[发明专利]一种基于多时相雷达图像纹理特征的震害建筑物识别方法

说明书

技术领域

本发明属于SAR图像识别技术领域，特别涉及一种基于多时相雷达图像纹理特征的震害建筑物识别方法。

背景技术

现有的多时相SAR图像变化检测方法主要基于SAR图像的强度图像，不同时相强度图像之间易受SAR图像中固有的相干斑噪声及地表目标变化的影响，导致变化检测结果不能完全反映真实地物的变化。然而，城镇完好建筑物分布规则有序，通常在墙体、地面等形成大量的二面角结构，在SAR图像中表现为亮度高的点，纹理特征明显。在SAR图像形成过程中，二面角结构形成的角反射器受入射角与相干斑噪声的影响较少，纹理特征表现相对稳定，使得基于纹理特征的变化检测成为可能。同时在震区环境中，震害建筑物更为突出的是空间纹理特征而不是孤立的像素灰度信息，这也使得纹理分析成为震害图像解译的重要方法，现有基于纹理特征变化检测的方法通常是利用单一通道的纹理特征，未实现多特征融合与有效利用，造成潜在可利用信息的丢失。

发明内容

本发明的目的在于：针对强度图像易受斑点噪声、地物变化影响，利用单一纹理特征易造成可利用信息丢失的问题，提供一种基于多时相雷达图像纹理特征的震害建筑物识别方法。

本发明的技术方案为：基于多时相雷达图像纹理特征的震害建筑物识别方法，它包括如下步骤：

A.将震前SAR图像I与震后SAR图像I'建立对应关系；

B.主纹理分量计算；

主纹理分量计算是指通过主成分分析来分析获取多纹理特征的主成分分量；SAR图像具有丰富的纹理特征，纹理特征可描述地表的空间分布状态和粗糙程度，是判别地物的重要标志；纹理特征分析方法中，统计分析法是最常用的方法之一，该方法通过设置一定窗口大小单元统计分析获取纹理特征，获取的特征参数主要描述纹理局部模式的随机和空间统计特征，从而可用来表达区域的一致性和区域间的相异性；

B1.统计分析得到震前SAR图像I与震后SAR图像I'中的两个或两个以上的纹理特征参数；

B2.将震前SAR图像I中的两个或两个以上的纹理特征参数组合，计算多纹理特征的主分量；在获取主分量过程中，不需要参数设定和优化需求，所获得的主分量中的信息逐级减少，第一主分量包含了包含原纹理特征参数90％以上信息量，因此选择主成分分析后的第一主分量作为震前SAR图像I的第一主分量；

同理得到震后SAR图像I'的第一主分量；

C.特征相关性计算与分析；

设定计算窗口计算得到震前SAR图像I与震后SAR图像I'第一主分量的相关性图，相关性图中每一个像元的值代表相关系数r；相关系数的大小可反映地震前后图像中建筑物区域发生变化的程度；

式中：m和n为震前SAR图像I与震后SAR图像I'计算窗口大小，取奇数；Y_ls，X_ls分别为震前SAR图像I与震后SAR图像I'第一主分量对应像元的特征值；分别表示震前SAR图像I与震后SAR图像I'第一主分量的区域特征平均值；

相关系数r取绝对值，取值范围为[0,1]，r的绝对值越接近1，两幅图像的线性关系越密切，即相似程度越高，说明震后建筑物破坏程度越小；越接近0，两幅图像的线性关系越差，即差异性越大，说明震后建筑物的破坏程度越严重；

D.相关性阈值统计与分类；

D1.在震后获取的光学遥感影像中随机抽取不同破坏程度的建筑物样本点，根据样本点分布，统计在相关性图中样本点所在位置的相关系数的大小；

D2.基于光学影像解译中建筑物样本点的破坏程度，对统计获得的样本点相关系数大小按照破坏程度进行分类，从而获取不同震害程度建筑物相关系数分类阈值范围。

进一步的，利用步骤D2获得相关系数分类阈值范围后，将震后光学遥感影像第一分量的相关性图像分为三类，进而获得不同震害程度建筑物的空间分布。

有益效果：本发明融合相关性分析与纹理特征主成分分量，实现震害建筑物的自动化提取，提取精度满足地震应急与评估的需求，相较于传统的SAR图像震害建筑物识别方法，该方法识别精度有显著的提高，提取结果可为地震应急救援力量分配、灾损评估、灾后重建等提供科学的信息依据。

附图说明

图1是本发明的技术流程图；

图2为本发明实施例2中震害建筑物识别结果图。

具体实施方式

以下实施例用于说明本发明，但不用来限制本发明的范围。

实施例1