[發明專利]一種基于邊緣信息和反卷積的SAR圖像超分辨率方法

有效

• 申請號: 201410070411.8
• 申請日: 2014-02-28
• 公開（公告）號: CN103839238B
• 公開（公告）日: 2017-02-15
• 發明（設計）人: 侯彪;焦李成;牛志偉;王爽;張向榮;馬文萍;馬晶晶
• 申請（專利權）人: 西安電子科技大學
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 西安智萃知識產權代理有限公司61221
• 代理人: 李東京

1.一種基于邊緣信息和反卷積的SAR圖像超分辨率方法，其特征在于：包含以下步驟：?

步驟1：輸入低分辨SAR圖像為I_l，初始化超分辨結果高分辨率SAR圖像I_h＝0，其高度H_h與寬度W_h由H_h＝2×H_l－1和W_h＝2×W_l－1計算得出，其中H_l和W_l分別為I_l的高度和寬度；?

步驟2：利用步驟1中得到的初始高分辨率圖像I_h，將I_l的像素值通過一定的復制規則逐一隔點復制到I_h中；?

步驟3：利用步驟2中得到的I_h的已填充點的像素值，對I_h中部分未填充位置I_h2i+1,2j+1的邊緣方向性使用最大后驗概率進行估計；I_h2i+1,2j+1到四個最接近的已填充位置的像素值分別為I_h2i,2j、I_h2i,2j+2、I_h2i+2,2j和I_h2i+2,2j+2，計算I_h2i,2j、I_h2i,2j+2、I_h2i+2,2j和I_h2i+2,2j+2對應的四個方向灰度距離估計R₀、R₁、R₂和R₃；?

步驟4：計算I_h中未填充位置I_h2i+1,2j+1四個最接近已填充位置I_h2i,2j、I_h2i,2j+2、I_h2i+2,2j和I_h2i+2,2j+2灰度值的標準差σ；?

步驟5：利用R₀、R₁、R₂、R₃和σ計算對應I_h2i,2j、I_h2i,2j+2、I_h2i+2,2j和I_h2i+2,2j+2的權重系數w₀、w₁、w₂和w₃，其中k＝0，1，2，3;?

步驟6：利用步驟2中得到的I_h和權重系數w₀、w₁、w₂、w₃，計算在步驟二中未填充位置像素值I_h2i+1,2j+1＝w₀I_h2i,2j+w₁I_h2i,2j+2+w₂I_h2i+2,2j+w₃I_h2i+2,2j+2，并填充到對應位置；?

步驟7：將步驟6中得到的I_h，逆時針旋轉45°，擴展部分用零值填充，擴展后圖像矩陣為I'_h；對非擴展部分重復步驟3到步驟6過程；?

步驟8：將圖像I'_h順時針旋轉45°，保留I'_h中的非擴展部分，得到圖像I_h；?

步驟9：對圖像進行去卷積變換消除因邊緣插值引起的模糊；定義與圖像I_h同樣大小的全零矩陣P，用3×3的高斯點擴散函數填充矩陣P；?

步驟10：對矩陣P做快速二維傅里葉變換得到與I_h同樣大小的矩陣P_f，對I_h做快速二維傅里葉變換得到與I_h同樣大小的矩陣F_h，然后F_h中元素與P_f中對應位置元素做除法運算，重新賦值到矩陣F_h中，其計算公式為：F_hi，j＝F_hi，j/P_fi，j;?

步驟11：對矩陣F_h做快速二維逆傅里葉變換得到與I_h同樣大小的矩陣I_t，令矩陣I_h＝I_t，得到最終的超分辨結果I_h。?