CN101566688B - 基于邻域方向性信息的sar图像降斑方法 - Google Patents
基于邻域方向性信息的sar图像降斑方法 Download PDFInfo
- Publication number
- CN101566688B CN101566688B CN2009100228661A CN200910022866A CN101566688B CN 101566688 B CN101566688 B CN 101566688B CN 2009100228661 A CN2009100228661 A CN 2009100228661A CN 200910022866 A CN200910022866 A CN 200910022866A CN 101566688 B CN101566688 B CN 101566688B
- Authority
- CN
- China
- Prior art keywords
- coefficient
- mask
- image
- sar image
- dissociation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种基于邻域方向性信息的SAR图像降斑方法,主要解决现有SAR图像降斑方法不能很好的抑制斑点噪声和丢失部分边缘和细节信息的缺点。其过程为:对原图像做非下采样Contourlet分解,对每个高频子带系数估计二进掩码;然后根据所对应的不同的掩码用修正对数高斯分布和混合指数分布获得条件似然比的表达式;再根据邻域方向性模型获得的非下采样Contourlet变换的方向信息求先验比;最后根据似然比和先验比获得缩减因子对每个高频子带系数进行缩减,将缩减变化以后的系数再经过重构得到降斑后的输出图像。本方法具有有效去除SAR图像斑点噪声和边缘保持性好的优点,可用于对边缘信息丰富、细节丰富的SAR图像,特别是含有机场、跑道、道路的SAR图像降斑。
Description
技术领域
本发明属于图像处理技术领域,涉及图像噪声抑制,具体地说是一种基于非下采样contourlet域的SAR图像降斑方法,可用于合成孔径雷达图像的斑点噪声的抑制。
背景技术
合成孔径雷达(SAR)是一种高分辨率成像雷达。它具有全天候、多极化、多视角、多俯角数据获取能力及对一些地物的穿透能力,不仅广泛地被应用在军事上,在农业、气象、地形地貌、灾情监控等民用上也有大量的应用。但由于SAR发射的是相干波,这些相干波经过与地物的相干作用,特别是地物的后向散射作用,使目标回波信号产生了衰减,这种衰减表现在图像上就是相干斑点噪声。因此如何抑制图像中的相干斑噪声,提高图像的解读能力以及获得更多的信息成为一个重要的问题。
SAR图像降斑的首要目标是在滤除斑点噪声的同时,尽可能的保持图像的细节信息。斑点噪声是一种复杂的乘性噪声模型的信号。对于斑点噪声这种特殊的性质,在过去的近二十年中,人们提出了很多经典的SAR降斑方法,如Lee滤波,增强Lee滤波,Kuan滤波等等。这些方法是用一个已定义好的滤波器窗来估计局域斑点噪声的方差,并进行滤波处理,其结果通常过分的平滑了边缘细节信息,这些方法均在一定程度上收到了较好的效果。1995年,美国学者Donoho把小波理论引入到图像去噪中,提出了小波软阈值方法。小波软阈值方法是一种非线性的算法,依然存在破坏图像细节信息的问题,对图像的辐射特性保持也不好。
后来美国学者Do and Vetterli提出了一种图像的多方向多尺度表示方法,称作Contourlet变换。Contourlet变换利用金字塔滤波器获得多尺度分解,利用方向滤波器组获得方向的合成。由于使用了几何信息,Contourlet变换在例如降噪和纹理恢复中的图像分析应用中,得到了比离散小波变换更好的效果。Contourlet变换经过下采样和上采样,并不是平移不变的,而平移不变性在图像分析应用例如边缘检测、轮廓提取和图像去噪中是至关重要的。
最近美国学者Cunha等人提出了非下采样分解机制的Contourlet变换,这种变换具有完全的平移不变性,在图像去噪的应用中已初步显示出其明显的优势:与UWT相比,非下采样Contourlet变换能够提供更好的图像恢复结果,在细节结构的保持方面更优于前者;同样,非下采样Contourlet变换的平移不变性也解决了基于Contourlet变换的去噪处理中出现的块状效应问题。但是这种基于非下采样Contourlet变换域降斑方法并没有充分考虑图像的几何特征以及SAR图像在非下采样Contourlet变换域的统计特性和方向性特征,因而降斑后的SAR图像平滑区域的斑点噪声滤除不充分,同时图像的细节和边缘信息部分丢失。
发明内容
本发明的目的在于克服上述已有技术的不足,提出一种基于邻域方向性信息的SAR图像降斑方法,以充分滤除降斑后的SAR图像平滑区域的斑点噪声,完整的保留图像的细节和边缘信息。
实现本发明目的的技术思路是将修正的对数高斯分布和混合指数分布引入到缩减因子的表示中,并利用方向性邻域模型结合非下采样Contourlet变换的高度冗余性和多方向选择性对非下采样Contourlet变换分解高频系数进行自适应萎缩,得到滤除斑点噪声、保留细节信息的SAR图像。其具体实现步骤包括如下:
(1)对原始SAR图像y进行非下采样Contourlet分解,得到在尺度j上位置l处的分解系数nl,j;
(3)找出第j层系数掩码为1和为0的系数,分别用修正对数高斯分布和混合指数分布逼近,得到概率密度分布的条件似然比ξl;
(4)根据设定的方向性邻域模型计算分解系数的缩减因子先验比ηl,得到非下采样Contourlet分解系数缩减因子ql;
(5)应用缩减因子对非下采样Contourlet分解的系数nl,j进行缩减;
本发明与现有技术相比具有如下优点:
1)本发明由于采用改进对数高斯分布和混合指数分布对掩码为1和掩码为0的非下采样Contourlet分解系数的直方图进行逼近,更加准确的捕捉了SAR图像在非下采样Contourlet变换域的统计特性,充分滤除了SAR图像平滑区域的斑点噪声。
2)本发明由于将方向性邻域模型引入缩减因子的计算中,充分考虑到了SAR图像本身的方向特性,保留了更加丰富的图像边缘和细节信息。
3)仿真结果表明,本发明方法较其它几种现有的经典的SAR图像降斑方法,在平滑区域的平滑效果以及边缘保持能力方面都有显著的提高。
附图说明
图1是本发明的流程框图;
图2是本发明使用的邻域方向性模型示意图;
图3是本发明求得的SAR图像高频子带系数为1的掩码以及概率分布拟合示意图;
图4是本发明与现有两种降斑方法应用于X波段幅度SAR图像的降斑仿真结果对比图;
图5是本发明与现有两种降斑方法应用于Ku波段强度SAR图像的降斑仿真结果对比图;
图6是对图4和图5的局部放大图。
具体实施方式
参照图1,本发明的具体实现步骤如下:
步骤1,输入SAR图像记为y,该SAR图像本身就是被斑点噪声污染的图像,因此不需要像研究自然图像去噪一样,给原图加上一个随机噪声或者某种特性的噪声,可直接对该图像进行降斑处理。
步骤2,将图像y经过非下采样Contourlet分解,得到一个与原图像大小相同的低频图像和J个尺度,每个尺度有D个高频图像,高频图像在尺度j上位置l处的系数记为nl,j。
(2.1)将输入图像经过非下采样金字塔分解,分别得到与原图大小相同的一个低频图像和一个高频图像;
(2.2)将高频图像输入非下采样方向滤波器组,得到D个与原图像大小相等D个方向的高频图像;
(2.3)将在步骤(2.1)中输出的低频图像再重复以上步骤(2.1)和(2.2)J次,最终得到非下采样Contourlet分解的一个低频图像和J个尺度,每个尺度有D个与原图像大小相同的高频图像,该高频图像在尺度j上位置l处的系数记为nl,j,j∈1…J。
步骤3,估计尺度j上位置l处的系数nl,j二值掩码
其中,σj为噪声标准差,其表达式为:σj=Median(|nl,j-Median(nl,j)|)/0.6745,为系数nl,j的掩码:表明nl,j为重要系数,主要对应SAR图像中的细节和有用信息;表示nl,j为不重要系数,主要对应SAR图像中的斑点噪声信息。
步骤4,根据系数二值掩码计算概率密度分布的条件似然比。
概率密度分布的条件似然比ξl由掩码为1的系数概率分布函数p(ml|1)和掩码为0的系数概率分布函数p(ml|0)两部分组成。
第一部分掩码为1的系数的概率分布函数p(ml|1),采用改进对数高斯概率密度函数进行拟合:
其中,x为非下采样Contourlet分解的高频系数值,ml为nl,j的幅值,即ml=|nl,j|,μ是ln(ml)的均值,σ1是ln(ml))的标准差。
图3(a),(b)所示为图4(a)中SAR图像进行非下采样Contourlet分解的第三个方向粗尺度与次粗尺度的掩码估计结果,图3(c),(d)为掩码为1的概率分布拟合图,其中直线所示为掩码为1的系数的直方图,带星号的直线所示为改进对数高斯概率密度函数的拟合曲线,带叉号的直线为统计先验指导的降斑方法中应用的Gamma概率密度函数的拟合曲线,从图3中可以看出,本发明中的改进高斯概率密度函数比Gamma概率密度函数更好的逼近了掩码为1的系数直方图。
第二部分掩码为0的系数的概率分布函数p(ml|0),采用混合指数概率密度函数拟合:
其中,x为非下采样Contourlet分解的高频系数值,λ1为第一个指数分布尺度参数,λ2为第二个指数分布的尺度参数,w1和w2分别为这两个指数分布的权重参数,λ1,λ2,w1,w2这些参数的值均采用期望最大化算法估计;
这样,概率密度分布的条件似然比表示为掩码为1的系数概率分布函数与掩码为0的系数概率分布函数之比:
ξl=p(ml|1)/p(ml|0)。 (4)
步骤5,根据方向性邻域模型计算系数缩减因子先验比。
方向性邻域模型如图2所示,是以位置l为中心的各向同性邻域模型如图2(a)中所示,oD(l)是以位置l为中心的方向性邻域模型,图2(b)所示为4个方向的邻域模型o4(l),图2(c)所示为8个方向的邻域模型o8(l)。为方向性邻域模型内的系数掩码值,如果模型中的系数掩码值全部为0,则系数缩减因子先验比为ηl=1;否则根据各向异性邻域模型oD(l)判断,如果模型oD(l)中任意一个方向的系数掩码全部为1,则系数缩减因子先验比为ηl=exp{γ·(2·D-1)},γ为控制因子;若以上两个条件都不满足,则系数缩减因子先验比为
综上,系数缩减因子先验比表示为:
首先,根据步骤4和5中得到概率密度分布的似然比ξl和系数缩减因子先验比ηl,计算非下采样Contourlet分解系数的缩减因子ql为:
ql=ξlηl/(1+ξlηl), (6)
然后,通过系数缩减因子ql与非下采样Contourlet分解系数nl,j的点积得到缩减后的非下采样Contourlet分解系数为:
(7.1)将缩减后当前尺度j的D个方向的高频系数输入非下采样重构方向滤波器组,得到一个与原SAR图像大小相等的高频图像;
(7.2)将步骤(7.1)中得到的高频图像与步骤(2.3)中保留的低频图像一起输入到非下采样重构金字塔滤波器,得到重构的尺度j的低频图像;
本发明的效果可以通过以下仿真结果进一步说明。
1.仿真内容
应用本发明降斑方法、统计先验指导的降斑方法与Gamma-MAP降斑方法,分别对一幅X波段幅度SAR图像和一幅Ku波段强度SAR图像进行降斑,并给出了等效视数ENL,图像均值M和标准差V作为客观评价标准来说明本发明应用于SAR图像降斑性能。ENL越高,说明平滑区域的降斑性能越好,降斑后的均值越接近原图像均值越好,降斑后标准差越低,说明平滑效果越好。仿真结果中还给出了本发明方法与统计先验指导的降斑方法应用于SAR图像的降斑结果细节放大图,对比结果进一步说明了本发明在细节以及边缘保持方面的优越性。
仿真内容1,将现有的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)中所示矩形区域1、2为表1中计算ENL所需的同质区域,表1列出了仿真内容(1)的仿真结果所得的均值、方差和等效视数对比结果。
表1:不同降斑方法客观指标评价:等效视数ENL、均值M、标准差V
从表1中可以看出,本发明获得了比现有其他方法更高的等效视数,降斑后的图像均值非常接近原图像均值,降斑后的标准差最低。因此本发明取得了比其他降斑方法更优的平滑降斑效果。
图5是仿真内容2的仿真结果。其中图5(a)为原SAR图像,图5(b)为经典Gamma-MAP方法降斑结果,图5(c)为统计先验指导的方法降斑图,图5(d)为本发明降斑结果图。图5(a)中所示矩形区域3、4为表2中计算ENL所需的同质区域。表2列出了仿真内容(2)的仿真结果所得的均值、方差和等效视数对比结果。
表2:不同降斑方法客观指标评价:等效视数ENL、均值M、标准差V
从表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)的图像局部放大图。从细节放大对比图可以看出,本发明得到的降斑图像能更好的保持原图像的细节和边缘信息,获得了非常好的直观效果。
Claims (3)
1.一种基于邻域方向性信息的SAR图像降斑方法,包括如下步骤:
(1)对原始SAR图像y进行非下采样Contourlet分解,得到在尺度j上位置 l处的分解系数n l,j;
(3)找出尺度j的系数二值掩码为1和为0的系数,分别用修正对数高斯分布和混合指数分布逼近,得到概率密度分布的条件似然比ξ l ,该步骤具体实现如下:
(3b)找出尺度j的系数二值掩码为0所对应的系数的绝对值 掩码为0的系数的直方图用如下条件概率分布表示:
其中,λ1为第一个指数分布尺度参数,λ2为第二个指数分布的尺度参数,w1和w2分别为这两个指数分布的权重参数,这些参数λ1,λ2,w1,w2的值采用期望最大化算法估计;
(3c)根据掩码为1的系数的直方图的条件概率分布函数和掩码为0的系数的直方图的条件概率分布函数,计算概率密度分布的条件似然比为:
ξ l =p(m l |1)/p(m l |0);
(4)根据设定的方向性邻域模型计算分解系数的缩减因子先验比η l ,得到非下采样Contourlet分解系数缩减因子q l ,该步骤具体实现如下:
(4b)根据方向性邻域模型计算系数缩减因子先验比为:
(4c)根据概率密度分布的条件似然比ξ l 和系数缩减因子先验比η l ,计算系数缩减因子为:
q l =ξ l η l /(1+ξ l η l );
(5)应用缩减因子对非下采样Contourlet分解的系数n l,j进行缩减;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100228661A CN101566688B (zh) | 2009-06-05 | 2009-06-05 | 基于邻域方向性信息的sar图像降斑方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100228661A CN101566688B (zh) | 2009-06-05 | 2009-06-05 | 基于邻域方向性信息的sar图像降斑方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101566688A CN101566688A (zh) | 2009-10-28 |
CN101566688B true CN101566688B (zh) | 2012-02-08 |
Family
ID=41282942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100228661A Expired - Fee Related CN101566688B (zh) | 2009-06-05 | 2009-06-05 | 基于邻域方向性信息的sar图像降斑方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101566688B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105574832A (zh) * | 2015-12-11 | 2016-05-11 | 上海交通大学 | 基于迭代方向滤波器组可逆深度卷积网络结构 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887580A (zh) * | 2010-07-23 | 2010-11-17 | 扬州万方电子技术有限责任公司 | 非下采样轮廓波变换域的图像降噪方法 |
CN102073998B (zh) * | 2011-01-20 | 2012-07-04 | 山东大学 | 基于非下采样轮廓波的边缘检测方法 |
CN102129672B (zh) * | 2011-03-15 | 2013-06-12 | 西安电子科技大学 | 基于swt域改进粒子滤波的sar图像降斑方法 |
CN102368332B (zh) * | 2011-11-17 | 2013-04-17 | 西安电子科技大学 | 基于nsct域局部高斯模型的sar图像降斑方法 |
CN102903080B (zh) * | 2012-09-06 | 2015-10-28 | 西安工程大学 | 合成孔径雷达图像相干斑噪声抑制性能的非监督评估方法 |
CN107392861B (zh) * | 2017-06-29 | 2020-05-29 | 南京航空航天大学 | 一种基于高斯比例混合模型的稀疏表示sar图像降斑方法 |
CN108389215B (zh) * | 2018-03-15 | 2021-04-23 | 京东方科技集团股份有限公司 | 一种边缘检测方法、装置、计算机存储介质及终端 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303764A (zh) * | 2008-05-16 | 2008-11-12 | 西安电子科技大学 | 基于非下采样轮廓波的多传感器图像自适应融合方法 |
-
2009
- 2009-06-05 CN CN2009100228661A patent/CN101566688B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303764A (zh) * | 2008-05-16 | 2008-11-12 | 西安电子科技大学 | 基于非下采样轮廓波的多传感器图像自适应融合方法 |
Non-Patent Citations (4)
Title |
---|
Wilfried Philips等.DESPPECKLING SAR IMAGES USING WAVELETS AND A NEW CLASS OF ADAPTIVE SHRINKAGE ESTIMATORS.《Proc ICIP》.2001,第2卷233-236. * |
孙强等.统计先验指导的非下采样Contourlet变换域SAR图像降斑.《西安电子科技大学学报(自然科学版)》.2008,第35卷(第1期),14-21. * |
贾建等.利用方向特性实现非下采样Contourlet变换阈值去噪.《西安电子科技大学学报(自然科学版)》.2009,第36卷(第2期),269-273. * |
贾建等.基于概率模型的非下采样Contourlet变换图像去噪.《西北大学学报(自然科学版)》.2009,第39卷(第1期),13-18. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105574832A (zh) * | 2015-12-11 | 2016-05-11 | 上海交通大学 | 基于迭代方向滤波器组可逆深度卷积网络结构 |
CN105574832B (zh) * | 2015-12-11 | 2019-11-15 | 上海交通大学 | 基于迭代方向滤波器组可逆深度卷积网络结构 |
Also Published As
Publication number | Publication date |
---|---|
CN101566688A (zh) | 2009-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101566688B (zh) | 基于邻域方向性信息的sar图像降斑方法 | |
CN101482617B (zh) | 基于非下采样轮廓波的合成孔径雷达图像去噪方法 | |
CN101980286B (zh) | 结合双树复小波与双变量模型的sar图像降斑方法 | |
Xu et al. | Patch ordering-based SAR image despeckling via transform-domain filtering | |
CN103077508B (zh) | 基于变换域非局部和最小均方误差的sar图像去噪方法 | |
CN104809734A (zh) | 一种基于引导滤波的红外图像与可见光图像融合的方法 | |
Min et al. | Image denoising via bivariate shrinkage function based on a new structure of dual contourlet transform | |
CN103295204B (zh) | 一种基于非下采样轮廓波变换的图像自适应增强方法 | |
CN101483777B (zh) | 基于自适应多尺度Bandelet包的SAR图像去噪压缩方法 | |
CN101719268B (zh) | 基于改进Directionlet域的广义高斯模型图像去噪方法 | |
CN103020922A (zh) | 基于pca变换的sar图像相干斑抑制方法 | |
CN102509263A (zh) | 基于sar图像局部统计特性的k-svd相干斑抑制方法 | |
Gao et al. | Directionlet-based denoising of SAR images using a Cauchy model | |
Lu et al. | SAR speckle reduction using Laplace mixture model and spatial mutual information in the directionlet domain | |
CN103093432B (zh) | 基于极化分解和图像块相似性的极化sar图像降斑方法 | |
Wang et al. | An image NSCT-HMT model based on copula entropy multivariate Gaussian scale mixtures | |
CN103426145A (zh) | 一种基于多分辨分析的合成孔径声纳相干斑噪声抑制方法 | |
Gupta et al. | Homomorphic wavelet thresholding technique for denoising medical ultrasound images | |
CN102693530B (zh) | 基于目标提取和srad算法的sar图像去斑方法 | |
CN101984461A (zh) | 基于可操纵金字塔的统计模型图像去噪方法 | |
CN102722878A (zh) | 基于目标提取和ppb算法的sar图像去斑方法 | |
Wang et al. | An edge-preserving adaptive image denoising | |
CN103345739B (zh) | 一种基于纹理的高分辨率遥感影像建筑区指数计算方法 | |
Zhang et al. | Bayesian-based speckle suppression for SAR image using contourlet transform | |
Lu et al. | Directionlet-based method using the Gaussian mixture prior to SAR image despeckling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120208 Termination date: 20170605 |