基于广义概率假设密度的多目标运动估计方法研究.docx

基于广义概率假设密度的多目标运动估计方法研究 余萌 徐琰珂 胡茄乾摘 要： 面向空对空作业的发展需求， 提出了一种基于广义概率假设密度的多目标运动估计方法 在Faster-RCNN方法基础上引入多尺度分析， 并利用改进K-means方法对观测目标进行粗聚类， 以此为前置输入， 提出了基于广义泊松分布的概率假设密度滤波器， 将聚类信息纳入滤波估计的权重更新中， 增强了对变阵群目标的跟踪时效性 仿真结果表明， 本文方法在没有初始聚类信息的先验知识下， 依然能够完成对多目标的识别分类与跟踪， 且精度优于现有的集群目标运动估计方法关键词：多目标运动估计; K-means聚类; 目标识别; 概率假设密度; 态势感知：TJ765; TP391.4 ： A ：1673-5048（2021）04-0037-060 引 言为适应日趋复杂的作战环境， 多飞行器协同作战方案逐渐成为当前空中打击作业的热点方向 美国海军研究的低成本无人机蜂群（Drone Swarm）技术已经成功完成了多次实地测验， 初步实现了无人机集群编队飞行与快速机动 俄罗斯新的S-70“猎人”（Okhotnik-B）重型隐形无人机成功伴飞苏-57， 突出了战斗支援角色， 旨在扩大战斗机的能力。

这些无人机集群和无人机与有人机的编队， 大大加强了空中打击力量的突防能力， 为国土防空增加了新的压力 为实现对空中集群目标的有效防御， 需要发展具有集群目标针对性的跟踪应对技术通常情况下， 多飞行器集群目标都是按照一定的编队飞行方式进行空中作业， 通过在复杂背景中进行机动目标集群聚类分析， 可以实现对集群目标的阵群态势粗估计 传统的无监督聚类方法（如谱聚类、 K-means等）容错性差， 且对初始值非常敏感， 尤其对有噪声的机动目标聚类效果欠佳 本文在Faster-RCNN方法基础上引入多尺度分析， 并利用改进K-means方法对观测目标进行了粗聚类， 将集群目标的态势粗估计分为群目标识别以及类别数估计两部分， 旨在提升集群目标运动估计精度的同时， 增强估计方式的灵活性在群目标识别方面， 本文以Faster-RCNN为基础进行改进 经典R-CNN[1]采用了Selective Search算法来提取图像中可能存在的兴趣区域（Regions of Interest）， 之后对每个提取区域采用标准的卷积神经网络（CNN）进行分类 在R-CNN后， 2015年提出的Fast-RCNN[2]采用兴趣区域池化来共享计算量较大的部分， 以此提高模型的效率， 随后发展的Faster-RCNN[3-4]是建立在Fast-RCNN基础上的一个完全可微分模型， 在卷积层后使用Region Proposal Network（RPN）生成一系列靶块， 对其进行裁剪过滤后利用softmax判断靶块属于前景还是后景， 同时使用另一分支进行靶块修正， 凝练较为精确的提议（Proposal）。

Faster-RCNN与Fast-RCNN的明显区别就是使用了RPN代替了选择性搜索来产生提议， 并强调提议的CNN网络与目标识别的CNN共享 Faster-RCNN同时将特征抽取、 候选区域提取、 边界框回归、 分类整合到了一个网络当中， 使目标识别的综合性能得到了很大改进分类态势粗估计可以为多目标跟踪提供前置输入 传统的多目标跟踪方法有联合数据关联方法（JPDA）[5]、 假设目标跟踪方法（MHT）[6-7]等， 这些方法多利用先验假设， 将集群目标出现的几类运动情景综合考虑（如遮挡、 虚报或目标合并等） 考虑到目标以集群为单位运动时， 目标的数量、 阵群会发生变化， 而传统方法对阵群可变的目标运动估计解释性不足[8]， 本文基于随机有限集合思想的扩展目标概率假设密度（ET-PHD）滤波方法[9]， 面向阵群可变的多目标运动估计提出相应的方案 核心方法是利用先验的粗聚类信息为ET-PHD方法提供聚类初值， 并引入广义泊松分布来增强对概率强度估计的解释性， 进而提高对阵群可变的目标跟踪品质 航空兵器 2021年第28卷第4期余 萌， 等： 基于广义概率假设密度的多目标运动估计方法研究1 多目标编队的态势粗估计1.1 基于Faster-RCNN的多目标识别多目标识别是多目标运动估计的前置环节， 能够为后续态势感知、 运动估计提供基礎数据。

考虑应用背景， 本文选择深度学习方法——Faster-RCNN， 对组网结构进行量体裁衣， 发展了基于Faster-RCNN的集群目标识别方法， 使用的敏感器为机载单目相机1.1.1 数据库选择与组建选择的深度学习训练集为DOTA（A Large-Scale Dataset for Object Detection in Aerial Images）[10]， 是武汉大学与华中科技大学联合制作的遥感图像数据集， 共有2 806张遥感图像以及15个类别 经过适当修改， 该数据库最终包含了直升机、 无人机、 客机、 战斗机以及怀疑目标这五个类别， 其中在怀疑目标类别中归属了诸如热气球、 翼伞、 鸟类等空中目标1.1.2 基于改进Faster-RCNN的空中多目标识别考虑空中作业环境中前/后景差异相比地面场景中更为明显， 特征的颗粒度更高， 因此， 从计算代价角度考虑， 并不需要生成一系列稠密的靶块进行前/后景判断 为此， 本文对Faster-RCNN方法进行了改进， 在原图尺度上采用谱残差方法替换原方法中候选靶块生成， 利用多个尺度下的局部图像显著性表征， 通过非极大抑制实现候选框生成， 如图1所示。

这种方法在前/后景差异明显的情况下可以在保证识别精度的同时， 大幅降低训练所需计算代价 在获得Proposal的特征图后， 分类时采用全连接层与softmax计算每个Proposal的类别归属， 并输出51的概率向量 同时， 利用边框回归得到每个Proposal的位置偏移量， 实现更准确的目标检测 由于本文对数据库进行了重新设计， 因此对整体组网进行了重新训练， 并未使用流行的迁移学习方法利用Faster-RCNN能够实现基于单目相机的目标识别， 但考虑方法目的与实时性， 在此仅将其作为集群目标的初始信息， 其余集群运动估计及编队态势提取与此环节独立 图2为使用Faster-RCNN的飞行目标检测示意图利用改进Faster-RCNN进行算法测试， 结果表明， 基于显著性多尺度RPN表征的简化方法与原Faster-RCNN并无明显差异， 但改进方法的训练成本大幅降低， 在未来出现未编目飞行目标时， 重训练以及参数调整的灵活性优势明显1.2 多目标隶属估计通过Faster-RCNN可以输出目标的群类别以及在图像中的矩心位置， 其中， 群类别是以概率向量来描述， 一个典型的输出结果如下：Proposal={[k1， m1]， [k2， m2]， …， [kN， mN]}（1）式中： k=[p1， p2， …， pK]T为概率向量， [pi|i∈K]为属于其中某类的概率， K为类别总个数; m=[μ， ν]为群目标在图像中的矩心位置， 此处为目标所属检测框的矩心。

以Proposal提供的目标识别为基础， 结合K-means方法[11]提出了目标快速聚类方法， 其主要思想是在K-means方法的基础上引入分类评估置信率函数， 即wj=∑x∈Djln（Djx-mj）∑k∈K∑x∈Dkln（Dkx-mk）（2）式中： Dj为第j个类别中的目标聚类集; wj为相对应的权重参数 在K-means迭代过程中， 根据置信率评估结果剔除置信过小的分类， 算法的伪代码实现流程图如图3所示相比经典K-means方法， 本文方法的计算复杂度从O（itn）上升至O（kitn） 该方法能够通过置信率实时估计来增强对观测噪声的抗干扰能力， 相比经典方法在鲁棒性上有一定提升， 图4所示为观测噪声影响情况下的估计结果示例红色棱柱为在编队在相机像平面下的跟踪估计结果， 机载相机的观测结果以粒子点描述 将运动目标分为3个编队， 并在不同信噪比的条件下比较经典K-means与Proposal-K-means两种方法的估计结果Proposal-K-means群目标隶属估计方法 P-K-means（K， M， x， D）输入： （k1， m1）， （k1， m2）D1k2， m3D2， …， kK， mMDV∥Proposal 提供的分类及像素位置For （k， m）∈（K， M）针对每个Proposal：选择K个聚类点作为初始聚类中心： mi， i=1， …， D重复初始化sj=0， j=1， …， D， nj=0， j=1， …， DFor x∈Dj=argminfdist（x， mf）分配相应的x到类别j：sj=sj+x; nj=nj+1;Endmj=sj/nj， j=1， …， k直至mj迭代收敛。

对每一个分类评估其置信率：wj=∑x∈DjlnDjx-mj/∑k∈K∑x∈DklnDkx-mk剔除wj过小的分类输出聚类结果： K-， M-图4（a）中编队数目为3个， 由图可知， 经典K-means能够准确跟踪编队目标的像素位置信息， 仅在仿真末段编队1与编队2较为接近时， 出现6处估计精度欠佳的情况（规定估计结果与编队内任一成员的像素欧式距离大于60即为精度欠佳） 本文所提出的Proposal-K-means在无编队数目这一先验信息情况下， 与经典K-means的估计结果无显著差异 由图4（b）可以看出， 即便在已知编队数量为3个的前提下， 估计精度欠佳的比例仍达到了31%， 主要归因于当观测噪声与目标的像素位置比较接近时， 仅采用欧式距离作为判据法并不能做出准确的区分 相比经典K-means方法， Proposal-K-means在群矩心位置的估计精度上提升明显， 如图4（c）所示， 其估计精度欠佳的比例为18%（相比经典K-means降低了40%）， 其主要是因为Proposal-K-means将编队的隶属信息作为输入， 并通过实时监测分类的置信率来剔除隶属信息不可靠的估计结果， 即相比经典K-means引入了审查环节， 因此能够在一定程度上减少估计精度欠佳的分类。

图4（d）其群矩心位置估计与真实情况具有较为明显的差异， 聚类精度仅为60%， 此结果也符合仿真预期， 因为K-means本质上仍对噪声抗干扰能力不足， 需要使用集群运动估计方法来提升目标运动状态估计的鲁棒性由于Proposal能够提供给K-means分类以及矩心的初值， 故将经典K-means从有监督学习转变为无监督学习， 更加符合任务对自主性与灵活性的需求2 多目标运动估计ET-GM-PHD方法面向的是单个目标可产生多观测量的跟踪场景 在本文的应用背景下， 若将单阵群视为扩展目标（即能够产生多个观测量的目标）， 便可以使用ET-GM-PHD滤波方法进行集群目标運动估计 PHD滤波方法与传统滤波方法的主要区别是去除了计算量庞大的测量信息与运动目标进行匹配关联的步骤， 采用多目标随机有限集合的概率分布一阶矩进行状态估计， 将定义在多目标状态空间上的后验概率密度映射到单目标的状态空间上， 并将单目标状态空间分布定义为概率假设密度（Probability Hypothesis Density， PHD）， 又称目标强度[12] 由于在目标强度分布的迭代估计过程涉及到集合积分的运算， 其算力需求过大， Vo等人提出了GM-PHD方法[13]， 其利用加权高斯混合模型刻画该目标的强度分布， 并通过与Kalman滤波器相似的预测步与更新步操作， 实现时间上的传递， 显著降低了PHD算。