CN113189988A - 一种基于Harris算法与RRT算法复合的自主路径规划方法 - Google Patents
一种基于Harris算法与RRT算法复合的自主路径规划方法 Download PDFInfo
- Publication number
- CN113189988A CN113189988A CN202110428987.7A CN202110428987A CN113189988A CN 113189988 A CN113189988 A CN 113189988A CN 202110428987 A CN202110428987 A CN 202110428987A CN 113189988 A CN113189988 A CN 113189988A
- Authority
- CN
- China
- Prior art keywords
- node
- path
- algorithm
- topological
- grid map
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000000203 mixture Substances 0.000 title claims description 6
- 238000010408 sweeping Methods 0.000 claims abstract description 31
- 238000001514 detection method Methods 0.000 claims abstract description 18
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000000605 extraction Methods 0.000 claims abstract description 10
- 230000003044 adaptive effect Effects 0.000 claims description 22
- 239000004576 sand Substances 0.000 claims description 6
- 230000004888 barrier function Effects 0.000 claims description 5
- 239000002131 composite material Substances 0.000 claims 3
- 230000008859 change Effects 0.000 abstract description 4
- 238000005070 sampling Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 206010063385 Intellectualisation Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/02—Control of position or course in two dimensions
- G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
- G05D1/0212—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
- G05D1/0221—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory involving a learning process
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/02—Control of position or course in two dimensions
- G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
- G05D1/0257—Control of position or course in two dimensions specially adapted to land vehicles using a radar
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种基于Harris算法与RRT算法复合的自主路径规划方法,其步骤包括:1、利用骨架提取算法得到栅格地图中的可行路径的轮廓骨架;2、过角点检测算法得到轮廓骨架中变化较大的特征点,在对特征点进行自适应阈值处理后,得到拓扑节点;3、通过自适应阈值来设定拓扑节点之间的距离来进行分步优化;4、进行分步优化时,通过椭圆采样限制其采样空间来提高路径的优化效率以达到快速规划一条可行路径并优化出一条较短的路径。本发明能将Harris角点检测算法和RRT算法相结合来进行路径规划,从而能实现扫地机器人在室内大场景下能够快速的找到一条可行路径。
Description
技术领域
本发明涉及机器人自动导引路径规划方法领域,具体是一种基于Harris算法与RRT算法复合的自主路径规划方法。
背景技术
如今,智能机器人的应用已经逐渐成为主流,如智能停车场、智能家居等,其中室内家居逐渐转向智能化的发展趋势更为明显,用于室内家居的各种机器人如扫地机器人等已经成为众多学者研究的热点课题。路径规划作为机器人等智能化的核心部分,目前主要有两种规划方式:基于环境已知的全局路径规划和基于环境未知的局部路径规划方法。而针对基于环境已知的全局路径规划对环境变化的适应性较低,尤其是针对大场景下如果特定路径的起点和终点,往往机器人需要很长的时间才能到达所设定的终点,并且其规划的路径往往比较复杂,导致无法及时对室内突发状况进行处理,以及无法及时处理常规情景。
发明内容
本发明的目的是提供一种基于Harris算法与RRT算法复合的自主路径规划方法,以解决现有技术全局路径规划方法存在对环境变化适应性低的问题。
为了达到上述目的,本发明所采用的技术方案为:
一种基于Harris算法与RRT算法复合的自主路径规划方法,其特征在于:包括以下步骤:
步骤1、以目标场地的外接矩形的任意一个顶角作为原点o,以所述原点o相连的两条边分别为x轴和y轴,从而建立二维的直角坐标系oxy;
在所述直角坐标系oxy中,将所述目标场地划分成M×N的二维的栅格地图BM×N如公式(1)所示
其中,M表示所述栅格地图BM×N的横向栅格数,N表示栅格地图BM×N的纵向栅格数;GMn表示所述栅格地图BM×N上坐标为(nx,ny)的坐标点被占用的状态,M>nx>0,N>ny>0,若GMn=0表示所述的坐标点为空闲点,若GMn=1表示所述的坐标点为障碍点;
在栅格地图BM×N中设置坐标为(Sx,Sy)的起始点nodes,以及坐标为(gx,gy)的目标点nodeg;
步骤2、采用骨架提取算法对步骤1得到的栅格地图BM×N进行计算,经过骨架提取算法计算后,栅格地图BM×N中的连通区域分别被细化成细线,由此得到包含栅格地图中所有可行路径的多条细线;
步骤3、采用Harris角点检测算法对细化出来的各个可行路径分别进行计算,提取得到各个可行路径对应的特征点nodec,然后对各个特征点nodec进行自适应阈值处理,得到各自对应的拓扑节点nodet;
步骤4、根据拓扑节点nodet,建立拓扑结构地图mapt;
步骤5、采用随机快速搜索树算法,首先把拓扑节点nodet添加到随机快速搜索树算法的初始树Treei中,然后基于随机快速搜索树算法进行初始路径的搜寻与查找;
步骤6、根据自适应阈值公式(3)和公式(4)来对由拓扑节点所构成的初始路径进行局部优化,公式(3)、(4)如下:
Threshold2=μ*(M+N)/4 (3)
其中:Threshold2是一个自适应阈值,作用是为对由拓扑节点所构成的初始路径进行分段,当拓扑节点之间的距离大于自适应阈值Threshold2时,其下面的拓扑节点就变成下一段路径;
μ是根据栅格地图的长、宽和栅格地图大小的数量级来定义的自适应参数,M和N为栅格地图的长和宽;
a=min{M,N}为栅格地图中较短的边,b=max{M,N}为栅格地图中较长的边,ε为栅格地图大小的数量级,当栅格地图较短的边a和较长的边b相等时令μ为1,当栅格地图较短的边a和较长的边b不等时,令μ为公式(4);
通过Threshold2这个阈值把初始路径分段,这些分段的路径为局部路径;
步骤7、将由拓扑节点所构成的初始路径根据公式(2)进行自适应阈值处理后对其进行局部优化,每一段局部路径优化时都被限制在了椭圆内;
在每段局部路径规划中调用Informed RRT*算法,从所构成每段局部路径规划中的第一个拓扑节点nodet1向最后一个拓扑节点站点nodetl进行扩展搜索,利用公式(5)计算各个相邻拓扑节点之间的距离di,其中如果通过Informed RRT*算法扩展到的节点nodee所得到的路径比由拓扑节点所构成的路径短,那么则用nodee代替相应的拓扑节点,公式(5)如下:
公式(5)中,实际路径函数值Ld代表的是相邻两个节点之间的距离,其中每一段优化后的局部路径的长度可由公式(6)表示:
Li=Ld1+Ld2+…+Lde+…+Ldl (6)
其中,Ld1代表第一个拓扑节点nodet1和第二个拓扑节点nodet2;Lde代表通过Informed RRT*算法扩展到的节点nodee和相应的拓扑节点之间的距离;Ldl代表最后一个拓扑节点nodetl和其相邻的拓扑节点之间的距离;
步骤8、将所述每段优化后的路径Li按照拼接规则进行拼接后形成新的路径{L1,L2,L3,…,Lm,…,Ln-1,Ln},由此得到优化后的全局路径La。
所述的Harris算法与RRT算法复合的自主路径规划方法,其特征在于:步骤3中,基于公式(2)对特征点nodec进行自适应阈值处理,公式(2)如下:
公式(2)中,Threshold为自适应节点阈值,经过骨架提取后的栅格地图对其进行角点检测算法会得到很多特征点nodec,为了消除冗余的节点通过自适应阈值Threshold来消除冗余的节点;消除冗余的特征点后还剩下的特征点为nodet,ε为栅格地图大小的数量级。
所述的基于Harris算法与RRT算法复合的自主路径规划方法,其特征在于:步骤5中,初始树Treei按如下规则生成:
(5.1)、所有由拓扑节点nodet所组成的闭合环路生成为一颗初始树;
(5.2)、在交叉支路上的拓扑节点可以有多个父节点,定义为parent1,parent2,…,parentn,以这样的方式生成一颗初始树;
(5.3)、当新节点noden和最近的拓扑节点nodetn的子节点的连线穿过障碍物时,令起点为最近的拓扑节点的父节点,然后把起点插入到初始树Treei中;
(5.4)、当新节点noden和最近的拓扑节点nodetn的子节点的连线没有穿过障碍物时,令起点为最近的拓扑节点的子节点的父节点,同样把起点插入到初始树Treei中。
所述的一种基于Harris算法与RRT算法复合的自主路径规划方法,其特征在于:还包括步骤9:判断拼接后的全局路径La最后一个节点是否为目标点nodeg,若是,则表示路径规划完成,且所述全局路径La为扫地机器人行驶的完整路径,否则返回步骤7顺序执行。
本发明提出一种基于Harris算法与RRT算法复合的自主路径规划方法,通过RRT算法在全局路径规划上的速度优势并且根据初始路径进行优化让其更好的适应扫地机器人的动力学模型,从而提高了大场景下室内环境的路径规划的效率和对不同室内环境的适应性。
与现有技术相比,本发明的优点为:
1.本发明通过将Harris角点检测算法与RRT算法相结合的方式,通过Harris角点检测算法对已经进行骨架处理后的栅格地图提取出的拓扑节点,可以有效地引导RRT算法的扩展方向。防止其进行盲目无效的扩展,可以大大的缩短路径规划的时间,从而可以有效处理固定室内环境下的路径规划时间长和大场景下室内导航算法难以自适应的问题。
2、本发明通过使用Informed RRT*算法,并采取对整个路径分段局部进行优化的方式有效地提高了RRT算法所规划出的路径的曲折性,并且降低了规划路径的长度,从而提高了路径规划的柔性,更好的适应了扫地机器人的动力学模型,降低了室内环境对于处理设备高性能的要求,降低了能源的消耗。
3.本发明通过将Harris角点检测算法与RRT算法相结合的方式来进行路径规划,并对规划的路径进行分段优化,相比于现有的扫地机器人都是遇到障碍物进行旋转,没有进行有效地路径规划,本发明有效的提高了对于静态环境的路径规划的时效性、自适应性、高效性。
附图说明
图1为本发明实施例根据室内环境所仿真的真实场景。
图2为本发明实施例根据室内环境所建立直角坐标系oxy。
图3为本发明实施例根据Harris角点检测算法生成的拓扑结构图。
图4a为本发明实施例中通过Informed RRT*算法扩展新节点的过程。
图4b为本发明实施例中为新节点选择路径代价最小的父节点。
图5为本发明实施例中新节点寻找子节点的过程
图6为本发明实施例在静态环境时根据Harris角点检测算法与RRT算法搜索并生成的初始路径。
图7为本发明实施例在对初始路径进行分段局部路径优化后所得到的适应于机器人动力学模型的路径。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本实例中以室内扫地机器人为例,自重1.5KG,正常行驶速度为1.1m/s,急刹车以及制动时间t约为2s。总体按如下过程实现Harris角点检测算法与RRT算法复合的自主路径规划方法:
1.根据试验场地的外接矩形确定栅格地图的大小;
2.根据室内环境的栅格地图进行骨架提取,然后根据骨架提取后的栅格地图结果,对其进行Harris角点检测得到相应的特征点。通过自适应阈值消除冗余的特征点,其余剩下的称为拓扑节点,构成拓扑结构图。
3.根据Harris角点检测算法与RRT算法快速得到适合此时室内环境的初始路径。
4.采用自适应阈值的方式将由拓扑节点所组成的初始路径进行分段局部优化。
5.将局部路径以拼接规则拼接,并取代之前的全局路径。
具体地说,该自主路径规划方法,是应用于由室内环境较为复杂、若干个障碍物和扫地机器人组成的智能家居系统中。如图1所示,在整个室内场景中,有床柜子还有桌子垃圾桶分布在室内的栅格地图内。栅格地图的每个坐标点都包含坐标属性和状态属性。坐标属性值代表了该坐标点相较于整个室内环境的相对位置,状态属性值代表了室内环境中该位置的是否被障碍物所占用。
本实施例具体过程如下:
步骤1:如图2所示,以仿真的室内环境的外接矩形的任意一个顶角作为原点o,以原点o相连的两条边分别为x轴和y轴,从而建立直角坐标系oxy。
在所述直角坐标系oxy中,将所述目标场地划分成M×N的二维的栅格地图BM×N如公式(1)所示
其中,M表示所述栅格地图BM×N的横向栅格数,N表示栅格地图BM×N的纵向栅格数;GMn表示所述栅格地图BM×N上坐标为(nx,ny)的坐标点被占用的状态,M>nx>0,N>ny>0,若GMn=0表示所述的坐标点为空闲点,若GMn=1表示所述的坐标点为障碍点。
在栅格地图BM×N中设置坐标为(Sx,Sy)的起始点nodes,以及坐标为(gx,gy)的目标点nodeg。
扫地机器人接收到主人发送的清扫任务,任务内容:从扫地机器人的当前位置坐标为(Sx,Sy)的起始点nodes到坐标为(gx,gy)的终点nodeg对地上的垃圾进行清理。起始点nodes与目标点nodeg的相对距离大于两个RRT算法中扩展步长的距离,仅有两个扩展步长的路径规划空间太小,无法给扫地机器人提供足够的应急反应时间;以扫地机器人的当前车头方向作为起始点nodes的路径搜索方向。
步骤2:采用骨架提取算法对步骤1得到的栅格地图BM×N进行计算,经过骨架提取算法计算后,栅格地图BM×N中的连通区域分别被细化成细线,由此得到包含栅格地图中所有可行路径的多条细线。
步骤3:采用Harris角点检测算法对细化出来的各个可行路径分别进行计算,提取得到各个可行路径对应的特征点nodec,然后对各个特征点nodec进行自适应阈值处理,得到各自对应的拓扑节点nodet。
由于所提取出的特征点nodec在栅格地图上的分布不均,所以本实施例根据公式(2)对特征点nodec进行自适应阈值处理,减少冗余的特征点nodec。经过自适应阈值处理后的特征点为拓扑节点nodet,这样可以为接下来扫地机器人进行RRT算法进行路径搜索时,提供了有效地扩展方向,提高了扫地机器人的效率,大大的减少了寻路时间。公式(2)如下:
公式(2)中,Threshold为自适应节点阈值,经过骨架提取后的栅格地图对其进行角点检测算法会得到很多特征点nodec。为了消除冗余的节点通过自适应阈值Threshold来消除冗余的节点。消除冗余的特征点后还剩下的特征点为nodet,ε为栅格地图大小的数量级。
步骤4:在得到拓扑节点nodet后,根据拓扑节点nodet,建立拓扑结构地图mapt,如图3所示。
步骤5:根据扫地机器人得到的任务指令由起始点nodes到终点nodeg。然后其采用随机快速搜索树(RRT)算法进行路径的搜寻与查找,具体的首先把根据对室内环境的栅格地图进行Harris角点检测算法后所得到的拓扑节点添加到随机快速搜索树(RRT)算法的初始树Treei中,初始树Treei按如下规则生成:
(5.1)、所有由拓扑节点nodet所组成的闭合环路生成为一颗初始树
(5.2)、在交叉支路上的拓扑节点可以有多个父节点,定义为parent1,parent2,…,parentn,以这样的方式生成一颗初始树
(5.3)、当任务所设定的起点nodes或终点nodeg与其最近的拓扑节点nodetn的子节点的连线穿过障碍物时,令起点nodes或终点nodeg为最近的拓扑节点的父节点,然后把起点nodes和终点nodeg插入到初始树Treei中。
(5.4)、当任务所设定的起点nodes或终点nodeg与其最近的拓扑节点nodetn的子节点的连线没有穿过障碍物时,令起点为最近的拓扑节点的子节点的父节点,同样把起点nodes和终点nodeg插入到初始树Treei中。通过这样可以有效地把新节点插入到拓扑结构中。
传统的RRT算法中的每个节点只有一个父节点,这样节点不能和拓扑结构相结合,因此本发明提出的算法让拓扑节点可以同时具有多个父节点,这样可以有效的把拓扑结构以树的形式表示出来。
步骤5主要可以概括为首先在初始树中找到距离起点和终点最近的拓扑节点。经过步骤(5.3)和(5.4)的判断,通过本发明所设定的父子关系把起点和终点分别替代最近的拓扑节点,可以有效的把起点和终点加入到拓扑结构中来扩充初始树,这样可以使扫地机器人快速的得到初始路径,其初始路径的结果如图6所示,可以有效减少等待时间,防止需要清扫的东西对地面进行污染,为智能家居系统提高了运行的效率。
步骤6:扫地机器人根据自适应阈值公式(2)和公式(3)来对由拓扑节点所构成的初始路径进行局部优化,在进行局部优化时,可以让扫地机器人一边进行这一段优化好的路径进行行驶,一边对下一段路程进行优化,这样可以加快算法的收敛,同时减少了时间上的浪费,优化后的路径更为平滑,适应扫地机器人的运行,不会出现较大的曲折路径。公式(3)、(4)如下:
Threshold2=μ*(M+N)/4 (3)
Threshold2是一个自适应阈值,目的是为对由拓扑节点所构成的初始路径进行分段,当拓扑节点之间的距离大于自适应阈值Threshold2时,其下面的拓扑节点就变成下一段路径了,通过Threshold2这个阈值可以把初始路径分段,这些分段的路径称为局部路径,然后对每段局部路径通过步骤七进行优化。μ是根据栅格地图的长、宽和栅格地图大小的数量级来定义的自适应参数,M和N为栅格地图的长和宽。
a=min{M,N}为栅格地图中较短的边,b=max{M,N}为栅格地图中较长的边,ε为栅格地图大小的数量级,当栅格地图较短的边a和较长的边b相等时令μ为1,当栅格地图较短的边a和较长的边b不等时,令μ为公式(4);
步骤7:将由拓扑节点所构成的初始路径根据公式(3)和(4)进行自适应阈值处理后对其进行局部优化,每一段局部路径优化时都被限制在了椭圆内,并且优化时都分两个过程:1.为新节点选择代价最小的父节点,2.重新布线原则。具体如下:
1.为新节点选择代价最小的父节点。如图4a和图4b所示,首先根据扫地机器人通过随机采样得到的采样点noder(noderand)作为起点nodes(nodestart)寻找目前随机树中最近的节点noden(nodenear),然后通过Informed RRT*算法中的扩展方式得到了一个新节点nodene(nodenew)。接下来以nodene为圆心,一定半径画一个圆,逐个计算起点到圆内节点的距离与圆内节点到新节点nodene的距离之和,其中最小的距离和就是起点到新节点nodene的路径代价,相应的节点为nodene的父节点。nodes就是扫地机器人目前所在的位置,nodene就是其通过RRT算法扩展出来的一个新节点。扫地机器人需要对这个新扩展出来的节点寻找路径代价最小的父节点。
2.重新布线原则,也就是为扫地机器人所扩展的新节点寻找子节点的过程。在由新节点nodene所构成的圆中,对圆中的任意一个节点,计算以新节点为圆中节点的父节点时,其距离代价是否会降低,如果会,那么将该节点的父节点赋值为新节点,对应的代价更新为新的代价。从图4b中可以得到,节点e距离扫地机器人原来的路径代价为23,如果以节点c为父节点,那么新的路径代价为20,小于原来的路径代价,那么将节点e的父节点改为节点c,结果图5所示:
在每段局部路径规划(这个局部路径规划是通过步骤6的自适应阈值Threshold2对初始路径进行分段得到的)中调用Informed RRT*算法,从所构成每段局部路径规划中的第一个拓扑节点nodet1向最后一个拓扑节点站点nodetl进行扩展搜索,利用公式(5)计算各个相邻拓扑节点之间的距离di,其中如果通过Informed RRT*算法扩展到的节点nodee所得到的路径比由拓扑节点所构成的路径短,那么则用nodee代替相应的拓扑节点,公式(5)如下:
公式(5)中,实际路径函数值Ld代表的是相邻两个节点之间的距离,其中每一段优化后的局部路径的长度可由公式(6)表示
Li=Ld1+Ld2+…+Lde+…+Ldl (6)
其中,Ld1代表第一个拓扑节点nodet1和第二个拓扑节点nodet2。Lde代表通过Informed RRT*算法扩展到的节点nodee和相应的拓扑节点之间的距离。Ldl代表最后一个拓扑节点nodetl和其相邻的拓扑节点之间的距离;
步骤8:扫地机器人将所述每段优化后的路径Li按照拼接规则进行拼接后形成新的路径{L1,L2,L3,…,Lm,…,Ln-1,Ln}得到优化后的全局路径La。
对于距离远的全局路径规划一般耗时过长、效率低下,因此本发明采用Harris角点检测算法与RRT算法,减少了计算的复杂度,从而减少了对室内环境计算设备的高性能要求,降低了实施成本,也提高了计算的效率,对智能家居扫地机器人常发生的不必要的碰撞,以及无效的路径移动有了很大的改善与提高。比如刚刚清理完的区域会有二次重复清理的可能,间接的节约能源的消耗。
所述拼接规则为:按照构成初始路径的拓扑节点的顺序依次将所得到的将所得到的优化后的局部路径按顺序进行拼接,其拼接后的优化路径结果如图7所示。
由于室内环境现场可能出现一定的变化,但是其主要的拓扑节点不会出现较大的改变,如果扫地机器人移动的路线中遇到了障碍物,可以通过机器人上的雷达预知情况采取制动,一般民用较便宜的雷达对障碍物的精准感知距离为3m,发现障碍和系统处理时间约为0.5s,急刹车以及制动时间t约为2s,此时扫地机器人的制动距离为S=v×t=1×2.5=2.5m,距离与障碍物相撞还差0.5m,有效地避免了碰撞,如果家中有较为珍贵的物品,可以有效地减少财产损失,同时立即对其进行RRT算法对路径进行重新搜索,当RRT算法所扩展的节点与拓扑节点吻合时,继续进行步骤5、6、7。
步骤9:判断拼接后的全局路径La最后一个节点是否为目标点nodeg,若是,则表示路径规划完成,且所述全局路径La为扫地机器人行驶的完整路径,任务完成后上报控制中心已经完成,等待下一条任务指令,若是没有任务则在原地等待,否则返回步骤7顺序执行。
本发明所述的实施例仅仅是对本发明的优选实施方式进行的描述,并非对本发明构思和范围进行限定,在不脱离本发明设计思想的前提下,本领域中工程技术人员对本发明的技术方案做出的各种变型和改进,均应落入本发明的保护范围,本发明请求保护的技术内容,已经全部记载在权利要求书中。
Claims (4)
1.一种基于Harris算法与RRT算法复合的自主路径规划方法,其特征在于:包括以下步骤:
步骤1、以目标场地的外接矩形的任意一个顶角作为原点o,以所述原点o相连的两条边分别为x轴和y轴,从而建立二维的直角坐标系oxy;
在所述直角坐标系oxy中,将所述目标场地划分成M×N的二维的栅格地图BM×N如公式(1)所示
其中,M表示所述栅格地图BM×N的横向栅格数,N表示栅格地图BM×N的纵向栅格数;GMn表示所述栅格地图BM×N上坐标为(nx,ny)的坐标点被占用的状态,M>nx>0,N>ny>0,若GMn=0表示所述的坐标点为空闲点,若GMn=1表示所述的坐标点为障碍点;
在栅格地图BM×N中设置坐标为(Sx,Sy)的起始点nodes,以及坐标为(gx,gy)的目标点nodeg;
步骤2、采用骨架提取算法对步骤1得到的栅格地图BM×N进行计算,经过骨架提取算法计算后,栅格地图BM×N中的连通区域分别被细化成细线,由此得到包含栅格地图中所有可行路径的多条细线;
步骤3、采用Harris角点检测算法对细化出来的各个可行路径分别进行计算,提取得到各个可行路径对应的特征点nodec,然后对各个特征点nodec进行自适应阈值处理,得到各自对应的拓扑节点nodet;
步骤4、根据拓扑节点nodet,建立拓扑结构地图mapt;
步骤5、采用随机快速搜索树算法,首先把拓扑节点nodet添加到随机快速搜索树算法的初始树Treei中,然后基于随机快速搜索树算法进行初始路径的搜寻与查找;
步骤6、根据自适应阈值公式(3)和公式(4)来对由拓扑节点所构成的初始路径进行局部优化,公式(3)、(4)如下:
Threshold2=μ*(M+N)/4 (3)
其中:Threshold2是一个自适应阈值,作用是为对由拓扑节点所构成的初始路径进行分段,当拓扑节点之间的距离大于自适应阈值Threshold2时,其下面的拓扑节点就变成下一段路径;
μ是根据栅格地图的长、宽和栅格地图大小的数量级来定义的自适应参数,M和N为栅格地图的长和宽;
a=min{M,N}为栅格地图中较短的边,b=max{M,N}为栅格地图中较长的边,ε为栅格地图大小的数量级,当栅格地图较短的边a和较长的边b相等时令μ为1,当栅格地图较短的边a和较长的边b不等时,令μ为公式(4);
通过Threshold2这个阈值把初始路径分段,这些分段的路径为局部路径;
步骤7、将由拓扑节点所构成的初始路径根据公式(2)进行自适应阈值处理后对其进行局部优化,每一段局部路径优化时都被限制在了椭圆内;
在每段局部路径规划中调用Informed RRT*算法,从所构成每段局部路径规划中的第一个拓扑节点nodet1向最后一个拓扑节点站点nodetl进行扩展搜索,利用公式(5)计算各个相邻拓扑节点之间的距离di,其中如果通过Informed RRT*算法扩展到的节点nodee所得到的路径比由拓扑节点所构成的路径短,那么则用nodee代替相应的拓扑节点,公式(5)如下:
公式(5)中,实际路径函数值Ld代表的是相邻两个节点之间的距离,其中每一段优化后的局部路径的长度可由公式(6)表示:
Li=Ld1+Ld2+…+Lde+…+Ldl (6)
其中,Ld1代表第一个拓扑节点nodet1和第二个拓扑节点nodet2;Lde代表通过InformedRRT*算法扩展到的节点nodee和相应的拓扑节点之间的距离;Ldl代表最后一个拓扑节点nodetl和其相邻的拓扑节点之间的距离;
步骤8、将所述每段优化后的路径Li按照拼接规则进行拼接后形成新的路径{L1,L2,L3,...,Lm,...,Ln-1,Ln},由此得到优化后的全局路径La。
3.根据权利要求1所述的基于Harris算法与RRT算法复合的自主路径规划方法,其特征在于:步骤5中,初始树Treei按如下规则生成:
(5.1)、所有由拓扑节点nodet所组成的闭合环路生成为一颗初始树;
(5.2)、在交叉支路上的拓扑节点可以有多个父节点,定义为parent1,parent2,...,parentn,以这样的方式生成一颗初始树;
(5.3)、当新节点noden和最近的拓扑节点nodetn的子节点的连线穿过障碍物时,令起点为最近的拓扑节点的父节点,然后把起点插入到初始树Treei中;
(5.4)、当新节点noden和最近的拓扑节点nodetn的子节点的连线没有穿过障碍物时,令起点为最近的拓扑节点的子节点的父节点,同样把起点插入到初始树Treei中。
4.根据权利要求1所述的基于Harris算法与RRT算法复合的自主路径规划方法,其特征在于:还包括步骤9:判断拼接后的全局路径La最后一个节点是否为目标点nodeg,若是,则表示路径规划完成,且所述全局路径La为扫地机器人行驶的完整路径,否则返回步骤7顺序执行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110428987.7A CN113189988B (zh) | 2021-04-21 | 2021-04-21 | 一种基于Harris算法与RRT算法复合的自主路径规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110428987.7A CN113189988B (zh) | 2021-04-21 | 2021-04-21 | 一种基于Harris算法与RRT算法复合的自主路径规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113189988A true CN113189988A (zh) | 2021-07-30 |
CN113189988B CN113189988B (zh) | 2022-04-15 |
Family
ID=76977841
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110428987.7A Active CN113189988B (zh) | 2021-04-21 | 2021-04-21 | 一种基于Harris算法与RRT算法复合的自主路径规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113189988B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113587950A (zh) * | 2021-08-26 | 2021-11-02 | 清华大学 | 自动驾驶汽车静态路径规划方法、装置及存储介质 |
CN114355888A (zh) * | 2021-12-07 | 2022-04-15 | 西安交通大学 | 基于直骨架的室内多级拓扑地图的自动生成方法及系统 |
CN114764245A (zh) * | 2022-04-19 | 2022-07-19 | 北京能工荟智机器人有限责任公司 | 无人巡检机器人及其控制方法 |
CN114898013A (zh) * | 2022-07-15 | 2022-08-12 | 深圳市城市交通规划设计研究中心股份有限公司 | 一种交通等时圈生成方法、电子设备及存储介质 |
WO2023124035A1 (zh) * | 2021-12-28 | 2023-07-06 | 美的集团(上海)有限公司 | 机器人的建图方法、电子设备及非易失性可读存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100174435A1 (en) * | 2009-01-07 | 2010-07-08 | Samsung Electronics Co., Ltd. | Path planning apparatus of robot and method thereof |
CN104898660A (zh) * | 2015-03-27 | 2015-09-09 | 中国科学技术大学 | 一种提高机器人路径规划效率的室内地图构建方法 |
US20170120448A1 (en) * | 2015-10-29 | 2017-05-04 | Korea Institute Of Science And Technology | Robot control system and method for planning driving path of robot |
CN109541634A (zh) * | 2018-12-28 | 2019-03-29 | 歌尔股份有限公司 | 一种路径规划方法、装置和移动设备 |
CN110823241A (zh) * | 2019-11-19 | 2020-02-21 | 齐鲁工业大学 | 基于可通行区域骨架提取的机器人路径规划方法及系统 |
CN111141304A (zh) * | 2019-12-30 | 2020-05-12 | 福州大学 | 一种基于同心圆采样引导rrt算法的路径规划方法 |
CN111220153A (zh) * | 2020-01-15 | 2020-06-02 | 西安交通大学 | 基于视觉拓扑节点和惯性导航的定位方法 |
CN111323037A (zh) * | 2020-02-28 | 2020-06-23 | 武汉科技大学 | 一种移动机器人新型骨架提取的Voronoi路径规划算法 |
CN111857149A (zh) * | 2020-07-29 | 2020-10-30 | 合肥工业大学 | 一种a*算法与d*算法复合的自主路径规划方法 |
CN112577507A (zh) * | 2020-11-04 | 2021-03-30 | 杭州电子科技大学 | 基于哈里斯鹰优化算法的电动汽车路径规划方法 |
-
2021
- 2021-04-21 CN CN202110428987.7A patent/CN113189988B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100174435A1 (en) * | 2009-01-07 | 2010-07-08 | Samsung Electronics Co., Ltd. | Path planning apparatus of robot and method thereof |
CN104898660A (zh) * | 2015-03-27 | 2015-09-09 | 中国科学技术大学 | 一种提高机器人路径规划效率的室内地图构建方法 |
US20170120448A1 (en) * | 2015-10-29 | 2017-05-04 | Korea Institute Of Science And Technology | Robot control system and method for planning driving path of robot |
CN109541634A (zh) * | 2018-12-28 | 2019-03-29 | 歌尔股份有限公司 | 一种路径规划方法、装置和移动设备 |
CN110823241A (zh) * | 2019-11-19 | 2020-02-21 | 齐鲁工业大学 | 基于可通行区域骨架提取的机器人路径规划方法及系统 |
CN111141304A (zh) * | 2019-12-30 | 2020-05-12 | 福州大学 | 一种基于同心圆采样引导rrt算法的路径规划方法 |
CN111220153A (zh) * | 2020-01-15 | 2020-06-02 | 西安交通大学 | 基于视觉拓扑节点和惯性导航的定位方法 |
CN111323037A (zh) * | 2020-02-28 | 2020-06-23 | 武汉科技大学 | 一种移动机器人新型骨架提取的Voronoi路径规划算法 |
CN111857149A (zh) * | 2020-07-29 | 2020-10-30 | 合肥工业大学 | 一种a*算法与d*算法复合的自主路径规划方法 |
CN112577507A (zh) * | 2020-11-04 | 2021-03-30 | 杭州电子科技大学 | 基于哈里斯鹰优化算法的电动汽车路径规划方法 |
Non-Patent Citations (2)
Title |
---|
JONATHAND.GAMMELL: "Informed RRT*:Optimal Incremental Path Planning Focused through an Admissible Ellipsoidal Heuristic", 《2014 IEEE/RSJ INTERNATIONAL CONFERENCE ON INTELLIGENT ROBOTS AND SYSTEMS (IROS)》 * |
张玉伟等: "基于改进Informed-RRT*算法的路径规划研究", 《组合机床与自动化加工技术》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113587950A (zh) * | 2021-08-26 | 2021-11-02 | 清华大学 | 自动驾驶汽车静态路径规划方法、装置及存储介质 |
CN113587950B (zh) * | 2021-08-26 | 2024-04-09 | 清华大学 | 自动驾驶汽车静态路径规划方法、装置及存储介质 |
CN114355888A (zh) * | 2021-12-07 | 2022-04-15 | 西安交通大学 | 基于直骨架的室内多级拓扑地图的自动生成方法及系统 |
CN114355888B (zh) * | 2021-12-07 | 2023-09-26 | 西安交通大学 | 基于直骨架的室内多级拓扑地图的自动生成方法及系统 |
WO2023124035A1 (zh) * | 2021-12-28 | 2023-07-06 | 美的集团(上海)有限公司 | 机器人的建图方法、电子设备及非易失性可读存储介质 |
CN114764245A (zh) * | 2022-04-19 | 2022-07-19 | 北京能工荟智机器人有限责任公司 | 无人巡检机器人及其控制方法 |
CN114898013A (zh) * | 2022-07-15 | 2022-08-12 | 深圳市城市交通规划设计研究中心股份有限公司 | 一种交通等时圈生成方法、电子设备及存储介质 |
CN114898013B (zh) * | 2022-07-15 | 2022-11-22 | 深圳市城市交通规划设计研究中心股份有限公司 | 一种交通等时圈生成方法、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113189988B (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113189988A (zh) | 一种基于Harris算法与RRT算法复合的自主路径规划方法 | |
CN109945873B (zh) | 一种用于室内移动机器人运动控制的混合路径规划方法 | |
CN114510056B (zh) | 室内移动机器人的平稳移动全局路径规划方法 | |
CN110986953B (zh) | 路径规划方法、机器人及计算机可读存储介质 | |
CN110231824B (zh) | 基于直线偏离度方法的智能体路径规划方法 | |
CN112229419B (zh) | 一种动态路径规划导航方法及系统 | |
Sudhakara et al. | Trajectory planning of a mobile robot using enhanced A-star algorithm | |
CN108444490B (zh) | 基于可视图和a*算法深度融合的机器人路径规划方法 | |
CN112650256A (zh) | 一种基于改进双向rrt机器人路径规划方法 | |
CN106371445A (zh) | 一种基于拓扑地图的无人车规划控制方法 | |
CN110006429A (zh) | 一种基于深度优化的无人船航迹规划方法 | |
CN109947120A (zh) | 仓储系统中的路径规划方法 | |
CN114166235A (zh) | 一种基于优化a-star算法的全局动态平滑路径规划方法 | |
CN110231821A (zh) | 多机器人编队的改进自适应零空间行为融合方法 | |
CN115047880A (zh) | 一种未知动态环境下机器人智能路径规划方法 | |
Gu et al. | Path planning for mobile robot in a 2.5‐dimensional grid‐based map | |
Wang et al. | Path Planning for Mobile Robots Based on Improved A* Algorithm | |
CN114986501A (zh) | 一种机械臂路径规划方法、系统及机械臂 | |
CN117553804B (zh) | 路径规划方法、装置、计算机设备和存储介质 | |
CN111829526B (zh) | 一种基于防撞半径的距离地图重构与跳点路径规划方法 | |
CN117124335B (zh) | 一种基于路径标记回溯策略的改进式rrt路径规划方法 | |
Zeng et al. | An efficient path planning algorithm for mobile robots | |
CN113791610B (zh) | 一种移动机器人全局路径规划方法 | |
Xu et al. | Hybrid frontier detection strategy for autonomous exploration in multi-obstacles environment | |
CN114625815A (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 |