基于局部直方图均衡的灰度图像拟合增强方法
【技术领域】
[0001] 本发明涉及一种可应用于各种灰度图像，包括动态范围较高、灰度级位数等于或 大于14位的红外图像处理方法。具体涉及一种改进型的直方图增强方法，以实现对分辨率 低、对比度低、信噪比低和视觉效果模糊的灰度图像进行视觉效果提升的技术方案。
【背景技术】
[0002] 图像获取指的是将物体成像的过程，将模拟图像转换成数字图像的技术，其获取 的图像图像质量由于外界场景的气象因素与探测器成像原理的限制，存在对比度低、分辨 率低、信噪比低和视觉效果等不足。现有技术针对雾天降质图像对比度低，细节模糊的特 点，提出了基于分子块重叠直方图均衡算法0SHE的雾天图像清晰化方法。P0SHE是一种常用 的图像增强方法。P0SHE算法的雾天图像清晰化方法是在输入图像上定义矩形子块，利用该 子块图像的直方图信息对子块所有的像素进行均衡化，将子块以一定的步长移动，直至遍 历输入图像的所有像素。该方法虽然在对雾天降质图像处理时，可以达到局部细节清晰化 和全局对比度增强的效果，但P0SHE算法对红外图像进行增强处理时存在过增强等问题。由 于该算法存在过增强和块效应问题，不利于目标的识别和跟踪。P0SHE方法是在部分重叠子 块中作均衡处理，各个子块的灰度分布不同，在子块边缘会产生明显的灰度不连续现象，即 方块效应，简称块效应。如果图像中存在大面积的背景区域，小子块P0SHE处理会将这些背 景区域分成小的子块进均衡，这自然没考虑到该区域的整体灰度分布。
[0003] 图像增强的首要目标是处理图像，使其比原始图像更适合于特定应用，比如目标 识别、跟踪等。基于直方图均衡(HE)的图像增强方法以其简单快速的特点倍受青睐，全局直 方图均衡(GHE)根据图像直方图计算灰度变换函数，对整幅图像像素作单一映射，没有考虑 图像局部信息，造成部分高频信息丢失。GHE由于过度拉伸图像对比度产生过增强现象，即 累加直方图曲线过于陡峭。基于局部直方图均衡(LHE)算法根据图像局部特性计算灰度变 换函数，克服GHE难以适应局部度分布的缺陷，对比度增强效果较好，它的最初原型是Tom等 人于1982年提出的自适应直方图均衡算法(AHE )，过高的时间复杂度限制AHE算法的应用范 围；Kim等人于2001年提出的子块部分重叠直方图均衡P0SHE算法，具有良好对比度增强效 果，较AHE算法计算量明显减少，但是该算法存在明显块效应和局部过增强问题。采用P0SHE 算法对图像增强处理后，块效应程度、增强效果与子块大小和步长有着密切的关系。然而， 针对不同场景图像，要达到同样的增强效果和较弱的块效应，子块大小和步长的取法往往 不一样。在块大小一致的情况下，步长越大，块效应越突出；在步一致的情况下，块越小，块 效应越突出。大子块、大步长处理只存在轻微的块效应，小子快、小步长则会带来严重的块 效应。较小的步长可以有效平滑块效应，但如果步长太小，会大大增加算法时间复杂度，然 而块效应抑制效果不再明显提高。
[0004]直方图增强技术是最常用的图像增强方式之一。经典直方图均衡化的工作原理是 通过对输入图像直方图统计和累积，映射输入图像的灰度值，使输出图像在灰度区间内均 匀分布，但该方法存在明显的不足，如会过度增强背景和噪声，降低目标对比度，细节易丢 失。自适应直方图均衡化对利用滑动窗口局部区域进行直方图均衡化处理，避免了图像其 它区域的影响，但是该方法会产生块效应。对比度受限的自适应直方图均衡通过局部处理 的方法，限制局部平坦区域的对比度增强，再通过双线性插值的方法消除边界效应，但该方 法时间复杂度过高。子块部分重叠直方图均衡通过子块移动策略来使部分子块重叠并对每 个子块进行HE处理，再通过加权求和的方法来消除噪声和增大对比度，该方法在不降低局 部细节的条件下，虽然大大减少了计算量，但是该方法存在明显的块效应和过增强问题。传 统的灰度图像直方图均衡化算法，均衡处理不能很好的凸显细节，而固定阈值β偏大，又会 出现过增强现象，场景中的事物表现的不自然。小子块处理旨在锐化降质图像中模糊的边 缘，β应当较大，如果过度限制子块的单个灰度级像素的数目，局部对比度无法充分拉伸，不 能充分锐化图像。
[0005] 对比度限制在每层P0SHE处理前，先对其各子块的直方图截取然后均匀分配到各 灰度级上，通过限制其直方图高度降低累加直方图曲线的斜率，即降低对比度增强幅度，从 而限制噪声放大及局部过增强现象。视频仿真结果对比图图的增强效果，增强前图像近处 画面模糊不清，图像远处画面缺乏层次，丢失很多细节。
【发明内容】
[0006] 本发明的目的是针对传统的灰度图像直方图均衡化算法和灰度图像增强技术的 不足之处，提供一种能够充分抑制块效应、消除过增强现象，基于局部直方图均衡拟合的灰 度图像增强方法，以解决用P0SHE算法对红外图像进行增强处理时存在的过增强等问题。
[0007] 本发明的技术解决方案可通过以下步骤实现:一种基于局部直方图均衡的灰度图 像拟合增强方法，其特征在于包括如下步骤:采用基于FPGA架构或DSP架构或FPGA+DSP架构 的实时处理平台，首先对动态范围过宽的灰度图像(如图像数据宽度2 14位的红外图像）， 进行分段线性变换获得适合的动态范围内的灰度图像，再将图像灰度分布区间分隔成两段 乃至多段，然后对灰度图像进行分段线性变换;在每层P0SHE处理前，先对其各子块的直方 图截取然后均匀分配到各灰度级上，并调整图像灰度分布各区间的分割点和变换直线的斜 率，对任意灰度区间进行扩展或压缩;变换结果采用子块部分重叠直方图均衡化处理，将邻 域内的子块变换函数通过加权求和计算得到当前子块的变换函数，使用此变换函数对当前 子块进行直方图均衡化处理;最后对直方图均衡后的灰度图进行非线性拟合处理，然后对 子块部分重叠直方图均衡化处理后的灰度图像进行直方图分布修正，重新分配图像像素 值，实现图像灰度值在多段直方图区间的非线性拉伸，使动态范围内的直方图分布变得更 "均匀"的同时，在各个灰度分布区间具有更优的对比度。
[0008] 本发明具有如下有益效果： 能够充分抑制块效应。本发明利用对灰度图像进行分段线性变换，将图像灰度分布区 间分隔成两段乃至多段，分别对其进行线性变换，通过增加灰度区间分割的段数，以及调整 各区间的分割点和变换直线的斜率，可以对任意灰度区间进行扩展或压缩，从而使得直方 图灰度分布均匀，图像不会过亮或过暗；由于摒弃了对子块大小和步长的过度依赖，在子块 大小和步长一定的情况下，能取得较好的增强效果，充分抑制块效应、消除过增强现象。既 能保证一定的动态范围压缩，又能保留细节信息和颜色保真，从而获得较为令人满意的增 强效果。本发明通过MATLAB仿真对其基于局部直方图均衡的灰度图像拟合增强方法进行了 试验验证。结果表明，该方法不仅可以有效缓解过增强问题，消除经POSHE算法处理后产生 的块效应现象，而且还可以抑制制冷型红外探测器的〃冷反射〃图像。
[0009] 本发明在每层P0SHE处理前，先对其各子块的直方图截取然后均匀分配到各灰度 级上，通过限制其直方图高度降低累加直方图曲线的斜率，即降低对比度增强幅度，从而限 制噪声放大及局部过增强现象。变换结果经子块部分重叠直方图均衡化处理，将邻域内的 子块变换函数通过加权求和计算得到当前子块的变换函数，使用此变换函数对当前子块进 行直方图均衡化处理，对最后对直方图均衡后的灰度图进行非线性拟合处理，对子块部分 重叠直方图均衡化处理后的灰度图像进行直方图分布修正，解决其产生的块效应和过增强 问题。不会使得某些灰度级占有过多的像素，在增强灰度图像对比度和细节信息的同时，能 够消除块效应和过增强现象。在增强红外图像对比度的同时，能够有效保留红外图像的整 体信息，并且具有较好的视觉效果。增强后图像远处细节得到凸显，整体画面非常清晰，具 有真实自然的视觉效果。
[0010] 本发明针对雾天不同退化程度的低照度图像，将对比度限制引入P0SHE方法中，结 合多层融合处理消除了 P0SHE方法中存在的块效应和过增强现象，能够有效增强其对比度 和局部信息，在亮度、对比度和细节方面都能取得较好的视觉效果，具有自适应好、易于并 行实现以及鲁棒性强的优点，因此具有一定的通用性，在视频监控、目标识别和图像制导等 实时图像处理领域有较高的工程实用价值。
【附图说明】
[0011] 图1是本发明基于局部直方图均衡拟合灰度图像的增强流程图。
[0012] 图2是本发明三线段分段线性变换灰度范围示意图。
[0013] 图3是对动态范围大于或等于14位的红外图像进行分段线性变换的示意图。
【具体实施方式】
[0014] 本发明的技术解决方案的具体步骤结合【附图说明】如下： 参阅图1、图2。根据本发明，采用基于FPGA架构或DSP架构或FPGA+DSP架构的实时处理 平台，首先对灰度图像进行分