本发明涉及视频图像增强领域，特别是指一种分程度直方图均衡融合的图像增强方法。
背景技术：
直方图均衡(he)方法是一种经典的图像增强算法，它通过统计学原理，使得增强后图像的直方图趋于均匀分布，从而提高图像的对比度。然而，原图直方图的过小数据会导致增强后的细节丢失，过大数据会导致过度增强，不利于人眼观察。
因此，有人提出基于块的对比度受限直方图均衡算法(clahe)，此方法能有效的防止过增强现象，并对局部进行合理增强。但是由于基于子块统计直方图，并由双线性插值方式得到最终结果图，难免会出现块效应现象，同时由于分块局部处理，也产生了大量的计算复杂度。本发明提出一种不同程度直方图均衡融合的图像增强方法，对原图进行不同程度的全局处理再通过融合方式得到最终增强图，该方法具有局部增强的效果，符合人眼视觉感知，同时大大降低了局部处理的计算复杂度。
技术实现要素：
本发明的主要目的在于克服现有技术中的上述缺陷，提出一种分程度直方图均衡融合的图像增强算法，用来克服局部增强计算复杂度大、块效应现象明显等问题，使得局部增强效果简单有效。
本发明采用如下技术方案：
一种分程度直方图均衡融合的图像增强方法，包含如下步骤：
取输入图像的亮度图i，统计灰度级直方图并将所述灰度级直方图均值归一化得到hist；
分别取gamma1和gamma2对归一化得到的hist进行校正处理得到校正后的直方图hist1和hist2；
将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2；
计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o。
具体地，所述取输入图像的亮度图i，统计灰度级直方图并将所述灰度级直方图均值归一化得到hist中，所述取输入图像的亮度图i，具体为：
所述亮度图i由rgb转化为非rgb颜色空间后提取的亮度通道的灰度图。
具体地，所述取输入图像的亮度图i，统计灰度级直方图并将所述灰度级直方图均值归一化得到hist中，具体为：
hist＝hist_original/sum(hist_original)*l；
其中hist_original为原始统计的直方图，l为直方图非零灰度级数目。
具体地，所述分别取gamma1和gamma2对归一化得到的hist进行校正处理得到校正后的直方图hist1和hist2中，具体为：
所述gamma1取值在0～0.5，gamma2取值在0.5～1之间，hist1接近均匀直方图，hist2接近原始直方图。
具体地，所述将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2，具体为：
所述t1根据直方图均衡原理得到，所述t2根据直方图均衡原理得到。
具体地，所述将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2，具体为：
所述增强图i1由原图经过映射曲线t1映射得到，即i1＝t1(i)；
所述增强图i2由原图经过映射曲线t2映射得到，即i2＝t2(i)。
具体地，所述计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o，具体为：
直方图hist的灰度级非零占空比w1：w1＝l/drange；
其中l为直方图非零灰度级数目，drange表示整个直方图的动态范围灰度级数目。
具体地，所述计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o，具体为：
所述最终增强图像o的融合公式为：o＝w1*i1 w2*i2；
其中w2＝1-w1。
由上述对本发明的描述可知，与现有技术相比，本发明具有如下有益效果：
(1)本发明提出的一种分程度的直方图均衡融合增强方法针对经典方法的缺陷，分别取不同的校正系数对归一化得到的hist进行校正处理得到校正后的直方图；将经过校正处理后的直方图进行均衡处理，得到对应的映射曲线，并将原图灰度值分别按这两条映射曲线进行映射，得到不同程度的增强图像；计算均值归一化得到的直方图的灰度级非零占空比，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将不同程度增强的图像进行融合，得到最终增强图像o，即是通过全局不同程度的增强再融合的方式，达到局部细节纹理增强的效果，并降低了计算的复杂度，使得方法更具实时性，适合产品化。
附图说明
图1为本发明提供的非线性变换的图像增强方法实现的整体流程图；
图2为本发明方法实施的效果对比图，其中，图2(a)为输入图像，图2(b)为结果图像。
以下结合附图和具体实施例对本发明作进一步详述。
具体实施方式
如图1所示流程，本发明为一种不同程度直方图均衡融合的图像增强方法，包括：统计原图亮度成分直方图；对直方图进行两次不同程度的gamma校正；对两个校正的后的直方图进行均衡增强处理；将不同增强程度的两幅进行融合处理。具体实施方法包括如下步骤：
s1：取输入图像的亮度图i，统计灰度级直方图并将所述灰度级直方图均值归一化得到hist；
具体地，所述取输入图像的亮度图i，统计灰度级直方图并将所述灰度级直方图均值归一化得到hist中，所述取输入图像的亮度图i，具体为：
所述亮度图i由rgb转化为非rgb颜色空间后提取的亮度通道的灰度图。
取待处理图像的亮度图i，统计其灰度级直方图hist0，并将其均值归一化为hist，实施方法为：
hist＝hist0/sum(hist0)*l；
其中hist_original为原始统计的直方图，l为直方图中非零灰度级数目。
s2：分别取gamma1和gamma2对归一化得到的hist进行校正处理得到校正后的直方图hist1和hist2；
具体地，所述分别取gamma1和gamma2对归一化得到的hist进行校正处理得到校正后的直方图hist1和hist2中，具体为：
所述gamma1取值在0～0.5，gamma2取值在0.5～1之间，hist1接近均匀直方图，hist2接近原始直方图。
分别取范围在0～0.5的gamma1和0.5～1的gamma2对均值归一化后的直方图hist进行校正处理，依次得到校正后的直方图记为hist1和hist2，实施方法为：
hist1＝hist^gamma1；
hist2＝hist^gamma2；
s3：将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2；
具体地，所述将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2，具体为：
所述t1根据直方图均衡原理得到，所述t2根据直方图均衡原理得到。
具体地，所述将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2，具体为：
所述增强图i1由原图经过映射曲线t1映射得到，即i1＝t1(i)；
所述增强图i2由原图经过映射曲线t2映射得到，即i2＝t2(i)。
将校正后的直方图进行均衡处理，得到对应映射曲线t1和t2，并由原图灰度值经过这两条映射曲线映射得到较弱增强图i1和较强增强图i2，实施方法为：
t1＝cumsum(hist1/sum(hist1))*(lmax-lmin 1) lmin；
t2＝cumsum(hist2/sum(hist2))*(lmax-lmin 1) lmin；
i1＝t1(i)；
i2＝t2(i)；
其中cumsum(*)为数组累加操作，lmax为统计直方图范围的最大像素值，lmin为统计直方图范围的最小像素值。
s4：计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o。
具体地，所述计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o，具体为：
直方图hist的灰度级非零占空比w1：w1＝l/drange；
更具体的为：w1＝l/(lmax-lmin 1)；
其中l为直方图非零灰度级数目，drange表示整个直方图的动态范围灰度级数目；lmax为统计直方图范围的最大像素值，lmin为统计直方图范围的最小像素值。
具体地，所述计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o，具体为：
所述最终增强图像o的融合公式为：o＝w1*i1 w2*i2；
其中w2＝1-w1。
本发明实施例实验图像的增强效果见图2，其中图2(a)为原图，图2(b)为增强结果图。从增强前后对比不难看出，本发明提供的是一种不同程度直方图均衡融合的图像增强方法，通过不同程度的全局均衡增强进行融合，来达到局部增强的效果并且大大降低计算复杂度，既避免图像纹理细节丢失，又避免图像过度增强，更符合人眼的视觉特性。
上述仅为本发明的具体实施方式，但本发明的设计构思并不局限于此，凡利用此构思对本发明进行非实质性的改动，均应属于侵犯本发明保护范围的行为。
1.一种分程度直方图均衡融合的图像增强方法，其特征在于，包含如下步骤：
取输入图像的亮度图i，统计灰度级直方图并将所述灰度级直方图均值归一化得到hist；
采用两个不同的参数gamma1和gamma2对归一化得到的hist进行伽马校正处理得到校正后的直方图hist1和hist2；
将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2；
计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o。
2.如权利要求1所述的一种分程度直方图均衡融合的图像增强方法，其特征在于，所述取输入图像的亮度图i，统计灰度级直方图并将所述灰度级直方图均值归一化得到hist中，所述取输入图像的亮度图i，具体为：
所述亮度图i由rgb转化为非rgb颜色空间后提取的亮度通道的灰度图。
3.如权利要求1所述的一种分程度直方图均衡融合的图像增强方法，其特征在于，所述取输入图像的亮度图i，统计灰度级直方图并将所述灰度级直方图均值归一化得到hist中，具体为：
hist＝hist_original/sum(hist_original)*l；
其中hist_original为原始统计的直方图，l为直方图非零灰度级数目。
4.如权利要求1所述的一种分程度直方图均衡融合的图像增强方法，其特征在于，所述采用两个不同的参数gamma1和gamma2对归一化得到的hist进行伽马校正处理得到校正后的直方图hist1和hist2中，具体为：
所述gamma1取值在0～0.5，gamma2取值在0.5～1之间。
5.如权利要求1所述的一种分程度直方图均衡融合的图像增强方法，其特征在于，所述将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2，具体为：
所述t1根据直方图均衡原理得到，所述t2根据直方图均衡原理得到。
6.如权利要求1所述的一种分程度直方图均衡融合的图像增强方法，其特征在于，所述将经过校正处理后的直方图hist1和hist2进行均衡处理，得到对应的映射曲线t1和t2，并将原图灰度值分别按这两条映射曲线进行映射，得到增强图像i1和增强图像i2，具体为：
所述增强图i1由原图经过映射曲线t1映射得到，即i1＝t1(i)；
所述增强图i2由原图经过映射曲线t2映射得到，即i2＝t2(i)。
7.如权利要求1所述的一种分程度直方图均衡融合的图像增强方法，其特征在于，所述计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o，具体为：
直方图hist的灰度级非零占空比w1：w1＝l/drange；
其中l为直方图非零灰度级数目，drange表示整个直方图的动态范围灰度级数目。
8.如权利要求1所述的一种分程度直方图均衡融合的图像增强方法，其特征在于，所述计算均值归一化得到的hist的灰度级非零占空比，记为w1，并以所述灰度级非零占空比作为增强程度较弱图像的权重，将增强图像i1和增强图像i2进行融合，得到最终增强图像o，具体为：
所述最终增强图像o的融合公式为：o＝w1*i1 w2*i2；
其中w2＝1-w1。
技术总结