基于双特征聚类同质像素选择的极化SAR非局部降噪方法技术

本发明专利技术公开了一种基于双特征聚类同质像素选择的极化SAR非局部降噪方法，首先使用基于块的统计距离的mean shift聚类选择结构特征相似的像素，并且使用基于Freeman‑Durden散射模型的Wishart无监督分类器选择散射机制相似的像素。同质像素集合由同时含有相似结构特征和散射机制的像素构成；根据相似像素集合并结合原始极化SAR数据计算非局部权值加权求和实现极化SAR非局部均值滤波。本发明专利技术在有效抑制相干斑噪声的同时能进一步保持图像结构信息和极化散射信息，可用于极化SAR图像滤波处理。

Non local noise reduction method for polarimetric SAR based on dual feature clustering and homogeneous pixel selection

The invention discloses a non local denoising method of polarimetric SAR double feature clustering based on the selection of homogeneous pixels, first use pixel mean shift clustering statistical distance block selection based on similar structure features, and using the unsupervised classifier selection mechanism based on Freeman scattering pixel similarity Durden scattering model Wishart. Homogeneous pixel set pixel structure characteristics and scattering mechanism by simultaneously with similar composition; according to the similar pixel set and calculate the weighted sum of non local implementation of polarization SAR non local mean filter based on original polarimetric SAR data. The invention can further preserve the image structure information and polarization scattering information while effectively suppressing speckle noise, and can be used for polarization SAR image filtering processing.

【技术实现步骤摘要】
基于双特征聚类同质像素选择的极化SAR非局部降噪方法
本专利技术涉及图像滤波算法领域，具体是一种基于双特征聚类同质像素选择的极化SAR非局部降噪方法。
技术介绍
相干斑是极化电磁波和物理散射体之间内在的相互作用。这种相互影响引起了像素与像素之间的变化，表现为强度值上的斑点噪声。因此，相干斑抑制对于极化SAR数据的图像分割和特征分类尤其重要。过去几十年科研人员提出了多种滤波算法。简单常用的boxcar滤波器因为在窗口内不加区分的平均而效果不佳。极化白化滤波器只适用于单视极化SAR散射矩阵。Liu等人将其发展适合多视数据。Lee等人使用边界对齐窗和最小均方误差准则进一步减少了相干斑。强度驱动的自适应邻域滤波通过阈值来选择同质像素。基于散射模型的滤波根据相同的散射机制来选择同质像素。然而，这些滤波器都忽视了把结构特征和散射机制结合起来提高同质像素选择的准确性。近年来，非局部均值滤波成为研究热点。像素之间的相似性度量准则已由原来的像素本身扩展为像素块邻域。因此，可使用定义在图像块上的统计特性来选择同质像素。Deledalle等人基于块的统计距离使块的相似性度量适用于极化SAR数据并且提出了统一的非局部滤波框架。Chen等人和Torres等人分别使用Wishart分布和stochastic散度来选择同质像素。Liu和Zhong在相似性度量时将无噪信号和相干斑噪声信号区分开。Zhong等人基于分布式RefinedLee滤波器和非局部均值滤波器将局部统计特性和非局部统计特性结合在一起。然而，这些算法都有一个共同的问题，那就是同质像素的选择都是基于有监督的块匹配方法，该方法难以确定合适的相似性阈值。再者，对于结构特征和散射机制保持，只使用统计相似性也是有缺陷的。本文解决的是如何以无监督的方式来选择结构和散射相似的同质像素。聚类以无监督的方式比较像素与不同聚类中心的接近程度来自适应确定相似性，而不需要预先定义阈值。因此，就同质像素选择而言，聚类相比块匹配展现出了更好的结果。
技术实现思路
本专利技术的目的是提供一种基于双特征聚类同质像素选择的极化SAR非局部降噪方法，以提高同质像素选择的准确性，解决现有技术极化SAR图像结构信息、散射信息保持不足的问题。为了达到上述目的，本专利技术所采用的技术方案为：基于双特征聚类同质像素选择的极化SAR非局部降噪方法，其特征在于：包括以下步骤：(1)、对每个像素点选择同质像素，具体步骤如下：1a)、计算极化SAR协方差矩阵数据的Span图；1b)、对包含d维空间的第i个图像块向量，在Span图上计算含有块的统计距离(PPB)的meanshift向量Mh,g(xi):式中：P(x)＝cp,dp(||x||2)，g(x)＝-p′(x)，xi,k是第i个图像块矢量xi的第k个数据，xj,k是第j个图像块矢量xj的第k个数据，p(||x||2)是高斯核函数，cp,d是归一化常数，h是带宽；1c)、在核函数窗口依次计算所有特征空间数据点的meanshift向量，沿meanshift梯度方向移动到新位置，然后将这个点作为一个新点，继续计算meanshift向量并移动，如果mh,g(xi)<δ，则停止聚类，mh,g(xi)就是该点处的模点；否则xi＝mh,g(xi)且迭代计算Mh,g(xi)；1d)、每个极化SAR协方差矩阵进行Freeman-Durden分解：FS＝eS(1+|β|2)，FD＝eD(1+|α|2)，FV＝eV，式中eS,eD,eV分别代表各自散射类别的贡献；β、α是参数。1e)、使用Wishart无监督分类器对每个类别迭代合并成5个聚类：式中Ci和Cj分别表示各自聚类的平均协方差矩阵；类中心协方差矩阵Cc和像素点矩阵Cp之间的Wishart距离为1f)、在搜索邻域内选择同时含有相似结构特征和散射机制的像素构成同质像素集合；(2)、根据同质像素进行估计，具体步骤如下：2a)、根据第一步选择出的同质像素，计算原始极化SAR数据的权值w(i,j)：其中：式中H(k)是两个像素点的p×p协方差矩阵的相似性度量，S(i,j)是以像素点i，j为中心的图像块的协方差矩阵相似性度量；σ和Z分别是滤波参数和归一化因子，Q(i,k)意味着以i为中心的块的第k个像素的协方差矩阵，Q(j,k)同上，N是图像块的大小；2b)、根据步骤2a)得到的权值w(i,j)，对同质像素加权求和，得到中心像素的滤波估计值与已有技术相比，本专利技术的有益效果体现在：1)本专利技术推导出了基于Span图上PPB距离的meanshift聚类来选择结构特征相似像素并且通过基于散射模型的Wishart无监督分类器来自适应地将相似散射机制像素聚集在一起，解决散射机制保持不足的缺陷。克服了原有块匹配方法阈值选择带来的缺陷。2)本专利技术将结构特征和散射机制结合在一起共同选择同质像素参与滤波，能够有效抑制相干斑，没有出现人造假目标和模糊效应，且有效保持了边缘、点目标和散射机制。附图说明图1是本专利技术的基于双特征聚类同质像素选择的极化SAR非局部降噪方法流程图。图2是RefinedLee极化SAR滤波算法的结果图。图3是本专利技术中基于双特征聚类同质像素选择的极化SAR非局部降噪方法结果图。具体实施方式基于双特征聚类同质像素选择的极化SAR非局部降噪方法，包括以下步骤：(1)、对每个像素点选择同质像素，具体步骤如下：1a)、计算极化SAR协方差矩阵数据的Span图；1b)、对包含d维空间的第i个图像块向量，在Span图上计算含有块的统计距离(PPB)的meanshift向量Mh,g(xi):式中：P(x)＝cp,dp(||x||2)，g(x)＝-p′(x)，xi,k是第i个图像块矢量xi的第k个数据，xj,k是第j个图像块矢量xj的第k个数据，p(||x||2)是高斯核函数，cp,d是归一化常数，h是带宽；1c)、在核函数窗口依次计算所有特征空间数据点的meanshift向量，沿meanshift梯度方向移动到新位置，然后将这个点作为一个新点，继续计算meanshift向量并移动，如果mh,g(xi)<δ，则停止聚类，mh,g(xi)就是该点处的模点；否则xi＝mh,g(xi)且迭代计算Mh,g(xi)；1d)、每个极化SAR协方差矩阵进行Freeman-Durden分解：FS＝eS(1+|β|2)，FD＝eD(1+|α|2)，FV＝eV，式中eS,eD,eV分别代表各自散射类别的贡献；β、α是参数。1e)、使用Wishart无监督分类器对每个类别迭代合并成5个聚类：式中Ci和Cj分别表示各自聚类的平均协方差矩阵；类中心协方差矩阵Cc和像素点矩阵Cp之间的Wishart距离为1f)、在搜索邻域内选择同时含有相似结构特征和散射机制的像素构成同质像素集合；(2)、根据同质像素进行估计，具体步骤如下：2a)、根据第一步选择出的同质像素，计算原始极化SAR数据的权值w(i,j)：其中：式中H(k)是两个像素点的p×p协方差矩阵的相似性度量，S(i,j)是以像素点i，j为中心的图像块的协方差矩阵相似性度量；σ和Z分别是滤波参数和归一化因子，Q(i,k)意味着以i为中心的块的第k个像素的协方差矩阵，Q(j,k)同上，N是图像块的大小；2b)、根据步骤2a)本文档来自技高网...

【技术保护点】
基于双特征聚类同质像素选择的极化SAR非局部降噪方法，其特征在于：包括以下步骤：(1)、对每个像素点选择同质像素，具体步骤如下：1a)、计算极化SAR协方差矩阵数据的Span图；1b)、对包含d维空间的第i个图像块向量，在Span图上计算含有块的统计距离(PPB)的mean shift向量M

【技术特征摘要】
1.基于双特征聚类同质像素选择的极化SAR非局部降噪方法，其特征在于：包括以下步骤：(1)、对每个像素点选择同质像素，具体步骤如下：1a)、计算极化SAR协方差矩阵数据的Span图；1b)、对包含d维空间的第i个图像块向量，在Span图上计算含有块的统计距离(PPB)的meanshift向量Mh,g(xi):式中：P(x)＝cp,dp(||x||2)，g(x)＝-p′(x)，xi,k是第i个图像块矢量xi的第k个数据，xj,k是第j个图像块矢量xj的第k个数据，p(||x||2)是高斯核函数，cp,d是归一化常数，h是带宽；1c)、在核函数窗口依次计算所有特征空间数据点的meanshift向量，沿meanshift梯度方向移动到新位置，然后将这个点作为一个新点，继续计算meanshift向量并移动，如果mh,g(xi)<δ，则停止聚类，mh,g(xi)就是该点处的模点；否则xi＝mh,g(xi)且迭代计算Mh,g(xi)；1d)、每个极化SAR协方差矩阵进行Freeman-Durden分解：FS＝eS(1+|β|2)，FD＝eD(1+|α2)，FV＝...

【专利技术属性】
技术研发人员：周芳，刘留，杨学志，许开炜，
申请(专利权)人：合肥工业大学，
类型：发明
国别省市：安徽,34