基于视轴矢量序列的目标轨迹确认方法

薛永宏, 王铁兵, 乔凯, 樊士伟

薛永宏, 王铁兵, 乔凯, 樊士伟. 基于视轴矢量序列的目标轨迹确认方法[J]. 红外技术, 2022, 44(4): 357-363.
引用本文: 薛永宏, 王铁兵, 乔凯, 樊士伟. 基于视轴矢量序列的目标轨迹确认方法[J]. 红外技术, 2022, 44(4): 357-363.
XUE Yonghong, WANG Tiebing, QIAO Kai, FAN Shiwei. Target Confirmation Method Based on Line of Sight Vector Sequence[J]. Infrared Technology , 2022, 44(4): 357-363.
Citation: XUE Yonghong, WANG Tiebing, QIAO Kai, FAN Shiwei. Target Confirmation Method Based on Line of Sight Vector Sequence[J]. Infrared Technology , 2022, 44(4): 357-363.

基于视轴矢量序列的目标轨迹确认方法

详细信息
    作者简介:

    薛永宏(1985-),男,博士,助理研究员,主要研究方向为空间信息获取与处理、红外探测系统设计

    通讯作者:

    樊士伟(1961-),男,博士,研究员,主要研究方向为天基感知体系

  • 中图分类号: TP911.73

Target Confirmation Method Based on Line of Sight Vector Sequence

  • 摘要: 针对椭圆轨道红外监视系统图像畸变大,基于图像配准进行目标轨迹确认误差大的问题,提出基于视轴矢量序列的目标确认分析方法;建立了目标视轴矢量序列的状态转移模型和观测模型,利用GM-PHD(Gaussian mixture probability hypothesis density)滤波器对目标视轴矢量进行滤波实现目标确认分析;针对在目标状态未知条件下,GM-PHD滤波器易于丢失新目标的问题,采用反馈滤波方式对GM-PHD滤波器进行改进设计。仿真结果验证了基于视轴矢量序列进行目标轨迹确认分析的有效性以及改进GM-PHD滤波器对新目标状态估计的时效性。
    Abstract: Aiming at the serious image distortion of elliptical orbit infrared surveillance system and the low precision of target confirmation algorithm based on image registration, a novel target confirmation method based on line of sight (LOS) vector sequence is proposed. By establishing the state transition model and observation model of target LOS vector sequence, the GM-PHD filter is used to filter LOS vector sequence data and confirm target. Owing to that the GM-PHD filter will lose the new target when its states are unknown, the feedback GM-PHD filter is proposed. The simulation results verify the effectiveness of the proposed method based on LOS vector sequence and the timeliness of the improved GM-PHD feedback filter for new target states estimation.
  • 红外监视系统因其作用距离远、隐蔽性强而被广泛应用于各个领域。典型红外监视系统如美国天基红外系统(space based infrared system, SBIRS),包括地球静止轨道和大椭圆轨道两类卫星。目标轨迹确认是天基红外监视系统关键技术之一,参考文献[1-4]以地球静止轨道卫星为背景,采用“先图像配准、后确认分析”的策略,实现对目标像平面轨迹的确认。此类方法主要以地球静止轨道卫星为应用背景,相邻探测图像相对“静止”、帧间图像变化差异小、易于实现图像配准;当卫星运行于椭圆轨道时,探测图像随时间持续的缩放、旋转[5],帧间图像差异变大,图像配准残差增大,上述方法对目标轨迹确认的错误率将增加。

    以椭圆轨道探测条件下目标轨迹确认为研究背景,针对帧间探测图像畸变大、基于图像配准方法的轨迹确认误差大等问题,提出基于视轴矢量序列的目标轨迹确认方法,并利用GM-PHD(Gaussian mixture probability hypothesis density)滤波器对目标视轴矢量序列进行滤波,实现目标轨迹确认。

    目标视轴矢量是指三维空间中从光学传感器口面中心点到目标点连线所形成的矢量,如图 1所示。红外监视系统中,目标视轴矢量一般通过建立逆成像模型,利用目标在像平面位置、卫星轨道、姿态等参数计算得到。由于探测图像中的目标点既包括真实的目标点,也包括由杂波产生的虚假目标点;实际得到的视轴矢量集是真实目标视轴矢量和虚假目标(杂波)视轴矢量的合集。

    图  1  目标视轴矢量示意图
    Figure  1.  The sketch of target LOS vector

    目标视轴矢量ueci, Tar通常可表示为:

    $$ {{\boldsymbol{u}}_\text{eci,Tar}} = \frac{{{{\boldsymbol{r}}_\text{eci}} - {{\boldsymbol{r}}_\text{eci/sat}}}}{{\left\| {{{\boldsymbol{r}}_\text{eci}} - {{\boldsymbol{r}}_\text{eci/sat}}} \right\|}} $$ (1)

    将公式(1)对时间求导可得到目标视轴矢量运动速度$ {{\boldsymbol{\dot u}}_\text{eci,Tar}} $与目标运动速度$ {{\boldsymbol{\dot r}}_\text{eci}} $、卫星运动速度$ {{\boldsymbol{\dot r}}_\text{eci/sat}} $之间的关系:

    $$ {{\boldsymbol{\dot u}}_\text{eci,Tar}} = \frac{{{\boldsymbol{I}} - {{\boldsymbol{u}}_\text{eci,Tar}}{\boldsymbol{u}}_\text{eci,Tar}^{\text{T}}}}{{\left\| {{{\boldsymbol{r}}_\text{eci}} - {{\boldsymbol{r}}_\text{eci/sat}}} \right\|}}\left( {{{{\boldsymbol{\dot r}}}_\text{eci}} - {{{\boldsymbol{\dot r}}}_\text{eci/sat}}} \right) $$ (2)

    式中:I为3×3单位矩阵。由式(2)可知,目标视轴矢量的运动是卫星运动与目标运动的耦合。

    通常卫星在空间的运动主要受地球重力的影响,可采用J2或J4模型表示[6-7];目标则受推力、重力、空气摩擦力等多种力的作用[8],且不同的力随时间变化的大小不同,因而一般采用动力学模型对其运动特性进行建模[9]。由于目标视轴矢量运动是卫星运动和目标运动的耦合,其运动可看作多种力作用下复杂的变加速运动,如式(3):

    $$ \left\{ \begin{array}{l} {{\boldsymbol{u}}_k} = {{\boldsymbol{u}}_{k - 1}} + {{{\boldsymbol{\dot u}}}_{k - 1}} \cdot \Delta t + \frac{1}{2}{{{\boldsymbol{\ddot u}}}_{k - 1}} \cdot \Delta {t^2} + {{\boldsymbol{q}}_{\text{p}}} \hfill \\ {{{\boldsymbol{\dot u}}}_k} = {{{\boldsymbol{\dot u}}}_{k - 1}} + {{{\boldsymbol{\ddot u}}}_{k - 1}} \cdot \Delta t + {{\boldsymbol{q}}_{\text{v}}} \hfill \\ {{{\boldsymbol{\ddot u}}}_k} = {{{\boldsymbol{\ddot u}}}_{k - 1}} + {{\boldsymbol{q}}_{\text{a}}} \hfill \\ \end{array} \right. $$ (3)

    式中:qpqνqa分别为目标视轴矢量位置、速度、加速度扰动项;Δt为成像周期。同理,杂波在像平面表现为随机运动的位置点,其视轴矢量运动特性符合高斯分布,因此可利用目标与杂波视轴矢量运动特性差异,通过滤波进行目标轨迹确认。

    目标视轴矢量为一条指向特定方向的射线,视轴矢量长度的变化并不会对矢量所包含的信息产生影响。为降低计算复杂性,可对视轴矢量进行归一化降维处理,使得:

    $$ \left\| {{{\boldsymbol{u}}_\text{eci,Tar}}} \right\| = \sqrt {u_{\text{x}}^2 + u_{\text{y}}^2 + u_{\text{z}}^2} = 1 $$ (4)

    归一化后单位视轴矢量ueci, Tar可通过其任意两个坐标轴唯一确定。在椭圆轨道红外监视系统中,卫星轨迹与目标之间的距离较远,在ECI坐标系下,目标视轴矢量的Z轴分量可近似认为是一个常数,因此可利用XOY平面内X轴和Y轴分量的滤波结果唯一确定单位视轴矢量ueci, Tar。降维后公式(3)给出的简化变加速运动模型可表示为:

    $$ \left[ {\begin{array}{*{20}{c}} {{u_{{\text{x}},k}}} \\ {{u_{{\text{y}},k}}} \\ {{{\dot u}_{{\text{x}},k}}} \\ {{{\dot u}_{{\text{y}},k}}} \\ {{{\ddot u}_{{\text{x}},k}}} \\ {{{\ddot u}_{{\text{y}},k}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} 1&0&{\Delta t}&0&{\frac{{\Delta {t^2}}}{2}}&0 \\ 0&1&0&{\Delta t}&0&{\frac{{\Delta {t^2}}}{2}} \\ 0&0&1&0&{\Delta t}&0 \\ 0&0&0&1&0&{\Delta t} \\ 0&0&0&0&1&0 \\ 0&0&0&0&0&1 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{u_{{\text{x}},k - 1}}} \\ {{u_{{\text{y}},k - 1}}} \\ {{{\dot u}_{{\text{x}},k - 1}}} \\ {{{\dot u}_{{\text{y}},k - 1}}} \\ {{{\ddot u}_{{\text{x}},k - 1}}} \\ {{{\ddot u}_{{\text{y}},k - 1}}} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{q_{{\text{px}}}}} \\ {{q_{{\text{py}}}}} \\ {{q_{{\text{vx}}}}} \\ {{q_{{\text{vy}}}}} \\ {{q_{{\text{ax}}}}} \\ {{q_{{\text{ay}}}}} \end{array}} \right] $$ (5)

    目标视轴矢量状态转移模型即可表示为:

    $$ {{\boldsymbol{x}}_k} = {\boldsymbol{F}}{{\boldsymbol{x}}_{k - 1}} + {\boldsymbol{Q}} $$ (6)

    其中,$ {{\boldsymbol{x}}_k}{\text{ = }}{\left[ {\begin{array}{*{20}{c}} {{u_{{\text{x}},k}}}&{{u_{{\text{y}},k}}}&{{{\dot u}_{{\text{x}},k}}}&{{{\dot u}_{{\text{y}},k}}}&{{{\ddot u}_{{\text{x}},k}}}&{{{\ddot u}_{{\text{y}},k}}} \end{array}} \right]^{\text{T}}} $为简化的目标视轴状态矢量;F为状态转移矩阵;Q为过程噪声。

    基于视轴矢量进行目标轨迹确认其观测数据来源于单帧图像目标检测得到的疑似目标点,结合卫星位置、姿态角、传感器参数等解算得到视轴矢量,因此观测模型可构建为:

    $$ {{\boldsymbol{z}}_k} = {\boldsymbol{H}}{{\boldsymbol{x}}_k} + {\boldsymbol{R}} $$ (7)

    式中:$ {\boldsymbol{H}}{\text{ = }}\left( {\begin{array}{*{20}{c}} 1&0&0&0&0&0 \\ 0&1&0&0&0&0 \end{array}} \right) $为观测矩阵;R为观测噪声。目标视轴矢量序列的状态转移模型和观测模型都为线性模型,因此可采用高斯混合概率假设密度(GM-PHD)滤波器[10-12]对目标轨迹进行确认。

    假设目标状态转移模型和测量模型均为线性高斯模型,即:

    $$ \left\{ \begin{array}{l} {f_{k|k - 1}}\left( {{\boldsymbol{x}}|{\boldsymbol{\zeta}}} \right) = N\left( {{\boldsymbol{x}};{{\boldsymbol{F}}_{k - 1}}{\boldsymbol{\zeta}},{\boldsymbol{Q}_{k - 1}}} \right) \hfill \\ {g_k}\left( {{\boldsymbol{z}}|{\boldsymbol{x}}} \right) = N\left( {{\boldsymbol{z}};{\boldsymbol{H}_k}{\boldsymbol{x}},{\boldsymbol{R}_k}} \right) \hfill \\ \end{array} \right. $$ (8)

    式中:N(•;m, P)表示均值为m,方差为P的高斯分布函数;Fk-1为状态转移矩阵;Qk-1为过程噪声协方差矩阵;Hk为观测模型矩阵;Rk为观测噪声协方差阵。

    假设在k-1时刻,目标状态函数为:

    $$ {v_{k - 1}}\left( {\boldsymbol{x}} \right) = \sum\limits_{i = 1}^{{J_{k - 1}}} {\omega _{i,k - 1}^{}N\left( {{\boldsymbol{x}};{\boldsymbol{m}}_{i,k - 1}^{},{\boldsymbol{P}}_{i,k - 1}^{}} \right)} $$ (9)

    则在k时刻,预测的目标状态函数为:

    $$ {v_{k|k - 1}}\left( {\boldsymbol{x}} \right) = v_{k|k - 1}^{\left( {\text{S}} \right)}\left( {\boldsymbol{x}} \right) + v_{k|k - 1}^{\left( {\text{B}} \right)}\left( {\boldsymbol{x}} \right) + {\gamma _k}\left( {\boldsymbol{x}} \right) $$ (10)

    式中:$ v_{k|k - 1}^{\left( {\text{S}} \right)}\left( {\boldsymbol{x}} \right) $、$ v_{k|k - 1}^{\left( {\text{B}} \right)}\left( {\boldsymbol{x}} \right) $和$ {\gamma _k}\left( {\boldsymbol{x}} \right) $分别表示从上一时刻存留的目标、分裂的目标以及新生成的目标(即新目标),三者的状态函数为:

    $$ \left\{ \begin{array}{l} v_{k|k - 1}^{\left( {\text{S}} \right)}\left( {\boldsymbol{x}} \right) = p_k^{\left( {\text{S}} \right)}\sum\limits_{i = 1}^{{J_{k - 1}}} {\omega _{i,k - 1}^{}N\left( {{\boldsymbol{x}};{\boldsymbol{m}}_{i,k|k - 1}^{\left( {\text{S}} \right)},{\boldsymbol{P}}_{i,k|k - 1}^{\left( {\text{S}} \right)}} \right)} \hfill \\ v_{k|k - 1}^{\left( {\text{B}} \right)}\left( {\boldsymbol{x}} \right) = \sum\limits_{i = 1}^{{J_{k - 1}}} {\sum\limits_{l = 1}^{J_k^{\left( {\text{B}} \right)}} {\omega _{i,k - 1}^{}\omega _{l,k}^{\left( {\text{B}} \right)}N\left( {{\boldsymbol{x}};{\boldsymbol{m}}_{i,l,k|k - 1}^{\left( {\text{B}} \right)},{\boldsymbol{P}}_{i,l,k|k - 1}^{\left( {\text{B}} \right)}} \right)} } \hfill \\ {\gamma _k}\left( {\boldsymbol{x}} \right) = \sum\limits_{i = 1}^{J_k^{\left( {\rm{\Gamma }} \right)}} {\omega _{i,k}^{\left( {\rm{\Gamma }} \right)}N\left( {{\boldsymbol{x}};{\boldsymbol{m}}_{i,k}^{\left( {\rm{\Gamma }} \right)},{\boldsymbol{P}}_{i,k}^{\left( {\rm{\Gamma }} \right)}} \right)} \hfill \\ \end{array} \right. $$ (11)

    其中:

    $$ \begin{array}{l} \left\{ \begin{array}{l} {\boldsymbol{m}}_{i,k|k - 1}^{\left( {\text{S}} \right)} = {\boldsymbol{F}_{k - 1}}{\boldsymbol{m}}_{i,k - 1}^{\left( {\text{S}} \right)} \hfill \\ {\boldsymbol{P}}_{i,k|k - 1}^{\left( {\text{S}} \right)} = {\boldsymbol{Q}_{k - 1}} + {\boldsymbol{F}_{k - 1}}{\boldsymbol{P}}_{i,k - 1}^{\left( {\text{S}} \right)}{\left( {{\boldsymbol{F}_{k - 1}}} \right)^T} \hfill \\ \end{array} \right. \hfill \\ \left\{ \begin{array}{l} {\boldsymbol{m}}_{i,l,k|k - 1}^{\left( {\text{B}} \right)} = \boldsymbol{F}_{l,k - 1}^{\left( {\text{B}} \right)}{\boldsymbol{m}}_{i,k - 1}^{\left( {\text{B}} \right)} + {\boldsymbol{d}}_{l,k - 1}^{\left( {\text{B}} \right)} \hfill \\ {\boldsymbol{P}}_{i,l,k|k - 1}^{\left( {\text{B}} \right)} = \boldsymbol{Q}_{l,k - 1}^{\left( {\text{B}} \right)} + \boldsymbol{F}_{l,k - 1}^{\left( {\text{B}} \right)}{\boldsymbol{P}}_{i,k - 1}^{\left( {\text{B}} \right)}{\left( {\boldsymbol{F}_{l,k - 1}^{\left( {\text{B}} \right)}} \right)^T} \hfill \\ \end{array} \right. \hfill \\ \end{array} $$ (12)

    式中:Jk(B)Jk(Γ)分别表示分裂目标和新目标的个数;ωl, k(B)ωl, k(Γ)为对应高斯项的权重;pk(S)为目标的存留概率;dl, k-1(B)Pi, k-1(B)Ql, k-1(B)为分裂目标状态的先验信息。可以看出,预测的目标状态函数仍可写为高斯混合形式:

    $$ {v_{k|k - 1}}\left( {\boldsymbol{x}} \right) = \sum\limits_{i = 1}^{{J_{k|k - 1}}} {w_{i,k|k - 1}^{}N\left( {{\boldsymbol{x}};{\boldsymbol{m}}_{i,k|k - 1}^{},{\boldsymbol{P}}_{i,k|k - 1}^{}} \right)} $$ (13)

    k时刻,传感器的量测为zk,则更新后的目标状态函数为:

    $$ {v_k}\left( {\boldsymbol{x}} \right) = \left( {1 - {p_{\text{D}}}} \right){v_{k|k - 1}}\left( {\boldsymbol{x}} \right){\text{ + }}\sum\limits_{{\boldsymbol{z}} \in {{\boldsymbol{Z}}_k}} {v_k^{\left( {\text{D}} \right)}\left( {{\boldsymbol{x}},{\boldsymbol{z}}} \right)} $$ (14)

    其中,第1项表示漏检目标的PHD;第2项表示利用检测到的目标更新后的PHD;且:

    $$ \left\{ \begin{array}{l} v_k^{\left( {\text{D}} \right)}\left( {{\boldsymbol{x}},{\boldsymbol{z}}} \right) = \sum\limits_{i = 1}^{{J_{k|k - 1}}} {w_{i,k}^{}\left( {\boldsymbol{z}} \right)N\left( {{\boldsymbol{x}};{\boldsymbol{m}}_{i,k|k}^{}\left( {\boldsymbol{z}} \right),{\boldsymbol{P}}_{i,k|k}^{}} \right)} \hfill \\ w_{i,k}^{}\left( {\boldsymbol{z}} \right) = \frac{{{p_{\rm{D}}}w_{i,k|k - 1}^{}q_{i,k}^{}\left( {\boldsymbol{z}} \right)}}{{{\kappa _k}\left( {\boldsymbol{z}} \right) + {p_{\rm{D}}}\sum\nolimits_{j = 1}^{{J_{k|k - 1}}} {\omega _{j,k|k - 1}^{}q_{j,k}^{}\left( {\boldsymbol{z}} \right)} }} \hfill \\ \end{array} \right. $$ (15)

    GM-PHD滤波器目标状态更新滤波示意图如图 2所示。

    图  2  GM-PHD滤波器目标状态更新
    Figure  2.  Target state update of GM-PHD filter

    根据GM-PHD滤波器状态更新过程,通常新目标的状态信息仅包含在状态函数γk(x)中,将式(11)代入目标状态更新函数式(14)后,可将新目标状态更新公式简写为:

    $$ v_k^{\left( {\rm{\Gamma }} \right)}\left( {\boldsymbol{x}} \right) = \sum\limits_{{\boldsymbol{z}} \in {{\boldsymbol{Z}}_k}} {\sum\limits_{n = 1}^{J_{k - 1}^{\left( {\rm{\Gamma }} \right)}} {\omega _k^{\left( {\rm{\Gamma }} \right)}\left( {{\boldsymbol{z}},n} \right)N\left( {{\boldsymbol{x}};{\boldsymbol{m}}_{n,k}^{}\left( {\boldsymbol{z}} \right),{\boldsymbol{P}}_{n,k}^{}} \right)} } $$ (16)

    其中:

    $$ \omega _k^{\left( {\rm{\Gamma }} \right)}\left( {{\boldsymbol{z}},n} \right) = \frac{{{p_{\rm{D}}}\omega _{n,k}^{\left( {\rm{\Gamma }} \right)}N\left( {{\boldsymbol{z}};{g_k}\left( {{\boldsymbol{m}}_{n,k}^{\left( {\rm{\Gamma }} \right)}} \right),{\boldsymbol{S}}_{n,k|k - 1}^{}} \right)}}{{{\kappa _k}\left( {\boldsymbol{z}} \right) + {p_{\rm{D}}}\sum\nolimits_{j = 1}^{{J_{k|k - 1}}} {\omega _{j,k|k - 1}^{}q_{j,k}^{}\left( {\boldsymbol{z}} \right)} }} $$ (17)

    表征更新后新目标高斯项的权重。为减少计算资源的消耗,GM-PHD滤波器即依据该权重对高斯项进行裁剪。以图 3为例,式(17)中高斯项权重的更新,本质上为真实量测z和预测量测gk(mn, k(Γ))之间相关性的计算过程。当新目标位置等先验信息已知时,真实量测与估计的量测在测量空间的位置非常接近,相关性大,因而高斯项权重较大,目标不易于被剪除;反之当目标位置等先验信息未知时,二者在测量空间的距离一般较远,相关性小,高斯项权重也非常小,甚至接近于零,新目标易于被剪除而导致丢失。

    图  3  GM-PHD测量空间似然函数
    Figure  3.  Likelihood function of GM-PHD measurement

    可以看出,GM-PHD滤波器导致新目标丢失的一个重要原因就是假设的新目标状态与真实目标状态不符,进而导致在测量空间中,真实量测与估计量测距离较远所致。事实上,任何新目标都将会在传感器中产生量测,且由于目标运动速度有限,在相邻时间内测量空间中新目标量测距离将很近。因而,利用相邻时间历史测量状态直接假设为新目标量测状态,即将k-1时刻传感器量测作为k时刻新目标量测状态的反馈滤波方式,可有效解决新目标丢失问题。新目标高斯项权重更新方式如下:

    $$ \omega _k^{\left( {\rm{\Gamma }} \right)}\left( {{\boldsymbol{z}},n} \right) = \frac{{{p_\rm{D}}\omega _{n,k}^{\left( {\rm{\Gamma }} \right)}N\left( {{\boldsymbol{z}};{{\boldsymbol{z}}_{k{\text{ - }}1}},{\boldsymbol{S}}_{n,k|k - 1}^{}} \right)}}{{{\kappa _k}\left( {\boldsymbol{z}} \right) + {p_\rm{D}}\sum\nolimits_{j = 1}^{{J_{k|k - 1}}} {\omega _{j,k|k - 1}^{}q_{j,k}^{}\left( {\boldsymbol{z}} \right)} }} $$ (18)

    改进后的GM-PHD反馈滤波器目标状态更新滤波示意图如图 4所示。

    图  4  GM-PHD反馈滤波器目标状态更新
    Figure  4.  Target state update of GM-PHD feedback filter

    以一颗运行于Molniya轨道的大椭圆轨道红外监视卫星为例,假定多目标场景。场景中包含3个不同时刻,从不同位置发射的导弹目标,目标1从仿真开始即出现,目标2从仿真开始后40 s出现,目标3从仿真开始后80 s出现。仿真中假设杂波密度κ(z)为10-5;成像周期Δt为5 s,像元角分辨率为60 μrad;卫星姿态确定误差为10'';经视轴矢量校正后,视轴矢量偏差为60 μrad。利用目标在像平面运动轨迹,解算后即可得到目标视轴矢量序列,对目标视轴矢量进行归一化并降低其维数后,得到序列目标视轴矢量在ECI坐标系XOY平面的运动轨迹。目标运动轨迹如图 5所示。

    图  5  目标运动轨迹
    Figure  5.  Target trajectory

    分别采用GM-PHD滤波器(GM-PHDF)和改进的GM-PHD滤波器(IGM-PHDF)对目标轨迹进行滤波;通过比较估计目标个数与真实目标个数的差异、估计目标视轴矢量与真实目标视轴矢量的偏差以及OSPA(optimal subpattern assignment)距离[13-14]对两个算法的性能进行评价。试验中Monte-Carlo仿真次数为1000次,两种滤波器对目标视轴矢量XY的估计结果以及Z轴解算结果如图 6所示;两种滤波器对目标个数估计结果、平均视轴矢量偏差及OSPA距离如图 7所示。

    图  6  两种滤波器对目标视轴矢量的估计结果
    Figure  6.  LOS estimation results of two different filters
    图  7  两种滤波器对目标状态的估计结果
    Figure  7.  Target state estimation results of two filters

    仿真结果表明:①基于GM-PHD滤波器,利用序列目标视轴矢量进行目标确认分析是可行的,可实现对目标个数以及目标视轴矢量状态的估计;②改进后的GM-PHD滤波器可在新目标出现后5 s内得到所有目标的视轴矢量信息,从而及时地发现新目标;与之相比,传统的GM-PHD滤波器可以在5 s内得到第一个目标和第二个目标的视轴矢量信息,但对于第三个目标则需在15 s后才能得到其视轴矢量信息,发现第三个目标的所需时间多于改进的算法。

    进一步比较两种滤波器对目标数量估计结果、平均视轴矢量偏差以及OSPA距离可知:待发现目标后,两种滤波器对目标视轴矢量信息的估计精度相当;但改进后的GM-PHD滤波器对目标数量的估计结果更为准确。

    目标视轴矢量序列运动特性与目标像平面轨迹运动特性相同,利用目标视轴矢量序列对目标轨迹进行确认分析,可有效解决椭圆轨道监视系统图像畸变大、基于图像配准的确认分析方法误差大等问题;改进设计的GM-PHD滤波器可在目标状态未知条件下,快速、有效、准确地对新目标状态进行估计,避免了新目标丢失和确认时效性差等问题,从而实现目标快速确认。尽管论文所提反馈滤波器是在高斯混合模型下进行的,但是通过采用PHD的EKF、UKF或粒子滤波实现形式,也可推广到非线性、非高斯模型应用的场景之中。

  • 图  1   目标视轴矢量示意图

    Figure  1.   The sketch of target LOS vector

    图  2   GM-PHD滤波器目标状态更新

    Figure  2.   Target state update of GM-PHD filter

    图  3   GM-PHD测量空间似然函数

    Figure  3.   Likelihood function of GM-PHD measurement

    图  4   GM-PHD反馈滤波器目标状态更新

    Figure  4.   Target state update of GM-PHD feedback filter

    图  5   目标运动轨迹

    Figure  5.   Target trajectory

    图  6   两种滤波器对目标视轴矢量的估计结果

    Figure  6.   LOS estimation results of two different filters

    图  7   两种滤波器对目标状态的估计结果

    Figure  7.   Target state estimation results of two filters

  • [1] 司马端, 龙云利, 安玮等. 基于显著性与帧间差分的红外弱小目标检测[J]. 航天电子对抗, 2015, 31(4): 32-35. DOI: 10.3969/j.issn.1673-2421.2015.04.009

    SIMA duan, LONG Yunli, An Wei, et al. Infrared small target detection based on saliency and frame to frame difference[J]. Aerospace Electronic Warfare, 2015, 31(4): 32-35. DOI: 10.3969/j.issn.1673-2421.2015.04.009

    [2]

    SUN Yang, YANG Jungang, LI Miao, et al. Infrared small-faint target detection using non-i. i. d. mixture of Gaussians and flux density[J]. Remote Sensing, 2019, 11(23): 2831. DOI: 10.3390/rs11232831

    [3]

    GAO Jinyan, LIN Zaiping, AN Wei. Infrared small target detection using a temporal variance and spatial patch contrast filter[J]. IEEE Access, 2019, 7(99): 32217-32226.

    [4]

    ZHU Ran, LONG Yunli, AN Wei. Multitarget tracking based on PHD smoother with unknown clutter spatial density[J]. Progress in Electromagnetics Research C, 2018, 82: 123-133. DOI: 10.2528/PIERC17120408

    [5] 张涛, 薛永宏, 丁永和, 等. 背景在椭圆轨道红外监视平台的成像特性分析[J]. 红外与激光工程, 2016, 45(3): 282-287. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201603044.htm

    ZHANG Tao, XUE Yonghong, DING Yonghe, et al. Analysis of imaging characteristics of background in elliptical orbit infrared surveillance platform[J]. Infrared and Laser Engineering, 2016, 45(3): 282-287. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201603044.htm

    [6] 肖业伦. 航空航天器运动的建模: 飞行动力学的理论基础[M]. 北京: 北京航空航天大学出版社, 2003.

    XIAO Yelun. Modeling of Aerospace Motion: Theoretical Basis of Flight Dynamics[M]. Beijing: Beihang University Press, 2003.

    [7] 张雅声, 徐艳丽, 杨庆. 航天器轨道理论与应用[M]. 北京: 清华大学出版社, 2020.

    ZHANG Yasheng, XU Yanli, YANG Qing. Spacecraft Orbit Theory and Application[M]. Beijing: Tsinghua University Press, 2020.

    [8] 张毅, 肖龙旭, 王顺宏. 弹道导弹弹道学[M]. 长沙: 国防科技大学出版社, 2005.

    ZHANG Yi, XIAO Longxu, WANG Shunhong. Ballistic Missile Ballistics[M]. Changsha: National University of Defense Technology Press, 2005.

    [9] 盛卫东. 天基光学监视系统目标跟踪技术研究[D]. 长沙: 国防科技大学, 2011.

    SHENG Weidong. Research on Target Tracking Technologies for Space-based Optical Surveillance System[D]. Changsha: National University of Defense Technology, 2011.

    [10]

    LUO Qingqing, GAO Zhisheng, XIE Chunzhi. Improved GM-PHD filter based on threshold separation clusterer for space-based starry-sky background weak point target tracking[J]. Digital Signal Processing, 2020, 103: 102766. DOI: 10.1016/j.dsp.2020.102766

    [11]

    Kim Dohyeung, Kwon Cheolhyeon, Hwang Inseok. Gaussian mixture probability hypothesis density filter against measurement origin uncertainty[J]. Signal Processing, 2020, 171: 107448. DOI: 10.1016/j.sigpro.2019.107448

    [12]

    ZHANG Qian, Taek Lyul Song. Improved bearings-only multi-target tracking with GM-PHD filtering[J]. Sensors, 2016, 16(9): 1469. DOI: 10.3390/s16091469

    [13]

    Dominic Schuhmacher, Ba-Tuong Vo, Ba-Ngu Vo. A consistent metric for performance evaluation of multi-object filters[J]. IEEE Trans. on Signal Processing, 2008, 56(8): 3447-3457. DOI: 10.1109/TSP.2008.920469

    [14] 刘伟峰, 文成林. 随机集多目标跟踪性能评价指标比较与分析[J]. 光电工程, 2010, 37(9): 14-20. https://www.cnki.com.cn/Article/CJFDTOTAL-GDGC201009006.htm

    LIU Weifeng, WEN Chenglin. Performance comparison and analysis of multitarget tracking indices based on random set[J]. Opto-Electronic Engineering, 2010, 37(9): 14-20. https://www.cnki.com.cn/Article/CJFDTOTAL-GDGC201009006.htm

图(7)
计量
  • 文章访问数:  150
  • HTML全文浏览量:  47
  • PDF下载量:  24
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-09
  • 修回日期:  2021-01-09
  • 刊出日期:  2022-04-19

目录

/

返回文章
返回