[發明專利]一種用于城市區域合成孔徑雷達圖像的噪聲抑制方法

在審

• 申請號: 201611039420.6
• 申請日: 2016-11-21
• 公開（公告）號: CN106780361A
• 公開（公告）日: 2017-05-31
• 發明（設計）人: 肖輝;嵇瑋瑋;姚崇斌;陳占勝;金戈;李曙光;張寧;孫竹
• 申請（專利權）人: 上海航天測控通信研究所
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 上海漢聲知識產權代理有限公司31236
• 代理人: 胡晶

1.一種用于城市區域合成孔徑雷達圖像的噪聲抑制方法，其特征在于，包括以下步驟：

步驟S1：對SAR圖像實部/虛部數據上滑動窗口內的樣本數據估計Alpha-Stable分布參數，記錄所有滑動窗口的參數α_i,j、γ_i,j，其中，參數α和γ的下標i、j表示滑動窗口的行、列數；

步驟S2：對SAR圖像幅度數據進行對數變換，在對數變換后的圖像上劃分各對應于實部/虛部的滑動窗口，計算各對應滑動窗口內的Alpha-Stable分布表達式：

$P_A\left(X\right)=e^{2X}{\∫}_0^{\mathrm{\∞}}u\exp (-{\mathrm{\γu}}^{\mathrm{\α}})J_0\left({\mathrm{ue}}^X\right)du$

其中，隨機變量X服從阿爾法穩定分布的概率密度函數，X為SAR圖像實部/虛部幅度值的自然對數，u∈[0,∞)，參數α∈(0,2]，決定該分布脈沖特性的程度，參數γ為分散系數，J₀(x)為第一類0階貝賽爾函數；

和噪聲分布表達式：

$p_A\left(N\right)=\frac{2L^L\exp (2NL-{\mathrm{Le}}^{2N})}{\mathrm{\Γ}\left(F\right)}$

其中，隨機變量N服從Nakagami分布的概率密度函數，N為SAR圖像幅度值對數變換后加性斑點噪聲，L是SAR圖像的視數，Γ(F)是一個均值為1、方差為的具有Γ分布的的概率密度函數，F為SAR圖像幅度值的乘性噪聲幅度；

步驟S3：根據貝葉斯公式推出

$P_{X|Y}\left(X\right|Y)=\frac{P_{Y|X}\left(Y\right|X)P_X\left(X\right)}{P_Y\left(Y\right)}=P_N\left(N\right)P_X\left(X\right)$

并由最大后驗概率準則：

$\hat{X}\left(Y\right)=\arg {\mathrm{maxP}}_N\left(N\right)P_X\left(X\right)$

即在X動態范圍內使P_X(X)*P_N(N)值最大的X_i作為滑動窗口中心象素的值；

步驟S4：將所有經最大后驗估計后的滑動窗口中心象素值，作指數變換得到濾波后幅度影像。