一种对极化SAR数据相干斑噪声的抑制方法.pdf

(19)中华人民共和国国家知识产权局*CN101976433A*(12)发明专利申请(10)申请公布号CN101976433A(43)申请公布日2011.02.16(21)申请号201010278044.2(22)申请日2010.09.09(71)申请人西安电子科技大学地址710071陕西省西安市太白南路2号(72)发明人王爽焦李成李悦侯彪钟桦于昕凤宏晓沈威(74)专利代理机构陕西电子工业专利中心61205代理人王品华朱红星(51)Int.Cl.G06T5/00(2006.01)G01S13/90(2006.01)权利要求书3页说明书7页附图4页(54)发明名称一种对极化SAR数据相干斑噪声的抑制方法(57)摘要本发明公开了一种对极化SAR数据相干斑噪声的抑制方法，主要解决现有极化SAR数据相干斑抑制方法不能很好的抑制相干斑噪声和丢失部分边缘和细节信息的缺点。该方法是基于改进sigma滤波和非局部均值的思想相结合的极化SAR相干斑抑制方法，其实现过程为：(1)输入极化SAR数据的相干矩阵T；(2)对相干矩阵T进行亮目标检测；(3)对像素使用方向窗判定方向并做改进的sigma滤波和非局部均值滤波；(4)滤波结果进行结合并保留亮目标；(5)完成T矩阵的滤波并合成伪彩图。本发明与现有技术相比显著提高了极化SAR数据的相干斑噪声抑制的能力，能够有效的平滑同质区域和保持边缘细节信息，可用于极化SAR数据的预处理过程。CN1097643ACCNN110197643301976438A权利要求书1/3页1.一种对极化SAR数据相干斑噪声的抑制方法，包括如下步骤：(1)将一组极化SAR数据表示为含有9个元素的3x3相干矩阵T，并使用T矩阵中的元素T11和T22对T矩阵进行亮目标检测，得到T矩阵的亮目标位置；(2)对相干矩阵T元素的像素判定方向：2a)定义八个方向窗w1～w8：2b)取相干矩阵T元素的一个像素x，以像素x为中心确定3x3区域A，将区域A分别与八个方向窗w1～w8相乘，得到八个乘积结果；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的非局部均值滤波结果xNLmeans；(5)对像素x进行改进的sigma滤波，得到改进的sigma滤波结果xsigma；(6)将步骤(4)和(5)得到像素x的滤波结果xNLmeans和xsigma按比例相加，得到结合的滤波结果其中a为比例值，表示为a＝[1-m/(n×n)]2，(n×n)为搜索窗Ω的大小；(7)对亮目标像素进行还原处理：7a)如果像素x的位置为步骤(1)判定的亮目标位置，对得到的结合的滤波结果做如下处理：为像素x的最终滤波结果；7b)如果像素x的位置不为步骤(1)判定的亮目标位置，则像素x的最终滤波结果(8)对相干矩阵T各元素的逐个像素进行上述步骤(2)～(7)的滤波过程，完成对相干2CCNN110197643301976438A权利要求书2/3页矩阵T的滤波，得到滤波后的相干矩阵T；(9)用Pauli向量法将滤波后的相干矩阵T合成伪彩图。2.根据权利要求书1所述的相干斑抑制方法，其中步骤(1)所述的用T矩阵中的元素T11和T22对T矩阵进行亮目标检测和保留，按如下步骤进行：2a)将相干矩阵T表示为：其中[ABC]＝2[Shh+SvvShh-Svv2Shv]，这里Shh表示h向发射和h向接收的回波数据，**Svv表示v向发射和v向接收的回波数据，Shv表示h向发射v向接收的回波数据，AA到CC*2*2分别代表T11到T33九个元素，T11＝AA＝|Shh+Svv|T22＝BB＝|Shh-Svv|；2b)将T11的所有像素从小到大排列，取出第t个像素，得到该像素值k1，z为T11像素总数；2c)使用3x3滑窗对T11的逐个像素进行扫描，当滑窗中9个像素中大于k1的个数超过TK时将此3x3区域视为亮目标区域，TK通常取5或6；2d)对T22做与上述步骤2b)-2c)相同的处理，得到T22的亮目标区域；2e)将T11和T22得到的亮目标区域的位置一起作为整个极化SA