专利名称:基于非下采样轮廓波的合成孔径雷达图像去噪方法
技术领域:
本发明属于图像处理技术领域,涉及该技术在合成孔径雷达Synthetic Aperture Radar 图像,即SAR图像去噪领域中的应用,具体地说是一种基于非下采样轮廓波的合成孔径雷 达图像去噪方法。该方法可用于SAR图像的预处理阶段。
背景技术:
合成孔径雷达可以全天时、全天候成像,且具有高空间分辨率和强穿透能力。SAR图 像在军事和民用方面均获得了广泛的应用。然而,受其自身固有的相干斑噪声的影响,SAR 图像不能有效地反映地物目标的散射特性,严重的影响了图像的质量,对图像的自动解译 造成很大困难。SAR图像去噪是SAR图像处理和分析的关键步骤,也是后续图像分割与 目标识别的基础。SAR图像去噪需要对同质区域有效滤波的同时,尽可能多得保持SAR 图像的方向以及纹理等重要信息,以利于后续解译任务的顺利完成。
SAR图像去噪方法大体可以分为两类。 一类是传统的基于空域的滤波方法,包括Lee 滤波、Gamma Map滤波和维纳滤波等。这些基于空域的滤波方法,难以保持图像的细节 特征,其去噪性能的好坏,很大程度依赖于所选滤波窗口的大小。
另一类是基于变换域的SAR图像去噪方法。由于小波变换具有良好的时频局部化特 性,所以基于小波变换的去噪方法在SAR图像去噪领域中得到了广泛的应用。基于小波变 换的SAR图像去噪方法分别在具有下采样分解机制和非下采样分解机制的两类小波变换 域中展开。在小波变换族中,非下采样小波变换解决了下采样小波变换缺乏平移不变性的 问题。实践证实了基于非下采样小波变换的去噪方法在综合性能上要优于基于下采样小波 变换的去噪方法。而不论是在下采样还是在非下采样的小波变换域对SAR图像进行去噪, 所得到的图像往往不能很好地保留原始SAR图像的细节特征,且小波分析在二维空间不是 最优的函数表示方法,不能很好地刻画图像中具有线奇异的几何信息。现有的基于变换域 的SAR图像去噪方法,多是将乘性斑点噪声简单通过对数变换转换为加性噪声来估计去噪 的阈值,导致去噪不充分。
近年来,人们提出新的多尺度几何分析工具,用来解决二维或是更高维奇异性。其中, 轮廓波Contourlet变换因对二维图像的表示性能优异、实现快速方便而备受瞩目。非下采
5样轮廓波变换Nonsubsampled Contourlet Transform, NSCT是轮廓波变换的平移不变版本,拥有比轮廓波变换更加丰富的基函数集合,可以提供更好的频率选择性和正则性,有利于更好地捕捉图像中的细节信息。
发明内容
本发明的目的在于克服已有技术的不足,即无法在有效去除SAR图像相干斑的同时保持图像的细节,提出一种基于非下采样轮廓波的SAR图像去噪方法及装置,以提高现有的基于变换域的降斑方法的细节保持性能。
本发明的技术方案是采用非下采样轮廓波作为对SAR图像变换的工具,并通过估计到的斑点噪声方差对变换系数进行判定,保留信号系数,抑制斑点噪声系数,实现SAR图像去噪,具体实现步骤如下-
(1) 输入SAR图像X,并对其进行£层非下采样轮廓波变换,得到低频子带Z),(a,6)和不同尺度上的高频方向子带{<^,(>》),OS/SZ-1, 1&、W, q,(a,6)表示SAR图像在尺度2-'上的第/个高频方向子带的系数,丄为3 5;
(2) 计算不同尺度各个高频方向子带的斑点噪声方差0^ ;
(3) 利用高频方向子带系数<:/,(",6)的局部均值附^"[(:,,^,6)]和高频方向子带的斑
点噪声方差cr^,将高频方向子带系数q,(",6)区分为信号或是噪声的变换系数,判定准则如下
其中,k在1 5之间取值,mea"[q,(a,州是以高频方向子带系数G,0,6)为中心的5x5大小的局部区域内系数的均值;
(4)对判定后的高频方向子带系数C,,,(a》)进行如下去噪操作
得到去噪后的高频方向子带系数《,(",6);
(5)对步骤(1)中得到的低频子带A(a,6)和去噪后的高频方向子带系数《,.(",6),作非下采样轮廓波逆变换,得到去噪后的SAR图像Z'。
上述步骤(2)中计算不同尺度各个高频方向子带的斑点噪声方差cr^ ,按如下步骤进
行
1)计算对SAR图像Z进行非下采样轮廓波变换时,施加在不同尺度各个高频方向子
(7,,0,6)€信号,如果 wea"[C,,(a,6)]2A;cj(C,, (fl, 6) e噪声,如果 me朋[C,, (a, 6)] <
C,.,0,6), 如果 C,^,6)e信号0, 如果 <:,,(",6)£噪声带的等价滤波器操作算子的二阶累积量《(2),
其计算公式为
,ri(W)2)f[(W)2)
其中,A表示尺度2-'上的第&个高频方向子带,//厂表示第/级塔形滤波的非下采样
拉普拉斯滤波器系数,^/ ]表示第y级方向滤波的非下采样方向滤波器系数;
2) 计算不同尺度各个高频方向子带的均值/^和方差CT^ ;
3) 计算不同尺度各个高频方向子带的斑点噪声方差C^,其计算公式为-
0> 二-:;-Cr7
其中,G7是输入SAR图像X的规范化标准差,且对于强度SAR图像有C^ = VTT^ ,
对于幅度SAR图像有C^ =A/(4;r-1)/i , i 为已知输入SAR图像的视数,i 21 。本发明提供的基于非下采样轮廓波的合成孔径雷达图像去噪的虚拟装置,包括变换装置,用于输入SAR图像X,并对其进行丄层非下采样轮廓波变换,得到低频
子带Z^(a,6)和不同尺度上的高频方向子带(q,(a,6), 0S/S丄—1, 1S,^W , C,,(a,6) g
SAR图像在尺度2—'上的第/个高频方向子带的系数,Z为3~5;
计算装置,用于计算不同尺度各个高频方向子带的斑点噪声方差o^ ;
判定装置,用于利用高频方向子带系数<^,(",的的局部均值,""0,6)]和高频方向
子带的斑点噪声方差o^ ,将高频方向子带系数C,,,(o,ZO区分为信号或是噪声的变换系数,
判定准则如下 2 A:crCsLc,, (a, &) e噪声,如果 we朋[C,, (", 6)] < Actcij
其中,k在1~5之间取值,me朋[q,0,6)]是以高频方向子带系数C^(a,6)为中心的5><5大小的局部区域内系数的均值;
去噪装置,用于对判定后的高频方向子带系数C,,,("力)进行如下去噪操作-
,fC,,.("》),如果 <:,,(",6)£信号
" 丄0, 如果 C,^,。e噪声
得到去噪后的高频方向子带系数《,&,6);
逆变换装置,用于对低频子带",(",w和去噪后的高频方向子带系数《,(",&),作非下
7采样轮廓波逆变换,得到去噪后的SAR图像Z'。
上述的非下采样轮廓波的合成孔径雷达图像去噪虚拟装置中的计算装置,包括
二阶累积量计算子装置,用于计算对SAR图像I进行非下采样轮廓波变换时,施加
在不同尺度各个高频方向子带的等价滤波器操作算子的二阶累积量《(2),
s('(2)=:rt )fi(2>")2),
J=lW 户l V
其中,it,表示尺度2-'上的第A个高频方向子带,/^']表示第/级塔形滤波的非下采样拉普拉斯滤波器系数,W"表示第y级方向滤波的非下采样方向滤波器系数;
高频方向子带均值和方差计算子装置,用于计算不同尺度各个高频方向子带的均值//义和方差o^ ;
斑点噪声方差计算子装置,用于计算不同尺度各个高频方向子带的斑点噪声方差7&,
其中,C^是输入SAR图像X的规范化标准差,且对于强度SAR图像有C^二^T^,对于幅度SAR图像有Gz =^(4丌-1)/i ,及为已知输入SAR图像的视数,i 2 1 。本发明与现有技术相比具有以下优点
1、 本发明由于釆用非下釆样轮廓波变换,可以有效避免因变换工具缺乏平移不变性而在去噪图像中产生的抖动失真。
2、 本发明可以针对SAR图像的斑点噪声在不同高频方向子带的统计特性,对SAR图像进行自适应地去噪。
3、 仿真结果表明,本发明相对于现有的SAR图像去噪方法,在抑制斑点噪声的同时可以有效保持图像中的细节信息。
本发明的技术过程和效果可结合以下附图详细说明
图1是本发明的实现过程示意图;图2是本发明的虚拟装置示意图3是用本发明和已有方法对SAR图像iml去噪结果图;图4是用本发明和已有方法对SAR图像im2去噪结果图。
具体实施例方式
参照图l,本发明的具体步骤如下
8步骤l,输入SAR图像X,并对其进行丄层非下采样轮廓波变换。
非下采样轮廓波变换是一种新的具有平移不变性的多尺度、局域的、多方向的过完备图像表示方法。非下采样轮廓波变换的构造基于非下采样的塔形滤波器组和非下采样的方向滤波器组,两部分之间是相互独立的。对SAR图像X进行一层非下采样轮廓波变换,
其过程为
1) 将SAR图像Z输入非下采样的塔形滤波器组,得到SAR图像X的一层非下采样轮廓波变换的低频信号和带通信号;
2) 将SAR图像I的带通信号输入非下采样的方向滤波器组,得到SAR图像1的一层非下采样轮廓波变换的高频方向子带,高频方向子带的数目可以是2的任意次幂;
3) 将SAR图像I的非下采样轮廓波变换的低频信号作为新的输入原图像,重复上述步骤l)和2),得到SAR图像Z进行i:层非下采样轮廓波变换的低频子带D,(a,6)和各个尺度上的高频方向子带(q,(a,6), 1, 1&、W, (^,,(",6)表示原图像1在尺度2"上的第/个高频方向子带的系数,Z为3 5;
步骤2,计算不同尺度各个高频方向子带的斑点噪声方差0^ ;
1)计算对SAR图像X进行非下采样轮廓波变换时,施加在不同尺度各个高频方向子带的等价滤波器操作算子的二阶累积量《(2),
其计算公式为
<formula>formula see original document page 9</formula>
其中,^表示尺度2—'上的第^个高频方向子带,/^']表示第/级塔形滤波的非下采样拉普拉斯滤波器系数,t/^表示第/级方向滤波的非下采样方向滤波器系数;
2) 计算不同尺度各个高频方向子带的均值A^和方差CT^ ;
3) 计算不同尺度各个高频方向子带的斑点噪声方差o^,
其计算公式为
其中,(^是输入SAR图像X的规范化标准差,且对于强度SAR图像有(^二VTT^,对于幅度SAR图像有& =>/(4;r-, / 为已知输入SAR图像的视数,i 2l;
步骤3,禾拥高频方向子带系数0, )的局部均值^朋[(:,,,("》)]和高频方向子带的斑点噪声方差o^,将高频方向子带系数(^,,(",6)区分为信号或是噪声的变换系数,判定准则如下:
<:,,(",6)£信号,如果 mea"[C,,( )]2A:cJ( C,,. (a, 6) e噪声,如果 wea"[C,, (a, <
其中,k在1 5之间取值,縦""[C,,(a,6)]是以高频方向子带系数C, )为中心的5x5大 小的局部区域内系数的均值;
步骤4,对判定后的高频方向子带系数q,(",zo进行如下操作
得到去噪后的高频方向子带系数《,("》);
步骤5,对步骤1中得到的SAR图像X的低频子带Z^(a,6)和去噪后的高频方向子带 系数《,(a,6),作非下采样轮廓波逆变换,得到去噪后的SAR图像I'。
非下采样轮廓波逆变换是利用非下采样轮廓波分解系数对图像进行重构的过程。对 SAR图像,的低频子带和去噪后的高频方向子带系数{《,(",6), 0^/《£ —1, 1&、W, 作非下采样轮廓波逆变换,其过程为
1) 依次对去噪后的高频方向子带(C,:,(",外0S《S丄-1, BS/y作非下采样方向滤波 器组重构,得到第Z, £ — 1, ..., l层非下采样轮廓波分解的带通信号;
2) 对低频子带和第Z层的的带通信号作非下采样塔形滤波器组重构,得到第丄-l层非 下采样轮廓波分解的低通信号;
3) 对第N层非下采样轮廓波分解的低通信号和第N层非下采样轮廓波分解的带通信号 作非下采样塔形滤波器组重构,得到第N-1层非下采样轮廓波分解的低通信号,依次令 N二丄一1, ^一2,…,1;
最终得到精确重构的去噪后SAR图像y ,即第0层非下采样轮廓波分解的低通信号。 参照图2,本发明的基于非下采样轮廓波的合成孔径雷达图像去噪虚拟装置,包括
A. 变换装置,用于输入SAR图像Z,并对其进行丄层非下采样轮廓波变换,得到低
频子带1)£0,6)和不同尺度上的高频方向子带{(^,"6), 0《/^丄—1, 1&、^, q,o,zo表
示SAR图像在尺度2-;上的第;个高频方向子带的系数,£为3~5,然后将不同尺度上各个
高频方向子带分别输入到计算装置和判定装置,同时将低频子带输入到逆变换装置;
B. 计算装置,用于计算不同尺度各个高频方向子带的斑点噪声方差a^ 。该计算装置
包括
二阶累积量计算子装置,用于计算对SAR图像X进行非下采样轮廓波变换时,施加
c,,(",6), 如果 q,(fl,6)s信号 o, 如果 <:,,(0,6)6噪声在不同尺度各个高频方向子带的等价滤波器操作算子的二阶累积量《'(2),
其中,^表示尺度2—'上的第^个高频方向子带,i^']表示第/级塔形滤波的非下采杼 拉普拉斯滤波器系数,f/^表示第y'级方向滤波的非下采样方向滤波器系数,将计算后的 《(2)输入到斑点噪声方差计算子装置;
高频方向子带均值和方差计算子装置,用于计算不同尺度各个高频方向子带的均值 ^和方差CT^ ,将计算后的^^和方差0"^输入到斑点噪声方差计算子装置;
斑点噪声方差计算子装置,用于计算不同尺度各个高频方向子带的斑点噪声方差C^ ;
"1 + GZ2 2
其中,Gz是输入SAR图像Z的规范化标准差,且对于强度SAR图像有C^ = VT7i , 对于幅度SAR图像有Gz =, / 为已知输入SAR图像的视数,A 2 1,将计算 后的o^输入到判定装置;
C.判定装置,用于利用高频方向子带系数<^,.(",6)的局部均值柳""[(:/,0,6)]和高频方向 子带的斑点噪声方差o^ ,将高频方向子带系数(^,,(",6)区分为信号或是噪声的变换系数, 判定准则如下
其中,k在1~5之间取值,附^"0》)]是以高频方向子带系数(^,,0,6)为中心的5x5 大小的局部区域内系数的均值,将判定后的高频方向子带系数C"(",6)输入到去噪装置;
D.去噪装置,用于对判定后的高频方向子带系数<^,,(0》)进行如下去噪操作
得到去噪后的高频方向子带系数c;;,(a,&),将去噪后的高频方向子带系数c;;,(",6)输入到逆
变换装置;
E.逆变换装置,用于对低频子带"£(0,6)和去噪后的高频方向子带系数《,(",6),作非 下采样轮廓波逆变换,得到去噪后的SAR图像I',将去噪后的SAR图像Z'输出。 以下通过仿真实验验证本发明方法的有效性。
C,, 0, 6) e信号,如果 ,(",州2 *
C,.,. (a, 6) s噪声,如果 (a,< A;crCj
C,,(a,6), 如果 C,Xa,6)e信号 0, 如果 q,(a,6)e噪声
11仿真条件使用的原图均为256X256大小的真实SAR图像,图像中包含了多种场景, 如农田、城镇和自然景物,具有丰富的边缘轮廓和细节信息。
仿真内容①选择了传统的SAR图像Lee滤波方法与本发明的方法NSCT进行对比; ②选择了图像处理领域里广泛使用的基于变换域的SAR图像去噪方法,即基于小波变换的 图像去噪方法WT与本发明的方法NSCT进行对比;③选择了不具有平移不变性的轮廓波 变换对SAR图像去噪的方法Contourlet与本发明的方法NSCT进行对比;④选择了非下采 样小波变换NSWT的SAR图像去噪方法与本发明的方法NSCT进行对比。
实验中对图像均采用三层分解。小波变换和NSWT选择DB4小波基函数;Contourlet 变换采用经典的"9-7"塔型分解和"pkva"方向滤波器组;NSCT变换选择"maxflat"非 下采样塔型分解和"dmaxflat7,,非下采样方向滤波器组。由粗尺度到细尺度,Contourlet 变换和NSCT变换的高频方向子带数目分别为2、 4、 8。 Lee滤波采用5x5的窗口。
本发明采用的评价指标为-
(1) SAR图像中通常采用等效视数ENL来评价去噪方法的平滑效果。£AX = //2 /ct2 , 这里^和^分别是SAR图像同质区域的均值和方差,ENL值越大,噪声抑制效果越好。 在实验中,选取SAR原图中方形框标示的同质区域作为测试数据计算等效视数。
(2) 均值比用于衡量去噪图像对原始图像辐射特'性的保持情况,均值比越接近l,表 明辐射特性保持得越好。
(3) 标准差反映的是图像中像素的灰度值与灰度平均之间的分散程度,通过抑制性滤 波,相干斑噪声的强度减弱,图像的方差减小。
(4) 为了考察SAR图像滤波后边缘保持效果的好坏,给出边缘保持指数ESI评价结 果。ESI定义为
,=i去噪图像
z |'=1原始图像
这里,m表示检验象素的个数,DA^和DA^分别是沿着边缘交界处左右或上下相邻 像素的灰度值,ESI越接近1,降斑方法的边缘保持性越好。在表1和表2的边缘保持指 数评价的第一行和第二行,本发明分别给出了水平方向和垂直方向边缘保持情况的指标。
仿真结果
(1)按照所述仿真内容对SAR图像iml仿真的几个图像度量指标如表1。
12表1 SAR图像iml不同去噪方法实验结果比较
原图imlLee滤波WTContourletNSWTNSCT
等效视数10.988597.195673.508796.8548101.6516127.9139
均值比1.00001.06681.05811.06711.0635l扁O
标准差53.725034.684137.579534.311735.342838.8410
边缘保持指数l扁O0.30620,31050.30630.32470.3725
l細O0.27450.23930.30570.31850.3530
(2)按照所述仿真内容对SAR图像im2仿真的几个图像度量指标如表2。表2SAR图像im2不同去噪方法实验结果比较
原图im2Lee滤波WTContourletNSWTNSCT
等效视数8扁884.487956.064596.255678.6345118.7014
均值比l扁O1.20711.19021.21421.1877l細O
标准差38.240624.490928.036225.557625.855327.6203
边缘保持指数l扁O0.28160.24380,32280.33290.3363
l扁O0.29180.25910.35310,34770.3660
从表1和表2对SAR图像iml和im2的实验数据可见,本发明提出的方法NSCT在 等效视数,均值比和边缘保持指数上和其他去噪方法相比均具有优势,如表1中对于SAR 图像iml的实验数据,原图像iml的等效视数为10.9885,采用本发明的方法NSCT得到 的去噪图像的等效视数为127.9139,高于采用Lee滤波,基于小波变换的SAR图像去噪方 法WT,基于轮廓波变换Contourlet,非下采样小波变换NSWT的去噪方法得到的去噪图 像的等效视数,其分别为97.1956, 73.5087, 96.8548, 101.6516;本发明对应的均值比结 果为1.0000,说明去噪结果与原图辐射特性保持一致,而且本发明方法得到的边缘保持指 数的水平方向和垂直方向的指标分别为0.3725和0.3530,都高出其他方法。五种方法结果 的标准差都比原图有所下降,表示这些方法都达到了滤波的作用。
(3) 用本发明和已有方法对SAR图像iml的去噪结果如图3。其中图3 (a)为SAR 图像iml;图3(b)为采用传统Lee滤波的去噪图像;图3 (c)为基于小波变换的SAR 图像去噪方法WT得到的去噪图像;图3 (d)为采用轮廓波变换的图像去噪方法Contourlet 得到的去噪图像;图3 (e)为采用非下采样小波变换的SAR图像去噪方法NSWT得到的 去噪图像;图3 (f)为采用本发明的方法NSCT得到的去噪图像。
(4) 用本发明和已有方法对SAR图像im2的去噪结果如图4。其中图4 (a)为SAR 图像im2;图4 (b)为采用传统Lee滤波的去噪图像;图4 (c)为基于小波变换的SAR 图像去噪方法WT得到的去噪图像;图4 (d)为采用轮廓波变换的图像去噪方法Contourlet 得到的去噪图像;图4 (e)为采用非下采样小波变换的SAR图像去噪方法NSWT得到的去噪图像;图4 (f)为采用本发明的方法NSCT得到的去噪图像。
参照图3和图4可见,Lee滤波可以对斑点噪声进行很好的平滑,但是模糊了区域的 边界,丢失了一些细节信息。基于小波的方法和基于Contourlet的方法由于采用的是非平 移不变的信号表示,所以易在边缘附近产生抖动,使得一些细小的边缘无法得到保持。基 于NSWT的去噪结果整体清晰度较高。相比于以上几种算法,本发明的基于非下采样轮廓 波的SAR图像去斑方法NSCT在提高了同质区域平滑性的同时,很好地保持了边缘、点 目标和线目标,增强了图像的可解释性。
本发明相比于现有的SAR图像去噪方法,不管从客观参数的评价上,还是从图像的视 觉质量上来看都具有优越性,可以有效避免某些变换因缺乏平移不变性而产生的图像失 真,在相干斑抑制的同时,可以很好地保持SAR图像的细节信息,是一种可行有效的SAR 图像去噪方法。
权利要求
1. 一种基于非下采样轮廓波的合成孔径雷达图像去噪方法,包括如下步骤(1)输入SAR图像X,并对其进行L层非下采样轮廓波变换,得到低频子带DL(a,b)和不同尺度上的高频方向子带{Cl,i(a,b),0≤l≤L-1,1≤i≤kl},Cl,i(a,b)表示SAR图像在尺度2-l上的第i个高频方向子带的系数,L为3~5;(2)计算不同尺度各个高频方向子带的斑点噪声方差;(3)利用高频方向子带系数Cl,i(a,b)的局部均值mean[Cl,i(a,b)]和高频方向子带的斑点噪声方差,将高频方向子带系数Cl,i(a,b)区分为信号或是噪声的变换系数,判定准则如下其中,k在1~5之间取值,mean[Cl,i(a,b)]是以高频方向子带系数Cl,i(a,b)为中心的5×5大小的局部区域内系数的均值;(4)对判定后的高频方向子带系数Cl,i(a,b)进行如下去噪操作得到去噪后的高频方向子带系数(5)对步骤(1)中得到的低频子带DL(a,b)和去噪后的高频方向子带系数作非下采样轮廓波逆变换,得到去噪后的SAR图像X*。
2.根据权利1所述的非下采样轮廓波的合成孔径雷达图像去噪方法,其特征在于步骤 (2)所述"计算不同尺度各个高频方向子带的斑点噪声方差OV/',按如下步骤进行-2a)计算对SAR图像I进行非下采样轮廓波变换时,施加在不同尺度各个高频方向 子带的等价滤波器操作算子的二阶累积量《'(2),其中,A表示尺度2"上的第&个高频方向子带,i/^表示第/级塔形滤波的非下采样拉普拉斯滤波器系数,c/ 表示第y级方向滤波的非下采样方向滤波器系数;2b)计算不同尺度各个高频方向子带的均值^x和方差cr^ ; 2c)计算不同尺度各个高频方向子带的斑点噪声方差0^,'6>/,,0,6)£信号,如果 mea"[C,,,0,6)]2^r( C,, (a, 6) e噪声,如果 we朋[C,, (a, 6)] < A:cj(C"0,6), 如果 C,.,(fl,6)e信号0, 如果 <:,,(0,6;^噪声《;(2^n(zw])2)n(i;(" 〗)2)其中,A是输入SAR图像Z的规范化标准差,且对于强度SAR图像有Gz二VTT^, 对于幅度SAR图像有Gz =V(4" —iV丑,i 为己知输入SAR图像的视数,/ 21。
3. —种基于非下采样轮廓波的合成孔径雷达图像去噪虚拟装置,包括变换装置,用于输入SAR图像X,并对其进行Z层非下采样轮廓波变换,得到低频 子带/\(",6)和不同尺度上的高频方向子带{(^,0,6), 0^/《丄—1, l&、w, q,o》)表示 SAR图像在尺度2"上的第/个高频方向子带的系数,丄为3 5;计算装置,用于计算不同尺度各个高频方向子带的斑点噪声方差cr^ ;判定装置,用于利用高频方向子带系数C,^,6)的局部均值,a"[q,(",6)]和高频方向 子带的斑点噪声方差o^ ,将高频方向子带系数<^,,(",6)区分为信号或是噪声的变换系数, 判定准则如下其中,k在1 5之间取值,柳朋0,州是以高频方向子带系数(:,,0》)为中心的5x5 大小的局部区域内系数的均值;去噪装置,用于对判定后的高频方向子带系数<^,,(",6)进行如下去噪操作得到去噪后的高频方向子带系数C/,(fl,6);逆变换装置,用于对低频子带D,(a》)和去噪后的高频方向子带系数《,(a》),作非下 采样轮廓波逆变换,得到去噪后的SAR图像义'。
4.根据权利3所述的非下采样轮廓波的合成孔径雷达图像去噪虚拟装置,其特征在于计算装置,进一步包括二阶累积量计算子装置,用于计算对SAR图像Z进行非下采样轮廓波变换时,施加 在不同尺度各个高频方向子带的等价滤波器操作算子的二阶累积量《'(2),其中,A表示尺度2-'上的第^个高频方向子带,//^表示第/级塔形滤波的非下采样 拉普拉斯滤波器系数,《;]表示第7级方向滤波的非下采样方向滤波器系数;C,, 0,6) e噪声,如果 weaw[C,,0,6)] < A;cr(C,,0,6), 如果 q,(a,6)e信号 0, 如果 (7,,(",6)£噪声高频方向子带均值和方差计算子装置,用于计算不同尺度各个高频方向子带的均值斑点噪声方差计算子装置,用于计算不同尺度各个高频方向子带的斑点噪声方差cr^,其中,C^是输入SAR图像Z的规范化标准差,且对于强度SAR图像有(^^VT7i, 对于幅度SAR图像有Gz二^/(4;r-1)/i , i 为已知输入SAR图像的视数,i 2 1 。
全文摘要
本发明公开了一种基于非下采样轮廓波的合成孔径雷达图像去噪方法,主要解决现有合成孔径雷达图像去噪方法难以有效保持图像细节的问题,其步骤是(1)输入SAR图像X,并对其进行L层非下采样轮廓波变换;(2)计算不同尺度各个高频方向子带的斑点噪声方差σ<sub>C#-[B</sub>];(3)利用高频方向子带系数C<sub>l,i</sub>(a,b)的局部均值mean[C<sub>l,i</sub>(a,b)],将其区分为信号或是噪声的变换系数;(4)对判定后的高频方向子带系数C<sub>l,i</sub>(a,b)中的信号部分进行保留,得到去噪后的高频方向子带系数C<sub>l,i</sub><sup>*</sup>(a,b);(5)对低频子带和去噪后的高频方向子带系数C<sub>l,i</sub><sup>*</sup>(a,b),作非下采样轮廓波逆变换,得到去噪后的SAR图像X<sup>*</sup>。本发明在有效去除相干斑噪声同时能够有效保持图像的细节,去噪图像无抖动失真,可用于合成孔径雷达图像的预处理阶段。
文档编号G01S13/90GK101482617SQ20091002095
公开日2009年7月15日 申请日期2009年1月16日 优先权日2009年1月16日
发明者彪 侯, 公茂果, 芳 刘, 霞 常, 杨淑媛, 焦李成, 爽 王 申请人:西安电子科技大学