CN106249232B - 基于目标运动态势信息数据关联策略的目标跟踪方法 - Google Patents
基于目标运动态势信息数据关联策略的目标跟踪方法 Download PDFInfo
- Publication number
- CN106249232B CN106249232B CN201610721016.0A CN201610721016A CN106249232B CN 106249232 B CN106249232 B CN 106249232B CN 201610721016 A CN201610721016 A CN 201610721016A CN 106249232 B CN106249232 B CN 106249232B
- Authority
- CN
- China
- Prior art keywords
- target
- mark
- tracking
- point mark
- moment
- 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
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
-
- 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
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
本发明公开了基于目标运动态势信息数据关联策略的目标跟踪方法,包括:步骤1,创建待选目标航迹;步骤2,对待选目标航迹进行跟踪处理,并根据目标起始条件建立合理的正式目标航迹;步骤3,利用卡尔曼滤波模型对跟踪目标和正式目标航迹进行滤波估计处理,得到航迹状态估计值。按船舶避碰雷达量测模型计算跟踪目标和正式目标的最小会遇距离点时间TCPA,将TCPA值与目标分离时间门限、目标接近时间门限和目标重叠时间门限进行比较,得到目标间所处的运动状态信息;步骤4:根据目标间所处的运动状态信息决定跟踪目标点迹航迹的关联策略和关联算法,最终得到跟踪目标与点迹的关联关系。
Description
技术领域
本发明属于雷达目标跟踪技术,尤其涉及基于目标运动态势信息数据关联策略的目标跟踪方法。
背景技术
数据关联过程是确定传感器接收到的量测信息和目标源对应关系的过程,贯穿于目标起始、目标维持和目标终结过程。目前,工程上常用的数据关联算法包括两大类,第一类是对最新量测集合进行研究,如Singer和Sea提出了最优邻近算法,其基本原理选择与目标预测位置最近的点迹作为关联点。当目标信噪比较高或者目标稀疏时,该算法跟踪效果较好。而低分辨雷达探测范围广,跟踪目标数量多。加之雷达分辨率低、方位测量精度差,目标极易呈现密集趋势或者目标回波重叠、目标追越、目标交叉现象。当目标处于目标密集区或者目标回波重叠场景下,基于单一位置信息的最优邻近算法目标跟踪效果较差。第二类数据关联算法是对当前时刻以前的所有量测集合进行研究,如Reid提出的多假设算法(MHT)。在目标密集区域,多假设算法利用后续多个扫描周期的量测数据,实现目标的稳定跟踪。Donald B.Reid采用多假设算法模拟仿真了两个目标回波重叠时的跟踪效果(AnAlgorithm for Tracking Multiple Target.IEEE Transctions on Automatic Control,Vol.AC-24,NO.6,December 1979.)。从理论上,多假设算法采用假设航迹分支的后验概率来解决密集区或者杂波区目标数据关联的模糊性。Masamichi Kojima证明多假设算法是一种解决复杂环境下的多目标数据关联准确度最优的算法(A Study of Target TrackingUsing Track-oriented Multiple Hypothesis Tracking.SICE,1998.5:29-31.)。但MHT算法所产生的假设航迹分支数目与虚警率、目标数和所处理的多个扫描周期数呈指数上升关系,因此算法的计算量和存储量都非常大,要实时地应用到对海雷达目标跟踪工程中仍然有一定难度。
发明内容
本发明的目的在于利用目标运动态势信息识别出目标所处环境,根据目标环境信息自适应采用不同的数据关联策略和数据关联算法,提高目标点迹航迹关联正确率,特别是改善多目标会遇、交叉等复杂场景中的跟踪能力。
实现本发明目的的技术解决方案为:基于目标运动态势信息数据关联策略的目标跟踪方法,步骤如下:
包括如下步骤:
步骤1,创建待选目标航迹:以跟踪目标为中心、r为半径的圆区域范围内的搜索点迹信息,搜索半径r与雷达天线周期有关。雷达天线周期越长,搜索半径越大。一般情况下可根据跟踪目标位置设置1~7Km即可。若点迹属于该区域,则利用该点迹创建待选目标航迹;否则不创建待选目标航迹;
步骤2,对待选目标航迹进行跟踪处理:在雷达检测周期内,若待选目标航迹满足目标航迹起始约束条件,则将待选目标航迹转为正式目标航迹,否则删除待选目标航迹;
步骤3,利用卡尔曼滤波模型对跟踪目标和正式目标航迹进行滤波估计处理,得到航迹状态估计值。按船舶避碰雷达量测模型计算跟踪目标和正式目标的最小会遇距离点时间TCPA(time to closest point of approach),将TCPA值与目标分离时间门限、目标接近时间门限和目标重叠时间门限进行比较,得到跟踪目标和正式目标相互之间所处的运动状态信息,如分离、接近或重叠状态;
步骤4:根据跟踪目标和正式目标相互之间所处的运动状态信息决定跟踪目标点迹航迹的关联策略和关联算法,最终得到跟踪目标与点迹的关联关系,从而实现目标跟踪。
本发明步骤1中建立跟踪目标区域范围内的待选目标航迹步骤如下,假设跟踪目标A的距离、方位表示为sA,wA。以(sA,wA)为中心,r为半径确定区域范围。对于落入该区域范围内的点迹建立待选目标航迹。
本发明步骤2包括如下步骤:
步骤2-1,假设在t1时刻创建了待选目标起始点,从t2时刻开始,采用最优邻近算法对待选目标进行点迹航迹关联:
其中,z表示雷达探测值,S-1为待选目标滤波新息协方差,d2(z)为新息加权范数,表示k时刻对k+1时刻目标预测状态的估计;
若d2(z)≤γ成立,则判定待选目标关联到有效点迹,否则,判定待选目标点迹丢失,γ为距离门限值。若待选目标关联到有效点迹,则记录待选目标关联点迹信息,待选目标关联点迹信息包括点迹径向距离g1、点迹方位g2、点迹能量g3、点迹方位扩展g4、点迹距离扩展g5和点迹时间g6,将所述待选目标关联点迹信息组成关联点迹特征向量G,如下式所示:
G=[g1 g2 g3 g4 g5 g6];
步骤2-2,依据卡尔曼滤波公式组,对待选目标的滤波状态和预测状态进行估计:
K(k)=P(k|k-1)H'(k)[H(k)P(k|k-1)H'(k)+Rk]-1
P(k|k-1)=Φ(k|k-1)P(k-1|k-1)Φ'(k|k-1)+Q(k)
P(k|k)=[I-K(k)H(k)]P(k|k-1)
其中,Z(k)表示k时刻传感器探测信息,Φ(k|k-1)表示k-1时刻到k时刻的状态转移矩阵,H(k)表示量测矩阵,U(k)表示输入控制项矩阵,u(k)表示已知输入或者控制信号,R(k)表示零均值、白色高斯量测噪声的协方差,Q(k)为零均值、白色高斯过程噪声的协方差,表示目标k时刻的状态向量,表示目标在k-1时刻对k时刻的预测,K(k)表示k时刻滤波增益,P(k|k)表示k时刻的协方差矩阵,P(k|k-1)表示预测协方差矩阵,I为单位矩阵,H'(k)为量测矩阵的转置矩阵,Φ'(k|k-1)为状态转移矩阵的转置矩阵。
步骤2-3,重复步骤2-1~步骤2-2,直至tn时刻,N为雷达检测周期,假设待选目标在N个雷达检测周期内共有M个关联点迹,关联点迹序列表示为[G1,G2...,GM];
步骤2-4,根据如下公式计算待选目标相邻关联点迹径向距离变化率:
其中,gi1表示第i个雷达扫描周期待选目标关联点迹的第一个特征,即点迹径向距离,gj1表示第j个雷达扫描周期待选目标关联点迹的第一个特征,gi6表示第i个雷达扫描周期待选目标关联点迹的第六个特征,即点迹时间,gj6表示第j个雷达扫描周期待选目标关联点迹的第六个特征,i=j+1,j=1,2,...,M-1。如g11表示第一个扫描周期目标关联点迹的第1个特征,g21表示第二个扫描周期目标关联点迹的第1个特征。tk+1-tk表示vk+1与vk之间的时间差,vk,vk+1表示k,k+1时刻目标径向距离速度,k=1,2,...,n-2,ak表示径向距离变化率;
步骤2-5,根据如下公式计算待选目标相邻关联点迹的方位变化率:
Azik=gi2-gj2,
Azik'=Azik+1-Azik,
其中,gi2和gj2分别表示第i个雷达扫描周期目标关联点迹的方位值(即第二个特征)和第j个雷达扫描周期目标关联点迹的方位值,Azik表示第j个雷达扫描周期目标点迹方位变化值,Azik'表示方位变化率,
步骤2-6,根据步骤2-4和步骤2-5分别计算得到n个目标距离变化率值和方位变化率值,按如下公式计算径向距离变化率的均方根和方位变化率的均方根RMSE:
其中,Xi表示目标距离变化率或者方位变化率,是目标径向距离变化率均值或者方位变化率的均值;
步骤2-7,计算待选目标关联点迹的检测概率p:
步骤2-8,设目标径向距离变化率均方根门限为RdG,目标方位变化率均方根门限为RaG,目标检测概率门限为Pd。若待选目标相邻关联点迹的径向距离变化率的均方根RMSEd、方位变化率的均方根RMSEa和检测概率p均满足以下条件:
RMSEd≤RdG,
RMSEd≤RaG,
p≤Pd,
则将待选目标航迹转为正式目标航迹,并依据卡尔曼滤波模型维持目标跟踪,同时记录目标航迹跟踪过程中的关联点迹特征向量。
本发明步骤3包括如下步骤:
步骤3-1,根据步骤2-1~步骤2-7所述的方法,利用卡尔曼滤波模型对跟踪目标和正式目标航迹进行滤波估计处理,得到航迹状态距离、方位、航向和航速估计值。利用船舶避碰雷达量测模型、跟踪目标和正式目标的距离、方位、航向、航速,计算跟踪目标和正式目标的最小会遇距离点时间TCPA:建立平面坐标系XOY,Y轴为正北方向,设跟踪目标S1位于坐标原点O并以航向d1航速V1匀速运动,跟踪目标周围检测出的正式目标S2以航向d2航速V2匀速运动,T1时刻,正式目标S2相对跟踪目标S1的位置为A点,相对距离为R,相对方位为θ,AA'为两目标相对航向,AA'与Y轴的夹角为φr,相对航速为Vr,采用如下公式计算φr和Vr:
t时刻后,计算正式目标S2相对跟踪目标S1的位移:
由ΔX和ΔY计算出正式目标S2相对跟踪目标S1在t时刻后的距离R(t)和方位θ(t):
其中,Xt0和Yt0分别表示t时刻正式目标S2的起始位置横坐标和t时刻正式目标S2的起始位置纵坐标,ΔX和ΔY分别表示目标S2在t时刻内X方向运动的位置和Y方向运动的位置,在实际工程应用中,通过t时刻雷达对目标S2的探测,可获得位置信息R(t)和θ(t)。通过如下公式计算跟踪目标和正式目标的最小会遇距离点时间TCPA:
TCPA=R(t)|cos(φr-θ(t))|/Vr,
步骤3-2,设GT1为目标回波重叠时间,GT2为目标回波接近时间。若TCPA绝对值小于GT1,则判定跟踪目标与正式目标为会遇状态;若TCPA绝对值大于GT1小于GT2,则判定跟踪目标与正式目标为接近状态;若TCPA绝对值大于GT2,则判定跟踪目标与正式目标为分离状态。
本发明步骤4包括如下步骤:
步骤4-1,若跟踪目标和目标处于分离状态,则采用最优邻近算法对跟踪目标进行数据关联;
步骤4-2,若跟踪目标和正式目标处于接近或者重叠状态,判断跟踪目标和正式目标是否关联到公共点迹。若跟踪目标和正式目标分别关联到不同的点迹,则按正常跟踪处理;若跟踪目标和正式目标关联点迹为公共点迹,则判断下述条件是否同时成立:
其中,和分别表示跟踪目标关联点迹的能量均值和正式目标关联点迹的能量均值,和分别表示跟踪目标关联点迹的回波方位扩展均值和正式目标关联点迹的回波方位扩展均值,Ed为雷达可探测到的公共点迹能量值,ld为雷达可探测到的公共点迹方位扩展值。α,β为比例系数,取值范围分别为0.7≤α≤1和0.7≤β≤1。GE为能量误差门限,Gl为方位扩展误差门限,如果条件同时成立,则点迹是跟踪目标和正式目标回波重叠点迹,点迹不分配给跟踪目标和正式目标,目标分别进行外推处理;如果条件不成立,则点迹来源于跟踪目标或者正式目标,按如下公式计算ΔA和ΔB:
ΔA表示公共点迹的能量与跟踪目标能量均值的相似性,ΔB表示公共点迹的能量与正式目标能量均值的相似性,若ΔA≤ΔB,点迹为跟踪目标关联点迹,若ΔA≥ΔB,点迹为正式目标关联点迹。
本发明步骤4-1中,当跟踪目标和正式目标处于分离状态时,采用最优邻近算法进行数据关联。设跟踪目标关联点迹能量特征集合为{E1,E2,...,En},回波方位扩展特征集合为{l1,l2,...,ln},n为跟踪目标关联点迹的计数值。正式目标关联点迹能量特征集合为{E'1,E2',...,Em'},回波方位扩展特征集合为{l1',l2',...,lm'}。m为正式目标关联点迹的计数值。点迹能量特征和方位扩展特征可采用成熟的点迹提取技术实现,即利用目标回波在方位、距离上的连续性进行点迹凝聚,得到目标方位扩展特征。能量特征通过将参与积累的方位距离对应的能量进行求和获得。
通过如下公式计算跟踪目标关联点迹的能量均值和正式目标关联点迹的能量均值
通过如下公式计算跟踪目标关联点迹的回波方位扩展均值和正式目标关联点迹的回波方位扩展均值
有益效果:本发明与现有技术相比,其显著优点:(1)目标数据关联过程中,不仅考虑到跟踪目标本身的状态估计,还考虑到跟踪目标周围的环境信息,建立了跟踪目标与周围其它目标之间的运动态势关系,实现了多种数据关联算法的自适应切换。(2)采用目标间态势信息判断出目标处于会遇、交叉等复杂跟踪场景时,结合点迹多元特征信息辅助数据关联算法,降低了基于单一位置信息关联导致的关联错误率。(3)与多假设算法相比,基于目标运动态势信息的数据关联策略方法的计算量、存储量较小,适合工程应用。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
图1为创建待选目标示意图。
图2为待选目标转正式目标处理流程图。
图3为船舶避碰雷达量测模型。
图4为目标运动态势判断处理流程图。
图5为TCPA运动态势判断说明图。
图6为多目标追越跟踪示意图。
图7为多目标追越跟踪过程中跟踪目标态势变化曲线图。
图8为多目标交叉跟踪示意图。
图9为多目标交叉跟踪过程中跟踪目标态势变化曲线图。
具体实施方式
第一步,结合图1说明待选目标航迹创建过程。以跟踪目标A为中心,r为半径的圆区域范围内的搜索点迹信息。其中,搜索半径r与雷达天线周期有关。雷达天线周期越长,搜索半径越大。一般情况下可根据跟踪目标位置设置1~7Km即可。若点迹落入该区域,则创建待选目标航迹;否则,不创建待选目标航迹。(雷达探测目标过程中,易受到杂波、目标等因素影响,同一时刻会在跟踪目标波门内出现多个点迹。图1中两个t1时刻的意义是在相同时刻内,目标波门内同时出现了两个点迹,两个点迹时间一致。同理,t3时刻也出现两个点迹。)
第二步,对跟踪目标A区域范围内的待选目标航迹进行跟踪处理。在若干检测周期内,若待选目标满足目标航迹起始条件,则将待选目标转为正式目标航迹B,否则删除待选目标航迹,具体如下所述。
(2-1):t1时刻已创建了待选目标起始点。从t2时刻开始,利用最优邻近算法公式对待选目标进行点迹航迹关联。其中,z表示雷达探测值,S-1为待选目标滤波新息协方差,d2(z)为新息加权范数。
若d2(z)≤γ成立,则待选目标关联到有效点迹,否则,待选目标点迹丢失。γ为距离门限值。若待选目标关联到点迹,则记录待选目标关联点迹信息。待选目标关联点迹信息主要由六部分组成:点迹径向距离g1、点迹方位g2、点迹能量g3、点迹方位扩展g4、点迹距离扩展g5,点迹时间g6,组成关联点迹特征向量G,如式所示。
G=[g1 g2 g3 g4 g5 g6]
(2-2):待选目标数据关联处理完毕之后,依据卡尔曼滤波公式组,对待选目标的滤波状态、预测状态进行估计。
K(k)=P(k|k-1)H'(k)[H(k)P(k|k-1)H'(k)+Rk]-1
P(k|k-1)=Φ(k|k-1)P(k-1|k-1)Φ'(k|k-1)+Q(k)
P(k|k)=[I-K(k)H(k)]P(k|k-1)
(2-3):重复(2-1)~(2-2)步骤直至tN时刻,N为检测周期。考虑到雷达虚警概率和保证待选目标有效关联点迹样本数,N一般取5~7。
(2-4):结合图2说明待选目标判断为正式目标的过程。假设待选目标在N个雷达探测周期内关联点迹序列可表示为[G1,G2...,GM],M为待选目标关联点迹个数,Gi为点迹特征向量。
(2-4-1):按公式计算待选目标相邻点迹径向距离变化率。
式中gij表示第i个雷达扫描周期目标关联点迹的第j个特征,i=j+1,j=1,2,...,m-1。vk,vk+1表示k,k+1时刻目标径向距离速度,k=1,2,...,n-2,tk+1-tk为vk+1与vk之间的时间差,ak为目标径向距离变化率值。
(2-4-2):按公式计算待选目标相邻点迹的方位变化率。
Azik=gi2-gj2,i=j+1
Azik'=Azik+1-Azik
gij表示第i个雷达扫描周期目标关联点迹的第j个特征,Azik'为方位变化率。
(2-4-3):计算待选目标径向距离变化率和方位变化率的均方根值。
式中Xi表示目标径向距离变化率值或者方位变化率值,是目标径向距离变化率均值或者方位变化率的均值,
(2-4-4):计算待选目标点迹的检测概率。
若待选目标径向距离变化率均方根、方位变化率均方根和点迹检测概率均满足以下目标航迹起始约束条件:
RMSEd≤RdG
RMSEd≤RaG
p≤Pd
则待选目标转为正式目标航迹。并依据卡尔曼滤波模型维持目标跟踪,同时记录目标跟踪过程中的关联点迹特征向量。
第三步:结合图3说明,利用船舶避碰雷达量测模型计算跟踪目标A和目标B的最小会遇距离点时间(TCPA)方法。平面坐标系XOY中,Y轴为正北方向。设跟踪目标S1位于坐标原点,以航向d1航速V1匀速运动。跟踪目标周围检测出的目标S2以航向d2航速V2匀速运动。T1时刻,S2相对S1的位置为A点,相对距离为R,相对方位θ。AA'为两目标相对航向,与Y轴的夹角为φr,相对航速为Vr。由几何关系可计算出φr和Vr,
t时刻后,计算目标S2相对目标S1的位移,
由ΔX,ΔY推导出S2相对S1距离和方位,
式中,Xt0和Yt0为t时刻正式目标S2的起始位置,ΔX和ΔY分别表示目标S2在t时刻内X方向和Y方向运动的位置。在实际工程应用中,可通过t时刻雷达对目标S2的探测获得位置信息R(t)和θ(t)。由相对距离和方位,按公式计算TCPA,
TCPA=R(t)|cos(φr-θ(t))|/Vr
结合图4说明依据TCPA值判断跟踪目标A与其周围目标B之间的运动态势流程图。GT1为目标回波重叠时间,GT2为目标回波接近时间。结合图5具体说明目标态势判断过程。设GT1=20,GT2=40。若TCPA绝对值小于GT1,则跟踪目标A与目标B判断为会遇状态;若TCPA绝对值大于GT1小于GT2,则跟踪目标A与目标B判断为接近状态;若TCPA绝对值大于GT2,则跟踪目标A与目标B判断为分离状态。一般情况,依据雷达周期设置GT1和GT2。本发明在工程应用中GT1设置为1分钟,GT2设置为5分钟。结合图6说明多目标追越跟踪过程中跟踪目标态势曲线图的变化,图7为追越目标0009批的态势曲线变化图。当跟踪目标0009批的TCPA大于GT1时,跟踪目标0009态势曲线为分离状态。在雷达第10至35扫描周期内,跟踪目标0009批与正式目标回波未分离,0009批跟踪目标TCPA值小于GT1时,态势曲线保持会遇状态。当雷达在第36个扫描周期时,跟踪目标完成追越过程。0009批跟踪目标和正式目标态势曲线重新更新为分离状态。结合图8说明多目标交叉跟踪过程中跟踪目标态势曲线图的变化,图9为交叉目标0015批的态势曲线变化图。当跟踪目标0015批的TCPA大于GT1时,跟踪目标0015批态势曲线为分离状态。当跟踪目标0015批和正式目标0018批回波接近时,TCPA值计算结果介入GT1和GT2之间,跟踪目标0015批态势曲线显示为接近状态。雷达第57至77扫描周期内,跟踪目标0015和正式目标0018批回波重叠,TCPA值计算结果小于GT1,跟踪目标0015批态势曲线为会遇状态。随着雷达扫描和目标运动,跟踪目标和正式目标回波逐渐分离,跟踪目标0015的态势曲线显示由会遇,变为接近,最终变为分离状态。
第四步,根据目标运动态势信息决定跟踪目标A点迹航迹的关联策略和关联算法。若跟踪目标A与其周围目标B处于分离状态,则说明跟踪目标A和目标B不会出现会遇或者交叉场景。跟踪目标A按最优邻近算法进行数据关联。若跟踪目标A与其周围目标B处于接近或者重叠状态时,可能出现两种点迹关联情况。第一种,目标各自关联不同的点迹;第二种,目标关联相同的点迹。对于第一种情况,无需考虑点迹特征。第二种情况需考虑点迹特征。假设跟踪目标A点迹能量特征集合为{E1,E2,...,En},方位扩展特征集合为{l1,l2,...,ln},n为跟踪目标关联点迹的计数值。目标B点迹能量特征集合为{E'1,E2',...,Em'},方位扩展特征集合为{l1',l2',...,lm'},m为跟踪目标关联点迹的计数值。
(4-1):按公式分别计算目标A、目标B关联点迹的能量均值,
(4-2):按公式分别计算目标A、目标B关联点迹的回波方位扩展均值。
(4-3):将公共点迹能量Ed、方位扩展ld分别与跟踪目标A、目标B能量均值 方位扩展均值行比较。若满足下列条件,说明公共点迹是两目标回波重叠点迹,则公共点迹不分配给目标A、目标B,两目标外推处理。否则,按第四步处理。式中GE为能量门限,Gl为能量门限。一般情况,根据雷达探测的目标能量和方位扩展进行设置GE和Gl。
(4-4):计算下列等式。
若ΔA≤ΔB,点迹为跟踪目标A关联点迹。若ΔA≥ΔB,点迹为目标B关联点迹。多元特征数据关联方法处理完毕。
尽管本发明就优选实施方式进行了示意和描述,但本领域的技术人员应当理解,只要不超出本发明的权利要求所限定的范围,可以对本发明进行各种变化和修改。
Claims (4)
1.基于目标运动态势信息数据关联策略的目标跟踪方法,其特征在于,包括如下步骤:
步骤1,创建待选目标航迹:在以跟踪目标为中心、r为半径的圆区域范围内搜索点迹信息,若点迹属于该区域,则利用该点迹创建待选目标航迹;否则不创建待选目标航迹;
步骤2,对待选目标航迹进行跟踪处理:在雷达检测周期内,若待选目标航迹满足目标航迹起始约束条件,则将待选目标航迹转为正式目标航迹,否则删除待选目标航迹;
步骤3,利用卡尔曼滤波模型对跟踪目标和正式目标航迹进行滤波估计处理,按船舶避碰雷达量测模型计算跟踪目标和正式目标的最小会遇距离点时间TCPA,将TCPA值与目标分离时间门限、目标接近时间门限和目标重叠时间门限进行比较,得到跟踪目标和正式目标相互之间所处的运动状态信息;
步骤4:根据跟踪目标和正式目标相互之间所处的运动状态信息决定跟踪目标点迹航迹的关联策略和关联算法,最终得到跟踪目标与点迹的关联关系,从而实现目标跟踪;
步骤2包括如下步骤:
步骤2-1,假设在t1时刻创建了待选目标起始点,从t2时刻开始,采用最优邻近算法对待选目标进行点迹航迹关联:
其中,z表示雷达探测值,S-1为待选目标滤波新息协方差,d2(z)为新息加权范数,表示k时刻对k+1时刻目标预测状态的估计;
若d2(z)≤γ成立,则判定待选目标关联到有效点迹,否则,判定待选目标点迹丢失,γ为距离门限值,若待选目标关联到有效点迹,则记录待选目标关联点迹信息,待选目标关联点迹信息包括点迹径向距离g1、点迹方位g2、点迹能量g3、点迹方位扩展g4、点迹距离扩展g5和点迹时间g6,将所述待选目标关联点迹信息组成关联点迹特征向量G,如下式所示:
G=[g1 g2 g3 g4 g5 g6];
步骤2-2,依据卡尔曼滤波公式,对待选目标的滤波状态和预测状态进行估计:
K(k)=P(k|k-1)H'(k)[H(k)P(k|k-1)H'(k)+R(k)]-1
P(k|k-1)=Φ(k|k-1)P(k-1|k-1)Φ'(k|k-1)+Q(k)
P(k|k)=[I-K(k)H(k)]P(k|k-1)
其中,Z(k)表示k时刻传感器探测目标的坐标位置信息,Φ(k|k-1)表示k-1时刻到k时刻目标运动状态转移矩阵,H(k)表示目标量测矩阵,U(k)表示输入控制项矩阵,u(k)表示已知输入或者控制信号,R(k)表示零均值、白色高斯量测噪声的协方差,Q(k)为零均值、白色高斯过程噪声的协方差,表示目标k时刻的状态估计向量,表示k-1时刻对k时刻的目标状态预测,K(k)表示k时刻滤波增益,P(k|k)表示k时刻的协方差矩阵,P(k|k-1)表示预测协方差矩阵,I为单位矩阵,H'(k)为目标量测矩阵的转置矩阵,Φ'(k|k-1)为状态转移矩阵的转置矩阵;
步骤2-3,重复步骤2-1~步骤2-2,直至tN时刻,N为雷达检测周期,假设待选目标在N个雷达检测周期内共有M个关联点迹,关联点迹序列表示为[G1,G2...GM];
步骤2-4,根据如下公式计算待选目标相邻关联点迹径向距离变化率:
其中,gi1表示第i个雷达扫描周期待选目标关联点迹的第一个特征,即点迹径向距离,gj1表示第j个雷达扫描周期待选目标关联点迹的第一个特征,gi6表示第i个雷达扫描周期待选目标关联点迹的第六个特征,即点迹时间,gj6表示第j个雷达扫描周期待选目标关联点迹的第六个特征,i=j+1,j=1,2,...,M-1,tk+1-tk表示vk+1与vk之间的时间差,vk表示k时刻目标径向距离速度,k=1,2,...,n-2,ak表示目标径向距离变化率;
步骤2-5,根据如下公式计算待选目标相邻关联点迹的方位变化率:
Azik=gi2-gj2,
Azik'=Azik+1-Azik,
其中,gi2和gj2分别表示第i个雷达扫描周期目标关联点迹的方位值和第j个雷达扫描周期目标关联点迹的方位值,Azik表示第j个雷达扫描周期目标点迹方位变化值,Azik'表示方位变化率;
步骤2-6,根据步骤2-4和步骤2-5分别计算得到n个目标距离变化率值和方位变化率值,按如下公式计算径向距离变化率的均方根和方位变化率的均方根RMSE:
其中,Xi表示目标距离变化率或者方位变化率,是目标径向距离变化率均值或者方位变化率的均值;
步骤2-7,计算待选目标关联点迹的检测概率p:
步骤2-8,设目标径向距离变化率均方根门限为RdG,目标方位变化率均方根门限为RaG,目标检测概率门限为Pd,若待选目标相邻关联点迹的径向距离变化率的均方根RMSEd、方位变化率的均方根RMSEa和检测概率p均满足以下条件:
RMSEd≤RdG,
RMSEd≤RaG,
p≤Pd,
则将待选目标航迹转为正式目标航迹,并依据卡尔曼滤波模型维持目标跟踪,同时记录目标航迹跟踪过程中的关联点迹特征向量。
2.根据权利要求1所述的方法,其特征在于,步骤3包括如下步骤:
步骤3-1,利用船舶避碰雷达量测模型计算跟踪目标和正式目标的最小会遇距离点时间TCPA:建立平面坐标系XOY,Y轴为正北方向,设跟踪目标S1位于坐标原点O并以航向d1航速V1匀速运动,跟踪目标周围检测出的正式目标S2以航向d2航速V2匀速运动,T1时刻,正式目标S2相对跟踪目标S1的位置为A点,相对距离为R,相对方位为θ,AA'为两目标相对航向,AA'与Y轴的夹角为φr,相对航速为Vr,采用如下公式计算φr和Vr:
t时刻后,计算正式目标S2相对跟踪目标S1的位移:
由ΔX和ΔY计算出正式目标S2相对跟踪目标S1在t时刻后的距离R(t)和方位θ(t):
其中,Xt0和Yt0分别表示t时刻正式目标S2的起始位置横坐标和t时刻正式目标S2的起始位置纵坐标,ΔX和ΔY分别表示目标S2在t时刻内X方向运动的位置和Y方向运动的位置,通过如下公式计算跟踪目标和正式目标的最小会遇距离点时间TCPA:
TCPA=R(t)|cos(φr-θ(t))|/Vr;
步骤3-2,设GT1为目标回波重叠时间,GT2为目标回波接近时间,若TCPA绝对值小于GT1,则判定跟踪目标与正式目标为会遇状态;若TCPA绝对值大于GT1小于GT2,则判定跟踪目标与正式目标为接近状态;若TCPA绝对值大于GT2,则判定跟踪目标与正式目标为分离状态。
3.根据权利要求2所述的方法,其特征在于,步骤4包括如下步骤:
步骤4-1,若跟踪目标和目标处于分离状态,则采用最优邻近算法对跟踪目标进行数据关联;
步骤4-2,若跟踪目标和正式目标处于接近或者重叠状态,则需判断跟踪目标和正式目标是否关联到公共点迹,若跟踪目标和正式目标分别关联到不同的点迹,则按正常跟踪处理;若跟踪目标和正式目标关联点迹为公共点迹,则判断下述条件是否同时成立:
其中,和分别表示跟踪目标关联点迹的能量均值和正式目标关联点迹的能量均值,和分别表示跟踪目标关联点迹的回波方位扩展均值和正式目标关联点迹的回波方位扩展均值,Ed为雷达可探测到的公共点迹能量值,ld为雷达可探测到的公共点迹方位扩展值,α和β为比例系数,GE为能量误差门限,Gl为方位扩展误差门限,如果条件同时成立,则点迹是跟踪目标和正式目标回波重叠点迹,点迹不分配给跟踪目标和正式目标,目标分别进行外推处理;如果条件不成立,则点迹来源于跟踪目标或者正式目标,按如下公式计算ΔA和ΔB:
ΔA表示公共点迹的能量与跟踪目标能量均值的相似性,ΔB表示公共点迹的能量与正式目标能量均值的相似性,若ΔA≤ΔB,点迹为跟踪目标关联点迹,若ΔA>ΔB,点迹为正式目标关联点迹。
4.根据权利要求3所述的方法,其特征在于,步骤4-1中,当跟踪目标和正式目标处于分离状态时,采用最优邻近算法进行数据关联,设跟踪目标关联点迹能量特征集合为{E1,E2,...,En},回波方位扩展特征集合为{l1,l2,...,ln},n为跟踪目标关联点迹的计数值,正式目标关联点迹能量特征集合为{E'1,E2',...,Em'},回波方位扩展特征集合为{l1',l2',...,lm'},m为正式目标关联点迹的计数值,利用目标回波在方位、距离上的连续性进行点迹凝聚,得到目标方位扩展特征,能量特征通过将参与积累的方位距离对应的能量进行求和获得,
通过如下公式计算跟踪目标关联点迹的能量均值和正式目标关联点迹的能量均值
通过如下公式计算跟踪目标关联点迹的回波方位扩展均值和正式目标关联点迹的回波方位扩展均值
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610721016.0A CN106249232B (zh) | 2016-08-24 | 2016-08-24 | 基于目标运动态势信息数据关联策略的目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610721016.0A CN106249232B (zh) | 2016-08-24 | 2016-08-24 | 基于目标运动态势信息数据关联策略的目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106249232A CN106249232A (zh) | 2016-12-21 |
CN106249232B true CN106249232B (zh) | 2018-08-21 |
Family
ID=57595546
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610721016.0A Active CN106249232B (zh) | 2016-08-24 | 2016-08-24 | 基于目标运动态势信息数据关联策略的目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106249232B (zh) |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180275265A1 (en) * | 2017-03-23 | 2018-09-27 | GM Global Technology Operations LLC | Target tracking using region covariance |
CN106932771A (zh) * | 2017-03-30 | 2017-07-07 | 成都紫瑞青云航空宇航技术有限公司 | 一种雷达仿真目标航迹跟踪方法及*** |
CN107192998B (zh) * | 2017-04-06 | 2019-05-31 | 中国电子科技集团公司第二十八研究所 | 一种基于协方差指标函数的自适应分布式航迹数据融合方法 |
CN107911196B (zh) * | 2017-10-27 | 2020-07-14 | 南京莱斯电子设备有限公司 | 一种雷达航迹报文传输方法 |
CN107896136B (zh) * | 2017-10-27 | 2020-09-08 | 南京莱斯电子设备有限公司 | 一种雷达航迹报文编码方法 |
CN108344992B (zh) * | 2017-12-20 | 2020-03-27 | 北京华航无线电测量研究所 | 一种用于车载毫米波雷达的多目标跟踪方法 |
CN108226863B (zh) * | 2017-12-27 | 2021-12-03 | 武汉理工大学 | 一种单目卫星跟踪定位方法 |
CN110031834A (zh) * | 2018-01-12 | 2019-07-19 | 西安艾索信息技术有限公司 | 一种改进的多目标雷达航迹处理方法 |
CN108344997B (zh) * | 2018-01-16 | 2021-10-29 | 华域汽车***股份有限公司 | 一种基于点迹特征的道路护栏快速检测方法 |
CN108303684B (zh) * | 2018-01-31 | 2021-09-03 | 长沙深之瞳信息科技有限公司 | 基于径向速度信息的地基监视雷达多目标跟踪方法 |
CN108680918B (zh) * | 2018-05-18 | 2020-08-28 | 森思泰克河北科技有限公司 | 应用于雷达的测速方法、测速装置及电子设备 |
CN109001725B (zh) * | 2018-06-07 | 2020-11-10 | 中国人民解放军海军工程大学 | 一种海上无人艇海上多目标跟踪方法 |
CN109031269B (zh) * | 2018-06-08 | 2020-07-07 | 上海西井信息科技有限公司 | 基于毫米波雷达的定位方法、***、设备及存储介质 |
CN109509214B (zh) * | 2018-10-15 | 2021-07-16 | 杭州电子科技大学 | 一种基于深度学习的船舶目标跟踪方法 |
CN110400491A (zh) * | 2019-06-10 | 2019-11-01 | 北京海兰信数据科技股份有限公司 | 一种宽水域多目标辅助避碰决策方法及决策*** |
CN110646798B (zh) * | 2019-09-23 | 2021-09-28 | 森思泰克河北科技有限公司 | 目标航迹关联方法、雷达和终端设备 |
CN110907903B (zh) * | 2019-10-30 | 2023-08-01 | 中电科技(合肥)博微信息发展有限责任公司 | 一种基于多传感器的自适应跟踪处理方法 |
CN110940971B (zh) * | 2019-11-06 | 2022-05-31 | 四川川大智胜软件股份有限公司 | 一种雷达目标点迹录取方法、装置及存储介质 |
CN111323757B (zh) * | 2019-12-30 | 2022-04-05 | 北京海兰信数据科技股份有限公司 | 一种航海雷达目标检测方法及装置 |
CN111289954B (zh) * | 2020-03-31 | 2022-03-15 | 四川长虹电器股份有限公司 | 用于毫米波雷达目标跟踪的点云划分与航迹匹配方法 |
CN111707997B (zh) * | 2020-06-03 | 2022-05-27 | 南京慧尔视智能科技有限公司 | 雷达目标跟踪方法、装置、电子设备及存储介质 |
CN111751816B (zh) * | 2020-07-07 | 2023-03-21 | 北京无线电测量研究所 | 一种基于点迹特征信息的雷达航迹起始方法 |
CN112363167A (zh) * | 2020-11-02 | 2021-02-12 | 重庆邮电大学 | 一种基于毫米波雷达与单目相机融合的扩展目标跟踪方法 |
CN112505681B (zh) * | 2020-11-02 | 2022-11-25 | 南京理工大学 | 一种基于tas的四面二维相控阵雷达多目标跟踪处理方法 |
CN112800349B (zh) * | 2021-02-02 | 2022-02-18 | 中华人民共和国广东海事局 | 水上移动目标的运动状态获取方法、装置、设备和介质 |
CN113744309B (zh) * | 2021-08-11 | 2024-05-24 | 中国科学院空天信息创新研究院 | 一种基于运动状态变化感知的机动目标跟踪方法 |
CN113702940B (zh) * | 2021-09-18 | 2022-06-21 | 中国人民解放军63921部队 | 基于多元特征信息分层融合空间集群目标分辨方法及应用 |
CN113721237B (zh) * | 2021-11-02 | 2022-02-08 | 南京雷电信息技术有限公司 | 一种多隶属度的目标智能匹配算法 |
CN117554922B (zh) * | 2024-01-12 | 2024-03-26 | 航天宏图信息技术股份有限公司 | 目标航迹的关联方法、装置、电子设备和计算机存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7142981B2 (en) * | 2003-08-05 | 2006-11-28 | The Boeing Company | Laser range finder closed-loop pointing technology of relative navigation, attitude determination, pointing and tracking for spacecraft rendezvous |
US7417583B2 (en) * | 2006-10-30 | 2008-08-26 | Raytheon Company | Methods and apparatus for providing target altitude estimation in a two dimensional radar system |
CN103727931A (zh) * | 2013-12-31 | 2014-04-16 | 南京理工大学 | 一种基于改进逻辑的航迹起始方法 |
CN105388472B (zh) * | 2015-11-20 | 2018-02-13 | 中国电子科技集团公司第二十八研究所 | 一种对海雷达***中交错重叠目标跟踪方法 |
-
2016
- 2016-08-24 CN CN201610721016.0A patent/CN106249232B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106249232A (zh) | 2016-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106249232B (zh) | 基于目标运动态势信息数据关联策略的目标跟踪方法 | |
CN109100714B (zh) | 一种基于极坐标系的低慢小目标跟踪方法 | |
CN103729859B (zh) | 一种基于模糊聚类的概率最近邻域多目标跟踪方法 | |
Newman et al. | Towards constant-time SLAM on an autonomous underwater vehicle using synthetic aperture sonar | |
CN106407677B (zh) | 一种测量数据丢失情况下的多目标跟踪方法 | |
CN109633590B (zh) | 基于gp-vsmm-jpda的扩展目标跟踪方法 | |
CN106980114A (zh) | 无源雷达目标跟踪方法 | |
CN107192998A (zh) | 一种基于协方差指标函数的自适应分布式航迹数据融合方法 | |
KR101628154B1 (ko) | 수신 신호 세기를 이용한 다중 표적 추적 방법 | |
CN110132284B (zh) | 一种基于深度信息的全局定位方法 | |
CN108957437B (zh) | 一种基于环境态势分析的雷达目标跟踪方法 | |
CN107861123A (zh) | 一种穿墙雷达在复杂环境下对多运动目标实时跟踪的方法 | |
CN105301584B (zh) | 同时解距离模糊的ipphdf机动多目标跟踪方法 | |
CN108734725A (zh) | 基于高斯过程的概率数据关联滤波扩展目标跟踪方法 | |
CN107436434B (zh) | 基于双向多普勒估计的航迹起始方法 | |
CN108614580A (zh) | 一种无人机目标跟踪中的分层避障控制方法 | |
CN109858526A (zh) | 一种目标跟踪中基于传感器的多目标轨迹融合方法 | |
CN108645412A (zh) | 一种多传感器自适应航迹起始方法 | |
CN109521420A (zh) | 基于多特征匹配的多目标跟踪方法 | |
CN114114242A (zh) | 一种基于点选航迹及航迹外推的雷达目标追赶跟踪方法 | |
CN105353352B (zh) | 改进搜索策略的mm‑pphdf机动多目标跟踪方法 | |
CN109509207A (zh) | 一种对点目标和扩展目标进行无缝跟踪的方法 | |
Qin et al. | Research on information fusion structure of radar and AIS | |
CN106772357A (zh) | 信噪比未知条件下的ai‑phd滤波器 | |
Chen et al. | Research on AIS and radar information fusion method based on distributed Kalman |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |