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

有效

• 申請號: 201310423429.7
• 申請日: 2013-09-16
• 公開（公告）號: CN103440490B
• 公開（公告）日: 2016-10-19
• 發明（設計）人: 王亞飛;鐘禮山;李滿春;程亮;劉永學;陳振杰;杜陪軍;陳偉;諶頌;孫越凡
• 申請（專利權）人: 南京大學
• 主分類號: G06K9/46
• 分類號: G06K9/46
• 代理公司: 南京同澤專利事務所(特殊普通合伙) 32245
• 代理人: 石敏

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為絕對定標因子。