APP下载

基于自适应分数阶微分的红外目标增强算法

2020-05-06代少升李东阳聂合文

红外技术 2020年3期
关键词:微分红外局部

代少升,李东阳,聂合文,姚 俐

基于自适应分数阶微分的红外目标增强算法

代少升,李东阳,聂合文,姚 俐

(重庆邮电大学 通信与信息工程学院,重庆 400065)

针对红外图像存在灰度范围窄、图像细节不清晰、目标边缘模糊的问题,提出了一种基于自适应分数阶微分的红外目标增强方法。该方法首先利用图像的梯度、信息熵进行有效融合,并且自适应调整分数阶微分以增强图像中的目标边缘;然后采用图像像素灰度的标准差和均值进行融合去确定目标的分割阈值,以区分出图像中的背景和目标部分;通过对图像中的目标区域进行线性增强,以进一步突显目标。经过实验验证:本文提出的方法能够有效地区分红外图像中的目标和背景,局部目标背景比(Target-to-Background Ratio,TBR)平均提高了0.5,视觉效果比较理想。

红外图像;目标增强;自适应分数阶微分;线性变换;局部目标背景比

0 引言

近些年来,随着红外技术不断发展,红外热成像在军事和民用上得到广泛应用[1]。由于红外图像多呈现昏暗、模糊,信噪比一般较低[2],并且红外目标经常被淹没在背景之中,所以红外目标的增强作为红外目标检测的前期工作尤为重要,影响后续的红外目标的检测以及跟踪[3]。红外目标增强主要有点增强、空域增强、频域增强3种类型。最早采用的点增强法主要利用了图像的灰度变换和几何变换,通过对比度增强和对比度拉伸等方法达到增强图像目标的目的。比如Zhang等人[4]提出了一种对比度受限自适应直方图均衡(Contrast Limited Adaptive Histogram Equalization,CLAHE)算法,对局部直方图的高度进行限制,以限制局部对比度的增强幅度;虽然增强了图像目标的对比度和限制了背景噪声的增加;但是也造成了目标增强程度的下降。另一方面,红外目标空域增强利用了图像的空间信息;通过目标的位置、形状、大小等特征区分目标和背景,从而有效地增强图像目标。在这方面,Luan等人[5]提出了一种多分辨多尺度的红外图像增强算法,利用目标的边缘轮廓特征进行模糊融合聚类和多分辨小波分解;虽然提高了图像信噪比并使得目标更为突出,但是目标细节信息的增强效果不好。除了点增强和空域增强,人们也从频域的角度研究了图像目标的增强。Qi等人[6]提出一种基于四元傅里叶相位谱(Phase Spectrum ofQuaternion Fourier Transform,PQFT)的红外目标增强方法,利用四元傅里叶变换的相位谱和导数滤波器区分目标和背景;虽然增强了红外目标,但是只适应于小目标。

本文提出了一种基于自适应分数阶微分的红外目标增强方法。为了解决红外目标整体对比度不高和边缘不清晰的问题,首先对红外图像进行自适应分数阶微分滤波处理,然后对红外图像进行局部目标线性增强,最后将两种目标增强效果融合。该算法不仅增强了目标的边缘细节,而且增加了目标的对比度和避免了背景噪声的增加,具有很好的红外目标增强视觉效果。

1 自适应分数阶微分

1.1 自适应参数融合

为了更好地获得红外目标局部特征信息,利用图像的梯度和局部信息熵联合得到目标特征,从而获得整幅图像的动态参数。如图1,图像统一采用3×3的模板进行矩形局部区域特征计算。目标区域梯度采用图像的水平梯度、垂直梯度、主对角线梯度和副对角线梯度联合获得,所以目标包含8个方向的梯度特征。图像的梯度()计算如下:

式中:h()是图像第点像素的水平梯度;v()是图像第点像素的垂直梯度;p()是图像第点像素的主对角线梯度;c()是图像第点像素的副对角线梯度。图像的局部熵()计算如下:

式中:1≤≤256,是像素值;是模板大小;Hist是以第点像素为中心的模板图像经过局部直方图统计所得的大小像素值的个数。图像的自适应特征参数值大小对应于图像每个像素点的增强程度,主要用0~1范围的数值表示,值越大,增强程度越强。为了使特征参数值映射到0~1之内,对以上两个特征参数分别进行归一化处理:

式中:()是图像第个像素点的一个特征参数;Max(())是特征参数()的最大值;Min(())是特征参数()的最小值。最后进行两个特征参数的融合。由于两个特征所占权重大小差别不大,所以采用求和平均法融合,计算公式如下:

1.2 分数阶微分增强

分数阶微分是数学一个重要分支,是由整数阶微积分演化而来。为了准确地使用分数阶微分对图像处理的计算,采用Grumwald-Letnikov[7]定义的分数阶微分以増强图像细节。分数阶微分优于整数阶微分,更加有利于突出目标边界细节。分数阶微分使用由G-L推广得到的阶分数阶微分一元函数表达式:

进行分数阶微分主要目的是增强图像高频部分,保留图像的中低频部分。图像的高频部分具体指图像的纹理和边缘部分,其中也包括目标的边缘部分;中低频主要对应于图像背景平滑区域。当阶数增大,图像的高频部分被增强;然而当过大,则过量増强目标部分和过少保留中低频部分,所以利用图像特征参数自适应调节微分阶数的大小,并且的取值范围在0~1范围之间。图像每个像素的分数阶微分増强阶数()赋值如下:

()=SIG() (7)

式中:SIG()是图像每个像素点的局部融合特征值。最后对整个图像所有像素点进行自适应地模版卷积滤波处理,图像目标区域对应的像素点的SIG()较大,对应的()也就较大,增强效果较明显;并且图像背景区域对应的像素点的SIG()较小,()也就较小,较好地保留背景信息。红外图像经过自适应分数阶微分增强后的效果如图4。在图4(b)中,目标边缘清晰和很好地被增强。

2 局部目标线性增强

线性拉伸可以增强图像的对比度,但是如果对整幅图像拉伸会导致目标与背景同时增强,目标增强效果不好;因此采用局部目标线性拉伸处理红外图像。局部目标线性拉伸首先采用图1中的局部模板对红外图像进行目标提取,然后对所提取的目标区域再进行线性变换处理。本文采用阀值判断法区别目标和背景区域。根据图像的局部均值和局部标准差融合得到目标特征值CTR(),图像的局部均值()和局部标准差()计算如下:

CTR()=1()+2() (10)

图2 图像像素点的8个方向

图3 5×5分数阶微分模板

图4 自适应分数阶微分增强

式中:1和2是融合系数。由于目标的局部均值的权重大于局部标准差,所以目标大小分别选取为0.7和0.3。当CTR()高于一定阀值Th时,如果图像第点像素属于目标区域,进行线性变换处理,反之不进行处理。图像的分割阀值Th是由目标特征值CTR()采用迭代法求得最佳数值,适应于各种不同模糊程度的图像。阀值Th主要计算方法是: 首先取出红外图像目标特征值CTR()的最大值和最小值,进行求均值得到初步阀值Tn;然后将红外图像中目标特征值CTR()大于阀值Tn的区域分割为一部分,小于阀值Tn的区域分割为另一部分;分别求出两部分的图像区域像素的目标特征值CTR()的均值,再进行平均得到阀值Th;如果Th和Tn的差值大于0.005,将Th赋值给Tn继续重复以上步骤,迭代优化直到取得目标和背景的分割阀值Th的最佳值。最后对图像阀值分割后的目标区域进行局部线性变换处理,主要过程是首先对阀值分割后的图像目标区域计算局部模板最大值max()和局部模板最小值min(),计算如下:

max()=max(m()) (11)

min()=min(m()) (12)

式中:m()是以图像目标区域中第个像素点中心的模板区域。然后进行像素线性区间变换,目标区域的像素值区间拉伸到局部模板最小值到1之间,计算如下:

式中:S()是图像目标区域中第个像素点经过拉伸后的像素值;()是图像目标区域中第个像素点的像素值。目标区域得到合适增强,背景区域很好地被保留,所以只有目标区域被拉伸和呈现高亮。局部目标线性增强效果如图5。

3 目标增强算法融合

红外图像在经过分数阶微分增强目标后,虽然目标边缘清晰,但是目标内部对比度过低;然而经过局部目标线性增强后,虽然目标内部得到拉伸,但是目标边缘不清晰。所以将两种方法处理后的结果进行叠加融合,为了更好地突出目标内部像素,局部目标线性增强后的结果占较大比重。对红外图像进行算法处理后的对比图以及像素分布图如图6。在图6(b)中,目标物体被明显增强,亮度明显提升。分数阶微分增强使边缘明显更加清晰,局部线性目标拉伸使目标明亮较完整。通过对比图6(c)和(d),更能观察到目标整体的像素值得到提高,而目标周围背景几乎没有被增强,避免了大部分的背景干扰。

图5 局部目标线性增强

4 结果与分析

为了分析本文算法的性能且突出在红外目标增强的效果,首先将自适应直方图均衡化[8](Adaptive Histogram Equalization,AHE)、CLAHE、小波变换[9](Wavelet Transform,WT)和PQFT 4种不同的目标增强算法和本文算法处理相同的红外图像进行对比;然后通过局部TBR指标评价各种目标增强方法的优劣。本文采用以森林和大地为背景的红外图像数据集。图像的尺寸大小为240×320。图像包含两种目标,分别是4个人体目标和少量汽车目标,如图7(a)。本文所有算法仿真实验都使用同一台个人笔记本电脑,电脑硬件性能为Intel i3-3110M 2.3GHz CPU核心、12G运行内存,使用的系统为windows7 64位,仿真软件为Matlab。

在图7中,图7(a)是原始图像。图像昏暗,目标模糊且边界不清晰,特别右侧人体目标和背景相近。如图7(b),原图经过AHE算法处理后,人和车两种目标都得到增强,但是背景同样也得到增强;图像整体亮度提高,但是模糊不清晰。如图7(c),CLAHE算法相比于AHE算法,背景区域噪声没有被增强;虽然目标区域得到增强,但是目标增强程度过小。如图7(d),WT算法处理后的图像目标明显被增强,人和车的一部分明亮突出,但是图像背景过度增强导致亮度过高且明显失真。图7(e)是PQFT算法,虽然图像背景噪声得到很好地抑制并且背景保留完整和不失真;但是目标的边缘和整体增强效果不好。图7(f)是本文算法处理后的结果,相比于以上4种算法,图像目标不仅明亮,而且目标边缘清晰,人体轮廓和肢体部分明显区分,车的部分细节十分清晰,很好地保留目标的特征;图像背景也十分清晰,噪声抑制效果良好,并且得到很好地保留,没有出现失真。由以上分析得出:本文算法在一定程度上优于以上4种目标增强算法。

为了全面地从数值上分析本文算法性能的优劣性能,使用局部TBR对比分析本文算法和其它4种算法。局部TBR是目标区域和目标周围的局部背景区域的像素均值之比,从目标和目标周围背景的对比程度反映红外图像目标增强的程度,局部TBR计算公式如下:

式中:T是图像目标区域像素点灰度的均值;B是图像目标周围背景区域像素点灰度的均值。局部TBR的值越高,说明该算法对红外图像目标增强的程度越大。本文算法与4种目标增强算法分别在人和车的部分两种目标的局部TBR值进行对比,如表1。在人体目标的增强方面,可以明显地得出本文算法局部TBR的值最大,相比于其它4种算法,平均提高0.5,计算公式如下:

式中:TBRAHE、TBRCLAHE、TBRWT、TBRPQFT、TBROUR分别是AHE、CLAHE、WT、PQFT和本文算法的局部TBR值,avg是本文算法相比其它4种算法局部TBR提高的平均值。所以本文算法对图像中的人体目标增强程度最大。同样在车的部分增强方面,除了AHE算法TBR的值优于PQFT算法外,与图像人体目标增强一样,本文算法的TBR的值最大。最终说明本文算法在红外图像目标增强程度效果最佳。综合以上所有实验结果得出,本文算法有效地增强了红外图像目标,优于其它4种目标增强算法。

表1 不同算法对于不同目标的局部TBR

5 结论

本文提出了一种新的红外目标增强方法,主要解决红外图像的目标边缘细节不清晰和对比度低的问题。首先利用了图像的梯度和局部信息熵联合得到新的目标边缘特征,作为阶数调整分数阶微分,实现了分数阶微分自适应增强图像目标;目标的边缘细节得到很好地增强,并且目标本身和图像背景更加清晰。然后实现了一种局部线性目标增强,由图像的局部均值和局部标准差融合得到的目标整体特征再进行自适应迭代化处理;分割出图像中的目标区域后,采用线性变换增强目标内部,得到的图像目标整体增强效果较好。经过实验结果表明,本文算法能够有效地区分红外图像中的目标和背景,适用于红外目标的增强。

[1] 吴旭景, 杜斌. 红外热成像无损检测技术现状及发展的相关研究[J]. 化工管理, 2018(29): 183.

WU X J, DU B. Research on the status and development of Infrared Thermal Imaging Nondestructive Testing Technology[J]., 2018(29): 183.

[2] 王好贤, 董衡, 周志权. 红外单帧图像弱小目标检测技术综述[J]. 激光与光电子学进展, 2019, 56(8): 9-22.

WANG H X, DONG H, ZHOU Z Q. Review on dim small target detection technologies in infrared single frame image[J]., 2019, 56(8): 9-22.

[3] 易诗, 张洋溢, 聂焱, 等. 红外图像中快速运动目标的检测与跟踪方法[J].红外技术, 2019, 41(3): 268-272.

YI Shi, ZHANG Yangyi, NIE Yan, et al. Fast moving target detection and tracking method in infrared image[J]., 2019, 41(3): 268-272.

[4] 张丽. 对比度受限自适应直方图均衡方法[J]. 电脑知识与技术, 2010, 6(9): 2238-2241.

ZHANG L. Contrast limited adaptive histogram equalization[J]., 2010, 6(9): 2238-2241.

[5] 栾孟杰. 一种多分辨多尺度的红外图像增强算法[J]. 激光杂志, 2019, 40(8): 81-84.

LUAN M J. A multiresolution and multiscale infrared image enhancement algorithm[J]., 2019, 40(8): 81-84.

[6] QI S X, MA J, LI H, et al. Infrared small target enhancement via phase spectrum of Quaternion Fourier Transform[J]., 2014, 62: 50-58.

[7] 牛为华, 孟建良, 崔克彬, 等. 利用Grümwald-Letnikov分数阶方向导数的图像增强方法[J]. 计算机辅助设计与图形学学报, 2016, 28(1): 129-137.

NIU W H, MENG J L, CUI K B, et al. Image enhancement method using Grümwald-Letnikov fractional directional derivative[J]., 2016, 28(1): 129-137.

[8] 刘德全, 崔涛, 杨雅宁. 局部对比度自适应直方图均衡化图像增强的算法研究[J]. 信息与电脑(理论版), 2016(7): 79-80.

LIU D Q, CUI T, YANG Y N. Algorithm Research on local contrast adaptive histogram equalization image enhancement[J].(), 2016(7): 79-80.

[9] 王瑞. 小波变换在红外图像处理中的应用研究[D]. 淮南: 安徽理工大学, 2016.

WANG R. The application of wavelet transform in infrared image processing[D]. Huainan: Anhui University of technology, 2016.

Linear Enhancement Algorithm of Infrared Target Based on Adaptive Fractional Differentiation

DAI Shaosheng,LI Dongyang,NIE Hewen,YAO Li

(,Chongqing University of Posts and Telecommunications, Chongqing 400065, China)

To solve the problems associated with infrared images, such as narrow gray range, unclear image details and fuzzy target edge, an infrared target enhancement method based on adaptive fractional differentiation is proposed. In this method, first, the gradient and information entropy of image are used for effective fusion, and the fractional differentiation is adaptively adjusted to enhance the edge of the target in the image. Subsequently, the standard deviation and mean value of the image pixel gray are fused to determine the segmentation threshold of the target, to distinguish the background and target in the image. The target area of the image is linearly enhanced to better highlight the target. Experimental results show that the proposed method can effectively distinguish the target and background in the infrared image. The average local target-to-background ratio (TBR) increased by 0.5, and the visual effect was ideal.

infrared image, target enhancement, adaptive fractional differentiation, linear transformation, local target-to-background ratio

TP391

A

1001-8891(2020)05-0257-07

2019-10-31;

2020-03-06.

代少升(1974-),男,河南潢川人,教授,主要从事红外图像处理方向的研究。E-mail: daiss@cqupt.edu.com。

国家自然科学基金(61671094);重庆市科学技术委员会国家科学基金(CSTC2015JCYJA40032)。

猜你喜欢

微分红外局部
与由分数阶Laplace算子生成的热半群相关的微分变换算子的有界性
网红外卖
一类带有Slit-strips型积分边值条件的分数阶微分方程及微分包含解的存在性
爨体兰亭集序(局部)
闪亮的中国红外『芯』
凡·高《夜晚露天咖啡座》局部[荷兰]
8路红外遥控电路
TS系列红外传感器在嵌入式控制系统中的应用
基于跟踪微分器的高超声速飞行器减步控制
丁学军作品