CN103759671A - 一种牙模三维表面数据非接触扫描方法 - Google Patents

Abstract

本发明公开了一种牙模三维表面数据的非接触扫描方法，采用激光发射器、高精度平移台、和两台网络摄像机构成的三维扫描系统，每个摄像机独立完成三维扫描，然后以其中一台摄像机为坐标原点，通过标定点的三维坐标信息，完成摄像机位置的刚体变换过程求解，最终将扫描的数据经过这个刚体变换，得到在同一个坐标系下的牙模三维表面数据信息，解决了单目场景下的激光平面的标定、线条实时高效的提取、三维实时扫描等相关的技术难题。最终的扫描误差低于0.05mm。在分辨率为640*480下，能够实现20秒小尺寸物体的快速扫描。最终左摄像机的扫描数据经过刚体变换与右摄像机的扫描数据合成，可以得到比较完整的牙模三维扫描数据。

Description

一种牙模三维表面数据非接触扫描方法

技术领域

本发明属于基于线结构激光的三维扫描技术领域，涉及双目计算机视觉位姿标定提前配准的非接触三维扫描，具体涉及一种三维表面数据非接触扫描方法。

背景技术

物体几何尺寸的精密测量在光电检测技术中是一个应用非常广泛且实用价值非常高的课题。计算机视觉测量方法具有无接触、准确度高、数据便于处理、易于自动控制等优点，因此基于计算机视觉的光电测量方法作为未来生产生活不可缺少的技术被许多研究者研究。随着计算机图像处理和模式识别技术的迅速发展，在精密测量中采用计算机视觉方法相对于纯机械方法或者其他光学方法具有操作直观简便、适应范围广、精度高等特点，具有很好的应用前景，已成为国内外研究的热点。

物体的三维轮廓测量法的用途很广，可以应用在机器人导航、视觉监控、监控制造等行业。同时在最近几年兴起的计算机虚拟现实，基于图像的绘制，三维动画等技术中都有广泛应用。物体的三维轮廓重构的实现方法很多，主要包括接触式和非接触式两大类。在非接触中，光学方法应用最为广泛。物体三维重构的光学方法又可以分为主动视觉和被动视觉两类。前者由物体辐射信号或者物体表面反射信号生成图像，成像设备不发出信号；后者是由成像设备发出一束信号，通过接收从物体表面反射回来的或穿透物体的信号生成图像。基于主动视觉的非接触的结构光法以其固有的非接触、易于实现和较高的精度等优点，这种方法的特点是测量过程中侧头不接触被测表面，且测量速度快，适用于各种软、硬材料的各种复杂曲面模型的三维测量。在考古测量与虚拟博物馆，虚拟现实、机械加工等逆向工程越来越多的得以应用。尤其在进入九十年代以来，随着工业检测技术、逆向工程技术和快速成型技术的迅速发展，对三维物体进行表面轮廓重构和测量的需求越来越大，对测量速度和精度要求越来越高，使得三角法式的结构光法成为使用最为广泛的方法。

发明内容

为获取牙模三维表面数据，本发明的目的在于，提供一种牙模三维表面数据的非接触扫描方法，该方法利用一套摄像机位置刚体变换求解的立体式三维扫描系统，能够在20秒内完成牙模数据扫描过程，且误差在0.05毫米范围内。

为了实现上述任务，本发明采取如下的技术解决方案：

一种牙模三维表面数据的非接触扫描方法，其特征在于，该方法采用激光发射器、高精度平移台、和两台网络摄像机构成的三维扫描系统，每个摄像机独立完成三维扫描，然后以其中一台摄像机为坐标原点，通过标定点的三维坐标信息，完成摄像机位置的刚体变换过程求解，最终将扫描的数据经过这个刚体变换，得到在同一个坐标系下的牙模三维表面数据信息，具体按下列步骤操作：

1）制作标定板，标定板为三列水平方向等距的三个黑圆圈，标定板放置在高精度平移台的台面上，标定板内的标定点分别按从左到右，从上到下的顺序排列，获取所有标定点在图像中的位置坐标；

2）两个摄像机对称设置在高精度平移台两侧，分别作为独立的单目扫描仪工作，两个摄像机与高精度平移台的台座平面有一定偏角，便于求解摄像机的刚体变换；

3）根据三维空间中非共线的三点确定唯一平面，即可得到台面的平面方程；

4）借用激光平面与物体表面相交求解物体表面三维空间位置坐标值，利用平移滑台的匀速运动的特性，添加物体平移值，得到物体的整个表面的三维点云坐标；

5）在扫描过程中，将激光条中心一点在图像中的坐标点转换为成像平面上一点，该点与坐标原点构成的直线与激光平面相交于一点，这点为激光线打在物体上面的物理空间中的点；

6）高精度平移台左侧摄像机经过选择和平移变换，得到与右摄像机的扫描数据在同一坐标系下的完整的三维扫描数据。

所述的坐标点为9个，其中每三个坐标点均匀分布在一条直线上，形成三条相互平行的直线。

所述的三维扫描数据在处理阶段，将扫描得到的数据经过刚体变换合并左右摄像机的扫描数据。

在图像中激光光条上中心线上，在台座平面内找一个点，在平移台面内找两点；这三点的三维坐标可通过相应公式计算得到，由这三个空间点即可确定激光平面方程。

所述的标定点的获取方法包括以下步骤：

第一步，读入图像，对图形进行灰度化、高斯滤波；

第二步，运用Canny算子提取轮廓；

第三步，提取轮廓信息，用椭圆逼近算法，提取轮廓的最小逼近椭圆。

第四步，根据OpenCV提供的函数接口，可以取得轮廓的周长和面积，由此可以得出一个约束条件：

$\mathrm{Thred}=\frac{C^2}{\mathrm{Area}}=4*\mathrm{\π};$

第五步，取得逼近椭圆的中心点，判断该点出的像素值大小，给定该阈值<127，判定该标定点为黑色中心点。

所述的激光条中心采用阈值重心法，包括以下步骤：

首先在图像的空间域中首先给定一个阈值，当图像中某个位置的亮度值大于设定的阈值时提取该点，否则跳过该点，进行下一个；

然后再将提取到的点运用重心法进行计算就得到光条中心结构光中心点；

最后按照图像的高斯分布法则，在上述提取的中心点上，对每个点进行Hessian矩阵运算，得到中心点的法线，最终得到光条中心的亚像素精确点。

与现有技术相比，本发明的牙模三维表面数据非接触扫描方法，带来的有益效果如下：

1、基于单目的经过几何投影变换求解得到的滑动平移台面的平面方程的精度最大相对误差小于0.2%，系统能够快速实现场景的高精度建模。

2、巧妙的借用了扫描平台与台座这两个平面，实现了一种快速的激光平面标定方法，最终的扫描绝对误差低于0.05mm。

在分辨率为640*480下，能够实现20秒小尺寸物体的快速扫描。最终左摄像机的扫描数据经过刚体变换与右摄像机的扫描数据合成，可以得到比较完整的三维扫描数据。

本发明提供的牙模三维表面数据的非接触扫描方法，复杂环境下结构光中心提取效果好，结构光中心的有效提取也直接影响三维扫描的准确度，目前所用到的光心提取算法为重心法，大多数对环境光要求比较苛刻，本发明可以利用光区域阈值法重心法快速提取光条中心，对环境光的要求比较低，提取效果良好。

本发明提供的牙模三维表面数据的非接触扫描方法，能够对被测量对象的三维数据进行有效存储：目前的存储方案中三维数据的抖动性较大，实验采样的重复性不好；而本发明可以得到测量对象的精确三维点数据，并按照obj格式存储，数据量较小，便于结果的显示和分析；对被扫描对象进行数字化，可通过相应的图形软件对被测对象进行更加细致的分析与研究。

本发明提供的牙模三维表面数据的非接触扫描方法，构建设备成本非常低廉：用最简单的仪器设备来构建，成本低廉，不超过8千元，在进行图像摄取时仅需要二个高质量网络摄像头（800元左右*2），一个一字线激光笔（1000元左右），一个电动平移台（5400元左右），相对于市场上现有的各种三维扫描系统动辄10余万、几十万来说成本非常低廉。

本发明提供的方法操作技术要求低，对于非接触方式的三维扫描方式一般要求较高的使用技能，方法简单易操作，不需要具备较高的专业知识。且扫描速度快，非接触测量一般要进行数据分析，所以测量时间一般比较慢；本发明可以实现实时扫描，能够即时显示测量结果，并且操作极其便利。

附图说明

图1为本发明的牙模三维表面数据非接触三维扫描系统的硬件构造图；

图2为本发明的牙模三维表面数据非接触三维扫描系统的工作流程图；

图3为本发明的牙模三维表面数据非接触三维扫描系统功能模块图；

图4为标定纸的样式；

图5为激光平面标定图片；

图6为标定点的提取效果图片；

图7为结构光光条中心（激光线）提取效果图片；

图8为扫描获得的牙模三维点云数据图片。

以下结合附图和实施例对本发明作进一步详细说明。

具体实施方式

本发明的牙模三维表面数据非接触三维扫描方法，以单线式结构光测量为基本模型，展开三维扫描系统的研究，主要研究基于单目机器视觉为基本原理的三维扫描系统，进一步获取牙模三维表面数据。由此提出一套摄像机刚体变换理论，即在原有的单目三维扫描系统里面同时加载两个摄像机，每个摄像机独立完成三维扫描，然后以其中一台摄像机为坐标原点，通过标定点的三维坐标信息，完成摄像机位置的刚体变换过程求解，最终将扫描的数据经过这个刚体变换，得到在同一个坐标系下的数据信息，即牙模三维表面数据信息。

此外在实现牙模三维表面数据获取过程中，分别研究了以下三个主要的分课题，这些课题的成功实现是开发完成这套高精度扫描系统的基础。

1、在机器视觉的图像处理中，一般都需要进行摄像机的标定工作。本系统研究了基于OpenCV的摄像机标定，根据OpenCV提供的棋盘标定法获取摄像机的内外参数，进而完成了图像的校正工作，是后续的背景标定以及投影几何变换的基础。

2、在构造单目机器视觉几何背景中，用到了一个原创设计的标定板，用以获取场景三维信息。标定点在图像中的精确定位方法也是本研究的一个重点。针对所使用标定点为圆圈的特征，提出了多种约束条件，进而可以完成精确的标定点定位。

3、激光的光带中心的提取，是所有主动视觉三维测量系统中的难点问题，针对所测量物体的特性，提出一种区域阈值重心法的结构光中心提取方法。提取结构能够精确到亚像素级，提取速度能够达到实时需要。

整个三维扫描系统中，主要包括了视频流显示模块、参数标定模块、三维扫描模块这三个大的子模块。其中参数标定模块实现的功能比较多，在这个模块设计的操作交互式操作流程也比较多。

视频流显示模块，主要实现摄像机的采集和视频流的显示功能。视频的采集主要使用的是最新版的OpenCV提供的VideoCapture类，实现视频流从硬件设备的读取，Mat类存储图像数据，图像显示是用到了QT的Widget类实现。这部分主要设计了载入摄像机和播放、暂停，释放摄像机功能，还包括左右摄像机的位置对换功能。

参数标定模块，主要是交互式实现背景板标定点的提取，激光平面的提取，摄像机刚体变换的求解。这部分主要为实现三维扫描提供基本的标定数据信息。交互式操作主要是为了更好的提取标定点以及标定激光平面。此部分设计的功能操作简单，鲁棒性强。

三维扫描模块，主要实现了光条中心的实时提取与点云数据实时生成功能，同时将生成的点云数据实时写入文件。当扫描物体滑出扫描区域，自动结束当前扫描。此模块完全依赖于前面两个模块运行的结果。

此外该三维扫描系统还设计了标定数据导入导出、图像分辨率大小更改等功能。可以极大的方便用户自己定义需要扫描模式。

三维扫描系统的开发环境为VS2010搭载QT5.0、OpenCV2.4。硬件环境：1mw632nm的一字线激光发射器、高精度平移台、罗技9000网络摄像机。系统的开发全部为C\C++语言实现。

参见图1，充分利用针孔投影模型的几何关系，构建了这套高效的三维扫描系统。左右摄像机分别作为独立的单目扫描仪工作。左右摄像机为参数相同的摄像机。激光发射器为一字线激光发射器。平移台上放置扫描物体。这里使用的平移台的平移速度可以提前读出。此外还需要将平移台放置在一个水平面上，在激光平面标定中需要用到平移台的底座平面。在搭建本套扫描系统时，左右摄像机需要有一定偏角，这样便于求解摄像机的刚体变换。激光发射器与摄像机的位置间留有一段距离。平移台的台面到台座的垂直距离需要提前手动测量一次（作为固定参数保存）。

参见图2，在三维扫描系统开始运行之前需要完成摄像机的标定，获取摄像机的焦距、成像原点、畸变系数等信息。初始化系统阶段从文件中读入左右摄像机的校正数据，以及原来的平移台面和激光的标定数据。图像获取阶段，自动校正获取的图像。参数获取阶段，选择激光标定和平移台面标定，或则直接使用以前标定的数据，获取左右摄像机的刚体变换矩阵，同时获取平移台平移速度。开始扫描阶段，启动平移台，确定扫描数据点存入位置。扫描阶段，提取光条中心，完成三维投影变换得到点云数据。数据处理阶段，将扫描得到的数据经过刚体变换合并左右摄像机的扫描数据。当光条移除扫描区域则系统自动结束扫描，或者手动结束完成扫描。

参见图3，整个三维扫描系统主要包括了视频流显示模块、参数标定模块、三维扫描模块这三个大的子模块。其中参数标定模块实现的功能比较多，在这个模块设计的操作交互式操作流程也比较多。

参见图4，本发明的标定板为三列水平方向等距的三个黑圆圈，标定板放置在平移台的台面里。标定板内的标定点分别按从左到右，从上到下的顺序排列，设所有标定点在图像中的位置坐标为PX{P1,P2,P3,…P9}。

参见图5，本发明充分利用了扫描平移台面与台座两个平面，由求解得到台面的标定点坐标，根据三维空间中非共线的三点确定唯一平面，即可得到台面的平面方程。

标定点的提取，目前大多数标定点都采用角点检测的方法，本发明参照现有的资料，经过大量的实验得出了一种精确的标定点提取方法，该方法提取精度高，可靠性好。

方法描述如下：

1、首先读入图像，对图形进行灰度化、高斯滤波。

2、运用Canny算子提取轮廓。

3、提取轮廓信息，用椭圆逼近算法，提取轮廓的最小逼近椭圆。

4、根据OpenCV提供的函数接口，可以取得轮廓的周长和面积。

根据圆圈的周长和面积计算公式：

周长：C=2*π*r

面积：Area=π*r²

则可以得出一个约束条件：

5、取得逼近椭圆的中心点，判断该点出的像素值大小，给定该阈值<127。判定该标定点为黑色中心点。

通过上述的约束条件，可以大大的提高标定点的提取的准确度和精度。提取效果参见图7。

在实验阶段，在分辨率为640*480，扫描了牙模模型。为了方便扫描数据的可视化和全方位查看，设计了一个伪彩色点云数据查看器。扫描结果显示图都是通过这个查看器显示的结果。参见图8的扫描结果都以深度伪彩色显示，由近到远分别用红-绿-蓝显示深度信息。

Claims (6)

1.一种牙模三维表面数据非接触三维扫描方法，其特征在于，该方法采用激光发射器、高精度平移台、和两台网络摄像机构成的三维扫描系统，每个摄像机独立完成三维扫描，然后以其中一台摄像机为坐标原点，通过标定点的三维坐标信息，完成摄像机位置的刚体变换过程求解，最终将扫描的数据经过这个刚体变换，得到在同一个坐标系下的牙模三维表面数据信息，具体按下列步骤操作：

1）制作标定板，标定板为三列水平方向等距的三个黑圆圈，标定板放置在高精度平移台的台面上，标定板内的标定点分别按从左到右，从上到下的顺序排列，获取所有标定点在图像中的位置坐标；

2）两个摄像机对称设置在高精度平移台两侧，分别作为独立的单目扫描仪工作，两个摄像机与高精度平移台的台座平面有一定偏角，便于求解摄像机的刚体变换；

3）根据三维空间中非共线的三点确定唯一平面，即可得到台面的平面方程；

4）借用激光平面与物体表面相交求解物体表面三维空间位置坐标值，利用平移滑台的匀速运动的特性，添加物体平移值，得到物体的整个表面的三维点云坐标；

5）在扫描过程中，将激光条中心一点在图像中的坐标点转换为成像平面上一点，该点与坐标原点构成的直线与激光平面相交于一点，这点为激光线打在物体上面的物理空间中的点；

6）高精度平移台左侧摄像机经过选择和平移变换，得到与右摄像机的扫描数据在同一坐标系下的完整的三维扫描数据。

2.如权利要求1所述的方法，其特征在于，所述的坐标点为9个，其中每三个坐标点均匀分布在一条直线上，形成三条相互平行的直线。

3.如权利要求2所述的方法，其特征在于，所述的三维扫描数据在处理阶段，将扫描得到的数据经过刚体变换合并左右摄像机的扫描数据。

4.如权利要求1所述方法，其特征在于，在图像中激光光条上中心线上，在台座平面内找一个点，在平移台面内找两点；这三点的三维坐标可通过相应公式计算得到，由这三个空间点即可确定激光平面方程。

5.如权利要求1所述的方法，其特征在于，所述的标定点的获取方法包括以下步骤：

第一步，读入图像，对图形进行灰度化、高斯滤波；

第二步，运用Canny算子提取轮廓；

第三步，提取轮廓信息，用椭圆逼近算法，提取轮廓的最小逼近椭圆。

第四步，根据OpenCV提供的函数接口，可以取得轮廓的周长和面积，由此可以得出一个约束条件：

$\mathrm{Thred}=\frac{C^2}{\mathrm{Area}}=4*\mathrm{\&pi;};$

第五步，取得逼近椭圆的中心点，判断该点出的像素值大小，给定该阈值<127，判定该标定点为黑色中心点。

6.如权利要求1所述的方法，其特征在于，所述的激光条中心采用阈值重心法，包括以下步骤：

首先在图像的空间域中首先给定一个阈值，当图像中某个位置的亮度值大于设定的阈值时提取该点，否则跳过该点，进行下一个；

然后再将提取到的点运用重心法进行计算就得到光条中心结构光中心点；

最后按照图像的高斯分布法则，在上述提取的中心点上，对每个点进行Hessian矩阵运算，得到中心点的法线，最终得到光条中心的亚像素精确点。

Images