本发明涉及一种运动目标的联合跟踪与识别方法,尤其是基于雷达宽带测量的弱机动目标跟踪与识别一体化处理方法。
背景技术:
在传统的运动目标跟踪与识别方法中,目标的识别和跟踪往往被视为两个独立的环节,其方法是先对目标进行识别,在此基础上对感兴趣目标维持跟踪。这种方法将识别作为一次性决策,其判断失误可能导致整个过程产生错误;同时,由于识别和跟踪环节不存在交联和反馈,识别结果也无法给跟踪环节提供有效的信息。为了提高对目标识别和跟踪的一体化处理能力,文献1:a.donka,m.lyudmila,“jointtargettrackingandclassificationwithparticlefilteringandmixturekalmanfilteringusingkinematicradarinformation”,digitalsignalprocessing,2006,16:180-204(a.donka等人于2006年在《数字信号处理》期刊第16卷发表的论文,题目为“利用雷达运动学信息和粒子滤波/混合卡尔曼滤波进行联合目标跟踪与分类”,称为背景技术1)公开了一种基于目标运动学测量信息(包括目标的距离、角度)的联合目标跟踪与识别处理方法。该方法从不同类型的空中飞机目标具有不同的机动特性(如加速度大小、速度包线等)出发,通过目标运动参数(如:位置、速度和加速度)的估计实现目标跟踪与识别的联合处理。该方法不足表现为:1)该方法中目标的识别主要通过对高速运动目标的机动模式判断来实现,在雷达探测过程中,若目标不执行机动飞行动作,识别条件无法得到满足;2)对于弱机动目标,如慢速运动大型舰船等,由于其在机动能力和速度包线等方面不存在显著的差异性,该方法不适用;3)此方法基于传统的窄带雷达,只能对应点目标,由于受分辨率条件的限制,只能获取目标的信号幅度或雷达截面积(radarcrosssection,rcs)等简单观测,无法得到关于目标的更多细节信息;且rcs严重敏感于目标姿态,在不同的观测角度下其大小起伏变化不定,通过目标rcs进行目标识别非常困难。
宽带雷达可以实现目标在纵向距离维度上的分辨力,由此得到的观测结果称为高分辨距离像(highrangeresolutionprofile,hrrp),相应地,此时的目标被称为扩展目标(表示目标占据多个距离分辨单元)。hrrp比rcs提供了更多关于目标的细节信息,在进行目标分类识别时可以摆脱对目标机动特性的依赖;但hrrp本身同样面临姿态敏感问题,在有多种混淆目标的条件下,不同类型目标在不同姿态下的hrrp特征可能具有高度的相似性,造成模式分类决策的模糊性,从而影响最终的目标识别效果。
因此,如何充分利用目标的运动学测量与高分辨距离像测量信息,提高目标跟踪与识别精度,是雷达系统整体性能提升中面临的一个重要问题。
解决该问题的一种可行思路是对扩展目标的运动状态和类型参数进行联合估计(滤波),其中,对目标的类型估计通过目标距离像(识别特征)匹配来实现,这就要求系统在联合处理过程中具备对目标距离像的预测能力。文献2:jieminhu,weiwang,qinglinzhai,jianpingou,ronghuizhan,andjunzhang,“globalscatteringcenterextractionforradartargetsusingamodifiedransacmethod”,ieeetransactionsonantennasandpropagation,2016,64(8):3573-3586(jieminhu等人于2016年在《天线与传播汇刊》第64卷第8期发表的论文,题目为利用修正的ransac方法进行雷达目标的全局散射中心提取,称为背景技术2)公开了一种雷达三维散射中心模型(见第3581页图8),该模型是目标几何实体在电磁域的抽象和等效,但没有涉及如何用三维散射中心模型进行距离像预测和目标识别。
技术实现要素:
本发明提出一种基于模型的联合跟踪与识别方法,该方法以雷达三维散射中心模型为基础,通过对目标运动状态的预测得到目标方向角(即雷达视线与目标轴向的夹角)的估计,并以此目标方向角作为输入条件,产生距离像的预测值,将该距离像预测值与雷达实际获取的距离像测量进行匹配,得到识别结果。由于跟踪过程的连续性,随着目标运动状态估计误差的不断减小和目标类型估计不确定性的不断降低,目标的识别概率也在不断提高,最终得到稳定识别结果。
本发明提供的基于模型的联合跟踪与识别方法主要包括以下几个步骤:
第一步,对目标运动状态与传感器观测进行建模,方法是:
1.1建立目标的运动状态模型:
对在x-y二维平面内作近似匀速运动的目标建立如公式(1)所示的状态演化模型:
xk=fxk-1+wk(1)
式中,
其中,ts为目标采样(观测时间)间隔,通常在0.5s~5s范围内取值;q为噪声强度,对于地海面目标其取值通常不大于1m/s2。
1.2建立雷达传感器对目标的观测模型:
1.2.1建立目标的运动学观测模型,如公式(3)所示:
式中,h(xk)为xk的运动学观测方程;βk表示k时刻包含了观测噪声的目标方位,rk表示k时刻包含了观测噪声的距离信息;vk=[v1,kv2,k]t,
1.2.2建立目标的距离像观测模型,如公式(4)所示
dk=g(φk,s)+nk=idft{ei=e(fi,φk,γk,s),i=0,1,2,…,i}+nk(4)
式中,g(φk,s)表示距离像观测方程,φk为k时刻目标方向角(雷达视线方向与目标运动方向的夹角),γk表示k时刻目标的俯仰角,当目标在二维平面内运动时γk≈0,s={an,αn,(xn,yn,zn);n=1,...,n}为目标的三维散射中心模型,idft表示逆离散傅立叶变换,e(·)为频率响应函数,ei表示fi在e(·)中的响应结果,i为频点个数,nk为i维零均值高斯观测噪声矢量。e(fi,φk,γk,s)的具体形式可用公式(5)表示为
其中fi为第i个频点对应的信号频率,fc为雷达载频,c为光速,j为虚数单位,an表示第n个散射中心的幅度,(xn,yn,zn)表示第n个散射中心在三维空间中的坐标,αn为第n个散射中心的类型参数,n为散射中心数目。
对于特定的目标,其模型参数(xn,yn,zn)、an和αn在使用之前已完全确定,这些参数可通过文献2中的三维散射中心模型重构处理过程(见文献2的第3579页图6)得到。
fi在[fc-b/2,fc+b/2]区间范围内,通过公式(6)来选取
式中,b为信号带宽,一般在几十到几百mhz量级,i的取值范围通常在几十到几百之间。
方向角φk通过公式(7)来确定
由公式(3)和(4)可得到k时刻目标运动学和距离像的联合观测{zk,dk},以及从τ=0时刻到τ=k时刻为止目标运动学和距离像的联合观测集
第二步:初始化,方法为:
2.1对目标状态进行初始化,方法是:
若实际应用条件下可能出现
式中,
2.2对k=0时刻第
式中,
2.3取q为小于1的正数,令ts为小于5的正数,一般取q=(0.1~0.5)m/s2,ts=(0.5~5)s,对公式(2)中的q和f进行初始化。
2.4取b=(50~500)mhz,i=50~500,对公式(6)中的fi(i=0,1,2,…,i)进行初始化,得到i个频点,即f0,…,fi,…,fi均按公式(6)初始化为相同的值。
第三步,令k=1。
第四步,从雷达接收k时刻的观测数据zk,预测k时刻的目标运动状态模型和观测模型,方法是:
令k-1时刻第
4.1令
4.2预测第
4.3预测第
4.4预测第
4.5预测第
式中,
4.5计算第
4.6用
式中,
4.7计算第
4.8判定
4.9计算第
4.10对第
第五步,估计k时刻目标的状态和类别,实现联合跟踪与识别,方法是:
5.1计算k时刻
5.2计算k时刻
5.3获取k时刻最终的目标状态估计
第六步,找到
本发明的有益效果主要有:
1)采用本发明可利用连续跟踪信息有效提高目标的正确识别概率,即随着时间下标k的不断推进,目标状态估计的精度越来越高,由公式(12)得到的目标方向角预测精度也逐渐提高,通过公式(15)得到的与真实目标类型相对应的粒子滤波器似然也相应增大,从而提高公式(20)中的目标分类识别概率);
2)突破了传统点目标跟踪与识别联合处理方法中实现目标分类识别需依赖于目标机动动作的限制(即传统中,点目标的识别只能通过对目标的机动模式估计来实现),本发明第四步通过公式(14)、(15)和(16)实现了基于宽带雷达目标运动学和距离像观测的联合预测,并最终实现了目标状态与类型的联合跟踪与识别,大大扩展了联合目标跟踪与识别技术的适用范围。
附图说明
图1是本发明总体流程图;
图2是三类舰船目标几何模型及散射中心位置;图2(a)是目标a的几何模型与对应的三维散射中心模型,图2(b)是目标b的几何模型与对应的三维散射中心模型,图2(c)是目标c的几何模型与对应的三维散射中心模型
图3是图2(a)目标a测试结果;图3(a)是图2(a)目标a的跟踪结果,图3(b)是图2(a)目标a的分类识别结果
图4是图2(b)(目标b)测试结果;图4(a)是图2(b)目标b的跟踪结果,图4(b)是图2(b)目标b的分类识别结果
图5是图2(c)(目标c)测试结果。图5(a)是图2(c)目标c的跟踪结果,图5(b)是图2(c)目标c的分类识别结果
具体实施方式
下面结合附图和实施例,对本发明进行详细解释和说明。
图1为本发明总体流程图,假定雷达传感器固定在坐标原点,对海面慢速运动目标进行观测,将连续观测数据与粒子滤波算法相结合,实现对目标航迹的跟踪和对目标类型的识别,具体实施过程如下:
第一步,对目标运动状态与传感器观测进行建模,方法是:在二维平面中建立如公式(1)所示的目标运动模型,并分别根据公式(3)和(4)产生目标的运动学观测结果和距离像观测结果。由公式(3)和(4)得到k时刻目标运动学和距离像的联合观测{zk,dk},以及从τ=0时刻到τ=k时刻为止目标运动学和距离像的联合观测集
第二步,在k=0时刻,通过公式(8)对
目标初始位置为(x0,y0)=(1.2,1.5)km,初始速度为
第三步,令k=1;
第四步,从雷达接收k时刻的观测数据zk,预测k时刻的目标运动状态模型和观测模型,方法是:对于k-1时刻的第
第五步,估计k时刻目标的状态和类别,实现联合跟踪与识别,方法是:利用公式(19)和(20)分别得到3个滤波器的目标状态估计
第六步,判断k时刻3个粒子滤波器得到的目标识别概率估计值
图2给出了三类目标的几何模型及与几何模型对应的三维散射中心模型,其中黑色星号“*”表示散射中心所处的位置。三类目标分别称为目标a、目标b和目标c,通过这些三维散射中心模型可以得到不同观测视角(方向角)下的距离像。
由文献2可知,目标的三维散射中心模型是通过目标的全姿态、宽带电磁仿真数据(称为原始宽带数据)重构得到的,作为一种电磁等效手段,由三维散射中心模型根据特定姿态角预测得到的距离像与原始宽带数据直接成像得到的距离像存在一定的误差,实验中将这一误差建模为(4)中的观测噪声n。根据目标与雷达的实际场景(相对位置和视线方向),对每一类目标分别进行联合跟踪与识别测试(即实际出现的目标是三类目标中的一个,其它两类目标作为潜在的、可能出现的混淆目标),用
分别假定真实出现的目标类型为目标a、目标b和目标c,对联合跟踪与识别方法进行测试,得到三类不同目标的航迹跟踪结果分别如图3~图5中的(a)部分所示,相应的目标识别结果(分类概率)如图3~图5中(b)部分所示。
如图3所示,在该实验中真实存在的目标类型为目标a,想定可能的目标类型还有目标b和目标c。在初始时刻,三类目标的识别概率均设为1/3,随着跟踪与识别过程的推进,目标a的分类概率(见图3(b)中的菱形线)依次变为0.36,0.48,0.62,0.68,0.79,…。而其他两类目标(分别为圆形线和星型线)的识别概率却逐渐减低。由此可见,在初始化的基础上,经过5个时间周期以后(对应第6个时间点)即可达到目标类型判定门限要求(门限为0.75,实际识别概率为0.79),且识别结果(a类)与真实的目标类型(a类)完全相符。由图3(a)可见目标初始位置坐标为(1200m,1500m),经过5个时间周期,目标真实位置依次变为(1193m,1505m),(1185m,1509m),(1177m,1513m),(1170m,1517m),(1164m,1522m),而得到的位置估计依次为(1192m,1506m),(1185m,1510m),(1174m,1516m),(1169m,1517m),(1165m,1523m),说明对目标轨迹的跟踪也是非常精确的。由图3(b)还可看出,达到目标识别门限之后,随着跟踪与识别过程的持续进行,目标识别概率还可进一步提高,至第25个时间点(对应25s处),其分类概率可达到1(100%),说明识别准确率为100%。
图4给出的真实目标为图2(b)所示b类目标的联合跟踪与识别实验结果,图5给出的是真实目标为图2(c)所示c类目标的联合跟踪与识别实验结果,其过程和结论与图3类似,在此不再赘述。
作为对比,表1同时给出了直接利用三维散射中心模型恢复的距离像与原始宽带数据成像所得的距离像进行相关匹配得到的识别结果(即不考虑跟踪估计环节,只考虑单纯的目标分类识别环节)。该实验中,三类目标的原始宽带距离像数据作为距离像识别模板,方位向0~360°全覆盖,角度间隔为1°(对应360个训练数据样本),俯仰方向取为0°角;通过目标三维散射中心模型(预测)产生距离像数据时,在0~360°全方位姿态范围内均匀采样,角度间隔为0.2°,这样每类目标各产生1440个测试数据样本。表中,目标识别混淆矩阵主对角线上的元素表示每类目标正确识别的距离像样本数,单类目标的识别率(pcc)和三类目标的整体识别率(oa-pcc)分别定义为
其中,
表1
由表1中的结果可以看出,如果不采用目标跟踪与识别联合处理方法,仅考虑单一的目标识别环节(称为常规识别方法),得到a、b、c三类目标的识别概率分别为1294/1440=0.8986,1295/1440=0.8993,1212/1440=0.8417,平均识别率为0.8799。通过对比可以清楚地看到,当粒子滤波器达到稳定之后,本发明联合跟踪与识别方法(识别概率为1)较目前现有常规识别方法(识别概率为0.8799)的识别概率提高0.1(10%)以上。
综上,本发明中的联合跟踪与识别方法不仅有效可行,且在目标分类识别方面较现有常规方法所具有的明显性能优势。