CN109521420A - 基于多特征匹配的多目标跟踪方法 - Google Patents

基于多特征匹配的多目标跟踪方法 Download PDF

Info

Publication number
CN109521420A
CN109521420A CN201811577747.8A CN201811577747A CN109521420A CN 109521420 A CN109521420 A CN 109521420A CN 201811577747 A CN201811577747 A CN 201811577747A CN 109521420 A CN109521420 A CN 109521420A
Authority
CN
China
Prior art keywords
target
moment
tracking
matrix
currently
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201811577747.8A
Other languages
English (en)
Other versions
CN109521420B (zh
Inventor
水鹏朗
廖沫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201811577747.8A priority Critical patent/CN109521420B/zh
Publication of CN109521420A publication Critical patent/CN109521420A/zh
Application granted granted Critical
Publication of CN109521420B publication Critical patent/CN109521420B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/66—Radar-tracking systems; Analogous systems
    • G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
    • G01S13/726—Multiple target tracking

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

一种基于多特征匹配的多目标跟踪算法，其步骤是：首先，针对量测点与目标关联时，匹配特征数量过少的问题，建立了所有目标的特征矩阵，特征矩阵包含了所有目标的检验统计量特征和多普勒频率特征；其次，在寻找最佳关联点时，将所有关联点的特征与目标的特征矩阵进行比对，找出与目标匹配度最高的关联点；最后，确定目标在当前跟踪时刻的运动状态估计值、检验统计量和多普勒频率特征，用于对目标在当前跟踪时刻之后运动状态的估计。本发明通过建立特征矩阵，并寻找与目标特征矩阵最匹配的关联点，解决了当多目标轨迹交叉重叠时目标轨迹跟踪不准确的问题。

Description

基于多特征匹配的多目标跟踪方法
技术领域
本发明属于雷达技术领域，更进一步涉及雷达目标跟踪技术领域中的一种基于多特征匹配的多目标跟踪方法。本发明可用于复杂场景下地面雷达对于多个目标运动轨迹的跟踪。
背景技术
多目标跟踪一直是雷达技术领域的重要问题之一。在雷达技术领域中，多目标跟踪技术的主要任务是从雷达接收机接收到的连续时间的多帧数据中检测出多目标的运动轨迹，进而在每一帧数据中分别确定这些目标的位置。多目标跟踪的核心部分是数据关联和状态估计。由于实际观测场景和待观测目标较为复杂，存在很多不确定因素，易导致多目标数据关联出现模糊，使得跟踪航迹容易拉偏，因此复杂场景下的多目标跟踪技术具有重要的研究意义。传统的用于数据关联的方法为最近邻域法，该方法容易受到其他目标和强噪声、杂波的影响，造成关联精度的显著下降，从而导致目标跟踪的准确性大大降低。
西北工业大学在其申请的专利文献“一种基于模糊聚类的概率最近邻域多目标跟踪方法”(专利申请号201310751067.4，申请日2013.12.30，申请公布号CN103729859A，申请公布日2014.04.16)中公开了一种利用基于最近邻分类器和均值漂移的目标跟踪方法。该方法实现的具体步骤是，(1)根据k-1时刻目标t的状态估计及其协方差矩阵得到k时刻目标t的状态预测值；(2)计算有效量测与目标t之间的最大关联概率；(3)对于M个传感器，反复执行步骤1和步骤2，依次求得M个传感器的各有效量测j与目标t之间的最大关联概率，得到各传感器观测下的相对应的M个局部状态估计值及其相应的估计均方误差阵；(4)依据步骤3中求得的各传感器的局部状态估计及其相应的估计均方误差阵，对各传感器局部航迹信息进行加权平均求和，得到状态融合更新，估计均方误差融合更新，得到所有目标的全局状态融合估计，实现多目标跟踪。该方法存在两点不足之处：其一，当多个目标的航迹发生交叉或重叠时，当前正在跟踪目标的量测点、其它目标的量测点和虚警点同时落入待跟踪目标在当前时刻跟踪波门的概率大大提升，由于该方法仅根据跟踪波门内所有量测点的坐标来求出这些量测点与目标之间的关联概率，所以这种只通过单一特征进行关联的方法有较大几率导致属于其它目标的量测点或者虚警点与当前跟踪目标关联，从而使得对当前跟踪目标运动状态的估计和未来轨迹的跟踪出现偏差。其二，该方法利用隶属度计算公式计算所有有效量测点与目标轨迹的关联概率，计算方法较为复杂，计算量较大，降低了工程实现效率。
发明内容
本发明的目的在于针对上述现有技术的不足，提出一种基于多特征匹配的多目标跟踪方法，用于改善多个目标的航迹发生交叉时，当前正在跟踪目标的量测点、其它目标的量测点和虚警点同时落入待跟踪目标在当前时刻跟踪波门的情况下雷达量测点的关联精度，从而提升了对每个目标轨迹的跟踪准确性。
实现本发明目的的思路是：首先，针对量测点与当前跟踪目标关联时，匹配特征数量过少的问题，建立了所有目标的特征矩阵，特征矩阵包含了所有目标的历史检验统计量特征和历史多普勒频率特征；其次，在寻找最佳关联点时，将所有关联点的特征与当前跟踪目标在当前时刻的特征矩阵进行比对，从而找出与当前跟踪目标的匹配度最高的关联点；最后，确定每个目标在当前跟踪时刻的运动状态估计值、检验统计量和多普勒频率特征，用于对每个目标在当前跟踪时刻之后运动状态的估计。
本发明的具体步骤如下：
(1)从当前跟踪时刻存在的多个目标中，选择一个未跟踪的目标，作为当前跟踪目标；
(2)预测运动状态：
将当前跟踪目标在当前跟踪时刻上一时刻的运动状态的估计值，代入状态方程，得到当前跟踪目标在当前跟踪时刻运动状态的预测值；
(3)建立当前跟踪目标的特征矩阵：
(3a)寻找与当前跟踪时刻最近的5个时刻；
(3b)从当前跟踪目标的历史数据中，提取当前跟踪目标在这5个时刻的检验统计量和多普勒频率；
(3c)将5个时刻当前跟踪目标的检验统计量，作为2×5的特征矩阵的第1行元素；
(3d)将5个时刻当前跟踪目标的多普勒频率，作为2×5的特征矩阵的第2行元素；
(4)从雷达接收机在当前跟踪时刻接收到的所有量测点中，找出符合关联条件的所有量测点；
(5)判断符合关联条件的量测点的数量是否等于0，若是，则将当前跟踪目标在当前跟踪时刻的关联状态设定为“无关联点”后，执行步骤(6)，否则，将当前跟踪目标在当前跟踪时刻的关联状态设定为“有关联点”后，执行步骤(7)；
(6)确定“无关联点”状态下的运动状态估计值和特征：
(6a)将当前跟踪目标在当前跟踪时刻运动状态的预测值，作为当前跟踪目标在当前跟踪时刻的运动状态估计值；
(6b)利用特征提取方法，从特征矩阵中得到当前跟踪目标在当前跟踪时刻的检验统计量、多普勒频率；
(6c)将得到的当前跟踪目标在当前跟踪时刻的运动状态估计值、检验统计量和多普勒频率加入历史数据后，执行步骤(9)；
(7)确定“有关联点”状态下与当前跟踪目标特征最匹配的关联点：
(7a)从所有关联点中，任选一个未被选过的关联点；
(7b)利用匹配值计算方法，计算所选关联点特征与目标特征矩阵的匹配值；
(7c)判定是否选完所有关联点：若是，则执行步骤(7d)；否则，执行步骤(7a)；
(7d)找出所有匹配值中的最小值，将该最小值所对应的关联点，作为与当前跟踪目标最匹配的量测点后执行步骤(8)；
(8)确定“有关联点”状态下的运动状态估计值和特征：
(8a)利用卡尔曼滤波公式，计算当前跟踪目标在当前跟踪时刻运动状态的估计值；
(8b)将最匹配量测点数据中包含的检验统计量、多普勒频率信息，分别作为当前跟踪目标在当前跟踪时刻的检验统计量、多普勒频率；
(8c)将得到的当前跟踪目标在当前跟踪时刻的运动状态估计值、检验统计量和多普勒频率加入历史数据后，执行步骤(9)；
(9)判断是否已选完当前跟踪时刻存在的所有目标；若是，执行步骤(10)；否则，执行步骤(1)；
(10)完成对当前跟踪时刻存在的所有目标在当前跟踪时刻的跟踪。
本发明与现有的技术相比具有以下优点：
第一，由于本发明对每一个目标建立了包含多个特征的特征矩阵，并在有关联点的状态下，将所有关联点的多个特征与特征矩阵进行匹配，克服了现有技术将属于其它目标的量测点或虚警点与当前跟踪目标关联，导致对当前跟踪目标运动状态的估计和未来轨迹的跟踪出现偏差的问题，使得本发明提高了量测点与待跟踪目标的关联准确率，进而能够对多个目标的轨迹实现更精确的跟踪。
第二，由于本发明通过计算所有关联点与特征矩阵匹配程度的方式，选择与目标最匹配的关联点，克服了现有技术计算方法较为复杂、计算量较大的问题，使得本发明在实际工程中实现更加容易，提高了工程效率。
附图说明
图1为本发明的流程图；
图2为本发明的仿真图。
具体实施方式
下面结合附图对本发明做进一步的描述。
结合附图1中，对本发明的具体步骤进一步的描述。
步骤1，从当前跟踪时刻存在的多个非机动目标中，选择一个未跟踪的目标，作为当前跟踪目标。
步骤2，预测运动状态。
将当前跟踪目标在当前跟踪时刻上一时刻的运动状态的估计值，代入状态方程，得到当前跟踪目标在当前跟踪时刻运动状态的预测值。
所述的运动状态，包括目标所处位置的直角坐标、目标在x轴方向的速度、目标在y轴方向的速度。
所述的状态方程如下：
其中，表示根据当前跟踪目标在第k-1时刻运动状态的估计值得到的当前跟踪目标在第k时刻运动状态的预测值，k的取值范围为大于等于1、小于等于雷达接收机工作时间的整数，Ak表示状态方程在k时刻的状态转移矩阵，表示当前跟踪目标在第k-1时刻的运动状态的估计值，Wk-1表示第k-1时刻环境中产生的动态噪声。
步骤3，建立当前跟踪目标的特征矩阵。
寻找与当前跟踪时刻最近的5个时刻。
从当前跟踪目标的历史数据中，提取当前跟踪目标在这5个时刻的检验统计量和多普勒频率。
所述的当前跟踪目标的历史数据，包含了当前跟踪目标从开始运动到当前跟踪时刻上一时刻的所有运动状态估计值、检验统计量和多普勒频率。
将5个时刻当前跟踪目标的检验统计量，作为2×5的特征矩阵的第1行元素。
将5个时刻当前跟踪目标的多普勒频率，作为2×5的特征矩阵的第2行元素。
步骤4，从雷达接收机在当前跟踪时刻接收到的所有量测点中，找出符合关联条件的所有量测点。
所述的关联条件为：
其中，zk表示待判断是否符合关联条件量测点的直角坐标，表示根据当前跟踪目标在k-1时刻运动状态的估计值得到的目标在k时刻量测点的预测直角坐标，Sk表示当前跟踪目标在k时刻的新息协方差矩阵，-1表示矩阵求逆操作，T表示矩阵转置操作，g表示跟踪波门的大小。
所述的新息协方差矩阵的计算步骤如下：
第1步，按照下式，计算当前跟踪目标在k时刻运动状态的协方差矩阵预测值：
Pk|k-1＝AkPk-1Ak T+Qk-1
其中，Pk|k-1表示根据当前跟踪目标在k-1时刻运动状态的协方差矩阵估计值，得到的当前跟踪目标在k时刻运动状态的协方差矩阵预测值，Ak表示状态方程在k时刻的状态转移矩阵，Pk-1表示当前跟踪目标在k-1时刻运动状态的协方差矩阵估计值，Qk-1表示k-1时刻环境中产生的动态噪声协方差矩阵。
第2步，按照下式，计算当前跟踪目标在k时刻的新息协方差矩阵：
Sk＝CkPk|k-1Ck T+Rk
其中，Sk表示当前跟踪目标在k时刻的新息协方差矩阵，Ck表示量测方程在k时刻的状态转移矩阵，Pk|k-1表示根据当前跟踪目标在k-1时刻运动状态的协方差矩阵估计值，得到的当前跟踪目标在k时刻运动状态的协方差矩阵预测值，T表示矩阵转置操作，Rk表示k时刻量测点的量测误差矩阵。
步骤5，判断符合关联条件的量测点的数量是否等于0，若是，则将当前跟踪目标在当前跟踪时刻的关联状态设定为“无关联点”后执行步骤6，否则，将当前跟踪目标在当前跟踪时刻的关联状态设定为“有关联点”后，执行步骤7。
步骤6，确定“无关联点”状态下的运动状态估计值和特征。
将当前跟踪目标在当前跟踪时刻运动状态的预测值，作为当前跟踪目标在当前跟踪时刻的运动状态估计值。
利用特征提取方法，从特征矩阵中得到当前跟踪目标在当前跟踪时刻的检验统计量、多普勒频率。
将得到的当前跟踪目标在当前跟踪时刻的运动状态估计值、检验统计量和多普勒频率加入历史数据后，执行步骤9。
所述的特征提取方法的具体步骤如下：
第1步，提取当前跟踪目标特征矩阵的第1行元素，作为目标历史检验统计量矩阵。
第2步，提取当前跟踪目标特征矩阵的第2行元素，作为目标历史多普勒频率矩阵。
第3步，将目标历史检验统计量矩阵的期望值、目标历史多普勒频率矩阵的期望值，分别作为当前跟踪目标在当前跟踪时刻的检验统计量、多普勒频率。
步骤7，确定“有关联点”状态下与当前跟踪目标特征最匹配的关联点。
(7.1)从所有关联点中，任选一个未被选过的关联点。
(7.2)利用匹配值计算方法，计算所选关联点特征与目标特征矩阵的匹配值。
所述的匹配值计算方法的具体步骤如下：
第1步，提取当前跟踪目标特征矩阵的第1行元素，作为目标历史检验统计量矩阵。
第2步，提取当前跟踪目标特征矩阵的第2行元素，作为目标历史多普勒频率矩阵。
第3步，按照下式，分别计算所选关联点与目标历史检验统计量矩阵元素的相似度、所选关联点与目标历史多普勒频率矩阵元素的相似度：
其中，l表示所选关联点与目标历史检验统计量矩阵元素或者与目标历史检验统计量矩阵元素的相似度，||表示求绝对值操作，c表示所选关联点的检验统计量或者多普勒频率，E表示目标历史检验统计量矩阵或者目标历史多普勒频率矩阵的期望值，σ表示目标历史检验统计量矩阵或者目标历史多普勒频率矩阵的标准差值。
第4步，计算每个关联点到当前跟踪目标在当前跟踪时刻预测位置的距离，从所有的距离值中选出最大值，作为归一化标准值。
第5步，按照下式，计算所选关联点到当前跟踪目标在当前跟踪时刻预测位置的归一化距离：
其中，表示所选关联点到当前跟踪目标在当前跟踪时刻预测位置的归一化距离，d表示所选关联点到当前跟踪目标在当前跟踪时刻预测位置的距离，D表示归一化标准值。
第6步，按照下式，计算所选量测点与目标特征矩阵的匹配值：
其中，M表示所选量测点与目标历史特征的匹配值，pstd表示特征匹配时检验统计量的参考权值，lstd表示所选量测点与目标历史检验统计量矩阵元素的相似度，pdpl表示特征匹配时多普勒频率的参考权值，ldpl表示所选量测点与目标历史多普勒频率矩阵元素的相似度，pd表示特征匹配时归一化距离的参考权值，表示所选量测点到当前跟踪目标在当前跟踪时刻运动状态预测值包含的直角坐标点的归一化距离。
(7.3)判定是否选完所有关联点：若是，则执行本步骤的(7.4)；否则，执行本步骤的(7.1)。
(7.4)找出所有匹配值中的最小值，将该最小值所对应的关联点，作为与当前跟踪目标最匹配的量测点后执行步骤8。
步骤8，确定“有关联点”状态下的运动状态估计值和特征。
利用卡尔曼滤波公式，计算当前跟踪目标在当前跟踪时刻运动状态的估计值。
所述的卡尔曼滤波公式如下：
其中，表示目标在k时刻的运动状态，Ak表示状态方程在k时刻的状态转移矩阵，表示目标在k-1时刻的运动状态，Hk表示k时刻的卡尔曼滤波增益矩阵，yk表示k时刻与目标关联的量测点直角坐标，Ck表示量测方程在k时刻的状态转移矩阵。
所述的卡尔曼滤波增益矩阵的计算公式如下：
Hk＝Pk|k-1Ck T(CkPk|k-1CT+Rk)-1
其中，Hk表示k时刻的卡尔曼滤波增益矩阵，Pk|k-1表示根据当前跟踪目标在k-1时刻运动状态的协方差矩阵估计值，得到的当前跟踪目标在k时刻运动状态的协方差矩阵预测值，Ck表示量测方程在k时刻的状态转移矩阵，T表示矩阵转置操作，Rk表示k时刻量测点的量测误差矩阵，-1表示矩阵求逆操作。
按照下式，计算当前跟踪目标在k时刻运动状态的协方差矩阵估计值，给k+1时刻跟踪使用：
Pk＝Pk|k-1-HkSkHk T
其中，Pk表示当前跟踪目标在k时刻运动状态的协方差矩阵估计值，Pk|k-1表示根据当前跟踪目标在k-1时刻运动状态的协方差矩阵估计值，得到的当前跟踪目标在k时刻运动状态的协方差矩阵预测值，Hk表示k时刻的卡尔曼滤波增益矩阵，Sk表示当前跟踪目标在k时刻的新息协方差矩阵，T表示矩阵转置操作。
将最匹配量测点数据中包含的检验统计量、多普勒频率信息，分别作为当前跟踪目标在当前跟踪时刻的检验统计量、多普勒频率。
将得到的当前跟踪目标在当前跟踪时刻的运动状态估计值、检验统计量和多普勒频率加入历史数据后，执行步骤9。
步骤9，判断是否已选完当前跟踪时刻存在的所有目标；若是，将当前跟踪时刻的下一时刻设为当前跟踪时刻，并执行步骤10；否则，执行步骤1。
步骤10，完成对当前跟踪时刻存在的所有目标在当前跟踪时刻的跟踪。
下面结合仿真实验对本发明的效果做进一步说明。
1.仿真实验条件：
本发明仿真实验的硬件测试平台是：CPU为intel Core i7 8700，主频为3.2GHz，内存16GB，软件平台是：Windows10专业版，64位操作系统，matlab r2018b。
2.仿真内容与仿真结果分析：
本发明的仿真实验是分别采用现有技术中的最近邻关联方法和本发明方法，对两个运动轨迹交叉的目标进行跟踪。
两个目标在第1秒时的初始状态：目标1的初始位置直角坐标为(0,6)，初始x轴方向速度为1，初始y轴方向速度为-1，检验统计量为10，并加上分布范围为[-1,1]的均匀噪声；目标2的初始位置直角坐标为(0,3.5)，初始x轴方向速度为1，初始y轴方向速度为-0.8，检验统计量为15，并加上分布范围为[-1,1]的均匀噪声。两个目标均作30秒的匀速直线运动，并将这两个目标在第2到第30秒位置的直角坐标转换为极坐标、x和y方向的速度分别转换为多普勒速度和切向速度。
状态噪声为方差等于0.0012的高斯噪声，径向距离和方位角的量测噪声均为方差等于0.052的高斯噪声。
第2秒到第30秒的虚警点位置分布服从泊松分布，虚警点x轴坐标取值范围为(-30,30)，y轴坐标分布范围为(-30,30)，设定每一秒在该范围内出现的个数期望值为100个；每一个虚警点的检验统计量为5，并加上分布范围为[0,2]的均匀噪声；每一个虚警点的多普勒速度为1，并加上分布范围为[-0.5,0.5]的均匀噪声。
本发明的仿真实验的结果如图2所示。其中，图2(a)为利用现有技术的最近邻关联方法对目标1、目标2方位角的跟踪结果图，图2(b)为利用本发明方法对目标1、目标2方位角的跟踪结果图；图2(c)为利用现有技术的最近邻关联方法对目标1、目标2径向距离的跟踪结果图，图2(d)为利用本发明方法对目标1、目标2径向距离的跟踪结果图。其中，图2(a)、图2(b)中的横坐标表示运动时刻，单位为秒，纵坐标表示方位角弧度值。图2(c)、图2(d)中的横坐标表示运动时刻，单位为秒，纵坐标表示径向距离，单位为千米。图2(a)、图2(b)、图2(c)、图2(d)中以“---”标示的曲线表示目标1在每个运动时刻方位角弧度真实值的变化曲线，以圆圈标示的曲线表示目标1在每个运动时刻方位角弧度量测值的变化曲线，以五角星标示的曲线表示目标1在每个运动时刻方位角弧度跟踪估计值的变化曲线，以“……”标示的曲线表示目标2在每个运动时刻径向距离真实值的变化曲线，以正方形标示的曲线表示目标2在每个运动时刻径向距离量测值的变化曲线，以星号标示的曲线表示目标2在每个运动时刻径向距离跟踪估计值的变化曲线。
从图2(a)、图2(c)可以看出，当两个目标的轨迹发生交叉时，现有技术的最近邻关联方法不仅不能将两个目标的轨迹分离，反而会使两个目标的跟踪轨迹发生重合。
而从图2(b)、图2(d)可以看出，当两目标轨迹发生交叉时，本发明方法依然能有较好的跟踪效果，正确分离两条重合的轨迹。
综上所述，本发明的仿真实验证明了在目标轨迹交叉的情况下，本发明能够对目标运动轨迹进行有效、准确的跟踪。

Claims (8)

1.一种基于多特征匹配的多目标跟踪方法，其特征在于，建立当前跟踪目标的特征矩阵，确定“无关联点”状态下的运动状态估计值和特征，确定“有关联点”状态下与当前跟踪目标特征最匹配的关联点，确定“有关联点”状态下的运动状态估计值和特征；该方法的具体步骤包括如下：
(1)从当前跟踪时刻存在的多个目标中，选择一个未跟踪的目标，作为当前跟踪目标；
(2)预测运动状态：
将当前跟踪目标在当前跟踪时刻上一时刻的运动状态的估计值，代入状态方程，得到当前跟踪目标在当前跟踪时刻运动状态的预测值；
(3)建立当前跟踪目标的特征矩阵：
(3a)寻找与当前跟踪时刻最近的5个时刻；
(3b)从当前跟踪目标的历史数据中，提取当前跟踪目标在这5个时刻的检验统计量和多普勒频率；
(3c)将5个时刻当前跟踪目标的检验统计量，作为2×5的特征矩阵的第1行元素；
(3d)将5个时刻当前跟踪目标的多普勒频率，作为2×5的特征矩阵的第2行元素；
(4)从雷达接收机在当前跟踪时刻接收到的所有量测点中，找出符合关联条件的所有量测点；
(5)判断符合关联条件的量测点的数量是否等于0，若是，则将当前跟踪目标在当前跟踪时刻的关联状态设定为“无关联点”后，执行步骤(6)，否则，将当前跟踪目标在当前跟踪时刻的关联状态设定为“有关联点”后，执行步骤(7)；
(6)确定“无关联点”状态下的运动状态估计值和特征：
(6a)将当前跟踪目标在当前跟踪时刻运动状态的预测值，作为当前跟踪目标在当前跟踪时刻的运动状态估计值；
(6b)利用特征提取方法，从特征矩阵中得到当前跟踪目标在当前跟踪时刻的检验统计量、多普勒频率；
(6c)将得到的当前跟踪目标在当前跟踪时刻的运动状态估计值、检验统计量和多普勒频率加入历史数据后，执行步骤(9)；
(7)确定“有关联点”状态下与当前跟踪目标特征最匹配的关联点：
(7a)从所有关联点中，任选一个未被选过的关联点；
(7b)利用匹配值计算方法，计算所选关联点特征与目标特征矩阵的匹配值；
(7c)判定是否选完所有关联点：若是，则执行步骤(7d)；否则，执行步骤(7a)；
(7d)找出所有匹配值中的最小值，将该最小值所对应的关联点，作为与当前跟踪目标最匹配的量测点后执行步骤(8)；
(8)确定“有关联点”状态下的运动状态估计值和特征：
(8a)利用卡尔曼滤波公式，计算当前跟踪目标在当前跟踪时刻运动状态的估计值；
(8b)将最匹配量测点数据中包含的检验统计量、多普勒频率信息，分别作为当前跟踪目标在当前跟踪时刻的检验统计量、多普勒频率；
(8c)将得到的当前跟踪目标在当前跟踪时刻的运动状态估计值、检验统计量和多普勒频率加入历史数据后，执行步骤(9)；
(9)判断是否已选完当前跟踪时刻存在的所有目标；若是，执行步骤(10)；否则，执行步骤(1)；
(10)完成对当前跟踪时刻存在的所有目标在当前跟踪时刻的跟踪。
2.根据权利要求1所述的基于多特征匹配的多目标跟踪方法，其特征在于：步骤(2)中所述的运动状态，包括目标所处位置的直角坐标、目标在x轴方向的速度、目标在y轴方向的速度。
3.根据权利要求1所述的基于多特征匹配的多目标跟踪方法，其特征在于：步骤(2)中所述的状态方程如下：
其中，表示根据当前跟踪目标在第k-1时刻运动状态的估计值得到的当前跟踪目标在第k时刻运动状态的预测值，k的取值范围为大于等于1、小于等于雷达接收机工作时间的整数，Ak表示状态方程在k时刻的状态转移矩阵，表示当前跟踪目标在第k-1时刻的运动状态的估计值，Wk-1表示第k-1时刻环境中产生的动态噪声。
4.根据权利要求1所述的基于多特征匹配的多目标跟踪方法，其特征在于：步骤(3b)中所述的当前跟踪目标的历史数据包含当前跟踪目标从开始运动到当前跟踪时刻上一时刻的所有运动状态估计值、检验统计量和多普勒频率。
5.根据权利要求1所述的基于多特征匹配的多目标跟踪方法，其特征在于：步骤(4)中所述的关联条件为：
其中，zk表示待判断是否符合关联条件量测点的直角坐标，表示根据当前跟踪目标在k-1时刻运动状态的估计值得到的目标在k时刻位置的预测直角坐标，Sk表示当前跟踪目标在k时刻的新息协方差矩阵，-1表示矩阵求逆操作，T表示矩阵转置操作，g表示跟踪波门的大小。
6.根据权利要求1所述的基于多特征匹配的多目标跟踪方法，其特征在于：步骤(6b)中所述的特征提取方法的具体步骤如下：
第一步，提取当前跟踪目标特征矩阵的第1行元素，作为目标历史检验统计量矩阵；
第二步，提取当前跟踪目标特征矩阵的第2行元素，作为目标历史多普勒频率矩阵；
第三步，将目标历史检验统计量矩阵的期望值、目标历史多普勒频率矩阵的期望值，分别作为当前跟踪目标在当前跟踪时刻的检验统计量、多普勒频率。
7.根据权利要求1所述的基于多特征匹配的多目标跟踪方法，其特征在于：步骤(7b)中所述的匹配值计算方法的具体步骤如下：
第一步，提取当前跟踪目标特征矩阵的第1行元素，作为目标历史检验统计量矩阵；
第二步，提取当前跟踪目标特征矩阵的第2行元素，作为目标历史多普勒频率矩阵；
第三步，按照下式，分别计算所选关联点与目标历史检验统计量矩阵元素的相似度、所选关联点与目标历史多普勒频率矩阵元素的相似度：
其中，l表示所选关联点与目标历史检验统计量矩阵元素或者与目标历史检验统计量矩阵元素的相似度，||表示求绝对值操作，c表示所选关联点的检验统计量或者多普勒频率，E表示目标历史检验统计量矩阵或者目标历史多普勒频率矩阵的期望值，σ表示目标历史检验统计量矩阵或者目标历史多普勒频率矩阵的标准差值；
第四步，计算每个关联点到当前跟踪目标在当前跟踪时刻预测位置的距离，从所有的距离值中选出最大值，作为归一化标准值；
第五步，按照下式，计算所选关联点到当前跟踪目标在当前跟踪时刻预测位置的归一化距离：
其中，表示所选关联点到当前跟踪目标在当前跟踪时刻预测位置的归一化距离，d表示所选关联点到当前跟踪目标在当前跟踪时刻预测位置的距离，D表示归一化标准值；
第六步，按照下式，计算所选量测点与目标特征矩阵的匹配值：
其中，M表示所选量测点与目标历史特征的匹配值，pstd表示特征匹配时检验统计量的参考权值，lstd表示所选量测点与目标历史检验统计量矩阵元素的相似度，pdpl表示特征匹配时多普勒频率的参考权值，ldpl表示所选量测点与目标历史多普勒频率矩阵元素的相似度，pd表示特征匹配时归一化距离的参考权值，表示所选量测点到当前跟踪目标在当前跟踪时刻运动状态预测值包含的直角坐标点的归一化距离。
8.根据权利要求1所述的基于多特征匹配的多目标跟踪方法，其特征在于：步骤(8a)所述的卡尔曼滤波公式如下：
其中，表示目标在k时刻的运动状态，Ak表示状态方程的状态转移矩阵，表示目标在k-1时刻的运动状态，Hk表示k时刻的卡尔曼滤波增益矩阵，yk表示k时刻与目标关联的量测点直角坐标，Ck表示量测方程的状态转移矩阵。
CN201811577747.8A 2018-12-20 2018-12-20 基于多特征匹配的多目标跟踪方法 Active CN109521420B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811577747.8A CN109521420B (zh) 2018-12-20 2018-12-20 基于多特征匹配的多目标跟踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811577747.8A CN109521420B (zh) 2018-12-20 2018-12-20 基于多特征匹配的多目标跟踪方法

Publications (2)

Publication Number Publication Date
CN109521420A true CN109521420A (zh) 2019-03-26
CN109521420B CN109521420B (zh) 2022-10-04

Family

ID=65797021

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811577747.8A Active CN109521420B (zh) 2018-12-20 2018-12-20 基于多特征匹配的多目标跟踪方法

Country Status (1)

Country Link
CN (1) CN109521420B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110780289A (zh) * 2019-10-23 2020-02-11 北京信息科技大学 基于场景雷达的多目标车辆跟踪方法及装置
CN111007880A (zh) * 2019-12-24 2020-04-14 桂林电子科技大学 一种基于汽车雷达的扩展目标跟踪方法
CN111581320A (zh) * 2020-05-09 2020-08-25 海南中智信信息技术有限公司 一种基于行为轨迹的碰头地点及时长识别方法及系统
CN113311448A (zh) * 2020-02-27 2021-08-27 郑州宇通客车股份有限公司 一种基于多特征信息的动态目标跟踪方法与装置
CN113902777A (zh) * 2020-06-22 2022-01-07 东南大学 多目标跟踪方法、系统以及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002168951A (ja) * 2000-12-04 2002-06-14 Mitsubishi Electric Corp 目標追尾方法とその装置
US20030048224A1 (en) * 2001-05-04 2003-03-13 Lockheed Martin Corporation System and method for measurement domain data association in passive coherent location applications
CN106443662A (zh) * 2016-10-28 2017-02-22 上海无线电设备研究所 一种低重频体制下速度模糊时的目标稳健关联方法
US20170116753A1 (en) * 2014-04-30 2017-04-27 Institute Of Automation Chinese Academy Of Sciences Large-Range-First Cross-Camera Visual Target Re-identification Method
CN106872955A (zh) * 2017-01-24 2017-06-20 西安电子科技大学 基于联合概率数据关联算法的雷达多目标跟踪优化方法
CN108198206A (zh) * 2017-12-22 2018-06-22 湖南源信光电科技股份有限公司 基于多特征联合和Camshift算法结合的多目标跟踪方法
CN108921873A (zh) * 2018-05-29 2018-11-30 福州大学 基于核相关滤波优化的马尔科夫决策在线多目标跟踪方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002168951A (ja) * 2000-12-04 2002-06-14 Mitsubishi Electric Corp 目標追尾方法とその装置
US20030048224A1 (en) * 2001-05-04 2003-03-13 Lockheed Martin Corporation System and method for measurement domain data association in passive coherent location applications
US20170116753A1 (en) * 2014-04-30 2017-04-27 Institute Of Automation Chinese Academy Of Sciences Large-Range-First Cross-Camera Visual Target Re-identification Method
CN106443662A (zh) * 2016-10-28 2017-02-22 上海无线电设备研究所 一种低重频体制下速度模糊时的目标稳健关联方法
CN106872955A (zh) * 2017-01-24 2017-06-20 西安电子科技大学 基于联合概率数据关联算法的雷达多目标跟踪优化方法
CN108198206A (zh) * 2017-12-22 2018-06-22 湖南源信光电科技股份有限公司 基于多特征联合和Camshift算法结合的多目标跟踪方法
CN108921873A (zh) * 2018-05-29 2018-11-30 福州大学 基于核相关滤波优化的马尔科夫决策在线多目标跟踪方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CHANGYUN LIU ET AL.: "Modified unscented Kalman filter using modified filter gain and variance scale factor for highly maneuvering target tracking", 《JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS》 *
MINGGANG GAN ET AL.: "Hierarchical particle filter tracking algorithm based on multi-feature fusion", 《JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS》 *
YU MA ET AL.: "Multi-target tracking based on multi-feature associated matrix reasoning", 《2011 4TH INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING》 *
倪琦等: "基于多特征融合的Mean_Shift跟踪算法", 《强激光与粒子束》 *
高岚等: "基于DSmT的机载目标断续航迹关联", 《传感器与微系统》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110780289A (zh) * 2019-10-23 2020-02-11 北京信息科技大学 基于场景雷达的多目标车辆跟踪方法及装置
CN111007880A (zh) * 2019-12-24 2020-04-14 桂林电子科技大学 一种基于汽车雷达的扩展目标跟踪方法
CN111007880B (zh) * 2019-12-24 2022-09-02 桂林电子科技大学 一种基于汽车雷达的扩展目标跟踪方法
CN113311448A (zh) * 2020-02-27 2021-08-27 郑州宇通客车股份有限公司 一种基于多特征信息的动态目标跟踪方法与装置
CN111581320A (zh) * 2020-05-09 2020-08-25 海南中智信信息技术有限公司 一种基于行为轨迹的碰头地点及时长识别方法及系统
CN113902777A (zh) * 2020-06-22 2022-01-07 东南大学 多目标跟踪方法、系统以及存储介质
CN113902777B (zh) * 2020-06-22 2024-07-02 东南大学 多目标跟踪方法、系统以及存储介质

Also Published As

Publication number Publication date
CN109521420B (zh) 2022-10-04

Similar Documents

Publication Publication Date Title
CN109521420A (zh) 基于多特征匹配的多目标跟踪方法
CN103729859B (zh) 一种基于模糊聚类的概率最近邻域多目标跟踪方法
CN104297748B (zh) 一种基于轨迹增强的雷达目标检测前跟踪方法
CN103885057B (zh) 自适应变滑窗多目标跟踪方法
CN105510896B (zh) 一种集中式多雷达数据处理的加权最近邻域数据关联方法
CN107396322A (zh) 基于路径匹配与编码译码循环神经网络的室内定位方法
CN107300698B (zh) 一种基于支持向量机的雷达目标航迹起始方法
CN104035083B (zh) 一种基于量测转换的雷达目标跟踪方法
Niedfeldt et al. Recursive RANSAC: Multiple signal estimation with outliers
CN108734725B (zh) 基于高斯过程的概率数据关联滤波扩展目标跟踪方法
CN108664930A (zh) 一种智能的多目标检测跟踪方法
CN104182609B (zh) 基于去相关的无偏转换量测的三维目标跟踪方法
CN108168564A (zh) 一种基于lhd灰色关联度的航迹关联方法
Ma et al. Target tracking system for multi-sensor data fusion
CN107339992A (zh) 一种基于行为的室内定位和地标语义标识的方法
CN104050368B (zh) 系统误差下基于误差补偿的群航迹精细关联算法
CN109324315A (zh) 基于双层次块稀疏性的空时自适应处理雷达杂波抑制方法
CN110187335A (zh) 针对具有非连续特性目标的粒子滤波检测前跟踪方法
CN105654518A (zh) 一种基于方差估计的跟踪模板自适应方法
Lin et al. Mobile robot self-localizationusing visual odometry based on ceiling vision
CN106291530B (zh) 一种基于最近邻法的概率数据关联优化方法
CN104569963A (zh) 用于超宽带穿墙雷达运动目标一维检测与跟踪的方法
CN110334674A (zh) 一种平面自由体轨迹识别跟踪及预测方法
CN114994656A (zh) 一种基于毫米波雷达的室内人员跟踪方法
CN111274529B (zh) 一种鲁棒的高斯逆威沙特phd多扩展目标跟踪算法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant