适应于环境空间变化的激光雷达SLAM建图方法.pdf
《适应于环境空间变化的激光雷达SLAM建图方法.pdf》由会员分享,可在线阅读,更多相关《适应于环境空间变化的激光雷达SLAM建图方法.pdf(11页珍藏版)》请在咨信网上搜索。
1、系统仿真学报系统仿真学报Journal of System Simulation第 35 卷第 8 期2023 年 8 月Vol.35 No.8Aug.2023适应于环境空间变化的激光雷达适应于环境空间变化的激光雷达SLAM建图方法建图方法焦嵩鸣,姚鑫*,丁辉,钟宇飞(华北电力大学 自动化系,河北 保定 071003)摘要摘要:为解决在空间大小发生明显变化的环境中,现有算法建图容易发生漂移等问题,提出了一种适应于环境空间变化的Adp-lio-sam建图方法,以提高lio-sam算法的通用性。改进了点云去畸变方法,使用KF融合激光雷达帧间位姿插值和IMU插值后得到的位姿进行运动补偿;使用模糊算法
2、对不同大小的空间环境适配不同的点云滤波阈值;优化了回环检测约束条件。实验结果表明:在环境空间大小改变、特征类型改变、结构相似度改变的实验中,改进算法比已有算法的回环误差平均可降低47.6%,z轴平均误差平均可降低36.4%,地图质量明显提升,且满足实时性需求。关键词关键词:机器视觉;空间变化;运动补偿;自适应阈值;地图构建中图分类号:TP391.9;TP212 文献标志码:A 文章编号:1004-731X(2023)08-1788-11DOI:10.16182/j.issn1004731x.joss.22-0489引用格式引用格式:焦嵩鸣,姚鑫,丁辉,等.适应于环境空间变化的激光雷达SLAM建
3、图方法J.系统仿真学报,2023,35(8):1788-1798.Reference format:Jiao Songming,Yao Xin,Ding Hui,et al.Lidar SLAM Mapping Method Adapted to Environmental Spatial ChangesJ.Journal of System Simulation,2023,35(8):1788-1798.Lidar SLAM Mapping Method Adapted to Environmental Spatial ChangesJiao Songming,Yao Xin*,Ding H
4、ui,Zhong Yufei(Department of Automation,North China Electric Power University,Baoding 071003,China)Abstract:In the environment with obvious changes in space size,aiming at the drift and other problems of the existing algorithm,Adp-lio-sam mapping method is proposed to adapt to the environment space
5、changes,and improve the generality of lio-sam algorithm.Point cloud dewarping method is improved,and Kalman filter algorithm is used to carry out the motion compensation data by fusing lidar inter-frame pose interpolation and IMU interpolation.Fuzzy algorithm is used to adapt different points filter
6、ing thresholds for different spatial environments and the constraints of loop closure detection are optimized.Experimental results show that,compared with the existing algorithm in the experiments of changing the size of the environment space,the type of features,and the similarity of the structure,
7、the improved algorithm can reduce the loopback error by 47.6%and the z-axis average error by 36.4%with high map quality and good real-time performance.Keywords:machine vision;spatial change;motion compensation;adaptive threshold;map building收稿日期:2022-05-12 修回日期:2022-07-23第一作者:焦嵩鸣(1972-),男,副教授,博士,研究方
8、向为机器人智能感知及电力机器人自主巡检。E-mail:通讯作者:姚鑫(1998-),男,硕士生,研究方向为多传感器融合SLAM。E-mail:第 35 卷第 8 期2023 年 8 月Vol.35 No.8Aug.2023焦嵩鸣,等:适应于环境空间变化的激光雷达SLAM建图方法http:/www.china-0 引言引言SLAM技术所受到的关注越来越多1,地图构建作为SLAM技术中的重要一环,一直是研究热点之一。常见的地图构建方法主要使用的传感器有相机和激光雷达2种,而激光雷达相较于相机具有精度高、抗干扰能力强、实时性好等特点,因此,广泛应用在建图问题中2-3。3D激光雷达能够提供丰富的环境信
9、息,实现较高精度的定位与建图。文献4提出一种基于模糊逻辑的SLAM方法,通过建立模糊规则,提高系统数据关联性,使系统在建图过程中具有更好的鲁棒性;LOAM(lidar odometry and mapping in real-time)算法5利用线面结构有效提取环境信息,并使用匀速模型对点云运动补偿,但建图时发生较为快速的旋转问题时,地图极易发生漂移;文献6添加了闭环检测功能,但建图依赖地面特征,适合地面较为平坦的环境;LIO_mapping7(tightly coupled 3d lidar inertial odometry and mapping)是一种激光雷达和IMU的紧耦合方案,提升
10、了系统在剧烈运动、长走廊等环境下建图时的鲁棒性,但并没有在变化场景中进行验证;文献8采用了多阶段匀加速模型对点云畸变进行运动补偿,取得了较好的效果,但仍不能对真实环境充分拟合;文献9提出一种自适应体素映射的方法,使得系统在场景退化的情况下也可以保持建图效果;Fast-LOAM10(fast lidar odometry and mapping)是一种采用非迭代的两级失真补偿算法,该方法有效降低了建图过程中的计算量,提升了建图效率;文献11在图优化框架上提出 lio-sam 算法,利用 IMU 对点云去畸变,比其他激光雷达与 IMU 紧耦合建图方案更优;文献12结合语义信息促进环境映射过程,与纯
11、几何方法相比,在动态环境中更具有优势;文献13使用基于关键帧和普通帧的局部地图匹配方法,弥补了地图匹配过程中数据关联错误的不足,并且提高了系统在灌木丛、动态对象场景中的建图准确性,但增加了计算量。综上所述,大部分算法对建图环境本身关注度较少,对出现明显空间大小变化的场景适应性较差。借鉴文献14-15中视觉SLAM采用自适应阈值提取特征点来应对环境变化的思想,以lio-sam算法为基础,提出了一种适应于环境空间变化的激光雷达与 IMU 融合的 Adp-lio-sam 建图方法,重点对空间明显变化的环境进行建图研究。1算法流程算法流程Adp-lio-sam算法流程如图1所示,整体分为4个部分。(1
12、)数据预处理:首先,对IMU数据进行预积分处理,为后续点云去畸变和激光雷达里程计提供位姿变换依据;其次,对嘈杂的激光点云进行剔除。图1 Adp-lio-sam算法流程Fig.1 Adp-lio-sam algorithm flowchart 1789第 35 卷第 8 期2023 年 8 月Vol.35 No.8Aug.2023系统仿真学报Journal of System Simulationhttp:/www.china-(2)KF点云插值去畸变:使用IMU数据和处理后的点云数据对当前激光点进行帧间插值,并将两者KF融合后的插值结果用于点云去畸变,之后对去畸变的点云进行特征分类。(3)自适
13、应阈值选取策略:采用模糊算法对不同大小的空间环境实时适配不同大小的滤波阈值,从而保证算法在建图过程中点云密度分布均匀,之后结合IMU数据进行局部地图构建。(4)回环检测与后端优化:同时使用激光里程计和IMU状态预测结果对系统位姿进行约束,并采用优化后的回环条件进行检测,最后输出期望的地图与轨迹。2改进点云去畸变算法改进点云去畸变算法激光扫描伴随着机器人运动的整个周期,因此每帧点云中的每个点都会保存着不同的数据信息。当这些点以同一时间戳向外发布时,就会产生运动畸变,从而对建图质量产生严重影响,所以,需要对点云进行运动补偿。目前,主流方法有纯估计法和传感器辅助法,然而纯估计法匀速模型或分段线性插值
14、模型的假设较为苛刻,IMU传感器通过加速度2次积分得到的位移累计误差较大。为达到更好的运动补偿效果,提出一种KF融合方法来改善运动畸变影响。2.1 激光雷达帧间位姿插值激光雷达帧间位姿插值帧间位姿插值指在匀速模型假设前提下,利用角度偏差将当前帧所有激光点补偿到当前帧起始时刻,插值构思如图2所示。记上一帧到当前帧的帧间位姿为Tse,当前点相对于当前帧起始时刻的水平角度为c,则当前点相对于当前帧起始时刻的相对位姿Tsc(可分解为旋转矩阵Rsc和平移矩阵tsc)可由式(1)计算得到。假设激光雷达符合匀速模型,匀速模型指在点云去畸变过程中,k帧到k+1帧之间的位姿变换与k-1帧到k帧之间求得的位姿变换
15、保持相同。Tsc=cTse(1)2.2 IMU插值插值IMU插值指利用IMU的高频位姿对低频的激光雷达进行插值,从而完成运动补偿,插值构思如图3所示,其中,tp为当前激光点对应的时间;tb为tp之前相邻的IMU时间;tf为tp之后相邻的IMU时间;*为IMU和Lidar各自对应的数据;t为时间轴;虚线为数据对应的时间。通过插值得到当前激光点的位姿:Tp=tp-tbtf-tbTf+tf-tptf-tbTb(2)式中:Tb为Tp之前相邻时刻tb的IMU位姿;Tf为Tp之后相邻时刻tf的IMU位姿。通过同样的方法可以求得当前帧起始时刻激光点对应的位姿Ts,则当前点相对于当前帧起始时刻的位姿Tsp(可
16、分解为旋转矩阵Rsp和平移矩阵tsp)可由式(3)得到:Tsp=TsT-1p(3)2.3 KF插值数据融合插值数据融合KF可以充分考虑传感器测量误差的影响,从而得到最优估计结果。滤波器输入为帧间位姿插图2 帧间位姿插值Fig.2 Inter-frame pose interpolation图3 IMU位姿插值Fig.3 IMU pose interpolation 1790第 35 卷第 8 期2023 年 8 月Vol.35 No.8Aug.2023焦嵩鸣,等:适应于环境空间变化的激光雷达SLAM建图方法http:/www.china-值和 IMU 插值得到的位姿,输出为最优估计位姿:Rcp
17、=Rsp+e1e1+e3(Rsc-Rsp)tcp=tsp+e2e2+e4(tsc-tsp)(4)式中:最优位姿Tcp可分解为旋转矩阵Rcp和平移矩阵tcp;e1、e2、e3、e4分别为IMU和激光雷达对旋转和平移的置信度。通过最优位姿Tcp把当前点P补偿到当前帧起始时刻,得到运动补偿后当前点的坐标:Pnew=RcpP+tcp(5)3自适应阈值选取方法自适应阈值选取方法滤波阈值u表示每个边长为u的立方体中只保存一个激光点。阈值太大,会使有效点云较少甚至建图匹配失败;阈值太小,虽然能保留充分的点云,但会在地图匹配上消耗较多时间。而实际建图过程中往往会遇到环境空间大小变化明显的情况,如从空旷的操场进
18、入教学楼、从宽阔的马路进入狭窄的巷道,此时使用固定阈值滤波就会对建图产生一定的影响。针对环境空间大小的改变,提出一种激光雷达自适应阈值滤波方法,采用模糊算法对不同的空间大小适配不同的滤波阈值。保存当前帧点云中每个激光点的距离,并将当前帧激光点距离在1100 m范围内的点分到5个不同的距离区间内,对5个距离区间从小到大编号15,计算每个距离区间内点的数量,并选取点数最多的距离区间对应的编号作为算法输入(记为n);选取滤波阈值作为算法输出(记为u),范围为0.2 0.6,m。对算法的输入、输出进行模糊处理;取编号n的模糊集为负大(NBn)、负小(NSn)、零(ZOn)、正小(PSn)、正大(PBn
19、),将编号n的变化范围分为5个等级:-2、-1、0、1、2,通过三角形隶属度函数计算编号n对应模糊集的隶属度;取滤波阈值u的模糊集为负大(NBu)、负小(NSu)、零(ZOu)、正小(PSu)、正大(PBu),将滤波阈值u的变化范围分为5个等级:-2、-1、0、1、2,通过三角形隶属度函数计算滤波阈值u对应模糊集的隶属度。隶属度函数分布如图4所示。通过多次实验发现,滤波阈值的大小与空间环境的大小成近似正比关系,即空间越小所适配的滤波阈值越小,空间越大所适配的滤波阈值越大。根据上述经验,设计如下模糊规则:若n负大则u负大,若n负小则u负小,若n为零则u为零,若n正小则u正小,若n正大则u正大。确
20、定模糊关系:R=(NBnNBu)(NSnNSu)(ZOnZOu)(PSnPSu)(PBnPBu)(6)确定观测量后,采用最大隶属度法计算模糊值:U=max(ui)(7)计算滤波阈值:u=a+b2+b-a2nU(8)式中:a=0.2、b=0.6。获取匹配当前环境的滤波阈值后,激光里程计通过使用降采样后的点云,匹配得到激光雷达帧间的旋转和平移变化量。首先,用当前帧和上一帧的平面点构建点到面的距离约束,并计算梯度下降方向的残差。其次,用当前帧和上一帧的图4 隶属度函数分布Fig.4 Membership function distribution 1791第 35 卷第 8 期2023 年 8 月V
21、ol.35 No.8Aug.2023系统仿真学报Journal of System Simulationhttp:/www.china-边缘点构建点到线的距离约束,并计算梯度下降方向的残差。最后,使用高斯牛顿算法对上述条件联合优化,从而得到帧间匹配的位姿。IMU输出结果为线加速度和角加速度,以距离激光雷达当前帧起始时刻之前最近的IMU数据所对应的时间为起点,距离激光雷达下一帧起始时刻之后最近的IMU数据所对应的时间为终点,求取起点到终点的时间变化量,通过IMU预积分得到IMU在这段时间内的旋转和平移变化量16。将IMU和激光里程计得到的帧间变换矩阵送入 GTSAM(georgia tech s
22、moothing and mapping),进而对当前状态进行预测并将此状态更新到地图中。4回环检测与后端优化回环检测与后端优化回环检测模块17通过计算当前帧与历史关键帧之间的欧式距离,作为判断系统是否发生回环的先验,在满足距离约束的关键帧中选取一个满足时间约束和序列号约束的关键帧,作为待回环帧,判断标准为 dcur-dpur20cur-pur100(9)式中:cur为当前帧序列号;pur为待回环帧序列号。使用 ICP(iterative closest point)算法构建待回环帧与当前帧位姿的误差函数:f(R t)=1Ni=1NQi-(RPi+t)2(10)式中:R为旋转矩阵;t为平移矩阵
23、;Pi为当前帧点云中一点;Qi为待回环帧中与Pi对应的点;N为点云点数。回环若没有触发,则等待准备下一次检测;若触发,则需要对待回环帧与当前帧之间所有的关键帧进行位姿优化,即固定待回环帧,将回环检测求得的变换矩阵作为参数输入优化器中进行优化求解。最后用iSAM优化器18-19求解出优化后的带回环检测位姿,从而构建全局一致地图。5实验与分析实验与分析实验主要使用的器件及软件平台:北科天绘R-Fans-16导航雷达,垂直方向共16线,视角为30(-15+15),线间距为2,水平方向单线采集激光点2 000个,视角为360,线间距约为0.090.36,采集频率520 Hz;IMU为xsens九轴惯性
24、传感器,型号为MTi 1-s,采集频率最高100 Hz;搭 载 IntelCore-i5-7300HQ 处 理 器 的 计 算 机;Ubuntu18.04下的机器人操作系统;开源SLAM算法评价软件EVO,通过自建数据集,分别验证点云运动补偿效果、自适应阈值效果、场景变化适应实验效果。5.1 运动补偿效果对比运动补偿效果对比为验证KF融合后的点云去畸变效果,实验采用控制变量法,即保持算法其他部分不变,只改变运动补偿部分,实验在一个较为方正的环境中进行,如图5所示。图中用红色直线方框标注墙壁,并为判断墙壁是否笔直提供参考,由实验结果可知,通过KF融合后的改进运动补偿算法效果最优。图5 运动补偿结
25、果比较Fig.5 Comparison of motion compensation results 1792第 35 卷第 8 期2023 年 8 月Vol.35 No.8Aug.2023焦嵩鸣,等:适应于环境空间变化的激光雷达SLAM建图方法http:/www.china-图5(a)为一个较为方正的楼梯口,通过墙面的笔直程度来判断点云运动补偿的效果。图 5(b)为没有运动补偿的效果,畸变严重,墙壁出现倾斜、断层现象。图5(c)为激光雷达帧间插值补偿效果,匀速运动时,下面红色方框内墙壁成左高右低式倾斜,这是因为激光雷达会出现误匹配的情况。图5(d)为IMU插值补偿效果,上面红色方框内墙壁左低
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 适应 环境 空间 变化 激光雷达 SLAM 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。