CN102332155A - 一种图像增强方法 - Google Patents
一种图像增强方法 Download PDFInfo
- Publication number
- CN102332155A CN102332155A CN201110318806A CN201110318806A CN102332155A CN 102332155 A CN102332155 A CN 102332155A CN 201110318806 A CN201110318806 A CN 201110318806A CN 201110318806 A CN201110318806 A CN 201110318806A CN 102332155 A CN102332155 A CN 102332155A
- Authority
- CN
- China
- Prior art keywords
- image
- coefficient
- frequency sub
- equalization
- 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.)
- Pending
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及一种基于DWT的整数小波系数均衡的图像增强算法,其包括如下:原始图像进行单层离散二维小波分解,得到低频子带的分解系数cA,和分别对应于水平、垂直和对角三个方向的三个高频子带分解系数cH,cV,cD;对所述低频子带分解系数cA做均衡化处理,得到均衡化的新系数;利用所述均衡化的新系数和所述高频子带分解系数cH,cV,cD重构图像;本发明的图像增强方法既保留了小波分析能够将信号的细节放大、利于分析处理、能够有效的增强图像的同时减小图像噪声的优点,又克服了采用直方图均衡化算法在噪声图像增强过程中会大大增强噪声,以及在采用同态滤波算法增强图像的过程中损失大量细节的缺点。
Description
技术领域
本发明涉及一种图像增强技术,属于图像处理领域,特别涉及一种基于DWT的整数小波系数均衡的图像增强方法。
背景技术
图像增强是数字图像处理的重要内容之一,其目的是突出图像中的“有用”信息,扩大图像中不同物体特征之间的差别,为图像的信息提取及其他图像分析技术奠定良好的基础。在图像处理领域,图像增强问题主要通过空域和频域处理两种方法来解决。空域方法通过在图像点上作用算子或掩码来解决,频域方法通过修改傅立叶变换系数来解决。这两种方法的优劣很明显,空域方法方便快速但会丢失很多之间的相关信息,频域方法可以很详细的分离出点之间的相关,但要做两次数量级为nlogn的傅立叶变换和逆变换的操作,计算量大得多。
直方图均衡化的基本思想是将原始图像的不均衡的直方图变换为均匀分布的形式。即将输入图像转换为在每一灰度级上都有相同像素点数(即输出的直方图是平坦的,分布是均匀的)。这使得图像像素值动态范围变宽,对比度得到有效拉伸,从而达到图像增强的效果。直方图均衡增强的缺点是若在对噪声图像增强的同时,一点小小的噪声也会被放大。
同态滤波是一种在频域中同时将图像亮度范围压缩和对比度增强的方法。其基本思想是将非线性问题转化成线性问题处理,同态滤波增强的缺点是在噪声图像的增强过程中会损失了大量的图像细节。
发明内容
本发明要解决的技术问题是提供一种算法简单、效果明显,特别是对噪声图像有显著增强效果的图像增强算法。
为解决上述技术问题,本发明提供基于DWT的整数小波系数均衡的图像增强算法,其包括如下步骤:
a、原始图像进行单层离散二维小波分解,得到低频子带的分解系数cA,和分别对应于水平、垂直和对角三个方向的三个高频子带分解系数cH,cV,cD;
b、对所述低频子带分解系数cA做均衡化处理,得到均衡化的新系数;
c、利用所述均衡化的新系数和所述高频子带分解系数cH,cV,cD重构图像;
所述步骤b中的所述均衡化处理的方法,包括如下步骤:
b1、对所述低频子带分解系数cA向下取整,即对作为二维数组的低频子带分解系数cA中的各元素k的值按四舍五入的原则取整,得到新的低频子带分解系数cA’;
b2、统计上述新的低频子带分解系数cA’中的各元素k的对应于不同取值所具有的个数n(k),并且求出各元素k中的最大值Kmax
b4、求出所述均衡化的新系数g (i,j)。
所述的元素k≥0。
0≤k≤Kmax,其中m,n为图像的行数与列数,i,j分别为图像矩阵的行、列变量。
本发明具有如下优点:
1、本发明涉及的图像增强方法仅对小波单层分解后的低频子带部分的系数进行运算,相比现有技术提供的方案,大大减少了图像增强过程中的运算量;本发明的图像增强方法既保留了小波分析能够将信号的细节放大、利于分析处理、能够有效的增强图像的同时减小图像噪声的优点,又克服了采用直方图均衡化算法在噪声图像增强过程中会大大增强噪声,以及在采用同态滤波算法增强图像的过程中损失大量细节的缺点。
2、本发明涉及的算法在增强效果上,特别是在噪声图像的增强效果上优于现有技术提供的方案。
附图说明
图1为本发明的原始图像增强对比结果;
图2为本发明的噪声图像增强对比结果;
图3为本发明的原始图像增强对比直方图。
具体实施方式
为了使本发明的内容更容易的被理解,下面根据具体实施例并结合附图,对本发明作进一步详细的说明:
本实施例采用的图像增强对象为标准的lena图像(如图1-a)与加入方差为0.003高斯噪声的lena图像(如图2-a)。(注:Lena图像是图像处理领域广泛使用的标准测试图像)。
本实施例需要的小波是光滑的、正交的、对称的,这样的小波处理图像的好处是速度快、图像重构精确性高、避免图像处理中发生相移。满足上述条件的有db小波系与sym小波系,同db小波系相比,其它性质同db小波一致的sym小波系有更好的对称性,更适合于图像处理,因此本实施例选取syml小波,对图像作单层小波分解。
本实施例具有以下具体步骤:
(1)利用单层离散二维小波分解函数[cA,cH,cV,cD]=dwt2(I,’syml’)(注:本函数为matlab中提供的标准二维单尺度小波变换函数,具体可参看国防工业出版社2006-6-1出版的《MATLAB 7.0图形图像处理》(ISBN:9787118045437)对标准的lena图像(图1-a)进行syml小波单层分解,得到一个低频子带分解系数cA和分别对应于水平、垂直、对角方向的三个高频子带的分解系数cA、cV、cD;
(2)由于图像的轮廓信息主要集中在低频部分,而细节和噪声则在高频部分,因此仅统计低频子带系数cA向下取整,得到新的低频子带系数cA’,统计上述新的低频子带分解系数cA’中的各元素k的对应于不同取值所具有的个数n(k),并且求出各元素k的最大值Kmax;
(3)对所述元素k在不同取值下的统计个数n(k)进行累积求和:
(5)利用二维小波分解系数的直接重构函数Y=upcoef2(’a’,g(k),’syml’,1)(注:本函数为matlab中提供的标准二维单尺度小波变换函数,具体可参看国防工业出版社2006-6-1出版的《MATLAB 7.0图形图像处理》(ISBN:9787118045437)中的13.3.1图像的二维小波分解与重构)重构基于所述均衡化的新系数和所述高频子带分解系数cH,cV,cD的图像(图1-d);
(6)重复上述步骤对加入方差为0.003高斯噪声的lena图像(图2-a)做图像增强,得到新的图像(图2-d)。
为检验本发明涉及的算法在图像增强与去噪方面的能力,本实施例中引入了常用的同态滤波增强算法、直方图均衡增强算法与本发明中所述的算法作对比仿真并且选用均方误差(MSE)、峰值信噪比(PSNR)和对比度增量(Δ)对本发明算法进行客观评价。
(1)均方误差:
(2)峰值信噪比:PSNR=101g[f2 max(i,j)/MSE];
(3)对比度增量:
上述算法(3)中为增强后图像局部对比度均值,Cf为原始图像局部对比度均值,对比度增量Δ为原始图像与增强后图像局部对比度之比,局部对比度以3×3的滑动窗口,按照(xmax-xmin)/(xmax+xmin)计算每个窗口的局部对比度,然后取其平均值,对比度增量越大说明增强效果越好。
实验中三种增强算法都对原始图像(图1-a)、加噪后的图像(图2-a)处理;(图1-b)、(图2-b)为同态滤波增强系数分别取HH=2.0,HL=0.5的增强结果。
图1中三种算法都对lena图像(图1-a)进行了增强,增强效果可看出直方图均衡(图1-c)与本发明算法(图1-d)两者增强的图像从视觉上难以分辨优劣,但两者效果明显比同态滤波增强好,图像层次也比较清晰,其对应的直方图(图4-c)、(图4-d)对比度动态范围也较宽,而同态滤波增强的图像整体偏亮,图像细节丢失,不清晰,其直方图(图4-b)也显示图像灰度集中在高亮区,对比度范围较窄,低灰度值基本没有。
对于图2中加入0.003高斯噪声的lena图像(图2-a)增强,三种算法中本发明(图2-d)算法效果最好,图像清晰、对比度好、噪声抑制的也相当好;直方图均衡算法(图2-c)在增强图像的同时噪声也被放大,能明显感觉到噪声存在;同态滤波增强(图2-b)在噪声抑制上与本文算法差不多,但是在增强过程中损失了大量细节。
上面是从主观的角度对三算法的增强效果进行了定性分析,下面运用评价图像降噪与增强效果的三个性能指标来定量研究三种增强算法对噪声图像的降噪与增强情况。经过图像增强后像素灰度值会发生改变,故用原始图像经过加噪的增强后图像与未加噪的原始图像增强后图像两者进行比较,分别计算这三个性能指标。计算结果如表1所示。
表1计算结果
从表1中可以看出对比度增量:直方图均衡>本发明算法>同态滤波,说明三算法的增强效果直方图均衡效果最好,其次是本发明算法,最后是同态滤波;对于噪声图像增强的噪声抑制,从表1的峰值信噪比可以看出:本发明算法>同态滤波>直方图均衡,即本文算法最优。
所述利用二维小波分解系数的直接重构函数Y=upcoef2(’a’,g(k),’syml’,1)重构基于所述均衡化的新系数和所述高频子带分解系数cH,cV,cD的图像时,也可参考《Matlab 6.X图形图像处理》,清华大学出版社,2002年5月,ISBN 7-302-05410-X。具体参看该文献第6.4.3章节——小波分析工具箱函数介绍
dwt2二维离散小波变换函数(第205页)
[cA,cH,cV,cD]=dwt2(X,’wname’)
[cA,cH,cV,cD]=dwt2(X,Lo_D,Hi_D)
[cA,cH,cV,cD]=dwt2(X,’wname’)使用指定的小波基函数’wname’对图像X进行二维离散小波变换,cA,cH,cV,cD分别为图像分解的近似分量、水平分量、垂直分量和细节分量。
[cA,cH,cV,cD]=dwt2(X,Lo_D,Hi_D)使用指定的低通和高通滤波器组Lo_D,Hi_D来对图像进行二维离散小波变换。
upcoef2函数用于利用多层小波分解重构近似分量或者细节分量(第207页)
Y=upcoef2(O,X,’wname’,N)
Y=upcoef2(O,X,Lo_R,Hi_R,N)
Y=upcoef2(O,X,’wname’,N,S)
Y=upcoef2(O,X,Lo_R,Hi_R,N,S)
Y=upcoef2(O,X,’wname’)
Y=upcoef2(O,X,Lo_R,Hi_R)
Y=upcoef2(O,X,’wname’,N)由二维离散小波变换得到的系数重构近似分量或者细节分量,’wname’为指定的小波基函数的名称,N指定重构次数,他必须是一正整数。若不指定,则N=i.
O指定细节信号的类型:
0=’a’重构近似信号,即X是第N层的近似系数。
0=’h’重构水平细节信号,即X是第N层的水平细节系数。
0=’v’重构垂直细节信号,即X是第N层的垂直细节系数。
0=’d’重构对焦细节信号,即X是第N层的对焦细节系数。
Y=upcoef2(O,X,Lo_R,Hi_R,N)利用指定的重构滤波器组Lo_R,Hi_R,N重构第N层的近似分量或者细节分量。
显然,上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而这些属于本发明的精神所引伸出的显而易见的变化或变动仍处于本发明的保护范围之中。
Claims (4)
1.一种图像增强方法,其特征在于,包括如下步骤:
a、对原始图像进行单层离散二维小波分解,得到低频子带分解系数cA,以及分别对应于水平、垂直和对角三个方向的三个高频子带分解系数cH,cV,cD;
b、对所述低频子带系数cA做均衡化处理,得到均衡化的新系数;
c、利用所述均衡化的新系数和所述高频子带分解系数cH,cV,cD重构图像。
3.根据权利要求2所述图像增强方法,其特征在于:所述的元素k≥0。
4.根据权利要求2所述图像增强方法,其特征在于:所述步骤b4中的得到均衡化的新系数的算法是:0≤k≤Kmax,其中m,n为图像的行数与列数,i,j分别为图像矩阵的行、列变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110318806A CN102332155A (zh) | 2011-10-19 | 2011-10-19 | 一种图像增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110318806A CN102332155A (zh) | 2011-10-19 | 2011-10-19 | 一种图像增强方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102332155A true CN102332155A (zh) | 2012-01-25 |
Family
ID=45483920
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110318806A Pending CN102332155A (zh) | 2011-10-19 | 2011-10-19 | 一种图像增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102332155A (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104036455A (zh) * | 2013-03-06 | 2014-09-10 | 江苏理工学院 | 基于二代小波的红外图像细节增强方法 |
CN104639800A (zh) * | 2013-11-08 | 2015-05-20 | 华为终端有限公司 | 一种用于图像降噪的方法及终端 |
CN105550998A (zh) * | 2015-12-09 | 2016-05-04 | 江苏理工学院 | 基于二代小波整数变换的图像增强方法及图像增强系统 |
CN105608685A (zh) * | 2015-11-17 | 2016-05-25 | 江苏理工学院 | 一种直方图修正的二次直方图均衡图像增强方法及系统 |
CN106040691A (zh) * | 2016-05-24 | 2016-10-26 | 扬州市东宇环保设备有限公司 | 一种环保清洁方法 |
CN106052759A (zh) * | 2016-05-24 | 2016-10-26 | 扬州市东宇环保设备有限公司 | 一种环保监控方法 |
CN106054693A (zh) * | 2016-05-24 | 2016-10-26 | 扬州市东宇环保设备有限公司 | 移动式环保清洁平台 |
CN106052758A (zh) * | 2016-05-24 | 2016-10-26 | 扬州市东宇环保设备有限公司 | 环保监控平台 |
CN106094607A (zh) * | 2016-05-24 | 2016-11-09 | 扬州市东宇环保设备有限公司 | 移动式环保监控平台 |
CN106296600A (zh) * | 2016-08-01 | 2017-01-04 | 大连理工大学 | 一种基于小波变换图像分解的对比度增强方法 |
CN107669425A (zh) * | 2016-05-24 | 2018-02-09 | 孙霞 | 心内科手术床 |
CN107811780A (zh) * | 2016-05-24 | 2018-03-20 | 杨琳 | 儿科手术床 |
CN108269246A (zh) * | 2018-01-30 | 2018-07-10 | 江苏理工学院 | 一种低频小波系数插值的图像均衡增强方法 |
CN108961179A (zh) * | 2018-06-19 | 2018-12-07 | 上海中和软件有限公司 | 一种医学图像后处理系统及其使用方法 |
CN109427047A (zh) * | 2017-08-28 | 2019-03-05 | 京东方科技集团股份有限公司 | 一种图像处理方法及装置 |
CN112016420A (zh) * | 2020-08-19 | 2020-12-01 | 云南电网有限责任公司临沧供电局 | 一种环网柜视频图像识别系统和方法 |
CN113837975A (zh) * | 2021-09-05 | 2021-12-24 | 桂林理工大学 | 基于双三次插值和奇异值分解的图像增强方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1172312A (zh) * | 1996-05-21 | 1998-02-04 | 三星电子株式会社 | 采用低通滤波和直方图均衡化的图像增强方法及其装置 |
-
2011
- 2011-10-19 CN CN201110318806A patent/CN102332155A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1172312A (zh) * | 1996-05-21 | 1998-02-04 | 三星电子株式会社 | 采用低通滤波和直方图均衡化的图像增强方法及其装置 |
Non-Patent Citations (4)
Title |
---|
温学兵 等: "基于小波去噪和直方图模板均衡化的手指静脉图像增强", 《吉林大学学报(理学版)》 * |
牛和明 等: "基于小波信息量分段调整的图像增强", 《计算机工程》 * |
王海峰 等: "一种小波变换的小波系数均衡图像增强算法", 《电视技术》 * |
郑辉 等: "基于小波分频与直方图均衡的图像增强算法", 《现代电子技术》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104036455B (zh) * | 2013-03-06 | 2017-04-12 | 江苏理工学院 | 基于二代小波的红外图像细节增强方法 |
CN104036455A (zh) * | 2013-03-06 | 2014-09-10 | 江苏理工学院 | 基于二代小波的红外图像细节增强方法 |
CN104639800A (zh) * | 2013-11-08 | 2015-05-20 | 华为终端有限公司 | 一种用于图像降噪的方法及终端 |
US9904986B2 (en) | 2013-11-08 | 2018-02-27 | Huawei Device (Dongguan) Co., Ltd. | Image denoising method and terminal |
CN104639800B (zh) * | 2013-11-08 | 2017-11-24 | 华为终端(东莞)有限公司 | 一种用于图像降噪的方法及终端 |
CN105608685A (zh) * | 2015-11-17 | 2016-05-25 | 江苏理工学院 | 一种直方图修正的二次直方图均衡图像增强方法及系统 |
CN105608685B (zh) * | 2015-11-17 | 2018-09-28 | 江苏理工学院 | 一种直方图修正的二次直方图均衡图像增强方法及系统 |
CN105550998A (zh) * | 2015-12-09 | 2016-05-04 | 江苏理工学院 | 基于二代小波整数变换的图像增强方法及图像增强系统 |
CN105550998B (zh) * | 2015-12-09 | 2019-01-08 | 江苏理工学院 | 基于二代小波整数变换的图像增强方法及图像增强系统 |
CN107669425A (zh) * | 2016-05-24 | 2018-02-09 | 孙霞 | 心内科手术床 |
CN106052759A (zh) * | 2016-05-24 | 2016-10-26 | 扬州市东宇环保设备有限公司 | 一种环保监控方法 |
CN106094607A (zh) * | 2016-05-24 | 2016-11-09 | 扬州市东宇环保设备有限公司 | 移动式环保监控平台 |
CN106052758A (zh) * | 2016-05-24 | 2016-10-26 | 扬州市东宇环保设备有限公司 | 环保监控平台 |
CN106054693A (zh) * | 2016-05-24 | 2016-10-26 | 扬州市东宇环保设备有限公司 | 移动式环保清洁平台 |
CN107811780A (zh) * | 2016-05-24 | 2018-03-20 | 杨琳 | 儿科手术床 |
CN106040691A (zh) * | 2016-05-24 | 2016-10-26 | 扬州市东宇环保设备有限公司 | 一种环保清洁方法 |
CN106296600A (zh) * | 2016-08-01 | 2017-01-04 | 大连理工大学 | 一种基于小波变换图像分解的对比度增强方法 |
CN106296600B (zh) * | 2016-08-01 | 2018-12-21 | 大连理工大学 | 一种基于小波变换图像分解的对比度增强方法 |
CN109427047A (zh) * | 2017-08-28 | 2019-03-05 | 京东方科技集团股份有限公司 | 一种图像处理方法及装置 |
CN109427047B (zh) * | 2017-08-28 | 2021-01-26 | 京东方科技集团股份有限公司 | 一种图像处理方法及装置 |
CN108269246A (zh) * | 2018-01-30 | 2018-07-10 | 江苏理工学院 | 一种低频小波系数插值的图像均衡增强方法 |
CN108961179A (zh) * | 2018-06-19 | 2018-12-07 | 上海中和软件有限公司 | 一种医学图像后处理系统及其使用方法 |
CN108961179B (zh) * | 2018-06-19 | 2022-10-18 | 上海中和软件有限公司 | 一种医学图像后处理系统及其使用方法 |
CN112016420A (zh) * | 2020-08-19 | 2020-12-01 | 云南电网有限责任公司临沧供电局 | 一种环网柜视频图像识别系统和方法 |
CN113837975A (zh) * | 2021-09-05 | 2021-12-24 | 桂林理工大学 | 基于双三次插值和奇异值分解的图像增强方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102332155A (zh) | 一种图像增强方法 | |
CN100550978C (zh) | 一种保持边缘的自适应图像滤波方法 | |
CN102509269B (zh) | 一种结合曲线波的基于图像子块相似性的图像去噪方法 | |
CN101882304B (zh) | 一种sar图像自适应去噪和特征增强方法 | |
CN102842120B (zh) | 基于超复数小波相位测量的图像模糊程度检测方法 | |
CN105515695B (zh) | 基于调制宽带转换器的压缩采样信号检测方法 | |
CN102509263B (zh) | 基于sar图像局部统计特性的k-svd相干斑抑制方法 | |
CN105761290B (zh) | 一种自适应的多尺度分块压缩感知采样方法 | |
CN104657948B (zh) | 一种用于海洋探测的激光水下成像图像去噪和增强方法 | |
CN102572499B (zh) | 基于小波变换多分辨率预测的无参考图像质量评价方法 | |
CN105550998B (zh) | 基于二代小波整数变换的图像增强方法及图像增强系统 | |
CN104680485A (zh) | 一种基于多分辨率的图像去噪方法及装置 | |
CN102426701A (zh) | 基于双树复小波变换和pca的水下声纳图像的去噪方法 | |
CN100433795C (zh) | 基于变换域数学形态学的图像降噪方法 | |
CN106530244A (zh) | 一种图像增强方法 | |
CN104881847A (zh) | 一种基于小波分析和伪彩色处理的比赛视频图像增强方法 | |
Ismael et al. | A new approach of image denoising based on discrete wavelet transform | |
CN102226917A (zh) | 一种基于非下采样Contourlet扩散的图像增强方法 | |
CN103700077B (zh) | 一种基于人眼视觉特性的自适应图像增强方法 | |
CN104574424B (zh) | 基于多分辨率dct边缘梯度统计的无参照图像模糊度评价方法 | |
CN109724693B (zh) | 一种基于平稳小波的融合光谱去噪方法 | |
Nha et al. | Consideration of a robust watermarking algorithm for color image using improved QR decomposition | |
Bhandari et al. | SVD based poor contrast improvement of blurred multispectral remote sensing satellite images | |
Ashish et al. | Satellite image processing using discrete cosine transform and singular value decomposition | |
CN104240208A (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 | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20120125 |