专利名称：：结合双树复小波与双变量模型的sar图像降斑方法
技术领域：
：本发明属于图像处理
技术领域：
，涉及图像噪声抑制，具体地说是一种复小波域的SAR图像降斑方法，可用于合成孔径雷达图像的斑点噪声的抑制。
背景技术：
：合成孔径雷达(SAR)是一种高分辨率成像雷达。它具有全天候、多极化、多视角、多俯角数据获取能力及对一些地物的穿透能力，不仅广泛地被应用在军事上，在农业、气象、地形地貌、灾情监控等民用上也有大量的应用。但由于SAR发射的是相干波，这些相干波经过与地物的相干作用，特别是地物的后向散射作用，使目标回波信号产生了衰减，这种衰减表现在图像上就是相干斑点噪声。因此如何抑制图像中的相干斑噪声，提高图像的解读能力以及获得更多的信息成为一个重要的问题。SAR图像降斑的首要目标是在滤除斑点噪声的同时，尽可能的保持图像的细节信息。斑点噪声是一种复杂的乘性噪声模型的信号。对于斑点噪声这种特殊的性质，在过去的近二十年中，人们提出了很多经典的SAR降斑方法，如Lee滤波，增强Lee滤波，Kuan滤波等等。这些方法是用一个已定义好的滤波器窗来估计局域斑点噪声的方差，并进行滤波处理，其结果通常过分的平滑了边缘细节信息，这些方法均在一定程度上收到了较好的效果。1995年，美国学者Donoho把小波理论引入到图像去噪中，提出了小波软阈值方法。小波软阈值方法是一种非线性的算法，依然存在破坏图像细节信息的问题，对图像的辐射特性保持也不好。离散小波变换由于具有缺少平移不变性与较差的方向选择性的缺点，最近英国学者Kingsbury等人提出了双树复小波变换，在图像去噪中的应用已初步显示出其明显的优势与离散小波变换相比，双树复小波变换由于其具有近似平移不变性和更多的方向选择性，有效的解决离散小波变换中出现的振铃效应。但是这种双树复小波变换降斑方法并没有充分考虑图像的几何特征以及SAR图像在复小波域的统计特性和系数之间的局部相关性，降斑后的SAR图像平滑区域的斑点噪声滤除不充分，同时图像的细节和边缘信息部分丢失。
发明内容本发明的目的在于克服上述已有技术的不足，提出一种结合双树复小波与双变量模型的SAR图像降斑方法，以充分滤除降斑后的SAR图像平滑区域的斑点噪声，完整的保留图像的细节和边缘信息。实现本发明目的的技术思路是结合双树复小波变换的平移不变性和多方向选择性，利用复小波系数实虚部双变量模型对双树复小波变换分解高频系数进行自适应萎缩，得到滤除斑点噪声、保留细节信息的SAR图像。其具体实现步骤包括如下(1)对原始SAR图像I进行双树复小波分解，得到在尺度j上的分解复系数yr其实部和虚部分别为Lpyi,(2)利用非对数加性噪声模型，求解各尺度上的噪声方差σ2；(3)利用局部邻域窗口，分别求解在尺度j上复小波系数的实部边缘标准差Onj和虚部边缘标准差Oy;(4)分别对在尺度j上的复系数的实部和虚部进行阈值收缩，求得估计的无噪小波系数火(5)对缩减后的系数火运行双树复小波重构，得到降斑后的图像/。本发明与现有技术相比具有如下优点1)本发明由于采用非对数加性噪声模型，利用这种模型可以避免在乘性模型转化为加性模型时，因取对数操作带来的对原图像的辐射特性保持的不足。因此能够更加充分的保留原图像的辐射特性。2)本发明由于利用复小波系数实虚部双变量模型，充分考虑到了SAR图像本身的方向特性与局部特性，保留了更加丰富的图像边缘和细节信息，充分滤除了SAR图像平滑区域的斑点噪声。3)仿真结果表明，本发明方法较其它几种现有的经典的SAR图像降斑方法，在平滑区域的平滑效果以及边缘保持能力方面都有显著的提高。图1是本发明的流程框图；图2是本发明与现有两种降斑方法应用于X波段幅度SAR图像的降斑仿真结果对比图；图3是本发明与现有两种降斑方法应用于Ku波段强度SAR图像的降斑仿真结果对比图。具体实施例方式参照图1，本发明的具体实现步骤如下步骤1，对输入SAR原始图像进行双树复小波分解。输入原始SAR图像记为I，该SAR图像本身就是被斑点噪声污染的图像，因此不需要像研究自然图像去噪一样，给原图加上一个随机噪声或者某种特性的噪声，可直接对该图像进行降斑处理，对输入的原始SAR图像I进行双树复小波分解，得到一个与一个低频图像和J个尺度，每个尺度有6个高频图像，在尺度j上高频图像的复小波系数记为LYj=yrjJ+i·Yijj(1)其中yy为复小波系数实部，Yi,j为复小波系数虚部。步骤2，利用非对数加性噪声模型，求解各尺度上的噪声方差σ2。(2a)利用非对数加性噪声模型，将输入的原始SAR图像表示为I=RX=X+(R-I)X=X+N(2)其中R代表相干斑，其均值为1，方差为σ〗，Χ代表地物的真实后向散射强度，N是将要滤除的加性噪声；从输入的原始SAR图像中的每一点取方形窗I(k)，窗口大小为kXk，计算每个局部方形窗的噪声方差权利要求一种结合双树复小波与双变量模型的的SAR图像降斑方法，包括如下步骤(1)对原始SAR图像I进行双树复小波分解，得到在尺度j上的分解复系数yj，其实部和虚部分别为yr，j，yi，j；(2)利用非对数加性噪声模型，求解各尺度上的噪声方差(3)利用局部邻域窗口，分别求解在尺度j上复小波系数的实部边缘标准差σr，j和虚部边缘标准差σi，j；(4)分别对在尺度j上的复系数的实部和虚部进行阈值收缩，求得估计的无噪小波系数(4a)设任一尺度的复小波系数的实部和虚部近似满足如下分布<mrow><msub><mi>p</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>3</mn><mrow><mn>2</mn><msup><mi>πσ</mi><mn>2</mn></msup></mrow></mfrac><mo>·</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msqrt><mn>3</mn></msqrt><mi>σ</mi></mfrac><msqrt><msubsup><mi>y</mi><mrow><mi>r</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup></msqrt><mo>)</mo></mrow></mrow>其中σ为复小波系数边缘标准差，yr，j与yi，j分别为在第j尺度上复小波系数的实部和虚部；(4b)求解尺度j上无噪信号的最大后验估计，得到第j尺度上复小波系数实部收缩函数和虚部的收缩函数分别为<mrow><msub><mover><mi>w</mi><mo>^</mo></mover><mrow><mi>r</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>soft</mi><mrow><mo>(</mo><msqrt><msubsup><mi>y</mi><mrow><mi>r</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mrow><mi>r</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup></msqrt><mo>-</mo><mfrac><mrow><msqrt><mn>3</mn></msqrt><msubsup><mi>σ</mi><mi>n</mi><mn>2</mn></msubsup></mrow><mi>σ</mi></mfrac><mo>)</mo></mrow></mrow><msqrt><msubsup><mi>y</mi><mrow><mi>r</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mrow><mi>r</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup></msqrt></mfrac><mo>·</mo><msub><mi>y</mi><mrow><mi>r</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow><mrow><msub><mover><mi>w</mi><mo>^</mo></mover><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>soft</mi><mrow><mo>(</mo><msqrt><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup></msqrt><mo>-</mo><mfrac><mrow><msqrt><mn>3</mn></msqrt><msubsup><mi>σ</mi><mi>n</mi><mn>2</mn></msubsup></mrow><mi>σ</mi></mfrac><mo>)</mo></mrow></mrow><msqrt><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup></msqrt></mfrac><mo>·</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow>其中soft(g)定义为<mrow><mi>soft</mi><mrow><mo>(</mo><mi>g</mi><mo>)</mo></mrow><mo>=</mo><mfencedopen='{'close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>g</mi><mo><</mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>g</mi></mtd><mtd><mi>g</mi><mo>≥</mo><mn>0</mn></mtd></mtr></mtable></mfenced></mrow>(4c)求解在尺度j上复小波实部阈值Tr，j和虚部阈值Ti，j中的较大值TjTj＝max(Tr，j，Ti，j)其中σr，j，σi，j分别为第j尺度上复小波系数的实部边缘标准差和虚部边缘标准差；(4d)在第j尺度上利用下式进行阈值收缩，计算得到缩减后的复小波系数<mrow><msub><mover><mi>w</mi><mo>^</mo></mover><mi>j</mi></msub><mo>=</mo><mi>soft</mi><mrow><mo>(</mo><mo>|</mo><msub><mi>y</mi><mi>j</mi></msub><mo>|</mo><mo>-</mo><msub><mi>T</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>·</mo><msup><mi>e</mi><mrow><mi>i</mi><mo>·</mo><mi>θ</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow></msup></mrow>其中θ(yj)表示yj方向的弧度值；(5)对缩减后的系数运行双树复小波重构，得到降斑后图像FDA0000032058360000011.tif,FDA0000032058360000012.tif,FDA0000032058360000014.tif,FDA0000032058360000015.tif,FDA0000032058360000021.tif,FDA0000032058360000022.tif,FDA0000032058360000024.tif,FDA0000032058360000025.tif2.根据权利要求1所述的SAR图像降斑方法，其中步骤(2)所述的求解各尺度上的噪声方差σ2，按如下步骤进行(2a)根据非对数加性噪声模型，对原图像中I的每一点取方形窗I(k)，窗口大小为kXk，求每个局部方形窗的噪声方差^l(k)=(^2Hk)+mIk))■丨(^2R+！)3.根据权利要求1所述的SAR图像降斑方法，其中步骤(3)所述的在尺度j上复小波系数的实部边缘标准差σ^和虚部边缘标准差σu，是分别对第j尺度上复小波系数的实部lt’j和虚部Ii,j中的每一个点取方形窗N(I)，窗口大小为1X1，分别求解在尺度j上复小波系数的实部边缘标准差Qy和虚部边缘标准差Oiij全文摘要本发明公开了一种结合双树复小波与双变量模型的的SAR图像降斑方法，主要解决现有SAR图像降斑方法不能很好的抑制斑点噪声和丢失部分边缘和细节信息的问题。其实现过程为对原始SAR图像进行双树复小波分解，得到在各个尺度上分解系数的实部和虚部；利用非对数加性噪声模型，求解得到噪声系数的方差；根据利用局部邻域窗口，求得复小波系数实部边缘方差和虚部边缘方差；利用最大后验估计求得阈值收缩函数，对双树复小波分解的系数进行阈值收缩；对缩减后的系数运行双树复小波重构，得到降斑后图像。本发明具有有效去除SAR图像斑点噪声和边缘保持性好的优点，可用于对边缘信息丰富、细节丰富的SAR图像，特别是含有机场、跑道、道路的SAR图像降斑。文档编号G06T5/00GK101980286SQ201010541989公开日2011年2月23日申请日期2010年11月12日优先权日2010年11月12日发明者侯彪,凤宏晓,李军,焦李成,王爽,田小林,缑水平,钟桦申请人:西安电子科技大学