[發(fā)明專利]一種SAR圖像建筑物墻面窗戶信息提取方法

有效

• 申請?zhí)枺?/td>: 201910241719.7
• 申請日: 2019-03-28
• 公開（公告）號(hào): CN110276240B
• 公開（公告）日: 2021-05-28
• 發(fā)明（設(shè)計(jì)）人: 朱廈;劉鋒;康利鴻;馬喜樂;孫通;李志欣;王俊;薛涵;趙亮;范鵬程
• 申請（專利權(quán)）人: 北京市遙感信息研究所
• 主分類號(hào): G06K9/00
• 分類號(hào): G06K9/00
• 代理公司: 北京方韜法業(yè)專利代理事務(wù)所(普通合伙) 11303
• 代理人: 遆俊臣

1.一種SAR圖像建筑物墻面窗戶信息提取方法，具體步驟為：1)利用點(diǎn)特征提取方法得到SAR圖像中建筑物亮點(diǎn)；2)利用Radon變換得到SAR圖像中建筑物亮點(diǎn)對應(yīng)的方位向間距；3)利用紋理相關(guān)性最大法得到SAR圖像建筑物亮點(diǎn)對應(yīng)的距離向間距；4)利用Radon變換投影方差法得到SAR圖像中建筑物墻面方位角；5)基于SAR圖像建筑物亮點(diǎn)的Delaunay三角網(wǎng)，利用三角面元匹配法來提取準(zhǔn)確建筑物墻面窗口數(shù)量及空間分布；

其中，利用三角面元匹配法來提取準(zhǔn)確建筑物墻面窗口數(shù)量及空間分布的方法為：對于TRI中的第i個(gè)三角形Δ_i，計(jì)算其與基準(zhǔn)三角形Δ₀之間的三角形相似性SIM_i，如果SIM_i大于0，則三角形Δ_i中三個(gè)頂點(diǎn)均為建筑物墻面窗戶；滿足條件的所有點(diǎn)組成的點(diǎn)集W{p_i＝(x_i,y_i)|i＝1,2,…N_win}，其中，N_win為窗戶個(gè)數(shù)，窗戶空間分布為點(diǎn)集W中各元素的坐標(biāo)值(x_i,y_i)；其中，SAR圖像中建筑物墻面方位角，為建筑物墻面與SAR圖像方位向夾角；

所述的基準(zhǔn)三角形Δ₀為：其中一條邊長度l＝r，該條邊對應(yīng)的高h(yuǎn)_l＝a，該條邊對應(yīng)的夾角其面積其中，r為距離向間距，a為方位向間距，為墻面方位角。

2.如權(quán)利要求1所述的SAR圖像建筑物墻面窗戶信息提取方法，其特征在于，步驟1)所述的點(diǎn)特征提取方法為：對SAR圖像I_h×w中的每一個(gè)像元，其值大于其8鄰域中所有像元值且像元值大于0.02Q_max的點(diǎn)為建筑物亮點(diǎn)，Q_max為SAR圖像I_h×w的最大像素值；則建筑物亮點(diǎn)分布圖像BW_h×w是與SAR圖像I_h×w高度和寬度相同的數(shù)據(jù)矩陣，像素值為1表明是亮點(diǎn)，否則不是建筑物亮點(diǎn)。

3.如權(quán)利要求1所述的SAR圖像建筑物墻面窗戶信息提取方法，其特征在于，步驟2)所述的方位向間距計(jì)算方法為：首先，對建筑物亮點(diǎn)分布圖像BW_h×w進(jìn)行Radon變換，變換公式為其中，BW_h×w(x,y)表示點(diǎn)(x,y)處的像素值，δ()是Direac函數(shù)，ρ是圖像中心到直線的垂直距離，取值范圍1～N_ρ，間隔為1，θ是圖像中心到直線的垂線與x軸的夾角，取值范圍為0～180°，間隔0.1°；然后，θ取值為90°，得到一條隨ρ變化的曲線C(ρ)＝f(ρ,θ＝90°)，滿足條件diff(sign(diff(C)))＝-2的點(diǎn)為曲線C(ρ)的極大值點(diǎn)，對應(yīng)的橫坐標(biāo)系列{ρ_i|i＝1…m}，其中diff(C)為差分函數(shù)，sign()為符號(hào)函數(shù)；最后，方位向間距其中，mean[]是求和函數(shù)，差分函數(shù)diff(ρ_i)＝ρ_i-ρ_i-1。