一种基于多特征融合的面向对象的建筑物变化检测方法
【技术领域】
[0001] 本发明属于遥感图像数据的信息提取技术领域，具体涉及一种基于多特征融合的 面向对象的建筑物变化检测方法。
【背景技术】
[0002] 随着经济发展和城市化进程的加快，人类的各种生产建设活动正在日益改变着城 市及其周边的自然环境和土地覆盖类型。因此快速有效的监测这些变化信息，分析变化原 因和影响结果，对我国的可持续发展具有重要意义。遥感图像变化检测技术就是根据同一 地区不同时相的两幅遥感图像，来获取地表地物变化的一门技术，该技术的快速发展为更 新地理空间数据发挥着重大的作用。遥感变化检测技术是数字图像处理方法、计算机视觉 技术和人工智能、模式识别理论的综合应用。
[0003] 遥感图像变化检测包括三个层次：像元级变化检测、特征级变化检测和目标级变 化检测。变化检测的三个层次各有优缺点，在具体的变化检测过程中要检测到哪个层次是 根据任务的需要来确定的。像元级变化检测保持了尽可能多的原始信息，具有特征级和目 标级层次上所不具备的细节信息，但像元级变化检测仅考虑像素属性的变化，而未考虑其 空间等特征属性的变化；特征级变化检测不仅考虑到空间形状的变化，而且还考虑特征属 性的变化，但特征级变化检测依赖于特征提取的结果，况且单独使用某一特征进行变化检 测有可能造成很大的漏检和错检；目标级变化检测最大的优点是它接近于用户的需求，检 测的结果可直接应用，它的不足之处在于目标提取的困难性。
[0004] 随着卫星技术的快速发展，使得高分辨率遥感图像的变化检测成为可能。分辨率 在10m以内的高分辨率遥感图像（Very High Resolution Imagery，VHR)已经广泛应用于 社会经济的很多领域，中低分辨率遥感图像中的点目标在高分辨率图像中变成了面，图像 包含更多、更丰富的地理和地形信息。高效利用高分辨率遥感图像丰富的地物细节变化信 息、像元之间以及像元属性之间的相互关系，能有效地抑制自然地物变化和不同成像条件 引起的信息干扰，因此对高分辨率遥感图像的变化检测研宄具有重要的理论意义和现实意 义。传统的遥感图像变化检测是基于像素的，如图像差值法、图像比值法。图像差值法和比 值法对图像的质量和预处理要求相对较高，不可避免的几何配准误差、相对辐射校正精度、 阴影等都是基于差异图像的高分辨率遥感图像变化检测中的典型问题。简单的单波段相减 由于没有考虑波段之间的统计相关性，使得绝对值不同的数值相减得到同样大小的差值， 不同像素灰度值之间相比得到同样大小的比值，忽略了不同地物在不同敏感波段存在的差 异，导致潜在可利用信息的丢失。
[0005] 在这背景下，面向对象的变化检测方法能够将像元-像元之间的差异推广到对 象-对象，把传统的以像素为单位的变化检测推广到以对象为单位的变化检测，最后提取 变化/未变化信息。
【发明内容】
[0006] 为了解决上述的技术问题，本发明提出了一种基于多特征融合的面向对象的建筑 物变化检测方法。
[0007] 本发明所采用的技术方案是：一种基于多特征融合的面向对象的建筑物变化检测 方法，其特征在于，包括以下步骤：
[0008] 步骤1 :对所选区域的两幅不同时相的高分辨率遥感图像A和B进行预处理，主要 包括几何纠正、辐射纠正、几何配准和图像裁剪；
[0009] 步骤2:计算A和B两幅图像的每个像素点的形态学建筑指数（MBI)、固定窗口的 纹理特征和慢特征分析（SFA)，得到图像A和B的MBI特征图像、纹理特征图像和SFA特征 图像；
[0010] 步骤3 :选取步骤1中的一幅图像采用基于改进的FNEA多尺度分割方法进行分 害J，得到一幅多尺度分割图像，并利用步骤2的计算结果来改进多尺度分割算法进行分割， 分割的结果是得到一个个对象；并由分割的结果得到每个对象对应像素点的坐标或索引， 称为索引矩阵；然后按照此索引矩阵对步骤1中的另一幅图像进行分割，分割的结果使两 幅图像具有同样大小的对象；
[0011] 步骤4 :因步骤2中得到的MBI特征图像、纹理特征图像和SFA特征图像具有尺度 不统一的特点，故采用单位标准差归一化方法对不同的特征图像进行特征优化；
[0012] 步骤5:对步骤4归一化后的结果，求每个对象的特征均值，以得到各个对象的特 征图像；
[0013] 步骤6 :对每个对象求其在不同时相的特征图像的差值，利用k均值聚类算法对三 幅不同的差值图像求其阈值，使其阈值自动化，以此得到三幅不同特征的变化图像；
[0014] 步骤7:对步骤6中得到的三幅不同特征的变化图像进行AC指数后处理；
[0015] 步骤8:利用熵值法，对步骤7中得到的三幅不同特征的变化图像加不同的权重， 再设定阈值，以得到变化检测结果；
[0016] 步骤9:利用基于投票法的多尺度融合方法对变化检测的结果进行处理，以得到 更尚的检测精度；
[0017] 步骤10:对步骤9得到的结果进行精度评定。
[0018] 作为优选，步骤1中所述的几何纠正采用基于多项式的遥感图像几何纠正，控制 点选取分布均匀，重采样采用双线性内插法，最后得到误差要求标准为RMSE〈0. 5像素；所 述的辐射纠正方法采用的是相对辐射归一化纠正。
[0019] 作为优选，步骤3中所述的选取步骤1中的一幅图像进行分割，选取方法为：当两 幅图像分辨率不一致时，用空间分辨率高的图像来进行分割；当分辨率相同时，按照获取时 相的时间顺序，选取后一时期的图像进行分割。
[0020] 作为优选，步骤3中所述的基于改进的FNEA多尺度分割方法，其具体实现过程为： 从一个像元起步，先将单像元合并为较小的对象，然后把具有异质性最小的较小对象合并 成较大的对象，这样不断合并，直到判断条件不成立，合并操作就终止，最终分割的结果中 所有图像对象的平均异质性最小；
[0021] 在判断两相邻区域是否能够合并时，用总异质性值和先前设定好的尺度阈值进行 比较，如果小于尺度阈值就合并，否则就结束合并操作；
[0022] 总的异质性h计算公式为：
[0023] h = wspectral*hspectral+wshape*hshape+w祖朴咖+界加咖朴如如；
[0024]其中，wspeetral, wshape, wmi，wtextme分别为光谱异质性hspeetral、形状异质性h shape、MBI 异 质性和纹理异质性h texture对应的权重；MBI异质性和纹理异质性的计算方法和光谱异 质性的计算方法原理相同，即计算每一波段的标准差与该波段权重的乘积，再把个波段的 值进行累加。
[0025] 作为优选，步骤6中所述的利用k均值聚类算法对三幅不同的差值图像求其阈值， 其具体实现包括以下子步骤：
[0026] 步骤6. 1:从数据集中随机选取K个图像单元作为初始聚类中心；
[0027] 步骤6.2 :计算各个图像单元到聚类中心的光谱距离，将它们一一归类到最近的 那个聚类中心所在的类；
[0028] 步骤6.3 :计算新形成的每一个聚类的图像单元的光谱均值，从而得到新的聚类 中心；
[0029] 步骤6. 4:迭代实施步骤6. 2和步骤6. 3,直至前后两次的聚类中心没有任何变化， 说明聚类调整结束，聚类准则函数已经收敛；
[0030] 当输入图像分别为归一化后的MBI特征图像、SFA特征图像和纹理特征图像时，用 K均值二值聚类所得的结果，即将这些特征分为两个分离度最高的图像类别，即低相似性与 高相似性类别；因此，由此得到的两个图像类别可分别对应于变化区域与未变化区域，即可 得到三个阈值X，Y，Z。
[0031] 作为优选，步骤7中所述的AC指数后处理，具体的公式如下：
[0032] A