CN111292267A - 一种基于拉普拉斯金字塔的图像主观视觉效果增强方法 - Google Patents
一种基于拉普拉斯金字塔的图像主观视觉效果增强方法 Download PDFInfo
- Publication number
- CN111292267A CN111292267A CN202010080143.3A CN202010080143A CN111292267A CN 111292267 A CN111292267 A CN 111292267A CN 202010080143 A CN202010080143 A CN 202010080143A CN 111292267 A CN111292267 A CN 111292267A
- Authority
- CN
- China
- Prior art keywords
- image
- sub
- histogram
- block
- mean
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000000007 visual effect Effects 0.000 title claims abstract description 19
- 230000003044 adaptive effect Effects 0.000 claims abstract description 15
- 230000002708 enhancing effect Effects 0.000 claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 12
- 238000013507 mapping Methods 0.000 claims description 12
- 230000000694 effects Effects 0.000 abstract description 7
- 230000003321 amplification Effects 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 1
- 230000002844 continuous effect Effects 0.000 description 1
Images
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
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- 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/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- 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/20021—Dividing image into blocks, subimages or windows
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种基于拉普拉斯金字塔的图像主观视觉效果增强方法,属于数字图像处理领域。该方法包括:一、对原始图像进行拉普拉斯金字塔分解得到近似图像和多尺度细节图像。二、基于改进对比度受限自适应直方图均衡化算法对近似图像进行增强。三、对多尺度细节图像进行增强;四、对增强后的近似图和多尺度细节图进行拉普拉斯逆变换得到增强后图像。通过以上步骤,可以有效地缓解传统对比度受限自适应直方图均衡化增强后图像中可能出现的边界效应、过增强和噪声过大等问题,显著提升图像主观视觉效果。
Description
技术领域
本发明提供一种基于拉普拉斯金字塔的图像主观视觉效果增强 方法,属于图像处理技术领域。
背景技术
图像增强是图像处理的重要处理步骤,它可以有效地改善图像 质量,提高图像主观视觉效果。其中对比度受限自适应直方图均衡 化是最为常用的局部图像增强算法之一,它首先把图像分为多个连 续不重叠的图像块,然后再分别进行直方图修正和均衡化得到对应 的灰度映射函数,最后进行双线性插值缓解不同图像块之间的不连 续效应。该算法可以有效地增强局部图像的细节信息,提高图像主 观视觉效果,但其存在如下问题:(1)图像块效益明显,如果相 邻图像块之间的映射函数差别比较大,经过双线性插值以后,其边界块效应依然会比较明显;(2)图像过增强问题,如果子图像存 在大面积灰度比较接近的区域,经过直方图修正后,其峰值依然明 显,容易造成过增强;(3)图像放大噪声问题,由于图像噪声和 细节都属于图像高频信息,对比度受限自适应直方图均衡化方法在 增强图像局部细节的时候会放大图像的噪声。
发明内容
本发明的目的是提供一种基于拉普拉斯金字塔的图像主观视觉 效果增强方法,它能够缓解传统对比度受限自适应直方图均衡化方法 增强后图像中可能出现的边界效应、过增强和噪声过大等问题,显著 提升图像主观视觉效果。
为实现上述目的,本发明的基于拉普拉斯金字塔的图像主观视 觉效果增强方法,其步骤如下:
步骤S1:图像拉普拉斯金字塔分解;
所述的拉普拉斯金字塔分解步骤中,是对输入图像进行拉普拉斯 金字塔分解得到近似图像LpK和多尺度细节图像Lpk(k=1,2,...,K-1); 然后将多尺度细节图像缩放至近似图像一样大小,求出缩放后图像均 值得到多尺度细节均值图像其中, K为拉普拉斯金字塔层数,Lp_resizedk(k=1,2,...,K-1)为缩放后多 尺度细节图像。
步骤S2:近似图像增强;
所述的近似图像增强步骤中,采用改进对比度受限自适应直方图 均衡化算法对步骤S1中得到的近似图像进行增强;其中,该步骤包 括:
步骤S21,将步骤S1中获取的图像LpK分割为 round(M/s)×round(N/s)个大小为s×s的连续不重叠子块,记为 其中M×N为LpK图像尺寸,round表示取整操作。 子块在图像Lp_mean中对应(位置一样)的子块记 为Lp_meani(i=1,2,...,s2)。对图像中的每个像素点 (x,y)(x=1,2,...,M,y=1,2,...,N)进行遍历,在像素点灰度值对应的直方图 个数加上log(Lp_meani(x,y)):
从步骤S21中子块图像直方图中的定义可以看出,其和像素的 数目和其对应的细节信息相关,而且包含细节信息越少的像素点对直 方图的影响越小。当图像中存在大量灰度级接近的背景区域时,其在 传统灰度直方图中会形成一个峰值,进而导致图像过增强;由于灰度 级接近的区域其包含的细节信息较小,其对我们定义的直方图影响较 小,可以有效避免峰值。
步骤S24,采用传统对比度受限自适应直方图均衡化增强算法中 的双线性插值方法,对步骤S23中获取的各子块图像映射函数进行插 值得到图像中每个像素的映射灰度级,并进行灰度替换,完成近似图 像增强,得到增强后近似图像
步骤S3:多尺度细节图像增强;
所述的多尺度细节图像增强步骤中,对步骤S1中得到的各尺度 细节图像Lpk(k=1,2,...,K-1)进行放大,得到增强后细节图像。
其中为增强后细节图像,mean(|Lpk|)为 细节图像Lpk绝对值均值,a(1<a<3)为增强系数。对多尺度细节图 像进行放大,可以增强图像细节、锐化图像。通过设定自适应阈值, 可以在一定程度上抑制噪声的放大。
步骤S4:拉普拉斯逆变换
与现有技术相比本发明的技术效果在于:传统对比度受限自适 应图像均衡化算法统计每个子块图像的灰度直方图,容易出现直方 图峰值,导致过增强。虽然其对直方图进行了修正,但当峰值比较 明显时,过增强问题依然存在。本发明提出的图像直方图是对图像 进行拉普拉斯分解,然后基于多尺度细节图像包含的细节信息来统 计子块图像直方图。基于该直方图,本发明可以有效避免传统对比 度受限自适应图像均衡化算法中可能出现的图像过增强现象。
尤其,传统对比度受限自适应图像均衡化算法可能出现边界效 应,当相邻子块图像直方图差异较大时,子块图像边界的不连续现 象尤为明显。本发明所提出的子块图像直方图加权算法,该算法对 子块图像直方图和所有子块图像均值直方图进行加权,自适应降低 子块图像直方图之间的差异,有效避免对比度受限自适应图像均衡 化中可能出现的边界效应。
尤其,传统对比度受限自适应图像均衡化算法对图像局部细节 增强的同时,容易放大图像的噪声。本发明是对图像进行拉普拉斯 分解,得到多尺度高频细节图和低频近似图像。图像噪声属于高频 信息,基本都包含在细节图中,只需要对低频近似图像进行增强,可有效地避免图像噪声过大。
附图说明
图1为本发明的流程图;
图2为本发明步骤S3的流程图
图3为本发明的一个实施例中处理前的原始图像;
图4为四个子块图像的灰度直方图和本发明的加权后细节直 方图;
图5为采用传统对比度受限自适应直方图均衡化算法增强后 图像;
图6为采用本发明的方法增强后。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚明白,以下结 合具体实施例,并参照附图,对本发明作进一步的详细说明。
本发明所述的一种基于拉普拉斯金字塔的图像主观视觉效果增 强方法,其流程图如图1所示,包括一下步骤:
步骤S1:图像拉普拉斯金字塔分解
采用现有技术对输入图像进行拉普拉斯金字塔分解得到近似图 像LpK和多尺度细节图像Lpk(k=1,2,...,K-1);然后将多尺度细节图 像缩放至近似图像一样大小,求出缩放后图像均值得到多尺度细节均 值图像其中,K为拉普拉斯金字塔 层数,Lp_resizedk(k=1,2,...,K-1)为缩放后多尺度细节图像。
步骤S2:近似图像增强;对步骤S1中获取的近似图像LpK进行 增强,其包括如下步骤:
步骤S21,将步骤S1获取的图像LpK分割为 round(M/s)×round(N/s)个大小为s×s的连续不重叠子块,记为 其中M×N为LpK图像尺寸,round表示取整操作。 子块在图像Lp_mean中对应(位置一样)的子块记 为Lp_meani(i=1,2,...,s2)。对图像中的每个像素点 (x,y)(x=1,2,...,M,y=1,2,...,N)进行遍历,在像素点灰度值对应的直方图 个数加上log(Lp_meani(x,y)):
步骤S24,采用传统对比度受限自适应直方图均衡化增强算法中 的双线性插值方法,对步骤S23中获取的各子块图像映射函数进行插 值得到图像中每个像素的映射灰度级,并进行灰度替换,完成近似图 像增强,得到增强后近似图像
步骤S3:多尺度细节图像增强
对步骤S1中得到的各尺度细节图像Lpk(k=1,2,...,K-1)进行如 下放大:
步骤S4:拉普拉斯逆变换
为验证本专利提出算法的有效性、合理性、可行性及科学性,对 图3中原始图像采用传统对比度受限自适应直方图均衡化和本发明的 方法进行增强。拉普拉斯分解层数K=2,子块图像数目为4,细节增 强系数a=1.5。图4是原始图像四个子块图像的灰度直方图和对应的 加权后细节直方图;图5是传统对比度受限自适应直方图均衡化算法 增强后图像;图6是应用本发明的方法增强后图像。从图4中可以看出, 原始图像中大面积的灰色背景区域灰度比较接近,因而在传统灰度直 方图中形成了个峰值。我们的直方图由像素点个数和其包含的细节信 息共同决定。由于图像中的背景区域比较平缓,包含的细节信息较小, 其对直方图的影响较小。因此我们的细节直方图中没有产生峰值。另 一方面,通过直方图加权后,子块图像之间的直方图差异明显变小。 从图5中可以看出,由于传统灰度直方图存在峰值,导致直方图均衡 化出现过增强现象;另外,由于子块直方图之间差异比较明显,子图像边界之间存在少量不连续现象。从图6中可以看出,我们算法能够 有效避免图像过增强、边界效应和噪声的放大,大幅提升图像主观视 觉效果。
Claims (7)
1.一种基于拉普拉斯金字塔的图像主观视觉效果增强方法,其特征在于,包括以下步骤:
步骤S1:图像拉普拉斯金字塔分解;
所述的拉普拉斯金字塔分解步骤中,是对输入图像进行拉普拉斯金字塔分解得到近似图像和多尺度细节图像;然后将多尺度细节图像缩放至近似图像一样大小,求出缩放后图像均值得到多尺度细节均值图像;
步骤S2:近似图像增强;
所述的近似图像增强步骤中,采用改进对比度受限自适应直方图均衡化算法对步骤S1中得到的近似图像进行增强;其中,该步骤包括:
步骤S21,将所述近似图像分割为连续不重叠子块图像,利用步骤S1中得到多尺度细节均值图像的信息,统计各子块图像直方图,并求出各子块直方图均值。
步骤S22,把所述各子块图像直方图和直方图均值进行加权操作得到各子块图像加权后直方图;
步骤S23,将所述各子块加权后直方图进行累加和归一化操作,得到各子块映射函数;
步骤S24,采用传统对比度受限自适应直方图均衡化增强算法中的双线性插值方法,对步骤S23中获取的各子块图像映射函数进行插值得到图像中每个像素的映射灰度级,并进行灰度替换,完成近似图像增强[1]。
步骤S3:多尺度细节图像增强;
所述的多尺度细节图像增强步骤中,对步骤S1中得到的各尺度细节图像进行放大,得到增强后细节图像。
步骤S4:拉普拉斯逆变换
所述拉普拉斯逆变换步骤中,是对步骤S2中增强后的近似图像和步骤S3中增强后的各尺度细节图像进行拉普拉斯逆变换得到增强后图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010080143.3A CN111292267B (zh) | 2020-02-04 | 2020-02-04 | 一种基于拉普拉斯金字塔的图像主观视觉效果增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010080143.3A CN111292267B (zh) | 2020-02-04 | 2020-02-04 | 一种基于拉普拉斯金字塔的图像主观视觉效果增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111292267A true CN111292267A (zh) | 2020-06-16 |
CN111292267B CN111292267B (zh) | 2020-10-23 |
Family
ID=71026722
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010080143.3A Active CN111292267B (zh) | 2020-02-04 | 2020-02-04 | 一种基于拉普拉斯金字塔的图像主观视觉效果增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111292267B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113129244A (zh) * | 2021-04-16 | 2021-07-16 | 重庆邮电大学 | 一种基于局部拉普拉斯分解的多尺度医学图像融合方法 |
CN114565535A (zh) * | 2022-03-02 | 2022-05-31 | 北京锐影医疗技术有限公司 | 一种基于自适应梯度gamma校正的图像增强方法及装置 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101030298A (zh) * | 2007-03-29 | 2007-09-05 | 杭州电子科技大学 | 多尺度自适应对比度变换的医学图像增强方法 |
US20100129003A1 (en) * | 2008-11-21 | 2010-05-27 | Pace Teresa L | Method and system for enhancing images using multi-resolution histogram shaping |
CN102413283A (zh) * | 2011-10-25 | 2012-04-11 | 广州飒特红外股份有限公司 | 红外热图数字信号处理系统及方法 |
CN102439629A (zh) * | 2010-06-08 | 2012-05-02 | 伊姆普斯封闭式股份有限公司 | 数字图像校正方法 |
US20120147225A1 (en) * | 2008-08-01 | 2012-06-14 | Nikon Corporation | Image processing method |
CN102637292A (zh) * | 2011-02-10 | 2012-08-15 | 西门子公司 | 一种图像的处理方法和装置 |
US8411938B2 (en) * | 2007-11-29 | 2013-04-02 | Sri International | Multi-scale multi-camera adaptive fusion with contrast normalization |
CN103426176A (zh) * | 2013-08-27 | 2013-12-04 | 重庆邮电大学 | 基于改进直方图和聚类算法的视频镜头检测方法 |
CN103700075A (zh) * | 2013-12-25 | 2014-04-02 | 浙江师范大学 | 基于Tetrolet变换的多通道卫星云图融合方法 |
CN103903227A (zh) * | 2012-12-29 | 2014-07-02 | 上海联影医疗科技有限公司 | 一种图像降噪的方法和装置 |
CN104574337A (zh) * | 2015-01-21 | 2015-04-29 | 山东科技大学 | 基于双边伽马校正和多尺度图像融合的图像增强方法 |
CN105678330A (zh) * | 2016-01-05 | 2016-06-15 | 北京环境特性研究所 | 一种基于高斯加权的直方图匹配方法 |
CN106023057A (zh) * | 2016-05-26 | 2016-10-12 | 深圳大学 | 一种用于皮下静脉显影仪的控制处理系统及成像方法 |
CN106910159A (zh) * | 2016-10-12 | 2017-06-30 | 王玲 | 视频拼接方法及装置 |
CN110175964A (zh) * | 2019-05-30 | 2019-08-27 | 大连海事大学 | 一种基于拉普拉斯金字塔的Retinex图像增强方法 |
-
2020
- 2020-02-04 CN CN202010080143.3A patent/CN111292267B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101030298A (zh) * | 2007-03-29 | 2007-09-05 | 杭州电子科技大学 | 多尺度自适应对比度变换的医学图像增强方法 |
US8411938B2 (en) * | 2007-11-29 | 2013-04-02 | Sri International | Multi-scale multi-camera adaptive fusion with contrast normalization |
US20120147225A1 (en) * | 2008-08-01 | 2012-06-14 | Nikon Corporation | Image processing method |
US20100129003A1 (en) * | 2008-11-21 | 2010-05-27 | Pace Teresa L | Method and system for enhancing images using multi-resolution histogram shaping |
CN102439629A (zh) * | 2010-06-08 | 2012-05-02 | 伊姆普斯封闭式股份有限公司 | 数字图像校正方法 |
CN102637292A (zh) * | 2011-02-10 | 2012-08-15 | 西门子公司 | 一种图像的处理方法和装置 |
CN102413283A (zh) * | 2011-10-25 | 2012-04-11 | 广州飒特红外股份有限公司 | 红外热图数字信号处理系统及方法 |
CN103903227A (zh) * | 2012-12-29 | 2014-07-02 | 上海联影医疗科技有限公司 | 一种图像降噪的方法和装置 |
CN103426176A (zh) * | 2013-08-27 | 2013-12-04 | 重庆邮电大学 | 基于改进直方图和聚类算法的视频镜头检测方法 |
CN103700075A (zh) * | 2013-12-25 | 2014-04-02 | 浙江师范大学 | 基于Tetrolet变换的多通道卫星云图融合方法 |
CN104574337A (zh) * | 2015-01-21 | 2015-04-29 | 山东科技大学 | 基于双边伽马校正和多尺度图像融合的图像增强方法 |
CN105678330A (zh) * | 2016-01-05 | 2016-06-15 | 北京环境特性研究所 | 一种基于高斯加权的直方图匹配方法 |
CN106023057A (zh) * | 2016-05-26 | 2016-10-12 | 深圳大学 | 一种用于皮下静脉显影仪的控制处理系统及成像方法 |
CN106910159A (zh) * | 2016-10-12 | 2017-06-30 | 王玲 | 视频拼接方法及装置 |
CN110175964A (zh) * | 2019-05-30 | 2019-08-27 | 大连海事大学 | 一种基于拉普拉斯金字塔的Retinex图像增强方法 |
Non-Patent Citations (1)
Title |
---|
刘轩等: ""基于对比度受限自适应直方图均衡的乳腺图像增强"", 《计算机工程和应用》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113129244A (zh) * | 2021-04-16 | 2021-07-16 | 重庆邮电大学 | 一种基于局部拉普拉斯分解的多尺度医学图像融合方法 |
CN113129244B (zh) * | 2021-04-16 | 2022-10-18 | 重庆邮电大学 | 一种基于局部拉普拉斯分解的多尺度医学图像融合方法 |
CN114565535A (zh) * | 2022-03-02 | 2022-05-31 | 北京锐影医疗技术有限公司 | 一种基于自适应梯度gamma校正的图像增强方法及装置 |
CN114565535B (zh) * | 2022-03-02 | 2022-09-13 | 北京锐影医疗技术有限公司 | 一种基于自适应梯度gamma校正的图像增强方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111292267B (zh) | 2020-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107527333B (zh) | 一种基于伽马变换的快速图像增强方法 | |
CN110827229B (zh) | 一种基于纹理加权直方图均衡化的红外图像增强方法 | |
Yun et al. | Image enhancement using a fusion framework of histogram equalization and Laplacian pyramid | |
JP4502303B2 (ja) | 画像処理装置 | |
EP1909227B1 (en) | Method of and apparatus for minimizing ringing artifacts in an input image | |
CN111292267B (zh) | 一种基于拉普拉斯金字塔的图像主观视觉效果增强方法 | |
CN107292834B (zh) | 红外图像细节增强方法 | |
CN108280836B (zh) | 一种图像处理方法及装置 | |
CN108961172B (zh) | 一种基于Gamma校正的图像对比度增强方法 | |
CN112911166B (zh) | 一种调整图像亮度的方法、装置、芯片、介质及摄像设备 | |
Hossain et al. | Medical image enhancement based on nonlinear technique and logarithmic transform coefficient histogram matching | |
CN104463814A (zh) | 基于局部纹理方向性的图像增强方法 | |
CN111724315B (zh) | 一种基于自适应加权中值滤波的红外图像噪声去除方法 | |
CN104715453A (zh) | 一种分区域处理的图像增强方法及电路 | |
CN112541869A (zh) | 一种基于matlab实现的Retinex图像去雾方法 | |
CN115965544A (zh) | 自适应亮度调节的图像增强方法及系统 | |
CN114565535B (zh) | 一种基于自适应梯度gamma校正的图像增强方法及装置 | |
CN116012273A (zh) | 一种基于局部灰度波动率的图像增强方法和装置 | |
CN111325685A (zh) | 一种基于多尺度相对梯度直方图均衡化的图像增强算法 | |
CN109242797B (zh) | 基于均质和异质区域融合的图像去噪方法、系统及介质 | |
Ramya et al. | A novel method for the contrast enhancement of fog degraded video sequences | |
CN116862809A (zh) | 一种低曝光条件下的图像增强方法 | |
CN109801246B (zh) | 一种自适应阈值的全局直方图均衡方法 | |
CN111325680A (zh) | 一种具有抑制局部过曝的图像增亮方法 | |
CN115511737A (zh) | 基于语义分割的红外图像增强方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |