本发明专利技术公开了一种结合双树复小波与双变量模型的的ＳＡＲ图像降斑方法，主要解决现有ＳＡＲ图像降斑方法不能很好的抑制斑点噪声和丢失部分边缘和细节信息的问题。其实现过程为：对原始ＳＡＲ图像进行双树复小波分解，得到在各个尺度上分解系数的实部和虚部；利用非对数加性噪声模型，求解得到噪声系数的方差；根据利用局部邻域窗口，求得复小波系数实部边缘方差和虚部边缘方差；利用最大后验估计求得阈值收缩函数，对双树复小波分解的系数进行阈值收缩；对缩减后的系数运行双树复小波重构，得到降斑后图像。本发明专利技术具有有效去除ＳＡＲ图像斑点噪声和边缘保持性好的优点，可用于对边缘信息丰富、细节丰富的ＳＡＲ图像，特别是含有机场、跑道、道路的ＳＡＲ图像降斑。
【技术实现步骤摘要】
本专利技术属于图像处理
，涉及图像噪声抑制，具体地说是一种复小波域的 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&g本文档来自技高网...
【技术保护点】
一种结合双树复小波与双变量模型的的ＳＡＲ图像降斑方法，包括如下步骤： （１）对原始ＳＡＲ图像Ｉ进行双树复小波分解，得到在尺度ｊ上的分解复系数ｙ↓［ｊ］，其实部和虚部分别为ｙ↓［ｒ，ｊ］，ｙ↓［ｉ，ｊ］； （２）利用非对数加性噪声模型，求解各尺度上的噪声方差σ↓［ｎ］↑［２］； （３）利用局部邻域窗口，分别求解在尺度ｊ上复小波系数的实部边缘标准差σ↓［ｒ，ｊ］和虚部边缘标准差σ↓［ｉ，ｊ］； （４）分别对在尺度ｊ上的复系数的实部和虚部进行阈值收缩，求得估计的无噪小波系数ｗ＾↓［ｊ］； （４ａ）设任一尺度的复小波系数的实部和虚部近似满足如下分布： ｐ↓［ｙ］（ｙ）＝３／２πσ↑［２］.ｅｘｐ（－√３／σ√（ｙ↓［ｒ，ｊ］↑［２］＋ｙ↓［ｉ，ｊ］↑［２］）） 其中σ为复小波系数边缘标准差，ｙ↓［ｒ，ｊ］与ｙ↓［ｉ，ｊ］分别为在第ｊ尺度上复小波系数的实部和虚部； （４ｂ）求解尺度ｊ上无噪信号的最大后验估计，得到第ｊ尺度上复小波系数实部收缩函数＊↓［ｒ，ｊ］和虚部的收缩函数＊↓［ｉ，ｊ］分别为： ＊↓［ｒ，ｊ］＝ｓｏｆｔ（√（ｙ↓［ｒ，ｊ］↑［２］＋ｙ↓［ｒ，ｊ］↑［２］）－√３σ↓［ｎ］↑［２］／σ）／√（ｙ↓［ｒ，ｊ］↑［２］＋ｙ↓［ｒ，ｊ］↑［２］）.ｙ↓［ｒ，ｊ］ ＊↓［ｉ，ｊ］＝ｓｏｆｔ（√ｙ↓［ｉ，ｊ］↑［２］＋ｙ↓［ｉ，ｊ］↑［２］－√３σ↓［ｎ］↑［２］／σ）√（ｙ↓［ｉ，ｊ］↑［２］＋ｙ↓［ｉ，ｊ］↑［２］）.ｙ↓［ｉ，ｊ］ 其中ｓｏｆｔ（ｇ）定义为： ｓｏｆｔ（ｇ）＝＊＊＊ （４ｃ）求解在尺度ｊ上复小波实部阈值Ｔ↓［ｒ，ｊ］和虚部阈值Ｔ↓［ｉ，ｊ］中的较大值Ｔ↓［ｊ］ Ｔ↓［ｊ］＝ｍａｘ（Ｔ↓［ｒ，ｊ］，Ｔ↓［ｉ，ｊ］） 其中Ｔ↓［ｒ，ｊ］＝√３σ↓［ｎ］↑［２］／σ↓［ｒ，ｊ］，Ｔ↓［ｉ，ｊ］＝√３σ↓［ｎ］↑［２］／σ↓［ｉ，ｊ］，σ↓［ｒ，ｊ］，σ↓［ｉ，ｊ］分别为第ｊ尺度上复小波系数的实部边缘标准差和虚部边缘标准差； （４ｄ）在第ｊ尺度上利用下式进行阈值收缩，计算得到缩减后的复小波系数： ＊↓［ｊ］＝ｓｏｆｔ（｜ｙ↓［ｊ］｜－Ｔ↓［ｊ］）.ｅ↑［ｉ.θ（ｙ↓［ｊ］）］ 其中θ（ｙ↓［ｊ］）表示ｙ↓［ｊ］方向的弧度值； （５）对缩减后的系数＊↓［ｊ］运行双树复小波重构，得到降斑后图像＊。...
【技术特征摘要】
【专利技术属性】
技术研发人员：王爽，焦李成，李军，凤宏晓，侯彪，钟桦，缑水平，田小林，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：87[]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。