[發明專利]一種SAR圖像相干斑抑制方法

有效

• 申請號: 201710156351.5
• 申請日: 2017-03-16
• 公開（公告）號: CN106971382B
• 公開（公告）日: 2019-11-26
• 發明（設計）人: 王澤龍;余奇;譚欣桐;朱炬波;劉吉英;嚴奉霞;王丹
• 申請（專利權）人: 中國人民解放軍國防科學技術大學
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 11350 北京科億知識產權代理事務所(普通合伙)
• 代理人: 湯東鳳<國際申請>=<國際公布>=<進入

1.一種SAR圖像相干斑抑制方法，其特征在于，所述方法包括如下步驟：

S1、建立稀疏偏微分方程模型：

其中，(u-w)/u²為所述模型的逼近項；為所述模型的擴散項；u＝Ψc為所述模型的稀疏約束；w為觀測SAR圖像強度，u為待估SAR圖像強度；Ψ為過完備字典，c為稀疏表示系數，t表示時間，λ表示擴散步長；

S2、對所述稀疏偏微分方程模型進行數值求解，輸出稀疏表示系數c和待估SAR圖像強度u。

2.如權利要求1所述的SAR圖像相干斑抑制方法，其特征在于，所述步驟S2采用交替優化方法求解，具體包括以下步驟：

S2.1、固定變量u，變量c通過迭代求解：

其中，τ為非負參數，δ₁與δ₂為非負正則化參數，S(a,b)＝sgn(a)·max(|a|-b,0)，sgn為符號函數，上標H表示復共軛，k表示迭代次數，k為大于等于0的整數；初始值u⁰為原始圖像，c⁰滿足u⁰＝Ψc⁰；

S2.2、固定變量c，變量u通過迭代求解：

u^k+1＝u^k+ε[F(u^k)+2δ₁(u^k-Ψc^k)]

其中，ε為迭代步長，

S2.3、當||u^k+1-u^k||_*≤η時，迭代終止，輸出稀疏表示系數c和待估SAR圖像強度u；否則繼續執行步驟S2.1；其中η為迭代終止常數，其中所述||·||*表示核范數。

3.如權利要求1所述的SAR圖像相干斑抑制方法，其特征在于，所述過完備字典Ψ經過K-SVD方法獲得。

4.如權利要求2所述的SAR圖像相干斑抑制方法，其特征在于，λ、τ、δ₁、δ₂的取值均在0～20之間，ε的取值在0～1之間，η的取值在0～0.01之間。