[发明专利]超高分辨率合成孔径雷达图像斑点噪声滤波器构建方法

说明书

技术领域

本发明涉及合成孔径雷达图像处理领域，具体涉及一种超高分辨率合成孔径雷达图像斑点噪声滤波器构建方法。

背景技术

合成孔径雷达(SAR)作为一种有源微波遥感技术，具有全天候、全天时的工作能力。它可以在不同频段和不同极化下得到目标的较高分辨率图像，在军事应用和国民经济的各个领域中表现出巨大的潜力和应用前景。斑点噪声的存在使得合成孔径雷达的遥感图像不能正确地反映地物目标的散射特性，不仅增加了使用者判读的复杂性，而且容易造成误判，降低了SAR对目标探测、分析和解释能力；也大大地增加了SAR图像定量化应用的难度。另一方面，随着SAR图像应用的不断拓展，对其空间、辐射分辨率要求不断提高，对斑点噪声的抑制水平要求也越来越高。有效地抑制或消除SAR图像斑点噪声是提高雷达图像质量、使SAR图像数据得到成功应用的基础。

SAR图像的斑点滤波器已经发展了三十多年。诸多方法被提出，包括传统的自适应滤波如增强Lee滤波、Kuan滤波等，基于小波变换的多分辨率滤波，基于压缩感知的字典学习滤波，非局部均值滤波等，但仅有一些使用简单却性能有限的方法被广泛采用。当涉及到具有分米级或厘米级分辨率的超高分辨率合成空间雷达图像时，传统的斑噪模型如乘积模型等假设受到挑战，因为它们不能解释所观察到的脉冲非高斯散射。受到上述影响，现有斑点噪声滤波器难以取得满意效果，当采用小窗口尺寸时，斑点难以除尽，严重影响后续分割、分类效果；而大窗口尺寸将导致图像细节的严重损失，使高分辨率的优势化为虚有。

发明内容

为了解决上述问题，本发明公开一种超高分辨率合成孔径雷达图像斑点噪声滤波器构建方法，包括以下步骤：分区步骤，采用复广义高斯分布对超高分辨率合成孔径雷达图像进行拟合，获取拟合区域的形状因子，将所述超高分辨率合成孔径雷达图像分为亚高斯区域和高斯区域；以及去噪步骤，针对所述亚高斯区域采用加性噪声滤波器，针对所述高斯区域采用乘性噪声滤波器，实现对超高分辨率合成孔径雷达图像的自适应噪声去除。采用该构建方法可构建多种具有足够平滑和细节保持能力的超高分辨率SAR斑点噪声去除混合滤波器。

优选为，所述亚高斯区域对应结构性区域，所述高斯区域对应分布式地物。

优选为，所述复广义高斯分布的概率密度分布函数表示为

其中，μ表示均值，β表示形状因子，Γ()表示伽马函数，C表示协方差矩阵定义为

其中，H表示哈密顿转置，z_r和z_i分别代表复数z的实部和虚部。

优选为，区分所述亚高斯区域和所述高斯区域的阈值为，当形状因子β≥1时表示高斯分布，当形状因子β<1时表示亚高斯区域。

优选为，采用最大似然法来获取待测数据的形状因子。

优选为，所述最大似然法获取形状因子具体包括以下步骤：

通过以下公式获取初始形状因子

接着通过一个递归的操作对形状因子进行更新

L是(1)中概率密度函数的对数值，

当β<1，通过固定点方法更新C

C＝(2β·Γ(2/β)/2Γ(1/β)/N)∑(z^HC^-1z)^β-1zz^H(5)

并继续更新β。

优选为，所述加性噪声滤波器采用增强Lee滤波器，所述乘性噪声滤波器采用高斯低通滤波器。

优选为，所述增强Lee滤波器滤波方法具体表示如下

w(t)＝exp(-K(C_I-C_u)/(C_max-C_I))(7)

C_u＝0.523(8)

其中，F(t)代表位置t处的滤波结果，A(t)代表原始的影像强度或幅度，代表划窗中的强度或幅度均值，K代表阻尼系数，C_I代变异系数。

优选为，高斯低通滤波器以二维高斯函数作为核函数

其中，x和y代表像素位置，σ代表标准差。

上述加性噪声滤波器和乘性噪声滤波器仅为采用本发明构建混合斑噪滤波器的一个例子，任何类似的滤波器都可以应用，并构建新的滤波器。

本发明能够自适应地对结构性细节和分布式面状地物分开滤波，达到有效滤除斑噪同时良好保持影像细节的目的。