[发明专利]SAR影像时空相似性分析下的城市不透水面提取方法

权利要求书

1.一种SAR影像时空相似性分析下的城市不透水面提取方法，包括以下步骤:

第一步、对SAR影像数据集进行预处理，所述SAR影像数据集包含有N幅不同时相的SAR影像；

第二步、对SAR影像数据集的N幅SAR影像进行几何配准，得到配准后的SAR影像数据集；

第三步、针对坐标值相同的SAR影像像元的后向散射系数，根据时间顺序建立与坐标值相关的像元后向散射系数时间序列；

第四步、在SAR影像数据集的画幅范围内目视选取M个城市不透水面纯净像元，相应的获得M个城市不透水面纯净像元坐标，对每幅SAR图像中M个城市不透水面纯净像元坐标位置的像元的后向散射系数求均值，得到SAR图像集的N个城市不透水面后向散射系数均值，根据时间顺序建立城市不透水面纯净像元的后向散射系数均值时间序列；

第五步、在SAR影像数据集的图幅范围内目视选取K个城市不透水面混合像元，相应的获得K个城市不透水面混合像元坐标，计算所述城市不透水面纯净像元的后向散射系数均值时间序列与每个城市不透水面混合像元坐标对应的所述像元后向散射系数时间序列之间的DTW距离值，并将计算得到的K个DTW距离值作平均，得到用于划分城市不透水面的最大阈值；

第六步、利用3*3的滑动窗口滑动扫描SAR影像数据集的整个图幅范围，当滑动窗口内所有像元对应坐标的像元后向散射系数时间序列与所述城市不透水面纯净像元的后向散射系数均值时间序列之间的DTW距离值均大于第五步得到的划分城市不透水面的最大阈值时，保持所述滑动窗口的中心坐标对应的像元后向散射系数时间序列不变，否则将所述滑动窗口的中心坐标对应的像元后向散射系数时间序列修改为所述城市不透水面纯净像元的后向散射系数均值时间序列；最终得到与坐标对应的修改后的像元后向散射系数时间序列；

第七步、再次利用3*3的滑动窗口滑动扫描SAR影像数据集的整个图幅范围，当滑动窗口内所有像元对应坐标的修改后的像元后向散射系数时间序列与所述城市不透水面纯净像元的后向散射系数均值时间序列之间的DTW距离值均小于第五步得到的划分城市不透水面的最大阈值时，所述滑动窗口的中心坐标对应的像元为城市不透水面，否则为非城市不透水面。

2.根据权利要求1所述的一种电力系统自动化实时多业务数据通信复用方法，其特征在于：第一步中对SAR影像数据集的预处理包括影像的裁剪和拼接、元数据更新、辐射校正、几何校正和地形校正。

3.根据权利要求1所述的一种电力系统自动化实时多业务数据通信复用方法，其特征在于，第二步中进行几何配准的方法为：选取SAR影像数据集中任一幅SAR影像作为主影像，其他影像为从影像进行几何配准。

4.根据权利要求1所述的一种电力系统自动化实时多业务数据通信复用方法，其特征在于：SAR影像中第i行第j列像元的后向散射系数

σij0=DNij2Ksin(αij)]]>

式中，DN_ij为第i行第j列像元的辐射强度值；α_ij为第i行第j列像元的的入射角；K为绝对定标因子。