[发明专利]基于双特征聚类同质像素选择的极化SAR非局部降噪方法

说明书

技术领域

本发明涉及图像滤波算法领域，具体是一种基于双特征聚类同质像素选择的极化SAR非局部降噪方法。

背景技术

相干斑是极化电磁波和物理散射体之间内在的相互作用。这种相互影响引起了像素与像素之间的变化，表现为强度值上的斑点噪声。因此，相干斑抑制对于极化SAR数据的图像分割和特征分类尤其重要。

过去几十年科研人员提出了多种滤波算法。简单常用的boxcar滤波器因为在窗口内不加区分的平均而效果不佳。极化白化滤波器只适用于单视极化SAR散射矩阵。Liu等人将其发展适合多视数据。Lee等人使用边界对齐窗和最小均方误差准则进一步减少了相干斑。强度驱动的自适应邻域滤波通过阈值来选择同质像素。基于散射模型的滤波根据相同的散射机制来选择同质像素。然而，这些滤波器都忽视了把结构特征和散射机制结合起来提高同质像素选择的准确性。

近年来，非局部均值滤波成为研究热点。像素之间的相似性度量准则已由原来的像素本身扩展为像素块邻域。因此，可使用定义在图像块上的统计特性来选择同质像素。Deledalle等人基于块的统计距离使块的相似性度量适用于极化SAR数据并且提出了统一的非局部滤波框架。Chen等人和Torres等人分别使用Wishart分布和stochastic散度来选择同质像素。Liu和Zhong在相似性度量时将无噪信号和相干斑噪声信号区分开。Zhong等人基于分布式Refined Lee滤波器和非局部均值滤波器将局部统计特性和非局部统计特性结合在一起。然而，这些算法都有一个共同的问题，那就是同质像素的选择都是基于有监督的块匹配方法，该方法难以确定合适的相似性阈值。再者，对于结构特征和散射机制保持，只使用统计相似性也是有缺陷的。本文解决的是如何以无监督的方式来选择结构和散射相似的同质像素。

聚类以无监督的方式比较像素与不同聚类中心的接近程度来自适应确定相似性，而不需要预先定义阈值。因此，就同质像素选择而言，聚类相比块匹配展现出了更好的结果。

发明内容

本发明的目的是提供一种基于双特征聚类同质像素选择的极化SAR非局部降噪方法，以提高同质像素选择的准确性，解决现有技术极化SAR图像结构信息、散射信息保持不足的问题。

为了达到上述目的，本发明所采用的技术方案为：

基于双特征聚类同质像素选择的极化SAR非局部降噪方法，其特征在于：包括以下步骤：

(1)、对每个像素点选择同质像素，具体步骤如下：

1a)、计算极化SAR协方差矩阵数据的Span图；

1b)、对包含d维空间的第i个图像块向量，在Span图上计算含有块的统计距离(PPB)的mean shift向量M_h,g(x_i):

式中：

P(x)＝c_p,dp(||x||²)，

g(x)＝-p′(x)，

x_i,k是第i个图像块矢量x_i的第k个数据，x_j,k是第j个图像块矢量x_j的第k个数据，p(||x||²)是高斯核函数，c_p,d是归一化常数，h是带宽；

1c)、在核函数窗口依次计算所有特征空间数据点的mean shift向量，沿mean shift梯度方向移动到新位置，然后将这个点作为一个新点，继续计算mean shift向量并移动，如果m_h,g(x_i)<δ，则停止聚类，m_h,g(x_i)就是该点处的模点；否则x_i＝m_h,g(x_i)且迭代计算M_h,g(x_i)；

1d)、每个极化SAR协方差矩阵进行Freeman-Durden分解：

F_S＝e_S(1+|β|²)，

F_D＝e_D(1+|α|²)，

F_V＝e_V，

式中e_S,e_D,e_V分别代表各自散射类别的贡献；β、α是参数。