一种图像对比度增强方法
【技术领域】
[0001]本发明属于图像增强技术领域，具体涉及一种基于临界可视偏差增强函数的图像对比度增强方法。
【背景技术】
[0002]图像增强技术是图像处理技术的一种，它可以显著改善图像质量，使得图像内容更有层次感并且主观观测效果更符合人们需求。现实生活中，图像往往存在各种问题，例如:拍照时光圈偏小，导致图像偏暗；场景的对比度较低，而使得图像重点不突出；曝光过度，导致影像失常，照片泛白等。
[0003]目前图像对比度增强的方法主要有:Chang_Hsing Lee等人在文献【Imageenhancementapproach using the just-noticeabIe-difference model of the humancisual system, Journal ofElectronic Imaging, 2012】中将临界可见偏差 JND模型运用到直方图均衡中，从一定程度上克服了直方图均衡带来的图像灰度级减少、细节缺失的弱点，但仍然不能突出图像中重要部分的信息。基于Retinex理论的图像增强算法【许欣，陈强，王平安，孙怀江，夏德深.消除光晕现象的快速Retinex图像增强，计算机辅助设计与图形学学报，2008】对高光区域敏感度小，图像整体颜色泛白，肤色纹理不清晰，高光区域细节信息并没有得到明显改善。基于图像梯度场增强的算法【许欣，陈强，孙怀江，夏德深，结合视觉感知特性的梯度域图像增强方法，计算机辅助设计与图形学学报，2009】结合Weber定律进行图像增强，但Weber定律不适合描述低照度图像的视觉感知特性。【刘琳，基于GPU的图像梯度场增强算法的研宄与实现，电子科技大学硕士论文，2014】直接将JND函数作用于图像梯度场，对图像灰度处于128附近的低对比度有较好的增强效果，但对处于其它灰度等级或图中同时存在有部分质量较好的图像区域时，增强质量下降，鲁棒性差。
【发明内容】
[0004]本发明的目的是提供一种基于临界可视偏差增强函数的图像对比度增强方法，将符合人眼视觉系统的临界可见偏差JND与图像梯度场相结合，设计增强和压扩函数，实现图像增强的方法，解决现有图像增强算法细节缺失、产生光晕现象和增强效果鲁棒性欠缺的问题。
[0005]本发明所要解决的问题是:如何提供一种细节清晰、增强效果良好、鲁棒性好、符合人眼视觉的图像对比度增强方法。
[0006]本发明所提出的技术问题是这样解决的:提供一种基于临界可视偏差增强函数的图像对比度增强方法，其流程如图1所示，主要包括以下步骤:
[0007]步骤101.读取图像R、G、B分量图像，对每个分量图像执行步骤102至步骤110 ；
[0008]步骤102.求取输入图像I (X，y)的梯度▽ I (x, y)；
[0009]步骤103.对输入图像的每个像素点(X，y)，计算其平均背景亮度得到y (x, y)，计算方法为:首先对I (X，y)进行行方向的IX 3均值滤波，均值滤波器模板为[1/31/31/3]的行向量；其次对I(x，y)进行列方向的3X1均值滤波，均值滤波器模板为[1/31/31/3]τ的列向量其中T表示转置，将行向量转置为列向量；
[0010]将所述计算方法总计执行三次，得到μ (X，y)；
[0011]步骤104.根据图像平均背景亮度计算每个像素点的临界可见偏差阈值JND (X，y)，其公式如下:
[0012]若μ (x，y)彡 127，则:
[0013]JND (x, y) = 21* (l—sqrt ( μ (x, y)/127))+4 ;
[0014]否则:
[0015]JND (x, y) = 3* ( μ (x, y)-128)/128+4 ;
[0016]其中，μ (X，y)是步骤103中求出的平均背景亮度值；
[0017]步骤105.将图像梯度与图像临界可见偏差阈值相结合，得到临界可视偏差增强函数ν.(χ，γ)的计算公式，公式如下:
[0018]Vjnd(XJ) = 2-logsig{[ V I (x, y) |-JND(x, y) ]/10}
[0019]其中logsig{.}是S型传输函数；
[0020]步骤106.根据图像梯度，构造梯度压扩函数Φ (X，y)，其具体形式如下:
[0021 ] Φ (X，y) = β X exp {_ β 2/ (2 σ 2)} / σ 2
[0022]其中β = I I V I (χ, y) I |，σ是可调节参数，对不同图像可以对σ进行调整以得到更好的增强效果，σ —般取值为[45，80];
[0023]步骤107.将临界可视偏差增强函数V.(X，y)和梯度压扩函数Φ (x，y)作用于临界可视对比度上，获得增强后的梯度G(x, y) = V I (χ, y) XVjnd(x, y) X Φ (x, y)；
[0024]步骤108.计算增强后梯度场G (χ, y)的散度div(G)；
[0025]步骤109.通过求解泊松方程▽ 2Ilte (X，y) = div (G)，从增强后的梯度重建亮度图像 IKe (x, y)；
[0026]步骤110.对重建的亮度图像像素取值IKe(x，y)线性拉伸至O吣255区间内，其计算方法如下:
[0027]求取图像IKe(x, y)的像素最大值ImaxKe和像素最小值Imin Ee,则线性拉伸后的图像为:
[0028]IEeunify (x, y) = 255* (IEe (x, y) -1minEe) / (ImaxEe-1minEe)
[0029]步骤111.合并经步骤102至步骤110操作后的R、G、B三个分量图像，得到最终增强后的图像。
[0030]进一步的，步骤102所述梯度▽ I (x，y)以前向差分表示，同时，步骤108所述增强后梯度场的散度div(G)用后向差分表示；或者，所述梯度Vl(x，y)以后向差分表示，同时，增强后梯度场的散度div(G)用前向差分表示。
[0031]本发明的有益效果是:
[0032]该图像的基于临界可视偏差增强函数的增强方法将图像梯度和人眼视觉系统的临界可见偏差结合起来，增强函数随着人眼可见偏差的不同而自适应改变；同时，根据梯度幅度值不同而进行不同程度的梯度值压扩，使得增强后的图像细节清晰，噪声抑制和增强效果好，更符合人眼视觉感知特性，增强效果良好。
【附图说明】
[0033]图1是基于临界可视偏差增强函数的图像对比度增强流程图；
[0034]图2是不同背景亮度下的临界可见偏差JND阈值曲线图；
[0035]图3是【具体实施方式】实施例1采用本发明方法对曝光不足的实际图像进行增强的前后对比图:其中3a为原始图像，3b为增强后的图像；
[0036]图4是【具体实施方式】实施例2采用本发明方法对强逆光的实际图像进行增强的前后对比图:其中4a为原始图像，4b为增强后的图像。
【具体实施方式】
[0037]为使本发明的技术方案更加清楚明白，下面结合实例对本发明做进一步的详细说明。
[0038]本发明方法【具体实施方式】包括求图像梯度、获取图像临界可视度、构造增强与压扩函数以及通过求解泊松方程重建增强后亮度图像，具体步骤如下:
[0039]实施例1:以实际拍摄的曝光不足图像为例:
[0040]本实施例选取实际拍摄的曝光不足导致画面偏暗的图像作为示例，如图3a所示，其图像分辨率为640*480 (像素)；增强函数的可调节参数α = 0.3。
[0041]本实施例的具体实现过程如下:
[0042]步骤101.读取原始图像I (X，y)的R、G、B三个颜色分量图像，对每一个分量执行步骤102至步骤110，对应于图像中位置于(320，240)的像素点，其R、G、B三个分量值分别为 Ie(320, 240) = 25、Ig(320，240) = 30、Ib(320，240) = 34 ；
[0043]步骤102.求图像I (X，y)的梯度▽ I (x，y)，以一阶前向差分近似求出:
[0044]V I (x, y) ^ (I (x+1, y) -1 (x, y)，I (x, y+1) -1 (χ, y))；
[0045]此时图像中位置于(320，240)处的R、G、B三个通道的χ方向、y方向梯度值分别为▽ Ik (320，240) = (32，7)、▽ Ig (320，240) = (41，I)、▽ Ib (320，240) = (50, 3)；
[0046]步骤103.求图像I (X，y)的平均背景亮度μ (x, y):对应于图像中(320，240)处的R、G、B三个通道的背景亮度值分别为μΕ(320, 240) = 12.65、yG(320, 240) = 11.88、μΒ(320, 240) = 10.99 ；
[0047]步骤104.根据图像平均背景亮度y (x，y)，求图像在(x，y)处的临界可见偏差JND阈值，计算公式如下: