-
如果设备性能较差,或照明条件不理想,数字图像的对比度往往很低。低对比度的图像看起来非常模糊,不能充分展示图像的细节,因此研究人员提出了各种算法来提升这些图像的视觉效果。
对比度增强技术大致可分为直接增强和间接增强两类。文献[1]是比较典型的直接增强算法,通过应用各种非线性函数[2]或求解优化问题[3]来改进图像对比度。这些增强效果是通过人眼视觉系统(human visual system, HVS)来衡量,如Weber-Fechner定律或Retinex理论[4]等。此方法在图像细节增强和动态范围压缩方面有一些优势,但计算复杂度较高,并导致了“晕”伪影,特别是在强边缘[5]周围,“晕”伪影尤其明显。尽管近几年提出的直接方法,缓解了这些问题,但仍然缺少能够兼顾高图像对比度和实时计算效率的算法。因此,基于全局变换函数的间接方法有着更为广泛的应用。
在间接方法方面,直方图均衡化(histogram equalization, HE)是最具代表性的方法之一。当图像直方图中出现高峰值时,会产生视觉伪影,如噪声放大、轮廓化或显著的亮度变化等。为了缓解这些问题,基于HE的方法对图像直方图进行修改,使其峰值衰减,然后由修改后的直方图得到变换函数。如文献[6]通过求解双准则优化问题,将图像直方图与均匀分布相结合;文献[7]使用归一化幂律函数,平滑图像直方图中的峰值。由于间接方法能够将增强后图像的平均亮度转移到中等灰度级别,因此将图像直方图分割成多个子直方图,并对每个子直方图分别进行HE处理[8]。为了保留细节,文献[9]提出了利用二维(2D)直方图均衡化的间接方法,模糊上下文信息来构建图像直方图,赋予纹理区域的像素更高的权重。与其他间接方法相比,此方法具有更好的对比度增强性能和更少的信息损失。
以上这些方法在低照度的情况下,对图像的处理能力较差,难以达到理想的效果。这主要是因为在低照度情况下,人眼对图像处理的方式与正常照度下的模式不同。为此,本文首先引入了一种新的基于人类亮度感知的评价对比敏感度的S型函数。由于人类视网膜的对比敏感度随着亮度(取对数)的增加呈指数级下降,因此将对比敏感度建模为对数参数的指数函数。利用对比灵敏度模型,通过修改Steven幂律的指数,导出S型函数。此外,本文还提出了一种参数优化方法,在保持输入图像的平均亮度和扩展图像的直方图的同时使信息损失最小化。
-
人类视网膜上的光感受器包括视杆细胞和视锥细胞两种。视杆细胞对光线很敏感,在低亮度水平下能够提供消色差视觉,这种视觉也称为暗视觉(10−6~10 cd/m2)。在亮度水平高于10−2 cd/m2时,视杆细胞开始饱和,视锥细胞提供名为明视觉的色觉。在亮度水平10−2~10 cd/m2之间,视杆细胞和视锥细胞都活跃,人类视网膜处于中光视觉的过渡模式。因为神经元只能传输一个动态范围大约为1:103的信号。所以人的视网膜工作原理是首先压缩真实场景亮度的动态范围,然后将光感受器适应到一定的亮度水平(即适应水平),在适应水平附近的一个相当小的动态范围内感知图像。为了描述视网膜对人类亮度感知的反应,提出了各种基于神经科学实验的响应模型。
Naka-Rushton方程[10]是比较有代表性的响应模型。此方程描述了视网膜响应度R与亮度水平L的关系,具体公式如下:
$$ R(L) = \frac{{{L^n}}}{{{L^n} + {\sigma ^n}}} $$ (1) 式中,n为决定视网膜反应函数陡度的参数;σ为适应水平。
式(1)表示人眼视觉系统中采用以适应水平为中心的S形曲线函数,此函数将亮度水平转换为视网膜响应度。近几年的视觉映射技术一般都采用类似的S型曲线响应函数,将现实世界的光亮度转换成像素亮度。采用基于人眼亮度感知的S型曲线函数,在避免视觉伪影的同时,提高了对比度。
-
许多间接的方法都是利用人的视网膜所能感知到的最小亮度差异(just noticeable difference, JND)来推导变换函数。根据Weber-Fechner定律,JND(表示为dL/dS)与背景亮度L成正比:
$$ \frac{{{\rm{d}}L}}{{{\rm{d}}S}} = kL $$ (2) 式中,k是一个感知常数,为Weber分数。求解式(2)的常微分方程可知,感知到的亮度S可表示为:
$$ S = \frac{1}{k}\ln L + {S_0} $$ (3) 式中,S0是积分常数。各种间接的方法使用对数函数[11]或GC曲线[12]作为变换函数。然而,由于Weber-Fechner定律仅适用于亮度大于100 cd/m2的亮度水平[13],对于低亮度水平下的视网膜响应并不匹配。因此,通常会在黑暗区域导致对比度过度增强。为了解决这一问题,文献[14]提出了一种利用图像直方图确定参数的自适应伽马校正(adaptive Gamma correction, AGC)方法。此方法中变换函数R(L)为:
$$ R(L) = {\left( {\frac{L}{{{L_{\max }}}}} \right)^{1 - {\rm{CDF}}(L)}} $$ (4) 式中,Lmax为最大像素亮度;CDF(L)为输入直方图的累积分布函数。与其他间接方法相比,AGC方法表现出更好的视觉匹配效果。对于模糊图像,这种效果尤为明显。然而,此方法提高了图像的平均亮度,在亮区损失了一定的对比度。为了解决上述局限性,文献[15]通过修改GC曲线来构建S型曲线,但很难同时达到令人满意的对比度增强效果和平均亮度。
为此,文献[16]提出了一种S型函数的变换函数R(L),定义为:
$$ R(L) = L \left( {1 + C \frac{1}{{1 + \exp \left( { - L} \right)}}} \right) $$ (5) 式中,C为尺度参数,用来确定增强程度。受此启发,文献[17]通过修改式(5)扩展了文献[14]中的工作,如下所示:
$$ R\left( L \right) = L + {K_1} \frac{L}{{1 - \exp \left( {{K_1} \left( {{K_2} + L} \right)} \right)}} $$ (6) 式中,K1和K2是调优参数。虽然这些方法证明了S型曲线增强对比度的可行性,但它们的性能仍然较差,图像失真严重。
-
Stevens幂次定律是一个著名的刺激−反应模型,与Weber-Fechner定律相比,它涵盖了更广泛的感觉范围。在Stevens幂次定律中,感知亮度R(L)函数定义为:
$$ R\left( L \right) = {L^k} $$ (7) 式中,指数k取决于刺激的类型。式(7)以微分的形式可写为:
$$ \frac{{\text{1}}}{R}{\rm{d}}R = k\frac{1}{L}{\rm{d}}L $$ (8) 式中,k用作灵敏度参数。此参数决定了随着刺激强度的增加,视网膜响应度增长的速度。对于不同类型的视网膜响应度,k假设为一个常数,即常规GC曲线。然而,由于人类视网膜响应度是由背景亮度自适应决定的。因此,本文根据亮度感知情况,建立k模型。文献[18]的实验表明,人类视网膜亮度的对比响应度(即Weber分数)随着对数级亮度的增加呈指数式递减。因此,Weber分数可以近似拟合为对数级亮度的指数函数。基于这些研究,本文将k建模为以下指数函数:
$$ k\left( L \right) = \alpha {\beta ^{ - \ln \left( L \right)}} $$ (9) 式中,α和β分别为确定k(L)的最大值和陡度的参数。变换函数R(L)为:
$$ R\left( L \right) = 255 {\tilde L^{\alpha {\beta ^{ - \ln \left( {\tilde L} \right)}}}} $$ (10) 式中,L为原始像素强度,取值范围[0, 255];
$ \tilde L $ 为归一化到[0,1]的像素强度,即将单个独立的像素强度处理为与整体图像相关的对比度,具体为:$$ \tilde L = \frac{{L - {L_{\min }}}}{{{L_{\max }} - {L_{\min }}}} $$ (11) 式中,Lmin和Lmax分别为输入图像的最小和最大亮度等级。式(10)中的变换函数覆盖数字图像的全动态范围
$ \left[ {R\left( {{L_{\min }}} \right) = 0,R\left( {{L_{\max }}} \right) = 255} \right] $ 。 -
本节分析了式(10)中变换函数的两个参数 α 和 β对合成图像的影响,并引入一种新的自动参数估计方法。首先,将α从0.3逐步变化到1.1,每一次的增量为0.2,同时保持β为2.0,得到的输出图像如图1所示,当α增大时,输出图像变暗。由于 α 限制了图像中各点亮度的最大值,可通过控制 α 来调节输出图像的平均亮度。但α需要约束来避免由平均亮度变化[19]引起的图像失真和闪烁伪影。为保持原始图像的平均强度水平,修正变换函数为:
$$ R\left( {{L_{{\rm{mean}}}}} \right) = 255 {\tilde L_{{\rm{mean}}}}^{\alpha {\beta ^{ - \ln \left( {{{\tilde L}_{{\rm{mean}}}}} \right)}}} = {L_{{\rm{mean}}}} $$ (12) 式中,Lmean为归一化输入图像的平均强度等级。即保证经过转换后的图像平均亮度与原始图像的平均亮度相一致,在此情况下,通过求解式(12),可得到保持平均亮度的最优α:
$$ \alpha {\text{ = }}\frac{{\ln \left( {{L_{{\rm{mean}}}}/255} \right)}}{{\ln \left( {{{\tilde L}_{{\rm{mean}}}}} \right)}} {\beta ^{\ln \left( {{{\tilde L}_{{\rm{mean}}}}} \right)}} $$ (13) β对输出图像的影响如图2所示。如图2g所示,通过式(13)得到的变换函数保持输入图像的平均亮度水平,因而不需要考虑β值。
通过确定变换函数的陡度来控制对比度增强程度,当值太小时,此方法的增强性能不理想;如果太大,变换函数会出现陡坡,导致明暗区域的细节丢失明显。因此,为了保持图像细节的同时提高图像对比度的最优值,本文定义:
$$ E\left( {\beta ,\lambda } \right) = {E_{\rm{d}}}\left( \beta \right) - \lambda {E_{\rm{e}}}\left( \beta \right) $$ (14) 式中,Ed(β)为数据保真项;Ee(β)为增强项;λ为控制数据保真项与增强项之间权衡的正则化参数。
为了防止输出像素值被截断而导致信息丢失,Ed定义为离散熵(DE)损失:
$$ {E_{\rm{d}}}\left( \beta \right) = - \sum\limits_{k = 1}^N {p\left( {{h_{i,k}}} \right){{\log }_2}p\left( {{h_{i,k}}} \right)} + \sum\limits_{k = 1}^N {p\left( {{h_{o,k}}} \right){{\log }_2}p\left( {{h_{o,k}}} \right)} $$ (15) 式中,p(·)为概率质量函数;hi为原始直方图的第k个分量;ho为使用式(10)在给定β的条件下得到的输出直方图;N为直方图的分量总数。
对于增强项,本文采用像素区方法来测量图像中所有像素的灰度差。具体为:
$$ {E_{\rm{e}}}\left( \beta \right) = \frac{1}{{N\left( {N - 1} \right)}}\sum\limits_{i = 1}^N {\sum\limits_{j = 1}^N {{h_{\rm{o}}}\left( i \right)} } {h_{\rm{o}}}\left( j \right)\left( {j - i} \right) $$ (16) 式中,N为直方图分量的总个数;i、j表示对应分量的索引。
当图像直方图分布均匀而不集中于特定灰度层时,该项得分较高。为了找到最佳的β,最小化式(14)中的函数E(β),本文采用黄金分割搜索算法评估算法,估值为三元组的点值形式的黄金比例,逐步缩小搜索范围。本文算法中使用初始区间[1.2, −1.4]来实现迭代过程,直到间隔
${\;\beta _u} - {\beta _l} $ <10−4或达到最大迭代次数Kmax=500时停止。 -
通过取3个数据集[20]的550个测试图像,对比本文方法——基于灵敏度的S型曲线模型(image contrast enhancement algorithm with S-type function based on human visual system, ICEA)和其他方法。如加权自适应直方图均衡化(weighted adaptive histogram equalizatio, WAHE)、上下文和变分对比(contextual and variational contrast, CVC)[14]、分层差异表示、自适应图像灰度校正(adaptive gamma correction, AGC)[18]、模糊背景对比度增强(fuzzy-contextual contrast enhancemen, FCCE)[21]、两个S型函数的方法[16-17]和其他一些有代表性的直接方法[22-24]。为保证验证的公平性,所有算法都使用Intel i5-3550 3.30 GHz CPU和16 GB RAM的台式机测试每个数据集的平均处理时间,使用python进行效果模拟,并且python的版本和其他依赖库的版本都统一。
为了确定参数α的最优值,首先在Berkeley图像数据集上测试了本文方法,将λ从0改变为0.05,每一步的变化为0.001,如图3所示。当λ较大时,本文方法在像素区分数上有很大的改善,但在DE分数上的细节损失较大。为了实现更好的增强像素区分数,同时防止细节丢失,设λ=0.015。首先将图像转换到HSV颜色空间,以防止颜色失真,对亮度通道V进行处理,H和S通道保持不变。
为了进行客观评估,本文采用了对比度增强(enhancement measure evaluation, EME)[2]、绝对平均亮度误差(absolute mean brightness error, AMBE)[25]、梯度大小相似度偏差(gradient magnitude similarity deviation, GMSD)[26]和离散熵(discrete entropy, DE)措施。
首先,使用EME测量基于Weber算法的图像平均对比度,越大说明图像整体的对比度越好。对于EME,将增强后的图像X划分为N个子块Xi,j,其大小固定相同,W为横向子块的个数,H为纵向子块的个数,计算每个子块的灰度最大值与最小值之比。然后计算平均比例作为最终得分。EME为:
$$ {\rm{EME}}\left( X \right) = \frac{1}{{W H}}\sum\limits_{i = 1}^W {\sum\limits_{j = 1}^H {20\ln \frac{{\max ( {{x_{ij}}} )}}{{\min ( {{x_{ij}} + \delta } )}}} } $$ (17) 式中,max(xi,j)和min(xi,j)分别是子块xi,j中像素强度的最大值和最小值。δ设定为0.000 1,目的是为了避免除零操作;每个子块的大小为8×8。
其次,测量输入图像和增强图像之间的平均像素强度的绝对差值为:
$$ {\rm{AMBE}}(X,Y) = \frac{1}{{W H}}\sum\limits_{i = 1}^W {\sum\limits_{j = 1}^H {\left| {Y\left( {i,j} \right) - X\left( {i,j} \right)} \right|} } $$ (18) 式中,Y是输入图像;X是增强后的图像;其他参数如之前所述。AMBE评分越低,说明相应的增强方法有效地保持了输入图像的平均亮度,增强后的图像更好地还原了原始图像的整体亮度特性,对人眼有更好的适应性。
再使用GMSD计算输入图像和增强图像之间像素的梯度相似度,然后测量整体梯度相似度的标准差作为最终得分。GMSD得分越低,说明输入图像与增强图像之间的图像失真越小:
$$ {\rm{GMSD}} = \sqrt {\frac{1}{N}\displaystyle\sum\limits_{i = 1}^N {{{({\rm{GMS}}(i) - {\rm{GMSM}})}^2}} } $$ (19) 式中,GMS(i)为第i个子块的梯度幅值相似性;GMSD是全部局域梯度场的平均处理结果。
最后,使用DE评估图像中的信息量。DE值高表示图像包含更多信息。由于信息处理不等式的存在,利用全局变换函数生成的输出图像的DE不能高于输入图像。
$$ H(X) = - \sum\limits_{i = 1}^n {{P_i}\lg ({P_i})} $$ (20) 总的来说,直接算法整体比间接算法在人眼观测的效果上更好,特别是在对细节信息的保留上有更好的表现。其中LDR等间接算法在极大地增强细节的同时会产生光晕的效果。在图像中对比度较强的区域,颜色过于鲜艳,影响图像的整体效果。
表1~表3为依据文献[20]的3个数据集,对常见图像处理算法进行分项指标计算的实验结果。结果表明,本文算法表现出了一定的优势。根据EME结果,LSCN和本文方法在各个类别中表现出了最好的性能,避免了AMBE和GMSD评分中显示的平均亮度变化和图像失真。但是,LSCN会导致GMSD分数中明显的图像失真,需要更高的计算复杂度。此外,RSIE的性能与本文方法相当,但计算复杂度高出约50倍。
表 1 Database 1
Categories Algorithms EME AMBE GMSD DE Processing time/ms Direct Methods NPEA[22] 11.50 27.34 0.0434 7.1029 9631.2 LSCN[23] 77.56 14.89 0.0987 6.4856 6069.3 RSIE[24] 24.10 5.87 0.0252 7.4942 302.4 Indirect Methods WAHE[27] 15.83 5.98 0.0385 7.1203 5.4 CVC[14] 22.32 8.15 0.0324 7.1104 103.5 LDR[28] 24.34 8.65 0.0285 7.1601 18.6 AGC[18] 15.43 41.33 0.0605 6.8903 1.9 FCCE[21] 21.34 9.56 0.0674 7.5392 128.4 CESF[16] 13.67 50.03 0.1093 5.8912 0.3 EACE[17] 24.77 10.92 0.1284 7.0498 35.8 SMSF
(proposed)26.32 4.88 0.0256 7.1953 5.3 表 2 Database 2
Categories Algorithms EME AMBE GMSD DE Processing time/ms Direct Methods NPEA[22] 13.24 21.73 0.0409 7.1945 6531.2 LSCN[23] 75.34 17.21 0.0895 6.7834 7469.3 RSIE[24] 21.34 11.35 0.0321 7.4654 392.4 Indirect Methods WAHE[27] 16.33 8.63 0.0402 7.0634 3.4 CVC[14] 21.35 7.99 0.0492 7.1671 124,7 LDR[28] 21.54 9.03 0.0231 8.0842 12.3 AGC[18] 16.89 39.04 0.0597 6.9081 1.2 FCCE[21] 20.56 7.88 0.0740 7.0909 124.7 CESF[16] 14.01 49.02 0.1205 5.7812 0.1 EACE[17] 25.81 13.75 0.1107 7.0938 33.8 SMSF
(proposed)28.30 5.84 0.0209 7.9912 5.1 表 3 Database 3
Categories Algorithms EME AMBE GMSD DE Processing time/ms Direct Methods NPEA[22] 12.43 22.56 0.0511 7.1923 4352.2 LSCN[23] 87.34 12.77 0.0985 6.3456 1345.3 RSIE[24] 31.01 13.46 0.0294 7.7653 1234 Indirect Methods WAHE[27] 15.88 8.34 0.0502 7.0543 3.3 CVC[14] 25.21 9.62 0.0496 7.2345 113,7 LDR[28] 27.49 11.23 0.0203 8.1347 14.5 AGC[18] 17.83 38.53 0.0653 6.9854 0.7 FCCE[21] 21.37 9.21 0.0842 7.1344 12.3 CESF[16] 16.72 47.28 0.1198 5.9654 0.2 EACE[17] 22.89 16.23 0.1209 7.1045 12.4 SMSF
(proposed)26.57 5.54 0.0365 7.9342 3.4 在传统的间接方法中,LDR、FCCE和EACE在EME评分方面体现了高对比度。然而,如AMBE评分所示,LDR在平均亮度保持方面表现较弱。对于EACE和FCCE,输出图像GMSD得分显示图像失真严重。相比之下,从EME评分来看,本文方法的增强性能比第二优的EACE算法提高了10.5%。此外,AMBE和GMSD评分显示,与EACE相比,此方法平均亮度变化减少60.2%,图像失真减少76.9%。如DE得分所示,在间接方法中,本文方法是保存细节方面效果第二的方法。虽然FCCE提供了最高的DE分数,但FCCE包含了局部处理,其中输出像素强度是结合输入像素强度和变换函数结果自适应获得。另一方面,此方法在不进行任何局部处理的情况下,仅获得了较一般的DE得分。
本文方法实现了更高的对比度增强性能,降低了平均亮度的变化,需要更少的计算复杂度,实现了较低的图像失真。
Image Contrast Enhancement Algorithm with S-Type Function Based on Human Visual System
-
摘要: 为了克服传统变换函数在低照度情况下的局限性,在间接对比度增强领域,提出了一种基于人眼亮度感知对比度灵敏性的S型函数。对于不同的图像亮度,存在不同的视网膜响应值,因此将人眼视网膜的对比度灵敏性建模为对数参数的指数函数。该方法以灵敏度模型作为Steven幂律的指数,推导出一个感知亮度的转换函数。同时还提出了一种参数优化方法,在保持输入图像的平均亮度和直方图的同时,保持信息损失最小。实验结果表明,该方法在保持输入图像平均亮度的情况下,具有更少的信息损失和更低的计算复杂度。在对比度增强、平均亮度保持和细节保持方面具有优势。Abstract: In the field of indirect contrast enhancement, an S-type function based on the contrast sensitivity of human visual system is proposed, in order to overcome the limitations of traditional transformation functions under low illumination. This method models the contrast sensitivity of the human retina as an exponential function of logarithmic intensity, since there are different retinal response values for different stimulus intensities. This method uses the sensitivity model as the exponent of Steven's power law to derive a transfer function for perceived brightness. At the same time, a parameter optimization method is proposed, which maintains the average brightness of the input image, and stretches the image histogram while ensuring minimal information loss. Experimental results show that this method has less information loss and better computational complexity performance while maintaining the average brightness of the input image. It has certain advantages in existing methods in terms of contrast enhancement performance, average brightness and preservation of details.
-
Key words:
- contrast enhancement /
- sensitivity model /
- S-type function /
- Steven power-law
-
表 1 Database 1
Categories Algorithms EME AMBE GMSD DE Processing time/ms Direct Methods NPEA[22] 11.50 27.34 0.0434 7.1029 9631.2 LSCN[23] 77.56 14.89 0.0987 6.4856 6069.3 RSIE[24] 24.10 5.87 0.0252 7.4942 302.4 Indirect Methods WAHE[27] 15.83 5.98 0.0385 7.1203 5.4 CVC[14] 22.32 8.15 0.0324 7.1104 103.5 LDR[28] 24.34 8.65 0.0285 7.1601 18.6 AGC[18] 15.43 41.33 0.0605 6.8903 1.9 FCCE[21] 21.34 9.56 0.0674 7.5392 128.4 CESF[16] 13.67 50.03 0.1093 5.8912 0.3 EACE[17] 24.77 10.92 0.1284 7.0498 35.8 SMSF
(proposed)26.32 4.88 0.0256 7.1953 5.3 表 2 Database 2
Categories Algorithms EME AMBE GMSD DE Processing time/ms Direct Methods NPEA[22] 13.24 21.73 0.0409 7.1945 6531.2 LSCN[23] 75.34 17.21 0.0895 6.7834 7469.3 RSIE[24] 21.34 11.35 0.0321 7.4654 392.4 Indirect Methods WAHE[27] 16.33 8.63 0.0402 7.0634 3.4 CVC[14] 21.35 7.99 0.0492 7.1671 124,7 LDR[28] 21.54 9.03 0.0231 8.0842 12.3 AGC[18] 16.89 39.04 0.0597 6.9081 1.2 FCCE[21] 20.56 7.88 0.0740 7.0909 124.7 CESF[16] 14.01 49.02 0.1205 5.7812 0.1 EACE[17] 25.81 13.75 0.1107 7.0938 33.8 SMSF
(proposed)28.30 5.84 0.0209 7.9912 5.1 表 3 Database 3
Categories Algorithms EME AMBE GMSD DE Processing time/ms Direct Methods NPEA[22] 12.43 22.56 0.0511 7.1923 4352.2 LSCN[23] 87.34 12.77 0.0985 6.3456 1345.3 RSIE[24] 31.01 13.46 0.0294 7.7653 1234 Indirect Methods WAHE[27] 15.88 8.34 0.0502 7.0543 3.3 CVC[14] 25.21 9.62 0.0496 7.2345 113,7 LDR[28] 27.49 11.23 0.0203 8.1347 14.5 AGC[18] 17.83 38.53 0.0653 6.9854 0.7 FCCE[21] 21.37 9.21 0.0842 7.1344 12.3 CESF[16] 16.72 47.28 0.1198 5.9654 0.2 EACE[17] 22.89 16.23 0.1209 7.1045 12.4 SMSF
(proposed)26.57 5.54 0.0365 7.9342 3.4 -
[1] 张新明, 沈兰荪. 基于小波的同态滤波器用于图像对比度增强[J]. 电子学报, 2001, 29(4): 531-533. doi: 10.3321/j.issn:0372-2112.2001.04.025 ZHANG X M, SHEN L S. Image contrast enhancement by wavelet based Homomorphic filtering[J]. Acta Electronica Sinica, 2001, 29(4): 531-533. doi: 10.3321/j.issn:0372-2112.2001.04.025 [2] 李冠章, 罗武胜, 李沛. 一种高效地修正Retinex图像自适应对比度增强算法[J]. 测试技术学报, 2009, 23(5): 445-451. doi: 10.3969/j.issn.1671-7449.2009.05.014 LI G Z, LUO W S, LI P. An effective modified Retinex algorithm for adaptive contrast enhancement of image[J]. Journal of Test and Measurement Technology, 2009, 23(5): 445-451. doi: 10.3969/j.issn.1671-7449.2009.05.014 [3] NERCESSIAN S C, PANETTA K A, AGAIAN S S. Non-Linear direct multi-scale image enhancement based on the luminance and contrast masking characteristics of the human visual system[J]. IEEE Transactions on Image Processing, 2013, 22(9): 3549-3561. doi: 10.1109/TIP.2013.2262287 [4] 王书民, 张爱武, 胡少兴, 等. 一种基于全局和局部光照估计的Retinex图像增强算法[J]. 科学技术与工程, 2015, 15(34): 100-104. doi: 10.3969/j.issn.1671-1815.2015.34.017 WANG S M, ZHANG A W, HU S X, et al. The Retinex algorithm of image enhancement based on the global and local illumination estimation[J]. Science Technology and Engineering, 2015, 15(34): 100-104. doi: 10.3969/j.issn.1671-1815.2015.34.017 [5] 程宇, 邓德祥, 颜佳, 等. 基于卷积神经网络的弱光照图像增强算法[J]. 计算机应用, 2019, 39(4): 1162-1169. doi: 10.11772/j.issn.1001-9081.2018091979 CHENG Y, DENG D X, YAN J, et al. Study on low-light image enhancement based on convolutional neural network[J]. Journal of Computer Applications, 2019, 39(4): 1162-1169. doi: 10.11772/j.issn.1001-9081.2018091979 [6] CELLIK T, TJAHJADI T. Automatic image equalization and contrast enhancement using gaussian mixture modeling[J]. IEEE Transactions on Image Processing, 2011, 21(1): 145-156. [7] KIM M, CHUNG M. Recursively separated and weighted histogram equalization for brightness preservation and contrast enhancement[J]. IEEE Transactions on Consumer Electronics, 2008, 54(3): 1389-1397. doi: 10.1109/TCE.2008.4637632 [8] WANG Y, CHEN Q. Image enhancement based on equal area dualistic sub-image histogram equalization method[J]. IEEE Transactions on Consumer Electronics, 1999, 45(1): 68-75. doi: 10.1109/30.754419 [9] 王冰雪, 刘广文, 刘美, 等. 轮廓波域内局部对比度增强的彩色图像灰度化算法[J]. 液晶与显示, 2020, 35(2): 151-160. WANG B X, LIU G W, LIU M, et al. Color image graying algorithm with local contrast Enhancement in contour wave domain[J]. Chinese Journal of Liquid Crystals and Displays, 2020, 35(2): 151-160. [10] YUKIHIRO M, KAZUO T, TOSHIHIDE K. Effect of midazolam, medetomidine, and butorphanol tartrate combination anesthetic on electroretinograms of mice[J]. Molecular Vision, 2019, 25(3): 645-653. [11] WANG C, YE Z. Brightness preserving histogram equalization with maximum entropy: A variational perspective[J]. IEEE Transactions on Consumer Electronics, 2005, 51(4): 1326-1334. doi: 10.1109/TCE.2005.1561863 [12] 吕宗伟, 唐治德, 周林, 等. 利用累积分布函数的亮度保持均衡算法[J]. 计算机辅助设计与图形学学报, 2011, 23(9): 1495-1503. LYU Z W, TANG Z D, ZHOU L. Flattest histogram specification with accurate brightness preservation[J]. Journal of Computer-Aided Design & Computer Graphics, 2011, 23(9): 1495-1503. [13] CELIK T, TJAHJADI T. Contextual and variational contrast enhancement[J]. IEEE Transactions on Image Processing, 2011, 20(12): 3431-3441. doi: 10.1109/TIP.2011.2157513 [14] HUANG S C, CHANG F C, CHIU Y S. Efficient contrast enhancement using adaptive gamma correction with weighting distribution[J]. IEEE Transactions on Image Processing, 2013, 22(3): 1032-1041. doi: 10.1109/TIP.2012.2226047 [15] CELIK T, LI H C. Residual spatial entropy-based image contrast enhancement and gradient-based relative contrast measurement[J]. Journal of Modern Optics, 2016, 63(16): 1600-1617. doi: 10.1080/09500340.2016.1163427 [16] HASSAN N, AKAMATSU N. A new approach for contrast enhancement using sigmoid function[J]. International Arab Journal of Information Technology, 2004, 1(2): 221-225. [17] LAL S, CHANDRA M. Efficient algorithm for contrast enhancement of natural images[J]. International Arab Journal of Information Technology, 2014, 11(1): 95-102. [18] WYSZECKI G W, STILES W S. Color science: Concepts and methods, quantitative data and formulae[J]. Physics Today, 1968, 21(6): 83-84. [19] KWOK N M, HA Q P, LIU D, et al. Contrast enhancement and intensity preservation for gray-level images using multi-objective particle swarm optimization[J]. IEEE International Conference on Automation Science & Engineering, 2009, 6(1): 145-155. [20] RAHMAN S, RAHMAN M M, ABDULLAH A W M, et al. An adaptive gamma correction for image enhancement[J]. Eurasip Journal on Image & Video Processing, 2016, 5(1): 35-43. [21] PARIHAR A S, VERMA O P, KHANNA C. Fuzzy-contextual contrast enhancement[J]. IEEE Transactions on Image Processing, 2017, 26(4): 1810-1819. doi: 10.1109/TIP.2017.2665975 [22] WANG S, ZHENG J, HU H M, et al. Naturalness preserved enhancement algorithm for non-uniform illumination images[J]. IEEE Transactions on Image Processing, 2013, 22(9): 3538-3548. doi: 10.1109/TIP.2013.2261309 [23] ZHAN K, SHI J, TENG J, et al. Linking synaptic computation for image enhancement[J]. Neurocomputing, 2017, 238(5): 1-12. [24] FU X, WANG J, ZENG D, et al. Remote sensing image enhancement using regularized-histogram equalization and DCT[J]. IEEE Geoence & Remote Sensing Letters, 2015, 12(11): 2301-2305. [25] CHEN S D, RAMLI A R. Minimum mean brightness error bihistogram equalization in contrast enhancement[J]. IEEE Transactions on Consumer Electronics, 2003, 49(4): 1310-1319. doi: 10.1109/TCE.2003.1261234 [26] XUE W, ZHANG L, MOU X, et al. Gradient magnitude similarity deviation: A highly effificient erceptual image quality index[J]. IEEE Transactions on Image Processing, 2014, 23(2): 684-695. doi: 10.1109/TIP.2013.2293423 [27] ARICI T, DIKBAS S, ALTUNBASAK Y. A histogram modification framework and its application for image contrast enhancement[J]. IEEE Trans Image Process, 2009, 18(9): 1921-1935. doi: 10.1109/TIP.2009.2021548 [28] LEE C, LEE C, KIM C S. Contrast enhancement based on layered difference representation of 2D histograms[C]//IEEE Transactions on Image Processing. [S.l.]: IEEE, 2013,