本发明涉及红外图像增强技术领域，具体地说，本发明涉及一种基于自适应直方图分段的红外图像增强技术方案。
背景技术：
红外成像系统抗干扰能力强，隐蔽性能好，大气穿透能力强，能够适应多种特殊场合。但是由于红外探测器自身特性，如灵敏度高、动态范围大等特性，以及复杂工作环境中的各种噪声干扰，导致红外图像呈现出高背景，低反差的特点。具体表现为成像场景在整个红外成像系统的动态范围中只占了很小一部分，图像对比度差，模糊不清。因此需要对原始红外图像进行图像增强，从而提高红外图像的对比度，改善图像的视觉效果。
目前，主流的红外图像增强方法有局部映射法和全局映射法。二者的本质区别是是否依据局部的图像信息来指导灰度值的映射。一般而言，局部映射法能够取得较好的增强效果，这是由于全局映射法只能根据较为宏观的信息来进行图像的增强，很容易产生过度增强的问题。但是，由于局部映射法会打乱原有灰度值的顺序，而红外图像中像素的明暗代表的是温度的高低，这会对一些温度敏感的应用产生困扰。如果在全局映射法的框架下，提高增强的效果，一直以来都是一个难题。
技术实现要素：
本发明的目的在于提供一种红外图像的增强技术方案，通过自适应的将直方图进行分段，并根据每个分段的属性判断其是属于背景还是目标，然后采取相应的不同的增强策略。经本发明处理后的图像对比度高，细节清晰，质量好。
本发明的技术方案包括一种基于自适应直方图分段的红外图像增强方法，包括以下步骤：
步骤1，将红外图像统计直方图自适应的分割为若干段，对于背景分段进行抑制，分配较小的动态范围，对于目标分段进行增强，分配较大的动态范围；
步骤2，根据直方图每个分段的平均直方图值ahv分配其在增强图像中的动态范围；
步骤3，在所述分配的动态范围内对各分段内的灰度进行灰度映射，得到增强后的红外图像。
而且，步骤1中，所述将红外图像统计直方图自适应的分割为若干段，包括以下过程，
步骤1.1，对输入的红外图像进行直方图统计，得到输入图像的直方图h(x)＝nx，x＝0,1,…,l-1；其中nx表示灰度值为x的灰度在图像中出现的次数，l为图像的灰度级数目；
步骤1.2，剔除直方图h(x)中为0的元素，得到非0直方图hv(t)，设共有t个非零元素，则t＝1,2,…,t，相应数学定义如下，
步骤1.3，剔除非0直方图hv(t)两端a×t个元素，得到新的直方图hr(s)，相应数学定义如下，
其中，a为预设参数，
步骤1.4，用方差参数为σ的高斯滤波器对hr(s)进行滤波，参数σ为预设值；
步骤1.5，以宽度为w的窗口遍历滤波后的hr(s)，当窗口中心的元素小于其它窗口内元素时，记为局部最小值，其中hr(s)的第一个元素和最后一个元素也视为局部最小值，经过窗口的遍历得到各局部最小值对应灰度的集合{s1,s2,…,sm}，m为局部最小值的数目，则直方图被自适应的分为[s1,s2],[s2,s3],…,[sm-1,sm]共计m-1个分段。
而且，步骤2中，所述根据直方图每个分段的平均直方图值ahv分配其在增强图像中的动态范围，包括以下过程，
步骤2.1，计算每个分段的平均直方图值ahv，相应定义如下，
其中，n_subi为第i个分段内的灰度级数，n为图像的总像素个数，min{,}表示取其中的最小值，cdfi为第i个分段内所有灰度的累计统计概率，定义如下，
步骤2.2，根据下式，为每个直方图分段分配其在增强图像中的动态范围d：
其中，dj是第j个分段在增强图像中的动态范围，
而且，th和gcdfi的计算方式如下，
th＝β×n/t,
gcdfi＝(cdfi)γ+0.1×η,
η＝(maxi-s1)/(sm-s1),
其中，β为预设参数，maxi为第i分段中最大的灰度值，η、γ是伽马变换中的参数。
而且，步骤3中，在所述分配的动态范围内对各分段内的灰度进行灰度映射时，映射关系定义如下，
其中，y(x)表示原图中灰度值x在增强后图像中所映射的灰度值。
本发明提供一种基于自适应直方图分段的红外图像增强系统，包括以下模块：
第一模块，用于将红外图像统计直方图自适应的分割为若干段，对于背景分段进行抑制，分配较小的动态范围，对于目标分段进行增强，分配较大的动态范围；
第二模块，用于根据直方图每个分段的平均直方图值ahv分配其在增强图像中的动态范围；
第三模块，用于在所述分配的动态范围内对各分段内的灰度进行灰度映射，得到增强后的红外图像。
而且，第一模块中包括以下单元，
第一单元，用于对输入的红外图像进行直方图统计，得到输入图像的直方图h(x)＝nx，x＝0,1,…,l-1；其中nx表示灰度值为x的灰度在图像中出现的次数，l为图像的灰度级数目；
第二单元，用于剔除直方图h(x)中为0的元素，得到非0直方图hv(t)，设共有t个非零元素，则t＝1,2,…,t，相应数学定义如下，
第三单元，用于剔除非0直方图hv(t)两端a×t个元素，得到新的直方图hr(s)，相应数学定义如下，
其中，a为预设参数，
第四单元，用于用方差参数为σ的高斯滤波器对hr(s)进行滤波，参数σ为预设值；
第五单元，用于以宽度为w的窗口遍历滤波后的hr(s)，当窗口中心的元素小于其它窗口内元素时，记为局部最小值，其中hr(s)的第一个元素和最后一个元素也视为局部最小值，经过窗口的遍历得到各局部最小值对应灰度的集合{s1,s2,…,sm}，m为局部最小值的数目，则直方图被自适应的分为[s1,s2],[s2,s3],…,[sm-1,sm]共计m-1个分段。
而且，第二模块中包括以下单元，
第一单元，用于计算每个分段的平均直方图值ahv，相应定义如下，
其中，n_subi为第i个分段内的灰度级数，n为图像的总像素个数，min{,}表示取其中的最小值，cdfi为第i个分段内所有灰度的累计统计概率，定义如下，
第二单元，用于根据下式，为每个直方图分段分配其在增强图像中的动态范围d：
其中，dj是第j个分段在增强图像中的动态范围，
而且，th和gcdfi的计算方式如下，
th＝β×n/t,
gcdfi＝(cdfi)γ+0.1×η,
η＝(maxi-s1)/(sm-s1),
其中，β为预设参数，maxi为第i分段中最大的灰度值，η、γ是伽马变换中的参数。
而且，第三模块中，在所述分配的动态范围内对各分段内的灰度进行灰度映射时，映射关系定义如下，
其中，y(x)表示原图中灰度值x在增强后图像中所映射的灰度值。
本发明的增强方法通过对直方图的自适应分段，从而能够分析每个分段的具体属性，对于背景分段进行抑制，分配较小的动态范围，对于目标分段进行增强，分配较大的动态范围，实现了有限动态范围下的高效利用，增强后的图像对比度高，视觉效果好。
附图说明
图1为本发明实施例基于自适应直方图分段的红外图像增强方法流程图。
图2为本发明实施例自适应分段的结果示意图。
图3为本发明实施例动态范围分配以及灰度映射结果示意图。
具体实施方式
下面结合附图和实施例对本发明进行进一步说明。
本发明公开了一种基于自适应直方图分段的红外图像增强方法，包括以下步骤：自适应的将红外图像统计直方图分割为若干段；根据直方图每个分段的属性分配其在增强图像中的动态范围；在所述分配的动态范围内对各分段内的灰度进行灰度映射。本发明的增强方法增强对比度高，图像清晰，视觉效果好。
如图1所示，本发明实施例基于边缘提取的红外直方图增强方法包括以下步骤：
步骤1：将红外图像统计直方图自适应的分割为若干段；包含以下具体步骤：
步骤1.1：对输入的红外图像进行直方图统计，得到输入图像的直方图h(x)＝nx，x＝0,1,…,l-1；其中nx表示灰度值为x的灰度在图像中出现的次数，l为图像的灰度级数目，例如14bit的图像，l＝16384；
步骤1.2：剔除直方图h(x)中为0的元素，得到非0直方图hv(t)，假设共有t个非零元素，则t＝1,2,…,t，该步骤的数学定义如下：
步骤1.3：剔除非0直方图hv(t)两端a×t个元素，得到新的直方图hr(s)，该步骤的数学定义如下：
其中，a为预设参数，
步骤1.4：用方差参数为σ的高斯滤波器对hr(s)进行滤波，参数σ为预设值，设为t×0.1％；
步骤1.5：以宽度为w的窗口遍历滤波后的hr(s)，当窗口中心的元素小于其它窗口内元素时，记为局部最小值，其中hr(s)的第一个元素和最后一个元素也视为局部最小值，则经过窗口的遍历得到各局部最小值对应灰度的集合{s1,s2,…,sm}，m为局部最小值的数目，则直方图被自适应的分为[s1,s2],[s2,s3],…,[sm-1,sm]共计m-1个分段。
具体实施时，优选w＝t/50～t/100。
如图2所示，即为m＝7的一个分段结果示意图。直方图被分割为6个“峰”形状的分段，分别为s1,s2,…,s7。
步骤2：根据直方图每个分段的平均直方图值ahv分配其在增强图像中的动态范围；包含以下具体步骤：
步骤2.1：计算每个分段的平均直方图值ahv，其定义如下：
其中，n_subi为第i个分段内的灰度级数，n为图像的总像素个数，min{,}表示取其中的最小值，cdfi为第i个分段内所有灰度的累计统计概率，定义如下：
步骤2.2：根据下式，为每个直方图分段分配其在增强图像中的动态范围d：
其中，dj是第j个分段在增强图像中的动态范围，
th＝β×n/t,
gcdfi＝(cdfi)γ+0.1×η,
η＝(maxi-s1)/(sm-s1),
其中，β为预设参数，maxi为第i分段中最大的灰度值，η、γ是伽马变换中的参数，γ取0.55，β取2-6之间。
步骤3：在所述分配的动态范围内对各分段内的灰度进行灰度映射，得到增强后的红外图像。具体映射关系定义如下：
其中y(x)表示原图中灰度值x在增强后图像中所映射的灰度值。
具体实施时，本发明所提供方法可基于软件技术实现自动运行流程，也可采用模块化方式实现相应系统。
本发明提供一种基于自适应直方图分段的红外图像增强系统，包括以下模块：
第一模块，用于将红外图像统计直方图自适应的分割为若干段，对于背景分段进行抑制，分配较小的动态范围，对于目标分段进行增强，分配较大的动态范围；
第二模块，用于根据直方图每个分段的平均直方图值ahv分配其在增强图像中的动态范围；
第三模块，用于在所述分配的动态范围内对各分段内的灰度进行灰度映射，得到增强后的红外图像。
各模块具体实现可参见相应步骤，本发明不予赘述。
如图3所示，即为动态范围分配以及灰度映射结果示意图。可以看到形状最陡峭的分段被判定为背景分段，因此其分配的动态范围没有得到扩大，其余分段的动态范围均得到了不同程度的扩大。
经实验对比，扫描型红外系统输出的、分辨率大小为1024×1280的红外图像。该红外图像未经过图像增强处理，可以看到，图像对比度差，图像细节难以分辨。
经过本发明增强方法处理得到的红外图像。可以看到，图像对比度得到了增强，图像背景均匀。
经过传统直方图均衡化增强方法(he)处理得到的红外图像。可以看到，虽然图像对比度得到了增强，但是背景的噪声也增强了。
经过直方图均匀四分段增强方法(ehpe)处理得到的红外图像。可以看到，图像对比度没有得到有效的增强。
需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明并不限于具体实施方式中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，同样属于本发明保护的范围。