发明名称 基于非下采样轮廓波的合成孔径雷达图像去噪方法
摘要 本发明公开了一种基于非下采样轮廓波的合成孔径雷达图像去噪方法,主要解决现有合成孔径雷达图像去噪方法难以有效保持图像细节的问题,其步骤是:(1)输入SAR图像X,并对其进行L层非下采样轮廓波变换;(2)计算不同尺度各个高频方向子带的斑点噪声方差<img file="200910020952.9_ab_0.GIF" wi="44" he="26" />(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>。本发明在有效去除相干斑噪声同时能够有效保持图像的细节,去噪图像无抖动失真,可用于合成孔径雷达图像的预处理阶段。
申请公布号 CN101482617B 申请公布日期 2011.06.01
申请号 CN200910020952.9 申请日期 2009.01.16
申请人 西安电子科技大学 发明人 焦李成;常霞;王爽;侯彪;刘芳;杨淑媛;公茂果
分类号 G01S13/90(2006.01)I;G06T5/00(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于非下采样轮廓波的合成孔径雷达图像去噪方法,包括如下步骤:(1)输入SAR图像X,并对其进行L层非下采样轮廓波变换,得到低频子带D<sub>L</sub>(a,b)和不同尺度上的高频方向子带{C<sub>l,i</sub>(a,b),0≤l≤L-1,1≤i ≤k<sub>l</sub>},k<sub>l</sub>表示图像在尺度2<sup>-l</sup>上的高频方向子带的数目,C<sub>l,i</sub>(a,b)表示SAR图像在尺度2<sup>-l</sup>上的第i个高频方向子带的系数,L为3~5;(2)计算不同尺度各个高频方向子带的斑点噪声方差<img file="FSB00000465560600011.GIF" wi="105" he="51" />2a)计算对SAR图像X进行非下采样轮廓波变换时,施加在不同尺度各个高频方向子带的等价滤波器操作算子的二阶累积量<img file="FSB00000465560600012.GIF" wi="111" he="68" /><img file="FSB00000465560600013.GIF" wi="1066" he="135" />其中,<img file="FSB00000465560600014.GIF" wi="74" he="58" />表示第i级塔形滤波的非下采样拉普拉斯滤波器系数,<img file="FSB00000465560600015.GIF" wi="75" he="79" />表示第j级方向滤波的非下采样方向滤波器系数;2b)计算不同尺度各个高频方向子带的均值μ<sub>X</sub>和方差<img file="FSB00000465560600016.GIF" wi="168" he="65" />2c)计算不同尺度各个高频方向子带的斑点噪声方差<img file="FSB00000465560600017.GIF" wi="105" he="52" /><img file="FSB00000465560600018.GIF" wi="555" he="154" />其中,G<sub>Z</sub>是输入SAR图像X的规范化标准差,且对于强度SAR图像有<img file="FSB00000465560600019.GIF" wi="269" he="70" />对于幅度SAR图像有<img file="FSB000004655606000110.GIF" wi="420" he="86" />R为已知输入SAR图像的视数,R≥1;(3)利用高频方向子带系数C<sub>l,i</sub>(a,b)的局部均值mean[C<sub>l,i</sub>(a,b)]和高频方向子带的斑点噪声方差<img file="FSB000004655606000111.GIF" wi="105" he="59" />将高频方向子带系数C<sub>l,i</sub>(a,b)区分为信号或是噪声的变换系数,判定准则如下:<img file="FSB000004655606000112.GIF" wi="1308" he="197" />其中,k在1~5之间取值,mean[C<sub>l,i</sub>(a,b)]是以高频方向子带系数C<sub>l,i</sub>(a,b)为中心的5×5大小的局部区域内系数的均值;(4)对判定后的高频方向子带系数C<sub>l,i</sub>(a,b)进行如下去噪操作:<img file="FSB000004655606000113.GIF" wi="1129" he="205" />得到去噪后的高频方向子带系数<img file="FSB00000465560600021.GIF" wi="191" he="76" />(5)对步骤(1)中得到的低频子带D<sub>L</sub>(a,b)和去噪后的高频方向子带系数<img file="FSB00000465560600022.GIF" wi="190" he="63" />作非下采样轮廓波逆变换,得到去噪后的SAR图像X<sup>*</sup>。
地址 710071 陕西省西安市太白路2号