CN106067164B - 基于自适应小波域处理的彩色图像对比度增强算法 - Google Patents
基于自适应小波域处理的彩色图像对比度增强算法 Download PDFInfo
- Publication number
- CN106067164B CN106067164B CN201610362381.7A CN201610362381A CN106067164B CN 106067164 B CN106067164 B CN 106067164B CN 201610362381 A CN201610362381 A CN 201610362381A CN 106067164 B CN106067164 B CN 106067164B
- Authority
- CN
- China
- Prior art keywords
- wavelet
- color image
- image
- rgb
- function
- 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.)
- Active
Links
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 15
- 230000009466 transformation Effects 0.000 claims abstract description 22
- 230000002708 enhancing effect Effects 0.000 claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000001186 cumulative effect Effects 0.000 claims description 7
- 238000005315 distribution function Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 241001269238 Data Species 0.000 claims description 4
- 239000003086 colorant Substances 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000005428 wave function Effects 0.000 claims 1
- 230000003321 amplification Effects 0.000 abstract description 8
- 238000003199 nucleic acid amplification method Methods 0.000 abstract description 8
- 230000000007 visual effect Effects 0.000 abstract description 5
- 238000000034 method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 4
- 238000011426 transformation method Methods 0.000 description 3
- 230000002401 inhibitory effect Effects 0.000 description 2
- 239000002075 main ingredient Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004438 eyesight Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Compression Of Band Width Or Redundancy In Fax (AREA)
Abstract
本发明公开了一种基于自适应小波域处理的彩色图像对比度增强算法,其包括如下步骤:步骤1、彩色图像RGB三层数据分别进行小波分解;步骤2、根据RGB三通道低频小波系数分布特性对其进行非线性变换,实现自适应增强对比度;步骤3、根据RGB三通道高频小波系数不同尺度、方向以及噪声特性,构造高斯型阈值滤波函数进行去噪,实现抑制噪声和伪像;步骤4:小波逆变换,重构RGB图像。本发明适用于彩色图像的对比度增强,使图像更加清晰、色调更加鲜明,同时抑制噪声放大、提高信噪比,有效地改善图像的视觉效果。
Description
技术领域
本发明属于图像处理技术领域,涉及一种基于自适应小波域处理的彩色图像对比度增强算法。
背景技术
由于受到光照、天候等成像条件及光学传感器退化因素的影响,实际图像可能存在着对比度下降、动态范围偏窄、有效量化不足以及目标局部细节信息分辨力不强等问题,这将影响人眼对目标精细判读解译或机器自动识别的能力。在实际应用中,通常采用图像对比度增强技术来解决上述问题,并已被广泛应用于医学图像诊断、视频监视、故障检测、计算机视觉以及光学遥感图像增强与目标识别等领域。
目前的灰度图像对比度增强算法主要分为灰度变换法、直方图调整法和反锐化掩模。其中,灰度变换法可分为对数变换、指数变换等,该类方法仅是通过调整图像灰度动态范围来提高对比度,对视觉效果提升不够明显,而且会带来噪声放大、伪像等问题。直方图调整可分为直方图均衡化、直方图规定化等,该类方法是通过重新调整图像直方图分布的方式实现对比度的增强,也存在噪声放大的问题,特别是在处理彩色图像时可能导致色彩失真。反锐化掩模仅仅增强图像边缘能量,即增强高频信息,同样会导致噪声放大,而且计算繁杂,处理效率低。
发明内容
本发明的目的是针对现有彩色图像对比度增强算法存在的图像增强效果不佳、适应性不足、易出现噪声放大以及计算复杂等问题,提出一种基于自适应小波域处理的彩色图像对比度增强算法,适用于自适应地增强彩色图像对比度,使图像更加清晰、色调鲜明,并能够有效地抑制噪声放大,增强图像视觉效果。
本发明的目的是通过以下技术方案实现的:
一种基于自适应小波域处理的彩色图像对比度增强算法,包括如下步骤:
步骤1、RGB彩色图像三通道分别进行小波变换:
(1)将彩色图像动态范围调整到[0,1]区间:
式中,Imin为待处理彩色图像,Amax和Amin分别为彩色图像矩阵中最大值和最小值;
(2)RGB三通道分别进行小波分解:
式中,ψ为sym4小波基函数,j为小波分解尺度,R、G、B分别代表三基色通道,WR、WG和WB分别为RGB三通道小波系数,x、y和p、q分别表示图像域和小波域的坐标,Z为整数集合。
步骤2、根据RGB三通道低频小波系数分布特性对其进行非线性变换,实现自适应增强对比度:
(1)提取三通道低频小波系数,构成矩阵WA,并计算矩阵WA累积分布函数:
式中,m为WA中最小值,Nr(ω)表示数值等于ω的小波系数数量,N为小波系数长度;
在累积分布函数中找到T值为1%和99%时相应的小波系数ωm和ωM,即:
T(ωm)=1%,
T(ωM)=99%;
(2)根据低频小波系数整体分布特性,求解非线性变换函数f(ω):
f(ω)=aω2+bω+c;
式中,a、b、c为待定系数,利用该变换函数将ωm和ωM分别变换到WA中最小值m和最大值M处,同时WA均值ωz变换到2j-1附近,即:
解方程可确定待定系数:
(3)利用f(ω)对WA进行灰度变换得WAx,然后对WAx中含有的小于m或大于M的数据进行剔除:
式中,WAout为低频小波系数最终增强结果;
(4)从WAout中分离出三层数据,替换三通道低频小波系数。
步骤3、根据RGB三通道高频小波系数不同尺度、方向以及噪声特性,构造高斯型阈值滤波函数进行去噪,实现抑制噪声和伪像:
(1)根据RGB三通道高频小波系数的不同尺度、方向以及噪声特性,计算相应滤波阈值:
式中,Lj和σx分别为j尺度子带图像小波系数的长度和标准差,为噪声方差;
(2)构造高斯型阈值函数F(ω),对各通道、尺度、方向的高频小波系数分别进行阈值滤波,达到抑制噪声和伪像的目的:
步骤4、通过小波逆变换重构RGB图像:
对经过步骤2和步骤3处理后得到的小波系数进行逆变换,重新构成RGB图像,即得到最终增强结果。
本发明相比于已有算法,具有以下优点:
(1)本发明提出的基于自适应小波域处理的彩色图像对比度增强算法既能够自适应地增强图像对比度,又能够有效地保持彩色图像的色调并抑制噪声和伪像,避免色彩失真从而使得彩色图像的色彩更加鲜明,显著增强图像的视觉效果。同时,该算法可应用于灰度图像的对比度增强处理。
(2)本发明基于灰度变换法,根据RGB三通道低频小波系数整体分布特性选取变换函数,通过对三通道低频小波系数全局化非线性变换,自适应调整小波系数分布,从而实现了图像的自动对比度增强,同时有效保证了图像色调不变、亮度适中。而且该方法具有算法简单、适应性强、处理效率高等优点。
(3)本发明通过对RGB三通道各尺度、方向的高频小波系数进行不同阈值滤波,达到抑制噪声放大的目的。同时,针对硬阈值滤波易产生振铃和软阈值滤波易出现能量损失严重、重构误差大等问题,构建了高斯型阈值滤波函数,通过该函数能够实现抑制噪声的同时,保持景物纹理细节信息的目的,并且能够有效地抑制振铃并降低重构误差,从而使得算法具有较强的适应性。
附图说明
图1为基于小波变换的彩色图像自适应对比度增强算法流程;
图2为图像累积分布函数示意图;
图3为非线性灰度变换函数曲线示意图;
图4为软阈值函数和硬阈值函数示意图;
图5为高斯型阈值函数曲线示意图;
图6为噪声图像小波分解结果;
图7为小波阈值滤波结果;
图8为彩色图像(一);
图9为彩色图像(一)处理结果;
图10为彩色图像(二);
图11为彩色图像(二)处理结果;
图12为灰度图像;
图13为灰度图像处理结果。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
本发明提供了一种基于自适应小波域处理的彩色图像对比度增强算法,用于彩色图像对比度增强,提高视觉效果。以RGB图像为例,如图1所示,其具体实施步骤如下:
步骤1、RGB彩色图像三通道分别进行小波分解:
先将彩色图像动态范围调整到[0,1]区间:
式中,Imin为待处理彩色图像,Amax和Amin分别为彩色图像矩阵中最大值和最小值。
然后,RGB三通道分别进行小波分解:
式中,ψ为sym4小波基函数,j为小波分解尺度,R、G、B分别代表三基色通道,WR、WG和WB分别为RGB三通道小波系数,x、y和p、q分别表示图像域和小波域的坐标,Z为整数集合。
步骤2、根据RGB三通道低频小波系数分布特性对其进行非线性变换,实现自适应增强对比度:
为保证图像色调不变,对三层数据按相同法则进行处理。
(1)分别提取WR、WG和WB的低频部分,构成矩阵WA,并计算矩阵WA累积分布函数(CDF):
式中,m为WA中最小值,Nr(ω)表示数值等于ω的小波系数的数量,N为小波系数长度。
(2)如图2所示,在累积分布函数中找到T值为1%和99%时(此数值可根据实际要求适当调整)相应的小波系数ωm和ωM,即:
T(ωm)=1%,
T(ωM)=99%。
(3)图像的对比度是指图像中最亮部分与最暗部分的密度之比,图像纹理的沟纹越深、脊线越突出,对比度越高。因此,将ωm和ωM分别向两侧拉伸即可实现对比度增强;另外,为将小波系数分布密集区域宽度扩展、分布稀疏区域宽度压缩,同时保证图像亮度适中,如图3所示,本发明选取非线性变换函数为:
f(ω)=aω2+bω+c,
式中,a、b、c为待定系数。利用该变换函数将ωm和ωM分别变换到WA中最小值m和最大值M处,同时WA均值ωz变换到2j-1附近,即:
解方程可确定待定系数:
(4)利用f(ω)对WA进行灰度变换得WAx,然后对WAx中含有的小于m或大于M的数据进行剔除:
式中,WAout为低频小波系数最终增强结果。
(5)从WAout中分离出三层数据,替换WR、WG和WB的低频小波系数。
步骤3:根据RGB三通道高频小波系数不同尺度、方向以及噪声特性,构造高斯型阈值滤波函数进行去噪,实现抑制噪声和伪像。
图像和噪声在经小波变换后具有不同的统计特性,图像本身的能量对应着幅值较大的小波系数,而噪声能量则对应着幅值较小的小波系数。根据这一特征,设置一个阈值门限Th,对WR、WG和WB的高频子带进行滤波。认为大于该阈值的小波系数的主要成份为有用信号,给予保留;小于该阈值的小波系数,主要成份为噪声,予以剔除,达到降噪目的。
在小波变换中,不同尺度的高频小波系数差别很大,因此阈值需要根据不同尺度来进行相应的调整。
式中,Lj和σx分别为j尺度子带图像小波系数的长度和标准差,为噪声方差。
常用的阈值滤波函数有软阈值函数和硬阈值函数(如图4所示),但是硬阈值函数具有阶跃特性,容易引起振铃现象;软阈值函数虽然连续性好,但会造成一定信息损失,给重构带来误差。实际应用中,常常需要对阈值函数进行改进。
此处目的在于抑制噪声放大,而非去除噪声,针对该需求,本发明采用高斯函数构造阈值滤波函数(如图5所示),将阈值处小波系数幅值削弱1/e,幅值很小的小波系数大幅度削弱,而幅值较大的小波系数几乎完全保留,滤波公式如下:
利用该滤波函数分别对WR、WG和WB中高频小波系数进行滤波。如图6、图7所示,降噪效果良好。
步骤4:小波逆变换,重构RGB图像。
对经过步骤2和步骤3处理后得到的小波系数进行逆变换,重新构成RGB图像,即得到最终增强结果。
处理结果展示如图8~图13所示。
Claims (1)
1.一种基于自适应小波域处理的彩色图像对比度增强算法,其特征在于所述算法步骤如下:
步骤1、彩色图像RGB三层数据分别进行小波分解;
所述步骤1的具体步骤如下:
(1)将彩色图像动态范围调整到[0,1]区间:
式中,Imin为待处理彩色图像,Amax和Amin分别为彩色图像矩阵中最大值和最小值;
(2)RGB三通道分别进行小波分解:
式中,ψ为sym4小波基函数,j为小波分解尺度,R、G、B分别代表三基色通道,WR、WG和WB分别为RGB三通道小波系数,x、y和p、q分别表示图像域和小波域的坐标,Z为整数集合;
步骤2、根据RGB三通道低频小波系数分布特性对其进行非线性变换,实现自适应增强对比度;
所述步骤2的具体步骤如下:
(1)提取三通道低频小波系数,构成矩阵WA,并计算矩阵WA累积分布函数:
式中,m为WA中最小值,Nr(ω)表示数值等于ω的小波系数数量,N为小波系数长度;
在累积分布函数中找到T值为1%和99%时相应的小波系数ωm和ωM,即:
T(ωm)=1%,
T(ωM)=99%;
(2)根据低频小波系数整体分布特性,求解非线性变换函数f(ω):
f(ω)=aω2+bω+c;
式中,a、b、c为待定系数,利用该变换函数将ωm和ωM分别变换到WA中最小值m和最大值M处,同时WA均值ωz变换到2j-1附近,即:
解方程可确定待定系数:
(3)利用f(ω)对WA进行灰度变换得WAx,然后对WAx中含有的小于m或大于M的数据进行剔除:
式中,WAout为低频小波系数最终增强结果;
(4)从WAout中分离出三层数据,替换三通道低频小波系数;
步骤3、根据RGB三通道高频小波系数不同尺度、方向以及噪声特性,构造高斯型阈值滤波函数进行去噪,实现抑制噪声和伪像;
所述步骤3的具体步骤如下:
(1)根据RGB三通道高频小波系数的不同尺度、方向以及噪声特性,计算相应滤波阈值:
式中,Lj和σx分别为j尺度子带图像小波系数的长度和标准差,为噪声方差;
(2)构造高斯型阈值函数F(ω),对各通道、尺度、方向的高频小波系数分别进行阈值滤波,达到抑制噪声和伪像的目的:
步骤4:小波逆变换,重构RGB图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362381.7A CN106067164B (zh) | 2016-05-26 | 2016-05-26 | 基于自适应小波域处理的彩色图像对比度增强算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362381.7A CN106067164B (zh) | 2016-05-26 | 2016-05-26 | 基于自适应小波域处理的彩色图像对比度增强算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106067164A CN106067164A (zh) | 2016-11-02 |
CN106067164B true CN106067164B (zh) | 2018-11-20 |
Family
ID=57420085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610362381.7A Active CN106067164B (zh) | 2016-05-26 | 2016-05-26 | 基于自适应小波域处理的彩色图像对比度增强算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106067164B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106556940B (zh) * | 2016-11-10 | 2019-11-19 | 武汉精测电子集团股份有限公司 | 一种tft-lcd屏自动光学检测中的背景抑制方法 |
CN106709891A (zh) * | 2016-11-15 | 2017-05-24 | 哈尔滨理工大学 | 基于小波变换与自适应变换结合的图像处理方法 |
CN107861367A (zh) * | 2017-12-26 | 2018-03-30 | 南通尚力机电工程设备有限公司 | 一种钟表性能在线诊断系统 |
CN108961179B (zh) * | 2018-06-19 | 2022-10-18 | 上海中和软件有限公司 | 一种医学图像后处理系统及其使用方法 |
CN110365914B (zh) * | 2019-07-24 | 2020-03-20 | 中国人民解放军国防科技大学 | 一种图像动态范围展宽方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103310423A (zh) * | 2013-07-08 | 2013-09-18 | 中国矿业大学(北京) | 一种矿井图像增强方法 |
CN103440623A (zh) * | 2013-08-02 | 2013-12-11 | 中北大学 | 基于成像模型提高雾天图像清晰度的方法 |
CN104657948A (zh) * | 2015-02-09 | 2015-05-27 | 重庆大学 | 一种用于海洋探测的激光水下成像图像去噪和增强方法 |
CN104915936A (zh) * | 2015-06-23 | 2015-09-16 | 香港应用科技研究院有限公司 | 基于小波变换的图像去色和增强 |
CN105303532A (zh) * | 2015-10-21 | 2016-02-03 | 北京工业大学 | 一种小波域Retinex图像去雾方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5668105B2 (ja) * | 2013-06-25 | 2015-02-12 | アキュートロジック株式会社 | 画像処理装置、画像処理方法及び画像処理プログラム |
-
2016
- 2016-05-26 CN CN201610362381.7A patent/CN106067164B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103310423A (zh) * | 2013-07-08 | 2013-09-18 | 中国矿业大学(北京) | 一种矿井图像增强方法 |
CN103440623A (zh) * | 2013-08-02 | 2013-12-11 | 中北大学 | 基于成像模型提高雾天图像清晰度的方法 |
CN104657948A (zh) * | 2015-02-09 | 2015-05-27 | 重庆大学 | 一种用于海洋探测的激光水下成像图像去噪和增强方法 |
CN104915936A (zh) * | 2015-06-23 | 2015-09-16 | 香港应用科技研究院有限公司 | 基于小波变换的图像去色和增强 |
CN105303532A (zh) * | 2015-10-21 | 2016-02-03 | 北京工业大学 | 一种小波域Retinex图像去雾方法 |
Non-Patent Citations (1)
Title |
---|
基于Shearlet变换和多尺度Retinex的遥感图像增强算法;王静静;《计算机应用》;20150110;第35卷(第1期);第202-205页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106067164A (zh) | 2016-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106067164B (zh) | 基于自适应小波域处理的彩色图像对比度增强算法 | |
CN105303532B (zh) | 一种小波域Retinex图像去雾方法 | |
US11127122B2 (en) | Image enhancement method and system | |
Bhonsle et al. | Medical image denoising using bilateral filter | |
Suganya et al. | Survey on image enhancement techniques | |
CN107862666A (zh) | 基于nsst域的混合遥感图像增强方法 | |
CN108961172B (zh) | 一种基于Gamma校正的图像对比度增强方法 | |
CN103295204A (zh) | 一种基于非下采样轮廓波变换的图像自适应增强方法 | |
CN106570843A (zh) | 一种自适应小波阈值函数图像噪声抑制方法 | |
Zhan et al. | Infrared image enhancement based on wavelet transformation and retinex | |
CN109064413B (zh) | 图像对比度增强方法及采用其的图像采集医疗设备 | |
CN104766285A (zh) | 一种水下降质图像的自适应增强方法 | |
CN107784639B (zh) | 一种无人机遥感图像改进多边滤波去噪方法 | |
Hongqiao et al. | A new image denoising method using wavelet transform | |
Mohapatra et al. | Histogram equalization and noise removal process for enhancement of image | |
Wu et al. | Low-Light image enhancement algorithm based on HSI color space | |
Jagatheeswari et al. | A novel approach for contrast enhancement based on histogram equalization followed by median filter | |
CN113850744A (zh) | 基于自适应Retinex和小波融合的图像增强算法 | |
CN116012273B (zh) | 一种基于局部灰度波动率的图像增强方法和装置 | |
Kusnik et al. | On the robust technique of mixed Gaussian and impulsive noise reduction in color digital images | |
Panchaxri et al. | Image denoising using adaptive NL means filtering with method noise thresholding | |
Chang et al. | An image enhancement algorithm based on Gaussian weighted bilateral filtering and retinex theory | |
Subah et al. | Underwater image enhancement based on fusion technique via color correction and illumination adjustment | |
Qingrong et al. | Contrast enhancement of medical images using fuzzy set theory and nonsubsampled shearlet transform | |
CN113487505A (zh) | 一种基于噪声识别的红外图像混合降噪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |