[发明专利]结合变化信息和建筑特征的遥感影像新增建筑提取方法

权利要求书

1.一种结合变化信息和建筑特征的遥感影像新增建筑提取方法，其特征在于，包括以下步骤：

步骤1，采用SEEDS算法将后时相遥感影像分割，获得超像素图IM_seg；

步骤2，采用LSD算法对后时相遥感影像提取线段，得到线段集合L_set；

步骤3，基于线段集合L_set计算超像素S_i的边缘线段垂直性

步骤4，利用边缘线段垂直性和相交线段数计算建筑物线段指数

步骤5，利用形态学建筑指数MBI提取影像基于对象的形态学建筑特征

步骤6，融合形态学建筑特征和建筑线段指数计算影像的建筑物强度

步骤7，分别计算前后时相遥感影像的纹理特征和

步骤8，分别计算前后时相遥感影像的纹理特征差DT和光谱特征差DS；

步骤9，分别求纹理和光谱特征的变化强度IT和IS，通过直接求和获得变化强度CI；

步骤10，依据分割结果IM_seg求对象级的变化强度CI_seg，并与建筑物强度加权融合得到新增建筑强度图NBI；

步骤11，计算新增建筑强度图NBI的分割阈值T，对NBI二值化后获取新增建筑对象斑块NB；

步骤12，对每一新增建筑对象斑块NB分别计算形状指数GI，过滤出GI超过阈值的斑块NB，得到最终提取的新增建筑对象NB_post。

2.根据权利要求1所述的结合变化信息和建筑特征的遥感影像新增建筑提取方法，其特征在于，步骤3中，当有n条线段与超像素S_i相交时，边缘线段垂直性的表达式为：

式中，F(Δβ_jk)表示当两条线段l_j、l_k与超像素S_i相交时的垂直性值，Δβ_jk表示线段之间的夹角，θ表示倾斜角差容许误差。

3.根据权利要求2所述的结合变化信息和建筑特征的遥感影像新增建筑提取方法，其特征在于，步骤4中建筑物线段指数采用如下公式计算：

其中，N()表示线段包含的像素个数，N(S_i∩l_m)表示线段l_m与超像素S_i相交的像素个数。

4.根据权利要求3所述的结合变化信息和建筑特征的遥感影像新增建筑提取方法，其特征在于，步骤5中，基于对象的形态学建筑特征计算方式如下：

式中，MBI(x,y)表示形态学建筑指数MBI图的(x,y)处的象元值，p和q分别表示与超像素S_i等面积的矩形的长和宽。

5.根据权利要求4所述的结合变化信息和建筑特征的遥感影像新增建筑提取方法，其特征在于，步骤6中建筑物强度的计算方法如下：

式中，为可调节的加权因子，Nor()表示对建筑指数进行最大值最小值归一化处理。

6.根据权利要求5所述的结合变化信息和建筑特征的遥感影像新增建筑提取方法，其特征在于，步骤7中，采用GLCM的方差测度作为纹理特征，其计算方法如下：

其中，d∈{1,2}，P(i,j)代表灰度级别i和j在给定方向和距离情况下的出现的可能性，μ是GLCM的均值，N代表灰度级。

7.根据权利要求6所述的结合变化信息和建筑特征的遥感影像新增建筑提取方法，其特征在于，步骤8中，纹理特征差DT和光谱特征差DS的计算方法如下：

式中，S表示光谱波段，T₁和T₂分别表示遥感影像的两个时间，abs()表示求绝对值。