A High Dynamic Range Compression Technique Based on Infrared Contrast Enhancement
-
摘要: 如何将红外探测器采集的高动态范围的数据压缩为低动态范围图像数据的同时,能尽可能地保留图像的信息,提高图像的对比度一直是一个技术难点。针对这一问题,本文提出了一种新的红外图像压缩方法。该方法引入了直方图信息,通过对直方图进行分割,区分背景区域像素和目标区域像素; 然后计算压缩映射模型; 最后结合分割后的直方图对图像的像素采用不同强度的对比度增强。本文算法利用直方图区分背景区域像素与目标区域像素,在增强图像对比度时,能有效抑制背景噪声。通过实验对比,结果表明,本文所提出的算法更能较好地突出图像的细节,增强图像对比度。Abstract: It has always been technically difficult to compress the high dynamic range data collected by an infrared detector to low dynamic range image data, while preserving the image information as much as possible and improving the contrast of the image. To solve this problem, a new infrared image compression method was proposed. In this method, histogram information is introduced, and the pixels of the background and target regions are distinguished by the segmentation of the histogram. Then, the compression model is established. Finally, enhancing the contrast of the image pixels using different coefficients combines the segmented histogram. The algorithm proposed in this paper uses histogram information to distinguish the pixels of the background region and the pixels of the target region and can effectively suppress background noise when enhancing the image contrast. The experimental results show that the proposed algorithm can better highlight details and improve the contrast.
-
0. 引言
大视场多目标监视是目前红外搜索跟踪领域研究的热点,该领域正朝着高分辨率、超大视场成像的方向发展[1-3]。对应大视场下多目标的监视与跟踪任务,视场的增大,红外图像的背景成分变得极其复杂,又加上来自地面建筑物、地杂波等的干扰,极其容易在图像区域产生伪目标,造成红外全景搜索跟踪系统的虚警[4-7]。在实施全景搜索跟踪的任务过程中,视场中难免出现虚假运动目标,其中既有静止物体,也可能包含运动方向与真实目标不一致的物体。如果不能对以上虚假目标及时进行识别并剔除,将会严重影响后续搜索跟踪任务的执行。
针对红外搜索跟踪容易遇到的上述问题,王卫华[8]等人提出了一种基于兴趣区(region of interest)提取的目标实时监测算法,根据目标运动特性与灰度特性快速提取目标可能存在的兴趣区,再针对兴趣区内的局部目标图像切片,进一步剔除虚假目标干扰;陈炳文[9]提出一种基于时空融合与粒子滤波的弱目标检测算法,首先利用像素点的帧间和邻域相关性,采用短程非参数核建模算法检测候选目标,接着使用粒子滤波跟踪候选目标,以剔除虚假目标得到真实航迹;王雪梅[10]等提出了基于背景检测的单帧检测弱小目标改进算法,结合多方向多极梯度算法,对背景预测初步检测出来的可能目标点进一步检测,从而排除虚假目标;罗举平[11]等针对红外目标检测任务,提出了变门限检测与跟踪的方法,减少了虚假量测的数量。张艳艳[12]等引入Lucas-Kanade(LK)光流法,对光流结果进行最大类间方差分割,视光流不连续区域为真实运动目标;西北工业大学[13]利用Lucas-Kanade光流法对连续帧红外图像检测疑似运动区域,再通过SVD(singular value decomposition)分解与目标背景重构检测弱小目标;Lei Wang[14]等利用先验概率降低视频中目标检测的误报率;Recchia Giuseppe[15]等将Lucas-Kanade光流法应用于无人机光电视避系统,有效降低了误警率与虚警率。然而上述方法并不专门针对全景图像,因此对于全景图像中同一运动目标帧间距离过大造成的影响,这些方法并不能有效解决。类似方法甚至会错把真实运动目标当作虚假目标加以剔除,因而并不适用于降低红外搜跟系统的虚警率。本文针对以上问题提出了一种基于Lucas-Kanade光流特征的红外全景系统虚假目标甄别与剔除的方法。
1. 建立图像金字塔
选取前后相邻两帧红外全景图像I和J,为图像I和J分别建立金字塔表达${\{ {I^L}\} _{L = 0,...{L_m}}}$和${\{ {J^L}\} _{L = 0,...{L_m}}}$,这里Lm通常取为3,也就是金字塔层数一般为4。把原始图像作为第0层,通过高斯滤波并令宽高缩小2L倍的图像作为第L层,多层图像组合类似金字塔,如图 1所示。
2. 计算前后帧图像光流
由于红外全景搜索跟踪系统的任务在于利用360°全景影像实现对周围环境态势感知的能力,因此往往将目标定为包括人、车辆、舰船以及飞机等相比环境温度较高的物体。在实际图像中该类目标一般呈现较高的灰度值,因此本文利用二值化将图像中的高亮目标都事先提取出来,在这些高亮目标区域中随机选择若干像素记为该目标对应的像素点u,接下来利用LK光流法思想依次从最高层开始计算目标像素点u对应的LK光流dL,并取LK光流dL对应层光流初始值gL的和为下一层光流的初始值,直至迭代到原始图像层(第0层),最终获取目标像素点u在原始图像层的光流值按照如下步骤计算图像光流d,接着利用光流d的模值与角度为主要特征,使用聚类算法将真实运动目标和虚假运动目标区分开来。首先计算前后帧图像间的光流:
初始化最高层光流的初始值${g^{{L_m}}} = {\left[ {0\;\;0} \right]^{\text{T}}}$,令L=Lm,每迭代一次,L减1,直到L=0:
由于图像金字塔建立过程中,图像的长宽依次缩小一倍,因此针对第L层图像,计算图像IL中像素点u对应的位置:$u{\text{ = [}}u_x^L\;\;u_y^L{\text{] = }}u{\text{/}}{{\text{2}}^L}$;
计算图像IL在x方向的梯度:
$$ I_x^L(x,y) = \frac{{{I^L}(x + 1,y) - {I^L}(x - 1,y)}}{2} $$ (1) 计算图像IL在y方向的梯度:
$$ I_y^L(x,y) = \frac{{{I^L}(x,y + 1) - {I^L}(x,y - 1)}}{2} $$ (2) 计算空间矩阵:
$$ \boldsymbol{G} = \sum\limits_{x = u_x^L - {w_x}}^{x = u_x^L + {w_x}} {\sum\limits_{y = u_y^L - {w_y}}^{y = u_y^L + {w_y}} {\left[ {\begin{array}{*{20}{c}} {I_x^{\text{L}}I_x^{\text{L}}}&{I_x^LI_y^L} \\ {I_x^LI_y^L}&{I_y^LI_y^L} \end{array}} \right]} } $$ (3) 式中:wx、wy表示u处具有相同光流邻域的大小。
初始化迭代LK光流初始值d0=[0 0]T;
令k=1,每迭代一次k加1,直到匹配差异向量$ \left| {{\eta ^k}} \right| $小于1:
定义此时(x, y)处前后两帧图像的差异δI(x, y),则有:
$$ \delta {I^k}(x,y) = {I^L}(x,y) - {J^L}(x + g_x^L + d_x^{k - 1},y + g_y^L + d_y^{k - 1}) $$ (4) 计算图像差异:
$$ {b_k} = \sum\limits_{x = u_x^L - {w_x}}^{x = u_x^L + {w_x}} {\sum\limits_{y = u_y^L - {w_y}}^{y = u_y^L + {w_y}} {\left[ {\begin{array}{*{20}{c}} {\delta {I^k}(x,y){I_x}(x,y)} \\ {\delta {I^k}(x,y){I_y}(x,y)} \end{array}} \right]} } $$ (5) 计算匹配差异向量:
$$ η^{k}=\boldsymbol{ G}^{−1}b_{k} $$ (6) 计算光流:
为下一次迭代提供初始值dk=dk−1+ηk;
结束k上的迭代,得到第L层图像上的光流优化值dL=dk;
为第L-1层图像提供光流初始值gL−1=2(gL+dL);
结束L上的迭代,最终计算得到光流d=g0+d0;
得到图像J中对应点坐标v=u+d。
3. 重新表达目标光流
上述步骤求得前帧图像I中目标区域像素点u在前后两帧图像对应的光流d,以及该像素点在后帧图像J中的位置v。其中光流d为一个二维向量,该向量标识图像平面内由位置u指向位置v的运动位移。而同一个运动目标其运动位移应该是一致的,这就给区分真实运动目标和虚假运动目标带来了先决条件。针对图像I中经过二值化后获得的目标区域,每个区域选择m个像素点(u1, u2, …, um),本方案根据图像视场范围以及处理器运算能力综合考量选择m=5。对于每一个像素点ui,按照第2章所述方法得到对应的光流d(ui),然后利用算术平均值来求取对应目标区域Φ的光流dΦ:
$$ {d_\varPhi } = \frac{{\sum\limits_{i = 1}^m {d({u_i}} )}}{m} $$ (7) 4. 目标分类
对于地面建筑或者地杂波等产生的虚假目标,由于全景成像系统镜头本身的运动,造成其通过上述计算也会产生光流值,但该光流值并不能反映像素点的实际运动。虚假目标光流值与真实目标光流值会有显著不同,这主要表现在光流向量的方向和大小两个方面。
1)虚假目标由于本身静止或运动范围很小,在叠加了红外摄像机自身旋转运动后,其光流方向会有高度趋同性,而真实运动目标的光流方向应与其实际位移方向一致,从而造成虚假目标与运动目标在运动方向上的显著区别;
2)虚假目标因为本身基本处于静止状态,表征运动位移的光流仅反映红外摄像机自身的运动,实际观察光流值大小应不大于10个像素,这也是区别于实际运动物体的显著特征。
设$\overrightarrow d = (x,y)$,令$\left| d \right| = \sqrt {{x^2} + {y^2}} $表示光流向量的大小,令θ表示光流向量的方向,其中:
$$ \theta = \left\{ \begin{array}{l} \arccos \frac{x}{{\sqrt {{x^2} + {y^2}} }},y \geqslant 0 \hfill \\ 2{\text{π }} - \arccos \frac{x}{{\sqrt {{x^2} + {y^2}} }},y < 0 \hfill \\ \end{array} \right. $$ (8) 创建向量$\vec L = (\left| d \right|,\theta )$用于重新表达光流,目的在于后续根据光流区分真实运动目标与静止虚假目标。根据前述第3章中的分析,虚假静止目标的表征向量$\vec L$会具有很小的$\left| d \right|$值,且θ值会比较接近,将其标注到二维直角坐标系中虚假静止目标的向量L坐标会高度集中,而真实运动目标则会大范围随机分布,如图 2所示。
最后利用聚类算法将得到的目标光流进行分类,由图 2可以发现真实运动目标与虚假运动目标的光流分布区域有较大区别,为了把二者快速区分开来,首先在离θ轴较近的光流表征向量$\vec L$中随机选择一个作为第一个初始聚类中心,然后选择距离该中心最远的一个光流表征向量$\vec L$作为第二个初始聚类中心,通过计算各光流表征向量$\vec L$与初始化聚类中心的距离,将这些光流表征向量$\vec L$重新划归至距离最近的类别,并以新类别中各光流表征向量$\vec L$的均值作为新的聚类中心,持续迭代直至光流表征向量的分配聚类情况没有更新,以实现最终把真实运动目标(类别1)与虚假静止目标(类别2)分开的目的。
5. 实验结果
为验证本文的虚警率降低方法是否具有普遍性,项目组利用研发成功的红外全景成像设备采集不同环境下(室内、室外)的中波红外全景图像,在Windows操作系统下使用Visual C++开发环境将本文算法实现,将不同场景下前后两帧全景图作为输入,以验证算法有效性。
观察图 3和图 4可以发现,无论室外还是室内环境,复杂背景下总会引入干扰,同时搜索跟踪系统自身旋转导致其中的高亮物体都会得到表征运动的光值,因此很多由于干扰产生的伪目标极易被搜索跟踪系统误判为真实运动目标。引入本文算法后,系统对输入的前后两帧全景图像计算LK光流,随后将重点区域光流进行包括大小与方向的比对,按照聚类思想将其实施二分类,重复迭代直至真实运动物体被区分出来。根据图 3与图 4中处理后图像可以发现,红外全景图像中的真实运动物体均以实线框标记出来,而由于干扰造成的虚假运动目标被虚线框标记出来,由此可见本文算法较好地将红外搜索跟踪系统所获全景场景中的虚假运动目标加以抑制,从而有效降低了该类系统的虚警率。
6. 结论
红外全景监控设备是未来红外热像系统化发展的主要趋势,而大视场的引入必然加大后续图像处理设备的工作量,对于嵌入式图像处理系统由于本身运算能力有限,因此在预处理阶段自动排除虚假目标以减轻后续处理运算的压力就显得尤为必要。本文所述方法通过对实际采集的红外全景序列图像加以分析,很好地剔除了造成干扰的虚假运动目标,实验结果也验证了该方法的有效性。未来将把该方法向以GPU/FPGA+DSP为主处理器的嵌入式平台移植,同时针对特征提取阶段做更为深入的研究融入更多特征以提升虚假目标检测效果,为红外搜索跟踪系统等产品提供有效的技术支持。
-
表 1 不同算法在相同场景图像中的信息熵的对比结果
Table 1 Information entropy comparison results that different algorithms in the same scene image
Algorithm Scene 1 Scene 2 Scene 3 Scene 4 Scene 5 Linear mapping 7.1260 6.2631 7.1262 7.4640 7.2761 Literature[6]'s algorithm 7.2788 6.7646 7.0064 7.4837 7.4046 Literature[7]'s algorithm 7.6616 6.1768 7.3264 7.6890 7.4696 Our algorithm 7.8303 6.7306 7.6266 7.7293 7.7107 表 2 不同算法在相同场景图像中的峰值信噪比的对比结果
Table 2 PSNR comparison results that different algorithms in the same scene image
Algorithm Scene 1 Scene 2 Scene 3 Scene 4 Scene 5 Literature[6]'s algorithm 19.7842 24.5278 23.7493 19.2759 21.3895 Literature[7]'s algorithm 18.2785 18.7230 22.8438 20.1421 23.9762 Our algorithm 22.2963 23.3482 23.7893 21.5421 24.4664 -
[1] Silverman J. Display and enhancement of infrared images[C]//Image Processing and its Applications, 1992, International Conference on. IET, 1992.
[2] Rafael C Gonzalez, Richard E Woods. Digital image processing[J]. Prentice Hall International, 2008, 28(4): 484 - 486.
[3] 王炳健, 刘上乾, 周慧鑫, 等. 基于平台直方图的红外图像自适应增强算法[J]. 光子学报, 2005, 34(2): 299-301. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB20050200Z.htm WANG Bingjian, LIU Shangqian, ZHOU Huixin, et al. Self-adaptive contrast enhancement algorithm for infrared images based on plateau histgrom[J]. Acta Photonica Sinica, 2005, 34(2): 484 - 486. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB20050200Z.htm
[4] 宋岩峰, 邵晓鹏, 徐军. 基于双平台直方图的红外图像增强算法[J]. 红外与激光工程, 2008(2): 125-128. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ200802029.htm SONG Yanfeng. SHAO Xiaopeng, XU Jun. Infrared image enhancement algorithm based on dual platform histogram[J]. Infrared And Laser Engineering, 2008(2): 125-128. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ200802029.htm
[5] ZUO C, CHEN Q, LIU N. Display and detail enhancement for the visualization of high dynamic range infrared images[J]. Opt. Eng., 2011, 50(12): 127401. DOI: 10.1117/1.3659698
[6] HUANG J, YONG M, YING Z, et al. Infrared image enhancement algorithm based on adaptive histogram segmentation[J]. Applied Optics, 2017, 56(35): 9686. DOI: 10.1364/AO.56.009686
[7] Branchitta F, Diani M, Corsini G, et al. Dynamic-range compression and contrast enhancement in infrared imaging systems[J]. Optical Engineering, 2008, 47(7): 076401.1-076401.14. DOI: 10.1117/1.2956655
[8] Monobe Y, Yamashita H, Kurosawa T, et al. Dynamic range compression preserving local image contrast for digital video camera[J]. IEEE Transactions on Consumer Electronics, 2005, 51(1): 1-10. http://ieeexplore.ieee.org/document/1405691
[9] 王园园, 赵耀宏, 罗海波, 等. 海面红外图像的动态范围压缩及细节增强[J]. 红外与激光工程, 2019, 48(1): 307-315. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201901045.htm WANG Yuanyuan, ZHAO Yaohong, LUO Haibo, et al. Dynamic range compression and detail enhancement of sea-surface infrared image[J]. Infrared and Laser Engineering, 2019, 48(1): 307-315. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201901045.htm
[10] 张菲菲. 梯度域处理框架下的图像视见度增强技术研究[D]. 武汉: 武汉大学, 2015. WANG Feifei. Study on image visibility enhancement in the framework of gradient domain processing[D]. Wuhan: Wuhan University, 2015.
[11] 张临临. 基于图像分层和动态压缩的图像细节增强算法研究[D]. 西安: 西安电子科技大学, 2012. ZHANG Linlin. Study on image detail enhancement algorithm based on image stratification and dynamic compression[D]. Xi'an: XIDIAN University, 2012.
[12] 单瑞卿, 李斌, 韩伟, 等. 高动态范围红外图像的显示与细节增强[J]. 光学技术, 2019, 45(4): 475-481. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJS201904016.htm SHAN Ruiqing, LI Bin, HAN Wei, et al. Display and detail enhancement for high-dynamic-range infrared images[J]. Optical Technique, 2019, 45(4): 475-481. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJS201904016.htm
-
期刊类型引用(1)
1. 花文波,赵城慷,高大友,花文涛. 一种温度补偿调焦红外系统丢包检测方法. 红外技术. 2024(09): 1081-1086 . 本站查看
其他类型引用(0)