专利名称：：基于邻域方向性信息的sar图像降斑方法
技术领域：
：本发明属于图像处理
技术领域：
，涉及图像噪声抑制，具体地说是一种基于非下采样contourlet域的SAR图像降斑方法，可用于合成孔径雷达图像的斑点噪声的抑制。
背景技术：
：合成孔径雷达(SAR)是一种高分辨率成像雷达。它具有全天候、多极化、多视角、多俯角数据获取能力及对一些地物的穿透能力，不仅广泛地被应用在军事上，在农业、气象、地形地貌、灾情监控等民用上也有大量的应用。但由于SAR发射的是相干波，这些相干波经过与地物的相干作用，特别是地物的后向散射作用，使目标回波信号产生了衰减，这种衰减表现在图像上就是相干斑点噪声。因此如何抑制图像中的相干斑噪声，提高图像的解读能力以及获得更多的信息成为一个重要的问题。SAR图像降斑的首要目标是在滤除斑点噪声的同时，尽可能的保持图像的细节信息。斑点噪声是一种复杂的乘性噪声模型的信号。对于斑点噪声这种特殊的性质，在过去的近二十年中，人们提出了很多经典的SAR降斑方法，如Lee滤波，增强Lee滤波，Kuan滤波等等。这些方法是用一个已定义好的滤波器窗来估计局域斑点噪声的方差，并进行滤波处理，其结果通常过分的平滑了边缘细节信息，这些方法均在一定程度上收到了较好的效果。1995年，美国学者Donoho把小波理论引入到图像去噪中，提出了小波软阈值方法。小波软阈值方法是一种非线性的算法，依然存在破坏图像细节信息的问题，对图像的辐射特性保持也不好。后来美国学者DoandVetterli提出了一种图像的多方向多尺度表示方法，称作Contourlet变换。Contourlet变换利用金字塔滤波器获得多尺度分解，利用方向滤波器组获得方向的合成。由于使用了几何信息，Contourlet变换在例如降噪和纹理恢复中的图像分析应用中，得到了比离散小波变换更好的效果。Contourlet变换经过下采样和上采样，并不是平移不变的，而平移不变性在图像分析应用例.如边缘检测、轮廓提取和图像去噪中是至关重要的。最近美国学者Omha等人提出了非下采样分解机制的Corvtourlet变换，这种变换具有完全的平移不变性，在图像去噪的应用中已初步显示出其明显的优势与UWT相比，非下采样Contourlet变换能够提供更好的图像恢复结果，在细节结构的保持方面更优于前者；同样，非下采样Contourlet变换的平移不变性也解决了基于Contourlet变换的去噪处理中出现的块状效应问题。但是这种基于非下采样Contourlet变换域降斑方法并没有充分考虑图像的几何特征以及SAR图像在非下采样Contourlet变换域的统计特性和方向性特征，因而降斑后的SAR图像平滑区域的斑点噪声滤除不充分，同时图像的细节和边缘信息部分丢失。
发明内容本发明的目的在于克服上述己有技术的不足，提出一种基于邻域方向性信息的SAR图像降斑方法，以充分滤除降斑后的SAR图像平滑区域的斑点噪声，完整的保留图像的细节和边缘信息。实现本发明目的的技术思路是将修正的对数高斯分布和混合指数分布引入到縮减因子的表示中，并利用方向性邻域模型结合非下采样Contourlet变换的高度冗余性和多方向选择性对非下采样Contourlet变换分解高频系数进行自适应萎縮，得到滤除斑点噪声、保留细节信息的SAR图像。其具体实现步骤包括如下(1)对原始SAR图像y进行非下采样Contourlet分解，得到在尺度_/上位置/处的分解系数",j;(2)将第y+l层的分解系数作为该层的无斑系数，求解第/层的系数二值掩码、；(3)找出第/层系数掩码为1和为0的系数，分别用修正对数高斯分布和混合指数分布逼近，得到概率密度分布的条件似然比^;(4)根据设定的方向性邻域模型计算分解系数的縮减因子先验比/7,，得到非下采样Contourlet分解系数縮减因子&;(5)应用縮减因子对非下采样Contourlet分解的系数进行縮减；(6)对縮减后的系数";j运行非下采样Contourlet重构，得到降斑后图像j)。本发明与现有技术相比具有如下优点1)本发明由于采用改进对数高斯分布和混合指数分布对掩码为1和掩码为0的非下采样Contourlet分解系数的直方图进行逼近，更加准确的捕捉了SAR图像在非下采样Contourlet变换域的统计特性，充分滤除了SAR图像平滑区域的斑点噪声。2)本发明由于将方向性邻域模型引入缩减因子的计算中，充分考虑到了SAR图像本身的方向特性，保留了更加丰富的图像边缘和细节信息。3)仿真结果表明，本发明方法较其它几种现有的经典的SAR图像降斑方法，在平滑区域的平滑效果以及边缘保持能力方面都有显著的提高。图l是本发明的流程框图2是本发明使用的邻域方向性模型示意图3是本发明求得的SAR图像高频子带系数为1的掩码以及概率分布拟合示意图4是本发明与现有两种降斑方法应用于X波段幅度SAR图像的降斑仿真结果对比图5是本发明与现有两种降斑方法应用于Ku波段强度SAR图像的降斑仿真结果对比图6是对图4和图5的局部放大图。具体实施例方式参照图1，本发明的具体实现步骤如下步骤l，输入SAR图像记为y，该SAR图像本身就是被斑点噪声污染的图像，因此不需要像研究自然图像去噪一样，给原图加上一个随机噪声或者某种特性的噪声，可直接对该图像进行降斑处理。步骤2，将图像j;经过非下采样Coritoiirlet分解，得到一个与原图像大小相同的低频图像和J个尺度，每个尺度有D个高频图像，高频图像在尺度y上位置/处的系数记为。(2.l)将输入图像经过非下采样金字塔分解，分别得到与原图大小相同的一个低频图像和一个高频图像；(2.2)将高频图像输入非下采样方向滤波器组，得到D个与原图像大小相等Z)个方向的高频图像；(2.3)将在步骤(2.1)中输出的低频图像再重复以上步骤(2.1)和(2.2)/次，最终得到非下采样Contourlet分解的一个低频图像和J个尺度，每个尺度有D个与原图像大小相同的高频图像，该高频图像在尺度y上位置/处的系数记为w,"步骤3，估计尺度_/上位置/处的系数、二值掩码^"。通过比较当前尺度/上位置/处的分解系数"^和较粗尺度_/+1上位置/处的分解系数《的点积与当前尺度层的噪声方差的值，得到二值掩码；"i"j.io;，(1)其中，(T/为噪声标准差，其表达式为-MW/a^nu-Me^fl^r^^l)/0.6745，、e{0,1}为系数",j的掩码S,，,=1表明",，;.为重要系数，主要对应SAR图像中的细节和有用信息；、=0表示"/;为不重要系数，主要对应SAR图像中的斑点噪声信息。步骤4，根据系数二值掩码A,.计算概率密度分布的条件似然比。概率密度分布的条件似然比^由掩码为1的系数概率分布函数/(^ll)和掩码为0的系数概率分布函数iO)两部分组成。第一部分掩码为1的系数的概率分布函数p(m,|1)，采用改进对数高斯概率密度函数进行拟合xo^V2;r(2)其中，x为非下采样Contourlet分解的高频系数值，^为",，,的幅值，即w,=|",」，p是In(m,)的均值，c^是ln(m,)的标准差。图3(a)，(b)所示为图4(a)中SAR图像进行非下采样Contourlet分解的第三个方向粗尺度与次粗尺度的掩码估计结果，图3(c)，(d)为掩码为1的概率分布拟合图，其中直线所示为掩码为1的系数的直方图，带星号的直线所示为改进对数高斯概率密度函数的拟合曲线，带叉号的直线为统计先验指导的降斑方法中应用的Gamma概率密度函数的拟合曲线，从图3中可以看出，本发明中的改进高斯概率密度函数比Gamma概率密度函数更好的逼近了掩码为1的系数直方图。第二部分掩码为0的系数的概率分布函数p(^|0)，采用混合指数概率密度函数拟合10)=，_v+w2e_v其中，；c为非下采样Contourlet分解的高频系数值，^为第一个指数分布尺度7参数，&为第二个指数分布的尺度参数，Wl和w2分别为这两个指数分布的权重参数，^,&，Wl,w2这些参数的值均采用期望最大化算法估计；这样，概率密度分布的条件似然比表示为掩码为1的系数概率分布函数与掩码为0的系数概率分布函数之比-^一,ll)/一,10)。(4)步骤5，根据方向性邻域模型计算系数縮减因子先验比。方向性邻域模型如图2所示，5(/)是以参考系数/为中心的各向同性邻域模型如图2(a)中所示，(/)是以参考系数/为中心的方向性邻域模型，图2(b)所示为4个方向的邻域模型04(/)，图2(c)所示为8个方向的邻域模型08(/)。《为方向性邻域模型内的系数掩码值，如果模型S(/)中的系数掩码值《全部为0，则系数縮减因子先验比为77,=1;否则根据各向异性邻域模型0/)(/)判断，如果模型(/)中任意一个方向的系数掩码全部为l，则系数縮减因子先验比为/^exp(r(2.Z)-l"，y为控制因子；若以上两个条件都不满足，则系数縮减因子先验比为化=exp{yZ(2;,-1)}。综上，系数縮减因子先验比表示为:1<eXp{y-(2.Z)-l)}ifUS;'=1(5)步骤6，计算非下采样Contourlet分解系数的縮减因子仏，进而得到縮减后的非下采样Contourlet分解系数";j。首先，根据步骤4和5中得到概率密度分布的似然比^和系数縮减因子先验比&,计算非下采样Contourlet分解系数的縮减因子^为&="/(1+")，(6)然后，通过系数縮减因子《,与非下采样Contourlet分解系数""的点积得到縮减后的非下采样Contourlet分解系数为(7)步骤7，将縮减后的系数";7进行非下采样Contourlet重构，最后得到降斑后的输出图像j)。(7.l)将縮减后当前尺度y的D个方向的高频系数输入非下采样重构方向滤波器组，得到一个与原SAR图像大小相等的高频图像；(7.2)将步骤(7.1)中得到的高频图像与步骤(2.3)中保留的低频图像一起输入到非下采样重构金字塔滤波器，得到重构的尺度/的低频图像；(7.3)将当前尺度变为_/-"，"-l…C/—1)，重复步骤(7.1)和(7.2)次，最终得到重构的降斑图像j)。本发明的效果可以通过以下仿真结果进一步说明。1.仿真内容应用本发明降斑方法、统计先验指导的降斑方法与Garama-MAP降斑方法，分别对一幅X波段幅度SAR图像和一幅Ku波段强度SAR图像进行降斑，并给出了等效视数ENL,图像均值M和标准差V作为客观评价标准来说明本发明应用于SAR图像降斑性能。ENL越高，说明平滑区域的降斑性能越好，降斑后的均值越接近原图像均值越好，降斑后标准差越低，说明平滑效果越好。仿真结果中还给出了本发明方法与统计先验指导的降斑方法应用于SAR图像的降斑结果细节放大图，对比结果进一步说明了本发明在细节以及边缘保持方面的优越性。仿真内容l，将现有的Gamma-MAP降斑方法、统计先验指导的降斑方法和本发明降斑方法应用于X波段的幅度SAR图像进行降斑；仿真内容2，将Gamma-MAP降斑方法、统计先验指导的降斑方法和本发明降斑方法应用于Ku波段的强度SAR图像进行降斑；仿真内容3，将统计先验指导的降斑方法和本发明降斑方法应用于X波段的幅度SAR图像与Ku波段的强度SAR图像进行降斑,得到降斑结果的局部放大对比图。2.仿真结果图4是仿真内容1的仿真结果。其中图4(a)为原SAR图像，图4(b)为Gamma-MAP方法降斑结果，图4(c)为统计先验指导的方法降斑图，图4(d)为本发明降斑结果图。图4(a)中所示矩形区域l、2为表1中计算ENL所需的同质区域，表l列出了仿真内容(1)的仿真结果所得的均值、方差和等效视数对比结果。表l:不同降斑方法客观指标评价等效视数ENL、均值M、标准差V9<table>tableseeoriginaldocumentpage10</column></row><table>从表l中可以看出，本发明获得了比现有其他方法更高的等效视数，降斑后的图像均值非常接近原图像均值，降斑后的标准差最低。因此本发明取得了比其他降斑方法更优的平滑降斑效果。图5是仿真内容2的仿真结果。其中图5(a)为原SAR图像，图5(b)为经典Garama-MAP方法降斑结果，图5(c)为统计先验指导的方法降斑图，图5(d)为本发明降斑结果图。图5(a)中所示矩形区域3、4为表2中计算ENL所需的同质区域。表2列出了仿真内容(2)的仿真结果所得的均值、方差和等效视数对比结果。表2:不同降斑方法客观指标评价等效视数ENL、均值M、标准差V<table>tableseeoriginaldocumentpage10</column></row><table>从表2中可以看出，本发明获得了比其他方法更高的等效视数，降斑后的图像均值非常接近原图像均值，降斑后的标准差最低。因此本发明取得了比其他降斑方法更优的平滑降斑效果。图6是仿真内容(3)的仿真结果。其中，图6(a)为图4(a)的图像局部放大图，图6(b)为图4(c)的图像局部放大图，图6(c)为图4(d)的图像局部放大图；图6(d)为图5(a)的图像局部放大图，图6(e)为图5(c)的图像局部放大图，图6(f)为图5(d)的图像局部放大图。从细节放大对比图可以看出，本发明得到的降斑图像能更好的保持原图像的细节和边缘信息，获得了非常好的直观效果。权利要求1.一种基于邻域方向性信息的SAR图像降斑方法，包括如下步骤(1)对原始SAR图像y进行非下采样Contourlet分解，得到在尺度j上位置l处的分解系数nl，j；(2)将第j+1层的分解系数作为该层的无斑系数，求解第j层的系数二值掩码(3)找出第j层系数掩码为1和为0的系数，分别用修正对数高斯分布和混合指数分布逼近，得到概率密度分布的条件似然比ξl；(4)根据设定的方向性邻域模型计算分解系数的缩减因子先验比ηl，得到非下采样Contourlet分解系数缩减因子ql；(5)应用缩减因子对非下采样Contourlet分解的系数nl，j进行缩减；(6)对缩减后的系数nl，jnew运行非下采样Contourlet重构，得到降斑后图像id="icf0002"file="A2009100228660002C2.tif"wi="4"he="4"top="114"left="167"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>2.根据权利要求1所述的SAR图像降斑方法，其中步骤(2)所述的求解第/层的系数二值掩码；/j，是根据当前层_/的分解系数",，;和较粗一层+1的分解系数","计算、=Jo，其中，o"乂为噪声方差，表达式为cr;:Me^m(ln,j-M^/fl"(n^.)|)/0.6745。3.根据权利要求1所述的SAR图像降斑方法，其中步骤(3)的实现如下(3a)找出第/层系数掩码为l所对应的系数附,"=1)，掩码为l的系数的直方图用如下条件概率分布表示—|1)=~~2。'2其中，x为非下采样Contourlet分解的高频系数值，p为掩码为l的系数绝对值取对数所得到的幅值ln(lm,"=1)|)的均值，ct,为幅值阔=1)|)的标准差；(3b)找出第/层系数掩码为0所对应的系数/,(;,=0)，掩码为0的系数的直方图用如下条件概率分布表示10)=，—v+w2e—一其中，A为第一个指数分布尺度参数，义2为第二个指数分布的尺度参数，W,和^分别为这两个指数分布的权重参数，这些参数4，;^，H^，K的值采用期望最大化算法估计；(3c)根据掩码为1的系数概率分布函数和掩码为0的系数概率分布函数，计算概率密度分布的条件似然比为4.根据权利要求1所述的SAR图像降斑方法，其中步骤(4)的实现如下(4a)设定以参考系数/为中心的邻域各向同性的邻域5(/)和各向异性的邻域(/)，构成方向性邻域模型，其中D为分解的总方向数；(4b)根据方向性邻域模型计算系数縮减因子先验比为<formula>formulaseeoriginaldocumentpage3</formula>exP{7Z(2;,_1)}e&e其中，g为方向性邻域模型内的系数掩码值；(4c)根据概率密度分布的条件似然比^和系数縮减因子先验比",，计算系数縮减因子为《/="/(l+")。5.根据权利要求1所述的SAR图像降斑方法，其中步骤(5)所述的用缩减因子对非下釆样Contourlet分解的系数w,，,进行自适应縮减，是通过系数縮减因子《与分解系数w,y的点积得到，艮P:<formula>formulaseeoriginaldocumentpage3</formula>全文摘要本发明公开了一种基于邻域方向性信息的SAR图像降斑方法，主要解决现有SAR图像降斑方法不能很好的抑制斑点噪声和丢失部分边缘和细节信息的缺点。其过程为对原图像做非下采样Contourlet分解，对每个高频子带系数估计二进掩码；然后根据所对应的不同的掩码用修正对数高斯分布和混合指数分布获得条件似然比的表达式；再根据邻域方向性模型获得的非下采样Contourlet变换的方向信息求先验比；最后根据似然比和先验比获得缩减因子对每个高频子带系数进行缩减，将缩减变化以后的系数再经过重构得到降斑后的输出图像。本方法具有有效去除SAR图像斑点噪声和边缘保持性好的优点，可用于对边缘信息丰富、细节丰富的SAR图像，特别是含有机场、跑道、道路的SAR图像降斑。文档编号G01S13/90GK101566688SQ20091002286公开日2009年10月28日申请日期2009年6月5日优先权日2009年6月5日发明者彪侯,卓刘,芳刘,朱虎明,焦李成,爽王,田小林,红赵,桦钟申请人:西安电子科技大学