一种基于多时相雷达图像纹理特征的震害建筑物识别方法技术

本发明专利技术属于SAR图像识别技术领域，特别涉及一种基于多时相雷达图像纹理特征的震害建筑物识别方法。基于多时相雷达图像纹理特征的震害建筑物识别方法，它包括如下步骤：A.将震前SAR图像I与震后SAR图像I'建立对应关系；B.主纹理分量计算；C.特征相关性计算与分析；D.相关性阈值统计与分类。本发明专利技术融合相关性分析与纹理特征主成分分量，实现震害建筑物的自动化提取，提取精度满足地震应急与评估的需求，提取结果可为地震应急救援力量分配、灾损评估、灾后重建等提供科学的信息依据。

A method for identifying earthquake damage buildings based on the texture features of multi phase radar images

【技术实现步骤摘要】
一种基于多时相雷达图像纹理特征的震害建筑物识别方法
本专利技术属于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'计算窗口大小，取奇数；Yls，Xls分别为震前SAR图像I与震后SAR图像I'第一主分量对应像元的特征值；分别表示震前SAR图像I与震后SAR图像I'第一主分量的区域特征平均值；相关系数r取绝对值，取值范围为[0,1]，r的绝对值越接近1，两幅图像的线性关系越密切，即相似程度越高，说明震后建筑物破坏程度越小；越接近0，两幅图像的线性关系越差，即差异性越大，说明震后建筑物的破坏程度越严重；D.相关性阈值统计与分类；D1.在震后获取的光学遥感影像中随机抽取不同破坏程度的建筑物样本点，根据样本点分布，统计在相关性图中样本点所在位置的相关系数的大小；D2.基于光学影像解译中建筑物样本点的破坏程度，对统计获得的样本点相关系数大小按照破坏程度进行分类，从而获取不同震害程度建筑物相关系数分类阈值范围。进一步的，利用步骤D2获得相关系数分类阈值范围后，将震后光学遥感影像第一分量的相关性图像分为三类，进而获得不同震害程度建筑物的空间分布。有益效果：本专利技术融合相关性分析与纹理特征主成分分量，实现震害建筑物的自动化提取，提取精度满足地震应急与评估的需求，相较于传统的SAR图像震害建筑物识别方法，该方法识别精度有显著的提高，提取结果可为地震应急救援力量分配、灾损评估、灾后重建等提供科学的信息依据。附图说明图1是本专利技术的技术流程图；图2为本专利技术实施例2中震害建筑物识别结果图。具体实施方式以下实施例用于说明本专利技术，但不用来限制本专利技术的范围。实施例1参见附图1，基于多时相雷达图像纹理特征的震害建筑物识别方法，它包括如下步骤：A.将震前SAR图像I与震后SAR图像I'建立对应关系；本例中，利用基于SIFT特征点的图像匹配方法，建立震前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'计算窗口大小，取奇数，m和n的取值与影像分辨率有关，一般取3、5或7；Yls，Xls分别为震前SAR图像I与震后SAR图像I'第一主分量对应像元的特征值；分别表示震前SAR图像I与震后SAR图像I'第一主分量的区域特征平均值；相关系数r取绝对值，取值范围为[0,1]，r的绝对值越接近1，两幅图像的线性关系越密切，即相似程度越高，说明震后建筑物破坏程度越小；越接近0，两幅图像的线性关系越差，即差异性越大，说明震后建筑物的破坏程度越严重；D.相关性阈值统计与分类；D1.对震后获取的光学遥感影像进行目视解译，在光学遥感影像中随机抽取不同破坏程度的建筑物样本点，根据样本点分布，统计在相关性图中样本点所在位置的相关系数的大小；D2.基于光学影像解译中建筑物样本点的破坏程度，对统计获得的样本点相关系数大小按照破坏程度进行分类，从而获取不同震害程度建筑物相关系数分类阈值范围。进一步的，利用步本文档来自技高网...

【技术保护点】
一种基于多时相雷达图像纹理特征的震害建筑物识别方法，其特征在于：所述识别方法包括如下步骤：A.将震前SAR图像I与震后SAR图像I'建立对应关系；B.主纹理分量计算；B1.统计分析得到震前SAR图像I与震后SAR图像I'中的两个或两个以上的纹理特征参数；B2.将震前SAR图像I中的两个或两个以上的纹理特征参数组合，计算多纹理特征的主分量，选择包含原纹理特征参数90％以上信息量的主分量作为震前SAR图像I的第一主分量；同理得到震后SAR图像I'的第一主分量；C.特征相关性计算与分析；设定计算窗口计算得到震前SAR图像I与震后SAR图像I'第一主分量的相关性图，相关性图中每一个像元的值代表相关系数r；

【技术特征摘要】
1.一种基于多时相雷达图像纹理特征的震害建筑物识别方法，其特征在于：所述识别方法包括如下步骤：A.将震前SAR图像I与震后SAR图像I'建立对应关系；B.主纹理分量计算；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'计算窗口大小，取奇数；Yls，Xls分别为震前SAR图像I与震后SAR图像I'第一主分量对应像元的特征值；分别表示震前SAR图像I与震后SAR图像I'第一主分量的区域特征平均值；相关系数r取绝对值，取值范围为[0,1]，r的绝对值越接近1，两幅图像的线性关系越密切，即相似程度越高，说明震后建筑物破坏程度越小；越接近0，两幅图像的线性关系越差，即差异性越大，说明震后建筑物的破坏程度越严重；D.相关性阈值统计...

【专利技术属性】
技术研发人员：李强，张景发，龚丽霞，薛腾飞，蒋洪波，罗毅，焦其松，
申请(专利权)人：中国地震局地壳应力研究所，
类型：发明
国别省市：北京,11