CN107255458B - 一种垂直投影光栅测量仿真系统的解算方法 - Google Patents

一种垂直投影光栅测量仿真系统的解算方法 Download PDF

Info

Publication number
CN107255458B
CN107255458B CN201710462258.7A CN201710462258A CN107255458B CN 107255458 B CN107255458 B CN 107255458B CN 201710462258 A CN201710462258 A CN 201710462258A CN 107255458 B CN107255458 B CN 107255458B
Authority
CN
China
Prior art keywords
point
grating
plane
iteration
line
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.)
Active
Application number
CN201710462258.7A
Other languages
English (en)
Other versions
CN107255458A (zh
Inventor
李文国
陈迎春
杨其乐
陈�田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN201710462258.7A priority Critical patent/CN107255458B/zh
Publication of CN107255458A publication Critical patent/CN107255458A/zh
Application granted granted Critical
Publication of CN107255458B publication Critical patent/CN107255458B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明涉及一种垂直投影光栅测量仿真系统及其实现方法,属于结构光测量技术领域。本发明以windows系统为开发平台,通过逆光线追迹法,运用计算机编程语言,编写光栅投影模块、工作台模块、图像采集模块,再模拟被测物体,实现具有光栅投影测量系统的仿真模型,其中模拟的投影机和参考平面相互垂直,并且投影机和照相机的光轴无交点。本发明仿真系统不需要投影机光轴和照相机光轴相交于一点,能够实现交叉光轴式光栅投影测量系统的仿真。

Description

一种垂直投影光栅测量仿真系统的解算方法
技术领域
本发明涉及一种垂直投影光栅测量仿真系统的解算方法,特别涉及一种利用计算机仿真现实中光栅调制过程的一种结构光投影仿真系统及其实现方法,属于结构光测量技术领域。
背景技术
随着科学技术的快速发展、社会生产率的持续提升,各个领域对测量技术的需求也在不断提高,特别是对其测量方式、测量效率和测量精度的要求变得更加苛刻。随着计算技术的。光学测量技术和数字图像处理技术快速发展,三维形貌测量方法层出不穷。结构光测量技术凭借其非接触和高精度的优点,在众多的测量方法中脱颖而出。
传统的接触式三维测量技术,需要通过机械探针接触被测物体,通过出点获得物体表面精确的坐标信息。接触式三维测量技术有很多优点如测量精度较高,对被测物体的色泽没有要求,但是也存在很大的缺点如测量过程耗时较长,机械探针在跟物体接触的过程中会造成物体表面接触变形。更重要的是有些物体表面形状太复杂导致机械探针无法抵达,对于这类物体接触式三维测量技术对其就无法进行测量。
光学测量技术是一种非接触式三维测量技术,可以弥补以上接触式三维测量技术存在的缺点。光栅投影测量技术就是一种很有效的光学三维测量技术,它是通过将光栅投影到物体表面,经过物体表面的调制光栅造成形变,接着利用图像采集器捕获变形后的光栅条纹,最后经过解调即可得出物体的表面信息。由于该技术信息量大,高效率和高精确度而备受关注,具有广阔的市场前景。
为了运用条纹光栅投影系统(Fringe projectionsystem,以下简称FPS)获取三维物体表面形貌信息,在实验之前最好采用仿真验证一下设计方案。仿真可以免除了搭建FPS平台的麻烦,而且仿真的速度更快,更灵活,而且可以方便他人使用。FPS理论验证需要建立一个高精度光栅投影系统,图像采集器,实验设备校准装置等设。因此有必要开发一个FPS仿真系统作为替换。FPS仿真系统相比搭建FPS平台,不仅具有高速、经济的有点,而且具有更高的灵活性和精度,除此之外它还能为系统参数的设置和误差的分析提供很大的便利。
许多研究人员已经提出了基于条纹投影技术的光学仿真系统。在大多文献中所使用的仿真算法都主要是运用光线追迹技术。目前光线追迹法被广泛应用与光学仿真领域,但是光线追迹法还存在许多问题:首先,由于布尔联合的问题,使三维信息减小一维,从而产生两个问题,对三维物体的操作应规范化,机器运算的有限精度带来了线面交点计算的误差。引入简单的接近正则化规则可以忽略一些小的误差,但是会产生一些副作用。其次,由于浮点运算的局限性,许多常用的求交点算法存在一些明显的数值问题。特别是光线经过自由曲面,要通过复杂的算法,对不同的情况进行讨论。
研究物体表面对光栅条纹的调制过程,对于研究光栅的调制方法、研究系统参数对测量的影响、评价算法的优劣、确定合理的技术方案和系统结构具有重要意义。但是目前存在的光栅投影仿真系统的仿真系统都是照相机光轴和参考平面垂直,而投影仪光轴相对于参考平面是倾斜的,并且要求两光轴相交于一点。本发明是垂直投影光栅投影测量仿真系统,对一些上方小下方大的物体,能够很好的避免阴影的产生,使在实际测量中能够减少处理阴影的麻烦。
发明内容
本发明提供了一种垂直投影光栅测量仿真系统的解算方法,能够很好的避免阴影的产生,使在实际测量中能够减少处理阴影的麻烦。
本发明的技术方案是:所述垂直投影光栅测量仿真系统的解算方法使用的垂直投影光栅测量仿真系统为:以windows系统为开发平台,通过逆光线追迹法,运用计算机编程语言,编写光栅投影模块、工作台模块、图像采集模块,再模拟被测物体,实现具有光栅投影测量系统的仿真模型,其中模拟的投影机和平面相互垂直,并且投影机和照相机的光轴无交点;
所述光栅投影模块用于模拟投影机生成光栅条纹图像,设置光栅调制频率和初始相位,从而将光栅条纹图像投射到被测物体表面;
所述被测物体和工作台模块用于仿真将要被测量的物体和物体所在的实验台;
所述图像采集模块用于模拟照相机采集调制过的光栅条纹,模拟光栅投影到物体表面后,光栅条纹会收到物体表面的调制,从而产生形变,通过模拟工业照相机采集形变光栅条纹图像。
一种垂直投影光栅测量仿真系统的解算方法,所述解算方法的具体步骤如下:
Step1、计算F点对应的长度和宽度;其中:F点为线段EC的延长线与参考平面R5R6R7R8的交点,E点为CCD象元,线段EC为象元E与相机光心C连线;
Step2、D点坐标的计算;其中D点为象元E与光心连线EC的延长线与平面R1R2R3R4相交点;
Step3、计算直线CD与物体表面的交点A的坐标;
Step4、计算平面R1R2R3R4上D点和物体表面上A点的相位值;
Step5、根据Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值;
Step6、根据上述计算出的内容,系统中即得到照相机拍到的图像,得到图像即仿真出了垂直投影测量系统。
进一步地,(1)、所述步骤Step1中运用相似三角形原理计算F点对应的长度和宽度,具体如下:
如图1,E点为CCD象元,每个象元的长和宽分别用Sx,Sy来表示,E点在CCD阵列上的行和列分别用i,j来表示,象元E与光心连线EC的延长线与物体表面相交于A点,与平面R1R2R3R4相交于D点,与参考平面R5R6R7R8相交于F点,F点与y轴的距离是F点对应的长度,F点与x轴距离是F点对应的宽度,根据相似三角形原理得出计算F点对应的长度和宽度的计算公式(1)、(2)
Figure GDA0002202411460000031
Figure GDA0002202411460000032
其中,i表示第i行的象元,f表示照相机焦距,lF表示F点对应的长度,bF表示F点对应的宽度;l'表示照相机距参考平面R5R6R7R8的距离;
(2)、所述步骤Step2中D点坐标的计算,具体如下:
由于D点在直线CF上,根据几何关系得出公式(3)、(4),由于D点总是在平面R1R2R3R4上,得到D(xD,yD,zD)坐标的计算公式(3)、(4)、(5):
Figure GDA0002202411460000033
Figure GDA0002202411460000034
zD=0 (5)
其中α、β分别是照相机光轴与yoz平面的夹角和xoz平面的夹角,O'为相机光轴和参考平面R5R6R7R8的交点;
(3)、所述步骤Step3中,计算直线CD与物体表面的交点A的坐标,具体步骤如下:
由于A点在直线CD上,根据相似三角形原理得出公式(6)、(7)、(8),
Figure GDA0002202411460000041
zA=f(xA,yA) (8)
其中,l为投影仪距平面R1R2R3R4距离,C(xC,yC,zC)为相机光心,其横纵坐标、Z坐标均通过人工设置,根据仿真的实际位置调动,在三维空间中线面交点不能用数学公式直接求出,因此在仿真系统的实现过程中要用数值迭代法计算A点的坐标,在三维空间中采用迭代法直接计算出交点A坐标,A点是直线CD和物体表面的交点,设想CD直线上有一点Q(xq,yq,zq),过Q点做辅助线GH垂直于平面,H为直线GH和平面的交点,G为直线GH和物体表面的交点令h1=QH,h2=GQ,其中h1,h2分别由公式(9)、公式(10)计算得到,
Figure GDA0002202411460000043
h2=f(xq,yq) (10)
其中,z=f(x,y)为曲面的函数表达式;
令D点为初始迭代点,初始步长s=1,向C点方向进行迭代搜索,每进行迭代一次得到的结果如下,
Figure GDA0002202411460000044
Figure GDA0002202411460000051
Figure GDA0002202411460000052
是n+1次迭代之后得到的Q点的横坐标和纵坐标,k表示Q点纵坐标和横坐标的比值,k=j/i,其中i,j表示E点像素的纵坐标和横坐标,
确定迭代停止点,定义三个变量d1、d2、kd,令第一次计算h1、h2的差值d1=h1-h2,对Q点迭代之后计算出h1、h2的差值d2=h1-h2,令kd=d1·d2,此过程每迭代一步就重新将d1、d2赋值一次,当Q点在DA线段上时h1>h2(QH>GH),有h1-h2<0;当Q点在AC线段上时h1>h2,有h1-h2>0;如果kd的符号在第n次迭代时发生变化,说明第n-1迭代时Q点在DA线段上,第n-1迭代时Q点在AC线段上,即说明Q点已经迭代到A点附近,Q点超过A点之后开始减小步长s的值并进行反方向迭代,令s=-s/2,然后运用上面的方法进行多次迭代,直到|s|<0.001迭代停止;
(4)、所述步骤Step4计算平面R1R2R3R4上D点和物体表面上A点的相位值,具体步骤如下:
在平面R1R2R3R4上投射的光栅栅距为λ,令O点的初始相位为0,则参考平面上的D点、被测物体表面上的A点处的光栅相位分别由公式(12)、(13)计算得到
Figure GDA0002202411460000054
Figure GDA0002202411460000055
其中,xD=OD,xA是A点X坐标值;
(5)、所述步骤Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值,具体步骤如下:
如果直线ED和物体表面之间不存在交点,则CCD阵列上的E(i,j)点表示参考平面上光栅条纹强度图像,E点强度和R1R2R3R4上的D点强度相对应,它的强度值根据公式(14)计算得到;如果直线ED和物体表面有交点,则CCD阵列上的E(i,j)点表示经过物体表面调制后的光栅条纹强度图像,E点强度和物体表面上A点强度相对应,它的强度值根据公式(15)计算得到,
Figure GDA0002202411460000061
其中,a表示背景光强,b表示投射光栅的最大强度。
本发明的有益效果是:
1、使用方便
为了运用条纹光栅投影系统(FPS)获取三维物体表面形貌信息,那么在实验之前最好采用仿真验证设计方案。仿真可以免除了搭建FPS平台的麻烦,而且仿真的速度更快,更灵活,而且可以方便他人使用。FPS理论验证需要建立一个高精度光栅投影系统,图像采集器,实验设备校准装置等。FPS仿真系统相比搭建FPS平台,不仅具有高速、经济的有点,而且具有更高的灵活性和精度,除此之外它还能为系统参数的设置和误差的分析提供很大的便利。
2、投影仪可以垂直投影
本发明是垂直投影光栅投影测量仿真系统,对一些上方小下方大的物体,能够很好的避免阴影的产生,使在实际测量中能够减少处理阴影的麻烦。
附图说明
图1为本发明系统原理图;
图2为本发明D点计算中的几何关系;
图3为本发明迭代法计算A点(直线和曲面交点)的二维示意图。
具体实施方式
实施例1:如图1-3所示,所述垂直投影光栅测量仿真系统的解算方法使用的垂直投影光栅测量仿真系统为:以windows系统为开发平台,通过逆光线追迹法,运用计算机编程语言,编写光栅投影模块、工作台模块、图像采集模块,再模拟被测物体,实现具有光栅投影测量系统的仿真模型,其中模拟的投影机和平面相互垂直,并且投影机和照相机的光轴无交点;
所述光栅投影模块用于模拟投影机生成光栅条纹图像,设置光栅调制频率和初始相位,从而将光栅条纹图像投射到被测物体表面;
所述被测物体和工作台模块用于仿真将要被测量的物体和物体所在的实验台;
所述图像采集模块用于模拟照相机采集调制过的光栅条纹,模拟光栅投影到物体表面后,光栅条纹会收到物体表面的调制,从而产生形变,通过模拟工业照相机采集形变光栅条纹图像。
采用64位win7操作系统,Intel(R)Core(TM)i5-2410M CPU@2.30处理器,4GB内存,选取C++Builder6.0为开发软件。通过逆光线追迹法,运用C++语言,通过编写光栅投影模块-projecter、模拟被测物体和工作台模块、图像采集模块-camera(其中C点是照相机的光心),实现具有光栅投影测量系统的仿真模型;其中模拟的投影机和照相机的光轴不需要有交点。实施例中仿真系统参数选取如下:CCD像素为800x600,摄像机焦距f=1.5mm,象元尺寸sx=0.0032mm,sy=0.0032mm,投影仪距平面R1R2R3R4距离l=100mm,照相机距参考平面R5R6R7R8的距离l`=150mm,投影仪安装位置:xp=0mm,yp=0mm,zp=100mm。
一种垂直投影光栅测量仿真系统的解算方法,所述解算方法的具体步骤如下:
Step1、计算F点对应的长度和宽度;其中:F点为线段EC的延长线与参考平面R5R6R7R8的交点,E点为CCD象元,线段EC为象元E与相机光心C连线;
Step2、D点坐标的计算;其中D点为象元E与光心连线EC的延长线与平面R1R2R3R4相交点;
Step3、计算直线CD与物体表面的交点A的坐标;
Step4、计算平面R1R2R3R4上D点和物体表面上A点的相位值;
Step5、根据Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值;
Step6、根据上述计算出的内容,系统中即得到照相机拍到的图像,得到图像即仿真出了垂直投影测量系统。
进一步地,(1)、所述步骤Step1中运用相似三角形原理计算F点对应的长度和宽度,具体如下:
如图1,E点为CCD象元,每个象元的长和宽分别用Sx,Sy来表示,E点在CCD阵列上的行和列分别用i,j来表示,象元E与光心连线EC的延长线与物体表面相交于A点,与平面R1R2R3R4相交于D点,与参考平面R5R6R7R8相交于F点,F点与y轴的距离是F点对应的长度,F点与x轴距离是F点对应的宽度,根据相似三角形原理得出计算F点对应的长度和宽度的计算公式(1)、(2)
Figure GDA0002202411460000081
Figure GDA0002202411460000082
其中,i表示第i行的象元,f表示照相机焦距,lF表示F点对应的长度,bF表示F点对应的宽度;l'表示照相机距参考平面R5R6R7R8的距离;
(2)、所述步骤Step2中D点坐标的计算,具体如下:
由于D点在直线CF上,根据几何关系得出公式(3)、(4),由于D点总是在平面R1R2R3R4上,得到D(xD,yD,zD)坐标的计算公式(3)、(4)、(5):
Figure GDA0002202411460000083
Figure GDA0002202411460000086
zD=0 (5)
其中α、β分别是照相机光轴与yoz平面的夹角和xoz平面的夹角,O'为相机光轴和参考平面R5R6R7R8的交点;
(3)、所述步骤Step3中,计算直线CD与物体表面的交点A的坐标,具体步骤如下:
由于A点在直线CD上,根据相似三角形原理得出公式(6)、(7)、(8),
Figure GDA0002202411460000084
Figure GDA0002202411460000085
zA=f(xA,yA) (8)
其中,l为投影仪距平面R1R2R3R4距离,C(xC,yC,zC)为相机光心,其横纵坐标、Z坐标均通过人工设置,根据仿真的实际位置调动,在三维空间中线面交点不能用数学公式直接求出,因此在仿真系统的实现过程中要用数值迭代法计算A点的坐标,在三维空间中采用迭代法直接计算出交点A坐标,A点是直线CD和物体表面的交点,设想CD直线上有一点Q(xq,yq,zq),过Q点做辅助线GH垂直于平面,H为直线GH和平面的交点,G为直线GH和物体表面的交点令h1=QH,h2=GQ,其中h1,h2分别由公式(9)、公式(10)计算得到,
Figure GDA0002202411460000091
h2=f(xq,yq) (10)
其中,z=f(x,y)为曲面的函数表达式;
令D点为初始迭代点,初始步长s=1,向C点方向进行迭代搜索,每进行迭代一次得到的结果如下,
Figure GDA0002202411460000092
Figure GDA0002202411460000093
Figure GDA0002202411460000094
Figure GDA0002202411460000095
是n+1次迭代之后得到的Q点的横坐标和纵坐标,k表示Q点纵坐标和横坐标的比值,k=j/i,其中i,j表示E点像素的纵坐标和横坐标,
确定迭代停止点,定义三个变量d1、d2、kd,令第一次计算h1、h2的差值d1=h1-h2,对Q点迭代之后计算出h1、h2的差值d2=h1-h2,令kd=d1·d2,此过程每迭代一步就重新将d1、d2赋值一次,当Q点在DA线段上时h1>h2(QH>GH),有h1-h2<0;当Q点在AC线段上时h1>h2,有h1-h2>0;如果kd的符号在第n次迭代时发生变化,说明第n-1迭代时Q点在DA线段上,第n-1迭代时Q点在AC线段上,即说明Q点已经迭代到A点附近,Q点超过A点之后开始减小步长s的值并进行反方向迭代,令s=-s/2,然后运用上面的方法进行多次迭代,直到|s|<0.001迭代停止;
(4)、所述步骤Step4计算平面R1R2R3R4上D点和物体表面上A点的相位值,具体步骤如下:
在平面R1R2R3R4上投射的光栅栅距为λ,令O点的初始相位为0,则平面上的D点、被测物体表面上的A点处的光栅相位分别由公式(12)、(13)计算得到
Figure GDA0002202411460000101
Figure GDA0002202411460000102
其中,xD=OD,xA是A点X坐标值;
(5)、所述步骤Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值,具体步骤如下:
如果直线ED和物体表面之间不存在交点,则CCD阵列上的E(i,j)点表示参考平面上光栅条纹强度图像,E点强度和平面R1R2R3R4上的D点强度相对应,它的强度值根据公式(14)计算得到;如果直线ED和物体表面有交点,则CCD阵列上的E(i,j)点表示经过物体表面调制后的光栅条纹强度图像,E点强度和物体表面上A点强度相对应,它的强度值根据公式(15)计算得到,
Figure GDA0002202411460000103
Figure GDA0002202411460000104
其中,a表示背景光强,b表示投射光栅的最大强度。
上面结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (2)

1.一种垂直投影光栅测量仿真系统的解算方法,其特征在于:所述解算方法的具体步骤如下:
Step1、计算F点对应的长度和宽度;其中:F点为线段EC的延长线与参考平面R5R6R7R8的交点,E点为CCD象元,线段EC为象元E与相机光心C连线;
Step2、D点坐标的计算;其中D点为象元E与光心连线EC的延长线与平面R1R2R3R4相交点;
Step3、计算直线CD与物体表面的交点A的坐标;
Step4、计算平面R1R2R3R4上D点和物体表面上A点的相位值;
Step5、根据Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值;
Step6、根据上述计算出的内容,系统中即得到照相机拍到的图像,得到图像即仿真出了垂直投影测量系统;
所述垂直投影光栅测量仿真系统的解算方法使用的垂直投影光栅测量仿真系统为:以windows系统为开发平台,通过逆光线追迹法,运用计算机编程语言,编写光栅投影模块、工作台模块、图像采集模块,再模拟被测物体,实现具有光栅投影测量系统的仿真模型,其中模拟的投影机和平面相互垂直,并且投影机和照相机的光轴无交点;
所述光栅投影模块用于模拟投影机生成光栅条纹图像,设置光栅调制频率和初始相位,从而将光栅条纹图像投射到被测物体表面;
所述被测物体和工作台模块用于仿真将要被测量的物体和物体所在的实验台;
所述图像采集模块用于模拟照相机采集调制过的光栅条纹,模拟光栅投影到物体表面后,光栅条纹会收到物体表面的调制,从而产生形变,通过模拟工业照相机采集形变光栅条纹图像。
2.根据权利要求1所述的垂直投影光栅测量仿真系统的解算方法,其特征在于:
(1)、所述步骤Step1中运用相似三角形原理计算F点对应的长度和宽度,具体如下:
E点为CCD象元,每个象元的长和宽分别用Sx,Sy来表示,E点在CCD阵列上的行和列分别用i,j来表示,象元E与光心连线EC的延长线与物体表面相交于A点,与平面R1R2R3R4相交于D点,与参考平面R5R6R7R8相交于F点,F点与y轴的距离是F点对应的长度,F点与x轴距离是F点对应的宽度,根据相似三角形原理得出计算F点对应的长度和宽度的计算公式(1)、(2)
Figure FDA0002202411450000022
其中,i表示第i行的象元,f表示照相机焦距,lF表示F点对应的长度,bF表示F点对应的宽度;l'表示照相机距参考平面R5R6R7R8的距离;
(2)、所述步骤Step2中D点坐标的计算,具体如下:
由于D点在直线CF上,根据几何关系得出公式(3)、(4),由于D点总是在平面R1R2R3R4上,得到D(xD,yD,zD)坐标的计算公式(3)、(4)、(5):
Figure FDA0002202411450000023
Figure FDA0002202411450000024
zD=0 (5)
其中α、β分别是照相机光轴与yoz平面的夹角和xoz平面的夹角,O'为相机光轴和参考平面R5R6R7R8的交点;
(3)、所述步骤Step3中,计算直线CD与物体表面的交点A的坐标,具体步骤如下:
由于A点在直线CD上,根据相似三角形原理得出公式(6)、(7)、(8),
Figure FDA0002202411450000025
Figure FDA0002202411450000026
zA=f(xA,yA) (8)
其中,l为投影仪距平面R1R2R3R4距离,C(xC,yC,zC)为相机光心,其横纵坐标、Z坐标均通过人工设置,根据仿真的实际位置调动,在三维空间中线面交点不能用数学公式直接求出,因此在仿真系统的实现过程中要用数值迭代法计算A点的坐标,在三维空间中采用迭代法直接计算出交点A坐标,A点是直线CD和物体表面的交点,设想CD直线上有一点Q(xq,yq,zq),过Q点做辅助线GH垂直于平面,H为直线GH和平面的交点,G为直线GH和物体表面的交点令h1=QH,h2=GQ,其中h1,h2分别由公式(9)、公式(10)计算得到,
Figure FDA0002202411450000031
h2=f(xq,yq) (10)
其中,z=f(x,y)为曲面的函数表达式;
令D点为初始迭代点,初始步长s=1,向C点方向进行迭代搜索,每进行迭代一次得到的结果如下,
Figure FDA0002202411450000032
Figure FDA0002202411450000033
Figure FDA0002202411450000034
是n+1次迭代之后得到的Q点的横坐标和纵坐标,k表示Q点纵坐标和横坐标的比值,k=j/i,其中i,j表示E点像素的纵坐标和横坐标,
确定迭代停止点,定义三个变量d1、d2、kd,令第一次计算h1、h2的差值d1=h1-h2,对Q点迭代之后计算出h1、h2的差值d2=h1-h2,令kd=d1·d2,此过程每迭代一步就重新将d1、d2赋值一次,当Q点在DA线段上时h1>h2(QH>GH),有h1-h2<0;当Q点在AC线段上时h1>h2,有h1-h2>0;如果kd的符号在第n次迭代时发生变化,说明第n-1迭代时Q点在DA线段上,第n-1迭代时Q点在AC线段上,即说明Q点已经迭代到A点附近,Q点超过A点之后开始减小步长s的值并进行反方向迭代,令s=-s/2,然后运用上面的方法进行多次迭代,直到|s|<0.001迭代停止;
(4)、所述步骤Step4计算平面R1R2R3R4或物体表面上任意点的相位值,具体步骤如下:
在平面R1R2R3R4上投射的光栅栅距为λ,令O点的初始相位为0,则平面上的D点、被测物体表面上的A点处的光栅相位分别由公式(12)、(13)计算得到
Figure FDA0002202411450000041
Figure FDA0002202411450000042
其中,xD=OD,xA是A点X坐标值;
(5)、所述步骤Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值,具体步骤如下:
如果直线ED和物体表面之间不存在交点,则CCD阵列上的E(i,j)点表示参考平面上光栅条纹强度图像,E点强度和平面R1R2R3R4上的D点强度相对应,它的强度值根据公式(14)计算得到;如果直线ED和物体表面有交点,则CCD阵列上的E(i,j)点表示经过物体表面调制后的光栅条纹强度图像,E点强度和物体表面上A点强度相对应,它的强度值根据公式(15)计算得到,
Figure FDA0002202411450000043
Figure FDA0002202411450000044
其中,a表示背景光强,b表示投射光栅的最大强度。
CN201710462258.7A 2017-06-19 2017-06-19 一种垂直投影光栅测量仿真系统的解算方法 Active CN107255458B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710462258.7A CN107255458B (zh) 2017-06-19 2017-06-19 一种垂直投影光栅测量仿真系统的解算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710462258.7A CN107255458B (zh) 2017-06-19 2017-06-19 一种垂直投影光栅测量仿真系统的解算方法

Publications (2)

Publication Number Publication Date
CN107255458A CN107255458A (zh) 2017-10-17
CN107255458B true CN107255458B (zh) 2020-02-07

Family

ID=60024048

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710462258.7A Active CN107255458B (zh) 2017-06-19 2017-06-19 一种垂直投影光栅测量仿真系统的解算方法

Country Status (1)

Country Link
CN (1) CN107255458B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109961455B (zh) * 2017-12-22 2022-03-04 杭州萤石软件有限公司 一种目标检测方法及装置
CN112116619B (zh) * 2020-09-16 2022-06-10 昆明理工大学 一种基于结构约束的多线结构光系统条纹中心线提取方法
CN112254681B (zh) * 2020-10-26 2022-06-07 昆明理工大学 一种发散型多线激光投影测量仿真系统及其实现方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1851618A (zh) * 2006-05-31 2006-10-25 北京航空航天大学 单目视觉半实物仿真系统及方法
CN103019011A (zh) * 2012-12-03 2013-04-03 哈尔滨飞羽科技有限公司 一种投影系统效果的仿真验证方法
CN104199313A (zh) * 2014-09-15 2014-12-10 中国农业科学院农业信息研究所 一种机器人智能仿真测试系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002516443A (ja) * 1998-05-15 2002-06-04 トリコーダー テクノロジー ピーエルシー 3次元表示のための方法および装置
CA2528791A1 (en) * 2005-12-01 2007-06-01 Peirong Jia Full-field three-dimensional measurement method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1851618A (zh) * 2006-05-31 2006-10-25 北京航空航天大学 单目视觉半实物仿真系统及方法
CN103019011A (zh) * 2012-12-03 2013-04-03 哈尔滨飞羽科技有限公司 一种投影系统效果的仿真验证方法
CN104199313A (zh) * 2014-09-15 2014-12-10 中国农业科学院农业信息研究所 一种机器人智能仿真测试系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于光栅投影的三维物体面形测量研究;邓萍萍等;《大众科技》;20160720;第18卷(第203期);第16-18页 *
空间相移法在双频光栅轮廓术中的应用;葛东东等;《微计算机信息》;20060910;第271-273页 *

Also Published As

Publication number Publication date
CN107255458A (zh) 2017-10-17

Similar Documents

Publication Publication Date Title
CN103530880B (zh) 基于投影高斯网格图案的摄像机标定方法
US11050983B2 (en) System and method for recalibrating a projector system
Luo et al. A simple calibration procedure for structured light system
CN107255458B (zh) 一种垂直投影光栅测量仿真系统的解算方法
CN102288131A (zh) 物体360°轮廓误差的自适应条纹测量装置和方法
Xiao et al. Optical fringe-reflection deflectometry with bundle adjustment
CN107036555B (zh) 一种交叉光轴光栅投影测量仿真系统及其实现方法
Zou et al. High-accuracy calibration of line-structured light vision sensors using a plane mirror
CN107167088A (zh) 玻璃变形量的测量方法及装置
Wei et al. Calibration approach for structured-light-stripe vision sensor based on the invariance of double cross-ratio
Zhao et al. Optical nanoscale positioning measurement with a feature-based method
CN116379965A (zh) 结构光系统标定方法、装置、结构光系统及存储介质
CN103559710B (zh) 一种用于三维重建系统的标定方法
Li et al. Monocular underwater measurement of structured light by scanning with vibrating mirrors
CN112985258B (zh) 一种三维测量系统的标定方法及测量方法
Wang et al. Positioning error calibration for two-dimensional precision stages via globally optimized image registration
CN113916128A (zh) 一种基于光笔式视觉测量系统的提高精度的方法
Jiang et al. Calibration method for binocular vision with large FOV based on normalized 1D homography
Fan et al. Near-field photometric stereo using a ring-light imaging device
CN110243308A (zh) 基于区域定位的实时相位测量轮廓术
CN112254681B (zh) 一种发散型多线激光投影测量仿真系统及其实现方法
Zhang et al. Real-time line-structured light measurement system based on Scheimpflug principle
Wang Theoretical and experimental study on a new reconstruction model based on geometric optics for a single-lens biprism-based stereovision system
Sun et al. Calibration Method for a Multi-line Structured Laser Light Vision System.
Liu et al. Comparison study of three camera calibration methods considering the calibration board quality and 3D measurement accuracy

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant