CN110852977B - 融合边缘灰度直方图与人眼视觉感知特性的图像增强方法 - Google Patents
融合边缘灰度直方图与人眼视觉感知特性的图像增强方法 Download PDFInfo
- Publication number
- CN110852977B CN110852977B CN201911036351.7A CN201911036351A CN110852977B CN 110852977 B CN110852977 B CN 110852977B CN 201911036351 A CN201911036351 A CN 201911036351A CN 110852977 B CN110852977 B CN 110852977B
- Authority
- CN
- China
- Prior art keywords
- gradient
- gray level
- pixel
- image
- edge
- 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
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000016776 visual perception Effects 0.000 title claims abstract description 16
- 230000035945 sensitivity Effects 0.000 claims abstract description 14
- 230000008447 perception Effects 0.000 claims abstract description 11
- 230000001186 cumulative effect Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000010606 normalization Methods 0.000 claims abstract description 5
- 230000004438 eyesight Effects 0.000 claims description 7
- 230000000007 visual effect Effects 0.000 abstract description 6
- 230000000694 effects Effects 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000010191 image analysis Methods 0.000 description 3
- 241001270131 Agaricus moelleri Species 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000010827 pathological analysis Methods 0.000 description 1
- 238000011160 research Methods 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/40—Image enhancement or restoration using histogram techniques
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
一种融合边缘灰度直方图与人眼视觉感知特性的图像增强方法：分别计算图像中每一个像素位置的行方向的梯度和列方向的梯度，以及该像素位置总的梯度；将所有像素位置总的梯度进行排序，找出其中最大的总的梯度，并设定梯度阈值；将所有像素位置总的梯度与梯度阈值进行比较，标记出像素位置总的梯度大于梯度阈值的所有图像边缘像素，并统计标记出的所有图像边缘像素中灰度级相同的边缘像素的数量，得到边缘灰度直方图；计算边缘灰度级小于设定值的累积概率分布；考虑人眼视觉非线性特性，计算人眼在不同亮度背景下的感知敏感度因子，并做归一化处理，得到最佳的灰度级调整参考值；从而得到增强图像。本发明显著提升图像的视觉效果。
Description
技术领域
本发明涉及一种图像增强方法。特别是涉及一种融合边缘灰度直方图与人眼视觉感知特性的图像增强方法。
背景技术
在现实生活中，由于受拍照环境光照不佳和成像设备性能不足等因素的影响，所获得的图像易出现细节模糊、颜色失真、对比度低、噪声较多等现象，这些降质图像对于图像分析与理解带来极大的困难。例如医学影像分析上，出于对患者安全检查的考虑，加之医疗设备图像信息采集能力的限制，所获取的各类医学影像的清晰度和对比度并不理想，直接影响医生病理诊断的准确性。显然，借助图像增强技术对大量降质图像进行增强处理具有重要的实际意义。一方面可以显著改善图像的视觉效果。另一方面增强的结果更有利于后续的图像分析与辨识处理。
在现有的图像增强算法中，直方图均衡方法(Global Histogram Equalization，GHE)是最经典且应用最广泛的增强算法，该算法将灰度累积概率分布作为灰度的调整函数，算法设计简单，易于编程实现，对于大多数降质图像的增强有一定的效果。后续很多研究人员针对GHE算法的灰度调整策略存在的问题做了一系列改进。改进的亮度保留双直方图均衡(BPBHE)算法首先根据输入图像亮度的均值将输入图像的直方图分解为两个子图，然后对这两个子直方图进行独立均衡。为了克服GHE这种全局算法对局部区域细节增强的不足，提出了改进的自适应直方图均衡(AHE)算法，该算法依据局部区域直方图计算灰度的映射关系，因此该算法的计算复杂度较大。另外，也有一些科研人员尝试将直方图均衡思想推广到频率域，提出了多尺度自适应直方图均衡(MAHE)，该算法首先使用小波变换将输入图像分解成子带，然后用自适应直方图均衡处理高频子带，最后将处理后的高频子带与低频子带合并得到增强图像。
尽管改进后的算法图像增强效果有所改善，但这些算法仍存在两方面的共性问题，影响了算法的增强性能。首先是分析数据存在的问题，GHE及其系列改进算法的分析数据都是图像的直方图，图像直方图仅仅是图像中各灰度简单的数量统计结果，并未真实反映各灰度对信息刻画的作用大小，因此依据灰度数量的概率分布进行灰度调整不能保证对信息刻画重要的灰度级动态范围进行有效拉伸。其次，人眼视觉感知研究发现，人眼对不同亮度背景的信息敏感程度存在明显差异，对亮度背景较暗的灰度偏差感知能力较弱，对亮度背景适中或较亮环境的灰度偏差感知能力较强。而GHE及其系列改进算法都是基于人眼不同亮度背景敏感度相同的假设，使得大量暗背景下的细节信息不能得到有效增强。
发明内容
本发明所要解决的技术问题是，提供一种能够显著改善图像的视觉增强效果的融合边缘灰度直方图与人眼视觉感知特性的图像增强方法。
本发明所采用的技术方案是：一种融合边缘灰度直方图与人眼视觉感知特性的图像增强方法，包括如下步骤：
1)分别计算图像中每一个像素I(x,y)位置的行方向的梯度和列方向的梯度，以及该像素位置总的梯度；
2)将所有像素位置总的梯度G进行排序，找出其中最大的总的梯度Gmax，并将0.15*Gmax设定为梯度阈值；
3)将所有像素位置总的梯度与梯度阈值进行比较，标记出像素位置总的梯度大于梯度阈值的所有图像边缘像素，并统计标记出的所有图像边缘像素中灰度级相同的边缘像素的数量，得到边缘灰度直方图；
4)依据边缘灰度直方图，使用如下公式计算边缘灰度级小于或等于k的累积概率分布Sk：
其中，nj为第j个边缘灰度级的像素个数，n为标记的所有像素的数量总和；
5)考虑人眼视觉非线性特性，计算人眼在不同亮度背景下的感知敏感度因子Em；
6)对人眼在不同亮度背景下的感知敏感度因子Em做归一化处理，得到最佳的灰度级调整参考值Am
7)找到与边缘灰度级小于或等于k的累积概率分布Sk最接近的灰度级调整参考值Am，将原始图像中像素灰度级数值为k的灰度级变换为灰度级m，从而得到增强图像。
步骤1)所述的：
行方向的梯度是采用如下公式得到：
列方向的梯度是采用如下公式得到：
其中，Gr为像素的行方向梯度，Gc为像素的列方向梯度，A0～A7是以像素I(x,y)的8邻域左上角的像素为起点顺时针遍历8个邻近像素得到的对应的8个灰度值。
步骤1)所述的总的梯度G是采用如下公式得到：
其中，Gr为像素的行方向梯度，Gc为像素的列方向梯度。
步骤5)是利用如下公式得到人眼在不同亮度背景下的感知敏感度因子Em
其中，J1是灰度级m＝1时的最小可见偏差Jm的值，最小可见偏差Jm由如下公式计算得到：
其中，Bm为灰度级为m的像素的背景区域的平均亮度值。
本发明的融合边缘灰度直方图与人眼视觉感知特性的图像增强方法，将边缘的灰度直接参与了图像信息的刻画，对边缘的灰度进行统计能更准确反映各灰度对图像信息刻画的重要性。另外，考虑到人眼对不同亮度背景的灰度偏差敏感的非线性特性，提出了新的灰度调整策略，确保在亮度适中或较亮的背景区域展示更多的图像细节信息，从而显著改善图像的视觉增强效果。
本发明相比经典的简单灰度级数量统计直方图，边缘灰度直方图的统计结果能更准确地反映不同灰度级对于信息刻画的重要程度，确保后续灰度级动态范围调整时，能对信息刻画重要的灰度级进行有效地拉伸，同时很好地抑制背景噪声。
本发明充分考虑了人眼视觉感知的非线性特点，在灰度级调整过程中将更多图像细节信息调整到人眼视觉比较敏感的较亮的灰度背景区域，显著提升图像的视觉效果。
本发明的方法可广泛应用于不同领域的影像处理，尤其在降质的医学影像处理方面，应用本发明算法可以更清晰地展示影像的细节信息且对噪声有很好的抑制效果，增强后的影像可以帮助医生做出更准确的诊断。
附图说明
图1是本发明融合边缘灰度直方图与人眼视觉感知特性的图像增强方法的流程图；
图2是本发明图像像素8邻域示意图；
图3是人眼视觉感知亮度阈值效应曲线图(横轴表示某一灰度级的像素的平均背景亮度值，纵轴表示最小可见偏差)；
图4a是原始工件图像；
图4b是GHE算法增强后的工件图像；
图4c是CLAHE算法增强后的工件图像；
图4d是本发明算法增强后的工件图像；
图5a是原始X光影像；
图5b是GHE算法增强后的X光影像；
图5c是CLAHE算法增强后的X光影像；
图5d是本发明算法增强后的X光影像；
图6a是原始自然场景图像；
图6b是GHE算法增强后的自然场景图像；
图6c是CLAHE算法增强后的自然场景图像；
图6d是本发明算法增强后的自然场景图像。
具体实施方式
下面结合实施例和附图对本发明的融合边缘灰度直方图与人眼视觉感知特性的图像增强方法做出详细说明。
由于图像的大部分主要的细节信息分布在图像的边缘位置，本发明的融合边缘灰度直方图与人眼视觉感知特性的图像增强方法，首先是计算各个像素位置总的梯度，然后找到所有像素位置总的梯度的最大值Gmax，设定0.15*Gmax为梯度阈值，然后标记出总的梯度大于梯度阈值的所有图像边缘像素，并统计标记出的所有边缘像素中灰度级相同像素的数量，得到边缘灰度直方图。由于考虑人眼视觉的非线性的影响，本发明的融合边缘灰度直方图与人眼视觉感知特性的图像增强方法，还要计算人眼在不同亮度背景下的感知敏感度因子Em，然后利用归一化之后的人眼视觉感知敏感度因子优化调整灰度级的动态范围。因而本发明的方法主要分为两个主要的部分，第一部分是计算得到图像的边缘灰度直方图，第二部分是利用人眼视觉感知特性进行灰度级动态范围的调整。
如图1所示，本发明的融合边缘灰度直方图与人眼视觉感知特性的图像增强方法，包括如下步骤：
1)分别计算图像中每一个像素I(x,y)位置的行方向的梯度和列方向的梯度，以及该像素位置总的梯度；其中，
所述的行方向的梯度是采用如下公式得到：
所述的列方向的梯度是采用如下公式得到：
其中，Gr为像素的行方向梯度，Gc为像素的列方向梯度，A0～A7是以像素I(x,y)的8邻域左上角的像素为起点顺时针遍历8个邻近像素得到的对应的8个灰度值，如图2所示。
所述的总的梯度G是采用如下公式得到：
其中，Gr为像素的行方向梯度，Gc为像素的列方向梯度。
2)将所有像素位置总的梯度G进行排序，找出其中最大的总的梯度Gmax，并将0.15*Gmax设定为梯度阈值；
3)将所有像素位置总的梯度与梯度阈值进行比较，标记出像素位置总的梯度大于梯度阈值的所有图像边缘像素，并统计标记出的所有图像边缘像素中灰度级相同的边缘像素的数量，得到边缘灰度直方图；
4)依据边缘灰度直方图，使用如下公式计算边缘灰度级小于或等于k的累积概率分布Sk：
其中，nj为第j个边缘灰度级的像素个数，n为标记的所有像素的数量总和；
5)考虑人眼视觉非线性特性，计算人眼在不同亮度背景下的感知敏感度因子，具体是利用如下公式得到人眼在不同亮度背景下的感知敏感度因子Em
其中，J1是灰度级m＝1时的最小可见偏差Jm的值，最小可见偏差Jm曲线如图3所示，并由如下公式计算得到：
其中，Bm为灰度级为m的像素的背景区域的平均亮度值。
6)对人眼在不同亮度背景下的感知敏感度因子Em做归一化处理，得到最佳的灰度级调整参考值Am
7)找到与边缘灰度级小于或等于k的累积概率分布Sk最接近的灰度级调整参考值Am，将原始图像中像素灰度级数值为k的灰度级变换为灰度级m，从而得到增强图像。
根据说明书附图，图4a、图5a与图6a是原始图像，图4b、图5b与图6b是GHE进行算法增强后的图像，图4c、图5c与图6c是CLAHE算法增强后的图像，图4d、图5d与图6d是本发明算法增强后的图像，通过说明书附图的图像增强算法的增强效果对比，可以看出，本发明算法可以更清晰地展示影像的细节信息且对噪声有很好的抑制效果。
Claims (1)
1.一种融合边缘灰度直方图与人眼视觉感知特性的图像增强方法，其特征在于，包括如下步骤：
1)分别计算图像中每一个像素I(x,y)位置的行方向的梯度和列方向的梯度，以及该像素位置总的梯度；其中所述的：
行方向的梯度是采用如下公式得到：
列方向的梯度是采用如下公式得到：
其中，Gr为像素的行方向梯度，Gc为像素的列方向梯度，A0～A7是以像素I(x,y)的8邻域左上角的像素为起点顺时针遍历8个邻近像素得到的对应的8个灰度值；
所述的总的梯度G是采用如下公式得到：
其中，Gr为像素的行方向梯度，Gc为像素的列方向梯度；
2)将所有像素位置总的梯度G进行排序，找出其中最大的总的梯度Gmax，并将0.15*Gmax设定为梯度阈值；
3)将所有像素位置总的梯度与梯度阈值进行比较，标记出像素位置总的梯度大于梯度阈值的所有图像边缘像素，并统计标记出的所有图像边缘像素中灰度级相同的边缘像素的数量，得到边缘灰度直方图；
4)依据边缘灰度直方图，使用如下公式计算边缘灰度级小于或等于k的累积概率分布Sk：
其中，nj为第j个边缘灰度级的像素个数，n为标记的所有像素的数量总和；
5)考虑人眼视觉非线性特性，计算人眼在不同亮度背景下的感知敏感度因子Em；是利用如下公式得到人眼在不同亮度背景下的感知敏感度因子Em
其中，J1是灰度级m＝1时的最小可见偏差Jm的值，最小可见偏差Jm由如下公式计算得到：
其中，Bm为灰度级为m的像素的背景区域的平均亮度值；
6)对人眼在不同亮度背景下的感知敏感度因子Em做归一化处理，得到最佳的灰度级调整参考值Am
7)找到与边缘灰度级小于或等于k的累积概率分布Sk最接近的灰度级调整参考值Am，将原始图像中像素灰度级数值为k的灰度级变换为灰度级m，从而得到增强图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911036351.7A CN110852977B (zh) | 2019-10-29 | 2019-10-29 | 融合边缘灰度直方图与人眼视觉感知特性的图像增强方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911036351.7A CN110852977B (zh) | 2019-10-29 | 2019-10-29 | 融合边缘灰度直方图与人眼视觉感知特性的图像增强方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110852977A CN110852977A (zh) | 2020-02-28 |
| CN110852977B true CN110852977B (zh) | 2023-04-11 |
Family
ID=69598164
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911036351.7A Expired - Fee Related CN110852977B (zh) | 2019-10-29 | 2019-10-29 | 融合边缘灰度直方图与人眼视觉感知特性的图像增强方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110852977B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114693682B (zh) * | 2022-06-01 | 2022-08-26 | 青岛大学附属医院 | 一种基于图像处理的脊椎特征识别方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2339797A1 (en) * | 1998-08-07 | 2000-03-02 | Arch Development Corporation | Method and system for the segmentation and classification of lesions |
| CN105608676A (zh) * | 2015-12-23 | 2016-05-25 | 浙江宇视科技有限公司 | 一种视频图像的增强方法及装置 |
| CN109801233A (zh) * | 2018-12-27 | 2019-05-24 | 中国科学院西安光学精密机械研究所 | 一种适用于真彩色遥感图像的增强方法 |
| WO2019144581A1 (zh) * | 2018-01-29 | 2019-08-01 | 江苏宇特光电科技股份有限公司 | 一种智能化红外图像场景增强方法 |
-
2019
- 2019-10-29 CN CN201911036351.7A patent/CN110852977B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2339797A1 (en) * | 1998-08-07 | 2000-03-02 | Arch Development Corporation | Method and system for the segmentation and classification of lesions |
| CN105608676A (zh) * | 2015-12-23 | 2016-05-25 | 浙江宇视科技有限公司 | 一种视频图像的增强方法及装置 |
| WO2019144581A1 (zh) * | 2018-01-29 | 2019-08-01 | 江苏宇特光电科技股份有限公司 | 一种智能化红外图像场景增强方法 |
| CN109801233A (zh) * | 2018-12-27 | 2019-05-24 | 中国科学院西安光学精密机械研究所 | 一种适用于真彩色遥感图像的增强方法 |
Non-Patent Citations (3)
| Title |
|---|
| 曾明等.视觉注意机制在图像增强中的应用研究.光子学报.2009,第38卷(第5期),1283-1287. * |
| 范晓鹏等.基于人眼亮度阈值特性的图像增强算法.光电子·激光.2014,第25卷(第8期),1606-1612. * |
| 赵瑞宇等.高铁闸机智能监控系统中的行人检测算法研究.计算机应用与软件.2013,第30卷(第3期),85-88. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110852977A (zh) | 2020-02-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Agrawal et al. | A novel joint histogram equalization based image contrast enhancement | |
| CN110505459B (zh) | 适用于内窥镜的图像颜色校正方法、装置和存储介质 | |
| Wang et al. | Simple low-light image enhancement based on Weber–Fechner law in logarithmic space | |
| Sun et al. | Fast single image haze removal via local atmospheric light veil estimation | |
| CN116205823B (zh) | 一种基于空域滤波的超声影像去噪方法 | |
| Acharya et al. | Particle swarm optimized texture based histogram equalization (PSOTHE) for MRI brain image enhancement | |
| CN110570360B (zh) | 一种基于Retinex的鲁棒和全面的低质量光照图像增强方法 | |
| CN116883537B (zh) | 基于图像增强的普通外科手术影像采集方法 | |
| CN108257099B (zh) | 一种基于视觉对比度分辨率的自适应红外图像增强方法 | |
| CN114298916B (zh) | 一种基于灰度拉伸和局部增强的X-Ray图像增强方法 | |
| CN117853386B (zh) | 一种肿瘤图像增强方法 | |
| Kumar et al. | Luminosity and contrast enhancement of retinal vessel images using weighted average histogram | |
| US20090074317A1 (en) | System and method for reducing halo effect in image enhancement | |
| CN114529475A (zh) | 基于二维伽马校正和色调映射的图像增强方法及系统 | |
| CN116245880A (zh) | 基于红外识别的电动车充电桩火灾风险检测方法 | |
| CN108447031A (zh) | 一种红外图像热背景噪声校正的方法 | |
| CN115578284A (zh) | 一种多场景图像增强方法及系统 | |
| CN114972067A (zh) | 一种x光小牙片图像增强方法 | |
| CN114359083B (zh) | 一种面向干扰环境的高动态热红外图像自适应预处理方法 | |
| CN110852977B (zh) | 融合边缘灰度直方图与人眼视觉感知特性的图像增强方法 | |
| Zhang et al. | A fundus image enhancer based on illumination-guided attention and optic disc perception GAN | |
| CN116630198A (zh) | 一种结合自适应伽马校正的多尺度融合水下图像增强方法 | |
| CN116452502A (zh) | 医学图像多阶段增强方法及系统 | |
| Weichao et al. | Research on color image defogging algorithm based on MSR and CLAHE | |
| Meng et al. | Improved multi-scale retinex algorithm for medical image enhancement |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20230411 |