本发明属于灰度图像质量评价技术领域，具体的说是一种基于人眼视觉特性的图像增强评价方法。
背景技术：
图像增强主要用于提高图像质量，使图像更有益于人眼观察。根据图像类型(雾霾天气成像、逆光条件成像、光照不足条件成像等)的不同，图像增强算法也是种类繁多。由于缺乏统一的评价标准，如何根据实际需求选择相应的增强算法往往缺乏理论依据。如何评价这些算法性能的优劣正是本发明要解决的问题。
图像增强算法的评价是建立在图像质量评价的基础之上，如果图像增强前后，图像质量有一定的提高，图像增强算法的使用是有益处的。根据实施评价的主体不同图像增强评价方法可以分为主观评价方法和客观评价方法，一般前者需要众多测试人员参与，分别对图像进行打分，用统计值作为图像质量的得分，十分耗时且人力投入较大；后者是通过软件实现图像质量的自动打分，效率高易操作，是目前的主流方法。人眼视觉系统是图像质量评价的主要依据，图像质量的评价应该是建立在人眼视觉特性的基础之上。而现有图像质量评价方法(如VIF算法、SSIM算法、EMEE算法等)常常避免使用人眼视觉模型，主要是因为该模型比较复杂。
技术实现要素：
针对现有技术中存在的上述不足之处，本发明考虑与人眼亮度掩盖效应和对比度掩盖效应有关的临界可见偏差，在临界可见偏差的基础上提供一种基于人眼视觉特性的图像增强评价方法。
本发明为实现上述目的所采用的技术方案是：一种基于人眼视觉特性的图像增强评价方法，包括以下步骤：
扫描原始图像与增强后图像，计算每一像素点(i0,j0)所在原始图像、增强图像的局部背景灰度差别、局部背景灰度，和原始图像与增强图像的相关系数：
所述局部背景灰度差别是指除去中心位置(i0,j0)像素，邻域内其他像素对应的灰度差别中最大的灰度差别；
针对原始图像与增强图像，分别根据局部背景灰度差别与局部背景灰度均值计算每一个像素位置上的临界可见偏差；
根据上述计算结果进行累计打分：在累计得分之前变量F与变量H的初始值为零；
针对每一个像素，当原始图像的局部背景灰度差别大于等于原始图像相应位置上的临界可见偏差时，变量F累计增加1分；
当该位置上原始图像与增强图像的相关系数与增强图像局部背景灰度差别的乘积大于等于增强图像相应位置的临界可见偏差时，变量H累计增加1分；
分数H/F即为增强算法的得分。
所述局部背景灰度差别的计算方法为：
mg0＝max{mg1,mg2,mg3,mg4}
其中，
(i0,j0)代表原始图像或者增强图像的像素位置；mg0代表像素点(i0,j0)的局部背景灰度差别，模板Gk的原点(0,0)在矩阵的左上角，I为原始图像或增强后图像。
所述局部背景灰度均值的计算方法为：
其中，
(i0,j0)代表原始图像或者增强图像的像素位置，I为原始图像或增强后图像，模板B为5*5的矩形邻域，并且原点(0,0)在矩阵的左上角。
所述原始图像与增强图像的相关系数的计算方法为：
其中，cov()代表协方差，Clocal,Dlocal分别代表原始图像与增强后图像以(i0,j0)为中心的局部区域，cov(Clocal,Clocal)代表Clocal的方差，cov(Dlocal,Dlocal)代表Dlocal的方差。
所述临界可见偏差的计算方法为：
Jnd＝max{f1(bg,mg),f2(bg)}
其中，
f1(bg,mg)＝mg·α(bg)+β(bg)
α(bg)＝bg·0.0001+0.115
β(bg)＝λ-bg·0.01
f1是人眼视觉特性中对比度掩盖效应带来的最小可分辨灰度差别的阈值；f2是人眼视觉特性中亮度掩盖效应带来的人眼最小可分辨灰度差别的阈值；f1与f2中的最大值代表了临界可见偏差参数；T0,T1,γ,λ与具体的显示器、观测环境有关，在实际应用中可事先测定，bg是局部背景灰度均值，mg代表局部背景灰度差别。
所述分数H/F：当H/F大于1时，说明通过图像增强算法处理，人眼可以分辨更多的图像细节；当H/F小于1时，可能是因为通过图像增强算法处理，人眼可以分辨的图像细节减少了，另外也可能是图像增强算法引入了很多的噪声，降低了原始图像与增强图像的相关系数；当H/F趋近1时，说明图像增强算法并没有带来很多人眼可感知的图像变化。
本发明具有以下优点及有益效果：
1.相比于现有的其他算法，本发明参考了人眼视觉特性对图像质量评价的影响；
2.该方法并非简单地对增强前后两幅图像进行打分，再对打分结果进行比较得到增强算法性能的评分，而是在对增强算法性能评分的同时，考虑前后两幅图像的相关性，当两幅图像不相关时，最终的评分应该是很低的。
附图说明
图1为本发明算法处理的流程图。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
参阅图1，一种基于人眼视觉特性的图像增强评价方法，包含以下步骤：
1.输入原始图像C与图像增强算法处理后的图像D，逐点扫描两幅图像，针对每一个位置(i0,j0)的像素使用公式(1)分别计算原始图像与增强图像的局部背景灰度差别；使用公式(2)分别计算原始图像与增强图像的局部背景灰度均值；使用公式(3)计算原始图像与增强图像的局部相关系数，局部大小选择以(i0,j0)为中心的3*3的矩形区域。
2.扫描原始图像与增强后图像，针对每一像素(i0,j0)计算该位置上的局部背景灰度差别。计算背景灰度差别选取的邻域大小为5*5。到目前为止，以上的计算过程中涉及三种用处不同的邻域大小：一个是用来计算局部背景灰度差别和局部背景灰度均值的邻域，该区域大小与所使用的模板大小一致都是5*5的；第二个是用来计算原始图像与增强图像的局部相关系数的邻域，该区域大小可以选择为3*3；第三个是用于计算局部背景灰度差别的邻域，选取为5*5大小。
3.将以上步骤中计算得到的局部背景灰度差别与局部背景灰度均值作为公式(4)的输入计算原始图像与增强图像在像素位置(i0,j0)上的临界可见偏差。其中，公式(4)中的某些参数具体值可以采用T0＝17,T1＝127,γ＝3/127,λ＝1/2，这些具体指来自于事先测定。
4.初始化变量F和变量H，初始值设置为零。扫描原始图像与增强图像，针对每一个像素判断：
当原始图像的灰度差别大于等于相应的临界可见偏差时，变量F累计增加1分；
当相关系数与增强图像的灰度差别的乘积大于等于相应的临界可见偏差时，变量H累计增加1分。
扫描结束后，分数H/F即为图像增强算法的最终得分。