一种基于汽车防撞雷达的道路多目标跟踪方法 |
|||||||
申请号 | CN201610113770.6 | 申请日 | 2016-02-29 | 公开(公告)号 | CN105699964A | 公开(公告)日 | 2016-06-22 |
申请人 | 无锡南理工科技发展有限公司; | 发明人 | 肖泽龙; 龚迅; 许建中; | ||||
摘要 | 本 发明 公开一种基于 汽车 防撞雷达的道路多目标 跟踪 方法,通过 有向图 切换的无迹卡尔曼概率 密度 假设(DS-UKGMPHD)方法实现。主要包括获取先验信息、建立观测与运动模型、建立有向图Φ、建立自适应有向图切换准则、初始化及计算混合参数、各目标模型的UKGMPHD滤波、估计融合以及输出估计轨迹等步骤。本发明不仅适用于非线性、复杂运动状态的多个目标,而且兼顾机动目标跟踪实时性和可实现性,进一步提升汽车防撞雷达系统的实时性和可靠性,保障道路安全,在汽车防撞雷达 信号 处理领域具有广阔的应用前景。 | ||||||
权利要求 | 1.一种基于汽车防撞雷达的道路多目标跟踪方法,具体步骤如下: |
||||||
说明书全文 | 一种基于汽车防撞雷达的道路多目标跟踪方法技术领域背景技术[0002] 随着我国汽车的保有量的增加、汽车的行驶速度的提高,汽车的碰撞事故越来越多。近年来,汽车预警技术成为研究开发的热点,汽车防撞雷达作为一项主动防撞措施,能够追踪汽车前方一定距离内的单个或多个目标,获取其与本车的距离、速度和方位角等信息,并可根据预警系统在危险情况下发出报警信号,以提示驾驶员采取措施,避免由于外后视镜视野不足或驾驶员注意力不集中发生的撞车事故。 [0003] 针对汽车防撞雷达道路多目标跟踪问题,现有技术的解决方法可归纳为两类。第一类是经典方法,主要基于滤波、数据关联及航迹管理;中国专利CN105000019A公开了一种用于检测、追踪和估计静止的路边物体的方法和系统,将来自几个前视雷达传感器的数据点用来确定是否存在静止的路边物体之前合并、归类和预分类数据点来减少数据量,通过数据拟合算法以估计物体的参数,但该专利仅针对路边静止的多目标进行跟踪与估计,未涉及机动目标的追踪问题;此外,随着道路车辆数目的增加,数据关联算法计算量将呈“组合爆炸”式增长,导致防撞系统实时性、准确度欠佳。第二类方法是基于随机有限集(RFS)理论,包括概率密度假设(PHD)方法及势概率假设密度(CPHD)方 法;上述理论均可基于高斯假设(GM)和基于非高斯假设的贯序蒙特卡罗(SMC)方法实现;将汽车防撞雷达道路目标跟踪置于目标检测步骤,如恒虚警检测、快速傅里叶变换等之后,虚警概率已得到一定程度的降低,故从实时性和实用性角度出发,可使用高斯混合概率假设密度(GMPHD)方法作为一种基于高斯假设的跟踪滤波方法。张昱、宋骊平、虎小龙在“基于概率假设密度的汽车防撞雷达多目标跟踪”(现代雷达,2014,36(6):82-87.)一文中提出使用多模型GMPHD方法进行道路目标跟踪,然而该法存在以下缺陷:首先,GMPHD法仅适用于线性条件,而雷达极坐标形式的量测方程导致直角坐标下的跟踪问题为非线性;其次,基于少量模型的跟踪方法,难以适应实际情况中复杂的、运动的、多个道路目标。 [0004] 针对上述问题,非汽车防撞雷达的多目标跟踪文献给出了解决方案,但并未全面解决这些问题。郝燕玲,孟凡彬,孙枫,沈锋等在“基于IMM的UK-GMPHDF算法在多机动目标跟踪中的应用”(系统工程理论与实践,2011,11:2225-2233)一文中,针对三维空中目标,基于交互多模型(IMM),在高斯混合概率假设密度滤波器(GMPHDF)中加入无迹卡尔曼滤波(UKF)解决非线性多目标跟踪问题,但仿真使用模型仅为匀速模型、“当前”模型,且在三维空间实施;中国专利CN104020466A公开了一种基于变结构多模型的机动目标跟踪方法,利用现有的模型扩展方法BMA的候选模型集与现有的期望模式扩展EMA中求期望模型的方法构成一种新的模型集合自适应策略,从而利用通用的变结构交互多模型算法实现单目标的多运动状态跟踪,但无多目标估计能力。中国专利105261036A公开了一种基于匹配的目标跟踪方法,利用FPGA电路进行去均值归一化积运算,从而实现对目标的跟踪,但该专利涉及到跟踪图像,且主要依赖硬件实施,无法解决实际情况中复杂的多个运动目标跟踪问题。 发明内容[0005] 本发明的目的在于克服现有技术的不足,提供一种基于汽车防撞雷达的道路多目标跟踪方法,实现防撞雷达对道路多目标的数目、位置的实时可靠估计与准确跟踪。 [0006] 本发明采用的技术方案是:以汽车防撞雷达对多目标探测后得到的数据为基础,基于变结构交互式多模型(IMM)方法,采用有向图切换(DS)法,在多个运动模型间实时选择切换合适的模型子集,并进行基于无迹变换的卡尔曼概率密度假设(UKGMPHD)滤波。 [0007] 本发明具体技术方案为: [0008] 一种基于汽车防撞雷达的道路多目标跟踪方法,其步骤如下: [0009] 步骤1、获取先验信息,建立第一观测与运动模型。所述先验信息为通过汽车防撞雷达获取的目标状态信息。 [0010] 确定状态方程 量测方程 及其道路目标的多种模型。依据实际情况选择不同模型,状态方程建立为: [0011] 其中, 是目标模型的状态转移矩阵,下标m代表不同目标模型,下标k表示时刻k,下标i表示不同目标。假设噪声服从高斯分布,有 [0012] 汽车防撞雷达观测向量 θk为观测方位角,rk为目标与观测点间距离。 [0013] 步骤2、建立第二观测与运动模型,初始化有向图模型集,建立有向图Φ。 [0014] 2(a)汽车防撞雷达车辆目标限于二维平面运动,目标i在k时刻的非线性量测模型如下: [0015] [0016] 其中, 为高斯分布的量测噪声,xs,k,ys,k为自车观测位置坐标。 [0017] 2(b)建立有向图Φ,及其有向子图Φ(l){βl-1,βl,βl+1}。 [0018] 步骤3、建立自适应有向图切换准则。设k-1时刻子图为Φ(k-1),切换规则如下: [0019] [0020] 式中,τDS是模型转换的阈值概率, 表示k-1时刻概率β(i±1); [0021] 步骤4、初始化及计算混合参数。 [0022] 4(a)时刻k=0时,每一向量子图Dn中含M个模型。对于当前模型q(q=1,...,M),初始化目标状态 目标协方差矩阵 模型概率 [0023] 4(b)设当前时刻k-1,下一时刻为k。则对于k>0时刻的第i个高斯分量,有[0024] [0025] 其中, 为k-1时刻的高斯分量i模型p的出现概率, [0026] 进一步地,当前有向量子图内模型数M,交互初始PHD函数为: [0027] [0028] 混合权值、状态、协方差分别为: [0029] [0030] [0031] [0032] 上式中, [0033] 步骤5、各目标模型的UKGMPHD滤波。 [0034] 以各高斯分量i在k-1时刻所得混合参数作为输入进行UKGMPHD滤波,得k时刻跟踪权值,状态及协方差矩阵 [0035] 进一步地对步骤5描述,包括以下过程: [0036] 5(a)对k=0时刻,初始化目标概率密度,有: [0037] [0038] 5(b)预测。对k>0时刻,k-1时刻,高斯分量i总数Jk-1,概率密度假设函数Dk-1(x)为: [0039] 进一步地,所述预测概率密度假设函数Dk|k-1(x)依据对象的不同分为存活目标、新生目标及衍生目标概率密度假设函数。本发明道路多目标跟踪通常不涉及衍生目标,故概率密度假设函数为:Dk|k-1(x)=DS,k|k-1(x)+γk(x) [0040] 存活目标、新生目标概率密度假设函数分别如下: [0041] [0042] [0043] 使用无迹卡尔曼滤波的思想滤波。依据状态方程,一步预测均值 方差 如下: [0044] [0045] [0046] 式中,l=0,……,2n为无迹滤波的Sigma点,L=2n+1为Sigma点总数。则k时刻预测概率密度假设函数为: [0047] 5(c)更新。 [0048] 对于雷达获取的非线性量测方程,一步预测均值、方差分别为: [0049] [0050] [0051] 结合步骤5(b)结果Dk|k-1(x)及检测概率PD,k,k时刻的更新概率密度假设函数Dk(x)为: [0052] [0053] DD,k(x;z)、估计均值 方差 可表示为: [0054] [0055] 其中, [0056] 高斯分量i的估计均值和方差如下: [0057] [0058] [0059] 上式中,状态增益矩阵、协方差分别为: [0060] [0061] [0062] 5(d)剪枝。 [0063] 对于更新得到的高斯分布进行剪枝,若高斯分布小于修剪阈值τ,直接予以删除。公式为: [0064] 5(e)合并。 [0065] 对于分布接近的若干高斯分布,合并为一个新的分布,即集合: [0066] [0067] 合并后,新的高斯分布为: [0068] [0069] [0070] 5(f)目标数目与状态估计。 [0071] 目标数目估计为融合估计的权值,即: [0072] 取整可得目标数目。按照上述步骤输出合并后状态记作 [0073] 进一步地,在步骤5完成后,存在模型概率更新步骤;完成模型概率更新后,返回步骤4。所述模型概率更新步骤如下: [0074] k时刻模型(q=1,...,M)下高斯分量i的似然函数为 模型q更新概率为: [0075] [0076] [0077] 步骤6、估计融合。估计融合PHD函数: [0078] [0079] 其中,Jk为高斯分量个数,高斯分量i的权值,均值,协方差分别为: [0080] [0081] [0082] [0083] 上式中, [0084] 步骤7、输出多目标的跟踪位置估计轨迹。 [0085] 本发明方法可以应用于具有自动驾驶模式的车辆驾驶系统,本发明方法将关于行驶路段的多目标跟踪信息输出形成多目标跟踪位置估计轨迹,所述多目标跟踪位置估计轨迹信息或数据提供到一个或多个车辆驾驶系统,所述车辆驾驶系统基于本发明提供的信息或数据发出自动驾驶模式的命令。更进一步地,所述多目标跟踪位置估计轨迹信息或数据与来自车辆驾驶系统或其它终端的车辆导航系统的地图数据相结合,提升自动驾驶的用户体验与安全性。 [0086] 本发明涉及的一种基于汽车防撞雷达的道路多目标跟踪方法,采用有向图切换的无迹卡尔曼概率密度假设(DS-UKGMPHD)方法,与现有技术相比,本发明达到的技术效果为: [0087] 1、本发明毋需数据关联、易于实施,适用于非线性、复杂运动状态下的多目标跟踪,与雷达极坐标形式的量测方程相匹配,保证了汽车防撞雷达多目标跟踪的实时性。 [0088] 2、在道路车辆数目与模型覆盖较多的情况下,本发明实时剔除不必要模型,同时兼顾机动目标跟踪实时性和可实现性,提高了跟踪准确度。 [0089] 3、本发明在杂波环境中实时估计多个目标的数目与位置等信息,进一步提高防撞雷达的可靠性,降低虚警率。 [0091] 图1为本发明有向图模型,其中l为模型序号,Φ(0)为默认有向子图。 [0092] 图2为本发明的实现流程图。 [0093] 图3为本发明实施例2中杂波湮没的道路多机环境下各目标真实轨迹及汽车防撞雷达对各目标的量测点迹。 [0094] 图4为本发明实施例2中使用DS-UKGMPHD方法对汽车防撞雷达多目标跟踪图。其中,“○”表示目标运动轨迹的起始端点,“☆”表示载有汽车防撞雷达的自车位置。 [0095] 图5为本发明实施例2中使用IMM-UKGMPHD方法对汽车防撞雷达多目标跟踪图。 [0096] 图6为本发明实施例2中使用DS-UKGMPHD方法,跟踪各目标在x方向随时间的目标位置-时间图。 [0097] 图7为本发明实施例2中使用DS-UKGMPHD方法,跟踪各目标在y方向随时间的目标位置-时间图。 [0098] 图8为本发明实施例2中使用IMM-UKGMPHD方法,跟踪各目标在x方向随时间的目标位置-时间图。 [0099] 图9为本发明实施例2中使用IMM-UKGMPHD方法,跟踪各目标在y方向随时间的目标位置-时间图。 [0100] 图10为本发明实施例2中使用DS-UKGMPHD方法与IMM-UKGMPHD方法目标个数跟踪对比图。 [0101] 图11为本发明实施例2中使用DS-UKGMPHD方法与IMM-UKGMPHD方法跟踪多目标位置-时间对比图。 具体实施方式[0102] 下面通过具体实施例,进一步对本发明的技术方案进行具体说明。应该理解,下面的实施例只是作为具体说明,而不限制本发明的范围,同时本领域的技术人员根据本发明所做的显而易见的改变和修饰也包含在本发明范围之内。 [0103] 实施例1 [0104] 一种基于汽车防撞雷达的道路多目标跟踪方法,具体步骤如下: [0105] 步骤1、获取先验信息,建立第一观测与运动模型。所述先验信息为通过汽车防撞雷达获取的目标状态信息。 [0106] 确定状态方程 量测方程 及其道路目标的多种模型。依据实际情况选择不同模型,状态方程建立为: [0107] 其中, 是目标模型的状态转移矩阵,下标m代表不同目标模型,下标k表示时刻k,下标i表示不同目标。假设噪声服从高斯分布,有 [0108] 汽车防撞雷达观测向量 θk为观测方位角,rk为目标与观测点间距离。 [0109] 步骤2、建立第二观测与运动模型,初始化有向图模型集,建立有向图Φ。 [0110] 2(a)汽车防撞雷达车辆目标限于二维平面运动,目标i在k时刻的非线性量测模型如下: [0111] [0112] 其中, 为高斯分布的量测噪声,xs,k,ys,k为自车观测位置坐标。 [0113] 2(b)建立有向图Φ,及其有向子图Φ(l){βl-1,βl,βl+1}。 [0114] 步骤3、建立自适应有向图切换准则。设k-1时刻子图为Φ(k-1),切换规则如 下: [0115] [0116] 式中,τDS是模型转换的阈值概率, 表示k-1时刻概率β(i±1); [0117] 步骤4、初始化及计算混合参数。 [0118] 4(a)时刻k=0时,每一向量子图Dn中含M个模型。对于当前模型q(q=1,...,M),初始化目标状态 目标协方差矩阵 模型概率 [0119] 4(b)设当前时刻k-1,下一时刻为k。则对于k>0时刻的第i个高斯分量,有[0120] [0121] 其中, 为k-1时刻的高斯分量k>0模型p的出现概率, [0122] 当前有向量子图内模型数M,交互初始PHD函数为: [0123] [0124] 混合权值、状态、协方差分别为: [0125] [0126] [0127] [0128] 上式中, [0129] 步骤5、各目标模型的UKGMPHD滤波。 [0130] 以各高斯分量i在k-1时刻所得混合参数作为输入进行UKGMPHD滤波,得k时刻跟踪权值,状态及协方差矩阵 包括以下过程: [0131] 5(a)对k=0时刻,初始化目标概率密度,有: [0132] [0133] 5(b)预测。对k>0时刻,k-1时刻,高斯分量i总数Jk-1,概率密度假设函数Dk-1(x)为: [0134] 所述预测概率密度假设函数Dk|k-1(x)依据对象的不同分为存活目标、新生目标及衍生目标概率密度假设函数。本发明道路多目标跟踪通常不涉及衍生目标,故概率密度假设函数为:Dk|k-1(x)=DS,k|k-1(x)+γk(x) [0135] 存活目标、新生目标概率密度假设函数分别如下: [0136] [0137] [0138] 使用无迹卡尔曼滤波的思想滤波。依据状态方程,一步预测均值 方差 如下: [0139] [0140] [0141] 式中,l=0,……,2n为无迹滤波的Sigma点,L=2n+1为Sigma点总数。则k时刻预测概率密度假设函数为: [0142] 5(c)更新。 [0143] 对于雷达获取的非线性量测方程,一步预测均值、方差分别为: [0144] [0145] [0146] 结合步骤5(b)结果Dk|k-1(x)及检测概率PD,k,k时刻的更新概率密度假设函数Dk(x)为: [0147] [0148] DD,k(x;z)、估计均值 方差 可表示为: [0149] [0150] 其中, [0151] 高斯分量i的估计均值和方差如下: [0152] [0153] [0154] 上式中,状态增益矩阵、协方差分别为: [0155] [0156] [0157] 5(d)剪枝。 [0158] 对于更新得到的高斯分布进行剪枝,若高斯分布小于修剪阈值τ,直接予以删除。公式为: [0159] 5(e)合并。 [0160] 对于分布接近的若干高斯分布,合并为一个新的分布,即集合: [0161] [0162] 合并后,新的高斯分布为: [0163] [0164] [0165] 5(f)目标数目与状态估计。 [0166] 目标数目估计为融合估计的权值,即: [0167] 取整可得目标数目。按照上述步骤输出合并后状态记作 [0168] 在步骤5完成后,存在模型概率更新步骤;完成模型概率更新后,返回步骤4。所述模型概率更新步骤如下: [0169] k时刻模型(q=1,...,M)下高斯分量i的似然函数为 模型q更新概率为: [0170] [0171] [0172] 步骤6、估计融合。估计融合PHD函数: [0173] [0174] 其中,Jk为高斯分量个数,高斯分量i的权值,均值,协方差分别为: [0175] [0176] [0177] [0178] 上式中, [0179] 步骤7、输出多目标的跟踪位置估计轨迹。 [0180] 实施例2 [0181] 本实施例基于计算机仿真对实施例1进行验证,验证环境为MatlabR2014a。具体实施方案如下: [0182] 1.仿真场景设置: [0183] 在二维平面区域[150m,150m]2内,模拟双向四车道道路环境,自左向右为第一车道至第四车道,以y轴正方向为第一、二车道流向。描述过程中,速度、加速度方向默认与流向一致。假设汽车防撞雷达探测范围在120m内,载有汽车防撞雷达的自车位置(40,0)m,并以25m/s的速度向y轴正方向做匀速运动。x坐标为各目标所处车道位置,y坐标为防撞雷达以自身为静止参考点探测到的目标相对位置。模拟道路目标运动情况如下:防撞雷达在0.1s于(62,120)m处探测到车辆目标1以24.8m/s的初速度与自车相向而行,做加速度 0.21m/s2的匀减速运动,1.2s起以加速度0.08m/s2做匀加速运动,至4s恢复匀速运动;在 2.1s时探测到车辆目标2在(70,106)m处以初速度22.3m/s从第四车道变道至第三车道,角速度为-0.1rad/s,3.6s起以角速度-0.04rad/s左转,驶回第四车道;探测到车辆目标3出现 2 于3.1s在(22,3)m处,以初速度26.7m/s做加速度0.5m/s的加速运动,在4.5s后做加速度为 0.02m/s2的减速运动;车辆目标4在6s时在(18,3)m处被探测到以27.5m/s做匀速运动,6.8s起变道至第二条车道,该过程角速度分别为0.05rad/s,-0.05rad/s,9.3s时驶于自车前方以27.5m/s匀速行驶。 [0184] 2.实施例验证 [0185] 对于状态模型,取目标状态向量 为 各分量分别代表道路目标在x轴方向的坐标、速度,y轴方向的坐标、速度,以及x轴方向、y轴方向加速度。 [0186] 使用五种模型进行跟踪,选取状态模型如下: [0187] [0188] [0190] 设噪声、杂波服从高斯分布,过程噪声协方差矩阵Q=diag([5,2,5,2,0.5,0.5]),量测噪声协方差矩阵R=diag([0.01,0.01])。目标存活概率PS=0.95,检测概率PD=0.99,合并门限为4,允许最大高斯分量个数Jmax为100,剪枝阈值τ为10-5,每次扫描平均杂波数为4。选取最优次模式分配(OSPA)距离为评价指标。状态模型中,取 对应 的角速度ω分别为-1rad/s,-0.05rad/s,0.05rad/s,1rad/s。仿真时长 10s,采样间隔0.1s。 [0191] 由仿真结果,防撞雷达在附图3的杂波道路及非线性条件下进行量测与跟踪。由附图4、附图6、附图7,DS-UKGMPHD算法跟踪位置与真实轨迹基本一致,而附图8、附图9表明IMM-UKGMPHD算法对目标的跟踪存在明显误差。由附图10,在使用5个跟踪模型的情况下,DS-UKGMPHD算法的目标误跟点个数较IMM-UKGMPHD算法更少。由附图11,在目标数目误跟时刻,两算法均出现峰值,而在其余时刻,DS-UKGMPHD算法的OSPA距离误差相对较小。 |