本发明公开了一种基于多模态复合编码和极线约束的高效三维图像获取方法，分别采取快速成像模式或高精度成像模式，在快速成像模式下，通过四幅条纹光栅获取两种不同频率相位图，利用极线约束以及左右一致性检验求得高频绝对相位并通过相位与三维坐标映射关系获取三维图像；在高精度测量模式下，通过N+2幅条纹光栅获取两种不同频率的相位，利用极线约束求得低频绝对相位并用低频绝对相位辅助高频相位展开从而获取高频绝对相位并最终通过相位与三维坐标映射关系获取三维图像。本发明不仅保证了成像效率的同时，还提高了成像精度。
1.一种基于多模态复合编码和极线约束的高效三维图像获取方法，其特征在于分别采取快速成像模式或高精度成像模式，在快速成像模式下，通过四幅条纹光栅获取两种不同频率相位图，利用极线约束以及左右一致性检验求得高频绝对相位并通过相位与三维坐标映射关系获取三维图像，其中快速成像方法的具体步骤如下：步骤一，成像系统标定；
步骤三，分析左相机、右相机采集到的光栅条纹分别得到一组高频相位和一组低频相位；
步骤四，利用极线约束搜索左相机高频相位图上每一点，即原始点在空间上的对应点，并通过深度约束剔除部分错误点；
步骤五，剩余的空间对应点被投射到右相机的高频以及低频相位图上，通过原始点和对应点之间相位差异确定最终对应点，并获取高频绝对相位；
步骤六，根据绝对相位获取三维图像，从而实现了对动态场景三维图像高效高精度获取；
对于高频相位图φH上任意一点p，即为原始点，其周期级次k有FH种可能性，这意味着其绝对相位ΦH＝φH+2kπ,k∈[0,FH-1]
可知p在三维空间中对应着FH个点，即FH个对应点，其中MZ,NZ以及CZ由步骤一中所获取的标定参数推导而得，W为投影仪的横向分辨率，xp为p在投影仪上的对应点；三维成像的关键就是要确认这FH个对应点中唯一正确的对应点，首先根据左相机、右相机和投影仪有效范围预先设定为(Zmin,Zmax)，其中Zmin，Zmax表示不同的实际空间深度值，并且Zmin＜Zmax；所有使得Z超出这个预先设定的范围的k及其对应的空间点将被确认为错误点并被排除掉，经过该步骤处理后，p的正确周期级次k及其空间对应点的范围将缩小至FH'种，其中FH'＜＜FH；
首先，将步骤四中在(Zmin,Zmax)内剩余的空间点投影到右相机的成像面上，得到FH'个位于右相机成像面上的二维对应点；
其次，通过φdiff＝φH(p)-φH'(p')进一步排除二维对应点中φdiff超过阈值0.5rad的点，其中φdiff表示原始点与对应点在包裹相位φH与φH'上的差异，此时将正确的对应点的范围从FH'缩小至FH”；
最后，在剩余的FH”个对应点中选取出使得φL(p)-φL'(p')最小的对应点，并认为该点为正确的对应点，那么该点所对应的周期级次k即为正确的周期级次，至此确认原始点唯一的ΦH；
在高精度测量模式下，通过N+2幅条纹光栅获取两种不同频率的相位，其中N≥3，利用极线约束求得低频绝对相位并用低频绝对相位辅助高频相位展开从而获取高频绝对相位并最终通过相位与三维坐标映射关系获取三维图像，高精度成像模式的步骤如下：步骤一，成像系统标定；
步骤二，生成、投射和采集N+2幅双频光栅条纹图案；
步骤三，分析左相机、右相机采集到的光栅条纹分别得到一组高频相位和一组低频相位；
步骤四，利用极线约束搜索左相机低频相位图上每一点，即原始点在空间上的对应点，并通过深度约束剔除部分错误点；
步骤五，剩余的空间对应点被投射到右相机的低频相位图上，通过原始点和对应点之间相位差异确定最终对应点，并通过低频绝对相位获取高频绝对相位；
步骤六，根据绝对相位获取三维图像，从而实现了对静态场景三维图像高效高精度获取；
对于低频相位图φL上任意一点q，即原始点，其周期级次k有NL种可能性，这意味着其绝对相位ΨL＝ψL+2lπ,l∈[0,NL-1]
可知q在三维空间中对应着NL个点，即对应点，其中MZ,NZ以及CZ由步骤一中所获取的标定参数推导而得，W为投影仪的横向分辨率，xq为q在投影仪上的对应点；三维成像的关键就是要确认这NL个对应点中唯一正确的对应点，首先根据左相机、右相机和投影仪有效范围可预先设定为(Zmin,Zmax)，其中Zmin，Zmax表示不同的实际空间深度值，并且Zmin＜Zmax，所有使得Zq超出这个预先设定的范围的l及其对应的空间点将被确认为错误点并被排除掉，经过该步骤q的正确周期级次l及其空间对应点的范围将缩小至NL'种，其中NL'＜＜NL；
首先，将步骤四中在(Zmin,Zmax)内剩余的空间点投影到右相机的成像面上，得到NL'个位于右相机成像面上的二维对应点；
其次，由于NL'＜＜NL从而可从NL'个对应点中选取出使得ψL(p)-ψL'(p')最小的对应点，并认为该点为正确的对应点，那么该点所对应的周期级次l即为正确的周期级次，至此可以确认原始点唯一的ΨL；
2.根据权利要求1所述的基于多模态复合编码和极线约束的高效三维图像获取方法，其特征在于快速成像模式步骤一中的成像系统包括计算机、左相机、右相机和投影仪，其中左相机、右相机与投影仪各自通过数据线与计算机连接，投影仪与左相机、右相机之间分别通过触发线连接，对成像系统进行标定，得到左相机、右相机与投影仪在一个世界坐标系下的标定参数，这些参数包含像素坐标系与世界坐标系之间的缩放参数、平移参数、旋转参数以及畸变参数。
3.根据权利要求1所述的基于多模态复合编码和极线约束的高效三维图像获取方法，其特征在于快速成像模式步骤二中，用计算机通过Matlab生成的四幅双频光栅条纹为两幅正弦条纹光栅加两幅三角波条纹光栅，条纹光栅如下：I1(x,y)＝A(x,y)+B(x,y)sin[πFH(2x/W-1)]I2(x,y)＝A(x,y)+B(x,y)cos[πFH(2x/W-1)]I3(x,y)＝A(x,y)+B(x,y)tri[(2FLx/W-1)]I4(x,y)＝A(x,y)-B(x,y)tri[(2FLx/W-1)]其中Ii(x,y)表示生成图像的像素坐标(x,y)处的光栅条纹强度，i＝1,2,3,4，表示第几幅光栅条纹图像，A为图像直流分量，B为振幅，tri为阈值区间为[-1,1]的三角波函数，FH,FL分别为I1,I2与I3,I4所包含的条纹周期数，W为整幅光栅条纹图像的像素宽度，A＝B＝
127.5，FH,FL分别取值为64和9，x取值范围是0至W-1；这些光栅经过投影仪投射后由左相机、右相机同步采集得到，其中左相机采集到的光栅条纹如下：其中， 为左相机实际拍到的光栅条纹图像，i＝1,2,3,4，(xc,yc)为摄像机拍摄图像的像素坐标，α为被测物表面反射率，β1为被反射的环境光，β2为直接进入摄像机的环境光，ΦH(xc,yc)为光栅条纹图 与光栅条纹图 所包含的相位，ΦL(xc,yc)为光栅条纹图与光栅条纹图 所包含的相位；假设Ac＝α(xc,yc)A(xc,yc)+α(xc,yc)β1(xc,yc)+β2(xc,yc)，Bc＝α(xc,yc)B(xc,yc)，并省去(xc,yc)则上面四式可化简为：右相机的采集过程与上述左相机的过程相同。
4.根据权利要求1所述的基于多模态复合编码和极线约束的高效三维图像获取方法，其特征在于快速成像模式步骤三中，分析左相机、右相机采集到的光栅条纹分别得到一组高频相位和一组低频相位，其中左相机的相位具体如下：其中φH为高频包裹相位，φL为低频包裹相位，利用与左相机相同的步骤可以求得右相机采集的条纹光栅对应的相位分别为φH',φL'。
5.根据权利要求1所述的基于多模态复合编码和极线约束的高效三维图像获取方法，其特征在于快速成像模式步骤六中的具体过程为：根据步骤五中获得的绝对相位ΦH结合下式求得三维图像坐标：
其中EX,FX,EY,FY由步骤一中的标定参数得到，Xp,Yp,Zp为被测物体空间三维坐标，至此可以获取被测物体三维图像。
6.根据权利要求1所述的基于多模态复合编码和极线约束的高效三维图像获取方法，其特征在于高精度成像模式步骤二中，生成、投射和采集N+2幅双频光栅条纹的具体过程如下：N+2幅条纹光栅为两幅低频正弦条纹光栅加N幅高频正弦条纹光栅，其中N≥3，以N＝3为例，条纹光栅如下:J1(x,y)＝A(x,y)+B(x,y)cos[πNH(2x/W-1)]J2(x,y)＝A(x,y)+B(x,y)cos[πNH(2x/W-1)+2π/3]J3(x,y)＝A(x,y)+B(x,y)cos[πNH(2x/W-1)+4π/3]J4(x,y)＝A(x,y)+B(x,y)sin[πNL(2x/W-1)]J5(x,y)＝A(x,y)+B(x,y)cos[πNL(2x/W-1)]其中Ji(x,y)表示生成图像的像素坐标(x,y)处的光栅条纹强度，i＝1,2,3,4，表示第几幅光栅条纹图像，A为图像直流分量，B为振幅，NH,NL分别为J1～J3与J4,J5所包含的条纹周期数，W为整幅光栅条纹图像的像素宽度，A＝B＝127.5，NH,NL分别取值为128和8，x取值范围是
0至W-1；这些光栅经过投影仪投射后由左相机、右相机同步采集得到，其中左相机采集到的光栅条纹如下：其中， 为左相机实际拍到的光栅条纹图像，i＝1,2,3,4，(xc,yc)为摄像机拍摄图像的像素坐标，α为被测物表面反射率，β1为被反射的环境光，β2为直接进入摄像机的环境光，ΨH(xc,yc)为光栅条纹图 到光栅条纹图 所包含的相位，ΨL(xc,yc)为光栅条纹图与光栅条纹图 所包含的相位；假设Ac＝α(xc,yc)A(xc,yc)+α(xc,yc)β1(xc,yc)+β2(xc,yc)，Bc＝α(xc,yc)B(xc,yc)，并省去(xc,yc)则上面四式可化简为：右相机的采集过程与上述左相机的过程相同。
7.根据权利要求1所述的基于多模态复合编码和极线约束的高效三维图像获取方法，其特征在于高精度成像模式步骤三中，分析左相机、右相机采集到的光栅条纹分别得到一组高频相位和一组低频相位，其中左相机的相位具体如下：其中ψH为高频包裹相位，ψL为低频包裹相位，利用与左相机相同的步骤可以求得右相机采集的条纹光栅对应的相位分别为ψH',ψL'。
8.根据权利要求1所述的基于多模态复合编码和极线约束的高效三维图像获取方法，其特征在于高精度成像模式步骤六中，根据绝对相位重构三维图像的具体过程为：根据步骤五中获得的绝对相位ΨH结合下式求得三维图像坐标：
其中EX,FX,EY,FY由步骤一中的标定参数得到，Xq,Yq,Zq为被测物体空间三维坐标，至此可以获取被测物体三维图像。
一种基于多模态复合编码和极线约束的高效三维图像获取
方法
技术领域
[0001] 本发明属于三维成像技术领域，特别是一种基于多模态复合编码和极线约束的高效三维图像获取方法。
背景技术
[0002] 在三维成像领域，快速获取目标物体的高精度三维数据一直以来都是一个重要的技术难点。早期人们通过机械的三坐标测量机来对目标物体逐点探测以获取三维图像，但是这种逐点接触式的图像获取技术一方面效率极其低下，另一方面对被测物体有所损坏，这些缺点使得该技术在诸如人体检测，文物保护等领域难以得到应用。相比于传统的机械三维图像获取技术，光学三维图像获取技术由于其非接触，高效率的优点在科学研究，工业检测等领域得到了广泛应用。近些年来由于数字投影设备的发展，光学三维成像方法中能够实现全场成像的条纹投影技术更是成为了研究热点(S.S.Gorthi and P.Rastogi,“Fringe projection techniques:whither we are？”Opt.Eng.48,133–140(2010).)。目前条纹投影领域中主流的两种技术分别是傅里叶轮廓术(M.Takeda and K.Mutoh“, Fourier transform profilometry for the automatic measurement of 3-d object shapes,”Applied optics 22,3977-3982(1983).)与相移轮廓术(V.Srinivasan,H.-C.Liu,and M.Halioua“, Automated phase-measuring profilometry of 3-d diffuse objects,”Applied optics 23,3105-3108(1984).)。
[0003] 相比于傅里叶轮廓术，相移轮廓术由于对环境光和噪声不敏感以及计算简单更适合自动化高精度三维图像的获取。该技术主要通过向被测物体投影多幅(至少三幅)相移光栅条纹并通过相机同步采集受到被测物体调制后的光栅条纹来获取被测物体的相位，最终通过相位到高度的映射关系来获取被测物体的三维图像。如何保证在较高成像效率下得到较高的成像精度一直以来都是相移轮廓术领域的一个重点与难点。一般而言，在相移轮廓术中，投影的光栅条纹越多越有利于相位精确获取，从而能得到越精确的被测物体三维图像。但是过多的光栅条纹会极大的影响相位获取的效率从而影响三维图像获取的速度(效率)(陈钱；冯世杰；顾国华；左超；孙佳嵩；喻士领；申国辰；李如斌.基于双频三灰阶正弦光栅条纹投影的时间相位去包裹方法:中国,201410027275.4[I].2013-04-30.)。针对动态(准静态)物体三维成像，为了保证成像的正确性，必须要求较高的成像效率。另一方面，尽管在测量静态物体时成像效率对最终成像正确性没有影响，但是面对流水线式的批量物体三维数据收集时，高的成像效率依旧有着绝对优势(龙佳乐；张建民；范智晖.一种基于三波长条纹投影的快速三维测量系统:中国,201620177719.7[I].2016-09-07.)。
[0004] 基于相移轮廓术的三维成像技术相比于其他方法在测量方式，成像效率与成像精度上有着很大优势并在诸如文物保护、人体检测等领域有着广泛应用，但是要想获得更多的应用必须进一步提升其成像效率与成像精度，克服二者之间的矛盾。
发明内容
[0005] 本发明的目的在于提供一种基于多模态复合编码和极线约束的高效三维图像获取方法，通过复合编码和极线约束提高了三维成像精度和效率。
[0006] 实现本发明目的的技术解决方案为：一种基于多模态复合编码和极线约束的高效三维图像获取方法，分别采取快速成像模式或高精度成像模式，在快速成像模式下，通过四幅条纹光栅获取两种不同频率相位图，利用极线约束以及左右一致性检验求得高频绝对相位并通过相位与三维坐标映射关系获取三维图像；在高精度测量模式下，通过N+2幅条纹光栅获取两种不同频率的相位，利用极线约束求得低频绝对相位并用低频绝对相位辅助高频相位展开从而获取高频绝对相位并最终通过相位与三维坐标映射关系获取三维图像。
[0007] 本发明与现有技术相比，其显著优点：(1)本发明的快速成像模式利用四幅光栅条纹获取被测物体三维图像，与现有技术相比，四幅复合编码光栅条纹组合保证了三维图像测量的高效性；于此同时，双目视觉中极线约束的引入使得该技术能用四幅复合光栅条纹获取高达64个条纹周期的高精度相位，保证了三维图像的高精度。最后，由于高频绝对相位直接通过极线约束求解，并没有借助于低频绝对相位，这避免了两种频率条纹光栅调制度差异带来的三维图像的不正确性。如图2所示，在均采用四幅条纹的前提下本发明的快速成像模式成像结果更精细即拥有更高的成像精度。因此本发明的快速成像模式在保证了成像效率的同时，提高了成像精度。(2)本发明的高精度成像模式通过N+2幅光栅条纹获取了被测物体三维图像，与现有技术相比，极线约束的引入使得8周期的低频绝对相位能够通过两幅低频光栅条纹直接获取，这极大的避免了现有技术中利用多帧光栅条纹(通常远大于2)获取低频绝对相位的冗余性，提升了测量效率。另一方面，N幅128条纹周期数的编码图案保证了最终获取的三维图像精度不低于现有的技术。如图4所示，本发明仅利用五幅条纹即可实现50um量级精度的三维成像。因此本发明的高精度成像模式在保证成像精度的同时，极高地提升了成像效率。
[0008] 下面结合附图对本发明作进一步详细描述。
附图说明
[0009] 图1为快速成像模式步骤流程示意图。
[0010] 图2为快速成像模式测试结果，图2(a)与图2(b)分别为本方法与现有技术获取的动态物体的三维图像，图2(c)与图2(f)为图2(a)局部细节，图2(d)与图2(g)为图2(b)局部细节，图2(e)与图2(h)为局部细节的对比。
[0011] 图3为高精度成像模式步骤流程示意图。
[0012] 图4为高精度成像模式测试结果，图4(a)为标准陶瓷平板三维成像结果，图4(b)为图4(a)中某一行的细节展示。
具体实施方式
[0013] 结合图1和图3，本发明基于多模态复合编码和极线约束的高效三维图像获取方法，分别采取快速成像模式或高精度成像模式，这两种成像模式共用相同的成像系统，以及相同的成像系统标定结果，并且都采用极线约束结合辅助相位信息作为相位展开的核心技术，其相比于传统多频时间相位展开方法在效率上大大提高，相比于传统的空间相位展开算法在效率和准确度上大大提高。所述的快速成像模式步骤如下：
[0014] 步骤一，成像系统标定；
[0015] 步骤二，生成、投射和采集四幅双频光栅条纹；
[0016] 步骤三，分析左相机、右相机采集到的光栅条纹分别得到一组高频相位和一组低频相位；
[0017] 步骤四，利用极限约束搜索左相机高频相位图上每一点，即原始点在空间上的对应点，并通过深度约束剔除部分错误点；
[0018] 步骤五，剩余的空间对应点被投射到右相机的高频以及低频相位图上，通过原始点和对应点之间相位差异确定最终对应点，并获取高频绝对相位；
[0019] 步骤六，根据绝对相位获取三维图像，从而实现了对动态场景三维图像高效高精度获取。
[0020] 所述的高精度成像模式步骤如下：
[0021] 步骤一，成像系统标定；
[0022] 步骤二，生成、投射和采集N+2幅双频光栅条纹图案；
[0023] 步骤三，分析左相机、右相机采集到的光栅条纹分别得到一组高频相位和一组低频相位；
[0024] 步骤四，利用极限约束搜索左相机低频相位图上每一点，即原始点在空间上的对应点，并通过深度约束剔除部分错误点；
[0025] 步骤五，剩余的空间对应点被投射到右相机的低频相位图上，通过原始点和对应点之间相位差异确定最终对应点，并通过低频绝对相位获取高频绝对相位；
[0026] 步骤六，根据绝对相位获取三维图像，从而实现了对静态场景三维图像高效高精度获取。
[0027] 下面详细说明两种成像模式的过程。
[0028] 本发明快速成像模式流程示意图如图3所示，具体实现步骤如下：
[0029] 步骤一，成像系统标定。
[0030] 成像系统包括计算机、左相机、右相机和投影仪，其中左相机、右相机与投影仪各自通过数据线与计算机连接，投影仪与左相机、右相机之间分别通过触发线连接。成像系统搭建完成后，用(Z.Zhang“,Aflexible new technique for camera calibration.”IEEE Transactions on pattern analysis and machine intelligence.22(11),1330-1334(2000).)中的标定方法进行系统标定，得到左相机、右相机与投影仪在一个世界坐标系下的标定参数，这些参数包含像素坐标系与世界坐标系之间的缩放参数、平移参数、旋转参数以及畸变参数。
[0031] 步骤二，生成、投射和采集四幅双频光栅条纹。
[0032] 计算机通过Matlab软件生成的四幅双频光栅条纹为两幅正弦条纹光栅加两幅三角波条纹光栅，条纹光栅如下：
[0033] I1(x,y)＝A(x,y)+B(x,y)sin[πFH(2x/W-1)]
[0034] I2(x,y)＝A(x,y)+B(x,y)cos[πFH(2x/W-1)]
[0035] I3(x,y)＝A(x,y)+B(x,y)tri[(2FLx/W-1)]
[0036] I4(x,y)＝A(x,y)-B(x,y)tri[(2FLx/W-1)]
[0037] 其中Ii(x,y)表示生成图像的像素坐标(x,y)处的光栅条纹强度，i＝1,2,3,4，表示第i幅光栅条纹图像，A为图像直流分量，B为振幅，tri为阈值区间为[-1,1]的三角波函数，FH,FL分别为I1,I2与I3,I4所包含的条纹周期数，W为整幅光栅条纹图像的像素宽度，A＝B＝127.5，FH,FL分别取值为64和9，x取值范围是0至W-1。这些光栅经过投影仪投出后由左相机、右相机同步采集得到。为了简单起见，这里只对左相机作具体分析，右相机的分析过程同左相机，其中左相机采集到的光栅条纹如下：
[0038]
[0039]
[0040]
[0041]
[0042] 其中， 为左相机实际拍到的光栅条纹图像，i＝1,2,3,4，(xc,yc)为摄像机拍摄图像的像素坐标，α为被测物表面反射率，β1为被反射的环境光，β2为直接进入摄像机的c c c c环境光，ΦH(x ,y)为光栅条纹图 与光栅条纹图 所包含的相位，ΦL(x ,y)为光栅条纹图 与光栅条纹图 所包含的相位。假设Ac＝α(xc,yc)A(xc,yc)+α(xc,yc)β1(xc,yc)+β2(xc,c c c c c c c c
y)，B＝α(x ,y)B(x ,y)，并省去(x ,y)则上面四式可化简为：
[0043]
[0044]
[0045]
[0046]
[0047] 右相机的采集过程与上述左相机的过程相同。
[0048] 步骤三，分析左相机、右相机采集到的光栅条纹分别得到一组高频相位和一组低频相位。
[0049] 根据步骤二中左相机所采集到的图像可以得到两组如下所示相位：
[0050]
[0051]
[0052] 其中φH为高频包裹相位，φL为低频包裹相位，利用与左相机相同的步骤可以求得右相机采集的条纹光栅对应的相位分别为φH',φL'。
[0053] 步骤四，利用极限约束搜索左相机高频相位图上每一点(原始点)在空间上的对应点，并通过深度约束剔除部分错误点。
[0054] 对于φH上任意一点p(原始点)，其周期级次k有FH种可能性，这意味着其绝对相位[0055] ΦH＝φH+2kπ,k∈[0,FH-1]
[0056] 有FH个不同值。再由下式
[0057]
[0058]
[0059] 可知p在三维空间中对应着FH个点，即FH个对应点，其中MZ,NZ以及CZ由步骤一中所获取的标定参数推导而得(K.Liu,Y.Wang,et al“Dual-frequency pattern scheme for high-speed 3-D shape measurement.”Optics express.18(5),5229-5244(2010).)。三维成像的关键就是要确认这FH个对应点中唯一正确的对应点。考虑到实际三维成像系统中测量空间范围有限，根据左相机、右相机和投影仪有效范围预先设定为(Zmin,Zmax)，如Zmin＝-200mm,Zmax＝200mm，所有使得Z超出这个预先设定的范围的k及其对应的空间点将被确认为错误点并被排除掉。经过此步骤，p的正确周期级次k及其空间对应点的范围将缩小至FH'种，其中FH'＜＜FH。
[0060] 步骤五，剩余的空间对应点被投射到右相机的高频以及低频相位图上，通过原始点和对应点之间相位差异确定最终对应点，并获取高频绝对相位。
[0061] 首先，将步骤四中在(Zmin,Zmax)内剩余的空间点投影到右相机的成像面上，得到FH'个位于右相机成像面上的二维对应点。
[0062] 实际上正确的对应点p'的包裹相位φH',φL'与原始点之间的包裹相位φH,φL应非常接近，从而可依据下式
[0063] φdiff＝φH(p)-φH'(p')
[0064] 进一步排除二维对应点中φdiff超过阈值0.5rad(rad为弧度单位，阈值大小先验确定)的点，其中φdiff表示原始点与对应点在包裹相位φH与φH'上的差异。此时可以将正确的对应点的范围从FH'缩小至FH”。
[0065] 然后，在剩余的FH”个对应点中选取出使得φL(p)-φL'(p')最小的对应点，并认为该点为正确的对应点，那么该点所对应的周期级次k即为正确的周期级次，至此可以确认原始点唯一的ΦH。
[0066] 步骤六，根据绝对相位获取三维图像。
[0067] 根据步骤五中获得的绝对相位ΦH结合下式求得三维图像坐标。
[0068]
[0069]
[0070] Xp＝EXZp+FX
[0071] Yp＝EYZp+FY
[0072] 其中EX,FX,EY,FY由步骤一中的标定参数得到(K.Liu,Y.Wang,et al “Dual-frequency pattern scheme for high-speed 3-D shape measurement.”Optics express.18(5),5229-5244(2010).)，Xp,Yp,Zp为被测物体空间三维坐标，至此可以获取被测物体三维图像。
[0073] 通过上述步骤可以看出，本发明在快速成像模式下，通过四幅光栅条纹获取被测物体三维图像，与现有技术相比，四幅复合编码光栅条纹组合保证了三维图像测量的高效性；于此同时，双目视觉中极线约束的引入使得该技术能用四幅复合光栅条纹获取高达64个条纹周期的高精度相位，保证了三维图像的高精度；最后，由于高频绝对相位直接通过极线约束求解，并没有借助于低频绝对相位，这避免了两种频率条纹光栅调制度差异带来的三维图像的不正确性。
[0074] 为了测试一种基于多模态复合编码和极线约束的高效三维图像获取技术效果，本发明给出了两组试验结果。图2(a)与图2(b)分别为本技术与现有技术获取的动态物体的三维图像，图2(c)与图2(f)为图2(a)局部细节，图2(d)与图2(g)为图2(b)局部细节，图2(e)与图2(h)为局部细节的对比。由这些实验结果可以看出本发明的快速成像模式与现有技术相比保留了更多的细节测量结果更加精细，在保证高效率三维成像的同时在成像精度上有了很大提升。
[0075] 本发明高精度成像模式流程示意图如图3所示，具体实现步骤如下：
[0076] 步骤一，成像系统标定。
[0077] 成像系统(包括计算机，左右两个相机、一个投影仪，其中左相机、右相机与投影仪各自通过数据线与计算机连接，投影仪与相机之间通过触发线连接)搭建完成后，用(Z.Zhang,“A flexible new technique for camera calibration.”IEEE Transactions on pattern analysis and machine intelligence.22(11),1330-1334(2000).)中的标定方法进行系统标定，得到两个相机与一个投影仪在一个世界坐标系下的标定参数。
[0078] 步骤二，生成、投射和采集四幅双频光栅条纹。
[0079] N+2幅条纹光栅为两幅低频正弦条纹光栅加N幅高频正弦条纹光栅，其中N≥3，为了简单起见，以N＝3为例，条纹光栅如下:
[0080] J1(x,y)＝A(x,y)+B(x,y)cos[πNH(2x/W-1)]
[0081] J2(x,y)＝A(x,y)+B(x,y)cos[πNH(2x/W-1)+2π/3]
[0082] J3(x,y)＝A(x,y)+B(x,y)cos[πNH(2x/W-1)+4π/3]
[0083] J4(x,y)＝A(x,y)+B(x,y)sin[πNL(2x/W-1)]
[0084] J5(x,y)＝A(x,y)+B(x,y)cos[πNL(2x/W-1)]
[0085] 其中Ji(x,y)表示生成图像的像素坐标(x,y)处的光栅条纹强度，i＝1,2,3,4，表示第几幅光栅条纹图像，A为图像直流分量，B为振幅，NH,NL分别为J1～J3与J4,J5所包含的条纹周期数，W为整幅光栅条纹图像的像素宽度，A＝B＝127.5，NH,NL分别取值为128和8，x取值范围是0至W-1。这些光栅经过投影仪投出后由左相机、右相机同步采集得到。为了简单起见，这里只对左相机作具体分析，右相机的分析过程同左相机，其中左相机采集到的光栅条纹如下：
[0086]
[0087]
[0088]
[0089]
[0090]
[0091] 其中， 为左相机实际拍到的光栅条纹图像，i＝1,2,3,4，(xc,yc)为摄像机拍摄图像的像素坐标，α为被测物表面反射率，β1为被反射的环境光，β2为直接进入摄像机的环境光，ΨH(xc,yc)为光栅条纹图 到光栅条纹图 所包含的相位，ΨL(xc,yc)为光栅条纹图 与光栅条纹图 所包含的相位。假设Ac＝α(xc,yc)A(xc,yc)+α(xc,yc)β1(xc,yc)+β2(xc,yc)，Bc＝α(xc,yc)B(xc,yc)，并省去(xc,yc)则上面四式可化简为：
[0092]
[0093]
[0094]
[0095]
[0096]
[0097] 右相机的采集过程与上述左相机的过程相同。
[0098] 步骤三，分析左相机、右相机采集到的光栅条纹分别得到一组高频相位和一组低频相位。其中左相机的相位具体如下：
[0099]
[0100]
[0101] 其中ψH为高频包裹相位，ψL为低频包裹相位，利用与左相机相同的步骤可以求得右相机采集的条纹光栅对应的相位分别为ψH',ψL'。
[0102] 步骤四，利用极限约束搜索左相机高频相位图上每一点(原始点)在空间上的对应点，并通过深度约束剔除部分错误点。
[0103] 对于φL上任意一点q(原始点)，其周期级次k有NL种可能性，这意味着其绝对相位[0104] ΨL＝ψL+2lπ,l∈[0,NL-1]
[0105] 有NL个不同值，再由下式
[0106]
[0107]
[0108] 可知q在三维空间中对应着NL个点(对应点)，其中MZ,NZ以及CZ由步骤一中所获取的标定参数推导而得(K.Liu,Y.Wang,et al“Dual-frequency pattern scheme for high-speed 3-D shape measurement.”Optics express.18(5),5229-5244(2010).)。三维成像的关键就是要确认这NL个对应点中唯一正确的对应点。考虑到实际三维成像系统中测量空间范围有限，根据左相机、右相机和投影仪有效范围可预先设定为(Zmin,Zmax)，如Zmin＝-200mm,Zmax＝200mm，所有使得Zq超出这个预先设定的范围的l及其对应的空间点将被确认为错误点并被排除掉。经过此步骤q的正确周期级次l及其空间对应点的范围将缩小至NL'种，其中NL'＜＜NL。
[0109] 步骤五，剩余的空间对应点被投射到右相机的低频相位图上，通过原始点和对应点之间相位差异确定最终对应点，并通过低频绝对相位获取高频绝对相位。
[0110] 首先，将步骤四中在(Zmin,Zmax)内剩余的空间点投影到右相机的成像面上，得到NL'个位于右相机成像面上的二维对应点。
[0111] 实际上正确的对应点q'的包裹相位ψL'与原始点之间的包裹相位ψL应非常接近，由于NL'＜＜NL从而可从NL'个对应点中选取出使得ψL(p)-ψL'(p')最小的对应点，并认为该点为正确的对应点，那么该点所对应的周期级次l即为正确的周期级次，至此可以确认原始点唯一的ΨL。
[0112] 然后，通过下式
[0113]
[0114] ψH＝ψH+2lπ
[0115] 求得最终的高频绝对相位ΨH。
[0116] 步骤六，根据绝对相位获取三维图像。
[0117] 根据步骤五中获得的绝对相位ΨH结合下式求得三维图像坐标。
[0118]
[0119]
[0120] Xq＝EXZq+FX
[0121] Yq＝EYZq+FY
[0122] 其中EX,FX,EY,FY由步骤一中的标定参数得到(K.Liu,Y.Wang,et al “Dual-frequency pattern scheme for high-speed 3-D shape measurement.”Optics express.18(5),5229-5244(2010).)，Xq,Yq,Zq为被测物体空间三维坐标，至此可以获取被测物体三维图像。
[0123] 通过上述步骤可以看出，本方法在高精度模式下，通过N+2幅光栅条纹获取了被测物体三维图像，与现有技术相比，极线约束的引入使得8周期的低频绝对相位能够通过两幅低频光栅条纹直接获取，这极大的避免了现有技术中利用多帧光栅条纹(通常远大于2)获取低频绝对相位的冗余性，提升了测量效率；另一方面，N幅128条纹周期数的编码图案保证了最终获取的三维图像精度不低于现有的技术。
[0124] 为了测试一种基于多模态复合编码和极线约束的高效三维图像获取方法的效果，我们给出了一组试验结果。图4(a)为静态物体在高精度成像模式下的三维成像结果，图4(b)为图4(a)中某一行的细节展示，由图4可知对静态物体本技术高精度成像模式能够实现高达50um的成像精度。
