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

摘要

本发明公开了一种用于城市区域合成孔径雷达图像的噪声抑制方法，该方法包括如下步骤对SAR图像实部/虚部数据上滑动窗口内样本数据估计Alpha‑Stable分布参数，记录所有滑动窗口的参数；对SAR图像幅度数据对数变换，在对数变换后图像上划分各对应于实部/虚部滑动窗口，计算各对应窗口内的Alpha‑Stable分布表达式和噪声分布表达式；计算Alpha‑Stable分布表达式与噪声分布表达式乘积的后验概率，取动态范围内最大后验概率的值作为滑动窗口中心象素的值；将所有经最大后验估计后的滑动窗口中心象素值作指数变换，得到滤波后幅度影像。本发明所述技术方案更有效地抑制城市区域SAR图像的斑点噪声，提高了SAR图像的质量。

主权项

一种用于城市区域合成孔径雷达图像的噪声抑制方法，其特征在于，包括以下步骤：步骤S1：对SAR图像实部/虚部数据上滑动窗口内的样本数据估计Alpha‑Stable分布参数，记录所有滑动窗口的参数αi,j、γi,j，其中，参数α和γ的下标i、j表示滑动窗口的行、列数；步骤S2：对SAR图像幅度数据进行对数变换，在对数变换后的图像上划分各对应于实部/虚部的滑动窗口，计算各对应滑动窗口内的Alpha‑Stable分布表达式：PA(X)=e2X∫0∞uexp(-γuα)J0(ueX)du]]>其中，随机变量X服从阿尔法稳定分布的概率密度函数，X为SAR图像实部/虚部幅度值的自然对数，u∈[0,∞)，参数α∈(0,2]，决定该分布脉冲特性的程度，参数γ为分散系数，J0(x)为第一类0阶贝赛尔函数；和噪声分布表达式：pA(N)=2LLexp(2NL-Le2N)Γ(F)]]>其中，随机变量N服从Nakagami分布的概率密度函数，N为SAR图像幅度值对数变换后加性斑点噪声，L是SAR图像的视数，Γ(F)是一个均值为1、方差为的具有Γ分布的的概率密度函数，F为SAR图像幅度值的乘性噪声幅度；步骤S3：根据贝叶斯公式推出PX|Y(X|Y)=PY|X(Y|X)PX(X)PY(Y)=PN(N)PX(X)]]>并由最大后验概率准则：X^(Y)=argmaxPN(N)PX(X)]]>即在X动态范围内使PX(X)*PN(N)值最大的Xi作为滑动窗口中心象素的值；步骤S4：将所有经最大后验估计后的滑动窗口中心象素值，作指数变换得到滤波后幅度影像。