CN109990796B - 基于双向扩展随机树的智能车路径规划方法 - Google Patents

基于双向扩展随机树的智能车路径规划方法 Download PDF

Info

Publication number
CN109990796B
CN109990796B CN201910328124.5A CN201910328124A CN109990796B CN 109990796 B CN109990796 B CN 109990796B CN 201910328124 A CN201910328124 A CN 201910328124A CN 109990796 B CN109990796 B CN 109990796B
Authority
CN
China
Prior art keywords
node
new2
new1
nodes
point
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
CN201910328124.5A
Other languages
English (en)
Other versions
CN109990796A (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.)
Chengdu University of Information Technology
Original Assignee
Chengdu University of Information Technology
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 Chengdu University of Information Technology filed Critical Chengdu University of Information Technology
Priority to CN201910328124.5A priority Critical patent/CN109990796B/zh
Publication of CN109990796A publication Critical patent/CN109990796A/zh
Application granted granted Critical
Publication of CN109990796B publication Critical patent/CN109990796B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/26Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/34Route searching; Route guidance
    • G01C21/3446Details of route searching algorithms, e.g. Dijkstra, A*, arc-flags, using precalculated routes

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Steering Control In Accordance With Driving Conditions (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种基于双向扩展随机树的智能车路径规划方法,其包括生成曲线,当曲线上有障碍物时构建起点相反的两棵扩展树;在地图上生成随机点,并查找扩展树距离随机点最近的节点;判断随机点与节点的连线上是否有障碍物;若有,则重新生成随机点,若无,在随机点与节点的连线上取一个新节点;判断最近的节点与新节点之间的连线上是否有障碍物,若有,则重新生成随机点,否则根据障碍物安全距离约束更新节点;之后根据智能车辆最大转向约束及节点,选取节点的父节点和子节点;判断两棵扩展树之间的距离是否小于设定阈值,若大于,则重新生成随机点,否则,根据两扩展树节点反向选取节点到起始点和目标点的路径构成智能车规划路径。

Description

基于双向扩展随机树的智能车路径规划方法
技术领域
本发明涉及物体运动路径的规划,具体涉及一种基于双向扩展随机树的智能车路径规划方法。
背景技术
在2001年Kuffner和LaValle提出并行生成随机树(RRT-Connect)的方法,以提高该算法的收敛速度。在2010年Karaman和Frazzoli提出RRT*算法,用于解决RRT算法产生的路径并非概率最优解的问题,但其运行效率较低。在2013 年,M.Jordan和A.Perez提出B-RRT*,用RRT*算法进行起始点和目标点的双向扩展,以减少RRT*算法所执行的时间,其应用于移动机器人中,B-RRT*虽然能够缩短一部分执行时间,但是其仍存在如下缺点:
(1)在运行过程中,使用全局均匀随机采样,使得算法运行时间增加,收敛的速度变慢;
(2)采用最近节点选择算法进行节点选择,容易导致其在复杂场景时所规划的路径失败;
(3)所规划的路径没有考虑车辆运动学约束,使得规划的路径不能运用于智能车的路径规划中。
发明内容
针对现有技术中的上述不足,本发明提供的基于双向扩展随机树的智能车路径规划方法,其规划的路径安全合理且更适用于车辆运动学特性。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种基于双向扩展随机树的智能车路径规划方法,其包括:
S1、采用Reeds-Shepp曲线连接智能车的起始点和目标点;
S2、当生成的曲线上存在障碍物时,从起始点向目标点构建扩展树tree1,从目标点向起始点构建扩展树tree2;
S3、在包含智能车起始点和目标点的地图上生成随机点qrand1/qrand2,并查找扩展树tree1/tree2距离随机点qrand1/qrand2最近的节点qnearrest1/qnearrest2
S4、判断随机点qrand1/qrand2与节点qnearrest1/qnearrest2的连线上是否有障碍物;若有,则返回步骤S3,否则进入步骤S5;
S5、在随机点qrand1/qrand2与节点qnearrest1/qnearrest2的连线上取位于一个扩展步长长度处的节点qnew1/qnew2
S6、判断节点qnearrest1/qnearrest2与节点qnew1/qnew2的连线上是否有障碍物;若有,则返回步骤S3,否则,进入步骤S7;
S7、选取以节点qnew1/qnew2所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新节点qnew1/qnew2
S8、根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/qnew2的父节点和子节点;
S9、判断扩展树tree1和扩展树tree2之间的距离是否小于设定阈值,若大于,则返回步骤S3,否则,进入步骤S10;
S10、根据扩展树tree1和扩展树tree2相遇的节点qa和qb,反向选取节点 qa和qb分别到起始点和目标点的路径构成智能车规划路径。
进一步地,在包含智能车起始点和目标点的地图上生成随机点qrand1/qrand2进一步包括:
S31、采用随机函数生成一个采样概率,并判断采样概率是否大于目标偏向概率,若大于,则进入步骤S32,否则,进入步骤S33;
S32、采用目标点/起始点作为扩展树tree1/扩展树tree2的随机点qrand1/随机点qrand2
S33、在扩展树tree1/扩展树tree2中寻找与目标点/起始点最近的节点qnr1/ 节点qnr2,以节点qnr1/节点qnr2为中心,以扩展步长为边长的方形区域,并判断方形区域内是否有障碍物;
S34、若没有障碍物,则采用节点qnr1/节点qnr2建立滑动窗口,并在滑动窗口中生成qrand1/随机点qrand2
S35、若有障碍物,则返回步骤S31。
进一步地,选取以节点qnew1/qnew2所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新节点qnew1/qnew2进一步包括:
S71、形成一块以节点qnew1/qnew2为中心,扩展步长为1/2边长的方形区域 A1/A2,并同时初始化规划路径与障碍物之间的最小距离dmin为无穷大;
S72、选取方形区域A1/A2中未选取过的节点poses1/poses2,并判断节点poses1/poses2是否为障碍物;若是,执行步骤S72,否则进入步骤S73;
S73、以方形区域A1/A2中的每个节点poses1/poses2为中心,以扩展步长为边长,形成区域B1/B2;
S74、选取区域B1/B2中未选取过的节点poses11/poses21,并判断节点 poses11/poses21是否为障碍物;若是,返回步骤S72,否则进入步骤S75;
S75、计算节点poses11/poses21到目标点的距离d,并判断距离d是否小于最小距离dmin;若是,进入步骤S76,否则返回步骤S72;
S76、采用距离d更新最小距离dmin,采用距离d对应的节点poses1/poses2 更新节点qnew1/qnew2
S77、判断方形区域A1/A2,是否存在未选取过得节点;若有,则执行S72;否则,输出节点qnew1/qnew2
进一步地,根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/ qnew2的父节点进一步包括:
S811、选取以节点qnew1/qnew2为中心,预设半径R为半径所覆盖区域内的所有节点为待选父节点qparent1/qparent2
S812、选取未选取的待选父节点qparent1/qparent2,计算待选父节点qparent1/ qparent2、待选父节点qparent1/qparent2的父节点和节点qnew1/qnew2所形成角度
Figure BDA0002036851610000045
S813、判断角度
Figure BDA0002036851610000046
是否小于车辆最大转角
Figure BDA0002036851610000047
若是,更新角度
Figure BDA0002036851610000048
对应待选父节点qparent1/qparent2的代价值为∞,并进入步骤S812;否则,进入步骤S814;
S814、计算角度
Figure BDA0002036851610000049
对应的待选父节点qparent1/qparent2的代价值:
Cparent=A*ed1
其中,Cparent为选取待选父节点的代价值;A为常数;d1为树根到qnew1/qnew2的距离,e为自然对数;
S815、判断所有的待选父节点qparent1/qparent2是否存在未被选取的待选父节点qparent1/qparent2;若存在,则返回步骤S812,否则进入步骤S816;
S816、选取代价值最小、且其所在区域内无障碍物的待选父节点qparent1/ qparent2作为节点qnew1/qnew2的父节点。
进一步地,根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/ qnew2的子节点进一步包括:
S821、选取以节点qnew1/qnew2为中心,预设半径R为半径所覆盖区域内的所有节点为待选子节点qchlid1/qchlid2
S822、选取未选取的待选子节点qchlid1/qchlid2,计算节点qnew1/qnew2、待选父节点qnew1/qnew2的父节点和待选子节点qchlid1/qchlid2所形成角度
Figure BDA0002036851610000041
S823、判断角度
Figure BDA0002036851610000042
是否小于车辆最大转角
Figure BDA0002036851610000043
若是,更新角度
Figure BDA0002036851610000044
对应待选子节点qchlid1/qchlid2的代价值为∞,并进入步骤S822;否则,进入步骤S824;
S824、计算角度
Figure BDA0002036851610000051
对应的待选子节点qchlid1/qchlid2的代价值:
Cchild=B*ed2
其中,Cchild为待选子节点的代价值;B为常数;d2为树根到待选子节点qchlid1/qchlid2的距离,e为自然对数;
S825、判断所有的待选子节点qchlid1/qchlid2是否存在未被选取的待选子节点qchlid1/qchlid2;若存在,则返回步骤S822,否则进入步骤S826;
S826、选取代价值最小、且其所在区域内无障碍物的待选子节点qchlid1/qchlid2作为节点qnew1/qnew2的子节点。
进一步地,基于双向扩展随机树的智能车路径规划方法还包括对所述智能车规划路径进行优化处理:
A1、在智能车规划路径上取出两个未同时选取过的路径点,并计算两个路径点间的距离d1
A2、判断离d1是否大于扩展步长,若是,则进入步骤A3;否则,返回步骤 A1;
A3、在两个路径点的连线上,取距离起点以扩展步长为长度的点作为路径点posenew
A4、选取以路径点posenew所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新路径点posenew,之后返回步骤A1。
进一步地,基于双向扩展随机树的智能车路径规划方法还包括采用贝塞尔曲线对智能车路径进行拟合以得到智能车最终的规划路径。
进一步地,障碍物判断时,采用地图上像素点的灰度值进行确定。
本发明的有益效果为:本方案提供的智能车路径规划方法克服了现有技术中基本的RRT、RRT*、B-RRT*算法存在转折点较多、路线不平滑、规划路径与障碍物过近和曲率不连续等问题,使得所规划的路径能更加满足车辆的运动学模型。采用Reeds-Shepp曲线连接智能车的起始点和目标点,能够解决目标点朝向的问题。
附图说明
图1为基于双向扩展随机树的智能车路径规划方法的流程图。
图2为路径规划示意图,其中(a)为起始点朝向与目标点朝向相同且存在转折的情况,(b)为起始点朝向与目标点朝向垂直的情况。
图3为采样概率小于等于目标偏向概率时,随机点的选取示意图。
图4为节点周围/所在区域内节点的选取示意图。
图5为父节点的选择对比示意图,其中(a)为选择父节点初的示意图;(b) 为B-RRT*算法选择父节点后的示意图;(c)为本方案选择父节点后的示意图。
图6为子节点的选择对比示意图,其中(a)为选择子节点初的示意图;(b) 为B-RRT*算法选择子节点后的示意图;(c)为本方案选择子节点后的示意图。
图7为使用贝塞尔处理过规划路径的示意图。
图8为posenew周围节点选取示意图。
图9为在某学校的校园地图。
图10为在图9内的区域一中生成的规划路径示意图,其中(a)为RRT算法在区域一中生成的规划路径示意图;(b)为RRT*算法在区域一中生成的规划路径示意图;(c)为B-RRT*算法在区域一中生成的规划路径示意图;(d)为采用本方案方法在区域一中生成的规划路径示意图。
图11为在图9内的区域二中生成的规划路径示意图,其中(a)为RRT算法在区域二中生成的规划路径示意图;(b)为RRT*算法在区域二中生成的规划路径示意图;(c)为B-RRT*算法在区域二中生成的规划路径示意图;(d)为采用本方案方法在区域二中生成的规划路径示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
参考图1,图1示出了基于双向扩展随机树的智能车路径规划方法的流程图;如图1所示,该方法S包括步骤S1至步骤S10。
在步骤S1中,采用Reeds-Shepp曲线连接智能车的起始点和目标点;具体参见图2,图中颜色较浅且较细的黑色路径为使用Reeds-Shepp曲线预处理生成的路径,颜色较深且较粗的黑色路径为B-RRT*算法规划的路径。从图2中可以看出使用Reeds-Shepp曲线能够有效解决车辆朝向的问题,并且所生成的轨迹也更加平滑。
在步骤S2中,当生成的曲线上存在障碍物时,从起始点向目标点构建扩展树tree1,从目标点向起始点构建扩展树tree2。当生成的曲线上不存在障碍物时,输出该曲线作为智能车的规划路径,并结束规划。
在步骤S3中,在包含智能车起始点和目标点的地图上生成随机点qrand1/ qrand2,并查找扩展树tree1/tree2距离随机点qrand1/qrand2最近的节点qnearrest1/ qnearrest2
在本发明的一个实施例中,在包含智能车起始点和目标点的地图上生成随机点qrand1/qrand2进一步包括步骤S31至步骤S35。
在步骤S31中,采用随机函数生成一个采样概率,并判断采样概率是否大于目标偏向概率,若大于,则进入步骤S32,否则,进入步骤S33;其中,采样概率小于等于目标偏向概率时,随机点的选取可以参考图3。
在步骤S2中,采用目标点/起始点作为扩展树tree1/扩展树tree2的随机点 qrand1/随机点qrand2
在步骤S33中,在扩展树tree1/扩展树tree2中寻找与目标点/起始点最近的节点qnr1/节点qnr2,以节点qnr1/节点qnr2为中心,以扩展步长为边长的方形区域,并判断方形区域内是否有障碍物;
在步骤S34中,若没有障碍物,则采用节点qnr1/节点qnr2建立滑动窗口,并在滑动窗口中生成qrand1/随机点qrand2;其中滑动窗口中的窗口长度优选为智能车的宽度。
在步骤S35中,若有障碍物,则返回步骤S31。
采用本方案的方式是对随机点qrand1/qrand2的选取进行约束,可以避免采样的随机性,从而避免了采样带来的不确定因素。
在步骤S4中,判断随机点qrand1/qrand2与节点qnearrest1/qnearrest2的连线上是否有障碍物;若有,则返回步骤S3,否则进入步骤S5;本方案中在进行障碍物判断时,均是采用地图上像素点的灰度值进行确定。
在步骤S5中,在随机点qrand1/qrand2与节点qnearrest1/qnearrest2的连线上取位于一个扩展步长长度处的节点qnew1/qnew2
在步骤S6中,判断节点qnearrest1/qnearrest2与节点qnew1/qnew2的连线上是否有障碍物;若有,则返回步骤S3,否则,进入步骤S7;
在步骤S7中,选取以节点qnew1/qnew2所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新节点qnew1/qnew2
在本发明的一个实施例中,选取以节点qnew1/qnew2所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新节点qnew1/qnew2进一步包括步骤 S71至步骤S76。
在步骤S71中,形成一块以节点qnew1/qnew2为中心,扩展步长为1/2边长(扩展步长的长度等于方形区域A1/A2的边长的一半)的方形区域A1/A2,并同时初始化规划路径与障碍物之间的最小距离dmin为无穷大;
在步骤S72中,选取方形区域A1/A2中未选取过的节点poses1/poses2,并判断节点poses1/poses2是否为障碍物;若是,执行步骤S72,否则进入步骤S73;
在步骤S73中,以方形区域A1/A2中的每个节点poses1/poses2为中心,以扩展步长为边长,形成区域B1/B2;
其中,节点poses1/poses2周围或所在区域(方形区域B1/B2)内节点的选取参考图4,poses1={A1,B1,C1,D1},poses11={A11,A12,A13,A14,B11,B12,B13,B14,C11,C12,C13,C14, D11,D12,D13,D14}),扩展步长为车辆宽度的一半,图4中的D为扩展步长。
在步骤S74中,选取区域B1/B2中未选取过的节点poses11/poses21,并判断节点poses11/poses21是否为障碍物;若是,返回步骤S72,否则进入步骤S75;
在步骤S75中,计算节点poses11/poses21到目标点的距离d,并判断距离d 是否小于最小距离dmin;若是,进入步骤S76,否则返回步骤S72;
在步骤S76中,采用距离d更新最小距离dmin,采用距离d对应的节点poses1/poses2更新节点qnew1/qnew2
在步骤S77中,判断方形区域A1/A2,是否存在未选取过得节点。若有,则执行S72。否则,输出节点qnew1/qnew2
在步骤S8中,根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/ qnew2的父节点和子节点。
在本发明的一个实施例中,根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/qnew2的父节点进一步包括:
S811、选取以节点qnew1/qnew2为中心,预设半径R为半径所覆盖区域内的所有节点为待选父节点qparent1/qparent2
S812、选取未选取的待选父节点qparent1/qparent2,计算待选父节点qparent1/ qparent2、待选父节点qparent1/qparent2的父节点和节点qnew1/qnew2所形成角度
Figure BDA0002036851610000103
S813、判断角度
Figure BDA0002036851610000104
是否小于车辆最大转角
Figure BDA0002036851610000105
若是,更新角度
Figure BDA0002036851610000106
对应待选父节点qparent1/qparent2的代价值为∞,并进入步骤S812;否则,进入步骤S814;
S814、计算角度
Figure BDA0002036851610000107
对应的待选父节点qparent1/qparent2的代价值:
Cparent=A*ed1
其中,Cparent为选取待选父节点的代价值;A为常数;d1为树根到qnew1/qnew2的距离,e为自然对数;
S815、判断所有的待选父节点qparent1/qparent2是否存在未被选取的待选父节点qparent1/qparent2;若存在,则返回步骤S812,否则进入步骤S816;
S816、选取代价值最小、且其所在区域内无障碍物的待选父节点qparent1/ qparent2作为节点qnew1/qnew2的父节点。
本方案的父节点的选取可以参考图5中的(a)和(c),图5中的(a)是以R 为半径⑨为圆心的重选选择父节点的过程,图5的(b)是使用基本B-RRT*算法执行重选父节点的结果(
Figure BDA0002036851610000108
→⑤→⑧→⑨其路径的长度为11),图5的(c)是采用本方案选择的父节点的结果(
Figure BDA0002036851610000109
→④→⑨其路径的长度为14)。
通过图5的(b)和(c)中可知,虽然(b)的路径比(c)短,但是
Figure BDA0002036851610000101
(
Figure BDA0002036851610000102
为路径⑤→⑧→⑨的角度)不满足路径转角小于智能车最大转角的要求,所以放弃节点⑧作为节点⑨的父点。
实施时,本方案优选根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/qnew2的子节点进一步包括:
S821、选取以节点qnew1/qnew2为中心,预设半径R为半径所覆盖区域内的所有节点为待选子节点qchlid1/qchlid2
S822、选取未选取的待选子节点qchlid1/qchlid2,计算节点qnew1/qnew2、待选父节点qnew1/qnew2的父节点和待选子节点qchlid1/qchlid2所形成角度
Figure BDA0002036851610000111
S823、判断角度
Figure BDA0002036851610000112
是否小于车辆最大转角
Figure BDA0002036851610000113
若是,更新角度
Figure BDA0002036851610000114
对应待选子节点qchlid1/qchlid2的代价值为∞,并进入步骤S822;否则,进入步骤S824;
S824、计算角度
Figure BDA0002036851610000115
对应的待选子节点qchlid1/qchlid2的代价值:
Cchild=B*ed2
其中,Cchild为待选子节点的代价值;B为常数;d2为树根到待选子节点qchlid1/qchlid2的距离,e为自然对数;
S825、判断所有的待选子节点qchlid1/qchlid2是否存在未被选取的待选子节点qchlid1/qchlid2;若存在,则返回步骤S822,否则进入步骤S826;
S826、选取代价值最小、且其所在区域内无障碍物的待选子节点qchlid1/qchlid2作为节点qnew1/qnew2的子节点。
本方案的子节点选取过程(重布随机树过程)可以参考图6,其中(a)为新生成的节点qnew⑨为圆心,R为半径并以④为父节点,圆内其他的节点为子节点(qnew的父节点除外)。图6的(a)表示基本B-RRT*算法重新选择子节点的结果, (b)表为基本B-RRT*算法重新布线的结果(
Figure BDA0002036851610000118
→④→⑨→⑥其路径为15);(c)采用本方案的子节点选取方法重新布线的结果(
Figure BDA0002036851610000119
④→⑨→⑩其路径为16)。
通过图6的(b)和(c)中可知虽然(b)的路径比(c)短,但是
Figure BDA0002036851610000116
(
Figure BDA0002036851610000117
为路径④→⑨→⑥的角度)不满足路径转角小于智能车最大转角的要求,所以放弃节点⑥作为节点⑨的子节点。
在步骤S9中,判断扩展树tree1和扩展树tree2之间的距离是否小于设定阈值,若大于,则返回步骤S3,否则,进入步骤S10;
S10、根据扩展树tree1和扩展树tree2相遇的节点qa和qb,反向选取节点 qa和qb分别到起始点和目标点的路径构成智能车规划路径。
在本发明的一个实施例中,基于双向扩展随机树的智能车路径规划方法还包括对所述智能车规划路径进行优化处理:
A1、在智能车规划路径上取出两个未同时选取过的路径点,并计算两个路径点间的距离d1
A2、判断离d1是否大于扩展步长,若是,则进入步骤A3;否则,返回步骤 A1;
A3、在两个路径点的连线上,取距离起点以扩展步长为长度的点作为路径点posenew
A4、选取以路径点posenew所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新路径点posenew,之后返回步骤A1。
实施时,本方案优选选取以路径点posenew所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新路径点posenew进一步包括:
1)首先获得节点posenew,初始化dmin(初始化定义dmin等于无穷大)。
2)获得posenew四周的节点poses111,依次从中选取pose并判断是否为障碍物。如果是,则从poses111中从新选取pose。否则在pose节点四周添加poses22,如图8所示,poses111={A11,B11,C11,D11},poses22={A11 1,A11 2,A11 3,A11 4,B11 1,B11 2,B11 3,B11 4, C11 1,C11 2,C11 3,C11 4,D11 1,D11 2,D11 3,D11 4};
3)在poses22中依次选择pose1,并判断poses1是否为障碍物。如果是,则执行b)。否则,计算pose到目标点距离d。
4)判断d和dmin大小,如果dmin>d,更新posenew和dmin。否则执行b)。
5)输出posenew
实施时,本方案优选基于双向扩展随机树的智能车路径规划方法还包括采用贝塞尔曲线对智能车路径进行拟合以得到智能车最终的规划路径。采用贝塞尔曲线进行优化可以参考图7,优化过程的公式为:
B(t)=(1-t)2*P0+2*t(1-t)*P1+t2*P2t∈[0,1]
其中,P0,P1,P2为二阶贝塞尔曲线的控制点。
本方案在对规划过程中的步骤S3至步骤S8中采用的符号“/”表示或的意思,且在对应过程中“/”前面部分对应所有“/”前面部分内容,“/”后面部分对应所有“/”后面部分内容。
下面以步骤S3至步骤S8为例对“/”前面和后面的内容进行分开说明:
对于扩展树tree1上节点实现步骤S3至S8的处理过程如下:
S3、在包含智能车起始点和目标点的地图上生成随机点qrand1,并查找扩展树tree1距离随机点qrand1最近的节点qnearrest1
S4、判断随机点qrand1与节点qnearrest1的连线上是否有障碍物;若有,则返回步骤S3,否则进入步骤S5;
S5、在随机点qrand1与节点qnearrest1的连线上取位于一个扩展步长长度处的节点qnew1
S6、判断节点qnearrest1与节点qnew1的连线上是否有障碍物;若有,则返回步骤S3,否则,进入步骤S7;
S7、选取以节点qnew1所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新节点qnew1
S8、根据智能车辆最大转向约束及节点qnew1,选取节点qnew1的父节点和子节点。
对于扩展树tree2上节点实现步骤S3至S8的处理过程如下:
S3、在包含智能车起始点和目标点的地图上生成随机点qrand2,并查找扩展树tree2距离随机点qrand2最近的节点qnearrest2
S4、判断随机点qrand2与节点qnearrest2的连线上是否有障碍物;若有,则返回步骤S3,否则进入步骤S5;
S5、在随机点qrand2与节点qnearrest2的连线上取位于一个扩展步长长度处的节点qnew2
S6、判断节点qnearrest2与节点qnew2的连线上是否有障碍物;若有,则返回步骤S3,否则,进入步骤S7;
S7、选取以节点qnew2所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新节点qnew2
S8、根据智能车辆最大转向约束及节点qnew2,选取节点qnew2的父节点和子节点。
本方案的扩展树tree1和扩展树tree2在节点是同步实现步骤S3至步骤S8 的,本方案每个步骤中“/”前面后面内容跟上面步骤S3至步骤S8是类似的,本方案就不在对每个步骤进行单独说明。
下面结合具体的实例,对本方案提供的智能车路径规划方法的效果进行说明:
分别将RRT、RRT*、B-RRT*和本方案算法应用于某大学学校园的栅格地图 (4846×2816)上(如图9所示),栅格大小为0.298米/pixel,选择不同的地点(此处的不同地点指代图9中的区域一和区域二)进行试验。实验基于ROS(开源机器人操作***)平台,并使用Rviz对规划结果进行可视化显示,计算机配置为:ubuntu14.04LTS,处理器intel i5-6500,主频为3.2Ghz,运行内存为8G。
本次实验通过在栅格地图上选择起始点和终点,使用RRT、RRT*、B-RRT* 和本方案算法规划。通过路径的长度、到障碍物的距离和转折数量(
Figure BDA0002036851610000153
Figure BDA0002036851610000154
则认为是转折)等方式来评价路径的好坏。
对区域一和区域二进行路径规划,规划结果见图10和图11。其中,表1和表2分别为连续转折区域、直角转折曲线区域的指标对比。
图10和图11中的方形节点为起始点向目标点的采样节点,圆形节点为目标点向起始点的采样节点。
从图10和图11可以看出,本方案算法所规划出来的路径较为平滑,且具有连续的曲率,能够解决基本RRT、RRT*、B-RRT*算法在初始时刻路径不合理的问题,同时与障碍物形成一定的安全距离并且能够较好的满足车辆运动学特性。
表1区域一测试对比
Figure BDA0002036851610000151
表2区域二测试对比
Figure BDA0002036851610000152
通过表1和表2能够看到虽然基本RRT、RRT*、B-RRT*算法能够在较短的时间内规划出路径,但是所规划的路径转折角度过大、路径点贴合障碍物,导致路径不能被智能车使用。而本方案算法在牺牲较短时间的基础上获得了更加平滑、安全且更符合车辆运动学的路径。

Claims (7)

1.基于双向扩展随机树的智能车路径规划方法,其特征在于,包括:
S1、采用Reeds-Shepp曲线连接智能车的起始点和目标点;
S2、当生成的曲线上存在障碍物时,从起始点向目标点构建扩展树tree1,从目标点向起始点构建扩展树tree2;
S3、在包含智能车起始点和目标点的地图上生成随机点qrand1/qrand2,并查找扩展树tree1/tree2距离随机点qrand1/qrand2最近的节点qnearrest1/qnearrest2
S4、判断随机点qrand1/qrand2与节点qnearrest1/qnearrest2的连线上是否有障碍物;若有,则返回步骤S3,否则进入步骤S5;
S5、在随机点qrand1/qrand2与节点qnearrest1/qnearrest2的连线上取位于一个扩展步长长度处的节点qnew1/qnew2
S6、判断节点qnearrest1/qnearrest2与节点qnew1/qnew2的连线上是否有障碍物;若有,则返回步骤S3,否则,进入步骤S7;
S7、选取以节点qnew1/qnew2所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新节点qnew1/qnew2
S8、根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/qnew2的父节点和子节点;
S9、判断扩展树tree1和扩展树tree2之间的距离是否小于设定阈值,若大于,则返回步骤S3,否则,进入步骤S10;
S10、根据扩展树tree1和扩展树tree2相遇的节点qa和qb,反向选取节点qa和qb分别到起始点和目标点的路径构成智能车规划路径;
选取以节点qnew1/qnew2所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新节点qnew1/qnew2进一步包括:
S71、形成一块以节点qnew1/qnew2为中心,扩展步长为1/2边长的方形区域A1/A2,并同时初始化规划路径与障碍物之间的最小距离dmin为无穷大;
S72、选取方形区域A1/A2中未选取过的节点poses1/poses2,并判断节点poses1/poses2是否为障碍物;若是,执行步骤S72,否则进入步骤S73;
S73、以方形区域A1/A2中的每个节点poses1/poses2为中心,以扩展步长为边长,形成区域B1/B2;
S74、选取区域B1/B2中未选取过的节点poses11/poses21,并判断节点poses11/poses21是否为障碍物;若是,返回步骤S72,否则进入步骤S75;
S75、计算节点poses11/poses21到目标点的距离d,并判断距离d是否小于最小距离dmin;若是,进入步骤S76,否则返回步骤S72;
S76、采用距离d更新最小距离dmin,采用距离d对应的节点poses1/poses2更新节点qnew1/qnew2
S77、判断方形区域A1/A2,是否存在未选取过得节点;若有,则执行S72;否则,输出节点qnew1/qnew2
2.根据权利要求1所述的基于双向扩展随机树的智能车路径规划方法,其特征在于,在包含智能车起始点和目标点的地图上生成随机点qrand1/qrand2进一步包括:
S31、采用随机函数生成一个采样概率,并判断采样概率是否大于目标偏向概率,若大于,则进入步骤S32,否则,进入步骤S33;
S32、采用目标点/起始点作为扩展树tree1/扩展树tree2的随机点qrand1/随机点qrand2
S33、在扩展树tree1/扩展树tree2中寻找与目标点/起始点最近的节点qnr1/节点qnr2,以节点qnr1/节点qnr2为中心,以扩展步长为边长的方形区域,并判断方形区域内是否有障碍物;
S34、若没有障碍物,则采用节点qnr1/节点qnr2建立滑动窗口,并在滑动窗口中生成qrand1/随机点qrand2
S35、若有障碍物,则返回步骤S31。
3.根据权利要求1所述的基于双向扩展随机树的智能车路径规划方法,其特征在于,根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/qnew2的父节点进一步包括:
S811、选取以节点qnew1/qnew2为中心,预设半径R为半径所覆盖区域内的所有节点为待选父节点qparent1/qparent2
S812、选取未选取的待选父节点qparent1/qparent2,计算待选父节点qparent1/qparent2、待选父节点qparent1/qparent2的父节点和节点qnew1/qnew2所形成角度
Figure FDA0003649181020000021
S813、判断角度
Figure FDA0003649181020000022
是否小于车辆最大转角
Figure FDA0003649181020000023
若是,更新角度
Figure FDA0003649181020000024
对应待选父节点qparent1/qparent2的代价值为∞,并进入步骤S812;否则,进入步骤S814;
S814、计算角度
Figure FDA0003649181020000025
对应的待选父节点qparent1/qparent2的代价值:
Cparent=A*ed1
其中,Cparent为选取待选父节点的代价值;A为常数;d1为树根到qnew1/qnew2的距离,e为自然对数;
S815、判断所有的待选父节点qparent1/qparent2是否存在未被选取的待选父节点qparent1/qparent2;若存在,则返回步骤S812,否则进入步骤S816;
S816、选取代价值最小、且其所在区域内无障碍物的待选父节点qparent1/qparent2作为节点qnew1/qnew2的父节点。
4.根据权利要求3所述的基于双向扩展随机树的智能车路径规划方法,其特征在于,根据智能车辆最大转向约束及节点qnew1/qnew2,选取节点qnew1/qnew2的子节点进一步包括:
S821、选取以节点qnew1/qnew2为中心,预设半径R为半径所覆盖区域内的所有节点为待选子节点qchlid1/qchlid2
S822、选取未选取的待选子节点qchlid1/qchlid2,计算节点qnew1/qnew2、待选父节点qnew1/qnew2的父节点和待选子节点qchlid1/qchlid2所形成角度
Figure FDA0003649181020000031
S823、判断角度
Figure FDA0003649181020000032
是否小于车辆最大转角
Figure FDA0003649181020000033
若是,更新角度
Figure FDA0003649181020000034
对应待选子节点qchlid1/qchlid2的代价值为∞,并进入步骤S822;否则,进入步骤S824;
S824、计算角度
Figure FDA0003649181020000035
对应的待选子节点qchlid1/qchlid2的代价值:
Cchild=B*ed2
其中,Cchild为待选子节点的代价值;B为常数;d2为树根到待选子节点qchlid1/qchlid2的距离,e为自然对数;
S825、判断所有的待选子节点qchlid1/qchlid2是否存在未被选取的待选子节点qchlid1/qchlid2;若存在,则返回步骤S822,否则进入步骤S826;
S826、选取代价值最小、且其所在区域内无障碍物的待选子节点qchlid1/qchlid2作为节点qnew1/qnew2的子节点。
5.根据权利要求1所述的基于双向扩展随机树的智能车路径规划方法,其特征在于,还包括对所述智能车规划路径进行优化处理:
A1、在智能车规划路径上取出两个未同时选取过的路径点,并计算两个路径点间的距离d1
A2、判断离d1是否大于扩展步长,若是,则进入步骤A3;否则,返回步骤A1;
A3、在两个路径点的连线上,取距离起点以扩展步长为长度的点作为路径点posenew
A4、选取以路径点posenew所在区域内使得规划路径与障碍物之间的距离大于安全距离的节点更新路径点posenew,之后返回步骤A1。
6.根据权利要求1-5任一所述的基于双向扩展随机树的智能车路径规划方法,其特征在于,还包括采用贝塞尔曲线对智能车路径进行拟合以得到智能车最终的规划路径。
7.根据权利要求1-5任一所述的基于双向扩展随机树的智能车路径规划方法,其特征在于,障碍物判断时,采用地图上像素点的灰度值进行确定。
CN201910328124.5A 2019-04-23 2019-04-23 基于双向扩展随机树的智能车路径规划方法 Active CN109990796B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910328124.5A CN109990796B (zh) 2019-04-23 2019-04-23 基于双向扩展随机树的智能车路径规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910328124.5A CN109990796B (zh) 2019-04-23 2019-04-23 基于双向扩展随机树的智能车路径规划方法

Publications (2)

Publication Number Publication Date
CN109990796A CN109990796A (zh) 2019-07-09
CN109990796B true CN109990796B (zh) 2022-07-19

Family

ID=67135122

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910328124.5A Active CN109990796B (zh) 2019-04-23 2019-04-23 基于双向扩展随机树的智能车路径规划方法

Country Status (1)

Country Link
CN (1) CN109990796B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110986953B (zh) * 2019-12-13 2022-12-06 达闼机器人股份有限公司 路径规划方法、机器人及计算机可读存储介质
CN111060125B (zh) * 2019-12-30 2021-09-17 深圳一清创新科技有限公司 碰撞检测方法、装置、计算机设备和存储介质
CN111176288A (zh) * 2020-01-07 2020-05-19 深圳南方德尔汽车电子有限公司 基于Reedsshepp全局路径规划方法、装置、计算机设备及存储介质
CN111369066B (zh) * 2020-03-09 2022-06-24 广东南方数码科技股份有限公司 路径规划方法、装置、电子设备及可读存储介质
CN111397598B (zh) * 2020-04-16 2022-02-01 苏州大学 人机共融环境中移动型机器人路径规划采样方法及***
CN111678523B (zh) * 2020-06-30 2022-05-17 中南大学 一种基于star算法优化的快速bi_rrt避障轨迹规划方法
CN111982125A (zh) * 2020-08-31 2020-11-24 长春工业大学 一种基于改进蚁群算法的路径规划方法
CN113064426B (zh) * 2021-03-17 2022-03-15 安徽工程大学 一种改进双向快速搜索随机树算法的智能车路径规划方法
CN113119115A (zh) * 2021-03-22 2021-07-16 深圳市优必选科技股份有限公司 一种机械臂运动规划方法、装置、可读存储介质及机械臂
CN112884256B (zh) * 2021-04-28 2021-07-27 深圳大学 路径规划方法、装置、计算机设备和存储介质
CN113188562B (zh) * 2021-07-01 2022-03-01 新石器慧通(北京)科技有限公司 可行驶区域的路径规划方法、装置、电子设备及存储介质
CN113934206B (zh) * 2021-07-22 2024-01-16 浙江科技学院 移动机器人路径规划方法、装置、计算机设备和存储介质
CN114137991B (zh) * 2022-02-08 2022-04-26 河北工业大学 一种基于二阶贝塞尔曲线的机器人连续路径优化方法
CN114500359B (zh) * 2022-04-15 2022-07-12 深圳市永达电子信息股份有限公司 集群动态组网方法和集群动态***
CN116070807B (zh) * 2022-11-08 2024-03-26 国电湖北电力有限公司鄂坪水电厂 一种基于空间关系的电站巡检路径规划方法及装置
CN115950439B (zh) * 2023-03-15 2023-06-02 季华实验室 双向rrt路径规划方法、装置、电子设备及存储介质

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107883961A (zh) * 2017-11-06 2018-04-06 哈尔滨工程大学 一种基于Smooth‑RRT算法的水下机器人路径优化方法
CN108073176B (zh) * 2018-02-10 2020-08-18 西安交通大学 一种改进型D*Lite车辆动态路径规划方法
CN108444489A (zh) * 2018-03-07 2018-08-24 北京工业大学 一种改进rrt算法的路径规划方法
CN108896052A (zh) * 2018-09-20 2018-11-27 鲁东大学 一种基于动态复杂环境下的移动机器人平滑路径规划方法
CN109582024A (zh) * 2018-12-27 2019-04-05 济南大学 一种智能铲运机的路径规划方法

Also Published As

Publication number Publication date
CN109990796A (zh) 2019-07-09

Similar Documents

Publication Publication Date Title
CN109990796B (zh) 基于双向扩展随机树的智能车路径规划方法
CN108762270B (zh) 变概率双向快速搜索随机树改进路径规划算法
CN106125764B (zh) 基于a*搜索的无人机路径动态规划方法
Zips et al. Optimisation based path planning for car parking in narrow environments
Yang et al. Spline-based RRT path planner for non-holonomic robots
Li et al. Multivehicle cooperative local mapping: A methodology based on occupancy grid map merging
US20160375901A1 (en) System and Method for Controlling Semi-Autonomous Vehicles
CN109764886A (zh) 一种路径规划方法
CN112590775B (zh) 一种自动泊车方法、装置、车辆及存储介质
CN103529843A (zh) Lambda*路径规划算法
Yang Anytime synchronized-biased-greedy rapidly-exploring random tree path planning in two dimensional complex environments
Manav et al. A novel cascade path planning algorithm for autonomous truck-trailer parking
Kim et al. Cloud RRT*: Sampling cloud based RRT
Du et al. An improved RRT-based motion planner for autonomous vehicle in cluttered environments
Wu et al. Mobile robot path planning based on a generalized wavefront algorithm
CN109556623B (zh) 融合触须算法与切线图法的迭代式路径规划算法
CN112197783B (zh) 一种考虑车头指向的两阶段多抽样的rrt路径规划方法
Sudhakara et al. Probabilistic roadmaps-spline based trajectory planning for wheeled mobile robot
Jafarzadeh et al. An exact geometry–based algorithm for path planning
CN113358129B (zh) 基于Voronoi图的避障最短路径规划方法
CN115328111A (zh) 一种用于自动驾驶车辆的避障路径规划方法
CN115454106B (zh) 一种基于双向搜索rrt*的auv回坞路径规划方法
Primatesta et al. Dynamic trajectory planning for mobile robot navigation in crowded environments
Wang et al. Reference path correction for autonomous ground vehicles driving over rough terrain
Yu et al. Robot path planning based on improved A* algorithm

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