CN111857134B - 一种基于贝叶斯网络的目标障碍车轨迹预测方法 - Google Patents

一种基于贝叶斯网络的目标障碍车轨迹预测方法 Download PDF

Info

Publication number
CN111857134B
CN111857134B CN202010601327.XA CN202010601327A CN111857134B CN 111857134 B CN111857134 B CN 111857134B CN 202010601327 A CN202010601327 A CN 202010601327A CN 111857134 B CN111857134 B CN 111857134B
Authority
CN
China
Prior art keywords
probability
grid
target obstacle
speed
obstacle vehicle
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
Application number
CN202010601327.XA
Other languages
English (en)
Other versions
CN111857134A (zh
Inventor
江洪
蒋潇杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongheng Zhilian (Jiangsu) Technology Co.,Ltd.
Original Assignee
Jiangsu University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangsu University filed Critical Jiangsu University
Priority to CN202010601327.XA priority Critical patent/CN111857134B/zh
Publication of CN111857134A publication Critical patent/CN111857134A/zh
Application granted granted Critical
Publication of CN111857134B publication Critical patent/CN111857134B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/0088Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots characterized by the autonomous decision making process, e.g. artificial intelligence, predefined behaviours
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0231Control of position or course in two dimensions specially adapted to land vehicles using optical position detecting means
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0257Control of position or course in two dimensions specially adapted to land vehicles using a radar
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/29Graphical models, e.g. Bayesian networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Data Mining & Analysis (AREA)
  • Development Economics (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Medical Informatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Traffic Control Systems (AREA)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

本发明公开无人驾驶汽车领域中的一种基于贝叶斯网络的目标障碍车轨迹预测方法,环境模型分别计算出历史路径概率PA、航向角概率PB、目标障碍车与其他障碍物距离影响概率PC、历史速度概率PD、历史避让障碍物速度概率PE和目标障碍车速度概率PF并输入贝叶斯网络模型,贝叶斯网络模型将概率PA、PB、PC、PD、PE、PF作为子节点输入,下一步路径栅格概率、下一步栅格速度概率作为根节点构建贝叶斯网络结构,得到从当前所在栅格到周围八个栅格的不同的后验概率以及各栅格对应的速度的后验概率,比较出后验概率值最大的栅格及速度作为目标障碍车的下一轨迹点,本发明考虑动态环境下目标障碍车与其周边其他障碍物之间距离的影响,预测的未来轨迹更加准确。

Description

一种基于贝叶斯网络的目标障碍车轨迹预测方法
技术领域
本发明涉及无人驾驶汽车领域,具体是对行驶在无人驾驶汽车周围的目标障碍车的行驶轨迹进行预测的方法。
背景技术
在无人驾驶领域中,需要对无人驾驶车进行轨迹规划,从而得到安全平滑的运行轨迹。为了使无人驾驶车能够更准确的避开周围的障碍物,通常需要对障碍物未来的轨迹进行预测。一般来说,通常是根据障碍物在当前以及当前之前预设时段内的运动状态,预测障碍物在未来多个预设时刻的轨迹点,从而得到障碍物未来的轨迹。但是,通过这种方式预测的障碍物未来的轨迹准确度较差,使得无人驾驶车轨迹规划的结果缺乏准确性。例如,在中国专利申请号为201910034446.9的文献中公布一种基于障碍物预测与MPC算法的无人驾驶轨迹规划方法,其通过车载传感器设备及相应的图像识别算法得到周围车辆的运动速度与方向,其假设在预测时域内,周围车辆沿当前车道按照当前速度行驶,得到周围车辆运动轨迹,保证了其后期轨迹规划的安全性,但是其假定在预测时域内,按照当前车道以及当前速度行驶,忽略了周围其他障碍物对其行驶的影响,与实际动态的行车环境不符。
发明内容
针对现有技术中存在的不足,本发明提出一种基于贝叶斯网络的目标障碍车轨迹预测方法,构建基于贝叶斯网络的轨迹预测模型,通过考虑相似环境下历史轨迹、航向角、目标障碍车与周围其他障碍物距离、避让障碍物速度、目标障碍车速度等多种影响因素,提高预测轨迹的准确性。
本发明是通过以下技术手段实现上述技术目的的:信息采集模块采集当前地图以及周围道路环境信息S、当前地点历史轨迹数据T、目标障碍车当前位置Ps、目标障碍车速度v0、目标障碍车与其他障碍物之间的相对距离p、其他障碍物速度v1和目标障碍车的航向角B信息并输入到环境模型中,历史轨迹数据T包括历史路径A、历史速度D和历史避让障碍物速度E,还包括以下步骤:
步骤(1)环境模型根据当前地图以及周围道路环境信息S将整个地图空间划分成大小相同的栅格,构建栅格地图,分别计算出历史路径概率PA、航向角概率PB、目标障碍车与其他障碍物距离影响概率PC、历史速度概率PD、历史避让障碍物速度概率PE和目标障碍车速度概率PF并输入贝叶斯网络模型;
步骤(2)贝叶斯网络模型将所述的概率PA、PB、PC、PD、PE、PF作为子节点输入,下一步路径栅格概率G、下一步栅格速度概率H作为根节点,构建贝叶斯网络结构,得到各节点的条件概率分布以及下一步路径栅格概率G、下一步栅格速度概率H的先验概率,运用贝叶斯公式得到从当前所在栅格Mk到周围M1到M8号栅格的不同的后验概率以及各栅格对应的速度的后验概率;
步骤(3)贝叶斯网络模型比较出后验概率值最大的栅格及速度作为目标障碍车的下一轨迹点,所有轨迹点连接起来得到预测的目标障碍车轨迹。
步骤(1)中,所述的历史路径概率PA
Figure BDA0002559113170000021
n(k→u)为当前地点相似历史路径A中包括栅格号串{Mk,Mu}的路径数目,nk为所有包含栅格Mk号的路径数目,Mk为目标障碍车当前所在栅格,Mu为目标障碍车即将到达的下一栅格,u∈{1,2,3,…,8};
选取九个相邻的栅格组成九宫格,以九宫格中心为圆心画一个内切圆,与九宫格内线相交八个点,这八个交点分别连接圆心,内切圆被分成八个扇形,分别对应八个栅格号,由航向角B得到目标障碍车的车头朝向和车辆中心线,得到车辆前轮转角范围∠C和前轮转角范围与不同扇形重叠部分角度∠A,则所述的航向角概率PB
Figure BDA0002559113170000022
所述的目标障碍车与其他障碍物距离影响概率PC
Figure BDA0002559113170000023
p为目标障碍车与周围其他障碍物的相对距离,p0为障碍物影响范围,v为目标障碍车与其他障碍物的相对速度,f为地面摩擦系数,a为安全系数,半径为p0+a的圆经过的栅格,其概率值最大,为55%;
所述的历史速度概率
Figure BDA0002559113170000024
nvi为当前地点环境下统计各区间速度出现的频次,nv为历史速度D总数;
所述的历史避让障碍物速度概率
Figure BDA0002559113170000025
navi为避让障碍物时各区间速度出现的频次,nav历史避让障碍物速度E总数;
当目标障碍车速度v0属于划分的速度区间时,则所述的目标障碍车速度概率PF
Figure BDA0002559113170000031
t为时间,其余区间目标障碍车速度概率PF为10%。
本发明的有益效果是:
(1)本发明将朴素贝叶斯网络模型运用到轨迹预测中,不仅考虑了历史轨迹的影响,同时还考虑了当前的航向角及动态环境下目标障碍车与其周边其他障碍物距离等影响,使影响因素不再单一,避免出现特殊情况而不能准确预测。
(2)本发明通过考虑动态环境下目标障碍车与其周边其他障碍物之间距离的影响,并提出对应的叶节点概率计算公式,更加符合实际动态环境,并对目标障碍车的每一个路径点的速度进行预测,因此,预测得到的目标障碍车未来轨迹更加准确。
(3)本发明设立滚动周期进行滚动预测,使预测出的路径能够适应复杂多变的环境。
附图说明
图1为实现本发明一种基于贝叶斯网络的目标障碍车轨迹预测方法的预测系统的硬件图;
图2为目标障碍车、本车和其他障碍物关系示意图;
图3为本发明目标障碍车轨迹预测方法的流程图;
图4为本发明路径预测网络结构图;
图5为本发明速度预测网络结构图;
图6为本发明运用模型预测路径栅格号和速度区间的流程图;
图7为本发明中贝叶斯网络预测方法原理图;
图8为本发明设计的相似历史路径栅格示意图;
图9为本发明设计的航向角示意图;
图10为本发明设计的目标障碍物与其他障碍物距离影响示意图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步的说明,但本发明的保护范围并不限于此。
本发明首先是要构建如图1所示的预测系统,该预测系统基于朴素贝叶斯网络建立的贝叶斯网络模型,通过贝叶斯网络模型对目标障碍车轨迹进行预测估计。该预测系统主要由依次串联的信息采集模块、信息处理模块和决策模块组成。其中,信息采集模块包括雷达、CCD相机和历史轨迹数据库及地图下载部分,信息处理模块包括依次串联的环境模型和贝叶斯网络模型。
雷达部分负责采集目标障碍车当前位置Ps、目标障碍车速度v0、目标障碍车与其他障碍物之间的相对距离p、其他障碍物速度v1,并将其发给信息处理模块中的环境模型。
CCD相机部分负责采集目标障碍车航向角B,并将其发给信息处理模块中的环境模型。
历史轨迹数据库及地图下载部分用于获取当前地点以往经过车辆的历史轨迹数据T和当前地图以及周围道路环境信息S,并且将其发给信息处理模块中的环境模型。其中,当前地点以往经过车辆的历史轨迹数据T主要包括当前地点相似历史路径A、当前地点历史速度D、当前地点历史避让障碍物速度E。
信息处理模块中的环境模型根据接收到的当前地图以及周围道路环境信息S进行环境建模,然后结合目标障碍车当前位置Ps、目标障碍车与其他障碍物之间的相对距离p、目标障碍车航向角B、目标障碍车速度v0、其他障碍物速度v1以及当前地点历史轨迹数据T计算出历史路径概率PA、航向角概率PB、目标障碍车与其他障碍物距离影响概率PC、历史速度概率PD、历史避让障碍物速度概率PE和目标障碍车速度概率PF,并将计算出的PA、PB、PC、PD、PE、PF,输入贝叶斯网络模型,作为贝叶斯网络模型各叶节点概率,最后,由贝叶斯网络模型根据各叶节点概率预测出目标障碍车轨迹l并发送给决策模块。决策模块根据信息处理模块发送的目标障碍车轨迹l信息进行下一步决策处理。
如图2所示,本发明的研究对象分为两类:无人驾驶汽车和障碍物。图2中的黑色圆代表无人驾驶汽车,无人驾驶汽车称为本车,黑色方块代表目标障碍车,灰色方块代表其他障碍物。其中,障碍物又可分为目标障碍车和其他障碍物,目标障碍车为本发明轨迹预测方法的预测对象,其他障碍物为相对于目标障碍车来说的路径范围内的障碍物,包括无人驾驶汽车相对于目标障碍车也属于其他障碍物,p是目标障碍车与其他障碍物之间的相对距离。图1中的预测系统安装在本车上,用来采集和处理当前目标障碍车周围环境和目标障碍车的状态信息,然后预测出目标障碍车未来n秒内的目标障碍车轨迹l。
本发明所述的一种基于贝叶斯网络的目标障碍车轨迹预测方法在进行预测时,分为预测模型构建和目标障碍车未来轨迹预测两个阶段,流程图如图3所示。首先,信息采集模块将当前地图以及周围道路环境信息S、当前地点历史轨迹数据T、目标障碍车当前位置Ps、目标障碍车速度v0、目标障碍车与其他障碍物之间的相对距离p、其他障碍物速度v1和目标障碍车的航向角B信息发送给信息处理模块,其次,信息处理模块中环境模型根据当前地图以及周围道路环境信息S进行环境建模,将整个地图空间划分成大小相同的栅格,形成一个栅格地图,栅格长度为L,在同一个栅格内的对象都被看作是一个质点,则栅格地图上的一条轨迹就可以由一连串栅格号表示,如图7所示,然后信息处理模块通过信息采集模块发送的信息计算出目标障碍车周围相邻八个栅格被选中的后验概率值以及目标障碍车速度的概率值,比较出后验概率值最大的栅格及速度作为目标障碍车的下一轨迹点,所有轨迹点连接起来得到一条目标障碍车轨迹l,最后信息处理模块将目标障碍车轨迹l发送给决策模块。其中信息采集模块每经过一个滚动周期将重新采集信息传输给信息处理模块重新计算。具体步骤如下:
第一阶段,预测模型的构建
步骤1:信息采集模块将雷达采集到的目标障碍车当前位置Ps、目标障碍车速度v0、目标障碍车与其他障碍物之间相对距离p、其他障碍物速度v1、CCD相机采集的目标障碍车航向角B以及下载的当前地点历史轨迹数据T、当前地图以及周围道路环境信息S发送给信息处理模块,信息处理模块中的环境模型可计算出目标障碍车与其他障碍物相对速度v,v=v0-v1。历史轨迹数据T包括历史路径A、历史速度D和历史避让障碍物速度E。
步骤2:信息处理模块中的环境模型根据信息采集模块传回的信息构建栅格地图,并根据目标障碍车当前位置Ps、目标障碍车与其他障碍物之间的相对距离p、目标障碍车的航向角B、目标障碍车速度v0以及当前地点历史轨迹数据T分别设计历史路径概率PA、航向角概率PB、目标障碍车与其他障碍物距离影响概率PC、历史速度概率PD、历史避让障碍物速度概率PE和目标障碍车速度概率PF计算公式。
参见图7和图8,将目标障碍车周围8个栅格从1~8标号,目标障碍车当前所在栅格称为栅格Mk,则目标障碍车即将到达的下一栅格为Mu,u∈{1,2,3,…,8},并将速度划分成连续的12个区间,速度区间间隔10km/h,分别为[0,10]、[10,20]、[20,30]…[110,120],单位为km/h。栅格与速度区间划分得越细,轨迹之间区分的越明显。
1.计算历史路径概率PA:如图8所示,为了充分利用历史轨迹T的信息,基于每个栅格的状态值来设计历史路径转移概率公式。目标障碍车从栅格Mk中行驶到栅格Mu中的历史路径概率PA的值等于当前地点相似历史路径A中包括栅格号串{Mk,Mu}的路径数目除以所有包含栅格Mk号的路径数目。因此,对于栅格图中的每对相邻栅格,都可以预计算出转移概率,公式为:
Figure BDA0002559113170000051
其中n(k→u)为当前地点相似历史路径A中包括栅格号串{Mk,Mu}的路径数目,nk为所有包含栅格Mk号的路径数目。
2.计算航向角概率PB:如图9所示,选取九个相邻的栅格组成九宫格,以九宫格中心为圆心,画一个内切圆,与九宫格内线相交八个点,这八个交点分别连接圆心,内切圆被分成八个扇形,分别对应八个栅格号。由目标障碍车航向角B得到目标障碍车的车头朝向,图9中栅格M5中的粗虚线表示车辆中心线,∠C为车辆前轮转角范围,∠A为前轮转角范围与不同扇形重叠部分角度,其中在一个滚动周期Δt内认为目标障碍车行驶方向不变。目标障碍车从Mk中行驶到Mu中的航向角概率PB的值等于前轮转角范围与不同扇形重叠部分角度和车辆前轮转角之间的比值:
Figure BDA0002559113170000061
其中,∠A为前轮转角与不同扇形重叠部分角度,∠C为车辆前轮转角范围。
3.计算目标障碍车与其他障碍物距离影响概率PC:如图10所示,网格颜色越深代表其概率值越低,基于环境建模,将整个道路栅格概率定义为50%,当周围出现其他障碍物时,越靠近障碍物周围的栅格概率值越低,距离越远的概率值越高;半径为p0+a的圆经过的栅格,其概率值最大,为55%;当距离超过p0+a时,则恢复原道路栅格概率值50%。其中p0为障碍物影响范围,随着目标障碍车与其他障碍物之间相对速度变化而变化。目标障碍车从栅格Mk中行驶到栅格Mu中的目标障碍车与其他障碍物距离影响概率用PC表示:
Figure BDA0002559113170000062
Figure BDA0002559113170000063
其中,p为目标障碍车与周围其他障碍物的相对距离,p0为障碍物影响范围,v为目标障碍车与其他障碍物的相对速度,f为地面摩擦系数,a为安全系数。
4.计算历史速度概率PD:将当前地点环境下统计以往经过车辆各区间速度出现的频次记为nvi,所有统计的当前地点历史速度D总数记为nv,目标障碍车从栅格Mk中行驶到栅格Mu中的历史速度概率PD的值等于各区间速度出现的频次与历史速度总数的比值:
Figure BDA0002559113170000071
其中nvi为当前地点环境下统计各区间速度出现的频次,nv为所有统计的历史速度D总数。
5.计算历史避让障碍物速度概率PE:统计避让障碍物时各区间速度出现的频次记为navi,所有统计的当前地点历史避让速度E总数记为nav,目标障碍车从栅格Mk中行驶到栅格Mu中的避让障碍物速度概率PE的值等于各区间速度出现的频次与行驶速度总数的比值:
Figure BDA0002559113170000072
其中navi为避让障碍物时各区间速度出现的频次,nav所有统计的历史避让速度E总数。
6.计算目标障碍车速度概率PF:当雷达检测到目标障碍车速度v0属于划分的某个速度区间时,其区间速度概率值用Pt表示,其余区间车速概率值记为10%,一个预测周期内传感器检测到的目标障碍车速度v0随着预测时间的推移而保持不变的概率Pt为:
Figure BDA0002559113170000073
其中t为时间。
则目标障碍车从Mk中行驶到Mu中的目标障碍车速度概率PF为:
Figure BDA0002559113170000074
步骤3:贝叶斯网络模型的构建包括可视化的体现节点相关性的有向无环图以及量化的各个节点的条件概率。现将目标障碍物车下一步即将行驶到栅格Mu的概率称为下一步路径栅格概率G,行驶到下一步路径栅格时的目标障碍物车的速度在某一速度区间概率称为下一步栅格速度概率H。如图4和图5所示,贝叶斯网络模型将环境模型输出的历史路径概率PA、航向角概率PB、目标障碍车与其他障碍物距离影响概率PC、历史速度概率PD、历史避让障碍物速度概率PE、目标障碍车速度概率PF作为子节点输入,下一步路径栅格概率G、下一步栅格速度概率H作为根节点,构建贝叶斯网络结构,得到各网络节点的条件概率分布,具体包括当前地点相似历史路径概率PA和下一步路径栅格概率G的条件概率P(Au|Gu),目标障碍车航向角概率PB和下一步路径栅格概率G的条件概率P(Bu|Gu),目标障碍车与其他障碍物距离影响概率PC和下一步路径栅格概率G的条件概率P(Cu|Gu),相似环境历史速度概率PD和下一步栅格速度概率H的条件概率P(Dj|Hj),当前地点历史避让障碍物速度概率PE和下一步栅格速度概率H的条件概率P(Ej|Hj),目标障碍车速度概率PF和下一步栅格速度概率H的条件概率P(Fj|Hj),下一步路径栅格概率G的先验概率分布Pf(Gu),下一步栅格速度概率H的先验概率分布Pf(Hj),其中,下一步路径栅格概率G包含8种结果{G1,G2,G3,G4,G5,G6,G7,G8},则先验概率Pf(G)分别为Pf(G1),Pf(G2)…Pf(G8),下一步栅格速度概率H包含{H1,H2,H3,H4,H5,H6,H7,H8,H9,H10,H11,H12},分别对应12个速度区间,分别为Pf(H1),Pf(H2)…Pf(H12)。
运用贝叶斯公式可以得到从当前所在栅格Mk到周围M1到M8号栅格的不同的后验概率P(Gu|Au,Bu,Cu)以及各栅格对应的速度的后验概率P(Hj|Dj,Ej,Fj):
(1)
Figure BDA0002559113170000081
(2)
Figure BDA0002559113170000082
其中u∈(1,2,3…8)表示当前所在栅格周围1号栅格到8号栅格,j∈(1,2,3…12)表示12个不同速度区间。
若max{P(G1|A1,B1,C1),P(G2|A2,B2,C2)…P(G8|A8,B8,C8)}=P(Gu|Au,Bu,Cu),则选择栅格Mu为目标障碍车的下一步栅格,将每步预测出栅格号连接起来,则可得到未来n秒内的一系列路径点。
若max{P(H1|D1,E1,F1),P(H2|D2,E2,F2)…P(H12|D12,E12,F12)}=P(Hj|Dj,Ej,Fj),则选择vj为下一步路径栅格上目标障碍车的速度vj,然后将路径点与对应的速度结合起来形成目标障碍车轨迹点,如图6所示,其中,信息采集模块每经过一个滚动周期Δt将重新搜集信息并传输给信息处理模块计算后验概率,更新轨迹点。如此构建成贝叶斯网络模型,即预测模型。
第二阶段:基于预测模型对目标障碍车未来轨迹进行预测。
步骤1:信息采集模块通过雷达、CCD相机采集目标障碍车周围以及目标障碍车当前的状态信息,并下载历史信息和地图,具体包括当前道路地图及周围道路信息S、当前地点以往经过车辆的历史轨迹数据T、目标障碍车的当前位置Ps、目标障碍车与其他障碍物之间的相对距离p、目标障碍车的航向角B、目标障碍车速度v0、其他障碍物速度v1等信息。
步骤2:先将信息采集模块中获取的参数信息发送给信息处理模块作为环境模型的输入,目标障碍车的当前位置Ps为目标障碍车的路径起点,根据第一阶段中步骤2的概率计算公式分别得到各叶节点概率,再由贝叶斯网络模型得到各条件概率,并运用贝叶斯定理分别计算出目标障碍车路径栅格及速度后验概率分布,选定下一路径栅格和速度,其中信息采集模块每经过一个滚动周期Δt将重新搜集信息并传输给信息处理模块计算后验概率,更新轨迹点,最终得到未来n秒的目标障碍车轨迹l。
步骤3:信息处理模块将预测出的目标障碍车轨迹l发送到决策模块,进行下一步决策处理。
所述实施例为本发明的优选的实施方式,但本发明并不限于上述实施方式,在不背离本发明的实质内容的情况下,本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。

Claims (4)

1.一种基于贝叶斯网络的目标障碍车轨迹预测方法,信息采集模块采集当前地图以及周围道路环境信息S、当前地点历史轨迹数据T、目标障碍车当前位置Ps、目标障碍车速度v0、目标障碍车与其他障碍物之间的相对距离p、其他障碍物速度v1和目标障碍车的航向角B信息并输入到环境模型中,历史轨迹数据T包括历史路径A、历史速度D和历史避让障碍物速度E,其特征是包括以下步骤:
步骤(1)环境模型根据当前地图以及周围道路环境信息S将整个地图空间划分成大小相同的栅格,构建栅格地图,分别计算出历史路径概率PA、航向角概率PB、目标障碍车与其他障碍物距离影响概率PC、历史速度概率PD、历史避让障碍物速度概率PE和目标障碍车速度概率PF并输入贝叶斯网络模型;
所述的历史路径概率PA
Figure FDA0003668566840000011
n(k→u)为当前地点相似历史路径A中包括栅格号串{Mk,Mu}的路径数目,nk为所有包含栅格Mk号的路径数目,Mk为目标障碍车当前所在栅格,Mu为目标障碍车即将到达的下一栅格,u∈{1,2,3,…,8};
选取九个相邻的栅格组成九宫格,以九宫格中心为圆心画一个内切圆,与九宫格内线相交八个点,这八个交点分别连接圆心,内切圆被分成八个扇形,分别对应八个栅格号,由航向角B得到目标障碍车的车头朝向和车辆中心线,得到车辆前轮转角范围∠C和前轮转角范围与不同扇形重叠部分角度∠A,则所述的航向角概率PB
Figure FDA0003668566840000012
所述的目标障碍车与其他障碍物距离影响概率PC
Figure FDA0003668566840000013
p为目标障碍车与周围其他障碍物的相对距离,p0为障碍物影响范围,v为目标障碍车与其他障碍物的相对速度,f为地面摩擦系数,a为安全系数,半径为p0+a的圆经过的栅格,其概率值最大,为55%;
所述的历史速度概率
Figure FDA0003668566840000014
nvi为当前地点环境下统计各区间速度出现的频次,nv为历史速度D总数;
所述的历史避让障碍物速度概率
Figure FDA0003668566840000021
navi为避让障碍物时各区间速度出现的频次,nav历史避让障碍物速度E总数;
当目标障碍车速度v0属于划分的速度区间时,则所述的目标障碍车速度概率PF
Figure FDA0003668566840000022
t为时间,其余区间目标障碍车速度概率PF为10%;
步骤(2)贝叶斯网络模型将所述的概率PA、PB、PC、PD、PE、PF作为子节点输入,下一步路径栅格概率G、下一步栅格速度概率H作为根节点,构建贝叶斯网络结构,得到各节点的条件概率分布以及下一步路径栅格概率G、下一步栅格速度概率H的先验概率,运用贝叶斯公式得到从当前所在栅格Mk到周围M1到M8号八个栅格的不同的后验概率以及各栅格对应的速度的后验概率;
当前所在栅格Mk到周围M1到M8号栅格的不同的后验概率
Figure FDA0003668566840000023
各栅格对应的速度的后验概率
Figure FDA0003668566840000024
u∈(1,2,3…8)表示当前所在栅格周围1号栅格到8号栅格,j∈(1,2,3…12)表示12个不同速度区间,当前地点相似历史路径概率PA和下一步路径栅格概率G的条件概率为P(Au|Gu),目标障碍车航向角概率PB和下一步路径栅格概率G的条件概率为P(Bu|Gu),目标障碍车与其他障碍物距离影响概率PC和下一步路径栅格概率G的条件概率为P(Cu|Gu),相似环境历史速度概率PD和下一步栅格速度概率H的条件概率为P(Dj|Hj),当前地点历史避让障碍物速度概率PE和下一步栅格速度概率H的条件概率为P(Ej|Hj),目标障碍车速度概率PF和下一步栅格速度概率H的条件概率为P(Fj|Hj),下一步路径栅格概率G的先验概率分布为Pf(Gu),下一步栅格速度概率H的先验概率分布为Pf(Hj);
步骤(3)贝叶斯网络模型比较出后验概率值最大的栅格及速度作为目标障碍车的下一轨迹点,所有轨迹点连接起来得到预测的目标障碍车轨迹。
2.根据权利要求1所述的一种基于贝叶斯网络的目标障碍车轨迹预测方法,其特征是:若max{P(G1|A1,B1,C1),P(G2|A2,B2,C2)…P(G8|A8,B8,C8)}=P(Gu|Au,Bu,Cu),则选择栅格Mu为目标障碍车的下一步栅格,将每步预测出栅格号连接起来得到未来n秒内的一系列路径点,若max{P(H1|D1,E1,F1),P(H2|D2,E2,F2)…P(H12|D12,E12,F12)}=P(Hj|Dj,Ej,Fj),则选择vj为下一步路径栅格上目标障碍车的速度,将一系列路径点与对应的速度结合起来形成目标障碍车轨迹点。
3.根据权利要求1所述的一种基于贝叶斯网络的目标障碍车轨迹预测方法,其特征是:信息采集模块每经过一个滚动周期△t重新采集信息,贝叶斯网络模型更新目标障碍车轨迹点。
4.根据权利要求1所述的一种基于贝叶斯网络的目标障碍车轨迹预测方法,其特征是:贝叶斯网络模型将预测出的目标障碍车轨迹发送给决策模块进行下一步决策处理。
CN202010601327.XA 2020-06-29 2020-06-29 一种基于贝叶斯网络的目标障碍车轨迹预测方法 Active CN111857134B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010601327.XA CN111857134B (zh) 2020-06-29 2020-06-29 一种基于贝叶斯网络的目标障碍车轨迹预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010601327.XA CN111857134B (zh) 2020-06-29 2020-06-29 一种基于贝叶斯网络的目标障碍车轨迹预测方法

Publications (2)

Publication Number Publication Date
CN111857134A CN111857134A (zh) 2020-10-30
CN111857134B true CN111857134B (zh) 2022-09-16

Family

ID=72988730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010601327.XA Active CN111857134B (zh) 2020-06-29 2020-06-29 一种基于贝叶斯网络的目标障碍车轨迹预测方法

Country Status (1)

Country Link
CN (1) CN111857134B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112529110B (zh) * 2020-12-29 2023-04-07 中国科学院自动化研究所 对手策略反演方法、系统、装置
CN112327888B (zh) * 2021-01-07 2021-03-30 中智行科技有限公司 路径规划方法、装置、电子设备和存储介质
CN113018866A (zh) * 2021-03-31 2021-06-25 完美世界(北京)软件科技发展有限公司 地图资源的加载方法和装置、存储介质、电子装置
CN113359763A (zh) * 2021-07-05 2021-09-07 广州蓝胖子移动科技有限公司 轮式机器人及动态障碍物轨迹预测方法、装置、存储介质
CN113642644B (zh) * 2021-08-13 2024-05-10 北京赛目科技有限公司 车辆环境等级的确定方法、装置、电子设备及存储介质
CN114670870A (zh) * 2022-03-18 2022-06-28 北京智行者科技有限公司 障碍物slt空间风险场环境建模方法及其装置、相关产品
CN116152782A (zh) * 2023-04-18 2023-05-23 苏州魔视智能科技有限公司 一种障碍物轨迹预测方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106950956A (zh) * 2017-03-22 2017-07-14 合肥工业大学 融合运动学模型和行为认知模型的行车轨迹预测系统
CN109885066A (zh) * 2019-03-26 2019-06-14 北京经纬恒润科技有限公司 一种运动轨迹预测方法及装置
CN110146100A (zh) * 2018-02-13 2019-08-20 华为技术有限公司 轨迹预测方法、装置及存储介质
CN110989636A (zh) * 2020-02-26 2020-04-10 北京三快在线科技有限公司 一种障碍物的轨迹预测方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106950956A (zh) * 2017-03-22 2017-07-14 合肥工业大学 融合运动学模型和行为认知模型的行车轨迹预测系统
CN110146100A (zh) * 2018-02-13 2019-08-20 华为技术有限公司 轨迹预测方法、装置及存储介质
CN109885066A (zh) * 2019-03-26 2019-06-14 北京经纬恒润科技有限公司 一种运动轨迹预测方法及装置
CN110989636A (zh) * 2020-02-26 2020-04-10 北京三快在线科技有限公司 一种障碍物的轨迹预测方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种基于高斯混合模型的轨迹预测算法;乔少杰 等;《软件学报》;20150202;全文 *

Also Published As

Publication number Publication date
CN111857134A (zh) 2020-10-30

Similar Documents

Publication Publication Date Title
CN111857134B (zh) 一种基于贝叶斯网络的目标障碍车轨迹预测方法
US11726477B2 (en) Methods and systems for trajectory forecasting with recurrent neural networks using inertial behavioral rollout
JP7140849B2 (ja) 確率的オブジェクト追跡及び予測フレームワーク
CN107063280B (zh) 一种基于控制采样的智能车辆路径规划系统及方法
US20200307561A1 (en) System and method for radar cross traffic tracking and maneuver risk estimation
US20190145765A1 (en) Three Dimensional Object Detection
EP3745158B1 (en) Methods and systems for computer-based determining of presence of dynamic objects
US10732639B2 (en) Method and apparatus for automatically generated curriculum sequence based reinforcement learning for autonomous vehicles
CN112840350A (zh) 自动驾驶车辆规划和预测
US20200148215A1 (en) Method and system for generating predicted occupancy grid maps
US11827214B2 (en) Machine-learning based system for path and/or motion planning and method of training the same
CN116249947A (zh) 预测运动规划系统及方法
CN112534483B (zh) 预测车辆驶出口的方法和装置
EP1290467A1 (en) Track map generator
WO2019134110A1 (en) Autonomous driving methods and systems
Maddern et al. Continuous appearance-based trajectory SLAM
US20220355825A1 (en) Predicting agent trajectories
US20230347979A1 (en) Methods and processors for controlling steering of self-driving car
US20230150510A1 (en) Methods and systems for tracking a mover's lane over time
US20230202473A1 (en) Calculating vehicle speed for a road curve
Jeong Predictive lane change decision making using bidirectional long shot-term memory for autonomous driving on highways
Schüle et al. Augmenting night vision video images with longer distance road course information
Kuhnt et al. Lane-precise localization of intelligent vehicles using the surrounding object constellation
CN114212107A (zh) 一种基于价值图的智能网联驾驶规划方法
US20240132112A1 (en) Path-based trajectory prediction

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
TR01 Transfer of patent right

Effective date of registration: 20230801

Address after: 230000 Room 203, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee after: Hefei Jiuzhou Longteng scientific and technological achievement transformation Co.,Ltd.

Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301

Patentee before: JIANGSU University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230818

Address after: 407-A1, Jiangdong Building, No. 101 Jiangdong North Road, Gulou District, Nanjing City, Jiangsu Province, 210000

Patentee after: Zhongheng Zhilian (Jiangsu) Technology Co.,Ltd.

Address before: 230000 Room 203, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee before: Hefei Jiuzhou Longteng scientific and technological achievement transformation Co.,Ltd.

TR01 Transfer of patent right