第51卷 第2期 激光与红外Vol.51,No.2 2021年2月 LASER & INFRAREDFebruary,2021
文章编号:1001 5078(2021)02 0250 07·图像与信号处理·基于引导滤波的高动态红外图像压缩算法
汪 璇1,赵金博2,刘智嘉3,高旭辉2,夏寅辉2
(1 湖北大学知行学院计算机与信息工程学院,湖北武汉430011;2 北京波谱华光科技有限公司,北京100015;
3 华北光电技术研究所,北京100015)
摘 要:非制冷红外机芯的原始图像位宽为14Bit,虽然相对于8Bit红外图像来说,14Bit红
外原始图像具有灰度动态范围更广、灵敏度更高、包含细节信息更多的优点,但是由于普通显
示设备能够显示的最大灰度范围为8Bit,所以需要对14Bit原始图像进行压缩,以满足常规显
示设备显示及后端图像处理的需求。若压缩算法性能不佳,在图像压缩过程中可能会丢失大
量的细节信息,直接影响成像质量。本文提出了一种基于引导滤波的压缩和显示算法,该算法
首先利用引导滤波对图像进行分层,根据不同图层的特点和在合成图像中的作用,分别进行增
强和降噪处理,使合成图片具有良好的显示效果。通过与常用算法的仿真对比实验,以及对视
觉效果和定量评价参数两个方面的实验结果分析,本文所提出算法在压缩和图像处理的性能
表现上均获得一定程度的提升。
关键词:红外;高动态;引导滤波;分层
中图分类号:TP391 41 文献标识码:A DOI:10.3969/j.issn.1001 5078.2021.02.021Highdynamicinfraredimagecompressionalgorithmbased
onguidedfiltering
WANGXuan1,ZHAOJin bo2,LIUZhi jia2,GAOXu hui2,XIAYin hui2
(1.CollegeofComputerandInformationEngineeringofZhiXingCollege,HubeiUniversity,Wuhan430011,China;
2.BeijingBopOpto ElectronicsTechnologyCo.,Ltd.,Beijing100015,China;
3.NorthChinaResearchInstituteofElectro Optics,Beijing100015,China)
Abstract:TheoriginalimageofuncooledinfraredmovementhasaBitwidthof14Bit.Comparedwiththe8Bitinfra
redimage,the14Bitinfraredoriginalimagehastheadvantagesofwidergraydynamicrange,highersensitivityand
moredetailedinformation.However,becausethemax
imumgrayrangethatordinarydisplayequipmentcanmakeis
8Bit,itisnecessarytocompressthe14Bitoriginalimagetomeettherequirementsofconventionaldisplayequipment
andback endimageprocessingrequirements.Iftheperformanceofthecompressionalgorithmispoor,alotofdetails
maybelostintheprocessofimagecompression,whichdirectlyaffectstheimagequality.Inthispaper,acompression
anddisplayalgorithmbasedonguidedfilterisproposed.Firstly,guidedfilterisusedtolayertheimage.Accordingto
thecharacteristicsofdifferentlayersandtheroleinthesyntheticimage,theenhancementandnoisereductionarecar
riedoutrespectively,sothatthesyntheticimagehasagooddisplayeffect.Bycomparingthesimulationresultswith
commonalgorithms,andanalyzingtheexperimentalresultsofvisualeffectandquantitativeevaluationparameters,itis
provedthatthealgorithmproposedinthispaperhasacertainimprovementintheoverallperformance.
Keywords:infrared;highdynamic;guidedfiltering;delamination
作者简介:汪 璇(1984-),女,硕士,讲师,主要研究方向为信号处理。E mail:61341883@qq com
收稿日期:2020 04 11;修订日期:2020 05 14
1 引 言
红外成像技术[1]
在军事侦察、工业生产、民用
生活等领域具有重要地位。目前,作为应用最为广泛的高位宽红外探测器,其非制冷红外机芯可输出位宽达14Bit的原始红外图像。与普通8Bit图像相比,
14Bit红外图像具有灰度动态范围更广,灵敏度更高,图像细节信息更丰富的特点。灰度动态范
围[
2]
即红外图像中亮度最大区域与亮度最小区域的亮度之比,红外图像灰度动态范围越大,其亮度信息越丰富,细节信息越明显。由于普通显示设备仅能显示8B
it灰度等级图像,14Bit红外图像的灰度范围远超过普通显示设备的响应范围,所以这样的图像
被称为高动态范围图像。由于高动态范围图像比传统显示设备的动态范围更广,因此必须将14Bit
原始红外图像的动态范围进行压缩至8Bit范围。高动态范围图像处理的核心思想是对红外摄像设备采集到的原始高动态图像信号进行数模转换和量化处理。灰度变化平缓的景物,生成图像的灰度分布会集中在较狭窄的区间内,而灰度变化剧烈的景物所生成图像的灰度会散落在范围较宽的区间。对于高动态范围图像红外图像处理,关键问题在于将1
4Bit图像中的信息转化为人眼可观察到的8Bit图像中同时还要保持图像原有的细节信息,具备较好的对比度以供人眼观察。
在此过程中,使用错误的压缩算法可能导致原始图像中重要信息丢失。另外,由于红外成像技术和设备的特殊性,红外成像常常出现整体亮度较暗、对比度低、边缘模糊、噪声大、细节表现力较差等成像缺点。为了解决上述问题,本文以非制冷红外图像为背景,在现有算法基础上提出了一个改进算法,该方法可以对高动态范围图像进行压缩并可以较好的显示图像细节,并通过仿真实验结果证明算法的有效性。
2 常用算法及原理
在红外图像处理领域中,自适应增益控制(AGC)是调节图像动态范围和对比度的最常用方法之一,该算法利用梯度线性变换原理将14Bit数据线性映射为8B
it数据,随后自适应调节图像的平均亮度和对比度。为方便叙述进行如下定义:[rmin,rmax]为原始图像的灰度取值区间,[smin,smax]为处理后的红外图像的灰度取值区间,r为原始图像的灰
度值,s为处理之后红外图像的灰度值,则灰度转化公式为:
s=smin+smax-smin
rmax-rmin
×(r-rmin
)(1)
式中,smax-smin
rmax-rmin
为压缩映射函数中的对比度增益。
由于14Bit高动态范围红外图像的灰度范围为0至214,所以导致rmax-rmin的数值较
大的同时对比度增益的数值却非常小,使得图像压缩后的图像动态范围(
s)取值过小,造成图像像素灰度过于集中,细节信息丢失严重,所以必须适当调整对比度增益。本文提出的自适应增益控制方法的核心思想是通过去掉一些极端灰度并自适应调节增益值,使压缩后图像尽可能利用可显示灰度范围,令图像整体具有较好的对比度和细节显示效果。
该方法的优点是高动态红外图像直方图中灰度值较为极端且数量较少的像素被算法直接清除,降低了极端灰度值对对比度增益的影响。随着极端灰度值像素消失,对比度增益变大,图像显示效果得到提升。不过,由于自适应增益控制算法中自适应仅针对删除某灰度像素的选择而言,映射增益和亮度范围仍需要人工调整,且动态范围采用是的线性映射转换的方法,不能显示原始高动态范围的详细信息。所以,对于高动态范围红外图像该算法无法同时显示每个灰度间隔,导致原始图像中大量细节信息丢失。
图1展示了使用自动增益控制算法显示不同动态范围的图像。图像中有许多细节,故须人为调整到动态范围的低响应区域(a)和高响应区域(b)来
搜索不同的细节。
图1 AGC手动调节图像的动态范围
Fig 1AGCmanuallyadjuststhedynamicrangeoftheimage
从图1中可明显看出,虽然自适应增益控制在一定的动态范围内实现了画面对比度增强,然而该方法的弊端是会丢失某些灰度间隔中携带的细节信
1
52激光与红外 No.2 2021 汪 璇等 基于引导滤波的高动态红外图像压缩算法
息。而且,图像显示出明显的过亮或过暗区域,整体显示效果欠佳。
为了改变线性转换算法的不足,直方图均衡化首先利用非线性拉伸的思路处理原始高动态范围图像,然后重新分配红外图像的像素值,令某些灰度范围内的像素数量近似于均匀分布,从而增强原始图像直方图中间峰段的对比度并降低两侧谷段的对
比度。
图2 直方图均衡化示意图Fig 2Histogramequalizationdiagram
图2所展示为灰度图像经过直方图均衡化前后直方图的变化情况,左侧为均衡化处理之前,右侧为均衡化处理之后。可以看出图像的灰度直方图形状由较为集中转变为接近平坦。其基本原理是通过对原始图像的像素灰度进行映射转化,使转化后的图像灰度的概率密度近趋于平坦,使不同灰度值的像素数量近似一致来达到增加图像的对比度的目的,便于视觉观察。图3为使用直方图均衡化算法的效
果图。
图3 直方图均衡化后的图像Fig 3Imageafterhistogramequalization
由图3可以看出,直方图均衡化法同样解决了高动态范围的原始红外图像在普通显示设备上的不可视性。但是图像背景灰度之间的间隔和占比增大,分布更加密集,背景区域对比度增大、噪点增多,图像中观测目标模糊且出现刷白现象。此外,由于直方图中间灰度间隔的过度拉伸使两端灰度的像素堆积,图像中出现了过暗和过亮的问题。
3 基于引导滤波图像分层的红外图像细节增强算法
为了改善前述算法中出现的问题,本文提出了一种基于引导滤波的高动态范围图像处理方法,通过将原始图像分为基本层和细节层,然后分别对相应部分进行灰度调整和增强并将处理后结果进行叠加,以实现对原始景物中的大动态信息和8Bit图像中的小动态细节信息的突出和增强。整体流程如图4
所示。
图4 引导滤波算法整体流程图
Fig 4Overallflowchartofguidedfilteringalgorithm
3 1 引导波图像分层
引导滤波的基本原理是通过定义图像中任意像素与其相邻像素之间的随机关系,使图像中全部像素具有不相同线性关系,根据每个像素的线性关系
到与原始图像梯度一致的图像基层[
3]
,用引导滤波的方法把原始图像分离成表征景物轮廓的基本层与表征景物纹理细节、噪声信息的细节层,即:
P=q+e
(2)
式中,基本层q与引导图像I之间为线性关系,表达式为:
qi=akIi+bk, i∈ωk
(3)
式中,k为引导图像中某个像素;ωk表示以k为圆心半径为r的局域窗口;ak和bk为引导图像在窗口内的比例系数。对上式两边求导得到:
dq=adI
(4)
由式(
4)可得知基本层q与引导图像I的梯度线性相关,因此,导引滤波能够较好的保持图像的边缘特征。
2
52激光与红外 第51卷
引导滤波方法就是寻使原始图像p与处理后图像q差异最小ak和bk的最优解,通常选择线性回归算法使拟合函数的输出值与原始图像最为接近,在窗口ωk内,
其代价函数为:Eak,b()k
=∑((akIi+bk-pi)2+εa2
k)(
5)式中,ε的作用是过滤删除过大的ak,
计算可得:ak=
1w∑i∈wk
(akIi+bk-p()i
2+εa2
k)δ2
k+ε
(6)bk=pk-akμk
(7)
式(6)和(7)中窗口内总像素数量为w,引导图像I在窗口ωk内的均值和方差分别为μk和σ2
k
,pk为窗口输入原始图像p的均值。3 2 直方图均衡的基础层压缩
由于基本层图像涵盖了全局背景信息,动态范围大且对比度低,经引导波分层处理后,其基本层依然具有很大的动态范围,需进一步压缩基本层图像,控制其灰度范围和调节对比度,本文使用一种基于直方图均衡化的压缩方法:首先设置图像统计直方图T为上限阈值,当灰度频数P大于T,则将T值赋于P,否则保持P值不变,直方图累积直方图函数如下:
S=T()r=∑r
k=
0PT(ek
)
(8)
通过设定合适的阈值,使用该函数能够删除一部分像素,而且由于在利用累积直方图进行灰度重新分配的过程中对灰度重新归一化,所以该方法能够在一定的范围内重新调整灰度间距。之后进行压缩处理,设置基本层压缩后的动态范围为Pb且令Pb=255,便于与细节层的合成操作。压缩后的基本层为:
sb=Pb×∑
r
k=0
PTr
()k
(9)
3 3 细节层的增强与噪声抑制
细节层是通过从图像基层中减去原始图像而获得,其灰度范围较广且携带细节信息的同时也包含大量的噪声信息,降噪处理必不可少。通过对数据信息进行分析,发现信号和噪声具有不同的特征:通常情况下,数据中信号幅值明显大于噪声幅值,且信号间往往具有较高的关联性而噪声
间为随机分布。可根据上述特征对信号与噪声进行区别。
根据人类视觉机理的相关研究发现,人眼视觉对图像平坦区域呈现的细节较为敏感,且容易忽略灰度变化较强区域的细节信息。基于这种掩蔽效
应,
Anderson等人[4]
提出了噪声可见性函数,来表征噪声的可视性特点。该函数为:
fni
,()j=1
Mi,()j
×θ+1(10)
式中,Mi,()j
为表征图像噪声可视化的函数;θ是一个取值范围为[0,1]的调节函数。图像区域越平坦,该区域噪声越明显,且可视化函数中的
Mi,()j越小、fi,()j越接近1;灰度变化越剧烈的区域噪声越不明显,其可视化函数中Mi,()j越大,fi,()j
越接近于0。根据上述特点可以从两个角度限制噪声放大:第一,限制噪声所在范围数据的增益系数;第二,根据局域特征调节增益系数,即设置平坦区域增强系数较小,剧烈区域增益系数相对较大。
假设细节放大增益Gi,()j的最大值为Gmax,最小值为Gmin,由噪声可视化函数可知,当fni,()j趋近0时,Gi,()j趋近于Gmax;fni,()j趋近于1时,Gi,()j趋近于Gmax。为简化运算,
细节图的增强函数构造为线性函数:
Gi,()j=Gmin+[1-f(i,j)](Gmax-Gmin
)(11)
通过仿真实验发现当增益函数中最大和最小增益分别为Gmax=1 5、Gmin=1时,图像可满足大多数应用需要,增强效果较好。4 实验与评价
本次仿真实验所使用的计算机主要硬件参数为IntelI78850H、RAM16G;使用操作系统为Win dows10及软件MatlabR2016a;实验数据为非制冷型中波红外焦平面探测器输出640×512原始图像,ADC分辨率为14Bit。
由图5可看出,
经过引导滤波后的基本层图像中目标边缘比原图更加锐化。由于细节层为原始图像与基本层相减得到,所以图像所呈现的是图像的高频特征,其中包含了大部分强边缘及较多噪声,图像细节较为细腻,动态范围也较小,远小于高斯滤波后细节层的范围。
3
52激光与红外 No.2 2021 汪 璇等 基于引导滤波的高动态红外图像压缩算法
图5 原始图像与引导滤波分层后的基本层与细节层图像
Fig 5Thebasiclayerandthedetaillayerimageafter
theoriginalimageandtheguidedfilterlayering
利用直方图均衡化法将基本层压缩,处理前后
的基本层如图6
所示。
图6 直方图均衡化处理前后的基本层
Fig 6Basiclayerbeforeandafterhistogramequaliza
tion
从图6图中可以对比看出经过直方图均衡化处
理后的图像相比于处理前图像的对比度明显增加,
且没有出现大面积的过度增强现象。细节层利用自
适应增强的方法的处理效果如图7
所示。
图7 自适应增强前后的细节层
Fig 7Detaillayerbeforeandafteradaptiveenhancement
由图7可以看出经过权重模板函数的自适应增
强后,细节层在噪声被大量清除的同时尽可能的保
留了细节信息,并且增加了图像的对比度。证明当
G
max
=1 5、G
min
=1时,增强效果较好,该数值也可
以根据具体探测背景通过实验自定确定,以达到不
同背景的最佳效果。
图8为采用不同融合系数(α)融合后8Bit输出
图像效果图,图(b)的纹理细节较图(a)更为明显,
但是图(a)的对比度相对均衡,更有利于视觉观察。
可根据不同场景和需求选择合适的融合系数,以达
到使用者的相应目的。
图8 不同比例融合的细节增强图
Fig 8Detailenhancementwithdifferentscalefusion
为了更准确验证算法性能,对多个红外图像压
直方图均衡化的基本原理缩与显示算法进行对比试验。图9图中(a)、(b)、
(c)、(d)分别是原始高动态范围图像、传统高频增
强算法处理图像、非锐化掩模增强算法处理图像和
引导滤波增强处理图像,原始图像为非制冷型640
×512红外探测器所成图像,模数转换分辨率为
14Bit
。
图9 四种图像效果展示
Fig 9Fourimageeffectsdisplay4
5
2激光与红外 第51卷
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论