[0048] 图3是直接使用线性拉伸处理的结果；
[0049] 图4是使用2%截断拉伸处理的结果；
[0050] 图5是使用本发明方法处理的结果。
【具体实施方式】
[0051] 以下结合附图，针对具体问题的实施过程做进一步说明。
[0052] 本发明的直方图紧致性变换，是一种以压缩图像灰度分布范围为目的的直方图变 换方法，它将包含像素个数很少的灰度级合并到图像其它灰度级上，除去了图像最大值与 最小值之间零概率灰度级和小概率灰度级，使得变换后的图像在其最大值和最小值之间的 所有灰度级上是连续的，并且每个灰度级上的像素个数都不小于基准个数，从而能有效地 压缩图像灰度的分布区间。
[0053] 对于本发明，依据输入图像的灰度分布特征，自动计算紧致性变换过程需要的参 数是具体实施方案的关键步骤，主要涉及前述方案中的步骤1~步骤3。下面，结合附图1 中的流程图，对关键步骤进行说明：
[0054] 步骤1 :确定直方图紧致性变换的初始搜索位置M ;
[0055] 根据直方图紧致性变换原理，起始搜索点的位置可以任意确定。但本发明考虑输 入图像的直方图分布特征，故起始搜索点位置不同，将对压缩结果产生影响。由于图像灰度 均值对图像灰度分布具有代表性，故本发明采用图像灰度的加权均值作为直方图紧致性变 换的起始搜索位置，计算公式如下
[0057] 其中，N是根据输入图像的最大值推算出来的图像量化位数，L1表示第i个灰度 级，Pi表示第i个灰度级的概率。该结果对步骤3产生直接影响。
[0058] 步骤2 :计算图像灰度的理想频度；
[0059] 直方图均衡方法的核心原理是使得变换后的图像灰度值在其量化区间内的每个 灰度级上等概率分布，此时图像具有较好的视觉效果。对于一幅大小为wXh的灰度图像， 当采用N位量化时，每个灰度级等概率出现时每个灰度级上应包含的像素个数为
[0061] 本发明将ζ。作为输入图像的理想频度。
[0062] 步骤3 :计算左基准频度和右基准频度；
[0063] 应用目的不同，一般对直方图压缩程度的要求也不同。故在本发明中引入调节系 数P，用户可以根据具体的应用目的设定。同时，不同的图像，通常具有不同的灰度直方图 分布。因而在直方图紧致性变换过程中，采用与图像直方图分布特征相适应的基准频度具 有更好的适应性。本发明依据直方图紧致性变换原理，将全局性基准频度分解为以起始搜 索点为中心的左基准频度G和右基准频度ζ R。左基准频度与直方图左端分布形态相关， 右基准频度与直方图右端分布形态相关。左基准频度用于直方图紧致性变换时在搜索起始 位置M的左侧搜索时的判断基准，右基准频度用于直方图紧致性变换时在搜索起始位置M 的右侧搜索时的判断基准。它们的计算公式如下： CN 105184759 A 说明书 5/5 页
[0065] 其中，LMax、LMin、Lm分别对应图像的最大值、最小值和均值，N为输入图像的量化位 数，P为用户调节系数（通常在0.5~0.8范围内取值即可满足一般的应用要求）。P越 大，压缩效果越接近等概率直方图；P越小，压缩效果越接近直接线性变换结果。
[0066] 效果展示
[0067] 图2所示为一幅灰度图像，灰度分布区间[37, 255]。图3所示为直接采用线性拉伸 变换之后的结果图像。从该图像上很容易观察，图像的主观视觉效果并没有得到有效改善。 图4是利用截断拉伸算法采用2%截断拉伸增强后的结果。由于该算法对处于中部灰度级 不进行处理，仅对图像直方图上两端的灰度级进行处理，不但没有有效提升图像视觉效果， 反而造成低灰度区域对比度的恶化，丢失了许多图像细节。图5是采用本发明算法（P = 0.6)公开的方法处理的结果。该结果不但有效地增强了图像的对比度，同时暗区的细节纹 理也得到较好的保持。将四幅图综合对比，很显然采用本发明公开的算法处理后的效果是 最好。它为图像对比度、清晰度自动增强处理，提供了一种有效的方法。
[0068] 本发明同样适合于彩色（多通道）图像的增强处理。与灰度图像处理不同的是， 需要将彩色图像首先分解成单一的通道。如RGB三通道图像，首先将图像分解为R、G、B三 个独立通道图像，每个通道图像即为一幅灰度图像，然后对每个通道实施上述增强处理，然 后将增强后的图像再次合成RGB三通道图像。
【主权项】
1. 一种基于直方图紧致性变换的图像自适应增强方法，其特征在于：根据输入图像的 特点计算变换参数，结合用户调节系数对输入图像进行直方图紧致性变换，然后利用线性 拉伸算法对输入图像重新赋值，输出增强图像； 所述直方图紧致性变换如下：设定一个频度阈值作为基准频度，与直方图上的每个灰 度级按顺序进行比较；除特殊情况外，对频度小于基准频度的灰度级，与该灰度级的下一个 或多个灰度级合并，直至当前所有合并灰度级的累积频度不小于基准频度；所述特殊情况 是指：合并后的累积频度与基准频度之差大于原频度与基准频度之差，则不作合并；最终 使得合并后的直方图在其灰度最大值和最小值之间所有灰度级的频度均接近或不小于基 准频度，且灰度连续、不间断； 所述变换参数包括直方图紧致性变换的起始位置和所述频度阈值。2. 根据权利要求1所述的基于直方图紧致性变换的图像自适应增强方法，其特征在 于：所述起始位置位于直方图的中部，以该起始位置将直方图划分为两个部分，并自该起始 位置按照正向、反向依次判别并处理这两个部分的每个灰度级。3. 根据权利要求2所述的基于直方图紧致性变换的图像自适应增强方法，其特征在 于：根据直方图两个部分的特点分别计算各自的基准频度，即所述基准频度分为左基准频 度和右基准频度，作相应的判别并即时处理。4. 根据权利要求1所述的基于直方图紧致性变换的图像自适应增强方法，其特征在 于：所述用户调节系数用于参与计算基准频度，以更符合不同应用目的需要。5. 根据权利要求1所述的基于直方图紧致性变换的图像自适应增强方法，其特征在 于，具体步骤如下： 步骤1 :统计输入图像f(x，y)的灰度直方图及其最大值LMax、最小值LMin和均值LM，根 据最大值确定图像的量化位数N，并采用下式计算初始搜索起点M其中，L1表示第i个灰度级，p i表示第i个灰度级的频度，int表示取整； 步骤2 :根据图像的宽度w和高度h，以及图像的量化位数N，计算图像灰度级的理想频 度步骤3 :根据图像的统计特征、理想频度及调节系数P，利用如下公式计算左基准频度 ζ JP右基准频度ζ R其中，P为用户调节系数，P >〇,取0.5~0.8; 步骤4 :进行直方图紧致性变换，同时确定合并灰度级的映射集合和压缩后图像灰度 的最小值L'Min和最大值L' Max; 步骤5 :计算线性拉伸变换系数k，并将截距置零； 步骤6 :利用线性拉伸变换关系式g(x，y) = k[f (X，y)_LMin]和灰度映射关系，对输入图 像f(x, y)重新赋值，并输出给增强图像g(x, y)。
【专利摘要】本发明提供一种基于直方图紧致性变换原理的图像增强方法，其步骤包括：步骤1：统计输入图像的灰度频度直方图，并计算图像直方图的特征值；步骤2：计算搜索起点及左右基准频度；步骤3：直方图紧致性变换；步骤4：建立灰度映射关系；步骤5：输入图像灰度重构。由于本方法顾及了输入图像的自身特点，因而能有效克服经典线性拉伸变换、直方图截断拉伸等增强方法的不足。实验表明，该方法增强效果明显，且具有良好的适应性和实用性，可广泛适用于各种低对比度图像的增强和高位宽量化图像向低位宽量化图像转化处理，是一种快速、有效的图像自动增强方法。
【IPC分类】G06T5/40
【公开号】CN105184759
【申请号】CN201510607885
【发明人】柳稼航, 陈铁桥, 王一豪, 朱锋
【申请人】中国科学院西安光学精密机械研究所
【公开日】2015年12月23日
【申请日】2015年9月22日