CN109886305B - 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 - Google Patents
一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 Download PDFInfo
- Publication number
- CN109886305B CN109886305B CN201910063618.5A CN201910063618A CN109886305B CN 109886305 B CN109886305 B CN 109886305B CN 201910063618 A CN201910063618 A CN 201910063618A CN 109886305 B CN109886305 B CN 109886305B
- Authority
- CN
- China
- Prior art keywords
- target
- measurement
- time
- state
- targets
- 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
- 238000012883 sequential measurement Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 title claims abstract description 22
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 12
- 238000005259 measurement Methods 0.000 claims abstract description 103
- 230000004927 fusion Effects 0.000 claims abstract description 49
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 15
- 230000002441 reversible effect Effects 0.000 claims abstract description 13
- 239000000203 mixture Substances 0.000 claims description 21
- 238000013138 pruning Methods 0.000 claims description 11
- 239000000126 substance Substances 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 9
- 230000004083 survival effect Effects 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于GM‑PHD滤波的多传感器非顺序量测异步融合方法,可用于解决杂波环境下基于雷达、声呐等多传感器异步非顺序量测的多目标跟踪问题。本发明采取了集中式特征级融合策略,融合中心对实时接收到的量测进行判断,针对顺序量测和到达滞后量测两种异步量测分别设计了基于GM‑PHD滤波器的融合算法。特别的对于到达滞后量测,本发明合理改进了GM‑PHD滤波器,解决了随机集框架下的反向状态预测及负时间量测更新问题,实现了目标状态的二次估计。本发明借助随机集理论的优势,避免了异步融合问题中复杂的数据关联问题,方法结构简单,易于实现的迭代更新,对于解决实际的多传感器多目标目标跟踪问题具有重要的现实意义。
Description
技术领域
本发明属于多目标跟踪及多传感器融合领域,涉及一种多传感器非顺序量测的异步融合方法,可用于解决杂波环境下基于雷达、声呐等多传感器异步量测的多目标跟踪问题。
背景技术
多传感器多目标跟踪技术在军事和民用方面有着广泛的应用,如空中预警,海洋监测,自动驾驶等。随着传感器种类的不断丰富和精度的有效提高,包含数据融合的多目标跟踪方法也在逐步向前发展。
针对多目标跟踪问题,经典解决方法集中在解决数据关联问题上,主要有最近邻算法(Nearest Neighbor,NN)、全局最近邻算法(Global Nearest Neighbor,GNN)、联合概率数据关联(Joint Probability Data Association,JPDA),多假设跟踪(MultipleHypothesis Tracking,MHT)等。然而,最近邻算法受杂波影响较大,标准JPDA的计算复杂度为NP-hard,MHT也存在组合爆炸问题。近年来,基于随机有限集理论(Random Finite Set,RFS)的目标跟踪方法不断引起学术界的广泛关注,相对于传统的多目标跟踪方法,该方法不需要进行数据关联,多目标的状态和量测都以随机集的形式出现,目标的出现、消失、传感器量测受到环境噪声的影响统一用随机集的形式给出。R.Mahler先后提出概率假设密度(Probability Hypothesis Density,PHD)滤波器,带有势分布的PHD滤波器(CardinalizedPHD,CPHD),多目标多伯努利(Multi-target Multi-Bernoulli,MeMBer)滤波器,逐步将随机有限集理论应用于目标跟踪领域。B.-N.Vo给出了PHD、CPHD、MeMBer等滤波器的两种收敛实现,即序贯蒙特卡洛(SequentialMonteCarlo,SMC)和高斯混合(Gaussian mixture,GM)方法。特别是GM-PHD滤波器,使得基于RFS的多目标跟踪方法在工程应用方面成为可能。已有研究表明,GM-CPHD的计算量相比JPDAF的多目标跟踪算法计算复杂度要小的多,并且能够得到更好的多目标跟踪结果。
在多传感器融合理论中,研究更多的是同步融合,即各传感器同步测量,且同步传输数据到融合中心。而实际中由于各传感器采样频率不同以及通信滞后等问题,到达融合中心的数据是异步的。异步量测主要分为顺序量测(in-sequence measurement,ISM)和非顺序量测(out-of-sequence measurement,OOSM),两者的差别主要在于量测到达融合中心的顺序与其时间戳的顺序是否一致,如图1所示。针对顺序量测异步融合问题,只需进行序贯滤波即可,而非顺序量测的异步融合则较为复杂。针对一步滞后非顺序量测问题,有学者先后提出基于向后状态倒推的最优A1算法、次优的B1、C1算法,这些算法各有优劣,在不同场景下能达到有效估计。上述算法主要针对单目标下的异步信息融合,由于涉及数据关联问题,多目标异步信息融合时将更为复杂,目前关于多传感器多目标异步融合问题的研究较少,且尚未有在随机集框架下非顺序量测异步融合的应用研究。
发明内容
本发明针对多传感器非顺序量测(OOSM)异步融合问题,研究设计了一种基于高斯混合概率假设密度滤波器(GM-PHD)的多传感器融合方法。方法借助随机集理论的优势,避免了异步融合问题中较为棘手的数据关联问题。本发明的技术方案为:
首先,建立多传感器多目标跟踪场景及基础模型,包括目标运动模型,量测模型、杂波模型等,并初始化。
其次,建立多传感器非顺序量测异步融合框架。该框架对异步量测进行集中式特征级融合。假设有s个同质传感器对目标进行探测,并将量测传给融合中心。融合中心每接到一个传感器的量测,便将量测时刻与当前最新估计结果时刻进行对比,当量测为顺序量测时,直接应用高斯混合PHD滤波器进行剪枝合并、预测、更新,得到新的估计结果;当量测为到达滞后量测时,应用基于直接更新法的高斯混合PHD滤波器,先进行反向状态预测,然后借助滞后量测对当前估计结果二次滤波。
最后,在需要获得目标状态估计的时刻,融合中心对当前最新估计结果进行目标个数及状态提取并输出。
本发明的创新点在于:
提出了基于高斯混合概率假设密度滤波器(GM-PHD)的多传感器异步融合方法,能对顺序量测和到达滞后量测两种异步数据进行有效处理,得到最新时刻多目标状态的估计值;对GM-PHD进行合理改进,以用于反向状态预测及非同步量测的状态直接更新问题,实现了非顺序量测下的目标状态估计;借助随机集理论的优势,避免了异步融合问题中复杂的数据关联问题。方法结构简单,易于实现的迭代更新,并保证了目标跟踪误差的收敛性。
附图说明
图1为异步非顺序量测示意图。
图2为本发明所提出的基于GM-PHD滤波的多传感器非顺序量测异步融合方法的流程图。
图3为仿真实验中三种不同测量情形。
图4为一次蒙特卡洛仿真中目标的真实轨迹及测量值。
图5为三种场景下100次蒙特卡洛仿真的跟踪估计结果。
图6为三种情况下100次蒙特卡洛仿真跟踪结果的平均Hausdoff距离和平均Wasserstein距离对比图。
图7为三种情况下100次蒙特卡洛仿真目标数目估计平均值。
具体实施方式
下面结合附图对本发明的实施进行详细的说明,并给出具体的操作方式以及实施步骤。
步骤一:建立多传感器多目标跟踪模型,并初始化
在随机有限集框架下,多目标的状态和量测均可视为随机有限集,假设tk时刻有N(k)个状态和M(k)个量测,多目标的状态和量测可以表示为有限集形式:
Xk={xk,1,xk,2,...,xk,N(k)} (1)
Zk={zk,1,zk,2,...,zk,N(k)} (2)
其中,xk,1,xk,2,...,xk,N(k)分别为第1、2、...、N(k)个目标状态,zk,1,zk,2,...,zk,N(k)分别为第1、2、...、M(k)个量测。
对于单个目标来说,运动模型为:
xk=F(k,k-1)xk-1+v(k,k-1) (3)
其中,为目标状态变量,其中的元素分别代表x方向的位置,x方向的速度,y方向的位置,y方向的速度。F(k,k-1)是目标从tk-1时刻到k时刻的状态转移矩阵,v(k,k-1)是从k-1时刻到tk-1时刻的状态误差累积过程,服从N(0,Q(k,k-1))的分布。
对于单个量测来说,量测模型为:
zk=H(k)xk+w(k) (4)
其中,zk为tk-1时刻的一个量测,H(k)为tk-1时刻量测函数,w(k)是分布为N(0,R(k))的量测噪声。
对于tk-1时刻多目标状态Xk-1,每个目标的状态xk-1∈Xk-1或者以ps,k+1的概率继续存活,变成新的状态xk,或者以1-ps,k+1的概率消失。所以对于给定的Xk-1,tk时刻的多目标状态Xk由继续生存的目标Sk|k-1(x)、和新生的目标Γk构成:
量测方面,对于给定目标或者被检测到,或者漏检,故每个目标x将产生一个随机有限集量测Θk(x)。除了源于目标的量测,传感器也可能会接收到杂波或虚警量测集Kk。所以多目标的量测Zk由真实量测和杂波构成,即:
在高斯混合概率密度假设下,起始时刻目标表示为高斯混合的形式:
杂波模型是强度为κk(z)的泊松有限集,即
κk(z)=λcVu(z) (8)
其中,u(z)表示监视区域的均匀概率分布,λc为杂波强度,V为传感器监视区域的面积。
步骤二:建立多传感器非顺序量测异步融合框架
如图2,本发明设计了一种多传感器非顺序量测异步融合框架,该框架对异步量测进行集中式特征级融合。假设有s个同质传感器对目标进行探测,所有传感器采样频率相同,融合中心以相同的频率输出目标状态的估计结果。导致传感器量测异步的原因主要有两个:各传感器探测时间未同步;由于信号传递时延不同,即便相同时间戳的量测到达融合中心的时间也不同。假设融合中心依次接收到来自不同传感器的量测,相应的测量时刻(时间戳)分别为:t1,t2,...,tk-1,tk,...,ts(其中下标代表接收顺序),融合中心当前最新估计值对应的时刻为tk,将其作为下一帧量测的比较基准,用Tk表示(Tk=tk),每接收到一个传感器的量测集合,融合中心依据其时间戳判断该量测集合属于顺序量测或到达滞后量测,然后采取相应的融合策略。即针对到达融合中心的量测的真实时刻tk+1,通过与当前最新估计时刻tk比较分为两种情况,顺序量测(tk+1>Tk),进行下述的剪枝合并策略,并执行步骤三;到达滞后量测(tk+1<Tk),不进行剪枝合并,直接执行步骤四。
(1)剪枝合并
假设tk时刻目标的后验强度可以表示为高斯混合的形式:
在高斯混合PHD滤波过程中,随着时间的推移,高斯分量的不断增加,因此需要对权重较小的高斯分量进行减枝;另外,一些高斯分量非常接近,需要合并为一个高斯分量。通常剪枝合并步骤在滤波器预测更新之后,但是当下一帧量测为滞后量测时,需用到剪枝合并前的相关信息,故将剪枝合并放在第二步,即确定量测为顺序量测时进行,到达滞后量测时不进行。设Tw为减枝阈值,U为合并阈值,Jmax为最大允许的高斯项数目,k表示剪枝合并的时刻。
循环:
步骤三:顺序量测(tk+1>Tk),依照GM-PHDF处理
(1)预测
假设tk时刻目标的后验强度可以表示为高斯混合的形式:
式中x为目标状态空间的任意值,Jk表示tk时刻目标后验强度的高斯分量个数,分别为第i个高斯分量的权重、均值和协方差。tk+1时刻目标状态集合可表示为源于tk时刻的存活目标Sk+1|k(x)和当前时刻的新生目标Γk+1,
即:
所以,tk+1时刻目标的先验强度可以表示为:
vk+1|k(x)=vs,k+1|k(x)+γk+1(x) (12)
其中,γk+1(x)为新生目标的强度,此处假设为高斯混合形式,即:
vs,k+1|k(x)为存活部分的先验高斯分布,可表示为:
(2)更新
如(1)中所述,tk+1时刻目标的先验强度可以统一表示为高斯混合形式:
得到tk+1时刻多目标PHD的预测后,由下式对PHD进行更新:
式中等号右面第一项为未被检测到的目标,直接用先验高斯分量代替,PD,k+1表示目标的检测概率,第二项代表对于每一个量测z∈Zk+1,都有对应的更新高斯分量,即:
式中对单个高斯分量使用了卡尔曼滤波的更新方程。具体步骤如下:
(3).最后,由于获得了更新时刻目标的状态估计,故Tk需要赋予新值,即Tk=tk+1
步骤四:一步滞后量测(tk+1<Tk)处理
同步骤三,tk时刻目标的后验强度可以表示为高斯混合的形式:
当融合中心接到的下一个传感器在tk+1(tk+1<Tk)时刻的量测集合Zk+1时,新的滤波时刻仍为tk(为强调量测的滞后性,这里用Zd代替滞后量测Zk+1,用td代替滞后量测的真实时刻tk+1,滞后的相关量均带有下标d)。本发明主要基于直接更新法中的次优方法(B1法),利用Zd对tk时刻目标状态直接进行二次更新,得到新的高斯分量
由步骤三可知,tk时刻的目标主要分为两类:
第一项为未被检测到的目标,直接用先验高斯分量代替,并利用目标未被检测到的概率(1-PD,k)对其权重进行修正,第二项代表被检测到的目标,且每一个目标都用所有量测z∈Zk对其先验高斯分量进行了更新。本发明对于两种目标分别进行了反向状态预测及二次更新。
1.tk时刻未被检测到的目标(先验目标)
(1).反向状态预测
未被检测到的目标,用先验高斯分量作为后验高斯分量,假设量测一步滞后,td、tk两时刻相距不远,td时刻目标状态集合可简化表示为tk时刻的所有目标:
所以,td时刻目标的先验强度可以为:
(2).滞后量测对最新目标状态的二次更新:
如(1)中所述,td时刻目标的先验强度可以统一表示为高斯混合形式:
得到td时刻多目标PHD的预测后,依照卡尔曼滤波算法,对于每一个由下式直接对PHD进行二次更新:
式中等号右面第一项代表td时刻未被检测到的目标,直接用原有后验高斯分量代替,PD,d表示td时刻目标的检测概率,第二项代表对于每一个量测z∈Zd,都有对应的更新高斯分量,即:
式中对单个高斯分量使用了卡尔曼滤波的更新方程。具体步骤如下:
因此,更新后的高斯分量的均值和协方差为:
2.tk时刻被检测到的目标(后验目标)
(1).反向状态预测
与1中先验目标相似,td时刻后验目标状态集合可简化表示为tk时刻的所有目标:
所以,td时刻目标的先验强度可以为:
(2).滞后量测对最新目标状态的二次更新:
如(1)中所述,td时刻目标的先验强度可以统一表示为高斯混合形式:
得到td时刻多目标PHD的预测后,依照卡尔曼滤波算法,对于每一个由下式直接对PHD进行二次更新:
式中等号右面第一项代表td时刻未被检测到的目标,直接用原有后验高斯分量代替,PD,d表示td时刻目标的检测概率,第二项代表对于每一个量测z∈Zd,都有对应的更新高斯分量,即:
式中对单个高斯分量使用了卡尔曼滤波的更新方程。具体步骤如下:
因此,更新后的高斯分量的均值和协方差为:
步骤五:目标个数及状态提取
在需要获得目标状态估计的时刻,融合中心对当前最新估计结果进行目标个数及状态提取,并输出结果。也可以由融合中心依据一定的输出频率进行该步骤。这里选择权重大于0.5的高斯分量进行输出,具体方法如下:
仿真实验对比了三种不同测量情形下的融合方法,如图3所示。
(a)顺序量测序贯PHD滤波器(ISM-PHD):两个传感器的量测到达融合中心顺序与采样顺序一致,融合中心基于PHD滤波器对异步测量数据进行序贯融合。
(b)非顺序量测PHD滤波器(OOSM-PHD):两个传感器的量测到达融合中心顺序与采样顺序相反,融合中心采取基于PHD滤波器的非顺序量测融合方法。
(c)丢帧情况下PHD滤波器:与(b)相同,两个传感器的量测到达融合中心顺序与采样顺序相反,但融合中心丢掉滞后量测,只利用顺序量测进行目标状态估计。
采取100次蒙特卡洛仿真实验,每次仿真目标真实轨迹保持不变,量测重新生成。图4为一次蒙特卡洛仿真中目标的真实轨迹及测量值,其余图皆为100次蒙特卡洛实验结果。图5为三种场景下跟踪估计结果,可以看出非顺序量测PHD滤波器的跟踪结果与基于顺序量测效果相当,而(c)中丢掉非顺序量测使得估计结果散乱,精度较差。图6为三种情况下跟踪结果的Hausdoff距离和Wasserstein距离对比图,图7为三种情况下目标数目估计值,分别从不同角度进一步说明所提出的非顺序量测PHD滤波器较为稳定,达到了和基于顺序量测的算法相近的效果,而直接将滞后量测丢掉的做法由于量测数目缺失,造成跟踪精度较低,结果易发散。
Claims (3)
1.一种基于GM-PHD滤波的多传感器非顺序量测异步融合方法,其特征在于包括如下步骤:
1)建立多传感器多目标跟踪模型,包括目标运动模型、量测模型、杂波模型,并初始化;
2)建立多传感器非顺序量测异步融合框架,融合中心当前最新的目标估计值对应的时刻为tk,将其作为下一帧量测的比较基准,用Tk表示(Tk=tk),当接收到一帧传感器的量测,融合中心将量测时间戳tk+1与Tk比较,判断该量测属于顺序量测或到达滞后量测,然后采取如下融合策略:对于顺序量测,即tk+1>Tk,应用高斯混合PHD滤波器进行剪枝合并、预测、更新,得到新的估计结果;对于到达滞后量测,即tk+1<Tk,应用基于直接更新法的高斯混合PHD滤波器,先进行反向状态预测,然后借助滞后量测对当前估计结果二次滤波,具体为:
tk时刻的目标主要分为两类:
第一项为未被检测到的目标,直接用先验高斯分量代替,并利用目标未被检测到的概率(1-PD,k)对其权重进行修正,这里称为先验目标;第二项代表被检测到的目标,且每一个目标都用所有量测z∈Zk对其先验高斯分量进行了更新,这里称为后验目标,两种目标分别进行反向状态预测及二次更新;
所述的两种目标分别进行反向状态预测及二次更新,具体为:
(a).tk时刻未被检测到的目标,即先验目标
(a1).反向状态预测
未被检测到的目标,用先验值作为后验高斯分量,td时刻先验目标状态集合可简化表示为tk时刻的所有目标:
所以,td时刻目标的先验强度可以为:
(a2).滞后量测对最新目标状态的二次更新:
如(a1)中所述,td时刻目标的先验强度可以统一表示为高斯混合形式:
得到td时刻多目标PHD的预测后,依照卡尔曼滤波算法,对于每一个由下式直接对PHD进行二次更新:
式中等号右面第一项代表td时刻未被检测到的目标,直接用原有后验高斯分量代替,PD,d表示td时刻目标的检测概率,第二项代表对于每一个量测z∈Zd,都有对应的更新高斯分量,即:
式中对单个高斯分量使用了卡尔曼滤波的更新方程,具体步骤如下:
因此,更新后的高斯分量的均值和协方差为:
(b).tk时刻被检测到的目标,即后验目标
(b1).反向状态预测
与(a1)中先验目标相似,td时刻后验目标状态集合可简化表示为tk时刻的所有目标:
所以,td时刻目标的先验强度可以为:
(b2).滞后量测对最新目标状态的二次更新:
如(b1)中所述,td时刻目标的先验强度可以统一表示为高斯混合形式:
得到td时刻多目标PHD的预测后,依照卡尔曼滤波算法,对于每一个由下式直接对PHD进行二次更新:
式中等号右面第一项代表td时刻未被检测到的目标,直接用原有后验高斯分量代替,PD,d表示td时刻目标的检测概率,第二项代表对于每一个量测z∈Zd,都有对应的更新高斯分量,即:
式中对单个高斯分量使用了卡尔曼滤波的更新方程,具体步骤如下:
更新后的高斯分量的均值和协方差为:
3)在需要获得目标状态估计的时刻,融合中心对当前最新估计结果进行目标个数及状态提取并输出。
2.如权利要求1所述的基于GM-PHD滤波的多传感器非顺序量测异步融合方法,其特征在于所述的步骤1)具体为:
1.1)建立多传感器多目标跟踪模型,并初始化;
用Xk={xk,1,xk,2,...,xk,M(k)}表示tk时刻多目标状态构成的随机有限集,用Zk={zk,1,zk,2,...,zk,N(k)}表示tk时刻量测构成的随机有限集,其中,xk,1,xk,2,...,xk,M(k)分别为第1、2、...、M(k)个目标状态,zk,1,zk,2,...,zk,N(k)分别为第1、2、...、N(k)个量测;
对于单个目标来说,运动模型为:
xk=F(k,k-1)xk-1+v(k,k-1)
其中,为目标状态变量,其中的元素分别代表x方向的位置,x方向的速度,y方向的位置,y方向的速度;F(k,k-1)是目标从tk-1时刻到tk时刻的状态转移矩阵,v(k,k-1)是从tk-1时刻到tk时刻的状态误差累积过程,服从N(0,Q(k,k-1))的分布;
对于单个量测来说,量测模型为:
zk=H(k)xk+w(k)
其中,zk为tk时刻的一个量测,H(k)为tk时刻量测函数,w(k)是分布为N(0,R(k))的量测噪声;
给定tk-1时刻多目标状态Xk-1,tk时刻的多目标状态Xk由继续生存的目标Sk|k-1(x)、和新生的目标Γk构成:
量测方面,多目标的量测Zk由真实量测Θk(x)和杂波Kk构成,即:
起始时刻,目标表示为高斯混合的形式:
杂波模型是强度为κk(z)的泊松有限集,即:
κk(z)=λcVu(z)
其中,u(z)表示监视区域的均匀概率分布,λc为杂波强度,V为传感器监视区域的面积。
3.如权利要求1所述的基于GM-PHD滤波的多传感器非顺序量测异步融合方法,其特征在于所述的步骤2中的,剪枝合并步骤在预测、更新之前,且剪枝合并步骤仅在顺序量测时进行,到达滞后量测时不进行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910063618.5A CN109886305B (zh) | 2019-01-23 | 2019-01-23 | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910063618.5A CN109886305B (zh) | 2019-01-23 | 2019-01-23 | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109886305A CN109886305A (zh) | 2019-06-14 |
CN109886305B true CN109886305B (zh) | 2021-05-04 |
Family
ID=66926510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910063618.5A Active CN109886305B (zh) | 2019-01-23 | 2019-01-23 | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109886305B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110596693B (zh) * | 2019-08-26 | 2021-10-22 | 杭州电子科技大学 | 一种迭代更新的多传感器gmphd自适应融合方法 |
CN110807942B (zh) * | 2019-09-24 | 2021-11-02 | 上海汽车工业(集团)总公司 | 智能驾驶汽车航迹更新方法及其更新系统 |
CN110703272B (zh) * | 2019-09-27 | 2022-05-17 | 江苏大学 | 一种基于车车通信和gmphd滤波的周边目标车辆状态估计方法 |
CN110646790B (zh) * | 2019-09-30 | 2021-08-03 | 电子科技大学 | 一种用于雷达组网失序量测集中式融合的目标跟踪方法 |
CN111127523B (zh) * | 2019-12-04 | 2023-03-24 | 杭州电子科技大学 | 基于量测迭代更新的多传感器gmphd自适应融合方法 |
CN111325776A (zh) * | 2020-02-19 | 2020-06-23 | 江苏理工学院 | 一种基于变分贝叶斯t分布卡尔曼滤波的phd多目标跟踪方法 |
CN111563918B (zh) * | 2020-03-30 | 2022-03-04 | 西北工业大学 | 一种多卡尔曼滤波器数据融合的目标跟踪方法 |
CN111736145B (zh) * | 2020-06-28 | 2022-04-19 | 电子科技大学 | 一种基于高斯混合概率假设密度滤波的多机动目标多普勒雷达跟踪方法 |
US11954179B2 (en) * | 2020-09-11 | 2024-04-09 | Inceptio Hongkong Limited | Resolving out of sequence data in sensor tracking |
CN112766151B (zh) * | 2021-01-19 | 2022-07-12 | 北京深睿博联科技有限责任公司 | 一种用于导盲眼镜的双目目标检测方法和系统 |
CN113987980B (zh) * | 2021-09-23 | 2022-05-20 | 北京连山科技股份有限公司 | 一种图形物理phd流行仿真实现方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182291A (zh) * | 2015-08-26 | 2015-12-23 | 西安电子科技大学 | 自适应目标新生强度的phd平滑器的多目标跟踪方法 |
CN105716620A (zh) * | 2016-03-16 | 2016-06-29 | 沈阳建筑大学 | 一种基于云计算和大数据的导航方法 |
CN107103331A (zh) * | 2017-04-01 | 2017-08-29 | 中北大学 | 一种基于深度学习的图像融合方法 |
CN107783104A (zh) * | 2017-10-17 | 2018-03-09 | 杭州电子科技大学 | 一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 |
CN108333569A (zh) * | 2018-01-19 | 2018-07-27 | 杭州电子科技大学 | 一种基于phd滤波的异步多传感器融合多目标跟踪方法 |
CN108344981A (zh) * | 2018-01-19 | 2018-07-31 | 杭州电子科技大学 | 面向杂波的多传感器异步检测tsbf多目标跟踪方法 |
CN108632764A (zh) * | 2018-04-02 | 2018-10-09 | 浙江大学 | 一种多传感器选择性量测数据融合估计算法 |
CN109031229A (zh) * | 2017-11-27 | 2018-12-18 | 电子科技大学 | 一种杂波环境下目标跟踪的概率假设密度方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201004488D0 (en) * | 2010-03-17 | 2010-05-05 | Isis Innovation | A method of tracking targets in video data |
FR3017481B1 (fr) * | 2014-02-07 | 2016-02-26 | Sagem Defense Securite | Procede de detection et de pistage de cibles |
CN106372646B (zh) * | 2016-08-30 | 2020-07-14 | 上海交通大学 | 基于srck-gmcphd滤波的多目标跟踪方法 |
CN107423539B (zh) * | 2017-03-29 | 2020-07-24 | 中国电子科技集团公司第三十八研究所 | 一种gm-phd滤波器的标签分配方法 |
-
2019
- 2019-01-23 CN CN201910063618.5A patent/CN109886305B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182291A (zh) * | 2015-08-26 | 2015-12-23 | 西安电子科技大学 | 自适应目标新生强度的phd平滑器的多目标跟踪方法 |
CN105716620A (zh) * | 2016-03-16 | 2016-06-29 | 沈阳建筑大学 | 一种基于云计算和大数据的导航方法 |
CN107103331A (zh) * | 2017-04-01 | 2017-08-29 | 中北大学 | 一种基于深度学习的图像融合方法 |
CN107783104A (zh) * | 2017-10-17 | 2018-03-09 | 杭州电子科技大学 | 一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 |
CN109031229A (zh) * | 2017-11-27 | 2018-12-18 | 电子科技大学 | 一种杂波环境下目标跟踪的概率假设密度方法 |
CN108333569A (zh) * | 2018-01-19 | 2018-07-27 | 杭州电子科技大学 | 一种基于phd滤波的异步多传感器融合多目标跟踪方法 |
CN108344981A (zh) * | 2018-01-19 | 2018-07-31 | 杭州电子科技大学 | 面向杂波的多传感器异步检测tsbf多目标跟踪方法 |
CN108632764A (zh) * | 2018-04-02 | 2018-10-09 | 浙江大学 | 一种多传感器选择性量测数据融合估计算法 |
Non-Patent Citations (2)
Title |
---|
Iterative RANSAC based adaptive birth intensity estimation in GM-PHD filter for multi-target tracking;Jingjing Wu 等;《Signal Processing》;20170228;第412-421页 * |
多目标无序量测OOSM-GMPHD滤波算法;赵凯 等;《传感器与微系统》;20171231;第0012-0067段 * |
Also Published As
Publication number | Publication date |
---|---|
CN109886305A (zh) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109886305B (zh) | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 | |
CN105182291B (zh) | 自适应目标新生强度的phd平滑器的多目标跟踪方法 | |
CN108333569B (zh) | 一种基于phd滤波的异步多传感器融合多目标跟踪方法 | |
CN110542885B (zh) | 一种复杂交通环境下的毫米波雷达目标跟踪方法 | |
CN108344981B (zh) | 面向杂波的多传感器异步检测tsbf多目标跟踪方法 | |
CN110596693B (zh) | 一种迭代更新的多传感器gmphd自适应融合方法 | |
CN111127523B (zh) | 基于量测迭代更新的多传感器gmphd自适应融合方法 | |
CN111580084B (zh) | 一种面向多距离扩展目标的多伯努利检测前跟踪方法 | |
Ristic et al. | Particle filter for joint estimation of multi-object dynamic state and multi-sensor bias | |
CN107436434B (zh) | 基于双向多普勒估计的航迹起始方法 | |
CN108717702B (zh) | 基于分段rts的概率假设密度滤波平滑方法 | |
CN108717174A (zh) | 基于信息论的预测快速协方差交互融合无源协同定位方法 | |
Li et al. | A parallel filtering-communication-based cardinality consensus approach for real-time distributed PHD filtering | |
CN108732564A (zh) | 一种双雷达修正序贯高斯混合概率假设密度滤波方法 | |
CN117169818B (zh) | 一种面向海面机动目标跟踪的雷达波形设计方法 | |
CN109239704A (zh) | 一种基于序贯滤波交互式多模型的自适应采样方法 | |
CN111340853B (zh) | 基于ospa迭代的多传感器gmphd自适应融合方法 | |
Koch | On optimal distributed Kalman filtering and retrodiction at arbitrary communication rates for maneuvering targets | |
Khan et al. | A fixed lag smoothing IPDA tracking in clutter | |
CN107330468B (zh) | 基于切尔诺夫融合准则的多传感器分布式数据融合方法 | |
CN104467742A (zh) | 基于高斯混合模型的传感器网络分布式一致性粒子滤波器 | |
Chen | An algorithm of mobile sensors data fusion tracking for wireless sensor networks | |
CN115619825A (zh) | 地面多目标跟踪状态及轨迹确定方法 | |
Liu et al. | Asynchronous multi-sensor fusion multi-target tracking method | |
CN115220002A (zh) | 一种固定单站的多目标数据关联跟踪方法和相关装置 |
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 |