CN106530237B - 一种图像增强方法 - Google Patents
一种图像增强方法 Download PDFInfo
- Publication number
- CN106530237B CN106530237B CN201610833349.2A CN201610833349A CN106530237B CN 106530237 B CN106530237 B CN 106530237B CN 201610833349 A CN201610833349 A CN 201610833349A CN 106530237 B CN106530237 B CN 106530237B
- Authority
- CN
- China
- Prior art keywords
- image
- histogram
- edge
- denoising
- carried out
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000002708 enhancing effect Effects 0.000 claims abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 11
- 230000007797 corrosion Effects 0.000 claims description 13
- 238000005260 corrosion Methods 0.000 claims description 13
- 230000011218 segmentation Effects 0.000 claims description 11
- RSWGJHLUYNHPMX-UHFFFAOYSA-N Abietic-Saeure Natural products C12CCC(C(C)C)=CC2=CCC2C1(C)CCCC2(C)C(O)=O RSWGJHLUYNHPMX-UHFFFAOYSA-N 0.000 claims description 8
- KHPCPRHQVVSZAH-HUOMCSJISA-N Rosin Natural products O(C/C=C/c1ccccc1)[C@H]1[C@H](O)[C@@H](O)[C@@H](O)[C@@H](CO)O1 KHPCPRHQVVSZAH-HUOMCSJISA-N 0.000 claims description 8
- KHPCPRHQVVSZAH-UHFFFAOYSA-N trans-cinnamyl beta-D-glucopyranoside Natural products OC1C(O)C(O)C(CO)OC1OCC=CC1=CC=CC=C1 KHPCPRHQVVSZAH-UHFFFAOYSA-N 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 4
- 230000003628 erosive effect Effects 0.000 claims description 4
- 230000000877 morphologic effect Effects 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000000873 masking effect Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000000007 visual effect Effects 0.000 description 4
- 241000023320 Luma <angiosperm> Species 0.000 description 3
- OSWPMRLSEDHDFF-UHFFFAOYSA-N methyl salicylate Chemical compound COC(=O)C1=CC=CC=C1O OSWPMRLSEDHDFF-UHFFFAOYSA-N 0.000 description 3
- 241000251468 Actinopterygii Species 0.000 description 2
- 102100032404 Cholinesterase Human genes 0.000 description 2
- 101000943274 Homo sapiens Cholinesterase Proteins 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- YLCXGBZIZBEVPZ-UHFFFAOYSA-N Medazepam Chemical compound C12=CC(Cl)=CC=C2N(C)CCN=C1C1=CC=CC=C1 YLCXGBZIZBEVPZ-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 239000013039 cover film Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003707 image sharpening Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification 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/73—Deblurring; Sharpening
- G06T5/75—Unsharp masking
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
-
- 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/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- 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/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- 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)
- Facsimile Image Signal Circuits (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及图像处理领域，更具体地，涉及一种图像增强方法。其具体步骤包括：a)对输入图像进行去噪处理得到去噪图像；b)对去噪图像进行边缘提取得到边缘图像；c)对边缘图像进行图像增强处理得到去噪且边缘增强的图像；d)利用亮度可控的直方图均衡方法对去噪图像进行处理得到全局增强图像；e)对c和d步骤所得到的图像进行线性叠加，得到最终的输出图像。本发明通过亮度可控的直方图均衡方法结合UM(Unsharp Masking，反锐化掩膜)算法思想，可以实现输出亮度可以跟进用户需求自动调节，而且通过设定合适的亮度值可以得到一幅对比度明显提升全局增强的输出图像，从而达到图像增强的目的。
Description
技术领域
本发明涉及图像处理领域，更具体地，涉及一种图像增强方法。
背景技术
图像增强是数字图像处理的一项最基本的技术，也是许多图像处理的预处理技术，其基本思想是：通过采用一系列技术去改善图像的质量和视觉效果，突出图像中感兴趣的特征，获得图像中有价值的信息，从而将图像转换成一种更适合于人或机器进行分析和处理的形式，使得处理后的图像对某些特定的应用有更好的效果。图像增强理论广泛应用于生物医学领域、工业生产领域、公共安全领域以及航空航天领域等。现有的图像增强方法很多，最基本的图像增强有空间域增强和频率域增强，但是图像增强理论的应用一般都是具有针对性的，对不同的应用运用不同的方法，其通用的图像增强方法是不存在的。
图像增强技术的关键在于如何有效地改善增强图像质量和增强视觉效果的同时更好地保留图像的边缘和细节信息，其中，与本发明相对比用到的图像增强方法，有如下两类：
(1)直方图均衡：该方法包括全局直方图均衡、直方图规定化和局部直方图均衡，这些方法通常用来增加许多图像的全局对比度，尤其是当图像的有用数据的对比度相当接近的时候。通过这种方法，亮度可以更好地在直方图上分布，从而可以用于增强局部的对比度而不影响整体的对比度，直方图均衡化通过有效地扩展常用的亮度来实现这种功能。但是，这些方法会增加背景噪声的对比度并且降低有用信号的对比度，也可能导致过增强与细节丢失问题。
(2)反锐化掩模：在印刷和出版界已用了多年的图像锐化处理过程是从原图像中减去一幅非锐化(平滑过的)版本，即将图像的模糊部分从原图像减去而得到清晰的图像，此即为反锐化掩模算法的思路，其算法处理步骤如下：(1)模糊原图像；(2)从原图像中减去模糊图像(产生的差值图像称为模板)；(3)将模板线性加到原图像上。该方法能够提升图像的高频信息，增强图像轮廓，但同时也可能增强噪声和出现振铃效应。
发明内容
本发明为克服上述现有技术所述的至少一种缺陷(不足)，提供一种能有效增强图像对比度，还能抑制噪声和保持图像细节的图像增强方法。
为解决上述技术问题，本发明的技术方案如下：
一种图像增强方法，包括如下步骤：
a)对输入图像进行去噪处理得到去噪图像；
b)对去噪图像进行边缘提取得到边缘图像；
c)对边缘图像进行图像增强处理得到去噪且边缘增强的图像；
d)利用亮度可控的直方图均衡方法对去噪图像进行处理得到全局增强图像；
e)对c和d步骤所得到的图像进行线性叠加，得到最终的输出图像。
步骤a的具体步骤包括：
a1)将输入图像转变为灰度图像；
a2)对灰度图像进行均值滤波得到平滑后的去噪图像。
步骤b的具体步骤为利用自定义的拉普拉斯模板对平滑滤波后的图像进行边缘提取得到边缘图像。
步骤c的具体步骤包括：
c1)对边缘图像进行预处理，得到边缘预处理图像，对边缘预处理图像进行图像增强得到边缘增强图像；
c2)对边缘图像进行二值化处理得到二值化图像，对二值化图像进行腐蚀处理得到腐蚀图像；
c3)将c1的边缘增强图像和c2的腐蚀图像进行综合处理得到去噪且边缘增强的图像。
步骤c1的具体步骤包括：
c11)求输入图像I(x,y)的最大和最小灰度值分别为Imax和Imin；
c12)对边缘图像进行预处理，得到边缘预处理图像Ew(x,y)，即
其中，E(x,y)为边缘图像，I(x,y)为输入图像；
c13)求出边缘预处理图像Ew(x,y)的灰度最大值max、灰度最小值min以及亮度平均值μ0、标准差σ0；
c14)求出边缘预处理图像的直方图，然后求出该直方图灰度值小于0的区域的阈值T1和直方图灰度值大于0的区域的阈值T2；
c15)根据c13和c14中求得的灰度最大值max、灰度最小值min以及图像的阈值T1和T2，将边缘预处理图像Ew(x,y)的直方图分成3个区域(min，T1)、(T1，T2)和(T2，max)，然后对该直方图进行基于均值和标准差剪切的分段直方图均衡，得到边缘增强图像Ee(x,y)。
步骤c14中阈值是通过Rosin方法求取的。
步骤c2和c3的具体步骤包括:
c21)根据Rosin方法求得边缘图像的阈值Tt，然后运用二值化求得二值图像B(x,y)；
c22)对二值图像进行形态学腐蚀得到腐蚀图像R(x,y)；
c31)结合步骤c和c2，得到去噪且边缘增强的图像Ewe(x,y)：
步骤d的具体步骤包括：
d1)求出去噪图像F(x,y)的直方图，并求得其亮度最大值fmax和最小值fmin，其中，0≤[fmin，fmax]≤255以及平均亮度值μ和标准差σ；
d2)根据平均亮度值μ和标准差σ，求出去噪图像F(x,y)直方图的分割点th1和th2为
其中，w为权值，可调节分割点的大小，一般取w＝1，0≤[th1，th2]≤255；
d3)根据fmax、fmin、th1、th2将去噪图像F(x,y)的直方图分成低、中、高三段，如下
其中，h(i)为图像F(x,y)的直方图统计函数，i表示0到255的灰度值。
d4)对直方图进行分段剪切与补偿，得到经过裁剪和补偿后的直方图；
d5)对图像F(x,y)的直方图进行分段剪切与补偿后，各段直方图的像素总数所占图像F(x,y)的总像素的比例均未变化，即
其中，toal代表图像F(x,y)的总像素；r1、r2、r3分别表示去噪图像F(x,y)直方图的低、中、高区域的像素总数占图像F(x,y)的像素总数的比例；
d6)计算直方图分段后的低、中、高区域的累积密度函数分别为
其中，Sl、Sm、Su分别为低中高直方图区域内的像素总数，hl″、hm″、hu″为分段剪切与补偿后各个区域的直方图统计函数。
d7)假设全局增强的输出图像G(x，y)的直方图的分割点分别为th1′和th2′，亮度平均值为μm，标准差为σm，根据输出图像直方图平均亮度估算模型和输出图像直方图标准差估算模型可分别估算出亮度平均值μm标准差σm，即
μm＝0.5[r1(th1′-1)+r2(th1′+th2′-1)+r3(th2′+255)]
其中，th2′＝th1′+2σm；
d8)令模型求出的亮度平均值和设定的亮度平均值mυ相等，即σm＝mυ，其中，mυ是用户可自行设定的平均亮度值；d7步骤中中的三个方程组是关于th1 1、th2′和σm三个未知数的方程组，可通过迭代的方式，计算出输出图像的直方图灰度分割点th1′和th2′；
d9)定义输出图像的动态范围[0，255]的图像映射曲线函数T为：
其中，th1和th2为自定义的去噪图像F(x,y)的直方图的两个分割点，th1′和th2′为输出图像直方图平均亮度估算模型和标准差估算模型计算出来的输出图像的直方图灰度级分割点；
d9)根据以上可得，从而求出全局增强的输出图像G(x,y)为
G(x，y)＝T(F(x，y))。
所述步骤e包括:
e1)将c和d得到的图像进行线性叠加，得到最终的输出图像O(x,y)，
O(x,y)＝G(x,y)+λ×Ewe(x,y)
其中，λ为比例因子。
与现有技术相比，本发明技术方案的有益效果是：
(1)本发明通过亮度可控的直方图均衡方法结合UM(Unsharp Masking，反锐化掩膜)算法思想，可以实现输出亮度可以跟进用户需求自动调节，而且通过设定合适的亮度值可以得到一幅对比度明显提升全局增强的输出图像，从而达到图像增强的目的。
(2)本发明通过对边缘图像进行一系列的处理，如预处理，分段直方图均衡、二值化和腐蚀操作，从而得到细节保持的边缘增强图像。
(3)本发明最后还通过结合UM算法思想将图像进行线性叠加，获得一幅信息流丰富、对比度提升、动态范围较高，且适合用户主观视觉效果的输出图像。
附图说明
图1为本发明图像增强方法的流程图.
图2为本发明图像增强方法中边缘增强的流程图。
图3为本发明图像增强方法中BCHE方法流程图。
图4为本发明输出图像直方图平均亮度估算模型。
图5为本发明输出图像直方图平均标准差估算模型。
图6为本发明采样的原始输入图像butterfly。
图7为通过本发明方法对图6图像处理后得到的输出图像，其初始亮度m＝80。
图8为通过本发明方法对图6图像处理后得到的输出图像，其初始亮度m＝110。
图9为本发明采样的原始输入图像fish。
图10为通过本发明方法对图9图像处理后得到的输出图像，其初始亮度m＝80。
图11为通过本发明方法对图9图像处理后得到的输出图像，其初始亮度m＝110。
图12为本发明的原始图像和输出图像对比度的对比。
图13为本发明的原始图像和输出图像熵的对。
具体实施方式
附图仅用于示例性说明，不能理解为对本专利的限制；
为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；
对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。
实施例1
如图1所示，本发明一种图像增强方法具体实施例的具体步骤报包括：
a)对输入图像I(x,y)进行去噪处理得到去噪图像F(x,y)；在此实施例中，去噪处理采用平滑滤波实现，首先将输入图像I(x,y)转变为灰度图像，然后对灰度图像进行均值滤波得到平滑后的去噪图像F(x,y)。
b)对去噪图像F(x,y)进行边缘提取得到边缘图像E(x,y)；在此实施例中，采用自定义拉普拉斯模板对平滑滤波后的图像进行边缘提取，例如定义拉普拉斯模板为w,
根据该自定义的拉普拉斯模板，求出边缘图像E(x,y)。
c)对边缘图像E(x,y)进行图像增强处理得到去噪且边缘增强的图像；
d)利用亮度可控的直方图均衡方法(Brightness Controllable HistogramEqualization，BCHE)对去噪图像进行处理得到全局增强图像；
e)对c和d步骤所得到的图像进行线性叠加，得到最终的输出图像。
在具体实施过程中，首先，对输入图像I(x,y)通过平滑滤波进行平滑去噪得到去噪图像F(x,y)；其次，采用自定义的拉普拉斯模板对图像进行边缘提取得到边缘图像E(x,y)；然后，分两步对边缘图像E(x,y)进行处理：
第一步：边缘增强，先对边缘图像E(x,y)进行预处理得到边缘预处理图像Ew(x,y)，再采用分段直方图均衡对边缘预处理图像进行增强，从而得到边缘增强的图像Ee(x,y)；
第二步：先利用Rosin算法寻找边缘图像的阈值，再对其进行二值化处理，得到去噪的边缘图像B(x,y)，之后再进行形态学腐蚀处理，消除双边缘效应，得到腐蚀图像R(x,y)；最后，对前面两步进行综合处理，得到提升后的具有细节保持和去噪的边缘增强图像Ewe(x,y)。
具体地，如图2所示，上述所述的第一步的具体步骤为：
c1)对边缘图像进行预处理，得到边缘预处理图像，对边缘预处理图像进行图像增强得到边缘增强图像；具体地：
c11)求输入图像I(x,y)的最大和最小灰度值分别为Imax和Imin；
c12)对边缘图像进行预处理，得到边缘预处理图像Ew(x,y)，即
其中，E(x,y)为边缘图像，I(x,y)为输入图像；
c13)求出边缘预处理图像Ew(x,y)的灰度最大值max、灰度最小值min以及亮度平均值μ0、标准差σ0；
c14)求出边缘预处理图像的直方图，其直方图形状类似高斯分布函数，其灰度值主要集中分布于原点及两侧附近，由此，可根据Rosin方法求出其直方图的分段阈值T1(直方图灰度值小于0区域的阈值)和T2(直方图灰度值大于0区域的阈值)，其中Rosin方法主要步骤如下:
c141)找到直方图的峰值点和谷值点；
c142)求出连接峰值点和谷值点的直线y＝Ax+B；
c143)在峰值点和谷值点区间的直方图曲线上找到一点到直线y＝Ax+B的垂直距离最大的点(x',y')，即可求得对应的阈值T＝x'。
c15)根据c13和c14中求得的灰度最大值max、灰度最小值min以及图像的阈值T1和T2，将边缘预处理图像Ew(x,y)的直方图分成3个区域(min，T1)、(T1，T2)和(T2，max)，然后对该直方图进行基于均值和标准差剪切的分段直方图均衡，得到边缘增强图像Ee(x,y)。
具体地，如图2所示，上述所述的第二步的具体步骤为：
c2)对边缘图像进行二值化处理得到二值化图像，对二值化图像进行腐蚀处理得到腐蚀图像；具体为：
c21)根据Rosin方法求得边缘图像的阈值Tt，然后运用二值化求得二值图像B(x,y)；
c22)对二值图像B(x,y)进行形态学腐蚀得到腐蚀图像R(x,y)；
c3)最后，将c1的边缘增强图像Ee(x,y)和c2的腐蚀图像R(x,y)进行综合处理得到去噪且边缘增强的图像Ewe(x,y)，
在具体实施过程中，步骤d中，先求平滑滤波后的去噪图像F(x,y)的亮度最大值、亮度最小值、亮度平均值和标准差，根据这4个值将去噪图像F(x,y)的直方图分成3段，并进行直方图剪切和补偿；接着，根据直方图亮度估算模型和标准差估算模型预测图像亮度与标准差，再通过计算求出灰度级分割点，从而估算图像的平均亮度值和标准差；然后，求出相对误差，再设定一个预设误差值，如果相对误差大于预设误差，则令初始设定的标准差等于估算的标准差，并重新求解灰度分割点，获得估算标准差，直到相对误差小于预设误差，从而求得最终的灰度分割点；最后，根据以上所求，运用直方图剪切和分段直方图均衡方法求得全局增强的图像。
如图3、图4和图5所示，步骤d的具体步骤为：
d1)求出去噪图像F(x,y)的直方图，并求得其亮度最大值fmax和最小值fmin，，其中，0≤[fmin，fmax]≤255以及平均亮度值μ和标准差σ，
其中，去噪图像F(x,y)h(i)为直方图统计函数，i表示0到255的灰度值，M、N为图像的行和列，M×N为图像总像素；
d2)根据平均亮度值μ和标准差σ，求出去噪图像F(x,y)直方图的分割点th1和th2为，
其中，w为权值，可调节分割点的大小，一般取w＝1，0≤[th1，th2]≤255；从而可将去噪图像F(x,y)的直方图分成低、中、高三段，分别为hl、hm和hu，
d3)对直方图进行分段剪切与补偿，得到经过裁剪和补偿后的直方图，步骤如下：
d31)定义r1、r2和r3分别为各段直方图在整体直方图中的比例，
d32)对第一段直方图hl进行裁剪，定义裁剪阈值Tl，
定义裁剪后的直方图为hl＇，
为了不改变hl在整体的比例，需要将裁剪多出的部分线性补偿到直方图中，定义补偿后的直方图为hl″，
其中，resl为第一段直方图经过裁剪出来的数量之和，i＝0，1，...，th1-1；
d33)对第二段直方图hm进行裁剪，定义裁剪阈值Tm，
定义裁剪后的直方图为hm′，
为了不改变hm在整体的比例，需要将裁剪多出的部分均匀补偿到直方图中，定义补偿后的直方图为hm″，
其中，resm为第二段直方图经过裁剪出来的数量之和，i＝0，1，...，th2-th1；
d34)对第三段直方图hu进行裁剪，定义裁剪阈值Tu，
定义裁剪后的直方图为hu′，
为了不改变hu在整体的比例，需要将裁剪多出的部分线性补偿到直方图中，定义补偿后的直方图为hu″，
其中，resu为第三段直方图经过裁剪出来的数量之和，i＝0，1，...，fmax-th2；
d35)定义经过裁剪和补偿的直方图h″，
d4)假设直方图亮度估算模型和标准差估算模型的输出图像G(x，y)的灰度级分割点为th1′和th2′，定义估计输出图像模型的亮度平均值为μm和标准差σm，则通过迭代求解方程得到分割点th1′；根据估计模型可得，
μm＝0.5[r1(th1′-1)+r2(th1′+th2′-1)+r3(th2′+255)]
th2′＝th1′+2σm
d5)通过以上方程，可求得th1，k′为
计算标准差σm，k的方程，
d6)定义迭代计数器k＝1，定义最大迭代次数K，定义预设误差δ，定义初始标准差σm，0，用户定义输入图像平均亮度为μm，0，计算th′1，0＝f1(μm，0，σm，0)，迭代开始；
d7)计算输出图像的标准差σm,k＝f2(th’1,k-1,σm,k-1)；
d8)则更新th’1,k＝f1(μm,0,σm,k)；
d9)若满足或k>K，迭代结束，输出最终的分割点th1′＝t’1,k，th2′＝th1＇+2σm,k；否则令k＝k+1，转到步骤d7)。
d10)计算直方图分段后的低、中、高区域的累积密度函数分别为，
其中，Sl、Sm、Su分别为低中高直方图区域内的像素总数
d11)定义输出图像的动态范围[0,255]的图像映射曲线函数T为，
d12)根据以上可得，从而求出全局增强的输出图像G(x,y)，
G(x,y)＝T(F(x,y))
d9)根据以上可得，从而求出全局增强的输出图像G(x,y),
G(x,y)＝T(F(x,y))。
步骤e包括:
e1)将c和d得到的图像进行线性叠加，得到最终的输出图像O(x,y)，
O(x,y)＝G(x,y)+λ×Ewe(x,y)
其中，λ为比例因子，一般选0到1之间，本具体实施例选其为0.5。
采用上述具体实施方式对原始图像进行增强后，可获得一幅信息量丰富、动态范围较高且视觉效果良好的图像。
基于上述具体实施方式，下面结合具体的实验对本发明的效果加以验证。
如图6和9所示，采集到的两幅质量较差的尺寸为256×256的butterfly图像和尺寸为248×333的fish图像，通过本发明提供的方法对其进行图像增强处理，分别得到信息量丰富、动态范围较高、对比度得到明显提升的输出图像，其视觉效果良好，其中，图7和图10的初始亮度值设定为80，图8和图11的初始亮度值设定为110，用户可根据自己需求来设定初始亮度值。图12和图13给出了本发明方法的输出图像与原图像的对比度和熵的对比，由图可知，通过本发明的方法其对比度和熵均得到提高。
相同或相似的标号对应相同或相似的部件；
附图中描述位置关系的用于仅用于示例性说明，不能理解为对本专利的限制；
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。
Claims (4)
1.一种图像增强方法，其特征在于，包括如下步骤：
a)对输入图像进行去噪处理得到去噪图像；
b)对去噪图像进行边缘提取得到边缘图像；
c)对边缘图像进行图像增强处理得到去噪且边缘增强的图像；
d)利用亮度可控的直方图均衡方法对去噪图像进行处理得到全局增强图像；
e)对c和d步骤所得到的图像进行线性叠加，得到最终的输出图像；
步骤d的具体步骤包括：
d1)求出去噪图像F(x，y)的直方图，并求得其亮度最大值fmax和最小值fmin，其中，0≤[fmin，fmax]≤255以及平均亮度值μ和标准差σ；
d2)根据平均亮度值μ和标准差σ，求出去噪图像F(x，y)直方图的分割点th1和th2为
其中，w为权值，0≤[th1，th2]≤255；
d3)根据fmin、fmax、th1、th2将去噪图像F(x，y)的直方图分成低、中、高三段，如下
其中，h(i)为图像F(x，y)的直方图统计函数，i表示0到255的灰度值；
d4)对直方图进行分段剪切与补偿，得到经过裁剪和补偿后的直方图；
d5)对图像F(x，y)的直方图进行分段剪切与补偿后，各段直方图的像素总数所占图像F(x，y)的总像素的比例均未变化，即
其中，toal代表图像F(x，y)的总像素；r1、r2、r3分别表示去噪图像F(x，y)直方图的低、中、高区域的像素总数占图像F(x，y)的像素总数的比例；
d6)计算直方图分段后的低、中、高区域的累积密度函数分别为
其中，Sl、Sm、Su分别为低中高直方图区域内的像素总数，hl″、hm″、hu″为分段剪切与补偿后各个区域的直方图统计函数；
d7)假设全局增强的输出图像G(x，y)的直方图的分割点分别为th1′和th2′，亮度平均值为μm，标准差为σm，根据输出图像直方图平均亮度估算模型和输出图像直方图标准差估算模型可分别估算出亮度平均值μm标准差σm，即
μm＝0.5[r1(th1′-1)+r2(th1′+th2′-1)+r3(th2′+255)]
其中，th2′＝th1′+2σm；
d8)令模型求出的亮度平均值和设定的亮度平均值mv相等，即σm＝mv，其中，mv是用户可自行设定的平均亮度值；d7步骤中的三个方程组是关于th1 1、th2′和σm三个未知数的方程组，通过迭代的方式，计算出输出图像的直方图灰度分割点th1′和th2′；
d9)定义输出图像的动态范围[0，255]的图像映射曲线函数T为：
其中，th1和th2为自定义的去噪图像F(x，y)的直方图的两个分割点，th1′和th2′为输出图像直方图平均亮度估算模型和标准差估算模型计算出来的输出图像的直方图灰度级分割点；
d9)根据以上可得，从而求出全局增强的输出图像G(x，y)为
G(x，y)＝T(F(x，y))；
步骤c的具体步骤包括：
c1)对边缘图像进行预处理，得到边缘预处理图像，对边缘预处理图像进行图像增强得到边缘增强图像；
c2)对边缘图像进行二值化处理得到二值化图像，对二值化图像进行腐蚀处理得到腐蚀图像；
c3)将c1的边缘增强图像和c2的腐蚀图像进行综合处理得到去噪且边缘增强的图像；
步骤c1的具体步骤包括：
c11)求输入图像I(x，y)的最大和最小灰度值分别为Imax和Imin；
c12)对边缘图像进行预处理，得到边缘预处理图像Ew(x，y)，即
其中，E(x，y)为边缘图像，I(x，y)为输入图像；
c 13)求出边缘预处理图像Ew(x，y)的灰度最大值max、灰度最小值min以及亮度平均值μ0、标准差σ0；
c14)求出边缘预处理图像的直方图，然后求出该直方图灰度值小于0的区域的阈值T1和直方图灰度值大于0的区域的阈值T2；
c15)根据c13和c14中求得的灰度最大值max、灰度最小值min以及图像的阈值T1和T2，将边缘预处理图像Ew(x，y)的直方图分成3个区域(min，T1)、(T1，T2)和(T2，max)，然后对该直方图进行基于均值和标准差剪切的分段直方图均衡，得到边缘增强图像Ee(x，y)；
步骤c14中阈值是通过Rosin方法求取的；
步骤c2和c3的具体步骤包括：
c 21)根据Rosin方法求得边缘图像的阈值Tt，然后运用二值化求得二值图像B(x，y)；
c22)对二值图像进行形态学腐蚀得到腐蚀图像R(x，y)；
c31)结合步骤c1和c2，得到去噪且边缘增强的图像Ewe(x，y)：
2.根据权利要求1所述的图像增强方法，其特征在于，步骤a的具体步骤包括：
a1)将输入图像转变为灰度图像；
a2)对灰度图像进行均值滤波得到平滑后的去噪图像。
3.根据权利要求2所述的图像增强方法，其特征在于，步骤b的具体步骤为利用自定义的拉普拉斯模板对平滑滤波后的图像进行边缘提取得到边缘图像。
4.根据权利要求1至3中任意中任一权利要求所述的图像增强方法，其特征在于，所述步骤e包括：
e1)将c和d得到的图像进行线性叠加，得到最终的输出图像O(x，y)，
O(x，y)＝G(x，y)+λ×Ewe(x，y)
其中，λ为比例因子。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610833349.2A CN106530237B (zh) | 2016-09-19 | 2016-09-19 | 一种图像增强方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610833349.2A CN106530237B (zh) | 2016-09-19 | 2016-09-19 | 一种图像增强方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106530237A CN106530237A (zh) | 2017-03-22 |
| CN106530237B true CN106530237B (zh) | 2019-06-07 |
Family
ID=58343979
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610833349.2A Active CN106530237B (zh) | 2016-09-19 | 2016-09-19 | 一种图像增强方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106530237B (zh) |
Families Citing this family (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107644403B (zh) * | 2017-08-23 | 2021-02-09 | 天津大学 | 恶劣环境条件下图像的非一致颜色校正方法 |
| CN108171661B (zh) * | 2017-12-11 | 2021-10-15 | 南京理工大学 | 基于改进Tri边缘算子的红外目标检测方法 |
| CN108186051B (zh) * | 2017-12-26 | 2021-11-30 | 珠海艾博罗生物技术股份有限公司 | 一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统 |
| CN108681734A (zh) * | 2018-05-28 | 2018-10-19 | 杭州电子科技大学 | 一种基于样本选择的图像相机来源鉴别方法 |
| CN108961179B (zh) * | 2018-06-19 | 2022-10-18 | 上海中和软件有限公司 | 一种医学图像后处理系统及其使用方法 |
| CN109724151A (zh) * | 2018-06-28 | 2019-05-07 | 永康市胜时电机有限公司 | 供暖器进水管道电控机构 |
| CN109064427A (zh) * | 2018-08-01 | 2018-12-21 | 京东方科技集团股份有限公司 | 增强图像对比度的方法、装置、显示设备和存储介质 |
| CN110956587B (zh) * | 2018-09-27 | 2023-07-21 | 中科创达软件股份有限公司 | 一种显示画面的调节方法及装置 |
| CN109544490B (zh) * | 2018-10-17 | 2021-07-13 | 北京达佳互联信息技术有限公司 | 图像增强方法、装置和计算机可读存储介质 |
| CN109410147A (zh) * | 2018-11-06 | 2019-03-01 | 哈尔滨工程大学 | 一种超空泡图像增强方法 |
| CN109378279B (zh) * | 2018-11-12 | 2020-12-18 | 武汉新芯集成电路制造有限公司 | 晶圆检测方法及晶圆检测系统 |
| CN109829860B (zh) * | 2018-12-26 | 2021-02-02 | 武汉高德智感科技有限公司 | 全图与局部相结合的图像线性动态范围压缩方法及系统 |
| CN111626967A (zh) * | 2019-02-28 | 2020-09-04 | 富华科精密工业（深圳）有限公司 | 图像增强方法、装置、计算机装置及可读存储介质 |
| CN110070508B (zh) * | 2019-04-23 | 2020-11-10 | 西安交通大学 | 一种基于门限值和线性映射的反锐化增强方法 |
| CN110310274B (zh) * | 2019-07-02 | 2021-12-17 | 河北农业大学 | 一种植物花朵数量检测方法 |
| CN110544231B (zh) * | 2019-07-24 | 2021-05-11 | 华南理工大学 | 基于背景标准化和集中补偿算法的锂电池电极表面缺陷检测方法 |
| CN110634107B (zh) * | 2019-08-09 | 2022-11-25 | 辽宁蓝枭科技有限公司 | 一种针对t1增强脑部磁共振图像亮度标准化方法 |
| CN110620924B (zh) * | 2019-09-23 | 2022-05-20 | 广州虎牙科技有限公司 | 编码数据的处理方法、装置、计算机设备及存储介质 |
| CN110896499B (zh) * | 2019-10-12 | 2020-10-02 | 南京杰游信息技术有限公司 | 移动终端播放位置定位系统及方法 |
| CN110683240A (zh) * | 2019-10-30 | 2020-01-14 | 西安航空学院 | 一种基于图像处理的垃圾分类处理系统 |
| CN113425319A (zh) * | 2020-03-08 | 2021-09-24 | 广东毓秀科技有限公司 | 面向基层医疗机构的呼吸系统疾病智能辅助诊断决策方法 |
| CN111429373A (zh) * | 2020-03-24 | 2020-07-17 | 珠海嘉润医用影像科技有限公司 | 一种图像增强处理方法及系统 |
| CN115298688A (zh) * | 2020-03-24 | 2022-11-04 | 华为技术有限公司 | 图像增强处理方法和装置 |
| CN112818944A (zh) * | 2021-03-08 | 2021-05-18 | 北方工业大学 | 面向地铁车站场景的密集人群计数方法 |
| CN113112438B (zh) * | 2021-04-29 | 2024-05-07 | 中国船舶重工集团公司第七二三研究所 | 一种基于裁剪直方图的图像增强方法 |
| CN114240927A (zh) * | 2021-12-28 | 2022-03-25 | 湖南云箭智能科技有限公司 | 一种板卡质量检测方法、装置、设备及可读存储介质 |
| CN114936981B (zh) * | 2022-06-10 | 2023-07-07 | 重庆尚优科技有限公司 | 一种基于云平台的场所扫码登记系统 |
| CN115019160B (zh) * | 2022-08-09 | 2022-11-11 | 中工重科智能装备有限责任公司 | 一种工厂仓库入库货物图像质量的控制方法 |
| CN117237243B (zh) * | 2023-11-16 | 2024-01-26 | 大连壹致科技有限公司 | 一种心血管造影图像局部自适应增强方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1741068A (zh) * | 2005-09-22 | 2006-03-01 | 上海广电（集团）有限公司中央研究院 | 一种基于边界的直方图均衡方法 |
| CN101877123A (zh) * | 2009-12-03 | 2010-11-03 | 北京中星微电子有限公司 | 一种图像增强方法和装置 |
| CN104063848A (zh) * | 2014-06-19 | 2014-09-24 | 中安消技术有限公司 | 一种低照度图像增强方法和装置 |
| CN104700376A (zh) * | 2014-11-25 | 2015-06-10 | 桂林电子科技大学 | 基于伽马校正和平滑滤波的图像直方图均衡化增强方法 |
-
2016
- 2016-09-19 CN CN201610833349.2A patent/CN106530237B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1741068A (zh) * | 2005-09-22 | 2006-03-01 | 上海广电（集团）有限公司中央研究院 | 一种基于边界的直方图均衡方法 |
| CN101877123A (zh) * | 2009-12-03 | 2010-11-03 | 北京中星微电子有限公司 | 一种图像增强方法和装置 |
| CN104063848A (zh) * | 2014-06-19 | 2014-09-24 | 中安消技术有限公司 | 一种低照度图像增强方法和装置 |
| CN104700376A (zh) * | 2014-11-25 | 2015-06-10 | 桂林电子科技大学 | 基于伽马校正和平滑滤波的图像直方图均衡化增强方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106530237A (zh) | 2017-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106530237B (zh) | 一种图像增强方法 | |
| CN109816605B (zh) | 一种基于多通道卷积的msrcr图像去雾方法 | |
| CN110889812B (zh) | 一种多尺度融合图像特征信息的水下图像增强方法 | |
| WO2019205707A1 (zh) | 基于暗通道的线状自适应改进全局大气光的图像去雾方法 | |
| KR102104403B1 (ko) | 단일영상 내의 안개 제거 방법 및 장치 | |
| CN106846263B (zh) | 基于融合通道且对天空免疫的图像去雾方法 | |
| Zhou et al. | Retinex-based laplacian pyramid method for image defogging | |
| CN102831591B (zh) | 一种基于高斯滤波的单幅图像的实时去雾方法 | |
| CN111161360B (zh) | 基于Retinex理论的端到端网络的图像去雾方法 | |
| Zhao et al. | Single image fog removal based on local extrema | |
| CN103914820B (zh) | 一种基于图像分层增强的图像去雾方法及系统 | |
| CN107798670B (zh) | 一种利用图像引导滤波器的暗原色先验图像去雾方法 | |
| Pei et al. | Effective image haze removal using dark channel prior and post-processing | |
| Shu et al. | Alternating minimization algorithm for hybrid regularized variational image dehazing | |
| CN112435184B (zh) | 一种基于Retinex和四元数的雾霾天图像识别方法 | |
| CN111598814B (zh) | 基于极端散射通道的单图像去雾方法 | |
| CN117252773A (zh) | 基于自适应颜色校正和导向滤波的图像增强方法及系统 | |
| Kaur et al. | Investigation on quality enhancement of old and fragile artworks using non-linear filter and histogram equalization techniques | |
| CN105608683B (zh) | 一种单幅图像去雾方法 | |
| CN104978720A (zh) | 一种视频图像雨滴去除方法及装置 | |
| CN109003247B (zh) | 彩色图像混合噪声的去除方法 | |
| He et al. | Fast single image dehazing via multilevel wavelet transform based optimization | |
| CN104504659A (zh) | 一种基于提升小波变换的快速iso去噪方法及系统 | |
| CN108416739B (zh) | 基于轮廓波和马尔可夫随机场的交通图像去雾方法 | |
| CN106780398B (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 |