CN100417191C - 一种图像组合降噪方法 - Google Patents
一种图像组合降噪方法 Download PDFInfo
- Publication number
- CN100417191C CN100417191C CNB2006100307462A CN200610030746A CN100417191C CN 100417191 C CN100417191 C CN 100417191C CN B2006100307462 A CNB2006100307462 A CN B2006100307462A CN 200610030746 A CN200610030746 A CN 200610030746A CN 100417191 C CN100417191 C CN 100417191C
- Authority
- CN
- China
- Prior art keywords
- noise reduction
- image
- contourlet
- centerdot
- noise
- 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
Abstract
本发明涉及一种图像组合降噪方法。本方法先通过Contourlet变换对输入的带噪图像进行多尺度、多方向的稀疏分解,然后根据Contourlet变换域系数服从广义高斯分布,在Contourlet域进行Bayes收缩阈值法降噪,并通过Contourlet反变换得到预降噪图像,最后,采用Wiener滤波法对预降噪图像进行进一步降噪处理,得到最终的降噪图像,达到图像降噪目的。采用本方法除去图像中的噪声,可提高图像质量,本发明可广泛应用于军事领域或非军事领域的光学成像、目标检测和安全监控等系统中。
Description
技术领域
本发明涉及一种图像组合降噪方法,该方法采用基于移不变Contourlet(轮廓小波)变换域降噪和维纳(Wiener)滤波的图像组合降噪方法,除去图像中的噪声,以提高图像质量。在军事领域和非军事领域如光学成像、目标检测、安全监控等系统中均有广泛应用。
背景技术
通常,我们得到的图像都受到不同程度的噪声污染,为了后续的进一步处理,很有必要进行降噪处理,滤出噪声,并尽可能地保留图像的所有特征信息,以提高图像的恢复质量。目前,图像降噪方法主要有:空域滤波,如均值滤波和中值滤波等,以及变换域滤波,如低通滤波等。
近二十来,小波变换以其良好的时频特性和多分辨率思想,在信号和图像处理领域得到了广泛的应用,并成功地应用到了图像降噪领域。但是,由一维小波通过张量积形成的二维可分离小波变换只能有效地表示一维奇异信息即点奇异信息,而不能有效地描述图像中的二维或高维奇异信息,如线、轮廓等重要信息,从而制约了小波降噪方法的性能。
Contourlet变换作为一种新的信号分析工具,解决了小波变换不能有效表示二维或更高维奇异性的缺点,能准确地将图像中的边缘捕获到不同尺度、不同频率、不同方向的子带中。它不仅具有小波变换的多尺度特性,还具有小波变换不具有的方向性和各向异性,能更稀疏的表示图像中的轮廓特征,因此能很好地应用于图像处理中,包括图像降噪方面。由于,Contourlet变换缺乏平移不变性,在降噪过程中会产生伪吉布斯(Gibbs)现象,循环平移(Cycle Spinning)方法可以用来抑制Contourlet域阈值降噪方法产生的伪Gibbs现象,降噪性能明显优于小波阈值降噪方法,在很大程度上提高了降噪方法的性能。但是,实际情况表明,这种方法还不能完全去除噪声,降噪图像中仍然残留少量噪声,影响图像质量,需要采取进一步的降噪处理。
发明内容
本发明的目的在于针对现有图像降噪方法存在的不足,提出了一种图像组合降噪方法,用于去除图像中的噪声,以提高图像质量。
为了达到上述目的,本发明采用下述技术方案:
一种图像组合降噪方法。其特征在于采用基于移不变Contourlet变换域降噪和Wiener滤波的图像组合降噪方法,该方法是先在Contourlet域,根据系数服从广义高斯分布,选取贝叶斯(Bayes)收缩阈值,提高图像降噪的效果。然后,再对经过Contourlet域降噪后得到的预降噪图像,采用Wiener滤波方法进行进一步的降噪处理,达到图像降噪的目的。
假设观测到的噪声图像为
I=f+n (1)
其中f为原始图像,n为独立同分布的高斯白噪声信号N(0,σ2).
上述降噪方法的具体步骤如下:
1.始化设置。令i=0,j=0,设定行方向和列方向上的最大平移量N1和N2.同时设定Contourlet变换中的拉普拉斯塔式分解层数K和每层中的方向分解数Lk;
2.对输入的带噪图像I在行和列方向上进行有效位移量的循环平移,得到平移图像
Sij=Cij(I), (2)
其中i∈(0,N1)和j∈(0,N2)分别为行方向和列方向上的平移量;
3.对得到的平移图像Sij进行多尺度、多方向的Contourlet稀疏分解,即
其中T(·)为Contourlet变换。从而得到一幅低频子图像Sif和一系列具有不同分辨率的高频子图像Shf (kl),其中k∈(1,K)和l∈(1,Lk)标明子图像位于第k层拉普拉斯塔式分解的第l方向;
4.对Contourlet变换后的高频子图像Shf (k,l)进行阈值降噪处理,得到降噪子图像,
其中,Λ(·)为阈值函数,TB为阈值参数。根据Contourlet域系数服从广义高斯分布,选用Bayes阈值
5.对第4步中得到的所有降噪高频子图像SDhf (k,l)和第3步中得到的低频子图像Sif实施Contourlet逆变换,得到在行方向和列方向上分别平移i和j后的降噪图像,
其中,T-1(·)为Contourlet逆变换;
6.对第5步中得到的图像Sij nf进行相应平移量的逆向循环平移,有
7.重复步骤2至6,直到i=N1和j=N2为止,停止重复;
8.对得到的所有Ii,j nf(i=0,…,N1;j=0,…,N2)求平均,得到预降噪图像:
9.对上一步得到的预降噪图像进一步进行Wiener滤波处理,得到最终的降噪结果
上述的第4步骤中的TB值的具体估计步骤为:
①对于噪声标准差σn,采用鲁棒性的中值估计,
其中Shf (Ki)(i=1…LK)为最高频率系数;
②由 有
其中, Shf (k i)是所考虑的高频系数;
③因此可得,阈值参数
本发明方法与现有技术相比较,具有如下显而易见的突出实质性特点和显著优点:
该发明旨在提供一种图像组合降噪方法,首先通过Contourlet变换对输入的带噪图像进行多尺度、多方向的稀疏分解,然后根据Contourlet变换域系数服从广义高斯分布,在Contourlet域进行Bayes收缩阈值法降噪,并通过Contourlet反变换得到预降噪图像,最后,采用Wiener滤波法对预降噪图像进行进一步降噪处理,以提高图像的恢复精度。具体特点和优点为:
(1)针对现有最具有代表性的小波域阈值降噪方法中小波变换的缺点------不能有效地表示图像中的二位或高维奇异性,将Contourlet变换应用到图像降噪中,进行多尺度、多方向分解,为后续降噪过程提供稀疏的图像描述系数。
(2)对现有图像降噪技术存在的不足,提出了一种图像组合降噪方法,即基于移不变Contourlet变换域降噪和Wiener滤波的图像组合降噪方法。
(3)本发明方法在Contourlet变换域的图像降噪阶段,针对图像的Contourlet域系数服从广义高斯分布(GGD),满足Bayes估计方法的假定条件——信号服从广义高斯分布,采用基于Bayes估计阈值,进行降噪,提高了降噪性能。
(4)本发明方法针对Contourlet变换域的图像降噪方法不能完全去除噪声,降噪图像中仍然残留少量噪声,采取Wiener滤波法进一步降噪处理。
本发明提供的图像降噪方法能提高降噪图像的质量、提供更加全面、准确的目标和背景信息,达到较理想的降噪效果。在军事领域和非军事领域如光学成像、目标检测、安全监控等系统中均有广泛应用前景。
附图说明
图1为本发明一个实施例的图像降噪方法框图。
图2是图1示例降噪结果照片图。图中,(a)到(e)为输入图像受到不同噪声污染情况下的降噪结果,噪声强度分别为10、20、30、40和50。每一行中,第一幅图为受噪声污染的输入,第二幅图为采用一种Contourlet变换域的图像降噪方法降噪后的降噪图像,第三幅图为采用本发明方法降噪后的降噪图像。
具体实施方式
本发明的一个优选实施例结合附图祥述如下:
本发明旨在提供一种图像组合降噪方法,如图1所示。该方法先通过Contourlet变换对输入的带噪图像进行多尺度、多方向的稀疏分解,然后根据Contourlet变换域系数服从广义高斯分布,在Contourlet域进行Bayes收缩阈值法降噪,并通过Contourlet反变换得到预降噪图像,最后,采用Wiener滤波法对预降噪图像进行进一步降噪处理,得到最终的降噪图像,达到图像降噪的目的。
具体步骤为:
1.初始化设置。令i=0,j=0,设定行方向和列方向上的最大平移量N1和N2.同时设定Contourlet变换的中LP分解层数K和每层中的方向分解数Lk;
2.对输入的带噪图像I在行和列方向上进行有效位移量的循环平移,得到平移图像
Sij=Cij(I),
其中i∈(0,N1)和j∈(0,N2)分别为行方向和列方向上的平移量;
3.对得到的平移图像Sij进行多尺度、多方向的Contourlet稀疏分解,即
其中T(·)为Contourlet变换。从而得到一幅低频子图像Sif和一系列具有不同分辨率的高频子图像Shf (k,l),其中k∈(1,K)和l ∈(1,Lk)标明子图像位于第k层LP(拉普拉斯塔式分解)的第l方向;
4.对Contourlet变换后的高频子图像Shf (kl)进行阈值降噪处理,得到降噪子图像,
其中,Λ(·)为阈值函数,TB为阈值参数。根据Contourlet域系数服从广义高斯分布,因此,本实施例选用Bayes阈值 具体估计步骤为:
①对于噪声标准差σn,采用鲁棒性的中值估计,
其中Shf (Ki)(i=1…LK)为最高频率系数;
②由 有
其中, Shf (k i)是所考虑的高频系数;
③因此可得,阈值
5.对第4步中得到的所有降噪高频子图像SDhf (k,l)和第3步中得到的低频子图像Sif实施Contourlet逆变换,得到在行方向和列方向上分别平移i和j后的降噪图像,
其中,T-1(·)为Contourlet逆变换;
6.对第5步中得到的图像Si,j nf进行相应平移量的逆向循环平移,有
7.重复步骤2至6,直到i=N1和j=N2为止,停止重复;
8.对得到的所有Ii,j kf(i=0,…,N1;j=0,…,N2)求平均,得到预降噪图像:
从图2可以看出,从视觉效果上看,本发明方法明显优于基于移不变Contourlet变换域的Bayes收缩阈值法。在降噪的过程中,后者只是对经过Contourlet变换后得到的高频系数进行Bayes阈值去噪,而没有考虑低频部分,故降噪后的图像仍然残留少量噪声。而本发明方法在此基础之上进一步采用具有MSE意义下最优的估计器——Wiener滤波器实施滤波,从而在保护图像细节信息的同时,进一步提高了降噪图像的PSNR,进一步降低了降噪图像的MSE。
表1给出了本发明降噪方法降噪结果的客观评价指标。
表中,方法1指基于移不变Contourlet变换域的Bayes收缩阈值降噪方法采用了峰值信噪比(PSNR)和最小均方误差(MSE)来衡量降噪图像的质量,进而评价降噪方法的优劣。
从表中可以看出,本图像降噪方法无论是在PSNR方面,还是在MSE方面,均能取得很好的降噪效果,有效地降低图像中的噪声信号,提高图像质量。
另外,从表中不难发现,随着噪声水平的升高,经本发明方法降噪后得到的降噪图像指标(PSNR和MSE)较基于移不变Contourlet变换域的Bayes收缩阈值法提高幅度不断增大,表明本发明降噪方法的优势越发明显。在视觉效果方面,也能从图2得出相同的结论。
总之,无论是从人眼视觉效果,还是从客观评价指标,均表明本发明方法更好地降低图像中的噪声信号,保护了图像中的重要细节信息,提高了图像的质量。
表1不同噪声水平下降噪结果的评价指标比较
Claims (1)
1. 一种图像组合降噪方法,其特征在于采用基于移不变Contourlet变换域降噪和Wiener滤波的图像组合降噪方法,该方法是先通过Contourlet变换对输入的带噪图像进行多尺度、多方向的稀疏分解,然后根据Contourlet变换域系数服从广义高斯分布,在Contourlet域进行Bayes收缩阈值法降噪,并通过Contourlet反变换得到预降噪图像,最后,采用Wiener滤波法对预降噪图像进行进一步降噪处理,得到最终的降噪图像,达到图像降噪的目的;具体步骤为:
1)初始化设置:令i=0,j=0,设定行方向和列方向上的最大平移量N1和N2;同时设定Contourlet变换中的拉普拉斯塔式分解层数K和每层中的方向分解数Lk;
2)对输入的带噪图像I在行和列方向上进行有效位移量的循环平移,得到平移图像
Sij=Cij(I),
其中i∈(0,N1)和j∈(0,N2)分别为行方向和列方向上的平移量;
3)对得到的平移图像Sij进行多尺度、多方向的Contourlet稀疏分解,即
其中T(·)为Contourlet变换;从而得到一幅低频子图像Slf和一系列具有不同分辨率的高频子图像Shf (k l),其中k∈(1,K)和l∈(1,Lk)标明子图像位于第k层拉普拉斯塔式分解的第l方向;
4)对Contourlet变换后的高频子图像Shf (k l)进行阈值降噪处理,得到降噪子图像,
其中,Λ(·)为阈值函数,Thr为阈值,根据Contourlet变换域系数服从广义高斯分布,选用Bayes阈值
5)对第4)步骤中得到的所有降噪高频子图像SDhf (k l)和第3)步骤中得到的低频子图像Slf实施Contourlet逆变换,得到在行方向和列方向上分别平移i和j后的降噪图像,
其中,T-1(·)为Contourlet逆变换;
6)对第5)步骤中得到的图像Si,j nf进行相应平移量的逆向循环平移,有
7)重复步骤2)至6),直到i=N1和j=N2为止,停止重复;
8)对得到的所有Ii,j nf(i=0,…,N1;j=0,…,N2)求平均,得到预降噪图像:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100307462A CN100417191C (zh) | 2006-09-01 | 2006-09-01 | 一种图像组合降噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100307462A CN100417191C (zh) | 2006-09-01 | 2006-09-01 | 一种图像组合降噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1917577A CN1917577A (zh) | 2007-02-21 |
CN100417191C true CN100417191C (zh) | 2008-09-03 |
Family
ID=37738459
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100307462A Expired - Fee Related CN100417191C (zh) | 2006-09-01 | 2006-09-01 | 一种图像组合降噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100417191C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101639934B (zh) * | 2009-09-04 | 2012-07-04 | 西安电子科技大学 | 基于轮廓波域块隐马尔可夫模型sar图像去噪方法 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090154567A1 (en) * | 2007-12-13 | 2009-06-18 | Shaw-Min Lei | In-loop fidelity enhancement for video compression |
CN101477679B (zh) * | 2009-01-16 | 2012-07-25 | 西安电子科技大学 | 基于轮廓波Contourlet变换的图像去噪方法 |
CN101719268B (zh) * | 2009-12-04 | 2011-10-19 | 西安电子科技大学 | 基于改进Directionlet域的广义高斯模型图像去噪方法 |
CN101799916A (zh) * | 2010-03-16 | 2010-08-11 | 刘国传 | 基于贝叶斯估计的生物芯片图像小波去噪方法 |
CN101847257B (zh) * | 2010-06-10 | 2012-06-20 | 上海电力学院 | 基于非局部均值与多级定向图像的图像降噪方法 |
CN102316282B (zh) * | 2011-09-20 | 2014-05-07 | 中国科学院理化技术研究所 | 一种基于光学杜比的图像降噪装置 |
CN106161875B (zh) * | 2015-03-25 | 2019-02-15 | 瑞昱半导体股份有限公司 | 图像处理装置与方法 |
CN114037747A (zh) * | 2021-11-25 | 2022-02-11 | 佛山技研智联科技有限公司 | 图像特征提取方法、装置、计算机设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005032122A1 (en) * | 2003-09-29 | 2005-04-07 | Samsung Electronics Co., Ltd. | Denoising method and apparatus |
CN1731388A (zh) * | 2004-12-30 | 2006-02-08 | 中国科学院长春光学精密机械与物理研究所 | 超声信号的数字信号处理方法 |
-
2006
- 2006-09-01 CN CNB2006100307462A patent/CN100417191C/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005032122A1 (en) * | 2003-09-29 | 2005-04-07 | Samsung Electronics Co., Ltd. | Denoising method and apparatus |
CN1731388A (zh) * | 2004-12-30 | 2006-02-08 | 中国科学院长春光学精密机械与物理研究所 | 超声信号的数字信号处理方法 |
Non-Patent Citations (4)
Title |
---|
基于Contourlet变换的图像去噪算法. 邓承志,汪胜前,钟华,刘祝华,邹道文.电视技术,第10期. 2004 |
基于Contourlet变换的图像去噪算法. 邓承志,汪胜前,钟华,刘祝华,邹道文.电视技术,第10期. 2004 * |
小波分析与维纳滤波相结合的消噪方法研究. 胡耀斌,陈艾华,张春良.电力系统通信,第27卷第162期. 2006 |
小波分析与维纳滤波相结合的消噪方法研究. 胡耀斌,陈艾华,张春良.电力系统通信,第27卷第162期. 2006 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101639934B (zh) * | 2009-09-04 | 2012-07-04 | 西安电子科技大学 | 基于轮廓波域块隐马尔可夫模型sar图像去噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1917577A (zh) | 2007-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100417191C (zh) | 一种图像组合降噪方法 | |
CN101847257B (zh) | 基于非局部均值与多级定向图像的图像降噪方法 | |
CN105913393B (zh) | 一种自适应小波阈值图像去噪方法及装置 | |
CN100550978C (zh) | 一种保持边缘的自适应图像滤波方法 | |
CN103049892B (zh) | 基于相似块矩阵秩最小化的非局部图像去噪方法 | |
Naimi et al. | Medical image denoising using dual tree complex thresholding wavelet transform and Wiener filter | |
CN101950414B (zh) | 自然图像非局部均值去噪方法 | |
CN100433062C (zh) | 一种Contourlet变换域的图像降噪方法 | |
CN101944230B (zh) | 基于多尺度的自然图像非局部均值去噪方法 | |
CN100433795C (zh) | 基于变换域数学形态学的图像降噪方法 | |
CN101477679A (zh) | 基于轮廓波Contourlet变换的图像去噪方法 | |
CN103455991A (zh) | 一种多聚焦图像融合方法 | |
CN104103041B (zh) | 超声图像混合噪声自适应抑制方法 | |
CN102890820A (zh) | 基于shearlet变换和维纳滤波的图像去噪方法 | |
CN102663702A (zh) | 基于区域划分的自然图像去噪方法 | |
CN100433046C (zh) | 一种基于稀疏变换的图像盲源分离方法 | |
CN103077507B (zh) | 基于Beta算法的多尺度SAR图像降噪方法 | |
CN102722879A (zh) | 基于目标提取和三维块匹配去噪的sar图像去斑方法 | |
CN103971345A (zh) | 一种基于改进双边滤波的图像去噪方法 | |
CN110852963B (zh) | 一种基于相关最大性的湍流退化图像复原方法 | |
CN102314675A (zh) | 基于小波高频的贝叶斯去噪方法 | |
CN104754183A (zh) | 一种实时监控视频自适应滤波方法及其系统 | |
CN101561879A (zh) | 基于曲波表示的图像欠定盲分离方法 | |
CN102622731B (zh) | 基于二维Otsu的轮廓波域维纳滤波图像去噪方法 | |
CN104182944A (zh) | 基于曲波与小波变换相串联的光学图像去噪方法 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080903 Termination date: 20110901 |