[發(fā)明專利]基于稀疏表示的SAR圖像相干斑抑制方法

有效

• 申請?zhí)枺?/td>: 201110346349.7
• 申請日: 2011-11-04
• 公開（公告）號: CN102346908A
• 公開（公告）日: 2012-02-08
• 發(fā)明（設計）人: 侯彪;焦李成;蔣繼光;王爽;劉芳;尚榮華
• 申請（專利權）人: 西安電子科技大學
• 主分類號: G06T5/00
• 分類號: G06T5/00;G06T5/50;G06T7/00
• 代理公司: 陜西電子工業(yè)專利中心 61205
• 代理人: 王品華;朱紅星

1.一種基于稀疏表示的SAR圖像相干斑抑制方法，包括以下步驟：

(1)對要進行相干斑抑制的SAR圖像Y進行對數(shù)變換，得到對數(shù)變換后的圖像Y′；

(2)將對數(shù)變換后的圖像Y′分割為大小為的重疊圖像塊；

(3)對分割后的每個圖像塊進行稀疏編碼，在編碼過程中引入該圖像塊的自身像素信息作為控制因子控制編碼：

(3a)將待編碼的每個圖像塊拉成列向量y，所有的列向量構成一個矩陣X；

(3b)選擇初始字典D為離散余弦變換DCT字典；

(3c)根據(jù)下列公式應用正交匹配追蹤算法OMP，實現(xiàn)對y的稀疏編碼：

α^=argminα||Dα-y||22+f(y)||α||0,]]>

其中α是列向量y在字典D下的稀疏表示向量，為迭代后的稀疏表示向量，||α||₀代表向量α的非零個數(shù)，為設定的懲罰項，用來控制編碼過程，是列向量y的歸一化L₂范數(shù)，σ為噪聲標準差，其計算公式為：

σ＝median(|W_Y|-median(W_Y))/0.6745，

其中W_Y表示原始SAR圖像Y經(jīng)過小波分解之后的高頻子帶系數(shù)；

(3d)對矩陣X中所有的列向量進行步驟(3c)的編碼操作，得到稀疏表示系數(shù)矩陣Λ；

(4)應用近似KSVD算法進行字典更新，得到更新后的字典D′和稀疏表示系數(shù)矩陣Λ′；

(5)利用更新后的字典D′和稀疏表示系數(shù)矩陣Λ′，得到處理后的圖像W；

(6)對處理后的圖像W進行指數(shù)變換，得到指數(shù)變換后的圖像R；

(7)將原始SAR圖像Y與指數(shù)變換后的圖像R相減得到差值圖像：V＝Y-R，對差值圖像V進行非線性各向異性擴散，進行k次迭代得到點目標圖像v^k；

(8)將點目標圖像v^k加到指數(shù)變換后的圖像R上，得到相干斑抑制后的圖像I＝R+v^k。