CN109443363A - 分治确定性路径寻优算法 - Google Patents

分治确定性路径寻优算法 Download PDF

Info

Publication number
CN109443363A
CN109443363A CN201811332691.XA CN201811332691A CN109443363A CN 109443363 A CN109443363 A CN 109443363A CN 201811332691 A CN201811332691 A CN 201811332691A CN 109443363 A CN109443363 A CN 109443363A
Authority
CN
China
Prior art keywords
grid
obstacle
path
angle
rectangle
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.)
Granted
Application number
CN201811332691.XA
Other languages
English (en)
Other versions
CN109443363B (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.)
Xiamen University
Original Assignee
Xiamen 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 Xiamen University filed Critical Xiamen University
Priority to CN201811332691.XA priority Critical patent/CN109443363B/zh
Publication of CN109443363A publication Critical patent/CN109443363A/zh
Application granted granted Critical
Publication of CN109443363B publication Critical patent/CN109443363B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/20Instruments for performing navigational calculations

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Processing Or Creating Images (AREA)
  • Design And Manufacture Of Integrated Circuits (AREA)

Abstract

分治确定性路径寻优算法,涉及路径规划领域。将环境地图分割成栅格,每个栅格只有空白状态或赋值状态,空白状态表示该栅格可以通行,赋值状态表示该栅格为障碍区域,无法通行;在相隔并且中间存在赋值栅格的空白栅格中确定起点和终点位置;根据障碍区域的位置及数量,提出一种角点合力法求取栅格障碍的角点位置并标记;将所有的标记角点两两相连;将每段可行路径分割为M个线段,每一线段的线长范围用一个矩形包围,求出该矩形四个角的坐标值并取整处理,四个角的坐标值与栅格化处理生成的矩阵坐标比较,判断每个坐标点是否位于空白栅格中;以起点作为出发点,终点作为结束点,比较各搜索路径的距离大小,从中搜索出一条从起点到终点的最优路径。

Description

分治确定性路径寻优算法
技术领域
本发明涉及路径规划领域,具体应用于精确求取复杂障碍下的机器人的最优路径。
背景技术
栅格路径规划的优点在于地图经栅格化处理后,可以将障碍分割为多个固定尺度的小栅格,这样可以精确识别出障碍的位置、大小及形状,使机器人更易于实现避障。但是地图环境经栅格化后,可选取的路线只能为栅格的四个角或中心点,不能直接选择两个目标的矢量连线路径,因此最终得到的路径并非最短路径。与栅格路径规划不同,在矢量路径规划中,选取的路线则可以为两个目标点的矢量连线路径,并没有栅格路径中的必须经过栅格中心点的约束。
传统的路径规划算法有A*算法,Dijkstra算法,PRM路径规划算法等。当A*算法和Dijkstra算法应用于栅格路径规划中时,这两种算法在某种程度上属于一种贪婪算法,通过不停的比较周围可行路径的权值来选取最优路径。然而随着地图的面积的扩大,这两种算法的搜索范围也将增加,算法的时间复杂度和搜索复杂度随之增大。即使是在无障碍的空白区域中,从起点到终点,这两种算法都无法直接搜索出起点到终点的直线路径,仍需要对周围的路径进行搜索比较。PRM路径规划算法一般应用于矢量路径规划中。该算法通常在地图中的可行范围内随机生成一定数量的可行点,然后将这些可行点连接,在未穿过障碍区域的连线中通过A*算法求取出最优路径。该算法的缺点在于,只有生成大量的随机可行点才有可能得到有一条从起点到终点的路径,且算法的复杂度也会随着可行点数量的增加而增大,尤其在存在狭缝通道的障碍地图中,一旦没有随机点落在该狭缝中,PRM算法将无法搜索出起点到终点的可行路径。同时,由于这些可行点是随机生成的,因此得到的最优路径并不是最短路径。
发明内容
本发明的目的在于提供可通过从经障碍区域角点部分的可行连线路径集合中搜索出一条最优路径的分治确定性路径寻优算法。
本发明包括以下步骤:
1)将需要搜索的环境地图分割成栅格,每个栅格只有空白状态或赋值状态,空白状态表示该栅格可以通行,赋值状态表示该栅格为障碍区域,无法通行;根据实际已知信息,在相隔并且中间存在赋值栅格的空白栅格中确定起点和终点位置;
在步骤1)中,所述需要搜索的环境地图可进行栅格化处理,根据障碍复杂程度及要求精度,设置所要产生地图矩阵的大小,生成一个N×N的方阵,将环境地图分割成多个规则的正方形栅格,其中空白栅格赋值为“0”,表示该区域可以通行,而障碍栅格赋值为“1”,表示该区域禁止通行;所述空白状态或赋值状态可对应于计算机中的“0”和“1”两种状态。
2)根据障碍区域的位置及数量,提出一种角点合力法求取栅格障碍的角点位置,并进行标记;
在步骤2)中,所述栅格障碍的角点位置可使用角点合力法计算并标记,所述角点合力法的计算步骤可为:
(1)对单位栅格障碍周围的上、下、左、右四个方向的栅格可行区域加权赋值1;
(2)当单位障碍的左或右和上或下的可行区域块都具有加权值1时,单位障碍被定义为障碍连通区域的角点位置;
(3)当单位障碍的左方和上方空白栅格的权值都为1时,标记该单位障碍的左上角的空白栅格;同理,当单位障碍的右方和上方空白栅格的权值都为1时,标记该单位障碍的右上角的空白栅格;当单位障碍的右方和下方空白栅格的权值都为1时,标记该单位障碍的右下角的空白栅格;当单位障碍的左方和下方空白栅格的权值都为1时,标记该单位障碍的左下角的空白栅格,通过力的合成公式:
规定力的合成方向指向障碍的角点处,其中,F为分力F1和分力F2的合力,α为两个分力的夹角。
3)将所有的标记角点两两相连,形成多条可行路径;然后将每段可行路径分割为M个线段,每一线段的线长范围用一个矩形包围,接着求出该矩形四个角的坐标值,并对坐标值进行取整处理,把这个矩形的四个角的坐标值与步骤1)中经栅格化处理生成的矩阵坐标进行比较,判断每个坐标点是否位于空白栅格中;若是,则为可行路径,否则为穿过障碍的路径,将其剔除可行路径集合外;
在步骤3)中,所述将所有的标记角点两两相连,形成多条可行路径,可连接各标记角点、起点和终点,产生多条路径;接着剔除经过障碍的连线路径,生成可行路径集合;引入微元的思想,将每一段连线微元化成M个线段,每一线段用一个矩形来等同其线长的范围;然后求出该矩形四个角的坐标值,并对坐标值进行取整处理,把矩形的四个角的坐标值与步骤1)中经栅格化处理生成的矩阵坐标进行比较,若这四个对角坐标中有一个坐标点位于障碍区域中,则该连线被视为穿过了障碍区域,由此被剔除出可行连线路径的集合外;剔除经过障碍连线的具体方法可为:引入微元法的思想,将所有的线段路径分割成多个线段,则每个线段所经过的地图范围用一个矩形包围,该矩形的四个角点的坐标表示如下:
式中,(x3,1,y3,1),(x3,2,y3,2),(x3,3,y3,3),(x3,4,y3,4)分别为矩形的四个定角坐标;为向下取整函数,为向上取整函数;dir为连线与地图水平轴的夹角,且dir=arctan((y2-y1)/(x2-x1));d表示每段小线段的长度,且存在关系d=D/N,N=D/0.5;因此只需通过判断矩形的四个顶角坐标点是否存在位于障碍区域的坐标,便可进一步判断出该连线路径是否穿过障碍区域,即是否属于可行路径集合;若穿过,则应删除。
4)以起点作为搜索出发点,终点作为搜索结束点,比较各搜索路径的距离大小,从中搜索出一条从起点到终点的最优路径。
在步骤4)中,所述比较各搜索路径的距离大小可采用动态对比的方法。
本发明通过从经过障碍角点区域的可行连线路径中搜索出一条最优路径,搜索最优路径所消耗的时间比Dijkstra算法更少,搜索范围比A*算法和Dijkstra算法更小,搜索得到的最优路径距离为最短,本发明搜索最优路径所消耗的时间比Dijkstra算法更少,搜索范围比Dijkstra算法更小,搜索得到的最优路径距离为最短。与PRM算法相比,本发明提出的算法的标记角点是随独立的障碍连通区域的角点的数量而生成的,并不是随机生成的,对于复杂障碍尤其是狭缝障碍等均可搜索出一条最优路径,并且可以保证该路径为最短路程。
附图说明
图1为本发明的流程图。
图2为两种不同的栅格类型地图算法测试。在图2中,(a)为“迷宫型”栅格类型地图,(b)为“陷阱型”栅格类型地图。
图3为测试本发明所提角点合力法的两幅障碍地图的赋值结果。在图3中,(a)为“迷宫型”障碍地图,(b)为“陷阱型”障碍地图。
图4为四种障碍栅格区域周围的空白栅格区域的搜索及加权赋值示意图。在图4中,(a)为左上角障碍栅格区域,(b)为右上角障碍栅格区域,(c)为右下角障碍栅格区域,(d)为左下角障碍栅格区域。
图5为使用对角合力法将图4所示障碍类型的赋值和标记结果。在图5中,(a)为左上角障碍类型,(b)为右上角障碍类型,(c)为右下角障碍类型,(d)为左下角障碍类型。
图6为测试本发明的两幅障碍地图的角点标记结果。在图6中,(a)为“迷宫型”障碍地图,(b)为“陷阱型”障碍地图。
图7为经过栅格的最短路径的几何示意图。
图8为采用微元法判断连线路径是否穿过障碍的示意图。
图9为剔除图6穿过障碍区域的标记点间的连线路径的仿真结果。在图9中,(a)为“迷宫型”障碍地图,(b)为“陷阱型”障碍地图。
图10为本发明中举例说明的一个具有障碍的环境地图。
图11为可行路径记录表Openlist从搜索开始到结束的部分更新记录。
图12为图10中起点、终点和各角点之间可行路径的距离记录。在图12中,(a)为所有角点的可行路径距离的更新记录,(b)为起点的可行路径距离的更新记录。
图13为图10中起点到各角点及终点的可行路径距离的部分更新记录。
图14为测试本发明的两幅障碍地图搜索的全局最优路径与A*算法和Dijkstra算法搜索的全局最优路径的仿真结果。在图14中,实线“D&C”为本发明的仿真结果,虚线“A*”为A*算法的仿真结果,点线“Dij”为Dijkstra算法的仿真结果;(a)为测试“迷宫型”障碍地图的仿真结果,(b)为测试“陷阱型”障碍地图的仿真结果。
具体实施方式
以下实施例将结合附图对本发明作进一步的说明。
参见图1,本发明实施例包括以下步骤:
1)将需要搜索的地图分割成规则且均匀的小栅格,每个栅格只有空白或赋值两种状态,这对应于计算机中的“0”和“1”两种状态。空白状态表示该栅格可以通行,赋值状态表示该栅格为障碍区域,无法通行;根据实际已知信息,在相隔并且中间存在赋值栅格的空白栅格中确定起点和终点位置;
在步骤1)中,对环境地图进行栅格化处理,根据障碍复杂程度及要求精度,设置所要产生地图矩阵的大小,生成一个N×N的方阵,将环境地图分割成多个规则的正方形小栅格。其中空白栅格赋值为“0”,表示该区域可以通行,而障碍栅格赋值为“1”,表示该区域禁止通行。本发明以2个50×50方阵为仿真对象,赋值为0的栅格以空白的形式表示,赋值为1(即障碍区域)的栅格以涂黑表示,如图2(a)和图2(b)所示。
2)根据障碍区域的位置及数量,提出一种角点合力法求取栅格障碍的角点位置,并进行标记;
在步骤2)中,使用角点合力法计算并标记出障碍角点,其效果如图6(a)和图6(b)所示。角点合力法的详细计算步骤为:
(1)对单位栅格障碍周围的上、下、左、右四个方向的栅格可行区域加权赋值1;
(2)当单位障碍的左或右和上或下的可行区域块都具有加权值1时,该单位障碍被定义为障碍连通区域的角点位置,如图3(a)和图3(b)所示;
(3)当单位障碍的左方和上方空白栅格的权值都为1时,如图4(a)所示,标记该单位障碍的左上角的空白栅格;同理,当单位障碍的右方和上方空白栅格的权值都为1时,如图4(b)所示,标记该单位障碍的右上角的空白栅格;当单位障碍的右方和下方空白栅格的权值都为1时,如图4(c)所示,标记该单位障碍的右下角的空白栅格;当单位障碍的左方和下方空白栅格的权值都为1时,如图4(d)所示,标记该单位障碍的左下角的空白栅格。通过力的合成公式:
得到上述4种情况的力的合成方向如图5(a)~(d)所示,规定力的合成方向指向障碍的角点处。其中,F为分力F1和分力F2的合力,α为两个分力的夹角。以图5(a)为例,障碍栅格(黑色区域)对附近的具有权值为1的空白栅格(白色区域)根据其在单位障碍的方位产生向上或向下或向左或向右方向力的作用。因此,该例产生指向上方的力F1和指向左方的力F2,由上述力的合成原理公式,得到合力F,合力F指向的空白栅格位置,并标记为角点。
3)将所有的标记角点两两相连,形成的多条可行路径;然后将每段可行路径分割为M个小线段,每一小线段的线长范围用一个小矩形包围,如图8所示,接着求出该矩形四个角的坐标值,并对坐标值进行取整处理,把这个小矩形的四个角的坐标值与步骤1)中经栅格化处理生成的矩阵坐标进行比较,判断每个坐标点是否位于空白栅格中。若是,则该这条路径为可行路径,否则为穿过障碍的路径,应将其剔除可行路径集合外;
在步骤3)中,连接各标记角点,起点和终点,产生多条路径。接着剔除经过障碍的连线路径,生成可行路径集合。引入微元的思想,将每一段连线微元化成M个小线段,每一小线段用一个小矩形来等同其线长的范围。然后求出该矩形四个角的坐标值,并对坐标值进行取整处理,把这个小矩形的四个角的坐标值与步骤1)中经栅格化处理生成的矩阵坐标进行比较,若这四个对角坐标中有一个坐标点位于障碍区域中,则该连线被视为穿过了障碍区域,由此被剔除出可行连线路径的集合外。剔除经过障碍连线具体方法为:引入微元法的思想,将所有的线段路径分割成多个小线段,则每个小线段所经过的地图范围用一个小矩形包围,如图8所示,该矩形的四个角点的坐标表示分别如下式:
式中,(x3,1,y3,1),(x3,2,y3,2),(x3,3,y3,3),(x3,4,y3,4)分别为小矩形的四个定角坐标;为向下取整函数,为向上取整函数;dir为连线与地图水平轴的夹角,且dir=arctan((y2-y1)/(x2-x1));d表示每段小线段的长度,且存在关系d=D/N,N=D/0.5。因此只需通过判断小矩形的四个顶角坐标点是否存在位于障碍区域的坐标,便可进一步判断出该连线路径是否穿过了障碍区域,即是否属于可行路径集合。若穿过,则应删除,结果如图9(a)和图9(b)所示。
4)以起点作为搜索出发点,终点作为搜索结束点,采用动态对比的方法比较各搜索路径的距离大小,从中搜索出一条从起点到终点的最优路径。
在步骤4)中,以下内容以图10为例进行详细说明。将步骤3)中选取出的可行连线路径放入Openlist中,如图11中的“Openlist”所示。将起点到所有标记角点的连线路径的长度放入Distancelist中,如图12(a),使用动态对比的方法,比较标记角点t到达表记角点k的多条可行路径的距离,决策得到起点到标记角点k的最短距离,如下式所示:
其中,fs,k表示将要搜索的从起点到标记角点k的最短可行路径距离,fs,t表示已决策的从起点到标记点t的最短可行路径距离,D(t,v,k)表示从标记角点t经标记角点v到标记角点k的可行路径距离,D(t,u,k)表示从标记角点t经标记角点u到标记角点k的可行路径距离(从标记角点t到标记角点k的可行路径可能需要通过多个其他标记角点才能到达,此处举例为经过一个其他标记角点即可到达)。fs,q表示起点到标记角点q无可行路径,其中,k,t,v,u,q为包括终点及所有标记角点的集合U内任意一标记角点。
接着对Distancelist进行更新,同时更新记录到达每个角点最优路径的父节点表Flaglsit,并将已确认为局部最优的连线路径从Openlist中剔除,直至Openlist为空。以图10为例,起点到标记角点1距离为3.16,到标记角点2距离为6.08,到标记点角3距离为∞(即无法直接到达),到标记角点4距离为∞,到标记角点5距离为8.06,到标记角点6的距离为5.10,到终点的距离为∞,如图12(a)中的“起点的Distancelist”所示。通过路径距离的第一次比较,起点可以通过标记角点1到达标记角点2再到达标记角点3,其路径距离为“起点—1”、“1—2”与“2—3”连线路径的距离之和,即3.16+3+4.47=10.63;起点还可以通过标记角点2到达标记角点3,其路径距离为“起点—2”与“2—3”连线路径的距离之和,即6.08+4.47=10.55,而10.55<10.63<∞,故起点的Distancelist中到达标记角点3的距离更新为10.55,如图12(b)中“起点的Distancelist第一次更新”所示。标记角点3的Flaglist更新为“起点—2—3”,如图13中“标记角点3的Flaglist”所示,Openlist中将“起点—1”、“1—2”、“起点—2”与“2—3”剔除,因为这三条连线路径已经比较完毕,如图11中“Openlist第一次更新”所示。同理,通过比较标记角点i到标记点j的路径距离,并更新其Distancelist和Flaglist,直至Openlist为空,如图11中“Openlist更新完毕”所示。最终从终点的Flaglist中得到起点到终点的最短路径连接点记录,如图13的“终点的Flaglist”所示,由此得到起点到终点的最短路径为“起点—5”和“5—终点”的连线路径,其最短路径距离为8.06+4.12=12.18。
本发明通过从经过障碍角点区域的可行连线路径中搜索出一条最优路径,对于图2(a)的障碍地图,本发明搜索最优路径所消耗的时间比Dijkstra算法更少,搜索范围比A*算法和Dijkstra算法更小,搜索得到的最优路径距离为最短,测试“迷宫型”障碍地图的相关数据如表1所示;对于图2(b)的障碍地图,本发明搜索最优路径所消耗的时间比Dijkstra算法更少,搜索范围比Dijkstra算法更小,搜索得到的最优路径距离为最短,测试“陷阱型”障碍地图的相关数据如表2所示。与PRM算法相比,本发明提出的算法的标记角点是随独立的障碍连通区域的角点的数量而生成的,并不是随机生成的,对于复杂障碍尤其是狭缝障碍等均可搜索出一条最优路径,并且可以保证该路径为最短路程。
表1
表2
表1和2为本发明算法与A*算法和Dijkstra算法对于图2的障碍地图搜索最优路径的时间、最优路径的距离和搜索范围的相关数据。
下面通过仿真实验结果测试本发明所提算法的有效性,并对本发明的具体实施步骤及过程作进一步地说明,需要强调的是,以下的仿真结果图均为本发明的真实测试结果。
第一步:首先将要搜索的地图进行了栅格化处理,将障碍区域分割成大小相同的栅格,本发明以50×50方阵为例,使用两种不同的栅格类型地图进行算法测试,如图2(a)、图2(b)所示。
第二步:提出角点合力法用于求取地图中的栅格角点区域,该方法类似于物理中力的合成,把栅格分为4种类型,并对栅格区域周围的空白区域进行赋值,分别如图4(a)、图4(b)、图4(c)和图4(d)所示。接着根据合力方向确定标记点,分别如图5(a)、图5(b)、图5(c)和图5(d)所示。对本发明的两幅障碍测试地图,对图2(a)和图2(b)使用角点合力法,赋值和标记结果如图6(a)和图6(b)所示。
第三步:连接图6(a)和图6(b)中的所有的已标记点、起点和终点,引入微元思想,对各连接路线进行片段化处理,判断每一条连线路径是否穿过的障碍区域,剔除穿过障碍区域的连线路径,得到所有的可行路径如图9(a)和图9(b)。该操作可以极大的缩小路径规划的搜索范围,为路径规划节省了时间。
第四步:计算图9(a)和图9(b)中所有的连线路径的长度并存储,然后对路径距离进行动态对比,如图11、图12(a)和图12(b)、图13所示,记录起点到各个标记点的最优路径所要经过的连线路径。最终可以从以结束点为终点的最优路径记录中得到起点到终点的全局最短路径,实验仿真结果如图14(a)和图14(b)所示,相关数据如表1和2所示。
在无障碍情形下,显然两点之间直线路径最短。若障碍不位于这条直线路径上,此时这些障碍不影响最短路径的选取。若障碍位于该直线路径上,则需要绕过该障碍通往终点,且通过几何知识可以证明,如图7所示,只有“贴着”障碍角点的路径才有可能为最短路径。受此启发,本发明提出一种分治确定性路径寻优算法。栅格化处理地图后,采用“分而治之”的思想,将障碍的凸出角点分为四种情形,使用角点合力法来快速确定地图中的障碍区域凸出角点的位置,并标记障碍区域的角点位置,然后连接各标记角点,通过线段微元法,进一步判断连线是否穿过了障碍区域,从而删除穿过障碍的连线。在剩余的可行路线中,对路径距离采用动态对比的方法比较局部线段的长度值,从而搜索出一条从起点到终点的全局最优路径。

Claims (5)

1.分治确定性路径寻优算法,其特征在于包括以下步骤:
1)将需要搜索的环境地图分割成栅格,每个栅格只有空白状态或赋值状态,空白状态表示该栅格通行,赋值状态表示该栅格为障碍区域,无法通行;根据实际已知信息,在相隔并且中间存在赋值栅格的空白栅格中确定起点和终点位置;
2)根据障碍区域的位置及数量,提出一种角点合力法求取栅格障碍的角点位置,并进行标记;
3)将所有的标记角点两两相连,形成多条可行路径;然后将每段可行路径分割为M个线段,每一线段的线长范围用一个矩形包围,接着求出该矩形四个角的坐标值,并对坐标值进行取整处理,把这个矩形的四个角的坐标值与步骤1)中经栅格化处理生成的矩阵坐标进行比较,判断每个坐标点是否位于空白栅格中;若是,则为可行路径,否则为穿过障碍的路径,将其剔除可行路径集合外;
4)以起点作为搜索出发点,终点作为搜索结束点,比较各搜索路径的距离大小,从中搜索出一条从起点到终点的最优路径。
2.如权利要求1所述分治确定性路径寻优算法,其特征在于在步骤1)中,所述需要搜索的环境地图进行栅格化处理,根据障碍复杂程度及要求精度,设置所要产生地图矩阵的大小,生成一个N×N的方阵,将环境地图分割成多个规则的正方形栅格,其中空白栅格赋值为“0”,表示该区域通行,而障碍栅格赋值为“1”,表示该区域禁止通行;所述空白状态或赋值状态对应于计算机中的“0”和“1”两种状态。
3.如权利要求1所述分治确定性路径寻优算法,其特征在于在步骤2)中,所述栅格障碍的角点位置使用角点合力法计算并标记,所述角点合力法的计算步骤为:
(1)对单位栅格障碍周围的上、下、左、右四个方向的栅格可行区域加权赋值1;
(2)当单位障碍的左或右和上或下的可行区域块都具有加权值1时,单位障碍被定义为障碍连通区域的角点位置;
(3)当单位障碍的左方和上方空白栅格的权值都为1时,标记该单位障碍的左上角的空白栅格;同理,当单位障碍的右方和上方空白栅格的权值都为1时,标记该单位障碍的右上角的空白栅格;当单位障碍的右方和下方空白栅格的权值都为1时,标记该单位障碍的右下角的空白栅格;当单位障碍的左方和下方空白栅格的权值都为1时,标记该单位障碍的左下角的空白栅格,通过力的合成公式:
规定力的合成方向指向障碍的角点处,其中,F为分力F1和分力F2的合力,α为两个分力的夹角。
4.如权利要求1所述分治确定性路径寻优算法,其特征在于在步骤3)中,所述将所有的标记角点两两相连,形成多条可行路径,是连接各标记角点、起点和终点,产生多条路径;接着剔除经过障碍的连线路径,生成可行路径集合;引入微元的思想,将每一段连线微元化成M个线段,每一线段用一个矩形来等同其线长的范围;然后求出该矩形四个角的坐标值,并对坐标值进行取整处理,把矩形的四个角的坐标值与步骤1)中经栅格化处理生成的矩阵坐标进行比较,若这四个对角坐标中有一个坐标点位于障碍区域中,则该连线被视为穿过障碍区域,由此被剔除出可行连线路径的集合外;剔除经过障碍连线的具体方法为:引入微元法的思想,将所有的线段路径分割成多个线段,则每个线段所经过的地图范围用一个矩形包围,该矩形的四个角点的坐标表示如下:
式中,(x3,1,y3,1),(x3,2,y3,2),(x3,3,y3,3),(x3,4,y3,4)分别为矩形的四个定角坐标;为向下取整函数,为向上取整函数;dir为连线与地图水平轴的夹角,且dir=arctan((y2-y1)/(x2-x1));d表示每段小线段的长度,且存在关系d=D/N,N=D/0.5;因此只需通过判断矩形的四个顶角坐标点是否存在位于障碍区域的坐标,判断出该连线路径是否穿过障碍区域,即是否属于可行路径集合;若穿过,则应删除。
5.如权利要求1所述分治确定性路径寻优算法,其特征在于在步骤4)中,所述比较各搜索路径的距离大小采用动态对比的方法。
CN201811332691.XA 2018-11-09 2018-11-09 分治确定性路径寻优算法 Expired - Fee Related CN109443363B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811332691.XA CN109443363B (zh) 2018-11-09 2018-11-09 分治确定性路径寻优算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811332691.XA CN109443363B (zh) 2018-11-09 2018-11-09 分治确定性路径寻优算法

Publications (2)

Publication Number Publication Date
CN109443363A true CN109443363A (zh) 2019-03-08
CN109443363B CN109443363B (zh) 2020-10-16

Family

ID=65551574

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811332691.XA Expired - Fee Related CN109443363B (zh) 2018-11-09 2018-11-09 分治确定性路径寻优算法

Country Status (1)

Country Link
CN (1) CN109443363B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111906779A (zh) * 2020-06-30 2020-11-10 珠海市一微半导体有限公司 一种越障结束判断方法、越障控制方法、芯片及机器人
CN112734878A (zh) * 2020-12-31 2021-04-30 南昌工学院 大型栅格地图两点之间连通性检测方法、设备及存储介质
CN112860829A (zh) * 2021-01-27 2021-05-28 广州朗国电子科技有限公司 一种地图指引系统、地图指引方法及存储介质
CN113108803A (zh) * 2021-04-12 2021-07-13 北京佰能盈天科技股份有限公司 针对双轴定位系统的路径规划方法及设备
CN113219975A (zh) * 2021-05-08 2021-08-06 珠海市一微半导体有限公司 一种路线优化方法、路径规划方法、芯片及机器人
CN113552881A (zh) * 2021-07-15 2021-10-26 浙江工业大学 一种用于神经网络训练的多路径规划数据集生成方法
CN113720344A (zh) * 2021-08-30 2021-11-30 深圳市银星智能科技股份有限公司 路径搜寻方法、装置、智能设备及存储介质
CN113804209A (zh) * 2021-05-27 2021-12-17 中国人民解放军战略支援部队信息工程大学 一种四角格网高精度长距离越野路径规划方法
CN115238525A (zh) * 2022-09-16 2022-10-25 广东工业大学 一种用于行人仿真客流组织的可行路径搜索方法
CN115437368A (zh) * 2022-06-02 2022-12-06 珠海云洲智能科技股份有限公司 一种救援路径确定方法、装置、救援设备及可读存储介质
CN116303479A (zh) * 2023-05-22 2023-06-23 浙大宁波理工学院 一种用于轨迹数据的实时压缩方法
CN113720344B (zh) * 2021-08-30 2024-06-04 深圳银星智能集团股份有限公司 路径搜寻方法、装置、智能设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060149465A1 (en) * 2004-12-30 2006-07-06 Samsung Electronics Co., Ltd. Method and apparatus for moving in minimum cost path using grid map
WO2009005188A1 (en) * 2007-07-03 2009-01-08 Electronics And Telecommunications Research Institute Path search method
JP2011227807A (ja) * 2010-04-22 2011-11-10 Toyota Motor Corp 経路探索システム、経路探索方法、及び移動体
CN106444755A (zh) * 2016-09-22 2017-02-22 江苏理工学院 基于改进遗传算法的移动机器人路径规划方法及系统
CN106708054A (zh) * 2017-01-24 2017-05-24 贵州电网有限责任公司电力科学研究院 结合地图栅格与势场法避障的巡检机器人路径规划方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060149465A1 (en) * 2004-12-30 2006-07-06 Samsung Electronics Co., Ltd. Method and apparatus for moving in minimum cost path using grid map
WO2009005188A1 (en) * 2007-07-03 2009-01-08 Electronics And Telecommunications Research Institute Path search method
JP2011227807A (ja) * 2010-04-22 2011-11-10 Toyota Motor Corp 経路探索システム、経路探索方法、及び移動体
CN106444755A (zh) * 2016-09-22 2017-02-22 江苏理工学院 基于改进遗传算法的移动机器人路径规划方法及系统
CN106708054A (zh) * 2017-01-24 2017-05-24 贵州电网有限责任公司电力科学研究院 结合地图栅格与势场法避障的巡检机器人路径规划方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
C. SARANYA等: "Real Time Evaluation of Grid Based Path Planning Algorithms: A comparative study", 《IFAC PROCEEDINGS VOLUMES》 *
李建建 等: "基于Moore模型统计路径生成算法的实现", 《南昌航空大学学报(自然科学版)》 *
罗德林 等: "基于势场蚁群算法的机器人路径规划", 《系统工程与电子技术》 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111906779B (zh) * 2020-06-30 2022-05-10 珠海一微半导体股份有限公司 一种越障结束判断方法、越障控制方法、芯片及机器人
CN111906779A (zh) * 2020-06-30 2020-11-10 珠海市一微半导体有限公司 一种越障结束判断方法、越障控制方法、芯片及机器人
CN112734878B (zh) * 2020-12-31 2023-06-20 南昌工学院 大型栅格地图两点之间连通性检测方法、设备及存储介质
CN112734878A (zh) * 2020-12-31 2021-04-30 南昌工学院 大型栅格地图两点之间连通性检测方法、设备及存储介质
CN112860829A (zh) * 2021-01-27 2021-05-28 广州朗国电子科技有限公司 一种地图指引系统、地图指引方法及存储介质
CN113108803A (zh) * 2021-04-12 2021-07-13 北京佰能盈天科技股份有限公司 针对双轴定位系统的路径规划方法及设备
CN113219975A (zh) * 2021-05-08 2021-08-06 珠海市一微半导体有限公司 一种路线优化方法、路径规划方法、芯片及机器人
CN113219975B (zh) * 2021-05-08 2024-04-05 珠海一微半导体股份有限公司 一种路线优化方法、路径规划方法、芯片及机器人
CN113804209B (zh) * 2021-05-27 2024-01-09 中国人民解放军战略支援部队信息工程大学 一种四角格网高精度长距离越野路径规划方法
CN113804209A (zh) * 2021-05-27 2021-12-17 中国人民解放军战略支援部队信息工程大学 一种四角格网高精度长距离越野路径规划方法
CN113552881A (zh) * 2021-07-15 2021-10-26 浙江工业大学 一种用于神经网络训练的多路径规划数据集生成方法
CN113552881B (zh) * 2021-07-15 2024-03-26 浙江工业大学 一种用于神经网络训练的多路径规划数据集生成方法
CN113720344A (zh) * 2021-08-30 2021-11-30 深圳市银星智能科技股份有限公司 路径搜寻方法、装置、智能设备及存储介质
CN113720344B (zh) * 2021-08-30 2024-06-04 深圳银星智能集团股份有限公司 路径搜寻方法、装置、智能设备及存储介质
CN115437368A (zh) * 2022-06-02 2022-12-06 珠海云洲智能科技股份有限公司 一种救援路径确定方法、装置、救援设备及可读存储介质
CN115437368B (zh) * 2022-06-02 2023-08-29 珠海云洲智能科技股份有限公司 一种救援路径确定方法、装置、救援设备及可读存储介质
CN115238525B (zh) * 2022-09-16 2023-04-18 广东工业大学 一种用于行人仿真客流组织的可行路径搜索方法
CN115238525A (zh) * 2022-09-16 2022-10-25 广东工业大学 一种用于行人仿真客流组织的可行路径搜索方法
CN116303479B (zh) * 2023-05-22 2023-08-22 浙大宁波理工学院 一种用于轨迹数据的实时压缩方法
CN116303479A (zh) * 2023-05-22 2023-06-23 浙大宁波理工学院 一种用于轨迹数据的实时压缩方法

Also Published As

Publication number Publication date
CN109443363B (zh) 2020-10-16

Similar Documents

Publication Publication Date Title
CN109443363A (zh) 分治确定性路径寻优算法
US6604044B1 (en) Method for generating conflict resolutions for air traffic control of free flight operations
CN108871368B (zh) 一种高精度地图车道横向拓扑关系的构建方法、系统及存储器
CN109115226B (zh) 基于跳点搜索的多机器人冲突避免的路径规划方法
CN111024092B (zh) 一种多约束条件下智能飞行器航迹快速规划方法
CN103901892B (zh) 无人机的控制方法及系统
CN106525047B (zh) 一种基于floyd算法的无人机路径规划方法
CN107402018B (zh) 一种基于连续帧的导盲仪组合路径规划方法
CN108151751A (zh) 一种基于高精度地图和传统地图结合的路径规划方法及装置
CN106202744A (zh) 海上风电场集电系统避障路径优化方法和系统
CN110333659B (zh) 一种基于改进a星搜索的无人驾驶汽车局部路径规划方法
CN108549378A (zh) 一种基于栅格地图的混合路径规划方法和系统
CN104615138A (zh) 一种划分移动机器人室内区域动态覆盖方法及其装置
KR100336597B1 (ko) 경로 탐색 장치
US20150262481A1 (en) System and method to determine occurrence of platoon
CN106125764A (zh) 基于a*搜索的无人机路径动态规划方法
CN105869512B (zh) 多信息的混杂度量地图建图方法和装置
AU2007280774B2 (en) Method for determining a route for an underwater vehicle
CN104464379B (zh) 基于分段匹配的航行计划与雷达航迹关联方法及系统
CN108180911A (zh) 一种agv自动生成修正路径方法
Yang et al. An enhanced weight-based topological map matching algorithm for intricate urban road network
JP5670078B2 (ja) 経路案内システム
KR100967931B1 (ko) 지번 검색 네비게이션 시스템 및 지번 검색 방법
JP2007003826A (ja) 地図データ出力装置、地図データ出力プログラム、路線図出力装置、路線図出力システムおよび路線図出力プログラム
CN105698796A (zh) 一种多机器人调度系统的路径搜索方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201016

Termination date: 20211109

CF01 Termination of patent right due to non-payment of annual fee