基于GLMB滤波的复杂场景下红外弱小目标自适应跟踪算法

蔡如华, 周健斌, 吴孙勇, 郑翔飞

蔡如华, 周健斌, 吴孙勇, 郑翔飞. 基于GLMB滤波的复杂场景下红外弱小目标自适应跟踪算法[J]. 红外技术, 2024, 46(7): 743-753.
引用本文: 蔡如华, 周健斌, 吴孙勇, 郑翔飞. 基于GLMB滤波的复杂场景下红外弱小目标自适应跟踪算法[J]. 红外技术, 2024, 46(7): 743-753.
CAI Ruhua, ZHOU Jianbin, WU Sunyong, ZHENG Xiangfei. An Adaptive Tracking Algorithm for Infrared Dim Small Targets in Complex Scenes Based on GLMB Filter[J]. Infrared Technology , 2024, 46(7): 743-753.
Citation: CAI Ruhua, ZHOU Jianbin, WU Sunyong, ZHENG Xiangfei. An Adaptive Tracking Algorithm for Infrared Dim Small Targets in Complex Scenes Based on GLMB Filter[J]. Infrared Technology , 2024, 46(7): 743-753.

基于GLMB滤波的复杂场景下红外弱小目标自适应跟踪算法

基金项目: 

桂林电子科技大学研究生教育创新计划资助项目 2022YCXS145

国家自然科学基金(地区科学基金项目) 62263007

详细信息
    作者简介:

    蔡如华(1971-),男,副教授,硕士,主要从事目标跟踪方面的研究。E-mail: 934019492@qq.com

    通讯作者:

    吴孙勇(1981-),男,教授,博士,主要从事目标跟踪方面的研究。E-mail: 53505713@qq.com

  • 中图分类号: TN219

An Adaptive Tracking Algorithm for Infrared Dim Small Targets in Complex Scenes Based on GLMB Filter

  • 摘要:

    针对红外弱小目标在复杂场景下受到漏检和杂波影响,导致跟踪不连续甚至失效的问题,本文提出一种红外弱小目标自适应跟踪算法。在预处理阶段,为了减少不必要的计算,首先定义一种衡量图像复杂度的算法。然后该算法通过计算红外图像多个特征得到场景复杂度来确认场景类型,再根据场景类型选取对应的检测算法提取目标候选位置、灰度以及局部直方图等特征建立对应的量测模型与似然函数。在目标跟踪阶段,为了自适应地匹配广义标签多伯努利(Generalized Labeled Multi-Bernoulli, GLMB)滤波器的滤波参数,在GLMB的基础上提出一种适应视频图像的新生算法进行航迹起始;针对红外图像序列目标检测概率未知的情况,将未知检测概率的基数化概率假设密度(Cardinality Probability Hypothesis Density, CPHD)滤波器集成到GLMB中实时估计目标检测概率以提升跟踪精度。仿真结果表明,所提出算法能有效地排除量测漏检和虚警的干扰,跟踪不同红外复杂场景下的弱小目标。

    Abstract:

    In this study, we propose a robust adaptive tracking algorithm for infrared dim objects that addresses the problem of tracking discontinuities and failures caused by missed detections and clutter in complex scenes. In the pre-processing stage, an algorithm that measures the image complexity eliminates unnecessary calculations. This algorithm determines the scene type by calculating multiple features of the infrared image to obtain the scene complexity, and then selects the corresponding detection algorithm to extract the target candidate location, grayscale and local histogram features. Subsequently, a measurement model and likelihood function are established based on the scene type. In the tracking stage, to flexibly match the filtering parameters of the generalized labeled multi-Bernoulli (GLMB) filter, an adaptive algorithm suitable for video image distribution is proposed for track initiation. Aiming at the unknown detection probability of an infrared image sequence, a cardinality probability hypothesis density (CPHD) filter was integrated into the GLMB to estimate the detection probability of the target in real time, thereby improving the accuracy of the tracker. The simulation results show that the proposed algorithm can effectively track small infrared objects in different complex scenarios.

  • 对于红外焦平面探测器,噪声等效温差是一个很重要的指标参数,与信噪声比相关,直观地反映了器件的灵敏度,也是热像仪设计考虑的关键参数,影响系统的作用距离与目标识别[1-6]

    在较早期文献资料中,单元探测器与信噪比相关的主要指标参数包括探测率、噪声等效功率等,而NETD只是一个表征热像仪系统温度灵敏度的重要指标,直到20世纪八九十年代开始红外焦平面探测器快速发展后,NETD才逐渐出现在探测器的主要性能指标中[7-13]

    探测器的NETD与热像仪系统的NETD基本定义是相同的[14]:噪声电压Vn与目标温差ΔT产生的信号电压ΔVs相等时的温差称为噪声等效温差,即信噪比为1时,目标与背景的温差。定义公式如下:

    $$ NETD = \frac{{\Delta T}}{{\Delta {V_{\rm s}}/{V_n}}} $$ (1)

    实际应用中,红外焦平面探测器的NETD与多种因素有关,如像元面积、光谱响应范围、量子效率、冷屏F数、背景温度等。除此之外,从定义公式可以看到,NETD并未像探测率一样对积分时间进行归一化处理,实际中器件的信噪比并不随着积分时间的加长而按比例增加,因此同一器件在不同的积分时间条件下,NETD在理论上也不同。另一方面,器件的响应与黑体温度也并不是线性关系,选择的背景温度和ΔT不同,NETD在理论上也会不同,同一背景温度条件下ΔT越小越接近实际应用情况,但由于受黑体精度的影响,ΔT很难控制得很小,采集到的响应信号也很容易受到其他因素的干扰,因此一般取较大的温差进行测试,如中长波探测器测试时温差ΔT通常取10℃~15℃[15-17],这样带来的问题是,与实际情况可能存在的差异。理论上,NETD与电荷存储量无关,但在某些特定条件下,如应用的背景温度或固定的势阱填充量,需要采取一些措施来满足使用或测试要求,此时NETD会受到一定的限制。

    在热像仪系统应用中,有的文献还提出了由于焦平面探测器像元之间响应不均匀性引起的空间NETD[18-21],为了简化计算,本文不考虑空间NETD,仅从NETD的测试方法和流程出发,通过黑体辐射及探测器的理论公式,对HgCdTe、InSb等常用光伏型红外焦平面探测器NETD的理论计算公式进行推导,分析影响器件NETD的主要因素,通过典型640×512 HgCdTe中波红外焦平面探测器对不同条件下的背景限NETD进行理论计算及分析,并对固定背景温度但不同测试温差下的NETD进行理论推算及对比测试。

    与探测率等性能指标测试方法一样,红外焦平面探测器NETD在测试时也是测出一定条件下噪声和温差响应信号,再根据式(1)进行计算,其测试系统示意图如图 1,包括辐射源、驱动电路、数据采集模块、计算机及控制软件等。一般情况下,采用面源黑体作为红外辐射源,设置温度T0模拟背景温度,T1模拟目标温度;驱动电路提供器件正常工作所必须的电源、偏置和脉冲时序等,并对器件的输出信号进行调理;数据采集模块对调理后的器件输出信号进行自动采集,并送至计算机中;计算机及控制软件对驱动电路、数据采集模块等进行控制,并对采集数据进行处理计算。

    图  1  红外焦平面探测器NETD测试系统
    Figure  1.  NETD test system for infrared FPA detectors

    具体测试流程如下:

    ① 连接测试系统与被测器件,使面源黑体面充满器件视场,器件达到工作温度后,通过驱动电路给被测器件加电,在规定的测试条件下使器件正常工作;

    ② 分别在黑体温度为T0T1条件下采集器件中每个像元的输出信号电压V0V1,连续采集F帧(F≥100)取平均值,并计算T0条件下F帧输出信号电压的均方根偏差为该像元的噪声电压Vn

    ③ 计算每个像元的输出信号电压差ΔVs,即V1V0

    ④ 根据Vn、ΔVs、黑体温差ΔT及式(1)计算出器件中每个像元的NETD。

    对于红外焦平面探测器,一般规定的主要测试条件为积分时间或势阱填充量。势阱填充量是指探测器光电转换产生的电子数占读出电路总电荷存储量的百分比,也就是探测器输出信号电压占输出动态范围的百分比,可以通过调整积分时间达到。

    根据普朗克黑体辐射公式,黑体温度为T时在λon~λoff光谱范围内的光子辐亮度为:

    $$L_{\mathrm{p}}(\lambda, T)=\int_{\lambda_{\mathrm{on}}}^{\lambda_{\text {off }}} \frac{2 c}{\lambda^4\left(\mathrm{e}^{\frac{h c}{\lambda k T}}-1\right)} \mathrm{d} \lambda$$ (2)

    式中:Lp(λ, T)为单位面积黑体在单位时间、单位立体角的光子辐亮度,个/(s·sr·cm2);λonλoff为探测器的前后截止波长,μm;c为光速,3×108 m/s;h为普朗克常数,6.63×10-34 Js;k为玻耳兹曼常数,1.38×10-23 J/K。

    面源黑体可以看成一个均匀大面积的朗伯源,当面源黑体充满探测器的整个视场时,探测器光电转换产生的电子数为:

    $$N=\eta A_{\mathrm{d}} t_{\mathrm{int}} \varOmega L_{\mathrm{p}}(\lambda, T)$$ (3)

    式中:N为探测器光电转换产生的电子数,个;η为探测器总量子效率,无量纲;Ad为像元的光敏面积,cm2tint为积分时间,s;Ω为探测器的视场角,sr,Ω=π/(4F2+1),F为探测器冷屏所限制的F数。

    假设探测器的暗电流为Idark,通过读出电路后,探测器的输出信号电压:

    $$V=\frac{Q}{C_{\mathrm{R}}}=\frac{q N+I_{\mathrm{dark}} t_{\mathrm{int}}}{C_{\mathrm{R}}}$$ (4)

    式中:q为电子电量,1.6×10-19 C;CR为读出电路的积分电容,F。

    T0T1黑体温度下,探测器的暗电流基本不发生改变,可得到这两个温差下响应电压ΔVs的计算公式:

    $$\Delta V_{\mathrm{s}}=\frac{\eta q A_{\mathrm{d}} t_{\mathrm{int}} \varOmega\left[L_{\mathrm{p}}\left(\lambda, T_1\right)-L_{\mathrm{p}}\left(\lambda, T_0\right)\right]}{C_{\mathrm{R}}}$$ (5)

    假设在T0温度下其总噪声为Vn(T0),根据式(1)可得到NETD的理论计算公式:

    $$\mathrm{NETD}=\frac{C_{\mathrm{R}} V_n\left(T_0\right)}{\eta q A_{\mathrm{D}} t_{\mathrm{int}} \varOmega\left[L_{\mathrm{P}}\left(\lambda, T_1\right)-L_{\mathrm{P}}\left(\lambda, T_0\right)\right]}\left(T_1-T_0\right)$$ (6)

    为了便于分析,温度T0T1下探测器光电转换产生的电子数记为N0N1,则N0ηAdtintΩLP(λ, T0)、N1ηAdtintΩLP(λ, T1)。探测器的总噪声包括背景辐射、暗电流、读出电路、测试系统等引起的噪声。其中背景辐射到达探测器上光子数目的均方起伏等于其光子到达率的平均值[22],因此T0温度下背景辐射噪声电压Vn-B为:

    $$V_{\mathrm{n}-\mathrm{B}}=q \sqrt{N_0} / C_{\mathrm{R}}$$ (7)

    假设探测器除了背景辐射噪声以外其他噪声电子数为Nde,则:

    $$V_{\mathrm{n}}\left(T_0\right)=\frac{q}{C_{\mathrm{R}}} \sqrt{N_0+N_{\mathrm{de}}^2}$$ (8)

    式(6)可简化为文献[23]中NETD的表达形式:

    $$\mathrm{NETD}=\frac{1}{K \sqrt{N_0}} \sqrt{1+\frac{N_{\mathrm{de}}^2}{N_0}}$$ (9)

    式中:K定义为T0温度下单位温差的光子数变化百分比,$K = ({N_1}/{N_0} - 1)/({T_1} - {T_0})$。

    对式(9)分析可以得到以下结果:

    当暗电流、读出电路等引起的其他噪声电子数Nde相对于背景T0产生的噪声电子数$\sqrt {{N_0}} $较大,即(Nde2/N0)≫1时,NETD主要受限于Nde大小,即为探测器噪声所限制的NETD。

    当(Nde2/N0)≪1时,NETD受探测器的其他噪声影响可以忽略,此时NETD主要受限于N0,即达到背景限的NETD。

    对于焦平面探测器,由于受读出电路积分电容设计尺寸的限制,能处理的电荷数不可能无限制地大,在上述两种情况下,当读出电路最大电荷存储量小于(Nde2+N0)时,NETD则受限于读出电路的电荷存储量,此时需要通过一些措施来减少探测器产生的电子数,例如减小像元尺寸、缩短积分时间、增大F数、降低暗电流等。另一种情况是,读出电路电荷存储量足够但受应用或测试条件的限制,例如要求在固定的势阱填充量条件下进行NETD测试,此时读出电路的电荷存储量越大,NETD就可以越小。

    由于探测器受光子辐照产生的电子数与温度不成比例关系,即使T0相同,选择不同的T1进行测试,K值发生变化,NETD也会有所不同。

    当只考虑背景辐射噪声时,根据式(9)可得到探测器在T0温度下背景限NETD的理论公式:

    $$\mathrm{NETD}=\frac{1}{K \sqrt{N_0}}$$ (10)

    代入KN0,即:

    $$\mathrm{NETD}=\frac{\sqrt{L_{\mathrm{P}}\left(\lambda, T_0\right)}}{\left[L_P\left(\lambda, T_1\right)-L_{\mathrm{P}}\left(\lambda, T_0\right)\right] \sqrt{\eta A_{\mathrm{D}} t_{\mathrm{int}} \Omega}}\left(T_1-T_0\right)$$ (11)

    从式(10)、(11)可以看出,T0T1不变的情况下,背景限NETD与$\sqrt {{N_0}} $成反比例关系,即与量子效率、像元面积、积分时间的平方根成反比例关系,与视场角平方根$\sqrt \varOmega $成反比例关系,即正比于$\sqrt {\frac{{4{F^2} + 1}}{\pi }} $,当F数较大时,与F数成近似正比例关系。

    从式(11)得到的另一个结论是,背景限NETD与读出电路的积分电容无关。实际中,当测试系统或热像仪系统带来的噪声不可忽略时,NETD与积分电容相关,其在他条件不变的情况下,随着积分电容变大,NETD也变大,但此时如果再采取增加积分时间、减小F数等方法达到相同的势阱填充量,由于测试系统或热像仪系统噪声在总噪声中占的比例减小,那么NETD则相对是变小的。

    根据上面的理论公式,NETD与测试温度相关,根据NETD定义,温度T1越接近T0,NETD越接近实际情况,但通常红外焦平面探测器在NETD测试时采用的都是较大温差。理论上,只要积分时间及T0保持不变,通过大温差下NETD的测试结果可推算出小温差下NETD的结果。下面是推算过程:

    T0T1温差下NETD实测结果记为${\rm NETD}_{{({T_1} - {T_0})}}$,T0T2温差下NETD的推算结果记为${\rm NETD}_{{({T_2} - {T_0})}}$,背景T0下探测器的总噪声为${V_{n - {T_0}}}$,在这两个条件下,温差响应信号分别$\Delta {V_{\rm s}}_{({T_1} - {T_0})}$和$ \Delta {V_{\rm s}}_{({T_2} - {T_0})} $。由式(1)有:

    $$ {\rm NETD}_{{({T_1} - {T_0})}} = \frac{{{V_{n - {T_0}}}}}{{\varDelta {V_{\rm s}}_{({T_1} - {T_0})}}}{\text{(}}{T_1} - {T_0}{\text{)}} $$ (12)
    $$ {\rm NETD}_{{({T_2} - {T_0})}} = \frac{{{V_{n - {T_0}}}}}{{\Delta {V_{\rm s}}_{({T_2} - {T_0})}}}({T_2} - {T_0}) $$ (13)

    可以得到:

    $$ {\rm NETD}_{{({T_2} - {T_0})}} = \frac{{({T_2} - {T_0})}}{{({T_1} - {T_0})}}\frac{{\Delta {V_{\rm s}}_{({T_1} - {T_0})}}}{{\Delta {V_{\rm s}}_{({T_2} - {T_0})}}}{\rm NETD}_{{({T_1} - {T_0})}} $$ (14)

    将式(5)代入上式得到T0T2温差下的NETD:

    $$ \begin{array}{l} {\rm NETD}_{{({T_2} - {T_0})}} = \hfill \\ \frac{{({T_2} - {T_0})}}{{({T_1} - {T_0})}}\frac{{[{L_P}(\lambda , {T_1}) - {L_P}(\lambda , {T_0})]}}{{[{L_P}(\lambda , {T_2}) - {L_P}(\lambda , {T_0})]}}{\rm NETD}_{{({T_1} - {T_0})}} \hfill \\ \end{array} $$ (15)

    通过上式,只需计算出两个温差条件下理论光子辐亮度差的比值,即可推算出另一温差下的NETD值。

    选择目前典型640×512 HgCdTe中波红外焦平面探测器,像元面积2.25×10-6 cm2,响应波段范围为3.7~4.9 μm,器件在整个响应范围内的量子效率为0.7,背景黑体温度范围250~350 K,温差取1 K,对其常用条件下NETD进行理论计算及分析。

    1)不同积分时间条件下背景限NETD

    冷屏F数为4,探测器读出电路电荷存储量不受限制时,不同积分时间下背景限NETD的理论计算结果如图 2(a),若读出电路电荷存储量为6.4 Me,实际应用中考虑势阱填充量为10%~90%,背景限NETD实际能达到的值如图 2(b)

    图  2  背景限NETD与背景温度、积分时间的关系
    Figure  2.  BLIP NETD with the background temperature and integration time

    从理论计算结果可以看出,随着背景温度和积分时间的降低,NETD均是增大的,但其前提条件是积分时间或背景温度不变,在较低背景温度下,只要系统可接受的积分时间足够长,NETD理论上可达到较小的值,但事实上,由于系统帧频限制,且由于暗电流的影响,在低背景条件下积分时间也不可能达到很大的值,NETD也会受到一定程度的限制。而高背景温度条件下,由于读出电路电荷存储量的限制,积分时间也不可能很长,需根据实际情况确定,在某些条件下,低背景温度下的NETD可能小于高背景温度下的NETD。从图 2(b)明显看出,当读出电路电荷存储量固定且势阱填充量在一定范围时,NETD受到较大限制,有些条件不可能达到理想状态,图中的标记为“半阱输出”的圆点位置是指某一背景度下,势阱填充量为50%时对应的积分时间和NETD,随着背景温度的升高,达到半阱时需要的积分时间减小,半阱NETD有增大的趋势。在本例中积分时间为1 ms和2 ms时,即使在高温350 K背景条件下也仍达不到半阱输出,此时NETD的值较大。

    2)不同电荷存储量条件下背景限NETD

    仍然是在冷屏F数为4条件下,假设读出电路电荷存储量变化范围为1 Me~12.8 Me,为了使NETD的评价更为客观,在不同背景温度下、调整积分时间使输出均达到半阱状态,NETD的理论计算结果图 3,图中实线和虚线分别为不同条件下半阱NETD和半阱积分时间。背景温度不变的情况下,随着电荷存储量的增加,半阱积分时间增大,半阱NETD减小。图中250 K背景温度下,12.8 Me电荷存储量时,半阱积分时间已超过200 ms,除特殊情况外已不具备实用的条件。

    图  3  背景限半阱NETD与背景温度、电荷存储量、半阱积分时间关系
    Figure  3.  Half well fill BLIP NETD with the background temperature, charge capacity and integration time

    3)不同F数条件下背景限NETD

    冷屏的F数直接决定了入射光子数量,假设F数范围为2~5,背景温度和积分时间不变,只改变F数,根据式(11),F数=5的背景限NETD近似为F数=2背景限NETD的2.5倍,但此时探测器输出势阱填充量差别较大,若调整积分时间使输出均达到半阱状态,在电荷存储量为6.4 Me时,不同背景温度下NETD的理论计算结果图 4。从结果可以看到,随着F数的增加,半阱积分时间增大,但相同背景温度下的NETD保持不变,不同F数下的NETD曲线均重合在一起,只要电荷存储量不变,半阱输出的电子数就已经固定,不论F数是否相同,半阱积分时间下的NETD都保持不变。

    图  4  背景限半阱NETD与背景温度、F数、半阱积分时间关系
    Figure  4.  Half well fill BLIP NETD with the background temperature, F number and integration time

    对典型640×512 HgCdTe中波红外焦平面探测器固定背景温度下的NETD进行实测及推算。取背景温度T0为293 K、F数为4、电荷存储量6.4 Me、积分时间为8 ms时,选择不同的T1,温差从±0.0001~±15 K,背景限NETD的理论计算值如图 5。从图中可以看出,T1T0时,采用的温差越大NETD越小;T1T0时,采用的温差越大NETD越大。在这两种情况下,绝对温差小于1 K时,NETD的变化都很小,绝对温差小于0.1 K时,NETD基本都不再发生变化。

    图  5  不同温差下NETD的变化曲线
    Figure  5.  NETD variation curves under the different temperature difference

    选取一只上述典型器件,采用不同的温差进行NETD实测,按常规测试方法T1T0,积分时间仍为8 ms,并按式(15)进行NETD推算,结果如表 1图 6。由于未扣除器件暗电流、读出电路、测试系统等噪声,NETD实测结果略高于背景限理论计算结果,但不影响推算结果。

    表  1  640×512中波探测器NETD实测结果与推算结果比较
    Table  1.  Comparison between test results and calculation results of 640×512 MW detector
    T0
    /mK
    T1
    /mK
    Tested NETD
    (T0T1)/mK
    Calculation
    NETD/mK
    293 308 21.3 /
    293 303 23.1 23.0
    293 298 24.7 24.8
    293 295 25.8 26.0
    293 294 26.3 26.4
    293 293.5 26.4 26.6
    293 293.2 / 26.7
    293 293.1 / 26.7
    293 293.01 / 26.8
    293 293.001 / 26.8
    293 293.0001 / 26.8
    下载: 导出CSV 
    | 显示表格
    图  6  640×512中波探测器NETD实测结果与推算结果比较
    Figure  6.  Comparison between test results and calculation results of 640×512 MW detector

    当温差小于0.5 K时,由于黑体控温、采集精度等的影响,响应信号基本不再发生变化,NETD测试已无意义。从已测的这组数据来看,随着与293 K温差的逐渐缩小,NETD逐渐增加,通过293~308 K常规测试结果推算出的NETD也呈同样规律,且两者结果非常接近,温差等于1 K时,推算结果和实测结果均比293~308 K常规结果大了约24%,温差小于1 K时,不管是实测结果还是推算结果,NETD都变化很小。由此可以得到如下结论,通过推算可以达到与实测相近的结果,在不改变目前探测器NETD测试规则的前提下,采用1 K温差进行推算,可以得到固定背景温度下更接近实际应用情况的NETD。

    HgCdTe、InSb等光伏型红外焦平面探测器的NETD除了和器件本身的像元面积、光谱响应范围、F数、量子效率、暗电流、读出电路等因素相关以外,也受限于各种实际应用条件,如帧频、背景温度等要求的限制,还与测试条件,如测试温度温差、积分时间等各种因素相关,这些因素相互影响和制约,不同的应用条件和测试条件可能导致NETD不同,需要根据具体情况结合理论计算进行分析。由于目前焦平面探测器常规NETD测试都是采用较大温差进行,与实际情况存在差异,本文推导的理论公式可在固定背景温度下从常规测试的NETD推算出较小温差的NETD,通过典型器件验证,不同温差的推算结果与实测结果符合较好,采用1K温差的NETD推算结果已非常接近实际情况,可为热像仪系统应用提供参考。本文虽以中波探测器为例进行验证,其分析方法适用于所有波长的红外焦平面探测器。

  • 图  1   红外弱小目标跟踪算法流程图

    Figure  1.   Diagram of infrared small target tracking algorithm

    图  2   新生强度控制遮罩

    Figure  2.   New born intensity controlling mask

    图  3   部分红外弱小目标数据集

    Figure  3.   Partial infrared small target datase

    图  4   数据集17图像序列特征

    Figure  4.   Image sequence characteristics of dataset 17

    图  5   数据集17图像复杂度分布

    Figure  5.   Image complexity distribution of dataset 17

    图  6   目标遮罩、目标局部响应与局部峰值

    Figure  6.   Target mask, local response and peak

    图  7   不同场景下的跟踪效果

    Figure  7.   Results about different scenes

    图  8   OSPA误差对比

    Figure  8.   OSPA error comparison

    图  9   基数估计分布对比

    Figure  9.   Comparison of cardinality distribution estimation

    表  1   红外弱小目标数据集信息表

    Table  1   Infrared small target data set information sheet

    Dataset Average SCR SCR Variance Description Type
    4 3.75 3.646 Close range, two targets, sky background, intersecting flight Simple
    6 5.11 1.571 Moving near and far, single target, ground background Simple
    8 6.07 0.159 Moving far and near, single target, ground background Complex
    9 6.29 17.086 Moving near and far, single target, ground background Complex
    13 1.98 0.886 Target moving near and far, single target, weak target, ground background Complex
    15 3.42 0.965 Single target, target maneuvering, ground background Complex
    17 1.09 0.353 Target moving near and far, single target, weak target, ground background Complex
    18 3.32 0.165 Target moving far and near, single target, ground background Complex
    20 3.01 1.485 Single target, target maneuvering, sky background Complex
    22 2.20 0.150 Target moving far and near, single target, ground background Complex
    下载: 导出CSV

    表  2   本文算法在不同数据集下的跟踪结果

    Table  2   Tracking results of the proposed algorithm on different datasets

    Dataset AMD AFA AR PR Rstd, x Rstd, y
    6 8% 2% 95.2% 93.6% 0.43761 0.86444
    8 8% 9.5% 97.2% 89.2% 5.3831 2.6797
    9 10% 11% 98.8% 92.4% 0.8031 0.8671
    13 2% 0% 100% 98.8% 0.5980 0.6135
    15 2% 8% 97.2% 86.4% 1.1927 0.6202
    17 10.8 6% 100% 67.6% 1.2950 0.9641
    18 1.2% 0% 99.6% 99.6% 0.6654 0.5817
    20 0% 0% 100% 100% 0.2274 0.2453
    22 21.6% 28.4% 96.8% 84.8% 1.2324 0.71792
    下载: 导出CSV

    表  3   每个滤波器的平均耗时

    Table  3   Average time cost of every filter  s

    Dataset Filter
    JLMB JLMBAD JGLMB JGLMBAD
    4 3.8125 3.2813 7.625 8.8175
    8 2.2031 1.951 2.625 1.6563
    下载: 导出CSV
  • [1] 张雪, 梁晓庚. 红外探测器发展需求[J]. 电光与控制, 2013, 20(2): 41-45. https://www.cnki.com.cn/Article/CJFDTOTAL-DGKQ201302013.htm

    ZHANG Xue, LIANG Xiaogeng. Development of and demands for infrared detectors[J]. Electronics Optics & Control, 2013, 20(2): 41-45. https://www.cnki.com.cn/Article/CJFDTOTAL-DGKQ201302013.htm

    [2] 赵琰, 刘荻, 赵凌君, 等. 基于YOLOv3的复杂环境红外弱小目标检测[J]. 航空兵器, 2019, 26(6): 29-34. https://www.cnki.com.cn/Article/CJFDTOTAL-HKBQ201906006.htm

    ZHAO Yan, LIU Di, ZHAO Lingjun, et al. Infrared dim and small target detection based on YOLOv3 in complex environment[J]. Aero Weapony, 2019, 26(6): 29-34. https://www.cnki.com.cn/Article/CJFDTOTAL-HKBQ201906006.htm

    [3]

    CHEN C L, LI H, WEI Y, et al. A local contrast method for small infrared target detection[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 52(1): 574-581.

    [4]

    Aghaziyarati S, Moradi S, Talebi H. Small infrared target detection using absolute average difference weighted by cumulative directional derivatives[J]. Infrared Physics & Technology, 2019, 101: 78-87.

    [5]

    FU J, ZHANG H, WEI H, et al. Small bounding-box filter for small target detection[J]. Optical Engineering, 2021, 60(3): 1-14.

    [6] 薛秋条, 宁巧娇, 吴孙勇, 等. 基于JMS-SMC-PHD滤波的检测前跟踪算法[J]. 红外技术, 2020, 42(8): 783-788. http://hwjs.nvir.cn/cn/article/id/hwjs202008013

    XUE Qiutiao, NING Qiaojiao, WU Sunyong, et al. A track-before-detect algorithm based on a JMS-SMC-PHD filter[J]. Infrared Technoiogy, 2020, 42(8): 783-788. http://hwjs.nvir.cn/cn/article/id/hwjs202008013

    [7]

    Bocquel M. Random finite sets in multi-target tracking-efficient sequential MCMC implementation[D]. Enschede: Centre for Telematics and Information Technology, 2013: 21-22.

    [8]

    Kim D Y, Vo B N, Vo B T. A labeled random finite set online multi-object tracker for video data[J]. Pattern Recognition, 2019, 90: 377-389. DOI: 10.1016/j.patcog.2019.02.004

    [9]

    XIAO S, MA Y, FAN F, et al. Tracking small targets in infrared image sequences under complex environmental conditions[J]. Infrared Physics & Technology, 2020, 104: 103102.

    [10]

    TAO H, SHEN X, DENG Q. Infrared target tracking algorithm based on Bernoulli filter and support vector machine[C]//International Conference on Information Science and Education (ICISE-IE). IEEE, 2020: 277-281.

    [11] 董小虎, 傅瑞罡, 高颖慧, 等. 复杂背景下红外小目标自适应检测跟踪[J]. 航空兵器, 2019, 26(6): 22-28. https://www.cnki.com.cn/Article/CJFDTOTAL-HKBQ201906005.htm

    DONG Xiaohu, FU Ruigang, GAO Yinghui, et al. Detection and tracking of small infrared targets adaptively in complex background[J]. Aero Weaponry, 2019, 26(6): 22-28. https://www.cnki.com.cn/Article/CJFDTOTAL-HKBQ201906005.htm

    [12]

    SHU C, DING X, FANG C. Histogram of the oriented gradient for face recognition[J]. Tsinghua Science and Technology, 2011, 16(2): 216-224. DOI: 10.1016/S1007-0214(11)70032-3

    [13]

    Eysa R, Hamdulla A. Issues on infrared dim small target detection and tracking[C]//International Conference on Smart Grid and Electrical Automation (ICSGEA). IEEE, 2019: 452-456.

    [14] 施天俊, 鲍广震, 王福海, 等. 一种适用于多场景的红外弱小目标检测跟踪算法[J]. 航空兵器, 2020, 26(6): 35-42. https://www.cnki.com.cn/Article/CJFDTOTAL-HKBQ201906007.htm

    SHI Tianjun, BAO Guangzhen, WANG Fuhai, et al. An infrared small target detection and tracking algorithm applying for multiple scenarios[J]. Aero Weaponry, 2020, 26(6): 35-42. https://www.cnki.com.cn/Article/CJFDTOTAL-HKBQ201906007.htm

    [15]

    Brown C E. Coefficient of Variation[M]. Berlin: Springer, 1998: 155-157.

    [16]

    Horé A, Ziou D. Image quality metrics: PSNR vs. SSIM[C]//20th International Conference on Pattern Recognition, 2010: 2366-2369.

    [17]

    Dubuisson Séverine. Tracking with Particle Filter for High-dimensional Observation and State Spaces[M]. New York: John Wiley & Sons, Inc., 2015.

    [18]

    Derpanis K G. The bhattacharyya measure[J]. Mendeley Computer, 2008, 1(4): 1990-1992.

    [19]

    Mahler R P S, Vo B T, Vo B N. CPHD filtering with unknown clutter rate and detection profile[J]. IEEE Transactions on Signal Processing, 2011, 59(8): 3497-3513. DOI: 10.1109/TSP.2011.2128316

    [20]

    Reuter S, Vo B T, Vo B N, et al. The labeled multi-Bernoulli filter[J]. IEEE Transactions on Signal Processing, 2014, 62(12): 3246-3260. DOI: 10.1109/TSP.2014.2323064

    [21]

    Do C T, Nguyen T T D, et al. Multi-object tracking with an adaptive generalized labeled multi-Bernoulli filter[J]. Signal Processing, 2022, 196: 1-20.

    [22]

    HUI B W, SONG Z Y, FAN H Q, et al. A dataset for dim small target detection and tracking of aircraft in infrared image sequences [DB/OL]. Science Data Bank, 2022 [2022-09-01]. https://www.scidb.cn/en/.

    [23] 刘洁, 李小昱, 李培武, 等. 基于近红外光谱的板栗水分检测方法[J]. 农业工程学报, 2010, 26(2): 338-341. https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201002059.htm

    LIU Jie, LI Xiaoyu, LI Peiwu, et al. Determination of moisture in chestnuts using near infrared spectroscopy[J]. Transaction of the Chinese Society of Agricultural Engineering, 2010, 26(2): 338-341. https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201002059.htm

    [24]

    Schuhmacher D, Vo B T, Vo B N. A consistent metric for performance evaluation of multi-object filters[J]. IEEE Transactions on Signal Processing, 2007, 56(8): 3447-3457.

    [25]

    HAO Q, GAO H, JUN G. Variational Bayesian labeled multi-Bernoulli filter with unknown sensor noise statistics[J]. Chinese Journal of Aeronautics, 2016, 29(5): 1378-1384. DOI: 10.1016/j.cja.2016.05.002

图(9)  /  表(3)
计量
  • 文章访问数:  92
  • HTML全文浏览量:  13
  • PDF下载量:  50
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-02-22
  • 修回日期:  2023-05-24
  • 网络出版日期:  2024-07-24
  • 刊出日期:  2024-07-19

目录

/

返回文章
返回