[发明专利]一种用于城市区域合成孔径雷达图像的噪声抑制方法

权利要求书

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：将所有经最大后验估计后的滑动窗口中心象素值，作指数变换得到滤波后幅度影像。