[发明专利]一种新型纹理特征的SAR建筑物震害估计方法

权利要求书

1.一种新型纹理特征的SAR建筑物震害估计方法，其特征在于，其方法包括如下步骤：

(S1)首先将原始PolSAR数据导入进行预处理，利用改进的Yamaguchi四分量分解方法对全极化合成孔径雷达数据进行极化分解，提取二次散射分量和体散射分量；

(S2)将二次散射分量对应的全极化合成孔径雷达图像中的待识别地物，分类为平行完好建筑物；

(S3)根据公式计算提取的全极化合成孔径雷达数据的SPAN图像，并根据公式计算全极化合成孔径雷达数据的MSD值；

(S4)为MSD选择合适的阈值，根据公式所示的分类规则，将极化分解得到的体散射分量划分为倒塌建筑物和方位建筑物；

(S5)最后，将方位建筑物和平行完好建筑物归并为完好建筑物，将倒塌建筑物定义为损毁建筑物类别。

2.根据权利要求1所述的一种新型纹理特征的SAR建筑物震害估计方法，其特征在于：所述步骤(S1)中，全极化合成孔径雷达的总功率图像数据包含所有极化通道的强度信息。

3.根据权利要求1所述的一种新型纹理特征的SAR建筑物震害估计方法，其特征在于：所述步骤(S3)中，使用SPAN图像计算全极化合成孔径雷达数据的计算方法如下：SPAN＝HH²+HV²+VH²+VV²。

4.根据权利要求1所述的一种新型纹理特征的SAR建筑物震害估计方法，其特征在于：所述步骤(S3)中，全极化合成孔径雷达数据的MSD值的计算方法为灰度共生矩阵(GLCM)，基于其中的均值和方差，提出均值标准差纹理特征指数(MSD)，MSD的计算方法如下：

式中，mean(·)和variance(·)分别表示计算(·)的基于GLCM的平均值和方差值。

5.根据权利要求1所述的一种新型纹理特征的SAR建筑物震害估计方法，其特征在于：所述步骤(S4)中，MSD的分类规则表示为：x属于体散射占优的建筑物，如果MSD值大于等于ε，则判断x为倒塌建筑物，如果MSD值小于ε，则判断x为方位建筑物。

6.根据权利要求5所述的一种新型纹理特征的SAR建筑物震害估计方法，其特征在于：所述ε表示利用MSD参数划分倒塌建筑物和方位建筑物的阈值，ε的计算方法如下所示：

式中，mean(·)表示计算(·)的平均值，CB和OB分别代表倒塌建筑物样本集合和方位建筑物样本集合。