基于曲率滤波和视觉显著性的红外小目标检测

第50卷第4期
2023年
北京化工大学学报(自然科学版)
Journal of Beijing University of Chemical Technology (Natural Science)
Vol.50,No.4
2023
引用格式:王钰,赵明晶,谢晓明,等.基于曲率滤波和视觉显著性的红外小目标检测[J].北京化工大学学报(自然科学版),
2023,50(4):75-86.
WANG Yu,ZHAO MingJing,XIE XiaoMing,et al.Infrared small target detection based on curvature filtering and visual significance[J].Journal of Beijing University of Chemical Technology (Natural Science),2023,50(4):75-86.
基于曲率滤波和视觉显著性的红外小目标检测
王　钰1　赵明晶2*　谢晓明1*　李　伟2
(1.北京化工大学信息科学与技术学院,北京　100029;2.北京理工大学信息与电子学院,北京　100081)摘　要:红外成像因具有隐蔽性强㊁环境适应能力强和抗干扰能力强等优点,被广泛用于军事和民用领域㊂为了实现对红外小目标的高精度检测,提出了一种基于曲率滤波和视觉显著性相结合的红外小目标检测算法㊂首次将曲率滤波引入红外小目标检测中,对图像中的背景进行估计,然后将背景估计结果与原图像进行差分,使得图像中的大部分背景被显著抑制㊂为了减小残余的部分高强度杂波对目标的正确检测产生影响,利用形态学方法将高强度杂波移除㊂为了进一步提高目标的检测精度,提出了一种局部对比度算法进行杂波抑制及目标增强㊂最后,采用自适应阈值分割方法得到显著的小目标㊂对本文所提算法与其他算法在5个数据集上进行了比较分析,结果表明,本文算法的信杂比(SCR)与背景抑制因子(BSF)远高于其他算法,在检测率和误报率方面也明显优于其他算法㊂
关键词:红外小目标;曲率滤波;局部对比度;单帧检测中图分类号:TP751 DOI :10.13543/j.bhxbzr.2023.04.010
收稿日期:2022-08-18
基金项目:航空科学基金项目(ASFC -20200051072001)第一作者:女,1996年生,硕士生*通信联系人
赵明晶,E⁃mail:631240891@ 谢晓明,E⁃mail:xmxie@
引　言
红外成像因具有隐蔽性强㊁可全天时工作㊁可穿透云雾以及抗电磁干扰能力强等优点,在军事和民用领域得到了广泛应用[1]㊂同时,目标种类的多样性㊁环境的复杂性等对红外目标检测技术提出了更高的要求[2]㊂在实际获取图像的过程中,目标和背景通常会以非常快的速度变化,导致获得的红外图像整体比较模糊㊁缺乏纹理细节信息㊁目标与背景的对比度降低等问题
[3]
,这些问题都给小目标的精确
检测增加了难度㊂现有的红外小目标检测算法从单帧和多帧两个方面来解决以上难点㊂单帧检测算法主要利用帧内的信息,对于复杂场景图像中的小目标较难检测,但其检测速度较快,易于通过硬件实现;多帧检测算法需要依赖帧间的信息,而帧间信息
的连续性通常受到红外设备和目标快速变化的影响[4],其检测性能会下降㊂
近年来,人们对传统的单帧红外小目标检测算
法进行了很多相关研究㊂基于空域滤波的检测算法通过对背景估计来检测小目标,例如最大中值和最大均值滤波[5]㊁双边滤波[6]㊁形态学滤波Top -hat 变换[7]等算法㊂经典的Top -hat 变换广泛用于红外目标检测,但是它对噪声很敏感,在处理背景杂乱的图像时会出现大量的虚警,且依赖结构元素的选取㊂为了解决这些问题,Zhao 等[8]提出了一种新的形态特征提取算法(MMP)用于红外小目标检测,该算法能够更好地利用红外图像中的空间信息㊂此外,基于视觉显著性的算法也获得了很好的检测效果㊂Chen 等[9]在人类视觉对比度机制的启发下,提出了一种局部对比度算法(LCM),该算法是基于目标与周围背景的对比度差异来实现目标增强的,但是在实际应用中会出现很高的虚警率㊂为了克服上述缺点,人们提出了很多改进的算法㊂Shi 等[10]提出了高速的多尺度局部对比度算法(HB -MLCM),该算法在检测速度和检测能力方面表现很好;Han 等[11]根据红外图像中噪声的种类,提出了多尺度局部对比度算法(RLCM),该算法对复杂背景下小目标的
处理具有较好的鲁棒性;为了避免区域交叠造成的目标漏检,穆靖等[12]提出了三层模板局部差异度量算法(TTLDM),该算法具有很好的实时性并且避免了多尺度运算导致的算法复杂度提高,但此类算法中大多数不能兼顾实时性和检测性能,对背景杂波的抑制也不充分,导致虚警率较高㊂还有一类算法利用红外图像背景的非局部自相关特性和目标的稀疏性,将图像的背景和目标分别视为稀疏矩阵和低秩矩阵进行处理,取得较好的效果[13-14]㊂如Gao 等[13]提出的红外图像块模型(IPI)就是经典的非局部先验算法,该模型对于高度变化的复杂场景具有很好的适应性,但此类算法的复杂度高,实时效果差㊂此外,人们还将一类新的算法应用于红外小目标检测,这种算法通过将红外图像的三维灰度图看成空间曲面,利用目标与背景的空间曲率差异来最小化图像曲率,进而实现对目标的检测㊂例如Zhao 等[15]提出了主曲率函数滤波的检测算法,该算法对于复杂背景下的目标检测效果较好,但此类算法需要计算每个像素的主曲率,其计算复杂度高,且要求图像二次可导,这一假设较难满足㊂为了解决这一问题,Gong 等[16]提出了曲率滤波理论,通过将图像视为局部可展曲面来最小化曲率,目前该算法在图像去噪方面已得到广泛使用,但是还未应用于红外小目标检测㊂
针对以上算法中存在的问题,本文提出了一种基于曲率滤波和视觉显著性相结合的红外小目标检测算法㊂首先采用改进的曲率滤波算法和背景差分操作来抑制图像中的大部分背景,然后提出了一种局部对比度算法,以进一步提高目标和背景之间的
对比度㊂与其他算法相比,本文提出的算法能够更有效地检测小目标,并且具有较低的误报率㊂
1　基于曲率滤波和视觉显著性相结合
的红外小目标检测算法
1.1　红外图像特性分析
在远距离情况下,红外图像容易受到红外设备和复杂天气的影响,形成的图像具有如下特点㊂
(1)背景通常分布较为均匀且占据图片中大量的像素,具有一定的连续性,其内部对比度不突出;
(2)一般来说,小目标只占据图像的几个或者
十几个像素,其大小不确定,需要自适应地检测目标大小㊂此外,由于目标的高速运动和背景的快速变化,形成的小目标可能会模糊,缺乏结构纹理信息,帧间信息也可能不连续㊂受到距离和设备等因素的影响,小目标的灰度通常也不明显,容易被背景杂波所淹没,导致其内部对比度不突出㊂
基于以上分析,本文提出了一种将曲率滤波和视觉显著性相结合的红外小目标检测算法㊂该算法的检测流程如图1所示㊂针对红外图像中目标占据的像素很小,大量的背景及变化容易对目标检测结果造成干扰的问题,提出了改进曲率滤波的背景抑制方法㊂采用全变分曲率滤波算法估计图像的整体背景,将得到的背景估计结果与原图进行差分,并用形态学方法进一步处理差分后的残余杂波,以减弱背景对目标的干扰㊂然后,依据人类视觉显著注意力机制的特点,提出了一种局部对比度算法,以进一步提高目标和背景之间的对比度,从而提高目标的检测率㊂最后,采用自适应阈值分割方法进行目标检测
㊂
图1　本文所提算法的总流程图
Fig.1　General flowchart of the proposed algorithm
㊃67㊃北京化工大学学报(自然科学版) 2023年
1.2　基于改进曲率滤波的背景抑制
在红外小目标检测中,由于拍摄距离通常较远,得到的图像中目标呈现出 弱”和 小”的特性㊂背景变化通常也会对目标产生严重干扰,使得目标不易被正确检测㊂因此,首先对背景进行抑制是提高目标检测率的一种重要手段㊂全变分曲率滤波算法[16]能够较好地去除图像中的噪声,因此本文采用全变分曲率滤波算法对红外图像的背景部分进行估计,然后将背景估计结果与原图像进行差分,得到对红外图像中的背景进行初步抑制的结果,这样能够极大地减少背景对目标正确检测的干扰㊂
全变分曲率滤波算法以变分模型为理论基础,通过构造的局部滤波器快速求解变分模型㊂与传统的变分模型求解相比,它不需要计算曲率,因此不需要假设图像是可导的㊂该算法将图像视为分段常值可展曲面,采用局部近似图像进行滤波,通过一定次数的迭代使得图像尽可能可展,从而逼近原始图像㊂
全变分曲率滤波算法通过构造一个3×3的窗口,计算中心像素R i,j到其邻域内所有切平面的投影距离d i(i=1,2, ,8)㊂
d1=R i-1,j-1+R i-1,j+R i,j-1+R i+1,j-1+R i+1,j
5-R i,j d2=R i-1,j+R i-1,j+1+R i,j+1+R i+1,j+R i+1,j+1
5-R i,j d3=R i-1,j-1+R i-1,j+R i-1,j+1+R i,j-1+R i,j+1
5-R i,j d4=R i+1,j-1+R i+1,j+R i+1,j+1+R i,j-1+R i,j+1
5-R i,j d5=R i-1,j-1+R i-1,j+R i-1,j+1+R i,j-1+R i+1,j-1
5-R i,j d6=R i-1,j-1+R i-1,j+R i-1,j+1+R i,j+1+R i+1,j+1
5-R i,j d7=R i+1,j-1+R i+1,j+R i+1,j+1+R i-1,j-1+R i,j-1
5-R i,j d8=R i+1,j-1+R i+1,j+R i+1,j+1+R i-1,j+1+R i,j+1
5-R i,j
(1)
然后通过选择最小切平面距离来校正中心像素值,其计算式如下㊂
Q i,j=R i,j+d m
s.t.|d m|=min{|d i|},i=1,2, ,8(2)式中:Q i,j为校正后的中心像素值,d m为找到的最小切平面距离㊂通过式(1)和(2)可以有效地判断中心像素是噪声点还是边缘点:当中心像素是一个噪声点时,它将会被平滑;当中心像素是一个边缘点时,它将不会被平滑㊂当经过一定次数的迭代后,就可以得到背景图像Q㊂因此,可以通过从原始的红外图像中减去背景图像Q来获得背景抑制的结果,将得到的结果记为A,公式如下㊂
A i,j=max{I i,j-Q i,j,0}(3)式中:A i,j为背景差分的图像,I i,j为原始红外图像㊂由式(3)易知平滑背景的输出为零,目标区域会被保留,大部分背景区域被显著抑制㊂但是全变分曲率滤波算法在对背景进行估计时会使得图像中的纹理细节信息过于平滑,因此在进行背景差分之后,会留下一些高强度的残余杂波,对目标的正确检测产生影响㊂为此,本文采用5×5的圆盘算子进行先腐蚀后膨胀的形态学操作来进一步减弱这些高强度的残余杂波㊂
基于改进曲率滤波的背景抑制过程如图2所示,其中红色方框部分为结果放大的区域,绿色圆框部分为目标区域㊂由图2(b)可以明显看出,经过曲率滤波操作之后,图像中的背景被显著抑制,目标得到明显突出㊂由图2(c)可以看出,经过形态学操作后,高强度的残余杂波得到有效减弱㊂最终得到的图像更适合使用局部对比度算法来增强小目标并抑制背景,也就是说在真正的小目标区域内局部对比度会更突出,而其他区域的局部对比度经过上述操作后会变弱㊂
1.3　基于局部对比度的目标增强
经典的局部对比度算法依据目标与邻域背景的差异性,将每个目标区域的最大像素点与其周围区域像素点的均值比作为目标区域的增强因子,这样极易增强图像中原本存在的噪声,造成大量虚警,并且由于需要逐像素运算,计算成本较高,实时性也会受到影响㊂因此,本文提出了一种改进的局部对比度算法,通过差值局部对比描述和多尺度运算结果融合来增强目标及扩大局部差异㊂
从图2(c)中可以看出,在经过改进曲率滤波操作后,大部分背景和高强度残余杂波得到抑制㊂为了提高目标检测率㊁降低虚警率,进一步提升目标和背景之间的对比度是至关重要的㊂因此,本文提出了一种增强局部差的方法,所使用的窗口模型如图3所示,模型分为内外两个区域,内部为目标区域T,外部为其邻域背景区域B㊂
采用小目标与其邻域背景的灰度均值之差d来
㊃77㊃
第4期 王　钰等:基于曲率滤波和视觉显著性的红外小目标检测
图2　基于改进曲率滤波的背景抑制
Fig.2　Background suppression based on improved
curvature filtering
图3　窗口模型
Fig.3　Window model 描述局部对比度(式(4))㊂
d =m T -m B
(4)m T =
1N T
∑N T
i =1
v i (5)
m B =
1
N B
∑N B
j =1
v j (6)
式中:m T 为目标区域的灰度均值,m B 为背景区域的灰度均值,N T 为目标区域的像素数,N B 为背景区域的像素数,v i 为目标区域第i 个像素的灰度值,v j 为背景区域第j 个像素的灰度值㊂
在实际应用中,红外小目标的大小通常是不断变化的,当小目标进入红外搜索和跟踪系统时,随着距离的改变,其大小变化范围一般在2×2和12×
12像素之间[17],因此需要选取不同窗口大小的目标区域来估计真实的红外小目标㊂
L k =m k T -m k
B ,k =1,2,3,4
(7)
式中:k 为当前选取的第k 个目标尺度,m k T 和m k
B 分
别为当前尺度下目标区域和背景区域的像素均值,L k 为当前选取的第k 个目标尺度的局部对比度㊂将得到的不同尺度的局部对比图进行Had⁃
amard 乘积后,可以获得显著的小目标㊂为了进一步增强红外小目标与背景的差异,对获得的显著目标结果进行平方,计算如下㊂S =(L 1 L 2 L 3 L 4)2
(8)
式中:S 为最终获得的局部对比度图㊂
从以上定义中可以看出,当目前像素属于背景
像素时,有L k ≈0和S ≈0,当目前像素属于目标像素时,有L k >0和S ≫0㊂因此,使用Hadamard 乘积能够进一步扩大目标和背景之间的对比度,从而达到使亮度较大的目标更亮,亮度较小的目标更暗的
目的㊂本文将小目标的区域大小分别设置为3×3㊁
5×5㊁7×7和9×9,背景区域大小设为固定尺寸15×15,这样有利于提高检测速度㊂1.4　自适应阈值分割
经过以上处理之后得到最终的小目标显著图,为了使得到的小目标更加直观,采用一种阈值分割方法[18]分离目标㊂
Thr =μ+λ×σ
(9)
式中:Thr 为阈值分割后的结果图;μ为背景均值;σ为背景标准差;λ的取值根据最佳的图像分割效果来确定,对于不同的算法和数据集而言,λ的取值差
异较大,其取值范围一般在10~90之间㊂
2　算法验证
为了验证本文所提算法的有效性,将本文算法与
IPI 算法[13]㊁LCM 算法[9]㊁RLCM 算法[11]㊁Top -hat 算法[7]㊁HB -MLCM 算法[10]和TTLDM 算法[12]进行对
㊃87㊃北京化工大学学报(自然科学版) 2023年
比分析㊂选取5组不同场景下的红外图像(见表1)进行测试,测试环境在MATLAB R2021a 中编译㊂
表1　选取的红外图像数据集
Table 1　Selected infrared image datasets
数据集帧数像素背景描述
Dataset130256×200较强云杂波空天场景
Dataset280
450×400目标淹没在较强云杂波空天场景Dataset3100300×256目标淹没在较弱云杂波空天场景
Dataset440256×200较弱云杂波空天场景Dataset5
40300×300
低空多建筑场景2.1　评价指标
为了评价不同红外小目标检测算法的背景抑制和目标增强效果,采用信杂比(SCR)㊁信杂比增益(SCRG)㊁背景抑制因子(BSF)作为评价指标,其计算式如下㊂
R SCR
=
|μt -μb |σb (10)G SCRG
=
SCR out SCR in (11)F BSF =
C in C out
(12)
式中:R SCR 为信杂比,G SCRG 为信杂比增益,F BSF 为背景抑制因子,μt 为目标的平均像素大小,μb 为目标周围区域的像素值大小,σb 为目标周围的像素值标准差,SCR in 为输入图像的信杂比,SCR out 为输出图像的信杂比,C in 和C out 分别为原图像和经过处理后的输出图像的标准差㊂SCR 值越大,小目标越容易被检测到;SCRG 反映了目标的输入输出相对于背景的增强程度,也可以用来描述小目标检测的难度,SCRG 值越大,目标的增强程度越大;BSF 反映了背景的抑制效果,BSF 值越大,抑制效果越好㊂对于检测结果,通常利用检测率(P d )和虚警率
(F a )来评价目标的检测精度,其计算式如下㊂
P d =N true
N actual
×100%(13)F a =
N false
N total
×100%(14)
式中:N true 为检测到的真实目标数,N actual 为总的真实目标数,N false 为检测到的虚假目标数,N total 为检测到的所有目标数㊂当以下条件都符合时,认为检测到的小目标是正确的[19]:(1)检测到的目标和真实目标的中心像素之差小于5个像素;(2)真实目标和检测到的目标的像素有重叠㊂
2.2　参数选择
2.2.1　迭代次数
在实验过程中发现,全变分曲率滤波的迭代次数会对算法的检测性能造成一定影响㊂本文比较了不同迭代次数下所提算法的性能,结果见表2㊂可以看出,当迭代次数为10次时在Dataset5中出现了虚警,当迭代次数为20次时算法性能达到了最佳,当迭代次数继续增加时算法性能并未出现变化,但检测效率下降㊂因此,最终选择20次为全变分曲率滤波的迭代次数㊂
表2　迭代次数对所提算法性能的影响
Table 2　Effect of the number of iterations on the
performance of the proposed algorithm
迭代次数Dataset1
Dataset2
Dataset3
Dataset4
Dataset5
P d /%
F a /%P d /%F a /
%P d /%F a /
%P d /%F a /
%P d /%
F a /
%1010001000100010001004.762010001000100010001000301000100010001000100050
10001000100010001000100100
100
0100
0100
0100
2.2.2　目标尺度
表3为不同目标尺度对本文算法的SCR㊁SCRG 和BSF 的影响㊂可以看出,当选用的目标尺度大小
表3　本文算法在不同目标尺度下的SCR㊁SCRG 和BSF 平均值
Table 3　Average values of SCR,SCRG and BSF of the proposed algorithm for different target scales
目标尺度Dataset1Dataset2Dataset3Dataset4Dataset5SCR
SCRG
BSF
SCR
SCRG
BSF
SCR
SCRG
BSF
SCR
SCRG BSF
SCR
SCRG
BSF
3×3㊁5×5
154.2866.3969.29269.23179.33458.93287.63394.00123.2845.427.7862.57347.3046.23
96.90
3×3㊁5×5㊁7×7
163.4667.5772.49262.02174.38462.02287.17401.06125.6344.217.5671.16316.5442.60134.923×3㊁5×5㊁7×7㊁9×9
175.5771.75
75.10260.77173.71463.07291.05406.34126.0450.75
8.73
71.19348.4746.38365.97
粗体数字代表最优结果㊂
㊃
97㊃第4期 王　钰等:基于曲率滤波和视觉显著性的红外小目标检测
为3×3和5×5时,数据集2的SCR 和SCRG 值大于其他目标尺度的融合结果;当目标尺度大小为
3×3㊁5×5㊁7×7和9×9时,数据集1㊁3㊁4㊁5的SCR㊁SCRG 和BSF 值均大于其他目标尺度的融合结果㊂因此,选用3×3㊁5×5㊁7×7和9×9的目标尺度大小进行目标结果融合会对本文算法产生较好的效果㊂2.3　结果分析2.3.1　定性分析
为了直观地比较不同算法的检测效果,从5个数据集中分别选取1张示例图片来展示检测结果及相应的三维灰度图,结果如图4~8所示㊂其中,三维灰度图的x 轴㊁y 轴分别表示图像像素的横坐标和纵坐标,z 轴表示图像的像素值,图像中的小目标由红框框出㊂可以看出,对于示例图片1~4,LCM㊁RLCM 和Top -hat 算法的处理结果中均含有大量的背景杂波和噪声,而IPI 和HB -MLCM 算法的检测结果较好,仅存在少量的背景杂波和噪声㊂对于示例图片5,以上5种算法表现均不佳,均含有大量的背景杂波和噪声㊂TTLDM 算法对示例图片1㊁2和4的检测效果不佳,其结果中含有大量的背景杂波和噪声㊂而本文所提算法的检测结果在各个示例图片中均没有背景杂波,目标的亮度能够被很好地提高,表现出较好的检测性能
㊂
图4　不同算法对示例图片1的检测结果
Fig.4　Detection results of different algorithms for example picture 1
2.3.2　定量分析
表4为不同算法的SCR㊁SCRG 和BSF 值㊂可以看出,LCM㊁RLCM㊁Top -hat 算法的SCR㊁SCRG㊁BSF 值均较低,IPI㊁HB -MLCM 和TTLDM 算法在不
同数据集中的SCR㊁SCRG 和BSF 值仅次于本文的算法,但与本文算法的差距很大,本文所提算法的结果最优㊂
表5为不同算法的检测率和虚警率㊂可以看
㊃08㊃北京化工大学学报(自然科学版) 2023年
图5　不同算法对示例图片2的检测结果
Fig.5　Detection results of different algorithms for example picture 2
表4　不同算法的SCR㊁SCRG 和BSF 平均值
Table 4　Average values of SCR,SCRG and BSF for different algorithms
算法Dataset1Dataset2Dataset3Dataset4Dataset5SCR
SCRG BSF
SCR
SCRG
BSF
SCR
SCRG BSF
SCR
SCRG BSF
SCR
SCRG BSF
IPI
141.2651.2822.3526.9016.41216.04136.37207.2342.0122.613.9524.69111.5115.3811.83LCM
3.511.290.70
4.062.671.03
5.528.720.894.020.751.13
6.18
0.820.98
RLCM
6.48
2.482.107.59
4.8818.499.4112.793.374.32
0.816.6211.841.5816.54Top -hat
16.015.63
1.87
11.257.1020.807.34
8.37
2.11
16.822.996.16
17.752.41
4.53HB -MLCM 59.8021.5620.0061.8940.38290.3768.5187.35
49.7830.575.3747.1449.556.637.23TTLDM
42.83
13.2616.8348.81
31.1962.05219.09254.7861.8121.813.8422.4697.9413.5761.64本文算法175.5771.7575.10260.77173.71463.07291.05406.34126.0450.75
8.73
71.19348.4746.38
365.97
粗体数字代表最优结果,下划线数字代表次优结果㊂
出:在检测率为100%的情况下LCM 算法的表现最差,在不同的示例图片上均有不同程度的虚警率;
IPI㊁RLCM㊁Top -hat 和TTLDM 算法在处理数据集时均有误检的情况;HB -MLCM 算法在处理数据集1~
㊃
18㊃第4期 王　钰等:基于曲率滤波和视觉显著性的红外小目标检测
图6　不同算法对示例图片3的检测结果
Fig.6　Detection results of different algorithms for example picture 3
表5　不同算法的检测率和虚警率
Table 5　Detection rates and false alarm rates of different algorithms
算法Dataset1
Dataset2
Dataset3
Dataset4
Dataset5
P d /%F a /%P d /%F a /%P d /%F a /%P d /%F a /%P d /%F a /%
IPI
1000
701.751000
97.50
10088.06LCM
10028.5710052.11000.991004.7610088.44RLCM
10056.5210011.11
1006.981006.981000
Top -hat
10030.2310032.21000.991000100
80.6HB -MLCM 100010001000100010097.88
TTLDM 1006.2591.2558.87
100010001000本文算法
100
01000
100
0100
010004时检测率很高,表现出很好的检测性能,但在处理数据集5时出现了大量虚警;TTLDM 算法可以较好地处理示例图片3~5,但在处理示例图片1和2时
均有误检的情况;本文所提出的算法可以在保证检测率为100%的情况下同时保持较低的虚警率,其结果均优于其他算法㊂
㊃28㊃北京化工大学学报(自然科学版) 2023年
图7　不同算法对示例图片4的检测结果Fig.7　Detection results of different algorithms for example picture4
表6为采用不同算法处理数据集的单帧平均运行时间㊂对每个数据集分别进行实验,取数据集总时间的平均值作为单帧的运行时间㊂可以看出,IPI 和RLCM的运行速度较慢,LCM㊁HB-MLCM㊁Top-
hat㊁TTLDM和本文算法的运行速度较快㊂
以上结果表明,LCM㊁RLCM和Top-hat算法对背景的抑制能力较差,目标增强的效果也不如其他算法;IPI算法容易造成目标漏检,使得检测率下降,且该算法的实时性较差;HB-MLCM算法在处理数据集5时BSF值很小,对具有强边缘背景的图像(如示例图片5)的抑制能力较差,并且SCR㊁SCRG 和BSF值远小于本文提出的算法;从检测率和虚警率来看,TTLDM算法对于目标被背景淹没的图像(如示例图片1和2)的检测能力较差,容易造成大量虚警,该算法适合检测目标与背景有明显差异的图像,此外,TTLDM算法的SCR㊁SCRG和BSF
表6　不同算法的单帧平均运行时间
Table6　Average running time of single frame for
different algorithms
算法
单帧平均运行时间/s
Dataset1Dataset2Dataset3Dataset4Dataset5 IPI5.43153.9511.204.6722.96 LCM0.460.550.520.460.49 RLCM2.766.253.932.774.20 Top-hat0.470.470.470.420.44 HB-MLCM0.480.500.490.460.47 TTLDM0.430.500.500.430.40本文算法0.471.740.520.490.53
值与本文提出的算法存在较大差异,虽然TTLDM算法具有较高的实时性,但其检测结果不太理想㊂因此,与其他算法相比,本文算法的综合检测性能最优㊂
㊃38㊃
第4期 王　钰等:基于曲率滤波和视觉显著性的红外小目标检测
图8　不同算法对示例图片5的检测结果Fig.8　Detection results of different algorithms for example picture5
3　结论
本文提出了一种将曲率滤波和视觉显著性相结合的红外小目标检测算法㊂首次将曲率滤波引入红外小目标检测中,通过背景差分操作来实现对背景的显著抑制,并采用形态学方法将残留的部分高强度噪声移除㊂为了加强目标与背景间的差异,采用目标区域与邻域背景间的灰度均值差来描述这种差异,然后将多尺度运算结果进行Hadamard乘积,使目标得到显著增强,背景得到显著抑制㊂最后,采用自适应阈值分割方法提取到显著的小目标㊂相较于其他算法,本文所提的算法不仅具有较快的检测速率,而且在检测精度上也有明显的提升㊂参考文献:
[1]　刘征,杨德振,李江勇,等.红外单帧弱小目标检测
算法研究综述[J].激光与红外,2022,52(2):154-
162.
LIU Z,YANG D Z,LI J Y,et al.A review of infrared
single⁃frame dim small target detection algorithms[J].
Laser&Infrared,2022,52(2):154-162.(in Chi⁃
nese)
[2]　范晋祥,杨建宇.红外成像探测技术发展趋势分析
[J].红外与激光工程,2012,41(12):3145-3153.
FAN J X,YANG J Y.Development trends of infrared im⁃
aging detecting technology[J].Infrared and Laser Engi⁃
neering,2012,41(12):3145-3153.(in Chinese) [3]　ZHANG H,ZHANG L,YUAN D,et al.Infrared small
target detection based on local intensity and gradient
㊃48㊃北京化工大学学报(自然科学版) 2023年
properties[J].Infrared Physics&Technology,2018,
89:88-96.
[4]　ZHAO M J,LI W,LI L,et al.Single⁃frame infrared
small⁃target detection:a survey[J].IEEE Geoscience
and Remote Sensing Magazine,2022,10(2):87-119.
[5]　DESHPANDE S D,ER M H,RONDA V,et al.Max⁃
mean and max⁃median filters for detection of small targets
[C]//Proceedings of SPIE,Signal and Data Processing
of Small Targets.Denver,1999:74-83. [6]　TOMASI C,MANDUCHI R.Bilateral filtering for gray
and color images[C]//Proceedings of Sixth IEEE Inter⁃
national Conference on Computer Vision.Bombay,1998:
839-846.
[7]　BAI X Z,ZHOU F G.Analysis of new top⁃hat transforma⁃
tion and the application for infrared dim small target detec⁃
tion[J].Pattern Recognition,2010,43(6):2145-
2156.
[8]　ZHAO M J,LI L,LI W,et al.Infrared small⁃target de⁃
tection based on multiple morphological profiles[J].
IEEE Transactions on Geoscience and Remote Sensing,
2021,59(7):6077-6091.
[9]　CHEN C L P,LI H,WEI Y T,et al.A local contrast
method for small infrared target detection[J].IEEE
Transactions on Geoscience and Remote Sensing,2014,
52(1):574-581.
[10]SHI Y F,WEI Y T,YAO H,et al.High⁃boost⁃based
multiscale local contrast measure for infrared small target
detection[J].IEEE Geoscience and Remote Sensing Let⁃
ters,2018,15(1):33-37.
[11]HAN J H,LIANG K,ZHOU B,et al.Infrared small tar⁃
get detection utilizing the multiscale relative local contrast
measure[J].IEEE Geoscience and Remote Sensing Let⁃
ters,2018,15(4):612-616.
[12]穆靖,李伟华,饶俊民,等.采用三层模板局部差异
度量的红外弱小目标检测[J].光学精密工程,2022,
30(7):869-882.
MU J,LI W H,RAO J M,et al.Infrared small target
detection using tri⁃layer template local difference measure
[J].Optics and Precision Engineering,2022,30(7):
869-882.(in Chinese)
[13]GAO C Q,MENG D Y,YANG Y,et al.Infrared patch⁃
image model for small target detection in a single image
[J].IEEE Transactions on Image Processing,2013,
22(12):4996-5009.
[14]ZHAO M J,LI W,LI L,et al.Infrared small⁃target de⁃
tection based on three⁃order tensor creation and tucker de⁃
composition[C]//Proceedings of the2021IEEE Interna⁃
tional Geoscience and Remote Sensing Symposium.Brus⁃
sels,2021:3129-3132.
[15]ZHAO Y,PAN H B,DU C P,et al.Principal curvature
for infrared small target detection[J].Infrared Physics&
Technology,2015,69:36-43.
[16]GONG Y H,SBALZARINI I F.Curvature filters effi⁃
ciently reduce certain variational energies[J].IEEE
Transactions on Image Processing,2017,26(4):1786-
1798.
[17]KIM S,LEE J.Scale invariant small target detection by
optimizing signal⁃to⁃clutter ratio in heterogeneous back⁃
ground for infrared search and track[J].Pattern Recog⁃
nition,2012,45(1):393-406.
[18]汪奎伟.红外小目标的检测与跟踪[D].大连:大连
理工大学,2013.
WANG K W.Infrared small target detection and tracking
[D].Dalian:Dalian University of Technology,2013.
(in Chinese)
[19]LI W,ZHAO M J,DENG X Y,et al.Infrared small tar⁃
get detection using local and nonlocal spatial information
[J].IEEE Journal of Selected Topics in Applied Earth
Observations and Remote Sensing,2019,12(9):3677-
3689.
㊃58㊃
第4期 王　钰等:基于曲率滤波和视觉显著性的红外小目标检测。