基于随机有限集的视频SAR多目标跟踪方法_陈李田.pdf

第 卷 第期空军工程大学学报 年月 收稿日期：基金项目：国家自然科学基金（）作者简介：陈李田（），男，河南范县人，硕士生，研究方向为信息与通信工程。：引用格式：陈李田，张云，李宏博，等基于随机有限集的视频 多目标跟踪方法空军工程大学学报，（）：，（）：基于随机有限集的视频 多目标跟踪方法陈李田，张云，李宏博，王勇（哈尔滨工业大学电子与信息工程学院，哈尔滨，）摘要视频合成孔径雷达具有高分辨与高帧速率成像的特点，可以连续获取地面感兴趣区域目标近似视频的信息，为基于 图像的目标识别与跟踪技术的快速发展奠定了基础。为了满足日益复杂的应用需求，多目标跟踪技术逐渐发展成熟，针对多目标跟踪过程中每个运动目标的状态都具有空时变性，并且目标的数量具有随机性的难题，首先建立了基于随机有限集的多目标跟踪算法，在此基础上讨论了贝叶斯框架下的概率假设密度算法，并在高斯混合模型下研究并实现了高斯混合概率假设密度滤波算法，进而实现了基于 的多目标跟踪算法，在复杂环境背景下验证了该算法的有效性。关键词视频；随机有限集；多目标跟踪；高斯混合概率假设密度 中图分类号 文献标志码文章编号 （），（，），（），（），（），（），；随着合成孔经雷达（，）技术的不断发展，图像目标自动识别、目标自动跟踪相关研究逐渐成为焦点。近年来，高分辨合成孔径雷达相关技术得到快速发展，并衍生出了高分辨视频 技术。视频 在对地侦测领域具有不可比拟的优势，它能够以近似光学的方式获取地面感兴趣区域的持续观测，更容易获取位于监测场景中的各种目标的相关信息，对基于 图像的目标识别与跟踪的研究起到了积极的推动作用。为了满足更高的应用需求，多目标跟踪技术逐渐发展成熟。在多目标跟踪过程中，每个目标的运动状态都有可能变化，且目标数量具有随机性，因此多目标跟踪需要在对监测区域内的目标状态进行预测和更新的同时也获取目标数量的变化情况，这对跟踪算法的性能提出更高的要求。针对上述需求，基于随机有限集（，）的视频 多目标跟踪技术成为跟踪领域的研究热点。其基本原理是在视频 系统中，利用有限集统计学将 的概念引入视频 多目标跟踪算法中，用有限集表示目标的状态集与量测集，并通过贝叶斯框架估计下的高斯混合概率假设密度（，）等次优滤波方式实现状态估计和更新，并对目标数量做出了预测。该方法将目标的状态集合用 表示，本身就可以对状态不确定和数量不确定的目标集合进行表示，不涉及传统算法的关联步骤，在复杂背景环境下具备较好的跟踪性能。多目标 跟踪算法最早是由 所提出的，典型代表分别为矩近似的 滤波和 滤波算法以及密度近似的多伯努利滤波算法。基于随机有限集的多目标跟踪算法近年来得到了充分了发展。文献 设计了一种基于随机有限集的新的概率假设密度（，）滤波器，通过对目标航迹的调整重组，有效提高了提取多目标航迹的准确度，进而提高了雷达在复杂背景下的持续跟踪性能。文献 提出了一种适用于 系统的基于 的改进高斯混合 滤波器的动目标跟踪方法，可以实现准确估计运动目标数量和状态，并实现运动目标轨迹的提取。本文在随机有限集条件的背景下，通过对 等次优滤波方式的研究，提出了一种适用于视频 连续多帧序列的基于随机有限集的多目标跟踪方法。多目标贝叶斯框架估计理论 基于 的视频 多目标跟踪系统多目标跟踪不仅需要估计当前时刻各目标的运动状态，还需要关注目标数量。由 基本概念可知，有限集中包含的元素数量是随机值，与多目标跟踪中目标数量未知的特点十分契合，所以在 的基础上构建多目标跟踪系统。基于 的视频 多目标跟踪系统模型如图所示。图视频 多目标贝叶斯跟踪系统模型 多目标运动模型和观测模型在目标跟踪的过程，每一时刻目标都有可能发生变化，存在种变化的可能性，分别是目标运动，目标消失，目标新生以及目标衍生。对目标被观测到的情况也分为种：观测到目标信息（似然）、未观测到目标信息（漏检）和观测到错误信息（杂波）。目标运动模型和观测模型示意图如图所示。图目标运动模型和观测模型示意图在本文构建的基于 的多目标跟踪系统中，主要考虑对线性运动目标的监测跟踪，且该类目标自身状态大概率会出现存活、消失和新生种情况，但不排除因为目标交会分离等原因出现目标衍生的情况。此外，在观测过程中，由于相干斑噪声等影响，会出现漏警或者虚警的情况，即目标观测过程中存在似然、漏警和杂波种情况。上述的目标类型、目标状态集合以及目标观测集合共同组成了本文所搭建的系统模型。多目标贝叶斯滤波原理多目标贝叶斯滤波是基于 的多目标跟踪理论的重要部分之一，基本思路是需要在递推出当前时刻多目标的后验概率密度的基础上，通过对后空军工程大学学报 年验概率密度进行积分估计出当前时刻监测区域内目标的运动状态和数量。整个过程的预测和更新递推方程分别为：（：）（：）（）（）（：）（）（：）（）（：）（）式中：（）表示状态转移概率密度函数；（）表示目标似然函数；：，表示直到时刻的观测积累。从预测和更新递推方程式中可以发现，基于 的多目标贝叶斯滤波涉及到对不确定维数集合的积分，随着目标状态 包含元素数量的增长，递推计算量会呈现出数值爆炸的趋势。因此，该方法很难在工程中实现，需要通过对基于多目标 的贝叶斯滤波算法进行近似，在确保跟踪性能的前提下使算法变得较为容易实现。基于 的 滤波算法 基于随机有限集的 滤波算法 是后验概率分布的一阶统计矩，也称之为强度函数。假设（：）是多目标状态有限集的后验概率密度函数，则其 可以表示为：（：）（）（：）（）式中：（）为狄拉克函数。在监测区域内对（：）进行积分就可以得到区域内期望的目标数量的估计值（），即目标数量估计方程式为：（）（：）（）为了简化多目标状态空间结构中所有不确定维数的多维复杂积分，可以通过引入 ，利用多目标强度函数（：）替换滤波递推过程中的概率密度，将多维复杂积分简化为单目标状态空间结构上的单一积分，有效提高滤波效率。具体流程如下：首先在 分布的假设条件下由式（）和式（）可以得到 滤波算法的预测方程和更新方程分别为：（：），（）（）（：）（）（：）（）（）（：）（，（）（：），（）（）（：）（），（）（）（：）（）式中：，（）为目标运动的概率；，（）为目标检测概率；（）和（）分别为从衍生和新 生 目 标 随 机 集 的 ；为 杂 波 随 机 集的。在得 到 当 前 时 刻 多 目 标 的 强 度 函 数（：）之后，由式（）可以估计出当前时刻监测区域内目标的数量：（：）（）基于随机有限集的 滤波有很多优点，比如：可以不通过数据关联便对多个目标的运动状态进行估计；可以对数量未知的多目标集合进行跟踪，自适应监测区域内目标数量的变化等。文献 指出 滤波是建立在 分布假设上的，这就导致了在递推迭代方程式中仍然存在多维的向量积分，要得到 滤波闭合解的解析形式还是比较困难的，并且容易引起维数灾难，但可以通过 滤波解决上述问题。基于随机有限集的 滤波为了优化 的滤波效果，本文采用通过建立高斯混合（，）模型的方法，在 滤波递推过程中用 参数替代 进行递推迭代，可以得到闭合解，进而优化滤波效果。对于线性的多目标跟踪系统，每一个目标的运动模型的状态转移函数和传感器的观测模型的似然函数均满足模型：（）（；，）（）（）（；，）（）式中：（）表示 分布；是前一时刻目标的状态；是状态转移矩阵；是观测矩阵；和分别为内部噪声和外部噪声的协方差矩阵。新生目标和衍生目标的 也需要满足 模型：（），（），（；（），（），）（）（），（），（；（）（），（），）（）式中：，为 模型下新生目标的属性，共同描 述 了新生目标的 。同样的，第期陈李田，等：基于随机有限集的视频 多目标跟踪方法，和，描 述 了 衍 生 目 标的。此时，在模型下时刻多目标的预测强度函数为：（），（），（）（）（）式中：，（），（），（）分别是目标运动、衍生和新生的预测。对于，（），其预测公式为：，（）.（）；（），；（），（）（）（），（）（）（），（）（）同样，（）和（）的预测公式为：，（），（）（），（；（，），；（，），）（）（），（），（；（），；（），）（）经过对，（），（）和（）的计算之后，可以得到高斯混合模型下时刻多目标的预测 为：（）（）（；（）；（）（）然后利用时刻传感器获得的模型下的量测对多目标的（）进行更新。时刻多目标 后验概率密度的强度函数为：（）（，）（），（；）（），（；）是传感器量测 的强度函数，其 形式的表达式为：，（；）（）（）（；（）；（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）综上，可 得时 刻 多 目 标 后 验 概 率的：（）（）（；（）；（）（）式中：，实验结果与分析为验证 本 文 提 出 的 基 于 随 机 有 限 集 的 视 频 多目标滤波跟踪方法的有效性，通过实测组视频 图像序列，进行跟踪验证，其中关键参数如表所示。表视频 系统参数参数数值工作频率 机载速度（）波长 中心斜距 场景大小（）脉冲重复时间 带宽 采样频率 实验共选取连续的 帧 图像序列，其中图是该图像序列中间的某一帧（配准后），作为场景示意图说明场景中的目标运动信息。图场景及目标运动状态示意图如图所示，这是本文实验所用视频 连续多帧图像序列中的某一帧，可以看出这是一个十字路口，在图示时刻的附近几帧中共有个目标，分别标记为目标、目标、目标和目标，在该配准坐标系下目标、沿轴负半轴方向运动，目标沿着轴正半轴方向运动。其中，目标从第帧已经存在到第 帧左右消失；目标在这 帧中一直存在；目标从第帧左右出现后一直存在；目标从第 帧左右出现后一直存在。基于此，结合目标检测所得目标位置信息，可以先估计出目标的运动轨迹，如图所示。需要注意的是，本次实验是定性分析，使用的是像素点位表示的目标位置，约为 像素。图目标运动轨迹空军工程大学学报 年然后进行基于随机有限集的 滤波实现这个目标的跟踪，需要得到检测区域的量测信息以及目标的新生强度信息，其中量测是基于阴影检测算法提取的阴影信息；目标的新生强度则是通过目标阴影的位置信息以及大小所占量测总像素点数量的比值近似估计目标的强度信息。为了方便观测到每个目标的跟踪情况，本实验对各目标进行了标签处理，实验结果见图。（）目标轴跟踪结果（）目标轴跟踪结果图多目标跟踪结果在图中，不同颜色的直线代表的是图中所估计的不同目标的运动轨迹，相应颜色实心点表示跟踪到的相应目标，灰色叉点表示杂波或虚警。然后进行误差分析。第帧、第 帧还有第 帧存在误差的原因是因为上述讨论目标数量的起始和消失帧数有所偏差，与跟踪算法的性能无关；第 帧时丢失了对目标的跟踪，第 帧丢失了对目标的跟踪，但由于算法在不断的预测与更新，在下一帧又重新对目标实现了跟踪；第 帧、第 帧和第 帧出现了错误的新生或衍生判断，由于目标当时位置与杂波位置较近，出现了误判现象，随着算法的更新修正，自动矫正了该错误，说明该跟踪方法具有一定的自纠正能力；第 帧第 帧，由于目标到了图像边缘，杂波影响较大，所以出现了明显的误差，该跟踪方法在边缘处理上有待改进。图给出了每个目标的跟踪情况。图每个目标的跟踪情况如图所示，判别值为说明该帧正常跟踪到了该目标，判别值为说明该帧没有该目标或者该目标消失了，判别值为说明跟踪时判定该目标发生了衍生或出现了其他新生目标。根据图所示情况，可以进一步得到本次实验对多目标目标跟踪估计的准确率约为 。为了更好地展示该跟踪算法的性能，本实验通过最优子模式分配距离（，）对目标的 滤波估计结果集合的脱靶距离进行测量，结果如图所示。图跟踪结果的 距离从图可以看出，该算法的 距离均小于 像素点，即小于，一定程度上说明了本算法的可行性与有效性。最后在成像场景下展示该跟踪算法的效果，仍用红色表示目标，蓝色表示目标，绿色表示目标，由于图像背景较暗，目标改用白色表示，结果如图所示。图跟踪结果展示第期陈李田，等：基于随机有限集的视频 多目标跟踪方法结语本文探索了基于随机有限集的多目标跟踪算法在视频 系统中的应用。在基于 的多目标跟踪系统结构和目标模型下，分析了多目标贝叶斯的原理；然后为了实现基于 的多目标跟踪算法，讨论了贝叶斯框架下的 算法，并在 模型下得到了适用性较强的 算法。最后在上述讨论的基础上，结合一组视频 连续多帧图像序列，实现了对该组图像序列中个目标的跟踪，验证了基于 的多目标跟踪算法在视频 连续多帧图像序列中实现多目标跟踪的可行性。参考文献丁金闪视频 成像与动目标阴影检测技术雷达学报，（）：，（），：，：左峰视频合成孔径雷达成像算法研究成都：电子科技大学，：何盟基于随机有限集的雷达多目标跟踪方法研究及软件实现 西安：西安电子科技大学，：，：，（）：，（）：，：，（）：孙海鑫基于随机有限集的雷达多目标跟踪算法研究 西安：西安电子科技大学，：，（）：陈迪高斯混合 滤波器理论与应用研究哈尔滨：哈尔滨工程大学，：（编辑：徐敏）（上接第 页），（）：，（）：李应红，吴云，梁华，等等离子体激励气动力学探索与展望力学进展，（）：，（）：，：，：，（）：，：，：，（）：化为卓，李应红，牛中国，等低速三角翼纳秒脉冲等离子体激励实验航空动力学报，（）：白涛，王晋军 后掠角三角翼不同前缘形状对其气动 特 性 影 响 的 实 验 研 究 实 验 流 体 力 学，（）：，（）：，（）：，（）：龙玥霄等离子体激励下分离流变化特性及机理研究西安：西北工业大学，：（编辑：姚树峰）空军工程大学学报 年