CN113344039B - 一种基于时空关联的多扩展目标跟踪方法 - Google Patents
一种基于时空关联的多扩展目标跟踪方法 Download PDFInfo
- Publication number
- CN113344039B CN113344039B CN202110547305.4A CN202110547305A CN113344039B CN 113344039 B CN113344039 B CN 113344039B CN 202110547305 A CN202110547305 A CN 202110547305A CN 113344039 B CN113344039 B CN 113344039B
- Authority
- CN
- China
- Prior art keywords
- target
- matrix
- clustering
- gaussian
- extended
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000005259 measurement Methods 0.000 claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims description 34
- 238000001914 filtration Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 238000013138 pruning Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000000638 solvent extraction Methods 0.000 abstract description 3
- 101001086191 Borrelia burgdorferi Outer surface protein A Proteins 0.000 description 6
- 230000004083 survival effect Effects 0.000 description 5
- 238000005192 partition Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于时空关联的多扩展目标跟踪方法。当目标占据传感器的多个分变率单元时,单个目标便会产生多个量测值,即为扩展目标。在这种背景下,当扩展目标产生交叉,一般的基于距离的划分方法会将不同目标的量测值划入同一量测集合中,造成滤波器的精度下降,势估计出现错误。本发明基于ET‑GM‑PHD算法,采用时空关联思想,利用扩展目标在相邻时刻量测值的关联性,并在一种有向图SNN划分基础上,对多扩展目标进行跟踪。本发明方法大大降低了扩展目标在交叉处的跟踪误差,对目标的个数和目标的位置实现了精确估计。同时将扩展目标与点目标的跟踪过程分离,大大降低了计算的计算量。
Description
技术领域
本发明属于信息融合领域,涉及一种基于时空关联的多扩展目标跟踪方法。
背景技术
随着现代传感器技术的发展,传感器的分辨率越来越高,单个目标会占据传感器的多个分辨率单元,同时得到多个关于目标的量测信息,该目标类型变为扩展目标。ET-GM-PHD (Extended Target Gaussian Mixture PHD)是在随机有限集的理论下对多扩展目标进行跟踪的滤波算法,该方法使用高斯混合形式表示多目标的强度函数,并以此来近似多目标的后验分布。基于随机有限集的多扩展目标的方法避免了数据关联,并降低了计算的复杂度,其中不可或缺的部分是对当前时刻量测集进行合理的划分,进而找到正确的量测集划分。目前广泛使用的划分方法是距离划分,该方法利用属于同一目标的量测值在空间上相隔较近的特性,对于满足距离门限的量测值都划分到同一个划分单元中,实现对量测集合的划分。但当目标距离较近的时候,不同目标的量测值会落入同一个量测集合中,造成滤波错误。
发明内容
本发明目的在于克服现有方法的不足,提出一种基于时空关联的多扩展目标跟踪方法,解决多扩展目标距离较近和航迹交叉下目标状态和目标数估计问题,具有良好的性能,同时该方法也不受量测密度的影响。本发明主要是应用在杂波环境中单个传感器实现对多扩展目标跟踪的场景下。
本发明的技术方案为:
一种基于时空的多扩展目标跟踪方法,考虑了扩展目标相邻时刻量测值之间的关联性,具体包括:
通过单传感器完成对扩展目标的探测,令k时刻传感器获得的量测集合为 zi为表示二维空间的单个量测值,nk为量测值数目;k时刻传感器检测区域中扩展目标状态集合为其中Mk表示扩展目标的数目,表示单个扩展目标的运动状态,xi,yi表示目标信息,表示目标运动信息。
建立扩展目标的运动方程为:
Xk=FXk-1+vk
扩展目标的量测方程为:
Zk=HXk+wk
基于时空关联的多扩展目标跟踪方法,如图1所示,主要包括以下步骤:
S2、当k≥1时,遍历高斯分量集合,按照状态转移矩阵作一步预测,并添加新生目标的高斯分量,表示为:
S6、对于S5产生的聚类单元,计算聚类单元的并集利用有向图SNN聚类方法进行进一步处理,其聚类数目为K∈[Kl,Ku],其中Kl为DBSCAN算法中的聚类数目,n为集合中量测值的数目,β为扩展目标的量测率;
其中,zs为zi的第个最邻近点,为zj的第个最邻近点,Ni表示zi在kNN图中的k个最邻近点,Nj表示zj在kNN图中的k个最邻近点,zi、zj和zs分别为kNN图中的顶点; Ni∩Nj表示zi和zj之间的共享最临近点,定义为:
式中,zij为一个虚顶点,表示顶点zi和zj是相互的最邻近点。
S62、对于相似度矩阵W和聚类数目范围[Kl,Ku],采用谱聚类的方法进行聚类,形成聚类单元;谱聚类的步骤为:
(1)依据相似度矩阵W计算未归一化的拉普拉斯矩阵L,
式中,D为对角矩阵,为kNN图的度矩阵。
(2)计算归一化的拉普拉斯矩阵Ls=D-1/2LD-1/2;
(3)计算矩阵Ls的特征值λi和特征向量ui;
(4)选取K个最小特征值对应的特征向量构成矩阵U;
(5)归一化矩阵U,得到矩阵Y,并将矩阵Y的每一行当作K维空间新的数据点,利用K-means++算法进行聚类。
S9、选择权重大于0.5的高斯分量作为目标的滤波结果,实现目标跟踪。
本发明的有益效果为:本发明方法降低滤波的计算量,并且能够有效处理扩展目标在交叉处的目标位置估计和目标数目估计。将扩展目标跟踪过程转换为扩展目标跟踪和点目标跟踪的形式,在提高运算效率的同时也使滤波算法有较好的鲁棒性,减少了杂波对估计结果的干扰。
附图说明
图1本发明应用于多扩展目标跟踪的流程图;
图2实施例1下的航迹估计、真实航迹、量测值;
图3实施例1下的OSPA误差对比;
图4实施例1下目标数目估计对比;
图5实施例1下运行时间对比;
图6实施例2下的航迹估计、真实航迹、量测值;
图7实施例2下OSPA误差对比;
图8实施例2下的目标数目估计;
图9实施例2下运行时间对比;
具体实施方式
仿真参数设置如下:目标存活概率为Ps=0.99,检测概率为PD=0.99,扩展目标量测个数服从期望值为10的泊松分布;杂波量测个数服从期望值为20的泊松分布;最大高斯项分量个数为Jmax=100,修剪门限T=10-6,合并门限μ=4;距离门限τ=50。过程噪声的协方差为sv=2;观测噪声方差为sε=20,监视的区域大小为:[-1000,1000]×[-1000,1000]m。
实施例1、
本实施例的目的是验证所提方法在两个扩展目标距离较近且平行运动场景下目标的滤波性能。目标1初始状态为[-600m,-500m,0m/s,0m/s]T,存活时间为1~100s;目标2的初始状态为[-600m,-600m,0m/s,0m/s]T,存活时间为1~100s;此场景下新生目标的强度函数为:
图2是实施例1下本发明的目标状态估计和目标真实状态的对比。可以看出本发明在扩展目标相隔较近的情况下也能得到一个较好的跟踪结果。
图3为实施例1下本发明方法和距离划分、有向图SNN划分方法的OSPA误差对比。可以看出,本发明方法的OSPA误差明显最小。
图4为实施例1下本发明方法和距离划分、有向图SNN划分方法的目标数估计对比。可以看出,本发明方法的目标数估计最为准确。
图5为实施例1下本发明方法和距离划分、有向图SNN划分方法的算法运行时间对比,可以看出,本发明方法的算法运行时间最低。
实施例2、
仿真参数如例1相同。本实施例的目的是验证算法在目标相交情况下的滤波性能。场景中存在2个目标,目标1的初始状态为[250m,250m,0m/s,0m/s]T,存活时间为1~100s;目标2的初始状态为[-250m,-250m,0m/s,0m/s]T,存活时间为1~100s;两个目标在56s 的时候相交。此场景下新生目标的强度函数为:
图6为实施例2下本发明的目标状态估计和目标真实状态的对比。可以看出,即使在目标航迹交叉时刻,也能得到一个较好的估计结果。
图7为实施例2下本发明方法和距离划分、有向图SNN划分方法的OSPA误差对比。可以看出,即使是在扩展目标航迹交叉时刻,本发明方法的OSPA误差也最小。
图8为实施例2下本发明方法和距离划分、有向图SNN划分方法的目标数估计对比。可以看出,在航迹交叉时刻,本发明方法的目标数估计最为准确。
图9为实施例2下本发明方法和距离划分、有向图SNN划分方法的算法运行时间对比,可以看出,本发明方法的算法运行时间最低。
Claims (1)
1.一种基于时空的多扩展目标跟踪方法,其特征在于,包括:
通过单传感器完成对扩展目标的探测,令k时刻传感器获得的量测集合为zi为表示二维空间的单个量测值,nk为量测值数目;k时刻传感器检测区域中扩展目标状态集合为其中Mk表示扩展目标的数目,表示单个扩展目标的运动状态,xi,yi表示目标信息,表示目标运动信息;
建立扩展目标的运动方程为:
Xk=FXk-1+vk
扩展目标的量测方程为:
Zk=HXk+wk
对扩展目标进行跟踪,具体方法为:
S2、当k≥1时,遍历高斯分量集合,按照状态转移矩阵作一步预测,并添加新生目标的高斯分量,表示为:
S6、对于S5产生的聚类单元,计算聚类单元的并集利用有向图SNN聚类方法进行进一步处理,其聚类数目为K∈[Kl,Ku],其中Kl为DBSCAN算法中的聚类数目,n为集合中量测值的数目,β为扩展目标的量测率;
其中,zs为zi的第个最邻近点,为zj的第个最邻近点,Ni表示zi在kNN图中的k个最邻近点,Nj表示zj在kNN图中的k个最邻近点,zi、zj和zs分别为kNN图中的顶点;Ni∩Nj表示zi和zj之间的共享最临近点,定义为:
式中,zij为一个虚顶点,表示顶点zi和zj是相互的最邻近点;
S62、对于相似度矩阵W和聚类数目范围[Kl,Ku],采用谱聚类的方法进行聚类,形成聚类单元;谱聚类的步骤为:
(1)依据相似度矩阵W计算未归一化的拉普拉斯矩阵L,
式中,D为对角矩阵,表示kNN图的度矩阵;
(2)计算归一化的拉普拉斯矩阵Ls=D-1/2LD-1/2;
(3)计算矩阵Ls的特征值λi和特征向量ui;
(4)选取K个最小特征值对应的特征向量构成矩阵U;
(5)归一化矩阵U,得到矩阵Y,并将矩阵Y的每一行当作K维空间新的数据点,利用K-means++算法进行聚类
S9、选择权重大于0.5的高斯分量作为目标的滤波结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110547305.4A CN113344039B (zh) | 2021-05-19 | 2021-05-19 | 一种基于时空关联的多扩展目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110547305.4A CN113344039B (zh) | 2021-05-19 | 2021-05-19 | 一种基于时空关联的多扩展目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113344039A CN113344039A (zh) | 2021-09-03 |
CN113344039B true CN113344039B (zh) | 2022-04-26 |
Family
ID=77469540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110547305.4A Expired - Fee Related CN113344039B (zh) | 2021-05-19 | 2021-05-19 | 一种基于时空关联的多扩展目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113344039B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116256746A (zh) * | 2023-05-12 | 2023-06-13 | 江苏亿连通信技术有限公司 | 基于雷达对防范区域周界空域异物入侵监测***及方法 |
CN117214857B (zh) * | 2023-11-09 | 2024-02-02 | 中国人民解放军海军航空大学 | 一种三维场景下高斯多假设多扩展目标的跟踪方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110297221A (zh) * | 2019-06-19 | 2019-10-01 | 西安电子科技大学 | 一种基于高斯混合模型的数据关联方法 |
CN110298865A (zh) * | 2019-05-22 | 2019-10-01 | 西华大学 | 基于阈值分离聚类器的天基星空背景弱小点目标跟踪方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11408978B2 (en) * | 2015-07-17 | 2022-08-09 | Origin Wireless, Inc. | Method, apparatus, and system for vital signs monitoring using high frequency wireless signals |
CN109671096B (zh) * | 2017-10-13 | 2023-07-28 | 南京航空航天大学 | 一种时空近邻目标检测及网格聚类量测划分下的多扩展目标跟踪方法 |
CN109991570B (zh) * | 2019-04-05 | 2023-02-03 | 电子科技大学 | 一种基于拟牛顿法的运动目标直接定位方法 |
CN110992396A (zh) * | 2019-10-18 | 2020-04-10 | 河南科技大学 | 一种基于K-means++聚类算法的多扩展目标跟踪方法 |
-
2021
- 2021-05-19 CN CN202110547305.4A patent/CN113344039B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110298865A (zh) * | 2019-05-22 | 2019-10-01 | 西华大学 | 基于阈值分离聚类器的天基星空背景弱小点目标跟踪方法 |
CN110297221A (zh) * | 2019-06-19 | 2019-10-01 | 西安电子科技大学 | 一种基于高斯混合模型的数据关联方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113344039A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113344039B (zh) | 一种基于时空关联的多扩展目标跟踪方法 | |
CN107831490A (zh) | 一种改进的多扩展目标跟踪方法 | |
CN110503071B (zh) | 基于变分贝叶斯标签多伯努利叠加模型的多目标跟踪方法 | |
CN103678949B (zh) | 基于密度分析和谱聚类的多扩展目标跟踪量测集划分方法 | |
CN110049549B (zh) | 基于WiFi指纹的多融合室内定位方法及其*** | |
CN111488552B (zh) | 基于高斯混合概率假设密度的紧邻多目标跟踪方法 | |
Qin et al. | Application of an efficient graph-based partitioning algorithm for extended target tracking using GM-PHD filter | |
Cao et al. | Correlation-based tracking of multiple targets with hierarchical layered structure | |
CN107832778A (zh) | 一种基于空间综合相似度的相同目标识别方法 | |
CN111274529B (zh) | 一种鲁棒的高斯逆威沙特phd多扩展目标跟踪算法 | |
CN109375160B (zh) | 纯方位无源定位中一种测角误差估计方法 | |
Yazdani et al. | New similarity functions | |
Leung et al. | Evaluating set measurement likelihoods in random-finite-set slam | |
CN108009272B (zh) | 一种基于有向带权图的低空运动目标搜索方法 | |
CN115619825A (zh) | 地面多目标跟踪状态及轨迹确定方法 | |
CN111273277B (zh) | 一种预测引导的多扩展目标量测集划分方法 | |
Tian et al. | Partition of GB-InSAR deformation map based on dynamic time warping and $ k $-means | |
Yuan et al. | Research on face tracking algorithm based on detection and supervision tracking | |
Chen et al. | A mean shift algorithm based on modified Parzen window for small target tracking | |
Yu et al. | A hypergraph matching labeled multi-Bernoulli filter for group targets tracking | |
Xinglin et al. | Fast density peak-based clustering algorithm for multiple extended target tracking | |
Li et al. | Extended Target Measurement Set Partition Algorithm Based on KFCM Clustering | |
Cai et al. | Multi-sensor Multi-Bernoulli target tracking algorithm based on measurement fusion | |
Gong et al. | Adaptive Soft-Encoding: A General Unsupervised Feature Aggregation Method for Place Recognition | |
Jiang et al. | Marine Targets Spatiotemporal Trajectory Similarity Measurement Method Based on Multidimensional Features |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220426 |
|
CF01 | Termination of patent right due to non-payment of annual fee |