CN113345008B - 考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法 - Google Patents
考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法 Download PDFInfo
- Publication number
- CN113345008B CN113345008B CN202110601426.2A CN202110601426A CN113345008B CN 113345008 B CN113345008 B CN 113345008B CN 202110601426 A CN202110601426 A CN 202110601426A CN 113345008 B CN113345008 B CN 113345008B
- Authority
- CN
- China
- Prior art keywords
- point cloud
- point
- obstacle
- characteristic
- dynamic
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1656—Programme controls characterised by programming, planning systems for manipulators
- B25J9/1664—Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning
- B25J9/1666—Avoiding collision or forbidden zones
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Abstract
本发明公开了一种考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,该方法包括:步骤1,通过激光点云数据的几何特性,提取符合预设曲率特征的点云作为特征点,建立同一特征点在两帧相邻时刻点云数据中的匹配关系,并构建代价函数,以轮式机器人位姿为变量构造ICP问题,获得轮式机器人位姿信息;步骤2,检测候选动态障碍物;步骤3,估计动态障碍物状态。本发明不需要结合多传感器,在使用单一传感器的前提下进行动态障碍物检测目的,使系统更加高效,而且安全性更高。
Description
技术领域
本发明涉及激光雷达障碍物检测技术领域,特别是关于一种考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法。
背景技术
现如今,在人工智能高速发展的推动下,机器人技术正从传统工业制造向教育娱乐、医疗服务和勘探勘测方向迅速扩展,甚至也影响着自动驾驶等交叉领域。这些科技应用不断进步使得人们生活越发便捷。在众多驱动方式不同的机器人中,轮式机器人(WMR)由于其机械结构成熟,被广泛应用在了各行各业,例如物流机器人、仓储机器人和扫地机器人等。其中,智能机器人系统便为其核心组成部分,具有重要的研究意义。
智能机器人系统主要由三大技术部分组成,即环境感知、决策规划和运动控制,是一种综合性的智能控制系统。为保证安全,完成作业任务,机器人的避障功能必不可少,如扫地机器人智能避障功能和自动驾驶汽车主动防撞技术等,而这些需求都离不开障碍物检测这一关键技术。
障碍物检测技术通常基于视觉传感器或激光传感器所开展,但由于单目相机在尺度上的模糊性以及双目相机在深度估计方面的范围局限性,不适用于机器人的实时避撞过程。而激光雷达由于测量精度高、测距范围广的特性常被作为主要感知传感器,应用到移动机器人及无人驾驶汽车上。
随着机器人应用场景的复杂性不断提升,在机器人运动过程中不仅需要对障碍物进行准确检测,同时也需要对障碍物的运动状态做进一步预测,根据障碍物的预测轨迹,从而使机器人在路径规划过程中更新安全路线以避免碰撞。根据运动状态不同,可将障碍物分成两类,一是静态障碍物,如建筑物和绿化带等,机器人与静态障碍物的相对速度来源于机器人自身运动;而机器人与动态障碍物的相对速度不仅取决于自身速度,还取决于动态障碍物的运动状态。如何最大限度地估计障碍物状态,以实现最高效安全的避障,仍是一个待研究的课题。现有方案主要有两种:一是结合深度学习技术,在提取出环境语义信息的基础上,判断障碍物是否为可能运动的物体,如行人、车辆,从而进行避障优先级决策。但其无法在不知晓机器人自身运动状态前提下仅利用单帧信息估计出障碍物运动状态,此外,深度学习技术强依赖于其训练数据库,需要并行计算加速,硬件成本高,环境适应性低,并非是障碍物检测技术落地的首要之选。二是将多源传感器信息结合,如激光雷达结合轮式编码器或惯性导航单元在推算出机器人自身运动状态的同时,对障碍物进行有效跟踪,从而进行避障规划。
为降低系统硬件拓扑复杂度,提高系统可靠性,通常希望复用少量传感器实现多目标功能。目前基于单一传感器的动态障碍物识别方法仍不成熟,相关理论仍有待进一步完善。
发明内容
本发明的目的在于提供一种考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法来克服或至少减轻现有技术的上述缺陷中的至少一个。
为实现上述目的,本发明提供一种考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,该方法包括:
步骤1,通过激光点云数据的几何特性,提取符合预设曲率特征的点云作为特征点,建立同一特征点在两帧相邻时刻点云数据中的匹配关系,并构建代价函数,以轮式机器人位姿为变量构造ICP问题,获得轮式机器人位姿信息;
步骤2,根据步骤1获得的轮式机器人的位姿信息,检测候选动态障碍物;
步骤3,根据步骤2检测得到的候选动态障碍物,估计动态障碍物状态。
进一步地,步骤1中,所述提取符合预设曲率特征的点云作为特征点的方法具体包括:
步骤a1,对点云数据进行体素滤波;
步骤b1,利用下式(1)计算点云数据中各点云的曲率值r,并按照曲率值r的大小排序,曲率值最大的点云定义为特征角点,曲率值最小的点云定义为特征平面点:
式中,X(k,i)表示第k帧点云数据Pk中的第i个点云在雷达坐标系中的坐标,X(k,j)表示Pk中的第j个点云在雷达坐标系中的坐标,S表示X(k,i)的邻域。
进一步地,所述特征点为特征角点的情形下,步骤1中,所述建立同一特征点在两帧相邻时刻点云数据中的匹配关系的方法具体包括:
在第k-1帧点云数据Pk-1中寻找Pk中的第a个特征角点X(k,a)对应的特征线段,该线段上具有点云b和点云c,其中,点云b为Pk-1中与特征角点X(k,a)距离最近的点云,点云c为Pk-1中在特征角点X(k,a)相邻扫描线上与特征角点X(k,a)距离最近的点云,点云线距离表示为下式(2):
式中,X(k-1,b)表示Pk-1中的第b个特征角点在雷达坐标系中的坐标,X(k-1,c)表示Pk-1中的第c个特征角点在雷达坐标系中的坐标;
所述特征点为特征平面点的情形下,步骤1中,所述建立同一特征点在两帧相邻时刻点云数据中的匹配关系的方法具体包括:
在Pk-1中寻找第d个特征平面点X(k,d)所对应的特征面,该特征面具有点云e、点云f和点云g,点云e为Pk-1中与特征平面点X(k,d)距离最近的点云,点云f为Pk-1中在点云e同一相邻扫描线上与点云d距离最近的点云,点云g为在相邻扫描线上与特征平面点X(k,d)距离最近的点云,点云面距离表示为下式(3):
式中,X(k-1,e)表示Pk-1中的第e个特征平面点在雷达坐标系中的坐标,X(k-1,f)表示Pk-1中的第f个特征平面点在雷达坐标系中的坐标,X(k-1,g)表示Pk-1中的第g个特征平面点在雷达坐标系中的坐标。
进一步地,
步骤1中,构建代价函数,以轮式机器人位姿为变量构造ICP问题,获得轮式机器人位姿信息的方法具体包括:
W为一个2×2矩阵,对其进行SVD分解,可得式(8):
W=U∑VT (8)
其中,∑为奇异值对角矩阵,U和V为正交矩阵;
步骤21,根据步骤1输出的轮式机器人位姿信息,提取出可能属于动态障碍物点云;
步骤22,分割步骤21提取出的可能属于动态障碍物点云中的地面点云与非地面点云,去除地面点云,得到非地面点云;
步骤23,将整体的非地面点云按照特定策略聚类为多个点云簇,一个点云簇对应一个障碍物;
步骤24,估计步骤23获得的每个障碍物的位置、形状和尺寸信息;
其中,步骤21具体包括:
步骤212,根据最近邻关系,Pk与Pk-1的交集部分确定为静态点云部分;
步骤213,将Pk与Pk-1的非交集部分的点云进行障碍物检测操作,通过观测多帧数据的方式去除属于静态障碍物的目标。
进一步地,步骤23具体包括:
进一步地,步骤24具体包括:
步骤241,将步骤23得到的点云簇的质心定义为该障碍物的坐标;
步骤242,将移除z轴维度后的点云进行随机采样一致性直线分割,将分割出的直线方向作为该障碍物主方向;
步骤243,利用步骤242确定的主方向建立二维坐标系,再根据障碍物所对应点云簇中每个点云的空间信息,计算点云在该坐标系x′、y′两个维度上的极大、极小值,并将极大值与相应的极小值的差值作为障碍物的长和宽。
进一步地,步骤3具体包括:
步骤31,利用候选动态障碍物信息根据多帧数据中障碍物的差异程度计算差异度函数,并建立关联矩阵;
步骤32,根据动态障碍物的关联情况,结合激光里程计计算得到的帧间相对运动,判断障碍物运动状态,从而得到动态障碍物的运动方向及速度。
进一步地,步骤31具体包括:
步骤313,根据步骤312中建立的关联矩阵(11),关联当前帧和上一帧点云数据中的动态障碍物,若最小差异度所对应的Pk-1中元素的下标u>v,则表明发生障碍物丢失;若最小差异度所对应的Pk-1中元素的下标u<v,则表明有新动态障碍物出现。
进一步地,步骤32具体包括:
根据激光雷达传感器的帧率f,则动态障碍物的速度为V=Δd×f,方向为v=tk-tk-1。
本发明由于采取以上技术方案,具有以下优点:1)不需要结合多传感器,在使用单一传感器的前提下进行动态障碍物检测目的,使系统更加高效;2)激光里程计在短时内的位姿估计较为准确,且不会出现零漂现象,也不需额外的标定工作,所以使用激光里程计代替轮速计或惯导可以使系统更简单和稳定,因为轮速计偏差较大,容易导致误匹配,而惯导会出现零漂现象,且需要和激光雷达进行外参标定,同样会使得系统变复杂,不确定因素增多;3)基于轮式机器人的运动特性将障碍物检测问题转化为二维问题,在达到动态避障目的前提下使算法复杂度降低,从而可得到更高的数据处理频率;4)充分利用激光里程计信息,在提取出潜在动态点云基础上再进行动态障碍物检测,大大减少了由分割聚类带来的冗余计算;5)由于障碍物跟踪是在估计了自身位姿的基础上进行的,相比于不考虑自身位姿的跟踪算法,本发明的障碍物跟踪效果更好。
附图说明
图1为动态障碍物检测算法架构示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明实施例提供的考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法包括:
步骤1,通过激光点云数据的几何特性,提取符合预设曲率特征的点云作为特征点,建立同一特征点在两帧相邻时刻点云数据中的匹配关系,并构建代价函数,以轮式机器人位姿为变量构造ICP(Iterative Closest Point;迭代最近点)问题,并使用SVD(Singular Value Decomposition;奇异值分解)分解在雷达坐标系的x和y二维平面解空间中进行迭代,直至位姿收敛,获得轮式机器人的位姿信息。
步骤2,根据步骤1获得的轮式机器人的位姿信息,检测候选动态障碍物。
步骤3,根据步骤2检测得到的候选动态障碍物,估计动态障碍物的运动方向及速度。
本发明实施例利用多帧激光点云数据进行特征点配准,解算机器人位姿,以实现局部里程计功能,进而替代轮式编码器或惯性导航单元。由于轮式机器人仅具有三个自由度,故本发明将激光点云的三维数据经算法处理后转化为二维结果进行描述,进而在达到动态避障目的前提下进一步降低算法复杂度。
在一个实施例中,步骤1中,所述提取符合预设曲率特征的点云作为特征点的方法具体包括:
步骤a1,对点云数据进行体素滤波,以在保证形状特征不变的前提下大幅度降低点云的数据量,为降低算法复杂度提供有利条件。
步骤b1,计算点云数据中各点云的曲率值,比如利用下式(1),并按照曲率值大小排序,曲率值最大的点云定义为特征角点,曲率值最小的点云定义为特征平面点。选择特征角点和特征平面点作为特征点,为后续点云配准做准备。
式中,X(k,i)表示第k帧点云数据Pk中的第i个点云在雷达坐标系中的坐标,描述为:X(k,i)=[x y z]T;X(k,j)表示Pk中的第j个点云在雷达坐标系中的坐标;S表示X(k,i)的邻域,其可以通过KD树等现有方法快速获得。
进一步地,步骤1中,通过分区域提取所述特征点,这样可以保证特征点分布均匀,即将点云云数据划分为不同区域,在每个区域内提取一定阈值数量的特征点。此外,在选取特征点过程中,应当避免将符合以下三种情况的点云被选作为特征点:1)局部平面与激光束大致平行的点云,此类点云可靠性较低,其中的局部平面可以理解为周围(一定邻域内)的点云与该点云共面而形成的平面;2)一定半径范围内(根据经验而定)已存在特征点的点云;3)在该点云位于遮挡区域的边界上。
步骤1中的“建立同一特征点在两帧相邻时刻激光点云数据中的匹配关系”具体包括特征角点和特征线之间的点云线距离和/或特征平面点和特征面之间的点云面距离通过和得到同一特征点在两帧相邻时刻激光点云数据中的匹配关系,以便计算特征点的相对运动关系。
设定:当前第k帧点云数据表示为Pk,其上一帧(k-1)点云数据表示为Pk-1,对于两类特征点处理分别如下:
对于Pk中的每一个特征角点,下面以第a个特征角点X(k,a)为例,对步骤1中的建立同一特征点在两帧相邻时刻点云数据中的匹配关系的方法进行说明。
在Pk-1中寻找特征角点X(k,a)所对应的特征线段,该线段可以由两点云(b,c)表示,点云b为Pk-1中与特征角点X(k,a)距离最近的点云,点云c是Pk-1中在特征角点X(k,a)相邻扫描线上与特征角点X(k,a)距离最近的点云,相邻扫描线即为雷达扫描线中在径向上与该扫描线最为接近的一条扫描线,其中,表示为下式(2):
式中,X(k,a)表示Pk中的第a个特征角点在雷达坐标系中的坐标,X(k-1,b)表示Pk-1中的第b个特征角点在雷达坐标系中的坐标,X(k-1,c)表示Pk-1中的第c个特征角点在雷达坐标系中的坐标。
对于Pk中的每一个特征平面点,下面以Pk中的第d个特征平面点X(k,d)为例,对步骤1中的建立同一特征点在两帧相邻时刻点云数据中的匹配关系的方法进行说明。
在Pk-1中寻找该特征平面点所对应的特征面,特征面可由三点云(e,f,g)表示,点云e为Pk-1中与特征平面点X(k,d)距离最近的点云,点云f为Pk-1中在点云e同一相邻扫描线上与点云d距离最近的点云,点云g为在相邻扫描线上与特征平面点X(k,d)距离最近的点云,其中,表示为下式(3):
式中,X(k-1,e)表示Pk-1中的第e个特征平面点在雷达坐标系中的坐标,X(k-1,f)表示Pk-1中的第f个特征平面点在雷达坐标系中的坐标,X(k-1,g)表示Pk-1中的第g个特征平面点在雷达坐标系中的坐标。
在一个实施例中,步骤1中,轮式机器人位姿信息可以通过轮式机器人的水平方向平移及偏航方向旋转,来表示相对运动。由于轮式机器人仅有三个自由度,本实施例仅考虑其水平方向平移及偏航方向旋转对应的二维变换矩阵SE(2),其表示为式(4):
式中,qk表示位姿参数,tx表示轮式机器人在惯性坐标系下x方向的平移量,ty表示轮式机器人在惯性坐标系下y方向的平移量,θz表示轮式机器人在惯性坐标系下绕z轴旋转的角度,表示结合上述状态量的二维刚体变换矩阵,即SE(2)。
文中采用上标L表示原三维数据去除z轴维度后的二维数据,变换矩阵等同理。
因此,步骤1中,构建代价函数,即:特征角点和特征线以及特征平面点和特征面之间对应几何关系误差函数,以轮式机器人位姿为变量构造ICP问题,获得轮式机器人位姿信息,该方法具体包括:
步骤b1,根据对应点云的距离关系,如上式(2)和式(3)所示,建立代价函数,将位姿解算问题转化成一个非线性最小二乘的求解问题,如下式(5):
W为一个2×2矩阵,对其进行SVD分解,可得式(8):
W=U∑VT (8)
其中,∑为奇异值对角矩阵,U和V为正交矩阵;
在一个实施例中,步骤2具体包括:
步骤21,根据步骤1输出的轮式机器人位姿信息,提取出可能属于动态障碍物点云。其可以通过如下步骤实现,例如:
步骤212,根据最近邻关系,Pk与Pk-1的交集部分确定为静态点云部分,Pk与Pk-1的非交集部分主要由三部分点云组成:1)属于动态障碍物的点云;2)在上一帧数据中被动态障碍物遮挡的点云;3)新观测到的点云,即首次进入雷达捕获范围物体上的点云。
步骤213,将Pk与Pk-1的非交集部分的点云进行障碍物检测操作,通过观测多帧数据的方式去除属于静态障碍物的目标,从而达到降低算法计算量目的。
步骤22,分割步骤21提取出的可能属于动态障碍物点云中的地面点云与非地面点云,去除地面点云,得到非地面点云。例如,可以将步骤21提取出的可能属于动态障碍物点云放入俯视方向上预设的栅格中,进行点云结构化操作,同时使点云二维化处理,这样可以避免障碍物之间由于地面点云数据造成牵连和噪声干扰。其中,俯视方向可以理解为雷达坐标系Z轴的反方向,即从上到下看。点云结构化操作可以采用比如深度图、KD树或八叉树等现有方法。点云二维化处理可以理解为将原本x、y、z三维信息中的z轴信息归一或者丢掉,由立体三维结构的转化成二维平面点云。
步骤23,将整体的非地面点云按照特定策略聚类为多个点云簇,一个点云簇对应一个障碍物,以区分为各自独立的障碍物。例如:步骤23可以通过如下步骤实现,也可以通过本领域公知的DBSCAN、超体聚类或K-means聚类方法等获得:
步骤24,估计步骤23获得的每个障碍物的位置、形状和尺寸信息。本实施例采用一个二维包围框将每个障碍物框出来,而二维包围框的长宽由其主方向和点云簇中的点云分布共同决定,以描述障碍物的位置、大小和主方向等信息。
例如:步骤24可以通过如下步骤实现,也可以通过本领域公知的最小面积法等方法获得:
步骤242,将移除z轴维度后的障碍物进行随机采样一致性直线分割,将分割出的直线方向作为该障碍物的主方向。
步骤243,利用步骤242确定的主方向建立二维坐标系,比如,设置所述主方向为x′轴,然后设置垂直该x′轴、并且垂直雷达坐标系的z轴y′轴。再根据障碍物所对应点云簇中每个点云的空间信息,计算点云在该坐标系x′、y′两个维度上的极大、极小值。最后将极大值与相应的极小值的差值作为二维包围框的长和宽,即其中,表示二维包围框的长,即障碍物的长度;表示二维包围框的宽,即障碍物的宽度。
在一个实施例中,步骤3具体包括:
步骤31,利用候选动态障碍物信息根据多帧数据中障碍物的差异程度计算差异度函数,并建立关联矩阵,最后通过关联矩阵,可以将多帧点云数据中障碍物关联起来,并成功在多帧数据中关联起来的障碍物即为动态障碍物。
例如:步骤31可以通过如下步骤实现,也可以通过本领域公知的卡尔曼滤波跟踪或基于深度学习等障碍物跟踪方法获得:
式中,表示由式(10)计算的差异度,其中的每一行的第二个下标相同,这意味每一行对应着当前帧的同一个障碍物,那么,第一行的下标均为1,则表示第一障碍物;第二行的下标均为2,则表示第二障碍物;……;第v行的下标均为v,则表示第v障碍物。
步骤313,根据步骤312中建立的关联矩阵(11),关联当前帧Pk和上一帧Pk-1中的动态障碍物。具体地,关联某一障碍物时,将当前帧Pk的关联矩阵的选定行所对应的障碍物作为选定障碍物,再计算选定行中的各元素分别与上一帧Pk-1中关联矩阵中所有的元素的差异度,并以计算得到的最小差异度所对应的Pk-1中元素所对应的障碍物作为选定障碍物的关联障碍物,此时判定选定障碍物与关联障碍物是同一障碍物。
重复上述步骤,直至当前帧Pk的关联矩阵的每一行对应的障碍物关联完毕。
由于以上操作是在潜在动态点云基础上实现的,故关联成功的障碍物一定为动态障碍物。
进一步地,若最小差异度所对应的Pk-1中元素的下标u>v,则表明发生障碍物丢失;若最小差异度所对应的Pk-1中元素的下标u<v,则表明有新动态障碍物出现。
比如:Pk中有10个障碍物,上一帧Pk-1有10个障碍物,并且这10个障碍物在这两帧点云数据中能够一一对应,那么,该关联矩阵为10行10列的矩阵,将第一行的10个元素对应的是当前帧Pk的第一个障碍物作为选定障碍物。而步骤313是将第一行的10个元素与Pk-1的关联矩阵中的所有元素,利用式(10)进行差异度计算,其中,计算得到的最小差异度所对应的Pk-1中元素所对应的障碍物作为选定障碍物的关联障碍物,二者被认定为同一个障碍物。
例如:步骤32可以通过如下步骤实现,也可以通过本领域公知的其它方法获得:
若某障碍物在连续两帧点云数据中被关联,其在Pk-1中的坐标为tk-1=[xk-1,yk-1],Pk中的坐标为tk=[xk,yk],步骤1计算得到的两帧点云数据相对运动为则利用式(12)计算障碍物的运动距离:
经上述步骤可以得到点云数据中动态障碍物信息。
若要得到静态障碍物信息,可建立栅格地图,根据轮式机器人高度将相应范围点云存入栅格地图,若栅格被占据,即表示为轮式机器人不可通过,从而达到在动态环境下高效安全避障目的。
本发明通过障碍物检测及跟踪技术,结合两帧点云数据相对运动,达到检测动态障碍物、提高系统安全性目的。
最后需要指出的是:以上实施例仅用以说明本发明的技术方案,而非对其限制。本领域的普通技术人员应当理解:可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.一种考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,其特征在于,包括:
步骤1,通过激光点云数据的几何特性,提取符合预设曲率特征的点云作为特征点,建立同一特征点在两帧相邻时刻点云数据中的匹配关系,并构建代价函数,以轮式机器人位姿为变量构造ICP问题,获得轮式机器人位姿信息;
步骤2,根据步骤1获得的轮式机器人的位姿信息,检测候选动态障碍物;
步骤3,根据步骤2检测得到的候选动态障碍物,估计动态障碍物状态;
步骤1中,构建代价函数,以轮式机器人位姿为变量构造ICP问题,获得轮式机器人位姿信息的方法具体包括:
W为一个2×2矩阵,对其进行SVD分解,可得式(8):
W=UΣVT (8)
其中,Σ为奇异值对角矩阵,U和V为正交矩阵;
3.如权利要求2所述的考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,其特征在于,所述特征点为特征角点的情形下,步骤1中,所述建立同一特征点在两帧相邻时刻点云数据中的匹配关系的方法具体包括:
在第k-1帧点云数据Pk-1中寻找Pk中的第a个特征角点X(k,a)对应的特征线段,该线段上具有点云b和点云c,其中,点云b为Pk-1中与特征角点X(k,a)距离最近的点云,点云c为Pk-1中在特征角点X(k,a)相邻扫描线上与特征角点X(k,a)距离最近的点云,点云线距离表示为下式(2):
式中,X(k-1,b)表示Pk-1中的第b个特征角点在雷达坐标系中的坐标,X(k-1,c)表示Pk-1中的第c个特征角点在雷达坐标系中的坐标;
所述特征点为特征平面点的情形下,步骤1中,所述建立同一特征点在两帧相邻时刻点云数据中的匹配关系的方法具体包括:
在Pk-1中寻找第d个特征平面点X(k,d)所对应的特征面,该特征面具有点云e、点云f和点云g,点云e为Pk-1中与特征平面点X(k,d)距离最近的点云,点云f为Pk-1中在点云e同一相邻扫描线上与点云d距离最近的点云,点云g为在相邻扫描线上与特征平面点X(k,d)距离最近的点云,点云面距离表示为下式(3):
式中,X(k-1,e)表示Pk-1中的第e个特征平面点在雷达坐标系中的坐标,X(k-1,f)表示Pk-1中的第f个特征平面点在雷达坐标系中的坐标,X(k-1,g)表示Pk-1中的第g个特征平面点在雷达坐标系中的坐标。
4.如权利要求3所述的考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,其特征在于,步骤2具体包括:
步骤21,根据步骤1输出的轮式机器人位姿信息,提取出可能属于动态障碍物点云;
步骤22,分割步骤21提取出的可能属于动态障碍物点云中的地面点云与非地面点云,去除地面点云,得到非地面点云;
步骤23,将整体的非地面点云按照特定策略聚类为多个点云簇,一个点云簇对应一个障碍物;
步骤24,估计步骤23获得的每个障碍物的位置、形状和尺寸信息;
其中,步骤21具体包括:
步骤212,根据最近邻关系,Pk与Pk-1的交集部分确定为静态点云部分;
步骤213,将Pk与Pk-1的非交集部分的点云进行障碍物检测操作,通过观测多帧数据的方式去除属于静态障碍物的目标。
5.如权利要求4所述的考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,其特征在于,步骤23具体包括:
6.如权利要求4所述的考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,其特征在于,步骤24具体包括:
步骤241,将步骤23得到的点云簇的质心定义为该障碍物的坐标;
步骤242,将移除z轴维度后的点云进行随机采样一致性直线分割,将分割出的直线方向作为该障碍物主方向;
步骤243,利用步骤242确定的主方向建立二维坐标系,再根据障碍物所对应点云簇中每个点云的空间信息,计算点云在该坐标系x′、y′两个维度上的极大、极小值,并将极大值与相应的极小值的差值作为障碍物的长和宽。
7.如权利要求4所述的考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,其特征在于,步骤3具体包括:
步骤31,利用候选动态障碍物信息根据多帧数据中障碍物的差异程度计算差异度函数,并建立关联矩阵;
步骤32,根据动态障碍物的关联情况,结合激光里程计计算得到的帧间相对运动,判断障碍物运动状态,从而得到动态障碍物的运动方向及速度。
8.如权利要求7所述的考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法,其特征在于,步骤31具体包括:
步骤312,Pk-1中共检测到的障碍物的总数量为m个,记为
步骤313,根据步骤312中建立的关联矩阵(11),关联当前帧和上一帧点云数据中的动态障碍物,若最小差异度所对应的Pk-1中元素的下标u>v,则表明发生障碍物丢失;若最小差异度所对应的Pk-1中元素的下标u<v,则表明有新动态障碍物出现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110601426.2A CN113345008B (zh) | 2021-05-31 | 2021-05-31 | 考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110601426.2A CN113345008B (zh) | 2021-05-31 | 2021-05-31 | 考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113345008A CN113345008A (zh) | 2021-09-03 |
CN113345008B true CN113345008B (zh) | 2022-05-24 |
Family
ID=77472911
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110601426.2A Active CN113345008B (zh) | 2021-05-31 | 2021-05-31 | 考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113345008B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114442615A (zh) * | 2021-12-31 | 2022-05-06 | 重庆特斯联智慧科技股份有限公司 | 一种基于障碍物属性的机器人行进策略确定方法和系统 |
CN114565616B (zh) * | 2022-03-03 | 2022-12-09 | 湖南大学无锡智能控制研究院 | 一种非结构化道路状态参数估计方法及系统 |
CN114779794B (zh) * | 2022-06-21 | 2022-10-11 | 东风悦享科技有限公司 | 基于无人巡逻车系统的台风场景下的街道障碍识别方法 |
CN115063550B (zh) * | 2022-07-22 | 2022-11-11 | 合肥工业大学 | 一种语义点云地图构建方法、系统及智能机器人 |
CN115993089B (zh) * | 2022-11-10 | 2023-08-15 | 山东大学 | 基于pl-icp的在线四舵轮agv内外参标定方法 |
CN116295507B (zh) * | 2023-05-26 | 2023-08-15 | 南京师范大学 | 一种基于深度学习的激光惯性里程计优化方法、系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105843223B (zh) * | 2016-03-23 | 2018-11-20 | 东南大学 | 一种基于空间词袋模型的移动机器人三维建图与避障方法 |
CN108152831B (zh) * | 2017-12-06 | 2020-02-07 | 中国农业大学 | 一种激光雷达障碍物识别方法及系统 |
GB201803292D0 (en) * | 2018-02-28 | 2018-04-11 | Five Ai Ltd | Efficient computation of collision probabilities for safe motion planning |
CN110221603B (zh) * | 2019-05-13 | 2020-08-14 | 浙江大学 | 一种基于激光雷达多帧点云融合的远距离障碍物检测方法 |
CN111015656A (zh) * | 2019-12-19 | 2020-04-17 | 佛山科学技术学院 | 一种机器人主动避障的控制方法、装置及存储介质 |
-
2021
- 2021-05-31 CN CN202110601426.2A patent/CN113345008B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113345008A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113345008B (zh) | 考虑轮式机器人位姿估计的激光雷达动态障碍物检测方法 | |
CN108647646B (zh) | 基于低线束雷达的低矮障碍物的优化检测方法及装置 | |
CN114384920B (zh) | 一种基于局部栅格地图实时构建的动态避障方法 | |
CN113345018B (zh) | 一种动态场景下的激光单目视觉融合定位建图方法 | |
Weon et al. | Object Recognition based interpolation with 3d lidar and vision for autonomous driving of an intelligent vehicle | |
CN112184736B (zh) | 一种基于欧式聚类的多平面提取方法 | |
WO2022188663A1 (zh) | 一种目标检测方法及装置 | |
CN109255302A (zh) | 目标物识别方法及终端、移动装置控制方法及终端 | |
CN113345009B (zh) | 一种基于激光里程计的无人机动态障碍物检测方法 | |
Jaspers et al. | Multi-modal local terrain maps from vision and lidar | |
CN112346463B (zh) | 一种基于速度采样的无人车路径规划方法 | |
CN113674399A (zh) | 一种移动机器人室内三维点云地图构建方法及系统 | |
Liu et al. | Real-time 6d lidar slam in large scale natural terrains for ugv | |
CN114004869A (zh) | 一种基于3d点云配准的定位方法 | |
CN115861968A (zh) | 一种基于实时点云数据的动态障碍物剔除方法 | |
GB2610410A (en) | Incremental dense 3-D mapping with semantics | |
Ewen et al. | These maps are made for walking: Real-time terrain property estimation for mobile robots | |
CN114066773B (zh) | 一种基于点云特征与蒙特卡洛扩展法的动态物体去除 | |
Bansal et al. | A lidar streaming architecture for mobile robotics with application to 3d structure characterization | |
CN112578673B (zh) | 无人方程式赛车多传感器融合的感知决策与跟踪控制方法 | |
CN116048120B (zh) | 一种未知动态环境下小型四旋翼无人机自主导航系统及方法 | |
Vatavu et al. | Modeling and tracking of dynamic obstacles for logistic plants using omnidirectional stereo vision | |
CN116052099A (zh) | 一种面向非结构化道路的小目标检测方法 | |
CN113554705A (zh) | 一种变化场景下的激光雷达鲁棒定位方法 | |
Wang et al. | A new grid map construction method for autonomous vehicles |
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 |