CN107783104A - 一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 - Google Patents
一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 Download PDFInfo
- Publication number
- CN107783104A CN107783104A CN201710966437.4A CN201710966437A CN107783104A CN 107783104 A CN107783104 A CN 107783104A CN 201710966437 A CN201710966437 A CN 201710966437A CN 107783104 A CN107783104 A CN 107783104A
- Authority
- CN
- China
- Prior art keywords
- particle
- target
- sensor
- tracking
- state
- 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.)
- Pending
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
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- 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/87—Combinations of radar systems, e.g. primary radar and secondary radar
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
技术背景
本发明属于多传感器单目标跟踪领域,涉及一种基于粒子滤波的多异步传感器检测前跟踪方法,用于解决弱目标检测跟踪情况的多异步传感器信息融合,达到提高目标跟踪精度和目标正确发现概率。
背景技术
在低信噪比的复杂环境下,雷达探测系统对弱目标的检测跟踪问题一直都是一项技术难题。检测前跟踪(Tracking Before Detection,TBD)是一种可以在低信噪比环境下提高传感器检测概率的方法,它在进行检测的过程中加入跟踪的思想,并且降低对目标的检测门限,然后再通过对目标原始能量的积累来进行目标检测,常用于对微弱目标信号的检测。目前的TBD方法主要有Hough变换、动态规划(Dynamic Programming,DP)及粒子滤波(Particle Filter,PF)。
现在已有的关于粒子滤波检测前跟踪的算法主要研究的问题是关于单传感器方面的,在多雷达目标检测跟踪方面的研究非常少。但是在实际情况下,一般都是通过多个传感器对目标进行联合探测,然后根据各个传感器之间的探测信息的互补性来提高目标正确发现概率和跟踪效果。但是当融合多个传感器的信息的时候,因为雷达采样时刻不同,多部传感器对相同目标的量测并不是目标处于同一位置的量测,如果直接将这些量测信息进行融合,最后得到的结果会有很大的误差。因此提出一种基于粒子滤波的多异步传感器检测前跟踪算法,提高目标跟踪精度和目标正确发现概率。
发明内容
本发明针对弱目标低信噪比环境下,多异步传感器将获得的不同时刻目标量测信息进行融合问题,本发明提出了一种基于粒子滤波的多异步传感器单目标检测前跟踪方法。
为了实现上述目的,本发明的技术方案的步骤为:
步骤1:设置粒子集初始分布;
步骤2:对粒子集的权重和粒子数进行初始化;
步骤3:构建多传感器单目标跟踪场景,并对目标的运动模型进行初始化,设置目标运动的相关参数,包括目标运动的过程噪声和传感器的量测噪声;
步骤4:粒子状态转移。将粒子的出现概率跟马尔科夫概率转移设定参数比较判断是否状态转移。
步骤5:权值计算。将存在的粒子根据方位和距离信息计算出对应的权重值,通过内插外推的空时校准方法将各传感器同一时刻的权重值进行数据融合;
其中通过内插外推的空时校准方法将各传感器同一时刻的权重值进行数据融合;具体包括如下步骤:
将其中一个雷达作为基准雷达,将其他雷达量测到的某一时刻的粒子的状态信息通过内插外推法,将该时刻的粒子状态进行时间和空间的回推,回推至基准雷达相邻的采样时刻,然后通过将同一粒子的权重相融合来实现对目标的检测跟踪获得目标航迹。
步骤6:将融合权值进行归一化。
步骤7:计算有效样本数Neff,若Neff<NT,NT为阈值,则进行重采样,之后将粒子集的各个权重值重新设置为初始值,跳至步骤8,若Neff>NT,直接跳至步骤8;
步骤8:估计目标状态;将目标估计状态返回至每个传感器,用于下一时刻的传感器量测估计,重复步骤4-8。
本发明的益处是:针对低信噪比弱目标环境下多异步传感器单目标检测前跟踪问题,本发明提出了一套完整的处理方法流程,最主要的是提出了一种多异步雷达目标量测信息融合方法。本发明可以有效的提高目标正确发现概率和跟踪精度。
具体实施方式
一种基于粒子滤波的多异步传感器单目标检测前跟踪方法,该方法具体包括以下步骤:
步骤1:设置粒子集初始分布p(s0,r0)。
步骤2:粒子数初始化取为M,根据粒子集的初始分布s0 j~p(s0,r0),j=1,2,...,M,得到粒子集并且将i=1,...,M各粒子的初始化权重设置为1/M。
步骤3:构建多传感器单目标跟踪场景,传感器的量测可能来自目标,也可能来自杂波,构建目标的运动模型、传感器的量测模型和杂波模型,并对它们进行初始化。
建立目标的运动模型:Sk+1=fk(tk,sk,ck,wk)
Prob(ck+1=i|ck=j)=[Ψ(tk)]ij
式中,表示在k时刻的目标状态,其中(xk,yk)表示目标位置,表示目标速度;tk表示k时刻的采样时间;ck表示模型变量,通过这个模型变量来表示目标是否出现;Ψ(tk)是马尔可夫状态转移矩阵。wk是高斯过程噪声,prob(w)是分布函数。模型之间的状态转移概率通过Ψ(tk)来完成,
传感器量测建模:设有L个雷达同时检测某一个区域,这些雷达的采样周期相同但是系统误差各不相同。假设第i个雷达在k时刻产生一系列距离—方位图像,其中每帧图像包含个分辨单元,将每个分辨单元设为(m,n),m=1,2,...,Nr,n=1,2,...,Nb对应一个矩形区域。则雷达i在时刻k的观测数据可以表示为:其中表示分辨单元(m,n)在时刻k的观测,其具体形式为:
式中,用vk,i来表示雷达i的按高斯分布的观测噪声。用p(vi)来表示它的分布函数。用Δsi来表示传感器i的测量系统误差,用来表示雷达i的第(m,n)分辨单元所接受到的量测信息,其函数形式可表示为:其中,Ri是与距离有关的常数,Bi是与方位有关的常数,和分为为衰减常数,rk,i是目标相对应的雷达的距离信息,bk,i是相对应的雷达的方位信息。设雷达i位置坐标为(xpi,ypi),则 噪声采用复噪声形式,即:式中,和是均值为零,方差为的高斯白噪声,并且相互独立。信噪比为 k传感器i获取的观测可以表示为:Zk,i={z1,i,z2,i,...,zn,i}。式中用Zk,1:L来表示在k时刻全部的雷达观测数据。
步骤4:粒子状态转移。按照噪声分布函数p(w)获取噪声样本集并依据状态转移方程式预测得到时刻k的基于状态装异的粒子集
步骤5:权值计算。由于多个雷达获得的是目标在不同位置的回波信息,这样直接计算得到的粒子j的各个权重是不能直接相乘融合的,所以需要以其中一个雷达的时刻作为基准雷达。假设以雷达1的时刻作为基准时刻,然后将其余雷达i=2,...,N对粒子j相对应的时刻的状态通过内插外推法分别对x坐标和y坐标进行时间回推。以匀加速直线运动x轴方向回推为例,当雷达i的采样时刻tik(k=1...n)回推到雷达1的采样时刻t1c(c=1...m):
目标在tik时刻的加速度为:
目标在tik时刻的速度为:
设配准时刻为t1c(tik<t1c<ti(k+1)),将雷达i向基准时刻配准后:
假设目标量测信息满足指数分布,则分布函数为其中当目标存在和不存在时均值分别为则当目标存在和不存在时权重的计算公式分别为为其中表示回推到基准时刻的状态然后将粒子j对应不同雷达i的权重相乘得到融合后的权重
步骤6:将融合权值进行归一化。由于各个雷达的测量精度和参数并不相同,所以需要将通过空时校准得到的融合的j粒子做归一化处理:得到
步骤7:计算有效样本数NT为阈值,若Neff<NT,则对粒子集进行重采样得到之后将粒子集的各个权重值重新设置为初始值。
步骤8:估计目标状态。估计目标状态将目标估计状态返回至每个传感器,用于下一时刻的传感器量测估计返回步骤4进行重复。
Claims (1)
1.一种基于粒子滤波的多异步传感器单目标检测前跟踪方法,其特征在于,该方法具体包括以下步骤:
步骤1:设置粒子集初始分布;
步骤2:对粒子集的权重和粒子数进行初始化;
步骤3:构建多传感器单目标跟踪场景,并对目标的运动模型进行初始化,设置目标运动的相关参数,包括目标运动的过程噪声和传感器的量测噪声;
步骤4:粒子状态转移;将粒子的出现概率跟马尔科夫概率转移设定参数比较判断是否状态转移;
步骤5:权值计算;将存在的粒子根据方位和距离信息计算出对应的权重值,通过内插外推的空时校准方法将各传感器同一时刻的权重值进行数据融合;
其中通过内插外推的空时校准方法将各传感器同一时刻的权重值进行数据融合;具体包括如下步骤:
将其中一个雷达作为基准雷达,将其他雷达量测到的某一时刻的粒子的状态信息通过内插外推法,将该时刻的粒子状态进行时间和空间的回推,回推至基准雷达相邻的采样时刻,然后通过将同一粒子的权重相融合来实现对目标的检测跟踪获得目标航迹;
步骤6:将融合权值进行归一化;
步骤7:计算有效样本数Neff,若Neff<NT,NT为阈值,则进行重采样,之后将粒子集的各个权重值重新设置为初始值,跳至步骤8,若Neff>NT,直接跳至步骤8;
步骤8:估计目标状态;将目标估计状态返回至每个传感器,用于下一时刻的传感器量测估计,重复步骤4-8。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710966437.4A CN107783104A (zh) | 2017-10-17 | 2017-10-17 | 一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710966437.4A CN107783104A (zh) | 2017-10-17 | 2017-10-17 | 一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107783104A true CN107783104A (zh) | 2018-03-09 |
Family
ID=61433725
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710966437.4A Pending CN107783104A (zh) | 2017-10-17 | 2017-10-17 | 一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107783104A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188424A (zh) * | 2018-09-14 | 2019-01-11 | 中国人民解放军海军航空大学 | 基于量测一致性的分布式多传感器多目标跟踪方法 |
CN109886305A (zh) * | 2019-01-23 | 2019-06-14 | 浙江大学 | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 |
CN110058222A (zh) * | 2019-03-29 | 2019-07-26 | 杭州电子科技大学 | 一种基于传感器选择的双层粒子滤波检测前跟踪方法 |
CN110187335A (zh) * | 2019-06-25 | 2019-08-30 | 电子科技大学 | 针对具有非连续特性目标的粒子滤波检测前跟踪方法 |
CN110286383A (zh) * | 2019-06-27 | 2019-09-27 | 电子科技大学 | 一种应用于目标跟踪的雷达和红外传感器部署方法 |
CN110376556A (zh) * | 2019-06-11 | 2019-10-25 | 杭州电子科技大学 | 一种基于锦标赛选择的双层粒子滤波检测前跟踪方法 |
CN111812636A (zh) * | 2020-06-01 | 2020-10-23 | 杭州电子科技大学 | 一种基于权重融合选择的粒子滤波检测前跟踪方法 |
CN111856415A (zh) * | 2020-06-19 | 2020-10-30 | 中国电子科技集团公司第二十七研究所 | 一种雷达数据处理设备超前校准方法、装置和存储介质 |
CN112415503A (zh) * | 2020-10-15 | 2021-02-26 | 杭州电子科技大学 | 一种基于目标重跟踪的多目标粒子滤波检测前跟踪方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106842165A (zh) * | 2017-03-16 | 2017-06-13 | 电子科技大学 | 一种基于不同距离角度分辨率雷达集中式异步融合方法 |
-
2017
- 2017-10-17 CN CN201710966437.4A patent/CN107783104A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106842165A (zh) * | 2017-03-16 | 2017-06-13 | 电子科技大学 | 一种基于不同距离角度分辨率雷达集中式异步融合方法 |
Non-Patent Citations (3)
Title |
---|
权宏伟 等: "低空目标检测融合系统的结构与框架设计", 《火力与指挥控制》 * |
程小会 等: "基于误差校准的多雷达粒子滤波检测跟踪算法", 《火力与指挥控制》 * |
程小会: "基于粒子滤波的多雷达单目标检测前跟踪算法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188424A (zh) * | 2018-09-14 | 2019-01-11 | 中国人民解放军海军航空大学 | 基于量测一致性的分布式多传感器多目标跟踪方法 |
CN109188424B (zh) * | 2018-09-14 | 2020-09-04 | 中国人民解放军海军航空大学 | 基于量测一致性的分布式多传感器多目标跟踪方法 |
CN109886305A (zh) * | 2019-01-23 | 2019-06-14 | 浙江大学 | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 |
CN109886305B (zh) * | 2019-01-23 | 2021-05-04 | 浙江大学 | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 |
CN110058222B (zh) * | 2019-03-29 | 2020-11-24 | 杭州电子科技大学 | 一种基于传感器选择的双层粒子滤波检测前跟踪方法 |
CN110058222A (zh) * | 2019-03-29 | 2019-07-26 | 杭州电子科技大学 | 一种基于传感器选择的双层粒子滤波检测前跟踪方法 |
CN110376556A (zh) * | 2019-06-11 | 2019-10-25 | 杭州电子科技大学 | 一种基于锦标赛选择的双层粒子滤波检测前跟踪方法 |
CN110376556B (zh) * | 2019-06-11 | 2021-05-11 | 杭州电子科技大学 | 一种基于锦标赛选择的双层粒子滤波检测前跟踪方法 |
CN110187335A (zh) * | 2019-06-25 | 2019-08-30 | 电子科技大学 | 针对具有非连续特性目标的粒子滤波检测前跟踪方法 |
CN110286383A (zh) * | 2019-06-27 | 2019-09-27 | 电子科技大学 | 一种应用于目标跟踪的雷达和红外传感器部署方法 |
CN110286383B (zh) * | 2019-06-27 | 2023-09-01 | 电子科技大学 | 一种应用于目标跟踪的雷达和红外传感器部署方法 |
CN111812636A (zh) * | 2020-06-01 | 2020-10-23 | 杭州电子科技大学 | 一种基于权重融合选择的粒子滤波检测前跟踪方法 |
CN111812636B (zh) * | 2020-06-01 | 2023-06-13 | 杭州电子科技大学 | 一种基于权重融合选择的粒子滤波检测前跟踪方法 |
CN111856415A (zh) * | 2020-06-19 | 2020-10-30 | 中国电子科技集团公司第二十七研究所 | 一种雷达数据处理设备超前校准方法、装置和存储介质 |
CN112415503A (zh) * | 2020-10-15 | 2021-02-26 | 杭州电子科技大学 | 一种基于目标重跟踪的多目标粒子滤波检测前跟踪方法 |
CN112415503B (zh) * | 2020-10-15 | 2023-10-20 | 杭州电子科技大学 | 一种基于目标重跟踪的多目标粒子滤波检测前跟踪方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107783104A (zh) | 一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 | |
CN103729859B (zh) | 一种基于模糊聚类的概率最近邻域多目标跟踪方法 | |
CN101509972B (zh) | 基于高分辨目标距离像修正相关矩阵的宽带雷达检测方法 | |
CN102914776B (zh) | 基于模糊c均值算法的多通道SAR动目标定位方法 | |
CN108333569A (zh) | 一种基于phd滤波的异步多传感器融合多目标跟踪方法 | |
CN105487060B (zh) | 一种双通道四斜率调制的多目标提取方法 | |
CN103954939B (zh) | 基于雷达组网的抗灵巧噪声干扰方法 | |
KR101628154B1 (ko) | 수신 신호 세기를 이용한 다중 표적 추적 방법 | |
CN113267773B (zh) | 一种基于毫米波雷达的室内人员准确检测和精确定位方法 | |
CN107576959B (zh) | 一种基于区域映射解模糊的高重频雷达目标检测前跟踪方法 | |
CN105093215A (zh) | 基于多普勒信息的雷达对低空慢速小目标的跟踪方法 | |
Jia et al. | A novel approach to target localization through unknown walls for through-the-wall radar imaging | |
CN102288943A (zh) | 基于两视实图像处理的单通道sar-gmti方法 | |
CN109858526A (zh) | 一种目标跟踪中基于传感器的多目标轨迹融合方法 | |
CN113777600A (zh) | 一种多毫米波雷达协同定位跟踪方法 | |
CN109716157A (zh) | 轴偏移估计装置 | |
CN104181524B (zh) | 一种自适应粒子数的多目标粒子滤波检测前跟踪方法 | |
CN105467370A (zh) | 一种复合双基地雷达进动目标isar图像横向定标方法 | |
CN106526549A (zh) | 联合两坐标雷达和三坐标雷达统计量的假目标鉴别方法 | |
CN109143224A (zh) | 一种多目标关联方法和装置 | |
CN105116387A (zh) | 基于位置和多普勒速度信息的pd雷达抗速度拖引方法 | |
CN106468771A (zh) | 一种低可观测高杂波条件下的多目标检测与跟踪方法 | |
CN110673130B (zh) | 一种基于航迹关联的运动目标航迹跟踪方法 | |
CN108645412A (zh) | 一种多传感器自适应航迹起始方法 | |
CN103792515B (zh) | 一种异平台2维雷达与红外传感器量测数据合成方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180309 |
|
RJ01 | Rejection of invention patent application after publication |