SAR图像相干斑滤波算法及评价

1、SAR图像相干斑滤波算法及评价目前已有大量的雷达相干斑抑制算法，这些算法可分为成像前的多视平滑预处理和成像后的滤波两大类。而成像后的滤波又包括空域滤波和频域滤波两种。为了减少相干斑噪声，早期的方法是在SAR成像处理中，通过降低处理器带宽形成多视图子图像，然后对多视子图像进行非相干叠加来降低相干斑噪声。这种非相干叠加来降低斑点噪声的方法称为多视处理。多视处理通过牺牲SAR图像的空间分辨率为代价来对相干斑进行抑制，已不能满足空间高分辨率的要求。空域滤波方法是利用图像像素的空间相关性对相干斑进行滤波，一般是利用一个滑动窗口，然后对窗口内的像素进行加权得到窗口中心点的像素值。频域的方法主要是利用小波变

2、换，比较著名的有小波软阈值方法，基于小波变换和多尺度分析的滤波方法。本文研究SAR图像边缘检测，采用了局域统计自适应滤波算法，因为该方法考虑了图像的不均匀性，以局域的灰度统计特性为基础来决定参与滤波的邻域像素点及其权值，能在平滑噪声的同时较有效的保持明显的边缘，而且能通过参数控制来调整平滑效果和边缘保持效果之间的权衡。本文采用了增强Lee滤波算法，Kuan滤波算法，Frost滤波算法，最大后验概率(MAP)滤波算法，边缘保持最优化(EdgePreservingOptimizedSpeckle，EPOS)滤波算法等。1.传统滤波方法传统滤波算法包括均值滤波、中值滤波等。这类算法的特点是直接对图像

3、进行处理，没有考虑任何噪声模型，也没有考虑噪声的统计特性。这些算法实现起来比较简单，但效果不太理想。它们计算简单，速度快，均匀区域的斑点噪声去除效果较好。缺点是细节保持得不好，图像边缘变模糊，点目标损失大，随着处理窗口的增大，图像的整体模糊和分辨率下降更严重。正是由于这两种传统滤波算法不适合相干斑噪声的乘性特点，实际中较少采用。11均值滤波均值滤波是将平滑窗口内所有像元的灰度值进行平均计算，然后赋给平滑窗口的中心像元，其数学表达式为：TOC o 1-5 h zR二丄玄nDN(1)i，jn2i,ji=1j=1式中，R为滤波后中心元素灰度值，DN为滤波窗口内各个像元的灰度值，i,ji,j窗口大小为

4、nxn。1.2中值滤波中值滤波是一种非线性信号处理技术。它假设信号有极端的数值，即认为在平滑窗口内噪声是极大值或极小值。中值滤波将平滑区域内所有像素的中值作为平滑区域中心像元值。设DN为奇数项离散系列(i=l,2n-l,j=l,2n-l),DN为DN按TOC o 1-5 h zi,ji,ji,j大小重新排列的奇数项离散系列，则中值滤波的数学表达式为： HYPERLINK l bookmark2 R=DN(2)i,jn,n式中，R为滤波后的中心像元灰度值，DN为滤波前平滑模板内各个像元 HYPERLINK l bookmark0 i,ji,j的原始灰度值，DN为平滑模板内各个像元重新排列后的灰度

5、值，窗口大小为i,jnxn。2基于局域统计特性的自适应滤波算法自适应滤波是近30年以来发展起来的一种最佳滤波方法。它是在维纳滤波、Kalman滤波等线性滤波基础上发展起来的一种最佳滤波方法。由于它具有更强的适应性和更优的滤波性能，从而在工程实际中，尤其在信息处理技术中得到了广泛的应用。自适应滤波的研究对象是具有不确定的系统或信息过程。这里的“不确定性”是指所研究的处理信息过程及其环境的数学模型不是完全确定的。其中包含一些未知因素和随机因素。2.1自适应滤波原理自适应滤波器的原理如图1所示：图1自适应滤波原理图图中x(j)表示j时刻的输入信号值，y(j)表示j时刻的输出信号值，d(j)表示j时刻

6、的参考信号值或所期望响应信号值，误差信号e(j)为d(j)与y(j)之差。(4)自适应滤波器的滤波参数受误差信号e(j)控制，根据e(j)的值而自动调整，使之适合下一时刻的输入x(j+1),以便使得输出y(j+1)接近于所期望的参考信号d(j+1)。局域自适应滤波算法是在图像上取一个平滑窗口，以窗口内所有像素值作为滤波器的输入值进行处理，得到的结果作为窗口内中心元素的滤波值。而在平滑窗口内如何完成滤波运算，是这类相干斑抑制算法研究的核心内容。基于局域统计的自适应滤波算法，应满足以下条件才能适合于SAR图像的处理：不要求确知信号的统计模型；要达到保留边缘且加强细节；它可分为以局域统计特性(均值和

7、方差)为依据和以其统计分布为依据两大类：有很好的相干斑抑制效果；算法高效，使每一个像素都能在滤波窗口内独立进行。2.2Sigma滤波该算法建立在SAR图像的乘性噪声模型上，假设斑点噪声的分布为高斯分布,窗口内的像素灰度值与其中心像素的灰度值比较接近。其基本原理为：Sigma滤波器将2o范围内的像素进行平均，即可去除差别过大的象素的影响。我们知道，对于一维高斯分布，采样点落在2o区间的概率是93.5%。在窗口滤波过程中，只选取窗口内像素灰度值落在2o范围内的点，将它们的平均值作为中心像素灰度的估计，而其它变化显著的像素则被视作边缘而不做滤波处理。首先计算滤波窗口内各像元灰度的平均值g作为滤波中心

8、像元(i,j)的平均ij值；然后再求窗口内标准差b作为滤波中心像元点(i,j)的标准差，公式如下(设(/窗口为(2M+1)(2N+1)：gij=g(i,j)二1(2M+1)(2N+1)(3)二b(i,j)二(2M+1；(2N+1)送记Ij)-京k=j-Ml=i-NSigma滤波器的算法表达式如下:迥艺8gklkl(5)(6)R=k二iml二jn近乞8klk二iml二jng(1-2a)gg(1+2c)jFklijFotherwisea=/g(7)F孤立散射体不应受到斑点平滑的影响，为此设置阈值，如果范围内的象素数小于或等于K=(滤波窗口大小+1)/2，则以中心象素周围最近的四点象素平均值作为滤波

9、输出。Lee滤波及其增强算法Lee滤波基于完全发育的斑点乘性噪声模型，假定先验均值和方差可由均质区内计算局域的均值和方差来得到，它是使用滤波窗口内样本均值和方差的自适应滤波算法。该方法是以MMSE(最小均方误差)准则作为基础，是固定窗口中观察强度g和局部平均强度g的线性组合，是一个优化的线性滤波器。该方ij法是在图像上对每个像元逐个滤波移动的过程，局部统计量随着空间位置的改变而改变。首先计算窗口内各像元灰度的平均值g作为滤波中心像元(i,j)的平均值;ij然后再求窗口内标准差a作为滤波中心像元点(i,j)的标准差，公式如下(设窗口为(2M+1)(2N+1)：g=g(z,j)=另迟g(k,l)j

10、(2M+1)(2N+1)k=jMl=iN理現(k,l)才)(2M+1)(2N+1)ijk=jMl=iNa=a(i,j)=ij(8)(9)Lee滤波表达式为:g=gw+g(1w)Lee滤波算法是在均质区域的基础上推导得到的，但这一点事实上在真实的SAR图像中是不成立的。因此，Lee滤波方法对于在保持边缘等细节信息方面不是十分理想，但同质区则比较有效。针对Lee算法的缺陷，A.Lopes提出根据图像不同区域采用不同滤波器的方法。A.Lopes把一个图像分为三类区域：第一类是均匀区域，其中的相干斑噪声可以简单地用均值滤波平滑掉；第二类是不均匀区域，在去除噪声时应保留纹理信息，应用Lee滤波；第三类是

11、包含分离点目标的区域，滤波器应尽可能地保留原始值。具体思想如图2所示：图2增强滤波算法流程增强的Lee滤波采用以下准则：CC时，在滤波子窗口内取均值代替中心像素的值；TOC o 1-5 h zIuCCC时，保留该中心像素值。Imax以上区域的划分和准则同样适用于下文提到的增强的Kuan滤波和增强的Frost算法。增强的Lee滤波表达式为：厂g二gjijg=gw+g(1-w)jjg二g HYPERLINK l bookmark30 Jijijocijucocijmax(10)其中：w=1-C2/C2是Lee滤波的权函数;uIc,c为阈值;umaxC二0.5227/_bI为局部标准差；uvL;C=

12、3C；l为成像视数。maxuKuan滤波及其增强算法Kuan滤波算法假设噪声为与信号相关的加法噪声，然后运用最小方差估计获得固定窗口中观察强度g和局部平均强度g的线性组合。Kuan滤波器与Lee(/滤波器的区别在于用一个信号加上一个依赖于信号的噪声来表示乘性模型的相干斑噪声。该方法是在图像上对每个像元逐个滤波移动的过程，局部统计量随着空间位置的改变而改变，公式如下：gijbij=g(i,j)二=b(i,j)=1(2M+1)(2N+1)1(2M+1)(2N+1)艺迟g(k,l)k=j-Ml=i-N理为(伙,i)-)ijk=j-Ml=i-N(11)(12)bcijucbcijmax(13)为阈值;

13、b为局部标准差;IKuan滤波表达式为:g=gw+g(1-w)ijijKuan滤波算法与Lee滤波算法一样，存在着保持边缘等细节信息不佳的问题。因此，它也有对应的增强算法。A.Lopes提出的增强的Kuan滤波表达式为:g，=gijijg=gw+g(1w)TOC o 1-5 h zijijgg=g HYPERLINK l bookmark34 Jijij其中：g为平滑处理后的像元灰度值；g为平滑窗口中各像元的原始灰度ijii值；1C2/C2g为窗口内像元灰度平均值；w=u-是Kuan滤波的权函数;i1+C2uc,cumax=1+2/L；L为成像视数。umaxFrost滤波及其增强算法Frost

14、滤波算法假定斑点噪声是乘性噪声的前提下，并假设SAR影像是平稳过程，对影像进行滤波。Frost滤波器的冲激响应为一双边指数函数，近似为低通滤波器，其滤波器参数由图像局域方差系数决定。冲激响应的衰减快慢取决于局域方差系数的大小，与其成正比关系。Frost自适应滤波器是以权重M值为自适应调节参数的环形对称滤波器，其数学表达式为：gijIJA=bijA.Lopes提出的增强的Frost滤波表达式为:g=g；ijij(14)CCC时Imaxg=g；ijij(16)其中：g为平滑处理后的像元灰度值;ijg为平滑窗口中各像元的原始灰度值;ijg为窗口内像元灰度平均值;ijM为平滑窗口中各个对应像元的权重指

15、数;ijT为平滑窗口内中心像元到其邻像元的绝对距离;ijb为平滑窗口中像元值的方差；n2是平滑窗口的大小;ij;C=guij；C=、；1+2/L；L为成像视数。max3.2.6MAP滤波算法最大后验概率(MAP)滤波法是假设相干斑为乘性gamma分布，所以又称GammaMAP滤波器。在知道c的概率密度函数(ProbabilityDensityFunction，PDF)先验知识情况下，就能获取更多的信息。这就是根据先验分布和似然函数的MAP滤波方法。g=ijCC时g=g；ImaxijijCCC时uImaxgQL-l)+fg2QL-1)2+4aLggjjjj(LH1)2aICg(a-2)+tg2(

16、a-2)2+8aggg=_jjjjj2a(L=1)。(17)(18)(19)(20)其中：g为平滑处理后的像元灰度值;jg为平滑窗口中各像元的原始灰度值;jg为窗口内像元灰度平均值;ja=(112-C2)；Iucj为平滑窗口中像元值的方差;n2是平滑窗口的大小；C=&j;C=1;C=1+2/L；L为成像视数。I，-guLmaxj3.2.7EPOS滤波算法边缘保持最优化(EdgePreservingOptimizedSpeckle,EPOS)滤波算法，通过噪声方差的知识，区分均匀区域和含有边缘或单散射点的区域，提出了一种变大小滤波矩阵找到每个像素点的最大均匀区域，在此区域作平均，则达到最优滤波并

17、能保证边缘不变模糊。EPOS算法基于改变窗口形状的思想，它改变步骤是逐步剔除与中心像素不相同的邻域，然后应用最后剩余的邻域内点的统计量估算最后的结果。EPOS算法将窗口的邻域分为8个不相重合的部分，如图3所示。图3EPOS算法的中心象素邻域划分首先计算滤波窗口的相对标准差b和8个邻域的相对标准差G(i=0,1,i7),然后与SAR图像的相对标准差阈值g=0.52272/L比较(L为成像视数):如果gg，整个区域均值滤波；否则剔除8个邻域中b最大的那个邻域，重新计算剩余区域的相对标准差i进行比较，直至满足均值滤波的要求，即搜索出与中心象素处于相同区域的最大相同区域的象素，用搜索到的全部象素的平均

18、值作为滤波输出值。当8个邻域全部剔除时，将窗口的长度减小1,重新计算，直至窗口的长度缩小为3时，直接保留中心像元的灰度值。3.滤波效果评价对于滤波算法的效果评价，可以从两个方面进行，一是主观评价标准，即通过人眼视觉效果进行定性评价;二是客观评价标准，比较常用的有以下评价指标：均方误差、峰值信噪比、以及等效视数17等，这些指标从量化角度考虑滤波效果，能够辅助主观评价标准进行定量评价。主观评价标准：人眼视觉效果视觉是一种极为复杂和重要的感觉，人所感受的外界信息80%以上来自视觉。人眼对图像高度敏感，很多情况下需要目视判断图像的质量或目视解释处理结果。对于SAR图像，图像中目标的识别可以是机器识别，

19、不严格要求时效性的话，也可以是人眼判别，而且人眼判别具有更好的准确性。对于相干斑噪声的滤除效果，除了一些客观评价指标外，主观视觉判断是一条重要标准，毕竟人的知识、经验和判断能力都是计算机目前无法取代的。客观评价标准：均值(Mean)和方差(StandardDeviation，STD)图像均值是整个图像的平均强度，它反映了图像的平均灰度，即图像所包含目标的平均后向散射系数；图像方差代表了图像区域中所有点偏离均值的程度，反映了图像的不均匀性。图像的均值和方差是反映图像整体特征的指标，一般情况下，如果地形、含水量(复介电常数)和表面粗糙程度不同，则会有不同的后向散射系数，反映到SAR图像中有不同的图像均值。图像区域中的地形差异大,人工目标多，图像的灰值变化大，对应的图像的方差变化也就越大。所以应当尽量保持图像的均值，同时减少图像的方差。若图像区域大小为MxN，图像在(i,j)处的像素灰度值为I，则图像均值i,jMNi=1卩和图像方差b2分别定义如下：(21)(22)丄空(I-卩)2MNi，ji=1j=1均方误差(MeanSquareError，MSE)均方误差衡量滤波后图像和理想图像之间的差异程度。如果MSE值越小,则反映滤波后