本发明涉及图像处理技术领域，更具体地，涉及一种颜色保持和灰度级压缩的CLHAE(Contrast Limited Adaptive Histogram Equalization，限制对比度自适应直方图均衡)图像增强方法。
背景技术：
随着社会的发展，汽车作为交通工具在生活中越发显得不可或缺，相应的车载产品也不断涌现，其中以车载图像产品最为突出。在实际应用中，车载产品采集的图像往往因为设备、环境亮度、图像压缩方法等等，导致获得的图像的对比度不高，清晰度低，不利于后续的处理。为了改善这些缺陷，需要对这类图像进行图像增强。
传统的图像增强方法分为两类：空域相关法和频域相关法。其中空域相关法一般指对图像像素直接进行处理；常用的方法有直方图均衡化，直方图规定化(匹配)，伽马调节等等。这些方法提出得较早，且最为成熟，都是针对全局图像进行图像增强的。对于整幅图像而言，图像的灰度级得到了扩展，但是对于灰度级分布较“均匀”的原始图像，这些方法处理的图像亮度较平滑，图像显得“泛白”，即对比度没有得到提升。因此上述图像增强方法均无法获得清晰满意的图像增强效果。
技术实现要素：
本发明的目的在于克服现有技术存在的上述缺陷，提供一种基于灰度级压缩的CLHAE图像增强方法。
为实现上述目的，本发明的技术方案如下：
一种基于灰度级压缩的CLHAE图像增强方法，包括以下步骤：
步骤S101：将输入的当前RGB颜色空间图像转换为YUV颜色空间图像，并从YUV颜色空间图像中提取出图像亮度Y的每个初始像素值；
步骤S102：对图像亮度Y进行均匀且不重叠地分块；
步骤S103：对图像亮度Y中的每个块进行灰度级压缩；
步骤S104：对图像亮度Y中的每个块进行直方图统计；
步骤S105：对图像亮度Y中的每个块进行限制对比度；
步骤S106：对图像亮度Y中的每个块进行直方图均衡，得到每个块的灰度级映射表；
步骤S107：根据灰度级映射表，对图像亮度Y进行双线性插值算法计算，得到图像亮度Y的每个调整后像素值；
步骤S108：计算图像亮度Y的每个初始像素值与其调整后像素值之间的颜色换算比例；
步骤S109：将颜色换算比例应用到当前RGB颜色空间图像中，输出新的RGB颜色空间图像。
进一步地，步骤S105中，通过设置限制对比度的阈值，对图像亮度Y的每个块进行限制对比度。
进一步地，通过进行限制对比度，得到每个块中超出阈值的总像素数，根据超出阈值的总像素数，得到每个块中每个灰度级的预分配像素数，并分配至每个块中。
进一步地，按照预分配像素数进行分配时，如第一次分配未将超出阈值的总像素数分配完时，则以第一次分配的数据为基础，进行第二次分配。
进一步地，步骤S106中，对完成限制对比度后的每个块的每个灰度级计算得到累进直方图，并根据累进直方图进行直方图均衡。
进一步地，步骤S102中，将图像亮度Y均匀且不重叠地分为四个块。
进一步地，步骤S107中，先通过以过图像亮度Y中每相邻两个块中点的连线为划分线，将图像亮度Y划分为九个小块，并按照各小块所在位置，将各小块分别定义为角区域小块、边区域小块和中心区域小块，然后进行双线性插值算法计算，得到图像亮度Y的每个调整后像素值。
进一步地，针对某角区域小块中的各像素，通过直接查该角区域小块所在的块对应的灰度级映射表，获得该角区域小块中各像素的调整后像素值；针对某边区域小块中的各像素，先查该边区域小块所在的两个块对应的两个灰度级映射表，获得该边区域小块中各像素的对应两个虚拟定点的像素值，然后采用线性插值算法，获得该边区域小块中各像素的调整后像素值；针对中心区域小块中的各像素，先查中心区域小块所在的四个块对应的四个灰度级映射表，获得中心区域小块中各像素的对应四个虚拟定点的像素值，然后采用双线性插值算法，获得中心区域小块中各像素的调整后像素值。
进一步地，所述颜色换算比例是调整后像素值与初始像素值的比值。
从上述技术方案可以看出，本发明通过对当前图像的颜色空间进行转换，针对图像亮度Y进行图像处理，再通过比例换算计算得到新像素的值，可防止伪彩色出现或者图像原R、G、B三个通道颜色比例的改变；同时，通过对图像亮度Y的灰度级进行压缩，对灰度级进行量化统计，可减少存储空间，提高图像处理效率；再次，利用CLAHE算法可实现对图像的增强处理，提高图像的动态范围；因此，本发明实现了在保持原图像颜色比例的基础上，有效地提高了处理效率，并且扩展了图像的动态范围，增强了图像的对比度，从而实现快速且高效的图像增强效果，较之传统的空域图像增强方法有着更好的适应性。
附图说明
图1是本发明的一种基于灰度级压缩的CLHAE图像增强方法流程图；
图2是本发明一较佳实施例中对图像亮度Y分块进行直方图统计的示意图；
图3是本发明一较佳实施例中对图像亮度Y每个块的灰度级进行限制对比度的示意图；
图4是本发明一较佳实施例中对图像亮度Y的像素值进行双线性插值的示意图。
具体实施方式
本发明通过对RGB空间的原始图像进行YUV空间转换，并且对亮度空间Y(Y空间)的灰度级进行压缩，以及使用CLAHE(Contrast Limited Adaptive Histogram Equalization，限制对比度自适应直方图均衡)算法对低对比度图像作图像增强处理，从而实现高效快速的图像增强效果。
下面结合附图，对本发明的具体实施方式作进一步的详细说明。
需要说明的是，下述的具体实施方式，仅为了清楚地解释本发明，因此，应避免以此作为对本发明的限定来加以理解。
在以下本发明的具体实施方式中，请参阅图1，图1是本发明的一种基于灰度级压缩的CLHAE图像增强方法流程图。如图1所示，本发明的一种基于灰度级压缩的CLHAE图像增强方法，包括以下步骤：
步骤S101：将输入的当前RGB颜色空间图像转换为YUV颜色空间图像，并从YUV颜色空间图像中提取出图像亮度Y的每个初始像素值。
首先输入一幅当前图像，其为RGB颜色空间图像，定义其画幅大小为XSize×YSize。由于直接在RGB颜色空间做图像处理，会改变R通道、G通道、B通道三者之间的颜色比例关系，最后输出的图像可能会出现伪彩色。为避免这种情况，因此对当前图像进行RGB颜色空间到YUV颜色空间的转换，并从YUV颜色空间图像中提取出图像亮度Y的每个初始像素值，记为Y_Old，从而得到一幅反映图像亮度Y的灰度图。
在图像亮度Y上做图像处理，具有以下优点：
1)将常规的对R、G、B三个通道做处理，变为对图像亮度Y做处理，使得图像处理的效率得到提高；
2)由于图像亮度Y只代表亮度，不涉及图像颜色，因而更好地保护了R、G、B三个通道之间的颜色比例，从而可防止伪彩色的出现。
步骤S102：对图像亮度Y进行均匀且不重叠地分块。
然后，对上述得到的Y空间的图像进行均匀的而且不重叠的分块，图2示出了这个过程的实现。其中，定义在坐标X方向上的分块数为ClipGrid_X(本实施例中假设块数为2)，每个块的宽度为XGridSize；定义在Y方向的分块数为ClipGrid_Y(本实施例中假设块数为2，共计4个块)，每个块的高度为YGridSize。
步骤S103：对图像亮度Y中的每个块进行灰度级压缩。
考虑到在实际的车载图像产品中，硬件资源的分配非常重要，所以需要尽量减少对图像亮度Y处理时所占据的硬件空间。
如果记原始图像的位宽为Orig_BitWidth，按照现有的算法，需要统计2Orig_BitWidth-1个灰度级，则在后续的算法步骤中共需要ClipGrid_X*ClipGrid_Y*(2Orig_BitWidth-1)*2大小的空间，需要消耗很大的硬件资源空间，这在硬件实施中是不可取的。
因此，需要对每个块的灰度级进行压缩。记需要压缩的位宽为Comp_BitWidth，记最后得到的最大灰度级级数为Graylevel_Max，则实际灰度级级数为：
Real_Graylevel＝[0,Graylevel_Max]
步骤S104：对图像亮度Y中的每个块进行直方图统计。
接着，可以借助CLAHE算法，对当前的低对比度图像进行图像增强处理。
首先，对图像亮度Y的每个块进行直方图统计，图2示出了这个过程的实现。在本实施例中，记每个块的直方图统计结果为：
Hist[i],(i∈[0,Graylevel_Max])
步骤S105：对图像亮度Y中的每个块进行限制对比度。
接着，对已经完成直方图统计的每个块进行限制对比度处理，图3示出了这个过程的实现。在本实施例中，先设置限制对比度的阈值，记限制对比度的阈值为ClipThreshold，对上述统计出的直方图Hist(i)进行限制对比度，得到每个块中超出阈值的总像素数ExcessPixels：
然后，根据求出的ExcessPixels，得到每个块中每个灰度级的预分配像素数ProcPixels：
ProcPixels＝ExcessPixels/Graylevel_Max
接下来，需要将预分配像素数进一步分配至每个块中。
本实施例总共需要对直方图Hist(i)进行两次限制对比度处理，即需要对每个块中每个灰度级的像素数进行两次分配。
其中，进行第一次分配(第一次限制对比度)的实施过程的公式如下：
其中，Hist(i)'代表第一次限制对比度时每个块的直方图统计结果(分配结果)。
当按照预分配像素数进行第一次限制对比度(第一次分配)之后，若每个块中超出阈值的总像素数ExcessPixels未分配完，还有剩余时，则以第一次分配的数据为基础，对每个块中每个灰度级的像素数进行第二次限制对比度(第二次分配)；将剩余像素点数记为Res_ExcessPixels，具体实施公式如下：
Hist(i)″＝Hist(i)'+Res_ExcessPixels/Graylevel_Max
其中，Hist(i)″代表第二次限制对比度时每个块的直方图统计结果(分配结果)。
步骤S106：对图像亮度Y中的每个块进行直方图均衡，得到每个块的灰度级映射表。
接着，对已经完成限制对比度的经过重新分配像素数后的每个块的每个灰度级计算累进直方图CDFn，具体实现过程是这样的：
然后，根据累进直方图对每个块进行直方图均衡，并输出经过直方图均衡后的每个块的灰度级映射表Map，具体实现过程是这样的：
其中，每个块的灰度级映射表分别为：
Mapij(i∈ClipGrid_X,j∈ClipGrid_Y)
步骤S107：根据灰度级映射表，对图像亮度Y进行双线性插值算法计算，得到图像亮度Y的每个调整后像素值。
接下来，对直方图均衡后的图像亮度Y进行双线性插值算法计算，图4示出了这个过程的实现。双线性插值算法的目的在于消除由于分块产生的块效应，本实施例中需要用到由步骤S106得到的当前像素位置周围的四个灰度级映射表，假设分别为：Map11,Map12,Map21,Map22。
图4示出，Y空间的图像原被分为4个不重叠区域(4个块)。可对这4个不重叠区域再取各自中点，以过每相邻两个块中点的连线为划分线，按图示连接划分为9个不重叠的小区域(9个小块)，记编号1～9。其中，按照各小块所在位置，将编号1，3，7，9所在小块定义为角区域小块；将编号2，4，6，8所在小块定义为边区域小块；而将编号5定义为中心区域小块，然后再进行双线性插值算法计算。
对图像亮度Y进行双线性插值可分为三个部分实现：
针对角区域小块1，3，7，9中的各像素，可采用各角区域小块所在的块对应的灰度级映射表，直接查表获得各角区域小块中各像素的调整后像素值Y。比如A点，可带入YA到Map11，直接获得YA′。
针对各边区域小块2，4，6，8中的各像素，先查各边区域小块所在的两个块对应的两个灰度级映射表，获得各边区域小块中各像素的对应两个虚拟定点的像素值，然后采用线性插值算法，获得该边区域小块中各像素的调整后像素值。比如B点，查Map11获得点查Map12获得点以两个点作为B点的2个虚拟定点，然后通过线性插值，获得B点的调整后像素值YB′。其线性插值公式为：
其中，W代表边区域小块2的X向(横向)坐标宽度，x为B点的X向坐标值。
针对中心区域小块中的各像素，先查中心区域小块所在的四个块对应的四个灰度级映射表，获得中心区域小块中各像素的对应四个虚拟定点的像素值，然后采用双线性插值算法，获得中心区域小块中各像素的调整后像素值。比如C点，分别查表Map11,Map12,Map21,Map22，获得以四个点作为C点的4个虚拟定点，通过双线性插值，获得C点的调整后像素值Yc′。其双线性插值公式为：
Yc'＝[Map11×(W-x)×(H-y)+Map12×x×(H-y)+Map21×(W-x)×y+Map22×x×y]/W×H
其中，W代表中心区域小块5的X向(横向)坐标宽度，H代表中心区域小块5的Y向(纵向)坐标高度，x为C点的X向坐标值，y为C点的Y向坐标值。
根据上述计算，即可得到图像亮度Y每个像素经过双线性插值之后的新的像素值(调整后像素值)Y_New(例如YA'、YB'、Yc')。
步骤S108：计算图像亮度Y的每个初始像素值与其调整后像素值之间的颜色换算比例。
为了保持原图像(当前图像)的颜色不出现伪彩色，应该保持原图像的R通道、G通道、B通道的颜色比例关系，这里的颜色通道的换算比例(颜色换算比例)K为：
即颜色换算比例K是调整后像素值Y_New与初始像素值Y_Old的比值。
步骤S109：将颜色换算比例应用到当前RGB颜色空间图像中，输出新的RGB颜色空间图像。
根据上面得到的颜色换算比例K，可以得到本实施例的结果图像(新的RGB颜色空间图像)。得到原图像R、G、B三个通道的新像素值的实现过程是：
其中，R、G、B是原图像三个通道的像素值，R'、G'、B'是输出图像(新的RGB颜色空间图像)三个通道的新的像素值。
综上所述，本发明通过对当前图像的颜色空间进行转换，针对图像亮度Y进行图像处理，再通过比例换算计算得到新像素的值，可防止伪彩色出现或者图像原R、G、B三个通道颜色比例的改变；同时，通过对图像亮度Y的灰度级进行压缩，对灰度级进行量化统计，可减少存储空间，提高图像处理效率；再次，利用CLAHE算法可实现对图像的增强处理，提高图像的动态范围；因此，本发明实现了在保持原图像颜色比例的基础上，有效地提高了处理效率，并且扩展了图像的动态范围，增强了图像的对比度，从而实现快速且高效的图像增强效果，较之传统的空域图像增强方法有着更好的适应性。
以上所述的仅为本发明的优选实施例，所述实施例并非用以限制本发明的专利保护范围，因此凡是运用本发明的说明书及附图内容所作的等同结构变化，同理均应包含在本发明的保护范围内。