光学仪器  2017, Vol. 39 Issue (4): 1-6   PDF    
一种光学投影测量仿真方法
陈迎春, 李文国     
昆明理工大学 机电工程学院, 云南 昆明 650500
摘要: 提出一种新的光栅条纹投影轮廓术仿真系统,该系统通过坐标系平移、转换以及运用相似三角形等方法,使仿真系统不再要求投影仪光轴和照相机光轴相交于一点。克服了投影仪和照相机之间的位置约束,该系统能对更复杂的测量系统进行仿真和算法研究。不同系统参数下的仿真结果表明,该系统可以实现光轴无交点情况下的仿真,得出的阴影区域符合实际情况。
关键词: 光栅投影     测量     仿真    
A novel algorithm for fringe projection profilometry
CHEN Yingchun, LI Wenguo     
Faculty of Mechanical and Electrical Engineering, Kunming University of Science and Technology, Kunming 650500, China
Abstract: Previously, the algorithm for grating projection simulation is limited by the condition that the optical axes of camera and projector must have an intersection point.We put forward a new algorithm that is based on the translation and transformation of coordinate systems and the use of similar triangles.The algorithm does not need the same intersection point of the optical axes.The system overcomes the position constraint of the projector and the camera, and it can carry out simulation and algorithm study of more complex measurement systems.The simulation results with different system settings show that the system can realize simulation for the condition without an intersection point, and the obtained shadow is consistent with actual situation.
Key words: grating projection     measurement     simulation    
引言

三维物体表面轮廓测量在航空航天[1]、汽车制造[2]、造船[3]、精密制造[4]、逆向工程[5]、文物保护[6]等领域有着广泛的应用。光学传感器具有高精度、非接触性等优点[7], 因此被广泛应用于三维物体表面轮廓测量之中。其中相位轮廓测量术(PMP)是结构光照明型主动光学三维传感方法, 由于其量程大、精度高、实时性强等优点而被广泛应用。因此研究三维物体表面对结构光光场的空间调制过程, 对研究解调方法、研究系统参数对测量的影响、评价算法的优劣、确定合力的技术方案和系统结构具有重要的意义[8]

三维面型对结构光进行空间调制的计算机仿真系统, 可以对三维面型测量系统的性能、主要结构参数、环境条件、系统精度、算法评价进行仿真研究。近年来大量国内外学者对传统系统模型进行了研究, 许多研究人员提出了基于条纹投影技术的光学仿真系统, 所使用的仿真算法主要是运用光线追迹技术。如Liang等[9]根据相位轮廓测量术原理运用逆光线追迹技术提出的三维面型对结构光的空间调制仿真系统, 实现了对投影光场、观察光场和三维面型几何关系在不同结构参数下的仿真, 该系统运用斜率逼近的方法计算光线和物体表面的交点, 计算比较复杂, 对物体阴影的判断也只是简单的近似, 精度较低。Li等[10]提出了一种高效的系统仿真算法, 有效地提高了运算速度, 仿真系统中阴影区域的判断也从简单的近似变得更加精确, 但是依然要求照相机光轴和投影仪光轴相交于一点。由于此种仿真系统设计的局限性, 因而不能使用该系统对一些比较复杂的测量系统进行模拟仿真研究[11-12]

本文提出一种新的结构光投影仿真系统, 该系统可以在投影仪光轴和照相机光轴没有交点的情况下, 对更加复杂的光栅条纹投影测量系统进行仿真。

1 仿真原理

系统中照相机光轴和投影仪光轴没有交点, 系统原理如图 1所示, 投影仪坐标系下的光路如图 2所示。图中C点为照相机的光心, 平面R1R2R3R4是照相机的参考平面, O点为照相机光轴与其参考平面的交点, 并且令O点为世界坐标系的原点, 其中xOy面与R1R2R3R4面重合, 照相机光轴与坐标系z轴重合, 令OC=l; P点为投影仪光轴的轴心, 平面R5R6R7R8是投影仪的参考平面, O′点为投影仪光轴与其参考平面的交点, 并且令O′点为投影仪坐标系的原点, 其中xOy′面与R5R6R7R8面重合, 照相机光轴与投影仪坐标系z′轴重合, 令OP=l′; 在模拟结构光光场的调制过程时, 根据物体表面形状建立数学模型, 得到物面方程z=f(x, y)。

图 1 交叉光轴的投影仿真系统原理图 Figure 1 Principle of projection simulation system with intersected optical axes

图 2 投影仪坐标系下的光路图 Figure 2 Optical path in the projector coordinate system

该系统应用逆光线追迹的计算方法, 系统的算法构架如下:1) 计算出CCD像元中心E点在其参考平面R1R2R3R4上所对应的D点坐标; 2) 运用变步长迭代法[10]计算得到物体表面上A点坐标; 3) 对A点进行从世界坐标系到投影仪坐标系的坐标变换; 4) 计算物体表面和参考平面上的相位和光强。

1.1 D点坐标的计算

E点为CCD像元中心, 每个像元的长和宽分别用Sx, Sy来表示, E点在CCD阵列上的行数和列数分别用i, j来表示。E点与光心连线EC的延长线与物体表面相交于A点, 与参考平面R1R2R3R4相交于D点。根据相似三角形原理可以得出D点坐标计算公式为

(1)
(2)

式中:f表示照相机焦距。

1.2 A点坐标的计算

图 1所示, 由于A点在直线CD上, 根据相似三角形原理可得[10]

(3)
(4)
(5)

式中:(xA, yA, zA)为点A的三维坐标。

由于在三维空间中无法利用式(3)~(5) 直接计算交点坐标, 需通过数值迭代法求出A点的三维坐标。

1.3 AD点在投影仪坐标系中坐标的求解

由于投影仪投射出的光是发散光, 光心P与参考平面R1R2R3R4的距离l′为有限远, 所以投影仪在平面R1R2R3R4上投射的光栅周期不是线性变化的。参考平面R5R6R7R8与投影仪的数字微镜器件(DMD)阵列是平行的, 因此投影仪在参考平面R5R6R7R8上投射的光栅是等距的, 应该在平面R5R6R7R8上计算物体表面和工作台平面的光栅相位。

(6)

式中:T为坐标的四阶平移矩阵; R为世界坐标系到照相机坐标系的四阶旋转矩阵。

已知点P(xP, yP, zP), 投影仪光轴与世界坐标系xOz面、yOz面的夹角分别为θ1θ2, 投影仪光心到参考平面R5R6R7R8的距离为l′。根据世界坐标系到投影仪坐标系的转换公式[式(6)]可以得到AD两点坐标在投影仪坐标系的对应位置A′、D′。

1.4 B′、W′点在投影仪坐标系中的x坐标值

在投影仪坐标系下, 如图 2所示, A′、B′、D′、W′分别对应世界坐标系中的ABDW, 在投影仪坐标系下观察, B′点是直线PA与平面R5R6R7R8的交点, 直线BF′与y′轴垂直; W′点是直线PD′与平面R5R6R7R8的交点, 直线WG′与y′轴垂直。根据相似三角形原理可以得到式(7)、(8), 继而可求出B′、W′点的x坐标。

(7)
(8)

在平面xOy′上投射的光栅栅距为λ, 令O′点的初始相位为0, 则B′点、W′点处的光栅相位分别为

(9)
(10)

式中:BF′=xB; WG′=xW

1.5 计算CCD阵列上与物体表面和工作台表面相对应像素点的强度值

如果直线ED和物体表面有交点, 则CCD阵列上的E(i, j)点表示经过物体表面调制后的光栅条纹强度图像, E点强度和参考平面R5R6R7R8上的B′点强度相对应。如果直线ED和物体表面不存在交点, 则CCD阵列上的E(i, j)点表示参考平面上光栅条纹强度图像, E点强度和参考平面R5R6R7R8上的W′点强度相对应。B′点和W′点的强度值计算公式分别为

(11)
(12)
2 阴影区域的判断

由于在阴影区域部分相位解包裹的过程不能连续进行, 所以在进行相位解包裹时出错的地方大多是从阴影区域开始, 因此阴影对相位解包裹过程的影响不可忽略。

阴影区域是指投影仪光线照不到、但是照相机能看到的区域, 该区域的光强等于背景光强。如图 3(a)所示, 对于物体表面上的A点, 照相机可以看到, 但是投影仪的光线照不到, 所以照相机看到的A点的强度是黑色背景的光强, 此时把背景的光强赋值给A点所对应的像元。如图 3(b)所示, 对于A点照相机、投影仪都可以看到, 所以A点就是非阴影区域, A点的亮度就是条纹在该点的强度。

图 3 二维情况下的阴影判别图 Figure 3 Two-dimensional shadow diagram
3 仿真结果

仿真采用Intel(R)Core(TM)i5-2410M CPU@2.30处理器, 4 GB内存, 64位Windows 7操作系统, 选用的操作软件是C++ Builder6.0, 并在MATLAB软件中对计算结果进行重构。根据上述系统原理建立仿真程序, 并将本文采用的算法在程序中实现, 在不同系统参数下对相应算法进行一系列实验, 结果证实了该系统的正确性。图 4所示为不考虑阴影情况下对物体的仿真结果。其中图 4(a)为模拟的圆锥物体在MATLAB软件中得到的三维重构图, 其方程为

(13)
图 4 计算机模拟结果(不含阴影) Figure 4 The simulation results without shadow

图 4(b)为模拟的抛物面物体(顶部下凹)在MATLAB软件中得到的三维重构图, 其方程为

(14)

图 4(c)图 4(d)分别为该系统下对图 4(a)图 4(b)中物体的光栅投影模拟结果, 从结果可以看出, 对两个物体的仿真结构均有效, 没有出现形状畸变的现象。

图 5是投影仪在不同安装位置下对抛物面物体的阴影判断的模拟结果, 系统参数设置如下:CCD像素数为600×800, 照相机焦距f=1.5 mm, 像元尺寸Sx=0.003 2 mm, Sy=0.003 2 mm, 照相机与参考平面R1R2R3R4距离l=100 mm, 投影仪与参考平面R5R6R7R8的距离l′=150 mm。图 5(a)中投影仪安装位置:xP=-150 mm, yP=0 mm, zP=100 mm; 图 5(b)中投影仪安装位置:xP=-100 mm, yP=0 mm, zP=100 mm。从图 5(a)图 5(b)P点的x坐标减小, 即投影仪和照相机的距离减小, 其他条件不变, 结果使阴影面积减小, 符合实际情况。

图 5 抛物面物体的计算机模拟结果(含阴影) Figure 5 The simulation results of the paraboloid with shadow

从实验结果来看, 在照相机光轴和投影仪光轴没有交点的情况下, 该系统所应用的算法是有效的, 克服了投影仪和照相机之间的位置约束, 使PMP测量系统不再要求两个光轴相交于一点, 由原来的三个约束条件减少为两个约束条件。该系统能对更复杂的测量系统进行仿真和算法研究。

4 结论

完成了逆光线追迹的程序设计, 求出与CCD阵列中E点相对应的D点、A点在世界坐标系中的坐标。通过坐标转换原理, 求出D点、A点在投影仪坐标系中的坐标值, 简化了投影仪参考平面上相位的计算, 解决了传统光学仿真系统要求投影仪光轴与照相机光轴相交于一点的局限性问题, 使PMP测量系统不再要求两个光轴相交于一点, 由原来的三个约束条件减少为两个约束条件。该仿真系统可以模拟出空间位置更加复杂的测量系统。

参考文献
[1] 郭峰. 坐标测量技术在航空制造中的应用[J]. 航空制造技术, 2013(19): 26–29. DOI:10.3969/j.issn.1671-833X.2013.19.001
[2] 万军. 采用CMM与光学扫描仪的汽车覆盖件模具联合测量技术[J]. 现代制造工程, 2013(2): 101–105.
[3] 黄若波, 张杰. 基于全站仪和船舶3D设计系统的三维精度测量技术研究[J]. 造船技术, 2011(4): 14–16.
[4] 刘进, 浑明, 李文俊. 三维摄影测量技术在水轮机叶片制造中的应用[J]. 大电机技术, 2013(1): 40–43.
[5] 张宁, 吴嘉辉. 基于三维数字化测量技术的逆向工程应用研究[J]. 长春大学学报, 2013, 23(12): 1555–1560.
[6] 张序, 李兆堃, 罗小华, 等. 应用三维激光扫描测量进行文物逆向工程恢复[J]. 测绘科学, 2013, 38(6): 169–171.
[7] 安冬, 盖绍彦, 达飞鹏. 一种新的基于条纹投影的三维轮廓测量系统模型[J]. 光学学报, 2014, 34(5): 0512004.
[8] 梁晓萍, 苏显渝. 物体面形对结构光场的空间调制:计算机模拟[J]. 光电工程, 1994, 21(3): 25–32.
[9] LIANG X P, SU X Y. Computer simulation of a 3-D sensing system with structured illumination[J]. Optics and Lasers in Engineering, 1997, 27(4): 379–393. DOI:10.1016/S0143-8166(96)00038-3
[10] LI W G, FANG S P, DUAN S J. Accurate and efficient algorithm for simulation of fringe projection system[J]. Optik-International Journal for Light and Electron Optics, 2013, 124(18): 3542–3548. DOI:10.1016/j.ijleo.2012.10.057
[11] BÖTTNER T, KÄSTNER M.The virtual fringe projection system (VFPS) and neural networks[C]//OSTEN W.Fringe 2005:5th International Workshop on Automatic Processing of Fringe Patterns.Heidelberg:Springer, 2005:191-194. http://rd.springer.com/chapter/10.1007/3-540-29303-5_22
[12] HASKAMP K, KÄSTNER M, REITHMEIER E.A virtual telecentric fringe projection system[C]//OSTEN W, KUJAWINSKA M.Fringe 2009:6th International Workshop on Advanced Optical Metrology.Heidelberg:Springer, 2009:1-4. http://link.springer.com/10.1007/978-3-642-03051-2_98