发明名称 一种室内环境中移动机器人定位标志点的快速识别方法
摘要 本发明提出了一种室内环境中移动机器人定位标志点的快速图像识别方法,属于计算机视觉技术。该方法通过对摄像机采集图像的二值分割,形态学处理,候选特征点搜索与分类,以及求取候选特征点的权重质心,可以准确的识别天花板灯罩的角点。该方法不受环境光线强度变化的影响,识别精度高、速度快,满足移动机器人实时精确定位的要求。
申请公布号 CN103020632B 申请公布日期 2016.03.30
申请号 CN201210473351.5 申请日期 2012.11.20
申请人 北京航空航天大学 发明人 贾英民;陈霄汉;王龙辉
分类号 G06K9/60(2006.01)I 主分类号 G06K9/60(2006.01)I
代理机构 代理人
主权项 一种室内环境中移动机器人定位标志点的快速识别方法,通过对采集图像的二值图像分割、图像形态学处理、候选特征点的搜索与分类以及候选特征点的质心求取,实现对定位标志点的快速可靠识别,其特征在于该方法的步骤如下:步骤1,将摄像机安装在移动机器人上,调整光轴方位,使摄像机的成像平面平行于天花板;步骤2,颜色分割摄像机抓取的图像,将分割的结果保存为一个二值图像I<sub>A</sub>,其中高于阈值的像素点被认为是前景像素点,像素值为1,低于阈值的像素点被认为是背景像素点,像素值为0;步骤3,对二值图像进行形态学处理,先使用膨胀算子然后再使用腐蚀算子,从而去除图像前景像素中的“洞”和“凹”的部分,处理结果仍然保存在二值图像I<sub>A</sub>中;步骤4,创建方阵M<sub>(2r+1)×(2r+1)</sub>,它的元素满足<img file="FDA0000816919260000011.GIF" wi="829" he="166" />步骤5,遍历I<sub>A</sub>中的每一个像素点,判断该点是否为边界点,如果是边界点,标记此点为停留点,以该停留点为中心点从I<sub>A</sub>中提取维数为(2r+1)×(2r+1)的子图像<img file="FDA0000816919260000015.GIF" wi="70" he="65" />步骤6,对子图像<img file="FDA0000816919260000016.GIF" wi="53" he="65" />实施如下运算<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mfrac><mrow><mi>&Sigma;</mi><mrow><mo>(</mo><mi>M</mi><mo>&CircleTimes;</mo><msubsup><mi>I</mi><mi>A</mi><mi>S</mi></msubsup><mo>)</mo></mrow></mrow><mrow><mi>&Sigma;</mi><mi>M</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000816919260000012.GIF" wi="350" he="164" /></maths>M为步骤4中创建的方阵,算子<img file="FDA0000816919260000013.GIF" wi="47" he="54" />表示矩阵按元素相乘,算子Σ(·)表示对矩阵元素求和,如果S∈[S<sub>1</sub>,S<sub>2</sub>],其中S<sub>1</sub><1/4<S<sub>2</sub>,对应的停留点被认为是一个候选特征点,记录该停留点的像素坐标;步骤7,对所有候选特征点分类,被半径为d<sub>r</sub>的圆包围的候选特征点被认为属于同一类;步骤8,求解每一类候选特征点的权重质心,该权重质心对应的像素坐标即为识别到的标志点坐标,每一类候选特征点的权重质心(u<sub>c</sub>,v<sub>c</sub>)的求取公式为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>u</mi><mi>c</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>u</mi><mi>k</mi></msub><mo>/</mo><msub><mi>S</mi><mi>k</mi></msub></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mn>1</mn><mo>/</mo><msub><mi>S</mi><mi>k</mi></msub></mrow></mfrac><mo>,</mo><msub><mi>v</mi><mi>c</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>v</mi><mi>k</mi></msub><mo>/</mo><msub><mi>S</mi><mi>k</mi></msub></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mn>1</mn><mo>/</mo><msub><mi>S</mi><mi>k</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000816919260000014.GIF" wi="686" he="175" /></maths>其中n为每一类候选特征点包含的候选特征点的数目,(u<sub>k</sub>,v<sub>k</sub>)为第k个候选特征点的像素坐标,1/S<sub>k</sub>为第k个候选特征点的权重,S<sub>k</sub>由步骤6计算得到。
地址 100191 北京市海淀区学院路37号