[发明专利]一种对极化SAR数据相干斑噪声的抑制方法

权利要求书

1.一种对极化SAR数据相干斑噪声的抑制方法，包括如下步骤：

(1)将一组极化SAR数据表示为含有9个元素的3x3相干矩阵T，并使用T矩阵中的元素T11和T22对T矩阵进行亮目标检测，得到T矩阵的亮目标位置；

(2)对相干矩阵T元素的像素判定方向：

2a)定义八个方向窗w₁～w₈：

w1=-101-101-101]]>w2=011-101-1-10]]>w3=111000-1-1-1]]>w4=11010-10-1-1]]>

w5=10-110-110-1]]>w6=0-1-110-1110]]>w7=-1-1-1000111]]>w8=-1-10-101011]]>

2b)取相干矩阵T元素的一个像素x，以像素x为中心确定3x3区域A，将区域A分别与八个方向窗w₁～w₈相乘，得到八个乘积结果；

2c)在八个乘积结果中找出均值最大的乘积，此乘积所对应方向窗的方向就是像素x的方向d；

(3)以像素x为中心确定搜索窗Ω，在搜索窗内取一个像素y，判断y像素在方向d上是否有贡献，即像素y的区域更接近方向d对应的方向窗的分布：

3a)以像素y为中心确定3x3区域B，将区域B与方向d所对应的3x3方向窗相乘，得到一个乘积结果；

3b)对步骤3a)得到的乘积结果取均值，若该均值大于0，则像素y在方向d上有贡献，否则，像素y在方向d上无贡献；

(4)若像素y在方向d上无贡献，则在搜索窗Ω内取下一个像素重新判定，若像素y在方向d上有贡献，则将像素y保留，并使用像素y对像素x进行的非局部均值滤波，然后在搜索窗Ω内取下一个像素重新判定，最终记录搜索窗Ω内所有在方向d上有贡献的像素个数m，并得到像素x的非局部均值滤波结果x_NLmeans；

(5)对像素x进行改进的sigma滤波，得到改进的sigma滤波结果x_sigma；

(6)将步骤(4)和(5)得到像素x的滤波结果x_NLmeans和x_sigma按比例相加，得到结合的滤波结果

x^=(1-a)·xNLmeans+a·xsigma]]>

其中a为比例值，表示为a＝[1-m/(n×n)]²，(n×n)为搜索窗Ω的大小；

(7)对亮目标像素进行还原处理：

7a)如果像素x的位置为步骤(1)判定的亮目标位置，对得到的结合的滤波结果做如下处理：

x^′=0.7·x+0.3·x^]]>

为像素x的最终滤波结果；

7b)如果像素x的位置不为步骤(1)判定的亮目标位置，则像素x的最终滤波结果

(8)对相干矩阵T各元素的逐个像素进行上述步骤(2)～(7)的滤波过程，完成对相干矩阵T的滤波，得到滤波后的相干矩阵T；

(9)用Pauli向量法将滤波后的相干矩阵T合成伪彩图。