CN111580084A - 一种面向多距离扩展目标的多伯努利检测前跟踪方法 - Google Patents
一种面向多距离扩展目标的多伯努利检测前跟踪方法 Download PDFInfo
- Publication number
- CN111580084A CN111580084A CN202010403701.5A CN202010403701A CN111580084A CN 111580084 A CN111580084 A CN 111580084A CN 202010403701 A CN202010403701 A CN 202010403701A CN 111580084 A CN111580084 A CN 111580084A
- Authority
- CN
- China
- Prior art keywords
- target
- bernoulli
- distance
- component
- density
- 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
Links
Images
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
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
- G06F18/24155—Bayesian classification
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Probability & Statistics with Applications (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于目标检测和跟踪技术领域,公开了一种面向多距离扩展目标的多伯努利检测前跟踪方法。本发明通过:雷达系统处理接收回波,得到距离‑方位回波幅度测量图像;对多目标状态密度分布进行预测;通过时域滑窗对距离目标历史测量数据进行储存,构造目标信息矩阵,并借助主成分分析对幅度分布信息进行估计;根据目标幅度估计和测量图像,结合预测密度对多目标状态进行更新,得到多目标后验密度;设置最大伯努利成份数、合并和修剪阈值,对多目标航迹进行管理,实现多目标检测前跟踪。本发明可以在线估计多个目标的幅度信息,有效适应不同分布模型雷达图像观测,保证检测与跟踪的自适应能力,实现时变数目多距离扩展目标的鲁棒检测和稳定跟踪。
Description
技术领域
本发明涉及目标检测和跟踪技术领域,更具体地说,涉及一种面向多距离扩展目标的多伯努利检测前跟踪方法
背景技术
在民用监控与军事侦察需求日益深化的背景下,要求传感器具有对多个扩展目标的快速检测与稳健跟踪能力。在传统单目标跟踪方法中,场景内至多存在一个目标,不需要使用复杂的数据关联对源于目标或杂波的量测进行划分和关联。而实际应用场景中,通常需要跟踪多个目标。传统多目标跟踪数据关联方法中,基于单帧点迹的联合概率数据关联(JPDA)和基于多帧点迹的多假设数据关联(MHT)都能够在信噪比较高的条件下解决目标和点迹量测匹配的问题。在低信噪比情况下,为避免弱目标出现漏检,检测器降低检测阈值,导致虚警密度增高。由于MHT是一种基于穷举的关联方法,其复杂度会随着点迹量测数目的增加而爆炸性增长。基于随机有限集的滤波器方法将多目标状态建模为随机有限集,不需要进行数据关联,能够较好避免组合爆炸问题。在随机有限集(Random Finite Set,RFS)框架下,概率假设密度(PHD)滤波器和势均衡概率假设密度(CPHD)滤波器最早被提出来解决多目标联合估计数目与状态。与PHD和CPHD这些矩近似滤波器不同,多目标多伯努利(MeMBer)滤波器是一种概率密度近似方法,实现过程与传统滤波器相似,不需要设计额外的多目标状态提取过程。此外,CPHD滤波器的高斯混合实现尽管优于多目标多伯努利滤波器,但是多伯努利滤波器的粒子近似实现却能取得更优的结果。
传统跟踪方法的点迹量测仅提取目标的运动状态信息而不包含回波幅度信息,Lerro(Lerro,D.,Bar-Shalom,Y.:‘Automated tracking with target amplitudeinformation’.American Control Conf.,San Diego,CA,USA,1990,pp.2875–2880)引入了带幅度信息的点迹量测,提出了单目标概率数据关联滤波器,利用运动似然和幅度似然区分目标和杂波,Keuk(Van Keuk,G.:‘Multihypothesis tracking using incoherentsignal-strength information’,IEEE Trans.Aerosp.Electron.Syst.,1996,32,(3),pp.1164–1170)将该方法应用到多目标多假设跟踪滤波器。有学者将利用幅度特征辅助跟踪方法嵌入到随机有限集滤波器中,(Ehrman,L.M.,Blair,W.D.:‘Comparison of methodsfor using target amplitude to improve measurement-to-track association inmulti-target tracking’.Proc.Int.Conf.Information Fusion,Florence,Italy,2006,pp.1–8)对比了使用目标幅度信息辅助的多目标跟踪方法,但所讨论的方法只适用于目标幅度先验已知或者是固定的情况。在实际应用中,非合作目标的实际回波幅度信息是未知且变化的。Bae分别针对单帧观测数据提出序贯蒙特卡洛(SMC)(Bae,S.-H.,Kim,D.Y.,Yoon,J.H.,et al.:‘Automated multi-target tracking with kinematic and non-kinematic information’,IET Radar Sonar Navig.,2012,6,(4),pp.272–281)方法和针对多帧观测数据提出最大似然估计(MAP)(Bae,S.-H.,Park,J.Y.,Yoon,K.-J.:‘Jointestimation of multi-target SNR and dynamic states in cluttered environment’,IET Radar Sonar Navig.,2017,11,(3),pp.539–549)方法在线估计目标的回波平均幅度。而针对传感器图像观测数据,Rollason(M.Rollason,D.Salmond,A particle filter fortrack-before-detect of a target with unknown amplitude,IET,2002,pp.14–19)首次在检测前跟踪(TBD)框架下提出基于未知目标回波幅度的单目标贝叶斯粒子滤波器,将目标平均回波幅度建模为状态变量,联合估计目标状态及平均回波幅度,并由Boers(BoersY,Driessen J N.Multitarget particle filter track before detect application[J].Radar Sonar&Navigation IEE Proceedings,2004,151(6):351-357)引入到多目标场景。相较于简单的点扩散函数,针对具有纹理信息的结构化图像数据,Rollason(MalcolmR,David S.Particle filter for track-before-detect of a target with unknownamplitude viewed against a structured scene[J].IET Radar,Sonar&Navigation,2018,12(6):603-609)提出改进的TBD粒子滤波。Liu(Chao Liu,Jinping Sun ID,PengLei,et al.δ-Generalized labeled multi-Bernoulli filter using amplitudeinformation of neighboring cells[J].Sensors.2018)在将幅度特征辅助方法加入到广义标签多伯努利滤波器并取得良好的跟踪效果。然而,上述方法都是基于点目标假设的,在高分辨率传感器中显然不适用,目前鲜有研究报道在线估计多个距离扩展目标平均回波幅度的跟踪方法。本发明针对未知幅度分布的多个距离扩展目标提出了一种基于主成分分析的幅度估计方法,实现了多个距离扩展目标的联合检测与跟踪。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供了一种面向多距离扩展目标的多伯努利检测前跟踪方法,该方法利用一个滑动窗口,对多帧雷达二维距离-方位观测图像数据进行缓存操作,并引入主成分分析工具,对各个目标的幅度分布信息进行估计,同时结合多伯努利滤波器在多目标跟踪问题上的优势,可以在线估计多个目标的幅度分布信息,可以有效适应不同分布模型雷达图像观测数据,保证检测与跟踪的自适应能力,实现时变数目多距离扩展目标的鲁棒检测和稳定跟踪。
本发明所采用的技术方案是:
一种面向多距离扩展目标的多伯努利检测前跟踪方法,包括以下步骤:
1)雷达系统处理接收回波信号,得到距离-方位回波幅度测量图像;
2)根据目标运动模型,对多目标状态密度分布进行预测;
3)通过时域滑窗对距离扩展目标的历史测量数据进行储存,构造目标信息矩阵(Target Information Matrix,TIM),并借助主成分分析(Principal ComponentAnalysis,PCA)工具对其幅度分布信息进行估计;
4)根据目标幅度估计结果和测量图像,结合预测密度对多目标状态进行更新,得到多目标后验密度;
5)设置最大伯努利成份数以及合并和修剪阈值,对多目标航迹进行管理,实现多目标检测前跟踪。
作为本发明的进一步改进,所述步骤1)通过以下步骤完成:
1.1)雷达系统发射大时宽大带宽信号,通过距离脉冲压缩得到高分辨距离像;
1.2)雷达系统对监测区域的不同角度进行扫描得到不同方位角度的高分辨距离像,一次扫描得到一幅包含Nr×Na个分辨单元的雷达距离-方位回波幅度测量图像,其中Nr为距离分辨单元数,Na为方位角度分辨单元数。
作为本发明的进一步改进,所述步骤2)通过以下步骤完成:
2.1)选取目标的远端,即目标上距离雷达最远端的点作为目标的参考点,如图1所示,并根据先验知识给定系统模型,以下给定先验知识:将k-1时刻的多目标状态表示为Xk-1,其中的每一个目标对应的状态表示为xk-1,则对应的后验多目标密度为其中和分别为目标的存活概率和概率密度分布,Mk-1为k-1时刻的伯努利成分数;k时刻,目标将以pS,k(xk-1)的概率转移到xk,并且状态转移概率密度为ζk|k-1(xk|xk-1),或者以1-pS,k(xk-1)的概率从场景中消失,相应的数学模型表示为{pS,k(xk-1),ζk|k-1(xk|xk-1)}的伯努利RFSSk|k-1(xk-1),由此将预测过程分为两部分:新生(newbirth)部分和持续存活(survival,也称为幸存部分)部分;
2.4)以持续存活目标的伯努利RFSSk|k-1(xk-1)和新生目标的伯努利RFS Γk的并集来描述k时刻的多目标预测状态其中多伯努利RFS中的伯努利RFS成分满足相互独立假设;多伯努利贝叶斯预测过程通过预测成分和新生成分的并来表示:Mk|k-1=Mk-1+MΓ,k是预测多伯努利RFS的伯努利成分数。
作为本发明的进一步改进,所述步骤3)通过以下步骤完成:
3.1)目标TIM通过一个滑动窗口,积累连续多帧图像测量数据中距离扩展目标所影响的图像区域分辨单元的强度数据,其中目标的图像影响范围由T(xk)描述,xk为目标状态,则多帧滑窗缓存获取的k时刻TIM的数学描述符号为Ak={Aτ|τ=k-Wl,…,k-1},其中 是k时刻第tk个分辨单元的强度测量数据,τ为时间索引,Wl是时间滑窗的长度;针对状态为xk的距离扩展目标,其图像影响范围T(xk)由下式进行描述: 其中rk(xk)和ak(xk)分别为目标相对雷达的径向长度和角度,r0和a0分别为监测空域的径向距离和角度的起始值,Δr和Δa分别为距离和角度分辨单元的大小,Lmax为目标的先验最大可能距离扩展长度;
3.3)将最新一帧雷达图像观测的目标信息加入TIM中,完成滑窗操作:Ak+1={Aτ|τ=k-Wl+1,…,k}。
作为本发明的进一步改进,所述步骤4)通过以下步骤完成:
对每一个伯努利成分计算目标存在概率以及相应的密度分布函数为其中 gz为基于雷达图像数据的多目标似然函数, 其中g1(·|xk)和g0(·)分别为有目标时和无目标时像素分辨单元的似然函数,lk(xk)为目标的似然比函数,为监测空域仅存在噪声时的似然函数。
作为本发明的进一步改进,所述步骤5)通过以下步骤完成:
5.1)设置伯努利成分修剪门限T和允许的最大伯努利成分数Jmax,将多伯努利RFS中存活概率小于修剪门限的成分删除;若伯努利成分数大于Jmax,则只保留存活概率最大的Jmax个成分。
5.3)设置合并门限U,将多伯努利RFS中的伯努利RFS按照存活概率由大到小进行排序,形成一个成分列表。取最大存活概率的成分作为聚类中心,与剩余的伯努利RFS成分与中心的马氏距离是否小于设置的门限;对小于合并门限U的成分进行合并,假设第i个成分和第j个伯努利成分需要进行合并,合并后的伯努利成分对应的存活概率为合并的成分从列表中删去,重复选取列表中存活概率最大的成分作为聚类中心,直至所有伯努利成分都完成修剪和合并。
与现有技术相比,本发明所具有的有益效果为:
本发明提供了一种面向多距离扩展目标的多伯努利检测前跟踪方法,针对扩展目标在检测前跟踪框架下需要幅度先验信息的问题,通过积累目标潜在航迹的历史幅度观测数据,利用主成分分析的方法提取目标幅度分布信息,并结合多伯努利滤波器在时变数目的多目标跟踪上的优势,实现多个距离扩展非合作目标的检测前跟踪。本发明的技术效果如下:1)本发明方法利用TBD技术能够发挥联合检测与跟踪的优势,结合有限集统计学理论对时变数目目标的建模,在雷达视频观测数据上进行滤波,实现对多个距离扩展目标的检测前跟踪;2)本发明方法通过粒子对多帧测量的积累和主成分分析工具来近似估计目标幅度分布信息,能够有效提取观测数据中的目标信息;同时,该估计方法对测量数据满足的分布类型没有特定的要求,具备通用性。
附图说明
图1为扩展目标示意图;
图2为多帧滑窗缓存示意图;
图3为本发明实施例流程示意图;
图4为本发明的流程图
具体实施方式
下面结合附图和具体实施例对本发明进一步详细说明。
参照附图3,本发明提出一种面向多距离扩展目标的多伯努利检测前跟踪方法,本发明涉及非线性状态估计问题,提供一种粒子滤波近似实现方案,具体通过以下步骤实现:
步骤1,雷达系统处理接收回波信号,得到距离-方位回波幅度测量图像。本实施例中,
1.1)雷达系统发射大时宽大带宽信号,接收回波通过距离脉冲压缩算法得到高分辨距离像;
1.2)雷达系统对监测区域的不同角度进行扫描得到不同方位角度的高分辨距离像,由此,一次扫描得到一幅包含Nr×Na个分辨单元的雷达距离-方位回波幅度测量图像,其中Nr为距离分辨单元数,Na为方位角度分辨单元数。
步骤2,根据目标运动模型,对多目标状态密度分布进行预测。本实施例中,
2.1)选取距离扩展目标上离雷达最远端的点作为目标的参考点,并假设刚体目标占据至多Lmax个距离分辨单元。根据先验知识给定目标的运动模型,本实施例中,假设目标仅进行近似匀速直线运动,则对应的目标的运动模型为xk=Fxk-1+wk,式中是k时刻目标的状态矢量,其中(x,y)′和分别是目标的位置和速度,(·)′表示矢量的转置,F为目标状态转移矩阵,wk是服从pw分布的过程噪声,假设为高斯白噪声,并记目标靠近雷达系统为正方向。将k-1时刻的多目标状态表示为Xk-1,其中的每一个目标对应的状态表示为xk-1,则对应的后验多目标密度为其中和分别为目标的存活概率和概率密度分布,Mk-1为k-1时刻的目标数,在本实施例中,采用粒子近似方法,因此目标后验密度可以由带权重的样本粒子表示,即 其中为第i个伯努利成分的第j个粒子对应的TIM。k时刻,目标将pS,k(xk-1)的概率转移到xk,并且状态转移概率密度为ζk|k-1(xk|xk-1),或者以1-pS,k(xk-1)的概率从场景中消失。由此将预测过程分为两部分:新生(new birth)部分和持续存活(survival,也称为幸存部分)部分;
在本实施例中,采用粒子近似实现,对于持续存活目标部分,在k-1时刻存在的目标,由目标运动学模型可以得到单目标的状态转移概率密度函数根据粒子k-1时刻的状态可以预测k时刻的粒子状态:当将目标的状态转移概率密度函数作为粒子的提议分布密度时,粒子权重为:归一化后得到: 得到持续存活目标部分的伯努利RFS预测参数预测存在概率为:以及对应的预测密度为:
在本实施例中,采用粒子近似实现,对于新生目标,在k时刻自发新生的目标部分,通过利用先验知识建立的新生密度分布产生相应的粒子: 当将目标的新生密度分布作为粒子的提议分布密度时,粒子对应的权重为:归一化后得到:由此,得到新生目标部分的伯努利RFS预测参数 由新生模型先验给定,新生目标概率分布密度为:
2.4)以持续存活目标的伯努利RFSSk|k-1(xk-1)和新生目标的伯努利RFS Γk的并集来描述k时刻的多目标预测状态其中多伯努利RFS中的伯努利RFS成分满足相互独立假设;多伯努利贝叶斯预测过程通过预测成分和新生成分的并来表示:Mk|k-1=Mk-1+MΓ,k是预测多伯努利RFS的伯努利成分数。
在本实施例中,通过上述2.2和2.3两个步骤得到k时刻预测的多伯努利RFS为 式中Mk|k-1=Mk-1+MΓ,k是预测多伯努利RFS的伯努利成分数,且每一个伯努利成分的参数可以由一群带权重的粒子进行近似:i=1,…,Mk|k-1。
步骤3,通过时域滑窗对距离扩展目标的历史测量数据进行储存,构造目标信息矩阵,并借助主成分分析工具对其幅度分布信息进行估计。
3.1)目标TIM通过一个滑动窗口,积累连续多帧图像测量数据中距离扩展目标所影响的图像区域分辨单元的强度数据,其中目标的图像影响范围由T(xk)描述,xk为目标状态,则多帧滑窗缓存获取的k时刻TIM的数学描述符号为Ak={Aτ|τ=k-Wl,…,k-1},其中 是k时刻第tk个分辨单元的强度测量数据,τ为时间索引,Wl是时间滑窗的长度;针对状态为xk的距离扩展目标,其图像影响范围T(xk)由下式进行描述: 其中rk(xk)和ak(xk)分别为目标相对雷达的径向长度和角度,r0和a0分别为监测空域的径向距离和角度的起始值,Δr和Δa分别为距离和角度分辨单元的大小,Lmax为目标的先验最大可能距离扩展长度;
在本实施例中,为了使粒子输出有效的幅度估计值,每个粒子须首先判断TIM中是否已经累计了Wl个历史帧观测数据,若达到Wl帧,则计算k时刻TIM的协方差矩阵Sk,求解相应的特征方程:SkG=-λmaxG,得到最大特征值λmax及其对应的特征向量G,代入 得到目标的幅度分布估计否则式中0为零向量。
3.3)将最新一帧雷达图像观测的目标信息加入TIM中,完成滑窗操作:Ak+1={Aτ|τ=k-Wl+1,…,k};
在本实施例中,每个粒子将估计得到的输出到校正器后,滑窗后移,将当前帧的数据积累到粒子的TIM中,并移除最Wl帧以前的测量数据,即将当前粒子预测状态影响的测量数据加入到中,该过程可描述为: 可通过数据结构中的队列结构完成滑窗步骤。
步骤4,根据目标幅度估计结果和测量图像,结合预测密度对多目标状态进行更新,得到多目标后验密度。
对每一个伯努利成分计算目标存在概率以及相应的密度分布函数为其中 gz为基于雷达图像数据的多目标似然函数, 其中g1(·|xk)和g0(·)分别为有目标时和无目标时像素分辨单元的似然函数,lk(xk)为目标的似然比函数,为监测空域仅存在噪声时的似然函数。
本实施例中,采用粒子近似方法,根据贝叶斯基本原理和多伯努利更新解析表达式可以得到各个伯努利成分的目标存在概率以及相应的目标状态密度分布函数其中为利用雷达图像观测数据计算得到的似然值对粒子权重进行更新后的结果,且 为似然比函数。
步骤5,设置最大伯努利成份数以及合并和修剪阈值,对多目标航迹进行管理,实现多目标检测前跟踪。本实施例中,
5.1)根据系统性能需求设置伯努利成分修剪门限T和允许的最大伯努利成分数Jmax,对多伯努利RFS中逐个伯努利成分进行遍历,将多伯努利RFS中存活概率小于修剪门限的成分删除;若伯努利成分数大于Jmax,则只保留存活概率最大的Jmax个成分。
5.3)设置合并门限U,将多伯努利RFS中的伯努利RFS按照存活概率由大到小进行排序,形成一个成分列表。取最大存活概率的成分作为聚类中心,与剩余的伯努利RFS成分与中心的马氏距离是否小于设置的门限;对小于合并门限U的成分进行合并,假设第i个成分和第j个伯努利成分需要进行合并,合并后的伯努利成分对应的存活概率为 合并的成分从列表中删去,重复选取列表中存活概率最大的成分作为聚类中心,直至所有伯努利成分都完成修剪和合并。
以上仅是本发明的实施方式之一,本发明的保护范围并不仅限于上述实例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应当视为落入本发明的保护范围。
Claims (6)
1.一种面向多距离扩展目标的多伯努利检测前跟踪方法,其特征在于,包含以下步骤:
步骤1,雷达系统处理接收回波信号,得到距离-方位回波幅度测量图像;
步骤2,根据目标运动模型,对多目标状态密度分布进行预测;
步骤3,通过时域滑窗对距离扩展目标的历史测量数据进行储存,构造目标信息矩阵,并借助主成分分析工具对其幅度分布信息进行估计;
步骤4,根据目标幅度估计结果和测量图像,结合预测密度对多目标状态进行更新,得到多目标后验密度;
步骤5,设置最大伯努利成份数以及合并和修剪阈值,对多目标航迹进行管理,实现多目标检测前跟踪。
2.如权利要求1所述的一种面向多距离扩展目标的多伯努利检测前跟踪方法,其特征在于,所述步骤1中,通过以下步骤完成:
1.1)雷达系统发射大时宽大带宽信号,通过距离脉冲压缩得到高分辨距离像;
1.2)雷达系统对监测区域的不同角度进行扫描得到不同方位角度的高分辨距离像,一次扫描得到一幅包含Nr×Na个分辨单元的雷达距离-方位回波幅度测量图像,其中Nr为距离分辨单元数,Na为方位角度分辨单元数。
3.如权利要求1所述的一种面向多距离扩展目标的多伯努利检测前跟踪方法,其特征在于,所述步骤2中,通过以下步骤完成:
2.1)选取目标的远端,即目标上距离雷达最远端的点作为目标的参考点,并根据先验知识给定系统模型,以下给定先验知识:将k-1时刻的多目标状态表示为Xk-1,其中的每一个目标对应的状态表示为xk-1,则对应的后验多目标密度为其中和分别为目标的存活概率和概率密度分布,Mk-1为k-1时刻的伯努利成分数;k时刻,目标将以pS,k(xk-1)的概率转移到xk,并且状态转移概率密度为ζk|k-1(xk|xk-1),或者以1-pS,k(xk-1)的概率从场景中消失,相应的数学模型表示为[pS,k(xk-1),ζk|k-1(xk|xk-1){的伯努利RFSSk|k-1(xk-1),由此将预测过程分为两部分:新生部分和持续存活部分;
4.如权利要求1所述的一种面向多距离扩展目标的多伯努利检测前跟踪方法,其特征在于,所述步骤3中,通过以下步骤完成:
3.1)目标TIM通过一个滑动窗口,积累连续多帧图像测量数据中距离扩展目标所影响的图像区域分辨单元的强度数据,其中目标的图像影响范围由T(xk)描述,xk为目标状态,则多帧滑窗缓存获取的k时刻TIM的数学描述符号为Ak={Aτ|τ=k-Wl,…,k-1},其中 是k时刻第tk个分辨单元的强度测量数据,τ为时间索引,Wl是时间滑窗的长度;针对状态为xk的距离扩展目标,其图像影响范围T(xk)由下式进行描述: 其中rk(xk)和ak(xk)分别为目标相对雷达的径向长度和角度,r0和a0分别为监测空域的径向距离和角度的起始值,Δr和Δa分别为距离和角度分辨单元的大小,Lmax为目标的先验最大可能距离扩展长度;
3.3)将最新一帧雷达图像观测的目标信息加入TIM中,完成滑窗操作:Ak+1={Aτ|τ=k-Wl+1,…,k}。
6.如权利要求1所述的一种面向多距离扩展目标的多伯努利检测前跟踪方法,其特征在于,所述步骤5中,通过以下步骤完成:
5.1)设置伯努利成分修剪门限T和允许的最大伯努利成分数Jmax,将多伯努利RFS中存活概率小于修剪门限的成分删除;若伯努利成分数大于Jmax,则只保留存活概率最大的Jmax个成分;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010403701.5A CN111580084B (zh) | 2020-05-13 | 2020-05-13 | 一种面向多距离扩展目标的多伯努利检测前跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010403701.5A CN111580084B (zh) | 2020-05-13 | 2020-05-13 | 一种面向多距离扩展目标的多伯努利检测前跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111580084A true CN111580084A (zh) | 2020-08-25 |
CN111580084B CN111580084B (zh) | 2021-02-02 |
Family
ID=72124941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010403701.5A Active CN111580084B (zh) | 2020-05-13 | 2020-05-13 | 一种面向多距离扩展目标的多伯努利检测前跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111580084B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112285698A (zh) * | 2020-12-25 | 2021-01-29 | 四川写正智能科技有限公司 | 一种基于雷达传感器的多目标跟踪装置及方法 |
CN112415468A (zh) * | 2020-10-27 | 2021-02-26 | 桂林电子科技大学 | 一种基于多伯努利滤波的doa跟踪方法 |
CN113048884A (zh) * | 2021-03-17 | 2021-06-29 | 西安工业大学 | 一种扩展目标跟踪实验平台及其试验方法 |
CN114061584A (zh) * | 2021-11-02 | 2022-02-18 | 江苏科技大学 | 一种基于多机器人的势均衡多伯努利滤波slam方法 |
CN114325686A (zh) * | 2021-12-23 | 2022-04-12 | 中国人民解放军国防科技大学 | 基于smc-phd滤波器的多目标跟踪方法 |
CN114609634A (zh) * | 2022-03-21 | 2022-06-10 | 电子科技大学 | 一种基于阴影的交互式多模型下的视频sar多目标跟踪方法 |
CN114690175A (zh) * | 2022-04-12 | 2022-07-01 | 电子科技大学 | 一种基于被动外辐射源雷达的目标直接检测与跟踪方法 |
CN115097437A (zh) * | 2022-06-06 | 2022-09-23 | 哈尔滨工程大学 | 一种基于标签多伯努利检测前跟踪算法的水下目标跟踪轨迹临近交叉解决方法 |
CN114061584B (zh) * | 2021-11-02 | 2024-05-24 | 江苏科技大学 | 一种基于多机器人的势均衡多伯努利滤波slam方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102043150A (zh) * | 2010-12-06 | 2011-05-04 | 电子科技大学 | 用于微弱目标检测的改进型粒子滤波检测前跟踪方法 |
CN104597434A (zh) * | 2015-01-14 | 2015-05-06 | 中国民航大学 | 改进包络移位补偿和分数阶傅里叶变换的多帧相参tbd方法 |
CN104597435A (zh) * | 2015-01-14 | 2015-05-06 | 中国民航大学 | 基于修正频域补偿和分数阶傅里叶变换的多帧相参tbd方法 |
CN104766320A (zh) * | 2015-04-02 | 2015-07-08 | 西安电子科技大学 | 阈值化量测下的多伯努利滤波弱目标检测与跟踪方法 |
JP2015180858A (ja) * | 2014-03-03 | 2015-10-15 | 三菱電機株式会社 | レーダ装置 |
US20180005566A1 (en) * | 2016-06-29 | 2018-01-04 | Semiconductor Energy Laboratory Co., Ltd. | Electronic device, operation method of the electronic device, and moving vehicle |
CN107703496A (zh) * | 2017-10-12 | 2018-02-16 | 桂林电子科技大学 | 一种交互式多模伯努利滤波的机动弱目标检测前跟踪方法 |
CN109901154A (zh) * | 2019-03-29 | 2019-06-18 | 中国人民解放军海军航空大学 | 基于递推rtht-tbd的自适应调整方法 |
CN109991597A (zh) * | 2019-04-04 | 2019-07-09 | 中国人民解放军国防科技大学 | 一种面向弱扩展目标的检测前跟踪方法 |
-
2020
- 2020-05-13 CN CN202010403701.5A patent/CN111580084B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102043150A (zh) * | 2010-12-06 | 2011-05-04 | 电子科技大学 | 用于微弱目标检测的改进型粒子滤波检测前跟踪方法 |
JP2015180858A (ja) * | 2014-03-03 | 2015-10-15 | 三菱電機株式会社 | レーダ装置 |
CN104597434A (zh) * | 2015-01-14 | 2015-05-06 | 中国民航大学 | 改进包络移位补偿和分数阶傅里叶变换的多帧相参tbd方法 |
CN104597435A (zh) * | 2015-01-14 | 2015-05-06 | 中国民航大学 | 基于修正频域补偿和分数阶傅里叶变换的多帧相参tbd方法 |
CN104766320A (zh) * | 2015-04-02 | 2015-07-08 | 西安电子科技大学 | 阈值化量测下的多伯努利滤波弱目标检测与跟踪方法 |
US20180005566A1 (en) * | 2016-06-29 | 2018-01-04 | Semiconductor Energy Laboratory Co., Ltd. | Electronic device, operation method of the electronic device, and moving vehicle |
CN107703496A (zh) * | 2017-10-12 | 2018-02-16 | 桂林电子科技大学 | 一种交互式多模伯努利滤波的机动弱目标检测前跟踪方法 |
CN109901154A (zh) * | 2019-03-29 | 2019-06-18 | 中国人民解放军海军航空大学 | 基于递推rtht-tbd的自适应调整方法 |
CN109991597A (zh) * | 2019-04-04 | 2019-07-09 | 中国人民解放军国防科技大学 | 一种面向弱扩展目标的检测前跟踪方法 |
Non-Patent Citations (3)
Title |
---|
RONGHUI ZHAN等: ""Maneuvering targets track-before-detect using multiple-model multi-bernoulli filtering"", 《2013 INTERNATIONAL CONFERENCE ON INFORMATION TECHNOLOGY AND APPLICATION》 * |
THARINDU R等: ""Labeled multi-bernoulli track-before-detect for multi-target tracking in video"", 《18TH INTERNATIONAL CONFERENCE ON INFORMATION FUSION》 * |
杨威等: ""弱目标检测前跟踪技术研究综述"", 《电子学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112415468A (zh) * | 2020-10-27 | 2021-02-26 | 桂林电子科技大学 | 一种基于多伯努利滤波的doa跟踪方法 |
CN112285698A (zh) * | 2020-12-25 | 2021-01-29 | 四川写正智能科技有限公司 | 一种基于雷达传感器的多目标跟踪装置及方法 |
CN113048884B (zh) * | 2021-03-17 | 2022-12-27 | 西安工业大学 | 一种扩展目标跟踪实验平台及其试验方法 |
CN113048884A (zh) * | 2021-03-17 | 2021-06-29 | 西安工业大学 | 一种扩展目标跟踪实验平台及其试验方法 |
CN114061584A (zh) * | 2021-11-02 | 2022-02-18 | 江苏科技大学 | 一种基于多机器人的势均衡多伯努利滤波slam方法 |
CN114061584B (zh) * | 2021-11-02 | 2024-05-24 | 江苏科技大学 | 一种基于多机器人的势均衡多伯努利滤波slam方法 |
CN114325686A (zh) * | 2021-12-23 | 2022-04-12 | 中国人民解放军国防科技大学 | 基于smc-phd滤波器的多目标跟踪方法 |
CN114325686B (zh) * | 2021-12-23 | 2024-05-03 | 中国人民解放军国防科技大学 | 基于smc-phd滤波器的多目标跟踪方法 |
CN114609634B (zh) * | 2022-03-21 | 2023-04-25 | 电子科技大学 | 一种基于阴影的交互式多模型下的视频sar多目标跟踪方法 |
CN114609634A (zh) * | 2022-03-21 | 2022-06-10 | 电子科技大学 | 一种基于阴影的交互式多模型下的视频sar多目标跟踪方法 |
CN114690175B (zh) * | 2022-04-12 | 2023-10-27 | 电子科技大学 | 一种基于被动外辐射源雷达的目标直接检测与跟踪方法 |
CN114690175A (zh) * | 2022-04-12 | 2022-07-01 | 电子科技大学 | 一种基于被动外辐射源雷达的目标直接检测与跟踪方法 |
CN115097437A (zh) * | 2022-06-06 | 2022-09-23 | 哈尔滨工程大学 | 一种基于标签多伯努利检测前跟踪算法的水下目标跟踪轨迹临近交叉解决方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111580084B (zh) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111580084B (zh) | 一种面向多距离扩展目标的多伯努利检测前跟踪方法 | |
Panta et al. | Probability hypothesis density filter versus multiple hypothesis tracking | |
US8405540B2 (en) | Method for detecting small targets in radar images using needle based hypotheses verification | |
CN109886305B (zh) | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 | |
Clark et al. | GM-PHD filter multitarget tracking in sonar images | |
CN108107429B (zh) | 基于最大似然估计的前视超分辨成像方法 | |
CN111999735B (zh) | 一种基于径向速度和目标跟踪的动静目标分离方法 | |
Vo et al. | An implementation of the multi-sensor generalized labeled multi-Bernoulli filter via Gibbs sampling | |
Grossi et al. | A heuristic algorithm for track-before-detect with thresholded observations in radar systems | |
CN108010066B (zh) | 基于红外目标灰度互相关和角度信息的多假设跟踪方法 | |
CN108196238B (zh) | 高斯背景下基于自适应匹配滤波的杂波图检测方法 | |
CN113960587A (zh) | 基于类别信息反馈的毫米波雷达多目标跟踪方法 | |
Willett et al. | MLPDA and MLPMHT applied to some MSTWG data | |
Yi et al. | Multi-target tracking via dynamic-programming based track-before-detect | |
Cong et al. | Bayesian track-before-detect algorithm for nonstationary sea clutter | |
CN111007880B (zh) | 一种基于汽车雷达的扩展目标跟踪方法 | |
Horridge et al. | Using a probabilistic hypothesis density filter to confirm tracks in a multi-target environment. | |
CN108387879B (zh) | 基于自适应归一化匹配滤波的杂波图单元中值检测方法 | |
Coraluppi et al. | Recent advances in multi-INT track fusion | |
Nandakumaran et al. | Improved multitarget tracking using probability hypothesis density smoothing | |
CN113126086B (zh) | 一种基于状态预测积累的生命探测雷达弱目标检测方法 | |
Vu et al. | Track-before-detect for an active towed array sonar | |
Mansouri et al. | Channel estimation and multiple target tracking in wireless sensor networks based on quantised proximity sensors | |
Coraluppi | Seeing the forest through the trees: Multiple-hypothesis tracking and graph-based tracking extensions | |
Zhu et al. | Multitarget tracking based on PHD smoother with unknown clutter spatial density |
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 |