CN107255458A - 一种垂直投影光栅测量仿真系统及其实现方法 - Google Patents
一种垂直投影光栅测量仿真系统及其实现方法 Download PDFInfo
- Publication number
- CN107255458A CN107255458A CN201710462258.7A CN201710462258A CN107255458A CN 107255458 A CN107255458 A CN 107255458A CN 201710462258 A CN201710462258 A CN 201710462258A CN 107255458 A CN107255458 A CN 107255458A
- Authority
- CN
- China
- Prior art keywords
- mrow
- points
- msub
- point
- grating
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 claims abstract description 39
- 230000003287 optical effect Effects 0.000 claims abstract description 39
- 238000004088 simulation Methods 0.000 claims abstract description 10
- 230000014509 gene expression Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 241000216843 Ursus arctos horribilis Species 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000006870 function Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 21
- 238000005516 engineering process Methods 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000009434 installation Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 230000009897 systematic effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000013041 optical simulation Methods 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 239000002932 luster Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring 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
Abstract
本发明涉及一种垂直投影光栅测量仿真系统及其实现方法,属于结构光测量技术领域。本发明以windows系统为开发平台,通过逆光线追迹法,运用计算机编程语言,编写光栅投影模块、工作台模块、图像采集模块,再模拟被测物体,实现具有光栅投影测量系统的仿真模型,其中模拟的投影机和参考平面相互垂直,并且投影机和照相机的光轴无交点。本发明仿真系统不需要投影机光轴和照相机光轴相交于一点,能够实现交叉光轴式光栅投影测量系统的仿真。
Description
技术领域
本发明涉及一种垂直投影光栅测量仿真系统及其实现方法,特别涉及一种利用计算机仿真现实中光栅调制过程的一种结构光投影仿真系统及其实现方法,属于结构光测量技术领域。
背景技术
随着科学技术的快速发展、社会生产率的持续提升,各个领域对测量技术的需求也在不断提高,特别是对其测量方式、测量效率和测量精度的要求变得更加苛刻。随着计算技术的。光学测量技术和数字图像处理技术快速发展,三维形貌测量方法层出不穷。结构光测量技术凭借其非接触和高精度的优点,在众多的测量方法中脱颖而出。
传统的接触式三维测量技术,需要通过机械探针接触被测物体,通过出点获得物体表面精确的坐标信息。接触式三维测量技术有很多优点如测量精度较高,对被测物体的色泽没有要求,但是也存在很大的缺点如测量过程耗时较长,机械探针在跟物体接触的过程中会造成物体表面接触变形。更重要的是有些物体表面形状太复杂导致机械探针无法抵达,对于这类物体接触式三维测量技术对其就无法进行测量。
光学测量技术是一种非接触式三维测量技术,可以弥补以上接触式三维测量技术存在的缺点。光栅投影测量技术就是一种很有效的光学三维测量技术,它是通过将光栅投影到物体表面,经过物体表面的调制光栅造成形变,接着利用图像采集器捕获变形后的光栅条纹,最后经过解调即可得出物体的表面信息。由于该技术信息量大,高效率和高精确度而备受关注,具有广阔的市场前景。
为了运用条纹光栅投影系统(Fringe projection system,以下简称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、计算参考平面和物体表面上一点的相位值;
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)
其中,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):
zD=0 (5)
其中α、β分别是照相机光轴与yoz平面的夹角和xoz平面的夹角,O'为相机光轴和参考平面R5R6R7R8的交点;
(3)、所述步骤Step3中,计算直线CD与物体表面的交点A的坐标,具体步骤如下:
由于A点在直线CD上,根据相似三角形原理得出公式(6)、(7)、(8),
(6)
zA=f(xD,yD) (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)计算得到,
h2=f(xq,yq) (10)
其中,z=f(x,y)为曲面的函数表达式;
令D点为初始迭代点,初始步长s=1,向C点方向进行迭代搜索,每进行迭代一次得到的结果如下,
是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上投射的光栅栅距为λ,令O点的初始相位为0,则参考平面上的D点、被测物体表面上的A点处的光栅相位分别由公式(12)、(13)计算得到
其中,xD=OD,xA是A点X坐标值;
(5)、所述步骤Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值,具体步骤如下:
如果直线ED和物体表面之间不存在交点,则CCD阵列上的E(i,j)点表示参考平面上光栅条纹强度图像,E点强度和参考平面R5R6R7R8上的D点强度相对应,它的强度值根据公式(14)计算得到;如果直线ED和物体表面有交点,则CCD阵列上的E(i,j)点表示经过物体表面调制后的光栅条纹强度图像,E点强度和物体表面上A点强度相对应,它的强度值根据公式(15)计算得到,
其中,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、计算参考平面和物体表面上一点的相位值;
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)
其中,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):
zD=0 (5)
其中α、β分别是照相机光轴与yoz平面的夹角和xoz平面的夹角,O'为相机光轴和参考平面R5R6R7R8的交点;
(3)、所述步骤Step3中,计算直线CD与物体表面的交点A的坐标,具体步骤如下:
由于A点在直线CD上,根据相似三角形原理得出公式(6)、(7)、(8),
(6)
zA=f(xD,yD) (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)计算得到,
h2=f(xq,yq) (10)
其中,z=f(x,y)为曲面的函数表达式;
令D点为初始迭代点,初始步长s=1,向C点方向进行迭代搜索,每进行迭代一次得到的结果如下,
是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上投射的光栅栅距为λ,令O点的初始相位为0,则参考平面上的D点、被测物体表面上的A点处的光栅相位分别由公式(12)、(13)计算得到
其中,xD=OD,xA是A点X坐标值;
(5)、所述步骤Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值,具体步骤如下:
如果直线ED和物体表面之间不存在交点,则CCD阵列上的E(i,j)点表示参考平面上光栅条纹强度图像,E点强度和参考平面R5R6R7R8上的D点强度相对应,它的强度值根据公式(14)计算得到;如果直线ED和物体表面有交点,则CCD阵列上的E(i,j)点表示经过物体表面调制后的光栅条纹强度图像,E点强度和物体表面上A点强度相对应,它的强度值根据公式(15)计算得到,
其中,a表示背景光强,b表示投射光栅的最大强度。
上面结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (3)
1.一种垂直投影光栅测量仿真系统,其特征在于:以windows系统为开发平台,通过逆光线追迹法,运用计算机编程语言,编写光栅投影模块、工作台模块、图像采集模块,再模拟被测物体,实现具有光栅投影测量系统的仿真模型,其中模拟的投影机和参考平面相互垂直,并且投影机和照相机的光轴无交点;
所述光栅投影模块用于模拟投影机生成光栅条纹图像,设置光栅调制频率和初始相位,从而将光栅条纹图像投射到被测物体表面;
所述被测物体和工作台模块用于仿真将要被测量的物体和物体所在的实验台;
所述图像采集模块用于模拟照相机采集调制过的光栅条纹,模拟光栅投影到物体表面后,光栅条纹会收到物体表面的调制,从而产生形变,通过模拟工业照相机采集形变光栅条纹图像。
2.一种权利要求1所述的垂直投影光栅测量仿真系统的实现方法,其特征在于:一种所述的垂直投影光栅测量仿真系统的实现方法,所述方法的具体步骤如下:
Step1、计算F点对应的长度和宽度;其中:F点为线段EC的延长线与参考平面R5R6R7R8的交点,E点为CCD象元,线段EC为象元E与相机光心C连线;
Step2、D点坐标的计算;其中D点为象元E与光心连线EC的延长线与平面R1R2R3R4相交点;
Step3、计算直线CD与物体表面的交点A的坐标;
Step4、计算参考平面和物体表面上一点的相位值;
Step5、根据Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值;
Step6、根据上述计算出的内容,系统中即得到照相机拍到的图像,得到图像即仿真出了垂直投影测量系统。
3.根据权利要求2所述的垂直投影光栅测量仿真系统的实现方法,其特征在于:
(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)
<mrow>
<msub>
<mi>l</mi>
<mi>F</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>i</mi>
<mo>&times;</mo>
<msub>
<mi>S</mi>
<mi>x</mi>
</msub>
<mo>&times;</mo>
<mfrac>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mi>f</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mi>F</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>j</mi>
<mo>&times;</mo>
<msub>
<mi>S</mi>
<mi>y</mi>
</msub>
<mo>&times;</mo>
<mfrac>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mi>f</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,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):
<mrow>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>l</mi>
<mi>F</mi>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mi>F</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
1
<mrow>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>l</mi>
<mi>F</mi>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&beta;</mi>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mi>F</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&beta;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
zD=0 (5)
其中α、β分别是照相机光轴与yoz平面的夹角和xoz平面的夹角,O'为相机光轴和参考平面R5R6R7R8的交点;
(3)、所述步骤Step3中,计算直线CD与物体表面的交点A的坐标,具体步骤如下:
由于A点在直线CD上,根据相似三角形原理得出公式(6)、(7)、(8),
(6)
<mrow>
<msub>
<mi>y</mi>
<mi>A</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>l</mi>
</mfrac>
<mo>+</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
zA=f(xD,yD) (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)计算得到,
<mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>l</mi>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>x</mi>
<mi>q</mi>
</msub>
<msub>
<mi>x</mi>
<mi>d</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
h2=f(xq,yq) (10)
其中,z=f(x,y)为曲面的函数表达式;
令D点为初始迭代点,初始步长s=1,向C点方向进行迭代搜索,每进行迭代一次得到的结果如下,
<mrow>
<msubsup>
<mi>x</mi>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>x</mi>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<mi>s</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>y</mi>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mi>k</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mi>x</mi>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
是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上投射的光栅栅距为λ,令O点的初始相位为0,则参考平面上的D点、被测物体表面上的A点处的光栅相位分别由公式(12)、(13)计算得到
其中,xD=OD,xA是A点X坐标值;
(5)、所述步骤Step4中的相位值计算CCD阵列上与A点、D点所对应的像素的强度值,具体步骤如下:
如果直线ED和物体表面之间不存在交点,则CCD阵列上的E(i,j)点表示参考平面上光栅条纹强度图像,E点强度和参考平面R5R6R7R8上的D点强度相对应,它的强度值根据公式(14)计算得到;如果直线ED和物体表面有交点,则CCD阵列上的E(i,j)点表示经过物体表面调制后的光栅条纹强度图像,E点强度和物体表面上A点强度相对应,它的强度值根据公式(15)计算得到,
其中,a表示背景光强,b表示投射光栅的最大强度。
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 true CN107255458A (zh) | 2017-10-17 |
CN107255458B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109961455A (zh) * | 2017-12-22 | 2019-07-02 | 杭州萤石软件有限公司 | 一种目标检测方法及装置 |
CN112116619A (zh) * | 2020-09-16 | 2020-12-22 | 昆明理工大学 | 一种基于结构约束的多线结构光系统条纹中心线提取方法 |
CN112254681A (zh) * | 2020-10-26 | 2021-01-22 | 昆明理工大学 | 一种发散型多线激光投影测量仿真系统及其实现方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002516443A (ja) * | 1998-05-15 | 2002-06-04 | トリコーダー テクノロジー ピーエルシー | 3次元表示のための方法および装置 |
CN1851618A (zh) * | 2006-05-31 | 2006-10-25 | 北京航空航天大学 | 单目视觉半实物仿真系统及方法 |
US20070206204A1 (en) * | 2005-12-01 | 2007-09-06 | Peirong Jia | Full-field three-dimensional measurement method |
CN103019011A (zh) * | 2012-12-03 | 2013-04-03 | 哈尔滨飞羽科技有限公司 | 一种投影系统效果的仿真验证方法 |
CN104199313A (zh) * | 2014-09-15 | 2014-12-10 | 中国农业科学院农业信息研究所 | 一种机器人智能仿真测试系统 |
-
2017
- 2017-06-19 CN CN201710462258.7A patent/CN107255458B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002516443A (ja) * | 1998-05-15 | 2002-06-04 | トリコーダー テクノロジー ピーエルシー | 3次元表示のための方法および装置 |
US20070206204A1 (en) * | 2005-12-01 | 2007-09-06 | Peirong Jia | Full-field three-dimensional measurement method |
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)
Title |
---|
葛东东等: "空间相移法在双频光栅轮廓术中的应用", 《微计算机信息》 * |
邓萍萍等: "基于光栅投影的三维物体面形测量研究", 《大众科技》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109961455A (zh) * | 2017-12-22 | 2019-07-02 | 杭州萤石软件有限公司 | 一种目标检测方法及装置 |
US11367276B2 (en) | 2017-12-22 | 2022-06-21 | Hangzhou Ezviz Software Co., Ltd. | Target detection method and apparatus |
CN112116619A (zh) * | 2020-09-16 | 2020-12-22 | 昆明理工大学 | 一种基于结构约束的多线结构光系统条纹中心线提取方法 |
CN112254681A (zh) * | 2020-10-26 | 2021-01-22 | 昆明理工大学 | 一种发散型多线激光投影测量仿真系统及其实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107255458B (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104809755B (zh) | 基于单幅图像文化遗存遗址的三维重建方法 | |
CN100489446C (zh) | 基于相位法的三维轮廓测量方法 | |
CN101697233B (zh) | 一种基于结构光的三维物体表面重建方法 | |
CN203084734U (zh) | 用于再生虚拟物体的系统 | |
CN103322911B (zh) | 一种基于图像的特色建筑立面图测绘方法 | |
CN102506711B (zh) | 线激光视觉三维旋转扫描方法 | |
CN105930858A (zh) | 一种带旋转、缩放的快速高精度几何模板匹配方法 | |
CN103838437A (zh) | 基于投影图像的触控定位控制方法 | |
CN107255458A (zh) | 一种垂直投影光栅测量仿真系统及其实现方法 | |
CN102262508B (zh) | 通过鼠标轨迹生成度量工具的方法 | |
CN101509763A (zh) | 单目高精度大型物体三维数字化测量系统及其测量方法 | |
CN106027962A (zh) | 视频监控的覆盖率计算方法及装置、布点方法及系统 | |
CN104123747A (zh) | 多方式触控三维建模方法和系统 | |
CN106949851A (zh) | 一种基于支持向量机的线结构光视觉传感器标定方法 | |
CN103927748B (zh) | 一种基于多矩形图像距离转换模型的坐标标定方法 | |
CN104297516A (zh) | 一种流体表面二维流速场测量方法 | |
CN100385197C (zh) | 便携式无导轨结构光三维扫描测量系统及其测量方法 | |
CN101178811A (zh) | 一种三步法结构光直线光条的图像特征提取方法 | |
CN107036555B (zh) | 一种交叉光轴光栅投影测量仿真系统及其实现方法 | |
CN106839989A (zh) | 一种三维扫描仪的扫描方法及扫描辅助标记板 | |
CN107631700A (zh) | 三维扫描仪与全站仪相结合的三维数据测量方法 | |
Luo et al. | A novel static deformation measurement and visualization method for wind turbine blades using home-made LiDAR and processing program | |
CN102306382A (zh) | 基于道路环境固有特征的事故现场图像校正方法 | |
Song et al. | Flexible line-scan camera calibration method using a coded eight trigrams pattern | |
CN112254681B (zh) | 一种发散型多线激光投影测量仿真系统及其实现方法 |
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 |