CN101697233B - 一种基于结构光的三维物体表面重建方法 - Google Patents
一种基于结构光的三维物体表面重建方法 Download PDFInfo
- Publication number
- CN101697233B CN101697233B CN2009102177397A CN200910217739A CN101697233B CN 101697233 B CN101697233 B CN 101697233B CN 2009102177397 A CN2009102177397 A CN 2009102177397A CN 200910217739 A CN200910217739 A CN 200910217739A CN 101697233 B CN101697233 B CN 101697233B
- Authority
- CN
- China
- Prior art keywords
- camera
- projection screen
- projector
- point
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提出一种基于结构光的三维物体表面重建方法，此方法使用标定板对投影机和照相机进行标定，并获得相应参数矩阵；通过投影机向投影屏幕投影格雷编码图像，照相机摄取格雷编码图像；通过对图像的预处理，得到投影机上的每个像素点与其在照相机成像平面上的位置对应关系，然后再对被测物体的调制后的图像进行预处理，对调制前和调制后的图像中的扫描点进行匹配，计算出被测物体表面的世界坐标。最后根据所得到的待测物体的三维坐标信息绘制相应的点云图像。此方法广泛适用于与建模、检测技术相关的虚拟现实、军事仿真、工业设计、荧幕检测等领域中各类立体模型表面的三维重建中，具有高精度、高效率、低成本、易实施等优点。
Description
技术领域
本发明涉及一种基于结构光的三维物体表面重建方法，属于计算机科学中图像处理技术领域。此方法广泛适用于与建模、检测技术相关的虚拟现实、军事仿真、工业设计、荧幕检测等领域中各类立体模型表面的三维重建中，具有高精度、高效率、低成本、易实施等优点。
背景技术
基于结构光的三维测量系统不同于传统的立体视觉系统，它用一个投射编码的设备代替了立体视觉中的一个摄像机，避免了立体匹配问题。结构光法具有精度高、原理简单、适于实时测量等优点，被广泛应用于自动加工、高速在线检测、质量控制、CAD/CAM、医学诊断、航空航天、汽车制造实物仿型、服装加工、鞋模等领域，是反向工程和计算机视觉中的重要组成部分。
基于结构光的三维物体表面重构包括投影机与照相机的标定、格雷编解码方法和重构方法。Y.F.Li等和Huynh提出了对结构光系统自动校准算法的研究；Giovanna Sansoni等人提出了一种对非连续表面测量结果进行分析估计的校正算法，具有较好的线性性质和测量精度。王宝光等人讨论了结构光传感器的数学模型及成像公式，以及结构光传感器参量对检测性能的影响，同时提出了提高线结构光传感器测量分辨力的光学图像拆分与合成技术，成功地应用于BGA芯片管脚的三维测量；张广军等人利用了RBF网络的非线性映射能力建立了结构光三维视觉检测的RBF网络模型；张宏伟等人讨论了一种结合三角Bezier曲面造型和传统NURBS曲面造型技术的实用曲面重构算法，同时提出了多视图数据的曲面匹配算法；翟广泉等人对具有多解特性的平面立体三维重建提出了新算法；刘志刚等人在测量与CAD、CMM和I-DEAS建模方面开展了研究，提出了一种基于特征线框的线性视觉传感器测量数据的分割与建模方法，并在I-DEAS中建立了CAD模型。如何实现高精度、高效率、低成本、易实施的三维重构至关重要。
发明内容
为了实现高精度、高效率、低成本、易实施的三维物体表面重建，本发明提供了一种基于结构光的三维物体表面重建方法。此方法使用标定板对投影机和照相机进行标定，并获得相应参数矩阵；通过投影机向投影屏幕投影格雷编码图像，照相机摄取格雷编码图像；通过对图像的预处理，得到投影机上的每个像素点与其在照相机成像平面上的位置对应关系，然后再对被测物体的调制后的图像进行预处理，通过对调制前和调制后的图像中的扫描点进行匹配，计算出被测物体表面点的世界坐标，根据所得到的待测物体的三维坐标信息绘制相应的点云图像。
一种基于结构光的三维物体表面重建方法，其特征在于，该方法所需设备和实施步骤如下：
1.所需设备：如图1所示，所需设备由8部分组成：投影机(1)，待测物体(2)，投影屏幕(3)，照相机(4)，计算机(5)，格雷编码图像序列(6)，标定板(7)，标定网格图像(8)；其中格雷编码图像序列(6)由20张图像组成，包括10张横向编码图像和10张纵向编码图像，如图2所示；格雷编码图像序列(6)存储在计算机(5)中。
2.实现步骤：
步骤1：对投影机(1)和照相机(4)进行标定，得到投影机(1)的光心的世界坐标和照相机(4)的内参数矩阵、外参数矩阵；
步骤110：开始；
步骤120：将标定板(7)平行于投影屏幕(3)，位于投影屏幕(3)和照相机(4)之间，同时也位于投影屏幕(3)和投影机(1)之间；计算机(5)读取照相机(4)所拍摄的标定板(7)的图像；
定义照相机(4)的CCD所在平面为成像平面；
定义图像坐标系如下：以成像平面为图像坐标系，其中：CCD的几何中心为图像坐标系的原点，成像平面水平向右为X轴的正方向，成像平面竖直向下为Y轴的正方向，单位长度为1个像素；
定义世界坐标系如下：以投影屏幕(3)所在的平面为世界坐标系的XY平面，其中：左上角的第一个角点处为世界坐标系的原点，垂直于投影屏幕(3)所在的平面并且指向投影机(1)一侧为Z轴的负方向，在投影屏幕(3)所在的平面内水平向右为X轴的正方向，在投影屏幕(3)所在的平面内竖直向下为Y轴的正方向，单位长度为1毫米；
步骤130：计算机(5)计算标定板(7)的图像上的各个角点的图像坐标和相应的世界坐标；
步骤140：计算机(5)根据步骤130中得到的每个角点的图像坐标和世界坐标，提取其对应关系，利用张正友两步标定法求出照相机(4)的内参数矩阵和外参数矩阵；
步骤150：计算机(5)将标定网格图像(8)通过投影机(1)投放到投影屏幕(3)上显示；
步骤160：照相机(4)拍摄投影屏幕(3)上的投影显示图像；
步骤170：计算机(5)计算步骤160中照相机(4)拍摄得到的投影显示图像上的各个角点的图像坐标和相应的世界坐标；
步骤180：根据步骤170得到的每个角点的图像坐标和世界坐标，提取其对应关系，利用张正友两步标定法求解出投影机(1)的光心的世界坐标、外参数矩阵；
步骤2：对调制前与调制后的条纹编码进行匹配；
步骤210：使用投影机(1)将格雷编码图像序列(6)中的20张图像依次投影到投影屏幕(3)上，每次投影后，使用照相机(4)拍摄投影屏幕(3)上的投影图像；
步骤220：将待测物体(2)放到投影屏幕(3)的前面，保证待测物体(2)在投影机(1)和照相机(4)的视场范围内；
步骤230：使用投影机(1)将格雷编码图像序列(6)中的20张图像依次投影到投影屏幕(3)上，每次投影后，使用照相机(4)拍摄投影屏幕(3)上的经过待测物体(2)调制后的投影图像；
步骤240：计算机(5)使用高斯法对照相机(4)拍摄的图像依次进行去噪处理，灰度化，二值化，并提出条纹边界曲线；
步骤250：计算机(5)对照相机(4)拍摄的图像进行解码操作，根据调制后的格雷编码图像序列(6)中相邻码间的固定海明距离特性，消除解码过程中出现的错误点，并找到调制后的格雷编码图像序列(6)中每条横条纹与调制前相应的横条纹之间的匹配关系以及调制后的格雷编码图像序列(6)中每条竖条纹与调制前相应的竖条纹之间的匹配关系；
步骤3：根据步骤2中得到的调制后的格雷编码图像序列(6)中每条横条纹与调制前相应的横条纹之间的匹配关系以及调制后的格雷编码图像序列(6)中每条竖条纹与调制前相应的竖条纹之间的匹配关系来计算待测物体(2)的表面点的世界坐标；
步骤310：在世界坐标系中，过投影机(1)的光心A1做垂直于投影屏幕(3)的垂线、垂足为A2，设待测物体(2)的表面点A5调制前在投影屏幕(3)上的投影点为A6，过A5做垂直于投影屏幕(3)的垂线、垂足为A4；
在世界坐标系中，做照相机(4)的光心A9与成像平面的几何中心A8的连线，设待测物体(2)的表面点A5与照相机(4)的光心A9连线交成像平面上于点A7，点A9与点A5的连线的延长线交投影屏幕(3)于点A3；
步骤320：将线段A1A2、A1A6、A2A6、A4A5、A9A3、A9A8、A7A8投影到与投影屏幕(3)垂直的水平面上，过A5点做与线段A7A8平行的直线交线段A2A6于点A10；得到两个分别以A1、A2、A6和A5、A4、A6为顶点的相似三角形；还得到角度关系：角A4A3A5＝角A8A7A9-角A3A10A5；
步骤330：由步骤320中得到的两个分别以A1、A2、A6和A5、A4、A6为顶点的相似三角形的相似关系和以A3、A4、A5为顶点的直角三角形的三角关系，列出方程式：
(A1A2-A4A5)/A1A2＝A2A4/A2A6；
A4A5/A3A4＝tan(角A8A7A9-角A3A10A5)；
其中：
A1A2：投影机(1)的光心A1到投影屏幕(3)的垂直距离；
A2A6：投影屏幕(3)上条纹宽度与A6到A2条纹数的乘积；
A2A3：投影屏幕(3)上条纹宽度与A3到A2条纹数的乘积；
角A3A10A5：标定时得到；
角A8A7A9：A7点到A8点的距离与焦距A8A9的比值的反余切值；
A4A5：待测物体(2)的表面点A5到投影屏幕(3)的垂直距离，即待测物体(2)的表面点A5的Z轴世界坐标值ZA5(未知)；
A2A4：A6到A2条纹数与A5点所在面的条纹宽度的乘积(未知)；
A3A4：A2A4-A2A3；
通过以上算式，计算出待测物体(2)的表面点A5的Z轴世界坐标值ZA5；
步骤340：由照相机内参数矩阵和外参数矩阵的乘积得到3行4列的线性变换矩阵M(其元素用Mij表示，其中1≤i≤3，1≤j≤4)，设点A7的图像坐标为(u，v)，根据照相机标定原理中的针孔线性模型得到两个方程：
u＝XA5*M11+YA5*M12+ZA5*M13+M14-u*XA5*M31-u*YA5*M32-u*ZA5*M33；
v＝XA5*M21+YA5*M22+ZA5*M23+M24-v*XA5*M31-v*YA5*M32-v*ZA5*M33；
通过这两个方程求出待测物体(2)的表面点A5的X轴世界坐标XA5和Y轴世界坐标YA5，得到待测物体(2)的表面点A5的三维世界坐标(XA5，YA5，ZA5)；
步骤350：根据步骤340中得到的待测物体(2)的表面点A5的三维坐标信息绘制相应的点云图像。
有益效果
本发明的目的是提供一种方便、自动化的三维物体表面重构方法。其优点在于使用计算机软件的方法，避免了手工操作，提高了重构精度；对投影机和照相机的位置没有严格的限制，避免了使用昂贵和复杂的专业光学平台，操作简单，易于实施。
附图说明
图1是一种基于结构光的三维物体表面重建方法构成图。此图也是说明书摘要附图。其中：1为投影机，2为待测物体，3为投影屏幕，4为照相机，5为计算机，6为格雷编码图像序列，7为标定板，8为标定网格图像。
图2是格雷编码图像序列图。
图3是重构算法平面示意图。
图4是待测物体重构效果图。左图为实物图，右图为本方法构造效果图。
具体实施方式
一种基于结构光的三维物体表面重建方法，该方法所需设备和实施步骤如下：
1.所需设备：如图1所示，所需设备由8部分组成：投影机(1)，待测物体(2)，投影屏幕(3)，照相机(4)，计算机(5)，格雷编码图像序列(6)，标定板(7)，标定网格图像(8)；其中格雷编码图像序列(6)由20张图像组成，包括10张横向编码图像和10张纵向编码图像，如图2所示；格雷编码图像序列(6)存储在计算机(5)中。
2.实现步骤：
步骤1：对投影机(1)和照相机(4)进行标定，得到投影机(1)的光心的世界坐标和照相机(4)的内参数矩阵、外参数矩阵；
步骤110：开始；
步骤120：将标定板(7)平行于投影屏幕(3)，位于投影屏幕(3)和照相机(4)之间，同时也位于投影屏幕(3)和投影机(1)之间；计算机(5)读取照相机(4)所拍摄的标定板(7)的图像；
定义照相机(4)的CCD所在平面为成像平面；
定义图像坐标系如下：以成像平面为图像坐标系，其中：CCD的几何中心为图像坐标系的原点，成像平面水平向右为X轴的正方向，成像平面竖直向下为Y轴的正方向，单位长度为1个像素；
定义世界坐标系如下：以投影屏幕(3)所在的平面为世界坐标系的XY平面，其中：左上角的第一个角点处为世界坐标系的原点，垂直于投影屏幕(3)所在的平面并且指向投影机(1)一侧为Z轴的负方向，在投影屏幕(3)所在的平面内水平向右为X轴的正方向，在投影屏幕(3)所在的平面内竖直向下为Y轴的正方向，单位长度为1毫米；
步骤130：计算机(5)计算标定板(7)的图像上的各个角点的图像坐标和相应的世界坐标；
步骤140：计算机(5)根据步骤130中得到的每个角点的图像坐标和世界坐标，提取其对应关系，利用张正友两步标定法求出照相机(4)的内参数矩阵和外参数矩阵；
步骤150：计算机(5)将标定网格图像(8)通过投影机(1)投放到投影屏幕(3)上显示；
步骤160：照相机(4)拍摄投影屏幕(3)上的投影显示图像；
步骤170：计算机(5)计算步骤160中照相机(4)拍摄得到的投影显示图像上的各个角点的图像坐标和相应的世界坐标；
步骤180：根据步骤170得到的每个角点的图像坐标和世界坐标，提取其对应关系，利用张正友两步标定法求解出投影机(1)的光心的世界坐标、外参数矩阵；
步骤2：对调制前与调制后的条纹编码进行匹配；
步骤210：使用投影机(1)将格雷编码图像序列(6)中的20张图像依次投影到投影屏幕(3)上，每次投影后，使用照相机(4)拍摄投影屏幕(3)上的投影图像；
步骤220：将待测物体(2)放到投影屏幕(3)的前面，保证待测物体(2)在投影机(1)和照相机(4)的视场范围内；
步骤230：使用投影机(1)将格雷编码图像序列(6)中的20张图像依次投影到投影屏幕(3)上，每次投影后，使用照相机(4)拍摄投影屏幕(3)上的经过待测物体(2)调制后的投影图像；
步骤240：计算机(5)使用高斯法对照相机(4)拍摄的图像依次进行去噪处理，灰度化，二值化，并提出条纹边界曲线；
步骤250：计算机(5)对照相机(4)拍摄的图像进行解码操作，根据调制后的格雷编码图像序列(6)中相邻码间的固定海明距离特性，消除解码过程中出现的错误点，并找到调制后的格雷编码图像序列(6)中每条横条纹与调制前相应的横条纹之间的匹配关系以及调制后的格雷编码图像序列(6)中每条竖条纹与调制前相应的竖条纹之间的匹配关系；
步骤3：根据步骤2中得到的调制后的格雷编码图像序列(6)中每条横条纹与调制前相应的横条纹之间的匹配关系以及调制后的格雷编码图像序列(6)中每条竖条纹与调制前相应的竖条纹之间的匹配关系来计算待测物体(2)的表面点的世界坐标；
步骤310：在世界坐标系中，过投影机(1)的光心A1做垂直于投影屏幕(3)的垂线、垂足为A2，设待测物体(2)的表面点A5调制前在投影屏幕(3)上的投影点为A6，过A5做垂直于投影屏幕(3)的垂线、垂足为A4；
在世界坐标系中，做照相机(4)的光心A9与成像平面的几何中心A8的连线，设待测物体(2)的表面点A5与照相机(4)的光心A9连线交成像平面上于点A7，点A9与点A5的连线的延长线交投影屏幕(3)于点A3；
步骤320：将线段A1A2、A1A6、A2A6、A4A5、A9A3、A9A8、A7A8投影到与投影屏幕(3)垂直的水平面上，过A5点做与线段A7A8平行的直线交线段A2A6于点A10；得到两个分别以A1、A2、A6和A5、A4、A6为顶点的相似三角形；还得到角度关系：角A4A3A5＝角A8A7A9-角A3A10A5；
步骤330：由步骤320中得到的两个分别以A1、A2、A6和A5、A4、A6为顶点的相似三角形的相似关系和以A3、A4、A5为顶点的直角三角形的三角关系，列出方程式：
(A1A2-A4A5)/A1A2＝A2A4/A2A6；
A4A5/A3A4＝tan(角A8A7A9-角A3A10A5)；
其中：
A1A2：投影机(1)的光心A1到投影屏幕(3)的垂直距离；
A2A6：投影屏幕(3)上条纹宽度与A6到A2条纹数的乘积；
A2A3：投影屏幕(3)上条纹宽度与A3到A2条纹数的乘积；
角A3A10A5：标定时得到；
角A8A7A9：A7点到A8点的距离与焦距A8A9的比值的反余切值；
A4A5：待测物体(2)的表面点A5到投影屏幕(3)的垂直距离，即待测物体(2)的表面点A5的Z轴世界坐标值ZA5(未知)；
A2A4：A6到A2条纹数与A5点所在面的条纹宽度的乘积(未知)；
A3A4：A2A4-A2A3；
通过以上算式，计算出待测物体(2)的表面点A5的Z轴世界坐标值ZA5；
步骤340：由照相机内参数矩阵和外参数矩阵的乘积得到3行4列的线性变换矩阵M(其元素用Mij表示，其中1≤i≤3，1≤j≤4)，设点A7的图像坐标为(u，v)，根据照相机标定原理中的针孔线性模型得到两个方程：
u＝XA5*M11+YA5*M12+ZA5*M13+M14-u*XA5*M31-u*YA5*M32-u*ZA5*M33；
v＝XA5*M21+YA5*M22+ZA5*M23+M24-v*XA5*M31-v*YA5*M32-v*ZA5*M33；
通过这两个方程求出待测物体(2)的表面点A5的X轴世界坐标XA5和Y轴世界坐标YA5，得到待测物体(2)的表面点A5的三维世界坐标(XA5，YA5，ZA5)；
步骤350：根据步骤340中得到的待测物体(2)的表面点A5的三维坐标信息绘制相应的点云图像。
本方法已在工程试验中得到应用，取得了好的效果。
Claims (2)
1.一种基于结构光的三维物体表面重建方法，其特征在于，该方法所需设备和实施步骤如下：
所需设备由8部分组成：投影机(1)，待测物体(2)，投影屏幕(3)，照相机(4)，计算机(5)，格雷编码图像序列(6)，标定板(7)，标定网格图像(8)；其中格雷编码图像序列(6)由20张图像组成，包括10张横向编码图像和10张纵向编码图像；格雷编码图像序列(6)存储在计算机(5)中。
实现步骤：
步骤1：对投影机(1)和照相机(4)进行标定，得到投影机(1)的光心的世界坐标和照相机(4)的内参数矩阵、外参数矩阵；
步骤110：开始；
步骤120：将标定板(7)平行于投影屏幕(3)，位于投影屏幕(3)和照相机(4)之间，同时也位于投影屏幕(3)和投影机(1)之间；计算机(5)读取照相机(4)所拍摄的标定板(7)的图像；
定义照相机(4)的CCD所在平面为成像平面；
定义图像坐标系如下：以成像平面为图像坐标系，其中：CCD的几何中心为图像坐标系的原点，成像平面水平向右为X轴的正方向，成像平面竖直向下为Y轴的正方向，单位长度为1个像素；
定义世界坐标系如下：以投影屏幕(3)所在的平面为世界坐标系的XY平面，其中：左上角的第一个角点处为世界坐标系的原点，垂直于投影屏幕(3)所在的平面并且指向投影机(1)一侧为Z轴的负方向，在投影屏幕(3)所在的平面内水平向右为X轴的正方向，在投影屏幕(3)所在的平面内竖直向下为Y轴的正方向，单位长度为1毫米；
步骤130：计算机(5)计算标定板(7)的图像上的各个角点的图像坐标和相应的世界坐标；
步骤140：计算机(5)根据步骤130中得到的每个角点的图像坐标和世界坐标，提取其对应关系，利用张正友两步标定法求出照相机(4)的内参数矩阵和外参数矩阵；
步骤150：计算机(5)将标定网格图像(8)通过投影机(1)投放到投影屏幕(3)上显示；
步骤160：照相机(4)拍摄投影屏幕(3)上的投影显示图像；
步骤170：计算机(5)计算步骤160中照相机(4)拍摄得到的投影显示图像上的各个角点的图像坐标和相应的世界坐标；
步骤180：根据步骤170得到的每个角点的图像坐标和世界坐标，提取其对应关系，利用张正友两步标定法求解出投影机(1)的光心的世界坐标、外参数矩阵；
步骤2：对调制前与调制后的条纹编码进行匹配；
步骤210：使用投影机(1)将格雷编码图像序列(6)中的20张图像依次投影到投影屏幕(3)上，每次投影后，使用照相机(4)拍摄投影屏幕(3)上的投影图像；
步骤220：将待测物体(2)放到投影屏幕(3)的前面，保证待测物体(2)在投影机(1)和照相机(4)的视场范围内；
步骤230：使用投影机(1)将格雷编码图像序列(6)中的20张图像依次投影到投影屏幕(3)上，每次投影后，使用照相机(4)拍摄投影屏幕(3)上的经过待测物体(2)调制后的投影图像；
步骤240：计算机(5)使用高斯法对照相机(4)拍摄的图像依次进行去噪处理，灰度化，二值化，并提出条纹边界曲线；
步骤250：计算机(5)对照相机(4)拍摄的图像进行解码操作，根据调制后的格雷编码图像序列(6)中相邻码间的固定海明距离特性，消除解码过程中出现的错误点，并找到调制后的格雷编码图像序列(6)中每条横条纹与调制前相应的横条纹之间的匹配关系以及调制后的格雷编码图像序列(6)中每条竖条纹与调制前相应的竖条纹之间的匹配关系；
步骤3：根据步骤2中得到的调制后的格雷编码图像序列(6)中每条横条纹与调制前相应的横条纹之间的匹配关系以及调制后的格雷编码图像序列(6)中每条竖条纹与调制前相应的竖条纹之间的匹配关系来计算待测物体(2)的表面点的世界坐标；
步骤310：在世界坐标系中，过投影机(1)的光心A1做垂直于投影屏幕(3)的垂线、垂足为A2，设待测物体(2)的表面点A5调制前在投影屏幕(3)上的投影点为A6，过A5做垂直于投影屏幕(3)的垂线、垂足为A4；
在世界坐标系中，做照相机(4)的光心A9与成像平面的几何中心A8的连线，设待测物体(2)的表面点A5与照相机(4)的光心A9连线交成像平面上于点A7，点A9与点A5的连线的延长线交投影屏幕(3)于点A3；
步骤320：将线段A1A2、A1A6、A2A6、A4A5、A9A3、A9A8、A7A8投影到与投影屏幕(3)垂直的水平面上，过A5点做与线段A7A8平行的直线交线段A2A6于点A10；得到两个分别以A1、A2、A6和A5、A4、A6为顶点的相似三角形；还得到角度关系：角A4A3A5＝角A8A7A9-角A3A10A5；
步骤330：由步骤320中得到的两个分别以A1、A2、A6和A5、A4、A6为顶点的相似三角形的相似关系和以A3、A4、A5为顶点的直角三角形的三角关系，列出方程式：
(A1A2-A4A5)/A1A2＝A2A4/A2A6；
A4A5/A3A4＝tan(角A8A7A9-角A3A10A5)；
其中：
A1A2：投影机(1)的光心A1到投影屏幕(3)的垂直距离；
A2A6：投影屏幕(3)上条纹宽度与A6到A2条纹数的乘积；
A2A3：投影屏幕(3)上条纹宽度与A3到A2条纹数的乘积；
角A3A10A5：标定时得到；
角A8A7A9：A7点到A8点的距离与焦距A8A9的比值的反余切值；
A4A5：待测物体(2)的表面点A5到投影屏幕(3)的垂直距离，即待测物体(2)的表面点A5的Z轴世界坐标值ZA5；
A2A4：A6到A2条纹数与A5点所在面的条纹宽度的乘积；
A3A4：A2A4-A2A3；
通过以上算式，计算出待测物体(2)的表面点A5的Z轴世界坐标值ZA5；
步骤340：由照相机内参数矩阵和外参数矩阵的乘积得到3行4列的线性变换矩阵M，其元素用Mij表示，其中1≤i≤3，1≤j≤4，设点A7的图像坐标为(u，v)，根据照相机标定原理中的针孔线性模型得到两个方程：
u＝XA5*M11+YA5*M12+ZA5*M13+M14-u*XA5*M31-u*YA5*M32-u*ZA5*M33；
v＝XA5*M21+YA5*M22+ZA5*M23+M24-v*XA5*M31-v*YA5*M32-v*ZA5*M33；
通过这两个方程求出待测物体(2)的表面点A5的X轴世界坐标XA5和Y轴世界坐标YA5，得到待测物体(2)的表面点A5的三维世界坐标(XA5，YA5，ZA5)；
步骤350：根据步骤340中得到的待测物体(2)的表面点A5的三维坐标信息绘制相应的点云图像。
2.如权利要求1所述的一种基于结构光的三维物体表面重建方法，具有以下特征：(1)使用计算机软件的方法，避免了手工操作，提高了重构精度；(2)对投影机和照相机的位置没有严格的限制，避免了使用昂贵和复杂的专业光学平台，提高效率的同时也降低了成本；(3)操作简单，易于实施。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009102177397A CN101697233B (zh) | 2009-10-16 | 2009-10-16 | 一种基于结构光的三维物体表面重建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009102177397A CN101697233B (zh) | 2009-10-16 | 2009-10-16 | 一种基于结构光的三维物体表面重建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101697233A CN101697233A (zh) | 2010-04-21 |
| CN101697233B true CN101697233B (zh) | 2012-06-06 |
Family
ID=42142336
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009102177397A Expired - Fee Related CN101697233B (zh) | 2009-10-16 | 2009-10-16 | 一种基于结构光的三维物体表面重建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101697233B (zh) |
Families Citing this family (47)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101852980B (zh) * | 2010-06-09 | 2011-10-05 | 长春理工大学 | 一种在cave投影系统上交互式播放全景视频流的方法 |
| CN101865656B (zh) * | 2010-06-18 | 2011-12-14 | 浙江大学 | 一种使用少数共面点精确定位多摄像头系统位姿的方法 |
| EP2400261A1 (de) | 2010-06-21 | 2011-12-28 | Leica Geosystems AG | Optisches Messverfahren und Messsystem zum Bestimmen von 3D-Koordinaten auf einer Messobjekt-Oberfläche |
| CN101901501B (zh) * | 2010-07-22 | 2012-02-08 | 首都师范大学 | 一种生成激光彩色云图的方法 |
| CN102339463B (zh) * | 2010-07-22 | 2015-10-21 | 首都师范大学 | 基于激光扫描仪标定线阵相机的系统及方法 |
| JP5067477B2 (ja) * | 2010-12-28 | 2012-11-07 | カシオ計算機株式会社 | 撮像パラメータ取得装置、撮像パラメータ取得方法、及び、プログラム |
| CN102184566A (zh) * | 2011-04-28 | 2011-09-14 | 湘潭大学 | 基于微型投影仪手机平台的便携式三维扫描系统及方法 |
| US10869020B2 (en) * | 2012-06-29 | 2020-12-15 | Inb Vision Ag | Method for capturing images of a preferably structured surface of an object and device for image capture |
| US8937657B2 (en) * | 2012-07-15 | 2015-01-20 | Erik Klass | Portable three-dimensional metrology with data displayed on the measured surface |
| CN103020954B (zh) * | 2012-10-31 | 2015-04-29 | 长春理工大学 | 面向非规则表面的自适应投影系统 |
| CN103528543B (zh) * | 2013-11-05 | 2015-12-02 | 东南大学 | 一种光栅投影三维测量中的系统标定方法 |
| CN103914871B (zh) * | 2014-03-06 | 2016-06-08 | 河南农业大学 | 基于点云数据的交互式选取物体表面坐标点的方法 |
| CN104408704A (zh) * | 2014-08-25 | 2015-03-11 | 太仓中科信息技术研究院 | 一种自动化变焦镜头标定装置及标定方法 |
| CN104217429A (zh) * | 2014-08-25 | 2014-12-17 | 太仓中科信息技术研究院 | 一种摄像机标定板的设计与检测方法 |
| CN104315996B (zh) * | 2014-10-20 | 2018-04-13 | 四川大学 | 用二进制编码策略实现傅里叶变换轮廓术的方法 |
| CN104574509A (zh) * | 2015-01-26 | 2015-04-29 | 上海交通大学 | 一种由投影重建物体三维图像的方法 |
| CN104677308B (zh) * | 2015-01-30 | 2017-08-29 | 宋展 | 一种高频二值条纹的三维扫描方法 |
| CN104713497B (zh) * | 2015-03-13 | 2017-07-04 | 香港应用科技研究院有限公司 | 相位偏移校准方法、3d形状检测的方法、系统及投影系统 |
| CN105095894A (zh) * | 2015-08-06 | 2015-11-25 | 磐纹科技（上海）有限公司 | 非接触型图书扫描设备 |
| CN106500626A (zh) * | 2015-09-08 | 2017-03-15 | 东南大学 | 一种手机立体成像方法及三维成像手机 |
| CN105303609A (zh) * | 2015-11-18 | 2016-02-03 | 湖南拓视觉信息技术有限公司 | 一种三维成像和实时建模的装置和方法 |
| CN105551039B (zh) * | 2015-12-14 | 2017-12-08 | 深圳先进技术研究院 | 结构光三维扫描系统的标定方法及装置 |
| US11057608B2 (en) * | 2016-01-04 | 2021-07-06 | Qualcomm Incorporated | Depth map generation in structured light system |
| CN107036555B (zh) * | 2016-10-12 | 2019-06-11 | 昆明理工大学 | 一种交叉光轴光栅投影测量仿真系统及其实现方法 |
| KR101851303B1 (ko) | 2016-10-27 | 2018-04-23 | 주식회사 맥스트 | 3차원 공간 재구성 장치 및 방법 |
| CN108242064B (zh) * | 2016-12-27 | 2020-06-02 | 合肥美亚光电技术股份有限公司 | 基于面阵结构光系统的三维重建方法及系统 |
| US9756303B1 (en) * | 2016-12-30 | 2017-09-05 | Texas Instruments Incorporated | Camera-assisted automatic screen fitting |
| CN107170043B (zh) * | 2017-06-19 | 2019-06-18 | 电子科技大学 | 一种三维重建方法 |
| CN108269279B (zh) | 2017-07-17 | 2019-11-08 | 先临三维科技股份有限公司 | 基于单目三维扫描系统的三维重构方法和装置 |
| CN107462184B (zh) * | 2017-08-15 | 2019-01-22 | 东南大学 | 一种结构光三维测量系统的参数重标定方法及其设备 |
| CN109425292B (zh) * | 2017-08-29 | 2020-06-16 | 西安知微传感技术有限公司 | 基于一维线结构光的三维测量系统标定装置及方法 |
| CN108458671B (zh) * | 2018-03-08 | 2020-01-03 | 西安知微传感技术有限公司 | 一种线结构光三维测量系统的标定方法 |
| CN108931209B (zh) * | 2018-05-04 | 2019-12-31 | 长春理工大学 | 一种高适应性的彩色物体三维重建方法 |
| CN109491518A (zh) * | 2018-11-13 | 2019-03-19 | 宁波视睿迪光电有限公司 | 一种定位交互方法、交互装置及交互系统 |
| CN109856900B (zh) * | 2018-11-27 | 2021-05-14 | 玉晶光电(厦门)有限公司 | 光学透镜组 |
| CN109917601B (zh) * | 2019-02-13 | 2021-05-28 | 盎锐（上海）信息科技有限公司 | 基于滚动快门的相机设置方法及装置 |
| CN110037808A (zh) * | 2019-05-14 | 2019-07-23 | 苏州大学 | 基于结构光扫描的术中肝脏表面信息实时采集方法及系统 |
| CN110337674B (zh) * | 2019-05-28 | 2023-07-07 | 深圳市汇顶科技股份有限公司 | 三维重建方法、装置、设备及存储介质 |
| CN111833292A (zh) * | 2019-05-28 | 2020-10-27 | 北京伟景智能科技有限公司 | 一种基于结构光的钢筋计数方法 |
| CN110415331B (zh) * | 2019-06-14 | 2022-08-19 | 中国电建集团昆明勘测设计研究院有限公司 | 一种基于点云数据的轮廓或孔洞快速检测判别和孔洞修补方法 |
| CN110677634B (zh) | 2019-11-27 | 2021-06-29 | 成都极米科技股份有限公司 | 投影仪的梯形校正方法、装置、系统及可读存储介质 |
| CN111340891B (zh) * | 2020-02-25 | 2021-02-19 | 卡莱特（深圳）云科技有限公司 | 一种利用led屏幕进行相机标定的方法及系统 |
| CN111336949B (zh) * | 2020-04-11 | 2023-06-02 | 武汉玄景科技有限公司 | 一种空间编码结构光三维扫描方法与系统 |
| CN112710254A (zh) * | 2020-12-21 | 2021-04-27 | 珠海格力智能装备有限公司 | 对象的测量方法、系统、装置、存储介质和处理器 |
| CN114280075B (zh) * | 2021-12-28 | 2024-06-28 | 无锡维度机器视觉产业技术研究院有限公司 | 一种管类零件表面缺陷在线视觉检测系统及检测方法 |
| TWI834493B (zh) * | 2023-02-08 | 2024-03-01 | 國立成功大學 | 基於多組編碼圖案的三維重建系統及方法 |
| CN116258822B (zh) * | 2023-05-16 | 2023-08-11 | 山东捷瑞数字科技股份有限公司 | 基于元宇宙的三维引擎边界限定方法、设备及存储介质 |
-
2009
- 2009-10-16 CN CN2009102177397A patent/CN101697233B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN101697233A (zh) | 2010-04-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101697233B (zh) | 一种基于结构光的三维物体表面重建方法 | |
| US9829309B2 (en) | Depth sensing method, device and system based on symbols array plane structured light | |
| CN103400366B (zh) | 基于条纹结构光的动态场景深度获取方法 | |
| CN109118569B (zh) | 基于三维模型的渲染方法和装置 | |
| CN107223269B (zh) | 三维场景定位方法和装置 | |
| US10008005B2 (en) | Measurement system and method for measuring multi-dimensions | |
| Tian et al. | Handling occlusions in augmented reality based on 3D reconstruction method | |
| US11816829B1 (en) | Collaborative disparity decomposition | |
| CN106839977B (zh) | 基于光栅投影双目成像技术的盾构渣土体积实时测量方法 | |
| CN107993263B (zh) | 环视系统自动标定方法、汽车、标定装置及存储介质 | |
| Song et al. | DOE-based structured-light method for accurate 3D sensing | |
| Lv et al. | Build 3D Scanner System based on Binocular Stereo Vision. | |
| CN104021588A (zh) | 一种实时恢复车辆三维真实模型的系统及方法 | |
| CN106797458A (zh) | 真实对象的虚拟改变 | |
| CN104634242A (zh) | 探针补点系统及方法 | |
| CN109741382A (zh) | 一种基于Kinect V2的实时三维重建方法与系统 | |
| CN111275776A (zh) | 投影增强现实方法、装置和电子设备 | |
| CN113379815A (zh) | 基于rgb相机与激光传感器的三维重建方法、装置及服务器 | |
| CN114792345B (zh) | 一种基于单目结构光系统的标定方法 | |
| CN108895979B (zh) | 线段编码的结构光深度获取方法 | |
| CN112530020A (zh) | 三维数据的重建方法、装置、处理器和电子装置 | |
| Harjoko et al. | Low-cost 3D surface reconstruction using Stereo camera for small object | |
| CN112070844A (zh) | 结构光系统的校准方法及装置、校准工具图、设备及介质 | |
| CN111914790A (zh) | 基于双摄像头的不同场景下实时人体转动角度识别方法 | |
| CN107515844B (zh) | 字体设置方法、装置及移动设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120606 Termination date: 20131016 |