发明名称 一种室内环境下考虑地面约束的RGB‑D视觉里程计方法
摘要 本发明属于移动机器人自主定位领域,特别涉及一种室内环境下考虑地面约束的RGB‑D视觉里程计方法。本发明包括:(1)室内场景中,利用搭载在移动平台的RGB‑D相机在移动过程中采集环境中的彩色图像信息和深度图像信息;相邻帧图像分别记为I<sub>s</sub>和I<sub>t</sub>;(2)根据深度图像信息通过逆投影函数π<sup>‑1</sup>求解获得环境的三维点云数据V<sub>s</sub>和V<sub>t</sub>;(3)ORB特征提取与匹配:采用ORB算法实现RGB‑D图像特征点的提取与匹配。本发明采用ORB算法完成特征提取与匹配的图像预处理过程,较SIFT和SURF算法,快速性提高一个数量级;利用深度图像检测获取点云地面并利用地面信息增强点集对齐,提高运动变换的估计精度。
申请公布号 CN106556412A 申请公布日期 2017.04.05
申请号 CN201610940748.9 申请日期 2016.11.01
申请人 哈尔滨工程大学 发明人 赵玉新;李亚宾;刘厂
分类号 G01C22/00(2006.01)I;G06T7/20(2017.01)I;G06T7/11(2017.01)I;G06T7/136(2017.01)I;G06K9/46(2006.01)I 主分类号 G01C22/00(2006.01)I
代理机构 代理人
主权项 一种室内环境下考虑地面约束的RGB‑D视觉里程计方法,其特征在于,包括以下几个步骤:(1)室内场景中,利用搭载在移动平台的RGB‑D相机在移动过程中采集环境中的彩色图像信息和深度图像信息;相邻帧图像分别记为I<sub>s</sub>和I<sub>t</sub>;(2)根据深度图像信息通过逆投影函数π<sup>‑1</sup>求解获得环境的三维点云数据V<sub>s</sub>和V<sub>t</sub>:<maths num="0001"><math><![CDATA[<mrow><mi>p</mi><mo>=</mo><msup><mi>&pi;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>Z</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><mi>u</mi><mo>-</mo><msub><mi>c</mi><mi>x</mi></msub></mrow><msub><mi>f</mi><mi>x</mi></msub></mfrac><mi>Z</mi><mo>,</mo><mfrac><mrow><mi>v</mi><mo>-</mo><msub><mi>c</mi><mi>y</mi></msub></mrow><msub><mi>f</mi><mi>y</mi></msub></mfrac><mi>Z</mi><mo>,</mo><mi>Z</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi></msup></mrow>]]></math><img file="FDA0001138256620000011.GIF" wi="807" he="135" /></maths>其中,p=(X,Y,Z,1)<sup>T</sup>表示环境三维点在相机坐标系下的齐次坐标;m=(u,v)<sup>T</sup>表示环境三维点投影到图像像素坐标系下的坐标;Z=Z(m)表示该点的深度数据;RGB‑D相机横纵方向上以像素为单位的焦距分别为f<sub>x</sub>,f<sub>y</sub>,光心为(c<sub>x</sub>,c<sub>y</sub>);(3)ORB特征提取与匹配:采用ORB算法实现RGB‑D图像特征点的提取与匹配;(4)以RGB‑D图像源帧I<sub>s</sub>对应的摄像头位姿的高度、俯仰角信息为先验知识并基于图像源帧到目标帧的相机小位移假设,对源、目标帧RGB‑D图像地面最可能区域进行局部搜索,从而快速而准确地对地面区域进行检测;考虑到点云数据异常点的存在,采用引入异常点剔除机制的地面平面检测算法;(5)添加地面约束:利用地面信息,通过施加扭动T<sub>inc</sub>将摄像头变换至新的位姿:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msup><mi>X</mi><mo>&prime;</mo></msup><mo>=</mo><msub><mi>T</mi><mrow><mi>i</mi><mi>n</mi><mi>c</mi></mrow></msub><mo>&CenterDot;</mo><mi>X</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>i</mi><mi>n</mi><mi>c</mi></mrow></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>R</mi><mrow><mi>i</mi><mi>n</mi><mi>c</mi></mrow></msub></mtd><mtd><msub><mi>t</mi><mrow><mi>i</mi><mi>n</mi><mi>c</mi></mrow></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0001138256620000012.GIF" wi="365" he="229" /></maths>其中,平移增量t<sub>inc</sub>可利用下式求解得到:t<sub>inc</sub>=(t<sub>x</sub>,t<sub>y</sub>,t<sub>z</sub>)=(x′‑x,y′‑y,z′‑z)式中,(x′,y′,z′)和(x,y,z)分别是平面重心m′和m的空间坐标;采用欧拉角表达旋转增量<img file="FDA0001138256620000013.GIF" wi="328" he="63" /><img file="FDA0001138256620000014.GIF" wi="51" he="46" />β,ψ表示<img file="FDA0001138256620000015.GIF" wi="38" he="47" />经旋转变换至<img file="FDA0001138256620000016.GIF" wi="40" he="47" />绕x,y,z坐标轴分别转过的角度,将欧拉角式转换为四元数,转换关系如下:<img file="FDA0001138256620000017.GIF" wi="942" he="527" />式中,字符c,s分别为cos和sin的缩写形式;(6)根据地面检测是否成功,即所检测获取的地面区域在图像区域占比是否超过20%,设置标志位flag:超过20%,地面拟合参数可靠准确,检测成功,flag=true,并在地面重合约束下对三维对应点集Fs和Ft施加扭动T<sub>inc</sub>,得到F<sub>s</sub>′和F<sub>t</sub>′;反之,,视作检测失败,flag=false;(7)依据标志位不同对地面检测失败和地面检测成功的对应点集利用随机采样一致性算法求解相机运动变换,并对成功进行地面校正的变换矩阵经恢复得到原始变换,从而得到了使原图像F<sub>s</sub>、F<sub>t</sub>三维匹配点集对齐的刚体变换估计T<sub>1</sub>、T<sub>2</sub>;(8)计算T<sub>1</sub>、T<sub>2</sub>二者对应的配准变换误差e<sub>i</sub>,误差较小者即为最佳的里程计运动变换T<sup>*</sup>。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室