[发明专利]一种极化SAR图像的快速滤波处理方法

说明书

技术领域

本发明涉及图像处理技术领域，特别是一种极化SAR图像快速滤波处理方法，可用于极化SAR图像的场景分析和图像理解的预处理。 

背景技术

合成孔径雷达(Synthetic Aperture Radar,SAR)就是利用雷达与目标的相对运动把尺寸较小的真实天线孔径用数据处理的方法合成一较大的等效天线孔径的雷达，也称综合孔径雷达。SAR不受天气、光照等条件的限制，能够连续、实时地进行对地观测。极化SAR为遥感图像中的信息获取和处理提供了更加丰富的目标信息，使其在采集地表或地面覆盖物的物理和电磁结构信息的应用中起着越来越重要的作用。由于极化SAR受相干斑噪声严重的干扰，相干斑噪声不仅导致数据质量下降，同时给极化SAR数据分类、目标检测等后处理带来困难。因此相干斑噪声的滤波处理成为极化SAR信息处理中的首要问题，相干斑噪声如果能够有效的抑制，可以显著提高极化SAR数据的可读性，有利于图像的处理、解译和信息提取。 

非局部(non local means，NLM)算法在近年来受到广泛关注和应用，其基本原理是利用图像块的相似性作为均值滤波的加权权重。这种算法对于自然图像的去噪结果十分理想，尤其是在边缘纹理等细节信息的保持方面明显优于其它现有的方法。目前该去噪方法正逐步推广到极化SAR相干斑噪声滤波处理领域。但是非局部均值算法的一个重要问题是计算量问题，在NLM类算法中像素点间相似性度量的计算最为耗时，因为任意两个图像块之间的相似性，都要通过计算对应像素之间的差异，然后进行累积获得。可以想象，随着滤波窗口的滑动，像素之间相似性度量运算存在大量的重复操作，使现有非局部均值类算法很难应用于工程实践，适用面小，计算速率低，费用高，为此需要提供极化SAR的非局部均值滤波类快速算法，但是今未见有公开报导。 

发明内容

针对上述情况，为克服现有技术之缺陷，本发明之目的就是提供一种极化SAR图像快速滤波处理方法，以拓展非局部均值滤波方法的适用范围，有效提高极化SAR图像滤波处理速度，节约费用的问题 

本发明解决的技术方案是，首先使用不局限于高斯分布的相似性度量计算滤波权重，然后利用积分图在极化SAR功率图上实现相似性度量的快速计算，最后用该相似度分别对极化散射矩阵或者协方差矩阵进行滤波处理，包括以下步骤： 

1)载入极化SAR图像数据，并进行扩展，为滤波做准备； 

2)利用扩展后的极化SAR图像计算得到总功率图像； 

3)记Δt表示目标像素t所在的邻域窗，Δs表示中心像素s所在的邻域窗，基于非局部均值中搜索窗Ω，利用积分图像计算实现快速滤波处理，步骤如下： 

3a)初始化邻域窗Δt在搜索窗内的当前位置 

3b)在扩展的总功率图像上，分别记A为邻域窗Δs中心像素构成的全局图像，B为邻域窗Δt中心像素构成的全局图像，按照下式计算全局图像A和B对应像素的相似性度量DS，其中I_A为图像A内像素的灰度值，I_B为图像B内像素的灰度值，(x,y)为像素坐标：