专利名称：基于各向异性扩散和互信息匀质性测度的sar图像降斑方法
基于各向异性扩散和互信息匀质性测度的SAR图像降斑方法技术领域
本发明属于图像处理技术领域，涉及一种SAR图像相干斑抑制方法，可用于对SAR 图像分割、分类和目标识别。
背景技术：
针对SAR图像的抑斑算法的基本目标是在抑制图像均勻区域斑点噪声的前提下， 保留图像的边缘、纹理和强反射点目标。针对SAR图像中的斑点噪声为乘性的特点，人们发展了许多的降斑算法。这类算法大致可以分为空域滤波算法和变换域滤波算法。典型的空域滤波算法有Lee滤波器，Kuan滤波器和Frost滤波器；典型的变换域滤波算法包括小波域、Contourlet域、curvelet域滤波算法。然而，不管是空域滤波算法，还是小波域滤波算法，它们在抑制相干斑时是各向同性的，即对边缘附近各个方向平滑程度相同，因此必然会破坏图像中的边缘和细节信息。基于非线性偏微分方程的图像滤波方法能够在降噪的同时增强图像的边缘。Yu ^ngjian等提出了一种各种异性扩散相干斑抑制算法—— SRAD (speckle reduction anisotropic diffusion)。SRAD 能够在均勻区域进行各向同性扩散，而在边缘附近进行各向异性扩散，因此可以在去除均勻区域相干斑噪声的同时增强和保护图像的边缘，但经过对SRAD的深入研究，发现其在均勻区域的降斑效果却不佳。主要是因为SRAD将变差系数作为勻质性测度来辨别图像中均勻区域、异质区域以及边缘，而变差系数极易受到相干斑噪声的影响，因此在均勻区域进行扩散时，有可能因为受到噪声影响而停止扩散，这样就使得均勻区域内噪声平滑不够充分。虽然可以通过增加迭代次数来解决，但随着迭代次数的增加图像中的边缘和纹理信息会最终而变得模糊。发明内容
本发明的目的在于克服SRAD的不足，提出一种基于互信息勻质性测度的各向异性扩散降斑方法，以提升各向异性扩散降斑方法在均勻区域的平滑效果并保留图像中的边缘和细节信息。
为实现上述目的，本发明的实现步骤如下
步骤1 根据图像的格式(强度与幅度)求解相干斑变差系数τ 0,当输入图像为强度格式时，Τ。=1/71 ；当为幅度图像时，Tfl = A^-IVZ ；
步骤2 计算图像中每一个像素点的互信息勻质性测度τ μ
步骤3 利用互信息勻质性测度构建各向异性扩散偏微分方程的扩散系数，即对像素点(i，j)，其对应的扩散系数表达式为 (、)=1/丨1 + [ -%]/[小1 +《)；1};
步骤4:基于步骤(3)得到的扩散系数，利用各向异性扩散偏微分方程 dI"j/dn = &v[q,y(、_)V仏]对图像中每一个像素点的灰度值进行迭代更新；
步骤5 更新迭代次数η :η = η+1，并判断η是否超过设定的最大迭代次数，若是输4出处理后的降噪图像；否则，返回步骤2。
本发明与现有技术相比的优点在于本发明相对于经典的空域滤波器如Lee滤波器、Frost滤波器及变换域滤波方法，其优势在于不仅可以有效去除SAR图像中均勻区域的相干斑噪声，还能够保留图像中的重要纹理信息并增强图像的边缘；相对于以往的各向异性扩散抑斑方法，本发明的优势在于可以增强均勻区域内对噪声的平滑效果。本发明的特点是由于使用了更加精确且对噪声不敏感的互信息勻质性测度作为SAR图像勻质区域、异质区域及边缘区域的检测算子，因此可以在勻质区域内进行充分的各向同性扩散从而尽可能地平滑掉噪声，在边缘与纹理区域进行各向异性扩散从而尽可能地保留原始图像信息。
图1是本发明对SAR图像进行降噪的流程图2是实验用的真实SAR图像以及本发明对原始图像的降噪图3是实验用的真实SAR图像以及本发明对原始图像的降噪图.具体实施方式
如图1所示，本发明的具体实施步骤如下
步骤1根据图像的格式(强度与幅度)求解相干斑变差系数τ 0,当输入图像为强度格式时，^ =1/λ/Ζ ;当为幅度图像时，τ。=；
步骤2计算图像中每一个像素点的互信息勻质性测度依据Bruno Aiazzi 等人于 2004 年发表在 IEEE Transaction on Geoscience and Remote Sensing 上的文章 《Information-Theoretic Heterogeneity Measurement for SAR Imagery》，τ “ j 是按照如下步骤进行计算的
2a)在一个3X3的滑动窗口内计算出像素点(i，j)的局部均值/\与标准差
权利要求
1.利用各向异性扩散和互信息勻质性测度的SAR图像降斑方法，其特征在于如下步骤步骤1 根据图像的格式(强度与幅度)求解相干斑变差系数τ ^，当输入图像为强度格式时，r·。=VyfL ；当为幅度图像时，r。二 V(4;r-1)/Z ；步骤2 计算图像中每一个像素点的互信息勻质性测度τ i,」；步骤3 利用互信息勻质性测度构建各向异性扩散偏微分方程的扩散系数，即对像素点(i，j)，其对应的扩散系数表达式为Ι(Tlj) = 1/(1 + [Tlj -T0]/[T02(1 + r02)]}；步骤4 基于步骤⑶得到的扩散系数，利用各向异性扩散偏微分方程 KJ丨dn = divIClj (TKJ WIIJ ]对图像中每一个像素点的灰度值进行迭代更新；步骤5 更新迭代次数η :n = η+l，并判断η是否超过设定的最大迭代次数，若是，输出处理后的降噪图像；否则，返回步骤O)。
2.根据权利要求1所述的方法，其中步骤( 所述的计算图像的互信息勻质性测度 τ (i，j)，依据 Bruno Aiazzi 等人于 2004年发表在 IEEE Transaction on Geoscience and Remote Sensing上的文章《Information-Theoretic Heterogeneity Measurement for SAR Imagery》，是按照如下步骤进行计算的2a)在一个3X3的滑动窗口内计算出像素点(i，j)的局部均值/\与标准差~ :/^pAsJpw ·’2b)计算出所有像素点的么和弋后，以么和弋为坐标轴，绘出该图像所有像素点对应的八和之的二维分布图；2c)对么和弋的二维分布图进行256X256的网格划分，得到二维直方图h(i，j)，并对 h(i，j)进行规范化；2d)对二维的离散直方图h(i，j)进行高斯平滑得到式与&的离散联合概率密度2e)根据树式，之)计算^与之的二维条件概率密度函数树之Ih2f)根据之IA5)计算像素点(i，j)的互信息勻质性测度τ (i，j) j) = log {ρ(σ( , j) I y))}其中，α为一常数，对强度图像《 = 1/VZ，对幅度图像α = 0.5227/VZ。
3.根据权利要求1所述的方法，其中步骤(4)所述的利用各向异性扩散偏微分方程对图像中每一个像素点进行灰度值更新，是按照以下步骤进行的3a)首先对各向异性扩散偏微分方程3/；, Idn = dVK^hJV/^]进行离散化，得到如下离散表达式其中ns表示像素点S的邻域，I ns|表示邻域内点的个数，P表示ns中点，At为时间步长。3b)利用前向差分代替步骤3a)中的离散偏微分方程中的梯度算子，得到各个像素点灰度值的更新方程
全文摘要
本发明公开了一种基于各向异性扩散和互信息匀质性测度的SAR图像相干斑抑制算法，主要解决现有多数相干斑抑制算法在去除匀质区域中的噪声的同时无法有效保留图像中的结构特征的缺陷。其实现过程是(1)对输入图像中的每个像素点，计算该点的互信息匀质性测度；(2)利用互信息匀质性测度构造各向异性扩散方程的扩散系数；(3)利用新建的扩散系数对均匀区域进行各向同性扩散，对边缘区域进行各向异性扩散；(4)对图像进行一定次数的迭代扩散后，输出最终的降噪图像。与现有的各向异性扩散降斑方法、空域及小波域降斑方法相比，本发明在增强均匀区域的噪声平滑效果的同时，还能有效保留图像中结构和细节特征，各方面的降噪性能提升显著，可用于SAR图像的分割、分类与目标识别。
文档编号G06T5/00GK102521811SQ20111042913
公开日2012年6月27日 申请日期2011年12月9日 优先权日2011年12月9日
发明者何友, 蔡复青, 袁湛 申请人:中国人民解放军海军航空工程学院