CN104268865B - 基于模板匹配的集中式多传感器编队目标跟踪方法 - Google Patents
基于模板匹配的集中式多传感器编队目标跟踪方法 Download PDFInfo
- Publication number
- CN104268865B CN104268865B CN201410480203.5A CN201410480203A CN104268865B CN 104268865 B CN104268865 B CN 104268865B CN 201410480203 A CN201410480203 A CN 201410480203A CN 104268865 B CN104268865 B CN 104268865B
- Authority
- CN
- China
- Prior art keywords
- formation
- measurement
- matrix
- matching
- template
- 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.)
- Active
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 52
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 238000001514 detection method Methods 0.000 claims abstract description 9
- 238000005755 formation reaction Methods 0.000 claims description 90
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 238000012790 confirmation Methods 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种可以在复杂背景下对非机动编队目标进行精确跟踪的集中式数据互联跟踪方法。该方法基于非机动情况下各探测周期内编队内目标真实回波位置相对固定的特性，通过预关联成功的编队状态集合与编队量测集合分别建立模板形状矩阵和待匹配形状矩阵，利用匹配搜索模型和匹配矩阵确认规则选出代价最小的匹配矩阵，并基于模板和对应的匹配矩阵利用kalman滤波完成编队内各目标航迹的状态更新，该方法能最大限度的消除杂波，并保证编队内目标结构的稳定性，从而不会出现编队内目标交叉错误互联的情况；同时，该方法能充分利用多个传感器的探测信息，避免冗余航迹的产生。
Description
技术领域
本发明涉及数据处理领域中基于模板匹配的集中式多传感器编队目标跟踪方法，适用于各种集中式结构的信息融合系统。
背景技术
在现实环境中，经常因为不可控制或特定人为目的等因素，会在一个较小的空域范围内构成一个复杂的目标群，如空间碎片的分裂、弹道导弹突防过程中伴随的大量诱饵、导弹和飞机编队等，这些目标分布范围较小，运动特征差异不明显，相对运动速度较低且特性接近。目标跟踪领域将此类目标称为编队目标。传统的多传感器多目标跟踪算法对编队目标的跟踪效果十分有限。此类算法通常基于测量直接对编队内目标进行建航，但因编队内目标间距较小，各目标跟踪波门会严重交叠，数据互联难度大增；而且，因编队内目标行为模式相似，错误的航迹起始及维持能在后续时刻得以延续，易造成整体态势混乱。
为解决该问题，近年来国内外学者提出了一系列编队目标跟踪算法，基本思路大多为：利用各种技术建立编队的等效量测，基于等效量测实现目标编队的整体跟踪。其优点是避免了编队内目标的相互影响，降低了跟踪混乱和计算量爆炸的几率，提高了整个跟踪系统的稳定性，节省了大量的雷达资源。但随着传感器分辨率的提高，逐步表现出以下不足：第一，现有编队目标跟踪算法的推导环境大多比较简单，通常假设编队中个体目标完全可辨，然而在实际探测过程中，因目标的互相遮挡、传感器分辨率不够充分等因素，编队目标通常是部分可辨的；第二，在一些实际工程应用中，如低空编队突防目标的拦截、编队内具有特殊价值的目标跟踪等，在跟踪整个编队的同时，十分需要对编队内个体目标进行单独跟踪，然而现有编队目标跟踪算法通常只能得到编队整体的状态，没有考虑编队内目标的精确跟踪问题。特别地，当工程上利用组网传感器探测编队目标时，必然需要进行数据互联和融合等处理，尤其是在云雨杂波、带状干扰等复杂背景下对于集中式多传感器编队目标跟踪技术，目前尚没有文献进行研究。
针对上述问题，有必要研究如何在复杂背景下实现编队内目标的多维点-航互联及量测合并，完成复杂环境下集中式多传感器非机动编队内目标的精确跟踪。这也是本发明的思路来源。
发明内容
要解决的技术问题
本发明的目的在于提供一种可以在复杂背景下对非机动编队目标进行精确跟踪的集中式数据互联跟踪方法。为了达到上述目的，本发明基于非机动情况下各探测周期内编队内目标真实回波位置相对固定的特性，利用模板匹配的方式区别目标真实测量和杂波，提出了一种基于模板匹配的集中式多传感器编队目标跟踪方法。该方法最大限度的消除了杂波，并保证了编队内目标结构的稳定性，从而不会出现编队内目标交叉错误互联的情况；同时，该方法既充分利用了多个传感器的探测信息，又避免了冗余航迹的产生。
技术方案
本发明所述的基于模板匹配的集中式多传感器编队目标跟踪方法，包括以下技术措施：首先，进行编队整体的预互联，利用循环阈值法进行编队的预分割，将量测集分割为不同的编队量测，再将编队量测与编队航迹的互联转变为编队中心量测与对应编队中心航迹的互联，从而实现编队的预互联；然后，建立模板匹配模型，通过预关联成功的编队状态集合与编队量测集合分别建立模板形状矩阵和待匹配形状矩阵，利用匹配搜索模型和匹配矩阵确认规则选出代价最小的匹配矩阵；最后，实现编队内航迹的更新，基于模板和对应的匹配矩阵利用kalman滤波，从而完成编队内各目标航迹的状态更新。
有益效果
本发明利用非机动模式下各时刻同一编队内目标真实回波整体结构的相似性，基于模板匹配消除背景杂波，实现编队目标的精确跟踪。其主要有以下优点：
(1)该跟踪方法通过模板与待匹配区域的搜索匹配，最大限度的消除了杂波，并保证了编队内目标结构的稳定性，不会出现编队内目标交叉错误互联的情况；
(2)该跟踪方法基于匹配矩阵确认模型，在所有传感器的匹配矩阵中，获取代价最小的匹配矩形完成编队内各目标的状态更新，既充分利用了多个传感器的探测信息，又避免了冗余航迹的产生。
附图说明
图1为跟踪方法流程图；
图2为循环阈值模型流程图；
图3为矩阵A划分示意图；
图4为匹配搜索模型示意图；
具体实施方式
以下结合说明书附图，对技术方案中的具体步骤进行详细描述。
1.基于编队整体的预互联
设G(k-1)为k-1时刻编队航迹状态更新值和协方差更新值的集合，
式中，Tg(k-1)为k-1时刻编队的个数；为k-1时刻第t个编队中的目标个数。
针对第i个传感器上报的量测集合Zi(k)，利用附图2中的循环阈值法进行编队的预分割，其中d0为常数阈值。设为分割后获得的第l个编队量测， 为中的量测个数。
设为的中心点，为编队Gt(k-1)中心航迹的状态更新值和协方差更新值。此时，编队量测与编队航迹的互联，暂时转变为与之间的互联。与传统目标相同，以的一步预测值为中心，建立关联波门。若满足式(2)，则认为落入关联波门内。
式中，γ为常数阈值；为新息协方差；
通常情况下，编队与编队之间相距较远，落入同一关联门的概率较小；假如同一传感器有多个编队量测落入同一关联波门内，基于最近邻域的思想，选取最近的编队量测为关联量。设各传感器上报量测中与Gt(k-1)关联成功的编队量测集合为
式中，Ns′为存在关联编队量测的传感器个数。
2.模板匹配模型的建立
基于预关联成功的Gt(k-1)和，建立模板匹配模型，剔除杂波并实现编队内航迹与真实量测的互联，具体过程分为以下四步。
1)模板形状矩阵的建立
基于Gt(k-1)建立模板，并获取x、y方向量测最大值和最小值的集合基于Exy(k)构建附图3所示的矩形A。
可以描述为
式中，为服从高斯分布的随机噪声；ωi为k时刻传感器i的位置；ωt为目标t 的真实位置；在传感器探测区域内服从均匀分布。当编队未发生机动时，编队内各目标真实位置构成的图形在相临时刻是基本一致的；受量测误差的影响，量测值会发生一定的摆动，由式(4)可知，在x、y方向摆动的范围大致为 按附图3的方式将矩形A划分成多个小矩形，其中μx、μy分别为x、y方向上的划分依据，其定义为
其中，分别为是Gt(k-1)各目标在x、y方向的最小差值。
定义矩形A模板形状矩阵为
式中，分别为矩形A在x、y方向上的矩形个数；
2)待匹配形状矩阵的建立
以为待匹配对象，建立待匹配形状矩阵描述具体分为以下三步：
(1)基于找出x、y方向上的最大值和最小值，建立如附图4所示的矩形B；
(2)依据把矩形B划分成多个小矩形，其中和的定义同式(5)。
(3)基于有无量测落入小矩形中建立待匹配形状矩阵为
其中，的定义同式(7)。同时记录小矩形内量测集为
其中，为量测集合中落入小矩形的量测个数。
3)匹配搜索模型的建立
如附图4所示，以矩形B在x、y方向上最小的矩形为起点，将矩形A嵌入B内，基于与建立Bij与A的匹配度矩阵为
其中，分别为x、y方向上进行匹配的次数，且
其中，
若则认为Bij与A匹配，即Bij中各量测的内部结构与A中各量测的内部结构相似。
4)匹配矩阵的确认
矩形B中满足匹配条件的Bij可能有多个，在此通过建立匹配代价矩阵确定最终的匹配矩阵。基于传感器s的匹配度矩阵定义匹配代价矩阵为
其中，di′j′为Bij中第i′j′个方格中量测与模板A对应方格中所包含目标状态一步的预测值的差值，即
同时，储存Bij各方格中用于计算di′j′的量测集对s个传感器，均建立取代价最小的矩阵为匹配矩形，记中各方格的互联量测集为
3.编队内航迹的状态更新
基于Gt(k-1)和利用模板A和匹配矩形若落入A中的小矩形ai′j′中，则基于落入中的量测利用kalman滤波完成状态和协方差更新。
在此，需要注意的是，在完成编队内所有航迹的状态更新后，还需要利用当前时刻剩余的量测集与前三个时刻剩余的量测集进行航迹起始，并建立航迹终结原则，判断编队内各条航迹是否终结。
Claims (3)
1.基于模板匹配的集中式多传感器编队目标跟踪方法，是一种用于密集杂波背景下编队目标跟踪的集中式跟踪方法，其特征在于包括以下步骤：
步骤1：利用集中式多传感器的量测信息，基于编队整体的预关联；
步骤2：基于预关联结果建立模板匹配模型，具体包括模板形状矩阵、待匹配形状矩阵、匹配度矩阵和匹配代价矩阵，剔除杂波并实现编队内航迹与真实量测的关联；
具体为，基于预关联成功的Gt(k-1)建立模板，生成模板形状矩阵，其中Gt(k-1)为k-1时刻第t个编队航迹状态更新值和协方差更新值的集合；以为待匹配对象，建立待匹配形状矩阵，其中为各传感器上报量测中与Gt(k-1)关联成功的编队量测集合；将模板形状矩阵和待匹配形状矩阵通过匹配搜索模型进行匹配，建立匹配度矩阵；由于匹配度矩阵具有不唯一性，在此建立匹配代价矩阵确定代价最小的匹配矩阵，从而得到关联量测集为由此，便建立了模板匹配模型，通过剔除杂波并实现编队内航迹与真实量测的关联；
步骤3：编队内航迹的状态更新。
2.根据权利要求1所述的方法，其特征在于步骤1包括：由于密集杂波环境下易产生交叉关联，造成多跟、错跟情况，要实现多传感器探测时编队内目标的精确跟踪，编队内目标与多传感量测值间的对应关联是关键，为此需要对编队目标进行预关联；对于第i个传感器上报的量测集，首先进行编队的预分割，将量测集分割为不同的编队量测；然后将编队量测与编队航迹的关联转变为与之间的关联，其中为采用循环阈值法对编队进行预分割后获得的第l个编队量测的中心点， 为k-1时刻第t个编队中心航迹的状态更新值，为k-1时刻第t个编队中心航迹的协方差更新值；与传统目标相同，以的一步预测值为中心，建立关联波门；若满足下式，则认为落入关联波门内；
其中，γ为常数阈值；为新息协方差；若同一传感器有多个编队量测落入同一关联波门内，基于最近邻域的思想，选取最近的编队量测为关联量；设各传感器上报量测中与Gt(k-1)关联成功的编队量测集合为其中Ns′为存在关联编队量测的传感器个数。
3.根据权利要求1所述的方法，其特征在于步骤3包括：基于Gt(k-1)和利用kalman滤波完成编队状态和协方差的更新，在完成编队内所有航迹的状态更新后， 还需要利用当前时刻剩余的量测集与前三个时刻剩余的量测集进行航迹起始，并建立航迹终结原则，判断编队内各条航迹是否终结。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410480203.5A CN104268865B (zh) | 2014-09-18 | 2014-09-18 | 基于模板匹配的集中式多传感器编队目标跟踪方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410480203.5A CN104268865B (zh) | 2014-09-18 | 2014-09-18 | 基于模板匹配的集中式多传感器编队目标跟踪方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104268865A CN104268865A (zh) | 2015-01-07 |
| CN104268865B true CN104268865B (zh) | 2017-06-16 |
Family
ID=52160384
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410480203.5A Active CN104268865B (zh) | 2014-09-18 | 2014-09-18 | 基于模板匹配的集中式多传感器编队目标跟踪方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104268865B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107290731B (zh) * | 2017-06-21 | 2019-12-10 | 哈尔滨工业大学 | 一种基于图像域模板匹配的雷达航迹起始方法 |
| CN109765928A (zh) * | 2019-01-07 | 2019-05-17 | 杭州电子科技大学 | 基于移动多智能体编队的协作式控制与目标跟踪方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102646279A (zh) * | 2012-02-29 | 2012-08-22 | 北京航空航天大学 | 一种基于运动预测与多子块模板匹配相结合的抗遮挡跟踪方法 |
| CN102737386A (zh) * | 2012-05-15 | 2012-10-17 | 北京硅盾安全技术有限公司 | 一种运动目标抗融合遮挡跟踪算法 |
| CN103677734A (zh) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | 基于特征匹配矩阵的多目标数据关联算法 |
| CN103886196A (zh) * | 2014-03-14 | 2014-06-25 | 中国电子科技集团公司第二十八研究所 | 一种抗大系统误差的航迹关联方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101984362B (zh) * | 2010-08-12 | 2012-10-31 | 中国人民解放军海军航空工程学院 | 基于数据压缩的集中式多源广义相关跟踪方法 |
| CN103471591B (zh) * | 2013-04-15 | 2017-06-06 | 中国人民解放军海军航空工程学院 | 基于逻辑法、全局最近邻和目标航向信息的机动多目标数据互联方法 |
| CN104050368B (zh) * | 2014-06-09 | 2017-04-12 | 中国人民解放军海军航空工程学院 | 系统误差下基于误差补偿的群航迹精细关联算法 |
-
2014
- 2014-09-18 CN CN201410480203.5A patent/CN104268865B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102646279A (zh) * | 2012-02-29 | 2012-08-22 | 北京航空航天大学 | 一种基于运动预测与多子块模板匹配相结合的抗遮挡跟踪方法 |
| CN102737386A (zh) * | 2012-05-15 | 2012-10-17 | 北京硅盾安全技术有限公司 | 一种运动目标抗融合遮挡跟踪算法 |
| CN103677734A (zh) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | 基于特征匹配矩阵的多目标数据关联算法 |
| CN103886196A (zh) * | 2014-03-14 | 2014-06-25 | 中国电子科技集团公司第二十八研究所 | 一种抗大系统误差的航迹关联方法 |
Non-Patent Citations (4)
| Title |
|---|
| 一种基于模板匹配的运动目标跟踪方法;党晓军 等;《计算机工程与应用》;20101231;第46卷(第5期);全文 * |
| 基于多关联模板匹配的火箭目标跟踪算法;尹宏鹏 等;《计算机工程》;20090930;第35卷(第18期);全文 * |
| 基于局部模板匹配的运动目标跟踪;卢建国 等;《光电子·激光》;20110228;第22卷(第2期);全文 * |
| 基于模板匹配的多模式车辆跟踪算法;林明秀 等;《系统仿真学报》;20070430;第19卷(第7期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104268865A (zh) | 2015-01-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103729859B (zh) | 一种基于模糊聚类的概率最近邻域多目标跟踪方法 | |
| CN106443622B (zh) | 一种基于改进联合概率数据关联的分布式目标跟踪方法 | |
| CN106842165B (zh) | 一种基于不同距离角度分辨率雷达集中式异步融合方法 | |
| CN104237879B (zh) | 一种雷达系统中的多目标跟踪方法 | |
| CN104050368B (zh) | 系统误差下基于误差补偿的群航迹精细关联算法 | |
| Ma et al. | Target tracking system for multi-sensor data fusion | |
| CN107229033B (zh) | 基于高度维分段搜索的多目标到达时间差定位方法 | |
| CN106291533A (zh) | 一种基于amd的分布式多传感器融合算法 | |
| CN103471591A (zh) | 基于逻辑法、全局最近邻和目标航向信息的机动多目标数据互联算法 | |
| CN107340516B (zh) | 基于多普勒速度的联合逻辑快速航迹起始方法 | |
| CN103152826A (zh) | 一种基于nlos状态检测补偿的移动目标追踪方法 | |
| CN106443664B (zh) | 系统误差下基于拓扑信息的雷达与esm航迹关联方法 | |
| CN109239702B (zh) | 一种基于目标状态集的机场低空飞鸟数量统计方法 | |
| CN105467370A (zh) | 一种复合双基地雷达进动目标isar图像横向定标方法 | |
| CN110187336A (zh) | 一种基于分布式phd的多站雷达站址定位和联合跟踪方法 | |
| CN104569923B (zh) | 基于速度约束的Hough变换快速航迹起始方法 | |
| CN103809161B (zh) | 雷达网抗距离欺骗+soj复合干扰方法 | |
| CN108445505A (zh) | 线环境下基于激光雷达的特征显著性检测方法 | |
| CN104535993A (zh) | 一种机载多主动雷达测距的地面物体高精度定位方法 | |
| CN104050641B (zh) | 基于形状方位描述符的集中式多传感器编队目标粒子滤波算法 | |
| CN103913740A (zh) | 一种基于空域分布特征的飞鸟群目标跟踪方法 | |
| CN104268865B (zh) | 基于模板匹配的集中式多传感器编队目标跟踪方法 | |
| CN106443665B (zh) | 基于辅助距离的雷达与esm航迹关联方法 | |
| CN107238835A (zh) | 一种编队目标点航抗叉关联方法 | |
| CN104020468B (zh) | 一种基于二次雷达的近距离目标更新周期的计算方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Haipeng Inventor after: Xu Congan Inventor after: Lin Xueyuan Inventor after: Wang Cong Inventor after: Dong Kai Inventor after: Liu Jun Inventor after: Pan Xinlong Inventor after: Guo Qiang Inventor before: Wang Haipeng Inventor before: Xu Congan Inventor before: Lin Xueyuan Inventor before: Wang Cong Inventor before: Dong Kai Inventor before: Liu Jun Inventor before: Pan Xinlong |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant |