一种改进的PolSAR图像非局部均值滤波方法技术

一种改进的PolSAR图像非局部均值滤波方法。该方法结合图像子块的统计特性和目标点的极化散射特性筛选同质像素点，引入变异系数自适应选取平滑系数计算滤波所需的权重，对同质像素点进行非局部均值滤波，最后输出滤波后的PolSAR图像。本发明专利技术方法具有以下优点：1)可有效滤除相干斑噪声2)具有较强的自适应程度，本发明专利技术方法能够根据图像块的异质性来自适应选取平滑系数，在滤除噪声的前提下，更好地保持了图像的边缘特征3)利用极化筛选使滤波后图像点目标的极化散射特征得到了较好地保持。

【技术实现步骤摘要】
一种改进的PolSAR图像非局部均值滤波方法
本专利技术属于PolSAR图像解译
，特别是涉及一种PolSAR图像非局部均值滤波方法。
技术介绍
极化合成孔径雷达(polarimetricsyntheticapertureradar，PolSAR)具有全天时、全天候以及高穿透性的工作特性，在地物分类、积雪覆盖测量和灾情估计等遥感观测中起到了越来越广泛的作用。相干斑噪声是PolSAR成像系统所固有的，它的存在不仅会降低图像的视觉质量，而且也不利于分类、识别等对图像的后续处理。因此，相干斑抑制成为了PolSAR图像解译的前提条件。但目前尚缺少有效的方法。
技术实现思路
为了解决上述问题，本专利技术的目的在于提供一种在保证相干斑抑制效果的前提下，能够更有效地保持滤波后图像的纹理结构和目标点的极化散射特性的改进的PolSAR图像非局部均值滤波方法。为了达到上述目的，本专利技术提供的改进的PolSAR图像非局部均值滤波方法包括按顺序进行的下列步骤：1)输入PolSAR图像，确定该图像中待滤波像素点的搜索窗并在其中选取同质像素点，组成同质像素点集合；2)利用变异系数计算PolSAR图像中待滤波像素点的自适应平滑系数；3)对步骤1)获得的同质像素点集合进行非局部均值滤波，结合步骤1)得到的相似度和步骤2)得到的自适应平滑系数计算出每个像素点的权重，并将权重进行归一化处理，得到归一化权重；4)利用步骤3)得到的归一化权重对同质像素点进行加权平均，计算出滤波后的像素点的极化协方差矩阵，最后重构滤波后的PolSAR图像。在步骤1)中，所述的输入PolSAR图像，确定该图像中待滤波像素点的搜索窗并在其中选取同质像素点，组成同质像素点集合的方法是：首先输入PolSAR图像，根据所设搜索窗大小确定待滤波像素点的搜索窗，然后对搜索窗内的所有像素点进行极化筛选，选取满足极化筛选条件的像素点作为参考滤波像素点；计算每个参考滤波像素点与待滤波像素点之间的相似度，然后利用阈值判别法筛选像素点，满足筛选条件的像素点将组成同质像素点集合，不满足条件的像素点则舍弃。在步骤2)中，所述的利用变异系数计算PolSAR图像中待滤波像素点的自适应平滑系数的方法是：根据极化协方差矩阵C计算每个像素点的散射功率，由计算出的每个像素点的散射功率得到散射功率图，在散射功率图上计算出每个像素点的变异系数，利用变异系数计算出待滤波像素点的自适应平滑系数。在步骤3)中，所述的对步骤1)获得的同质像素点集合进行非局部均值滤波，结合步骤1)得到的相似度和步骤2)得到的自适应平滑系数计算出每个像素点的权重，并将权重进行归一化处理，得到归一化权重的方法是：在同质像素点集合中，利用相似度和自适应平滑系数计算每个像素点的权重，最后将这些权重累加并计算每个同质像素点的归一化权重。在步骤4)中，所述的利用步骤3)得到的归一化权重对同质像素点进行加权平均，计算出滤波后的像素点的极化协方差矩阵，最后重构滤波后的PolSAR图像的方法是：将待滤波像素点的同质像素点集合中每个像素点的极化协方差矩阵与其归一化权重相乘，最后再累加得到新的极化协方差矩阵；将待滤波像素点进行滑窗处理，直到遍历整幅PolSAR图像，得到PolSAR图像中每一个像素点滤波后的极化协方差矩阵，最后重构滤波后的PolSAR图像。本专利技术提供的改进的PolSAR图像非局部均值滤波方法与现有方法相比具有以下优点：(1)本专利技术方法用极化筛选优化相似像素点的选取规则，将PolSAR图像的统计特性和极化散射特性相结合，提高了相似像素点选取的准确性，同时也提高了滤波效率；(2)本专利技术方法利用图像块的变异系数得到可变的自适应平滑系数，使得其能够随着图像块的异质性自适应改变，在有效滤除噪声的基础上，可使PolSAR图像的边缘保持效果更佳。附图说明图1为本专利技术提供的改进的PolSAR图像非局部均值滤波方法流程图。图2为非局部均值滤波搜索窗示意图。图3(a)为SanFrancisco地区PolSAR数据PauliRGB图。图3(b)为SanFrancisco地区对应的光学图。图4(a)为HaywardFault地区的PauliRGB图。图4(b)为HaywardFault地区区域对应的光学图。图5(a)为SanFrancisco地区精致LEE滤波效果图。图5(b)为SanFrancisco地区NL-Pretest滤波效果图。图5(c)为SanFrancisco地区本专利技术方法滤波效果图。图6(a)为HaywardFault地区精致LEE滤波效果图。图6(b)为HaywardFault地区NL-Pretest滤波效果图。图6(c)为HaywardFault地区本专利技术方法滤波效果图。图7为SanFrancisco地区图像中点(500,495)的同极化响应图。其中，(i)原始点目标(ii)精致LEE滤波方法(iii)NL-Pretest滤波方法(iv)本方法。图8为SanFrancisco地区图像中点(500,495)的交叉极化响应图。其中，(i)原始点目标(ii)精致LEE滤波(iii)NL-Pretest滤波(iv)本方法。图9为HaywardFault地区图像中点(217,83)的同极化响应图。其中，(i)原始点目标(ii)精致LEE滤波(iii)NL-Pretest滤波(iv)本方法。图10为HaywardFault地区图像中点(217,83)的交叉极化响应图。其中，(i)原始点目标(ii)精致LEE滤波(iii)NL-Pretest滤波(iv)本方法。具体实施方式下面结合附图和具体实施例对本专利技术提供的改进的PolSAR图像非局部均值滤波方法进行详细说明。如图1所示，本专利技术提供的改进的PolSAR图像非局部均值滤波方法包括按顺序进行的下列步骤：1)输入PolSAR图像，确定该图像中待滤波像素点的搜索窗并在其中选取同质像素点，组成同质像素点集合；具体步骤如下：(a)首先输入PolSAR图像，如图2所示，设待滤波像素点为x，以该像素点为中心选取一定范围内的图像块作为待滤波像素点x的搜索窗S(x)；以待滤波像素点x为中心的3×3邻域像素块为待滤波像素点x的邻域窗，记为X；对搜索窗S(x)内像素点进行极化筛选：利用Yamaguchi分解方法将PolSAR图像的协方差矩阵C分解为以下四个分量：C＝[PSPDPVPC](1)其中，PS、PD、PV、PC分别表示表面散射功率、二面散射功率、体散射功率以及螺旋散射功率。选取这四个散射功率中的最大散射功率对应的散射方式作为该像素点的主导散射机制，判断搜索窗S(x)中的像素点和待滤波像素点x的主导散射机制是否相同，若相同，将该像素点视为待滤波像素点x的同质像素点并作为参考滤波像素点加以保留；(b)对于上一步保留的参考滤波像素点，如图2所示，以邻域窗为单位计算其与待滤波像素点x之间的相似度，公式如下：其中，C(Xi)、C(Yi)分别表示以待滤波像素点x和像素点y为中心的邻域窗中第i个像素点的协方差矩阵；k表示邻域窗中包含的像素点个数。然后利用阈值判别法对第i个像素点进行相似度筛选，公式如下：式中，H表示相似度阈值，计算公式为：式中，k表示邻域窗中包含的像素点个数，l表示视数，K表示调整参数。调整参数K较大时表本文档来自技高网...

【技术保护点】
1.一种改进的PolSAR图像非局部均值滤波方法，其特征在于：所述的改进的PolSAR图像滤波方法包括按顺序进行的下列步骤：1)输入PolSAR图像，确定该图像中待滤波像素点的搜索窗并在其中选取同质像素点，组成同质像素点集合；2)利用变异系数计算PolSAR图像中待滤波像素点的自适应平滑系数；3)对步骤1)获得的同质像素点集合进行非局部均值滤波，结合步骤1)得到的相似度和步骤2)得到的自适应平滑系数计算出每个像素点的权重，并将权重进行归一化处理，得到归一化权重；4)利用步骤3)得到的归一化权重对同质像素点进行加权平均，计算出滤波后的像素点的极化协方差矩阵，最后重构滤波后的PolSAR图像。

【技术特征摘要】
1.一种改进的PolSAR图像非局部均值滤波方法，其特征在于：所述的改进的PolSAR图像滤波方法包括按顺序进行的下列步骤：1)输入PolSAR图像，确定该图像中待滤波像素点的搜索窗并在其中选取同质像素点，组成同质像素点集合；2)利用变异系数计算PolSAR图像中待滤波像素点的自适应平滑系数；3)对步骤1)获得的同质像素点集合进行非局部均值滤波，结合步骤1)得到的相似度和步骤2)得到的自适应平滑系数计算出每个像素点的权重，并将权重进行归一化处理，得到归一化权重；4)利用步骤3)得到的归一化权重对同质像素点进行加权平均，计算出滤波后的像素点的极化协方差矩阵，最后重构滤波后的PolSAR图像。2.根据权利要求1所述的改进的PolSAR图像非局部均值滤波方法，其特征在于：在步骤1)中，所述的输入PolSAR图像，确定该图像中待滤波像素点的搜索窗并在其中选取同质像素点，组成同质像素点集合的方法是：首先输入PolSAR图像，根据所设搜索窗大小确定待滤波像素点的搜索窗，然后对搜索窗内的所有像素点进行极化筛选，选取满足极化筛选条件的像素点作为参考滤波像素点；计算每个参考滤波像素点与待滤波像素点之间的相似度，然后利用阈值判别法筛选像素点，满足筛选条件的像素点将组成同质像素点集合，不满足条件的像素点则舍弃。3.根据权利要求1所述的改进的PolSAR图像非局部均值滤波方法，...

【专利技术属性】
技术研发人员：韩萍，贾锟，卢晓光，张喆，韩斌斌，
申请(专利权)人：中国民航大学，
类型：发明
国别省市：天津,12