面向星球探测的岩石三维表面评估技术

  • 格式:pdf
  • 大小:189.48 KB
  • 文档页数:2

第38卷 

、,0I-38 第2期 

No.2 计算机工程 Computer Effgineering 2012年1月 

January 2012 

・工程应用技术与实现・ 文章编号:1o00-__3428(2o12)o2—_o240__一l2 文献标识码:A 中图分类号;TP391 

面向星球探测的岩石三维表面评估技术 

高宏伟,李斌,陈付国 

(沈阳理工大学信息科学与工程学院,沈阳110159) 

摘要:针对星球探测任务中的岩石采样问题,利用双目立体视觉装置采集不同岩石图像,采用立体视觉算法获得岩石表面的三维点云数 

据。以三角剖分后生成的三角片法向量夹角为基本元素,提出基于C均值聚类的岩石表面评估算法,对岩石平整三维表面进行评估计算, 

并用不同颜色加以区分。仿真结果证明了该算法的有效性。 关健诃:星球探测;立体视觉;C均值聚类;岩石评估 

Rock 3D Surface Evaluation Technology for Planetary Exploration 

GAO Hong・wei,LI Bin,CHEN Fu-guo 

(School of Information Science&Engineering,Shenyang Ligong University,Shenyang 1 10159,China) 

[Abstract]According to the rock sample problem in planetary exploration task,binocular vision system is applied to sample different rock images 

and the rocks’surface 3D point cloud data are acquired by stereo vision algorithms in this paper.The relative plane on the rock is automatically evaluated by the proposed method based on C—means clustering between the triangulation pieces’normal vector after triangulation.and the different 

plane is signed by diferent color.Simulation results show the validity of the algorithm. [Key wordsl planetary exploration;stereo vision;C—means clustering;rock evaluation 

D0I:10.3969/j.issn.1000—3428.2012.02.080 

1概述 

我国月球探测系统工程已正式立项并实施,月球探测二 

期工程将发射月球探测器登陆月球,而实施这一目标的重要 

手段之一就是月球漫游车。目前,国内已有一批高校和科研 院所进行了一些月球车方面的预研工作,并取得了一些成 

果…。实现月面漫游是月球车进行其他后续工作的基础,而 

立体视觉系统提供的地形信息则是实现月球车安全漫游的关 键。除此之外,立体视觉系统还可以为车载机械臂提供待考 

察岩石的三维信息。本文以双目立体视觉定位技术为依托, 

针对车载机械臂待考察的岩石为目标,研究岩石重建及表面 

评估技术。 

2双日立体视觉 

2.1双日传感器标定 采用考虑径向畸变的Tsai经典两步方法在高精度的标定 

台上进行标定。在Tsai两步法的第1步中要标定3x3的旋转 

矩阵 、尺度因子S 、t 和t、,是平移矩阵 的前2个元素。 

第2步要标定焦距 厂、f 和径向畸变kt。在七l=0时计算_厂 

和t 的初值,然后根据RAC(Radial Alignment Constraint)模 

型 迭代求取f、t,和k.的最优值。 

2.2立体匹配 本文采用双向灰度相关匹配技术,对于一个候选匹配点 

进行双向确认,这对于剔除奇异点会起到较好的作用 j。具 

体做法是读取2幅图像,以存储左图数据数组的第】个元素 

为基准,在算法设定的搜索空间内在右图中通过灰度相关法 

搜索候选匹配点,记录分数值最大的点。以刚才记录的右图 

中最大分数值的点为基准,反过来在左图中进行相关分数值 

最大的候选匹配点的搜索并记录,判断该点与左图中的基准 

点是否相同,相同则认为匹配正确并保存匹配点对。依次类 推完成剩余的交叉互相关匹配工作。 

2.3三维重建 最为常用的方法是利用标定得到的相机参数根据最小二 

乘法求解空间点三维坐标,然后再根据三角剖分算法和纹理 映射算法在计算机屏幕上形成重建物体的三维造型 l。 

3岩石三维表面评估技术 

3.1 C均值禀类 由于研磨对象由操作者选取指定,一般多选取人眼能分 

辨出来的有平整表面的岩石,且近似平整表面个数可以凭肉 

眼直接判断,这便使得岩石表面自动评估算法的实现成为可 

能。C均值算法只用一个聚类中心点作为一类的代表,一个 

点往往不能充分反映该类的分布结构,从而损失很多有用的 

信息。为此,可以对C均值算法的C类进行调整,采用类核 代替类心,即采用改进的C均值算法,以提高其分类性能 0。 

C均值聚类算法步骤如下: 

(1)设己标准化的待分类特征向量集x.,X ,…,X ),N为 

待分类的特征向量的总数,希望它们分为C类。令特征向量 

Xi=( ,x2 , ) ,定义sum(i)=∑ 。,且令MA=max[sum(i)】, 

MI=minl sum(i)I,计算: 

ai:—(c-1)[sum—(i)-M1]+】,j:1,2,…,Ⅳ (】) MA—MI 显然1≤q≤c,若a 最接近整数 ,则把x 分划至W 

中。对所有样本进行上述处理,就可实现初始分类,从而产 

基金项目:辽宁省教育厅基金资助项目(20080611) 作者筒介:高宏伟(1978一),男,副教授、博士,主研方向:计算机 

视觉,图像处理;李斌、陈付国,硕士研究生 收稿日期:2011-07—11 E-mail:ghw1978@sohu.

com 第38卷第2期 高宏伟,李斌,陈付国:面向星球探测的岩石三维表面评估技术 241 

生初始聚类中心z z ,…,z ,且k=0,k代表迭代次数。 位置之间存在一定的偏离。三角剖分图像如图3所示。 

(2)计算各聚类W,的均值矢量“,和协方差阵 ,。 

1 1 一 一 “,=— ∑ , 7=—I_∑( 一UfU ) (2) J J , (3)将各样本 按最小距离原则分划到某一类中。采用类 

核代替类心的方法来定义一个类核函数K。:K(x,Vi),表示 

类1.,,的分布情况。其中, ,是关于类 ,的一个参数集, 

是一个样本。为了刻画待识别的样本 到 i的接近程度,定 

义距离: 

. 一 d(x,w ) lnll ̄l+( —uj) ’ ( —llJ)一21“ (3) 

(4)计算新类的聚类中心: 

, 、 1 z =— 。 (4) iEl' ̄_ 如果没有样本改变其类别则停止算法,否则转步骤(2)。 

3.2基于c均值聚类的岩石表面评估算法 以各个三角片的法向量为基本元素,采用c均值聚类法 

对各三角片平面进行聚类,聚类算法步骤如下: 

(1)观察岩石图像判断岩石表面有几个近似的平整表面, 

以此作为聚类分析算法的输入,从而确定待分类个数 (大多 

数情况下,一般K=2)。 

(2)求出构成岩石表面的所有三角片的重心和法向量。 (3)任意选取2个三角片的重心和法向量作为初始聚类 

中心。 

(4)对各法向量进行聚类,以向量间的夹角作为优化准 

则。与聚类中心法向量夹角小于5。的聚为一类。 (5)重新计算每类的重心法向量(每类所有法向量与各轴 

夹角的均值),以此作为下一次聚类的中心法向量,如此循环 

聚类直到相邻2次的聚类中心法向量夹角基本不变为止。输 

出聚类中心法向量及相应的重心三维坐标,供机械臂定位 

使用。 

4实验结果与分析 

为了验证本文算法的可行性,分别采集了不同距离与尺 寸的岩石图片,每种岩石上都有2个近似平整的表面。采用 

立体视觉算法对其进行三维重建,如图1、图2所示。 

圈_ 

(a)大岩石 (b)小岩石 

图1原始岩石图像 

(a)大岩石 (b)小岩石 

图2三维点云图像 

在摄像机标定过程中,本文采用欧式空间下的强标定, 

采样点的位置和数量决定了该标定模型应用到实际场景中的 

重建精度。图中重建数据包括了岩石、地面和背景,由于背 景距离摄像机镜头过远,超过了标定时的采样范围,对于这 

一类点的重建是根据原有标定模型推算出来的,因此与实际 嗣黛 

(a)大岩石 (b)小岩石 

图3三角嗣分图像 

由三角剖分及纹理映射效果可以看出,距大岩石很远的 

背景树木精度并不高,但这并不会影响对岩石的重建和评估, 

只要在图像采集时使岩石在标定模型测量范围内即可,虽然 

远处背景的重建效果粗糙,但岩石的主体却十分清晰,并且 

三角剖分和纹理映射效果较好。 

图4是分别对2组岩石数据进行聚类后的结果。 

适 

(a)大岩石侧视图 (b)大岩石聚类结果 

■薅 

(c)小岩石侧视图 (d)小岩石聚类结果 

图4岩石聚类结果 

大岩石被分为2类:岩石正表面和岩石上表面,虽然岩 

石正面的右下部分也是一个过渡圆滑的平整面,可供机械臂 

研磨,但其较小的曲率仅限制于相邻三角片之间,对于该圆 

滑平整面的整体来说,其最大角度范围已经超出三角片法向 量聚类的交角阈值,并且该表面的点云数据在数量上相对其 

他2个平整面较少,因此,在对三角片法向量进行聚类时, 

该部分的点云数据被滤除。 

5结束语 为解决星球探测科学考察任务中漫游车立体视觉系统的 

安全漫游问题,本文研究了岩石三维表面评估技术,提出聚 

类评估算法,并以实验证明了算法的有效性。下一步将实现 

全自动的多个岩石的分割与评估,为漫游车进行科学考察提 

供技术支撑。 

参考文献 

[1】 胡群芳,陈永杰.中国掀起月球车研制热[JJ.中国航天,2004, 

(51:8-9. [2]Tsai R Y A Versatile Camera Calibration Technique for High— accuracy 3D Machine Vision Metrology Using Off-the—shelf TV 

Cameras and Lenses[J].IEEE Journal of Robotics and Automation, 

1987,3(4):323—344. 【3】高宏伟,李斌,梁英,等.立体视觉中误匹配滤波方法的 

研究[J].计算机工程,2008,34(20):210—212. 

[4] Ma Yongzhuang,Liu Weijun.Progressive Matching Based on Segmentation for 3D reconstruction[C]//Proc.of the 5th 

International Conf.on Computer and Information Technology. 

Shanghai,China:IEEE Press,2005. [5】李峰,刘顺利,熊平,等.基于改进C一均值聚类算法的空