CN102226917A - 一种基于非下采样Contourlet扩散的图像增强方法 - Google Patents
一种基于非下采样Contourlet扩散的图像增强方法 Download PDFInfo
- Publication number
- CN102226917A CN102226917A CN2011101331737A CN201110133173A CN102226917A CN 102226917 A CN102226917 A CN 102226917A CN 2011101331737 A CN2011101331737 A CN 2011101331737A CN 201110133173 A CN201110133173 A CN 201110133173A CN 102226917 A CN102226917 A CN 102226917A
- Authority
- CN
- China
- Prior art keywords
- matrix
- nsct
- obtains
- coefficients
- image
- 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.)
- Pending
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及一种基于非下采样Contourlet扩散的噪声图像增强方法,基于非下采样Contourlet变换对边缘/轮廓的良好表示特性,利用非线性各向异性扩散方法有效消减增强图像中存在的伪Gibbs现象,在对图像的边缘特征和对比度增强的同时,对其中的噪声进行抑制,进一步提高图像质量。实验结果表明,本发明与传统的基于wavelet和非下采样Contourlet的方法相比能得到较好的增强效果。
Description
技术领域
本发明涉及一种基于非下采样Contourlet扩散的图像增强方法,属于数字图像处理领域。
背景技术
图像增强是按特定的需要突出一幅图像中的某些信息,同时削弱或去除某些不需要的信息的方法。其目的是提高图像在某种特定运用下的清晰度和对比度,从而改善图像的质量,使得处理的结果更符合人的视觉感官***或更容易被机器所识别。传统的图像增强方法主要有直方图均衡化、反锐化掩模等空域法,以及通过傅里叶变换增强感兴趣频率成分的频域法。这些方法增强图像对比度的同时也会放大噪声,使得图像细节信息被噪声淹没。
近年来,基于小波变换的图像增强技术取得很好的效果,然而,由一维小波拓展为二维小波变换只具有有限的方向,用它来对图像进行分析时,是以“点”为单位来捕捉图像的特征,不能“最优”表示含“线”或“面”奇异的高维函数。为了解决小波变换在这些方面的局限性,2002年,M.N.Do、Martin Vetterli提出了Contourlet变换,它有效地有效的弥补了上述不足,是一种“真正”二维图像表示方法,这种方法可以很好地抓住图像的几何结构,“很好”刻画图像边缘,被广泛应用于去噪,增强等领域。同时,由于Contourlet变换不具有平移不变性,2006年,具有平移不变特性的Contourlet变换-非下采样Contourlet变换(Nonsubsampled Contourlet Transform,NSCT)被M.N.Do等人提出,逐渐表现出其自身在去噪、增强方面的优势。
现阶段,利用NSCT进行图像增强的方法主要是先对含噪声的图像进行NSCT变换,在变换域通过一个合适的阈值去除主要由噪声组成的小于阈值的系数,然后进行NSCT逆变换得到增强后的图像。这种方法虽然很好的去除了噪声,但是,在小系数置零的过程中,同时也去除了部分有用的系数,从而造成了去噪后的图像具有一定的伪Gibbs现象。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于非下采样Contourlet扩散的图像增强方法,在对图像的边缘特征和对比度增强的同时,对其中的噪声进行抑制,并能有效消减伪Gibbs现象,进一步提高图像质量。
技术方案
一种基于非下采样Contourlet扩散的图像增强方法,其特征在于步骤如下:
步骤1:对噪声图像灰度值矩阵进行给定尺度参数下的非下采样Contourlet正变换,得到不同分解尺度下各个方向的NSCT系数矩阵;所述的尺度参数为:分解尺度数和方向数,尺度数取3-4,方向数取2-4;
步骤2:利用增益函数与步骤1得到的NSCT系数矩阵中的系数相乘实现系数矩阵的非线性变换,得到非线性变换后不同分解尺度下各个方向的NSCT系数矩阵;所述的增益函数为:
其中,x为步骤1得到的NSCT系数矩阵中的系数,p∈(0,1],s∈(0,1];T1、T2、T3的计算公式分别为:
T2=s3T1
T3=s4T2
其中,j为分解尺度,σ为按σ=median(|SHH|)/0.6745计算得到的噪声图像灰度值矩阵的对数矩阵的噪声标准差,median(·)表示取中值,|·|表示取模,SHH是对噪声图像灰度值矩阵进行一级小波分解得到的对角方向高频子带小波系数矩阵;σx为输入x的二范数;s1∈[1,5],s2∈[1,5],s3∈[2,3],s4∈[1.5,2.5];
步骤3:对步骤2得到的NSCT系数矩阵进行NSCT逆变换,得到NSCT逆变换结果矩阵,逆变换的尺度参数与正变换的尺度参数一致;
步骤4:对步骤3得到的NSCT逆变换结果矩阵进行NSCT扩散,得到最终的增强图像灰度值矩阵;
所述的NSCT扩散过程的具体步骤为:
步骤a:对输入的步骤3得到的NSCT逆变换结果矩阵进行NSCT正变换,得到不同分解尺度下各个方向的NSCT系数矩阵;
步骤b:利用逆阈值化函数对步骤(1)得到的NSCT系数矩阵进行逆阈值化处理,得到不同分解尺度下不同方向的NSCT系数矩阵;所述的逆阈值化函数为
步骤c:对步骤b得到的NSCT系数矩阵进行NSCT逆变换,得到NSCT逆变换结果矩阵;
步骤d:将步骤a得到的NSCT变换结果矩阵与步骤b得到的逆阈值化NSCT矩阵相减,得到差图像灰度值矩阵;
步骤e:以步骤c得到的NSCT逆变换结果矩阵为初始条件进行非线性各向异性扩散,得到平滑后的图像灰度值矩阵;
步骤f:将步骤d得到的差图像灰度值矩阵与步骤e得到的平滑图像灰度值矩阵相叠加,得到最终的NSCT扩散后的图像灰度值矩阵。
有益效果
本发明提出的一种基于非下采样Contourlet扩散的图像增强方法,各向同性扩散过程源于自然界中的热扩散,而非线性各向异性扩散是Perona及Malik为了避免各向同性扩散中出现的边缘模糊和定位问题提出的一种基于偏微分方程(Partial Differential Equation,PDE)的方法。该方法可同步实现噪声消除和边缘保持,但对噪声非常敏感。
附图说明
图1:本发明基于非下采样Contourlet扩散的噪声图像增强方法的流程图。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明中所用图像为定义在R2上的含有高斯噪声灰度图像I0,每幅图像都对应一个矩阵,图像中的每一个像素对应矩阵中相应的一个元素;
步骤1:对输入的噪声图像I0进行NSCT正变换,得到不同分解尺度j下不同方向l的NSCT系数矩阵C0(j,l),其中j=1,2,…,J,J为最细分解尺度,l=1,…,Lj,Lj为第j个尺度下的方向数;
步骤2:对输入图像按σ=median(|SHH|)/0.6745进行噪声标准差估计,其中median(·)表示取矩阵中所有系数的中值,|·|表示取模,SHH是对图像I0进行一级小波分解得到的对角方向高频子带(HH子带)小波系数矩阵;
步骤4:利用增益函数对对步骤1中的C0(j,l)在除最粗尺度以外的不同分解尺度各个方向下的NSCT系数矩阵分别进行拉伸/收缩,得到非线性变换后不同分解尺度各个方向的NSCT系数矩阵;
步骤5:对步骤4得到的NSCT系数矩阵进行NSCT逆变换,得到NSCT逆变换矩阵,即初始增强图像Icon;
步骤6:对步骤5得到的初始增强图像Icon进行NSCT正变换,得到不同分解尺度下不同方向的NSCT系数矩阵Ccon(j,l);
步骤7:对NSCT系数矩阵Ccon(j,l)在除最粗尺度外的不同分解尺度不同方向下进行逆阈值化处理,得到扩散系数矩阵Cdif(j,l),逆阈值化函数为
步骤8:对扩散系数矩阵Cdif(j,l)进行NSCT逆变换,得到NSCT逆变换矩阵,即扩散图像Idif,并用算子P(·)表示该过程,记为Idif=P(Cdif(j,l)),进而可得到差图像
步骤9:对扩散图像Idif进行非线性各向异性扩散,得到平滑后图像Dif(Idif),这里Dif(·)表示表示一个迭代次数为t的扩散过程Dif(I)=I(t+1)=I(t)+div(g(|▽P(I(t))|)▽P(I(t)));
步骤10:将步骤9得到的平滑后图像Dif(Idif)与差图像相叠加,得到最终的增强图像
Claims (1)
1.一种基于非下采样Contourlet扩散的图像增强方法,其特征在于步骤如下:
步骤1:对噪声图像灰度值矩阵进行给定尺度参数下的非下采样Contourlet正变换,得到不同分解尺度下各个方向的NSCT系数矩阵;所述的尺度参数为:分解尺度数和方向数,尺度数取3-4,方向数取2-4;
步骤2:利用增益函数与步骤1得到的NSCT系数矩阵中的系数相乘实现系数矩阵的非线性变换,得到非线性变换后不同分解尺度下各个方向的NSCT系数矩阵;所述的增益函数为:
其中,x为步骤1得到的NSCT系数矩阵中的系数,p∈(0,1],s∈(0,1];T1、T2、T3的计算公式分别为:
T2=s3T1
T3=s4T2
其中,j为分解尺度,σ为按σ=median(|SHH|)/0.6745计算得到的噪声图像灰度值矩阵的对数矩阵的噪声标准差,median(·)表示取中值,|·|表示取模,SHH是对噪声图像灰度值矩阵进行一级小波分解得到的对角方向高频子带小波系数矩阵;σx为输入x的二范数;s1∈[1,5],s2∈[1,5],s3∈[2,3],s4∈[1.5,2.5];
步骤3:对步骤2得到的NSCT系数矩阵进行NSCT逆变换,得到NSCT逆变换结果矩阵,逆变换的尺度参数与正变换的尺度参数一致;
步骤4:对步骤3得到的NSCT逆变换结果矩阵进行NSCT扩散,得到最终的增强图像灰度值矩阵;
所述的NSCT扩散过程的具体步骤为:
步骤a:对输入的步骤3得到的NSCT逆变换结果矩阵进行NSCT正变换,得到不同分解尺度下各个方向的NSCT系数矩阵;
步骤b:利用逆阈值化函数对步骤(1)得到的NSCT系数矩阵进行逆阈值化处理,得到不同分解尺度下不同方向的NSCT系数矩阵;所述的逆阈值化函数为
步骤c:对步骤b得到的NSCT系数矩阵进行NSCT逆变换,得到NSCT逆变换结果矩阵;
步骤d:将步骤a得到的NSCT变换结果矩阵与步骤b得到的逆阈值化NSCT矩阵相减,得到差图像灰度值矩阵;
步骤e:以步骤c得到的NSCT逆变换结果矩阵为初始条件进行非线性各向异性扩散,得到平滑后的图像灰度值矩阵;
步骤f:将步骤d得到的差图像灰度值矩阵与步骤e得到的平滑图像灰度值矩阵相叠加,得到最终的NSCT扩散后的图像灰度值矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101331737A CN102226917A (zh) | 2011-05-20 | 2011-05-20 | 一种基于非下采样Contourlet扩散的图像增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101331737A CN102226917A (zh) | 2011-05-20 | 2011-05-20 | 一种基于非下采样Contourlet扩散的图像增强方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102226917A true CN102226917A (zh) | 2011-10-26 |
Family
ID=44807882
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011101331737A Pending CN102226917A (zh) | 2011-05-20 | 2011-05-20 | 一种基于非下采样Contourlet扩散的图像增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102226917A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102496144A (zh) * | 2011-11-15 | 2012-06-13 | 哈尔滨工程大学 | 基于hsv色彩空间的nsct水声图像增强方法 |
CN102637294A (zh) * | 2012-03-03 | 2012-08-15 | 西北工业大学 | 基于非下采样Contourlet变换和改进的总变分的图像增强方法 |
CN102708545A (zh) * | 2012-04-13 | 2012-10-03 | 泰山学院 | 一种非局部变换域图像增强方法 |
CN103295204A (zh) * | 2013-06-20 | 2013-09-11 | 河海大学常州校区 | 一种基于非下采样轮廓波变换的图像自适应增强方法 |
CN104392442A (zh) * | 2014-11-18 | 2015-03-04 | 西北工业大学 | 基于非下采样Contourlet变换与主动轮廓的遥感图像变化检测方法 |
US9311564B2 (en) | 2012-10-05 | 2016-04-12 | Carnegie Mellon University | Face age-estimation and methods, systems, and software therefor |
CN110246106A (zh) * | 2019-06-22 | 2019-09-17 | 福州大学 | 基于量子和声搜索模糊集的nsst域浮选泡沫图像增强及去噪方法 |
CN115187465A (zh) * | 2022-05-30 | 2022-10-14 | 中国矿业大学 | 基于nsct和图像匹配的浮选泡沫图像增强方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7315656B2 (en) * | 2003-05-22 | 2008-01-01 | The Boeing Company | Methods and apparatus for enhanced viewing of aerial refueling operations |
CN101833754A (zh) * | 2010-04-15 | 2010-09-15 | 青岛海信网络科技股份有限公司 | 图像增强方法及*** |
CN101882305A (zh) * | 2010-06-30 | 2010-11-10 | 中山大学 | 一种图像增强处理的方法 |
-
2011
- 2011-05-20 CN CN2011101331737A patent/CN102226917A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7315656B2 (en) * | 2003-05-22 | 2008-01-01 | The Boeing Company | Methods and apparatus for enhanced viewing of aerial refueling operations |
CN101833754A (zh) * | 2010-04-15 | 2010-09-15 | 青岛海信网络科技股份有限公司 | 图像增强方法及*** |
CN101882305A (zh) * | 2010-06-30 | 2010-11-10 | 中山大学 | 一种图像增强处理的方法 |
Non-Patent Citations (2)
Title |
---|
ARTHUR L. DA CUNHA ET AL: "The nonsampled contourlet transform:theory, design, and applications", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 * |
董鸿燕 等: "基于Contourlet变换的带噪图像增强方法", 《信号处理》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102496144A (zh) * | 2011-11-15 | 2012-06-13 | 哈尔滨工程大学 | 基于hsv色彩空间的nsct水声图像增强方法 |
CN102637294A (zh) * | 2012-03-03 | 2012-08-15 | 西北工业大学 | 基于非下采样Contourlet变换和改进的总变分的图像增强方法 |
CN102708545A (zh) * | 2012-04-13 | 2012-10-03 | 泰山学院 | 一种非局部变换域图像增强方法 |
US9311564B2 (en) | 2012-10-05 | 2016-04-12 | Carnegie Mellon University | Face age-estimation and methods, systems, and software therefor |
CN103295204A (zh) * | 2013-06-20 | 2013-09-11 | 河海大学常州校区 | 一种基于非下采样轮廓波变换的图像自适应增强方法 |
CN103295204B (zh) * | 2013-06-20 | 2015-12-02 | 河海大学常州校区 | 一种基于非下采样轮廓波变换的图像自适应增强方法 |
CN104392442A (zh) * | 2014-11-18 | 2015-03-04 | 西北工业大学 | 基于非下采样Contourlet变换与主动轮廓的遥感图像变化检测方法 |
CN110246106A (zh) * | 2019-06-22 | 2019-09-17 | 福州大学 | 基于量子和声搜索模糊集的nsst域浮选泡沫图像增强及去噪方法 |
CN115187465A (zh) * | 2022-05-30 | 2022-10-14 | 中国矿业大学 | 基于nsct和图像匹配的浮选泡沫图像增强方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102226917A (zh) | 一种基于非下采样Contourlet扩散的图像增强方法 | |
Jain et al. | A survey of edge-preserving image denoising methods | |
CN106530237A (zh) | 一种图像增强方法 | |
CN101777179B (zh) | 一种基于偏微分方程的双树复小波图像去噪方法 | |
CN102509263B (zh) | 基于sar图像局部统计特性的k-svd相干斑抑制方法 | |
CN101901476A (zh) | 基于NSCT域边缘检测和Bishrink模型的SAR图像去噪方法 | |
CN104103040B (zh) | 一种图像增强的方法 | |
Shi et al. | An image denoising method based on multiscale wavelet thresholding and bilateral filtering | |
CN102637294A (zh) | 基于非下采样Contourlet变换和改进的总变分的图像增强方法 | |
Gopinathan et al. | Wavelet and FFT Based Image Denoising Using Non-Linear Filters. | |
CN102663679A (zh) | 基于Shearlet收缩和改进TV模型的图像去噪方法 | |
CN104182939A (zh) | 一种医疗影像图像细节增强方法 | |
CN108492268A (zh) | 基于小波系数融合的低照度图像增强算法 | |
CN102081790B (zh) | 一种基于非线性Curvelet扩散的噪声图像增强方法 | |
CN104616259A (zh) | 一种噪声强度自适应的非局部均值图像去噪方法 | |
CN111192204A (zh) | 一种图像增强方法、***和计算机可读存储介质 | |
CN103745442B (zh) | 基于非局部小波系数收缩的图像去噪方法 | |
CN103310423B (zh) | 一种矿井图像增强方法 | |
Jumiawi et al. | Image spectrum segmentation for lowpass and highpass filters | |
CN102184530B (zh) | 基于灰关联阈值的图像去噪方法 | |
Li et al. | Image de-noising of locust sections based on adaptive wavelet and partial differential equation method | |
CN104182944A (zh) | 基于曲波与小波变换相串联的光学图像去噪方法 | |
Thriveni | Satellite image enhancement using discrete wavelet transform and threshold decomposition driven morphological filter | |
Zhang et al. | Diffusion scheme using mean filter and wavelet coefficient magnitude for image denoising | |
Sharma et al. | Image enhancement using fractional partial differential equation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20111026 |