CN103940370B - 基于周期互质混合编码的目标物体三维信息获取方法 - Google Patents
基于周期互质混合编码的目标物体三维信息获取方法 Download PDFInfo
- Publication number
- CN103940370B CN103940370B CN201410178890.5A CN201410178890A CN103940370B CN 103940370 B CN103940370 B CN 103940370B CN 201410178890 A CN201410178890 A CN 201410178890A CN 103940370 B CN103940370 B CN 103940370B
- Authority
- CN
- China
- Prior art keywords
- target object
- phase
- image
- stripe
- template
- 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
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000003287 optical effect Effects 0.000 claims abstract description 14
- 230000010363 phase shift Effects 0.000 claims description 17
- 230000000737 periodic effect Effects 0.000 claims description 14
- 238000003384 imaging method Methods 0.000 claims description 11
- 238000012012 milestone trend analyses Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 19
- 230000003993 interaction Effects 0.000 abstract description 2
- 230000000903 blocking effect Effects 0.000 abstract 1
- 238000005070 sampling Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000002366 time-of-flight method Methods 0.000 description 2
- 238000010146 3D printing Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于周期互质混合编码的目标物体三维信息获取方法,主要解决现有相移方法在截断相位延拓时存在的二义性问题。其实现步骤为:1)设计三幅模板图像T1、T2、T3;2)将投影仪P与摄像机V水平放置,并使两者光轴平行,用投影仪P将模板图像T1、T2、T3依次投射到目标物体上;3)通过摄像机V同步拍摄三幅模板投影到目标物体后发生的对应变形图像 并传回计算机对其进行解码,得到其正弦条纹截断相位φps和方波条纹截断相位φsq;4)根据正弦条纹截断相位φps和方波条纹截断相位φsq,求解得到目标物体的三维信息;本发明具有抗干扰能力强、目标物体三维信息测量精度高、分辨率高的优点,可用于逆向工程、人机互动、文物三维重建等领域。
Description
技术领域
本发明属于计算机视觉测量技术领域,主要涉及一种目标物体三维信息的获取方法,可用于人机交互、3D打印、逆向工程、文物三维重建。
背景技术
随着当今计算机视觉技术的发展和人们生活水平的提高,传统的二维信息已经不能满足人们的需求,快速、高效的从场景中获取三维信息成为当前研究的热点。从场景中获取三维信息的过程叫做深度获取。当前深度获取方法主要分为接触式和非接触式两大类。
接触式是通过实际触碰物体表面来得到三维信息,如坐标测量机。接触式测量法虽可以获得高精度数据,但其硬件要求高、速率低、使用不方便,对物体存在一定损害性,且往往只能得到少数特征点的深度信息。
非接触式则不需要与目标物体接触就能获取目标物体的三维信息,如现有的激光雷达、CT、声呐等方法。非接触式的测量速度快,采样密集,已成为现阶段获取三维信息一个重要手段。在非接触式三维信息获取方法里面又主要分为被动方法和主动方法两大类。
被动式三维信息获取方法是通过测量由目标物体表面反射的周遭辐射性,如激光、可见光,来计算目标物体的三维信息,常见的计算方法有立体视觉法、阴影测量法、聚焦法等方法。其中应用最为广泛的立体视觉法是通过模拟人的视觉方式,采用两个或多个摄像机获得的多个角度下的被测物体的数字图像,根据特定采样点在两幅图像中的匹配像素点及摄像机的相对空间位置、光学参数,依据三角测量原理计算得出该点的三维坐标。虽然被动式测量方法对成像设备要求不高,不需要额外光源,且操作简便,容易实现,但是对于不存在明显特征的图像,如无明显角点,该方法计算量大,匹配精度低,难以获得精细的测量结果。
主动式三维信息获取是通过将额外的能量投射到目标物体,借助能量的反射来计算目标物体的三维信息,常见的投射能量有可见光、高能量光束、超声波和X射线等,主动式三维信息获取方法主要有激光扫描法、飞行时间法和结构光法。
激光扫描法采用激光扫描仪对物体逐点进行扫描从而获取深度信息,虽然精度较高,但是速度很慢;飞行时间法即time of flight,是新兴的三维成像技术之一,它的测量速度快,但采集得到的深度图像分辨率和精度都较低。
结构光法,其原理是利用投射设备将具有一定规律的结构光模板投射到被测目标表面,利用图像采集设备记录被测目标表面的图像,将采集到的图像数据结合三角测距原理以及图形处理技术计算出物体表面的深度数据信息。这种方法既利用了图像作为信息载体,又利用可控光源形成明显纹理特征,能大大减小对诸如缺乏纹理、光滑、无明显灰度的表面区域进行匹配的难度。鉴于结构光三角化技术结构简单、测量速度快、精度高等特点,已得到广泛应用。
结构光法根据投射光图案又可分为单点法、单线法和编码法。单点结构光法具有准确度高、简单可靠的优点,但其测量速度慢;单线结构光法虽然比单点法测量速度有所提高,但是相应确定对应关系的难度也增大,影响了测量准确度。相比之下,编码结构光法通过编码解码使识别简化,有效提高了测量效率,因此成为目前结构光法最有前途的发展方向。
编码结构光法可分为空间、直接和时间三种编码方法。空间编码向景物投射一副编码图案、得到一副对应的编码图像,将编码图像和编码方式对照进行解码,从而解决两者的对应关系,它具有适合动态场景测量的优点,但是由于受其空间领域内特征点的影响,解码困难,测量误差增大,存在分辨率较低、受景物表面发射率不一致及颜色的影响等缺点。直接编码是指每个被编码像素都由其本身的灰度或颜色来识别,由于编码一般被浓缩成唯一的一副图案,噪声值增加,所以所用灰度或颜色的频谱非常宽,必须投射附加的参考图案,以便区分所用的投射灰度或颜色。该方法投射的图案数目较少,理论上适合动态测量且可以达到高分辨率,但编码图案识别困难,降低了测量准确度。
时间编码是将多个不同的编码图案按时序先后投射到景物表面、得到相应的编码图像序列,将编码图像序列组合起来进行解码,它具有易于实现、3D测量准确度高、空间分辨率高等优点。时间编码结构光方法又可分为两灰度级编码方法、多灰度级编码方法和组合编码方法。
两灰度级编码方法中,最具代表性应用最为广泛的是格雷码法.其优点是鲁棒性好,可实现突变表面及不连续表面的测量,测量准确度高,但分辨率相对较低。
多灰度级编码方法中,最具代表性应用最为广泛的是相移法.其优点是测量准确度高,分辨率高,但在对相位进行相位解缠绕时存在二义性。为了解决这个问题,哈尔滨理工大学发明了一种三编码周期灰度梯形相移结构光三维信息获取方法,专利号为2011104314023。此发明设计了一种三编码周期灰度梯形相移结构光三维信息获取方法,它通过选择三个不同的编码周期,对同一个编码周期,投射三幅分别移动1/3周期的编码图案,总计9幅编码图案,利用同一个采样点在不同编码周期内的同余方程确定采样点对应的投影角度,从而确定采样点的三维信息。由于此发明需要额外的投射9幅编码图案,其实现相对复杂。
发明内容
本发明的目的在于针对上述已有技术的不足,提出一种结合相移编码和方波编码的周期互质混合编码方法,以在不增加投影模板数量的情况下,减小相位解缠绕时所存在的二义性,保证对三维信息测量的准确度和高分辨率。
实现本发明目的技术方案是:对相移法中的三幅光栅相位模版进行改进,将改进后的周期互质混合编码模版投影到目标物体上,通过摄像机拍摄目标物体上变形的条纹图案;根据拍摄到的条纹图案求解出每个像素点的截断相位,并进行基于中国剩余定理的匹配,进而得到每个像素点在编码模板中的位置坐标。具体步骤包括如下:
1)设计周期互质混合编码的三维扫描系统所需投影的三幅模板图像T1、T2、T3:
1a)设定三幅结合相移条纹和方波条纹的周期互质混合编码图像模板中正弦条纹的直流分量为A,正弦条纹幅度为B;
1b)设定三幅模板图像T1、T2、T3的初始相移δi分别为δ1=0,δ2=2π/3,δ3=-2π/3;
1c)分别定义三幅模板图像T1、T2、T3中像素点(x,y)处的灰度值:
其中,表示模板图像中像素点(x,y)处方波条纹的灰度值,C为模板图像中方波条纹幅度,mod(a,b)是指a除以b得到的余数;Imix(x,y,δi)为第i幅模板图像中(x,y)处像素点的灰度值,i=1,2,3;Tps和Tsq分别为三幅模板图像T1、T2、T3中正弦条纹周期和方波条纹周期,它们是互为质数的整数;
(2)将投影仪P与摄像机V水平放置,并使两者光轴平行,用投影仪P将模板图像T1、T2、T3依次投射到目标物体上;
(3)通过摄像机V同步拍摄三幅编码模板投影到目标物体后发生的对应变形图像并传回至计算机;
(4)计算机对摄像机捕获到的对应变形图像进行解码,得到这些变形图像中正弦条纹截断相位φps和方波条纹截断相位φsq;
(5)根据上步所解码出的正弦条纹截断相位φps和方波条纹截断相位φsq,进行基于互质原理的匹配,求解出变形图像中像素点(x,y)对应在模板图像中的位置(x0,y0),再利用三角测距原理求解得到目标物体在世界坐标系下的三维信息:
其中,Xw,Yw,Zw分别为目标物体在世界坐标系下沿x轴、y轴、z轴的三维坐标值,世界坐标系的原点选在摄像机的光心,x轴沿摄像机成像平面的水平方向,y轴沿摄像机成像平面的垂直方向,z轴与摄像机成像平面垂直,x、y分别为变形图像中像素点(x,y)的行、列坐标,x0为变形图像中像素点(x,y)对应在模板图像中的位置(x0,y0)的行坐标,f为摄像机的焦距,b为摄像机V光心和投影仪P光心的水平距离,M为变形图像的总行数,N为变形条纹图像总列数,MT为模板图像T1、T2、T3的总行数,hcp为摄像机单个像素所代表的实际高度,其值为θV为摄像机V的垂直视场角,wcp为摄像机单个像素所代表的实际宽度,hpp为投影仪模版中单个像素所代表的实际高度。
本发明与现有技术相比具有以下优点:
第一,本发明对传统相移法中的三幅光栅相位模版进行改进,不需增加额外的物理设备和测量步骤即可提高相位解缠绕的计算速度,实现简单。
第二,本发明所述编码方法结合了相移编码和方波编码,避开了传统相移法在相位解缠绕时所存在的二义性,保证对三维信息测量的准确度和高分辨率。
第三,本发明对变形图像中任意像素点均可实现匹配,且求解每个点的三维信息时互相独立,可实现并行计算。
附图说明
图1为本发明使用的系统框图;
图2为本发明的实现流程图;
图3为本发明设计的三幅周期互质混合编码模板;
图4为本发明获得的正弦条纹截断相位图和方波条纹截断相位图;
图5为本发明获得的目标物体的三维重构图。
具体实施方式
本发明是对传统相移法的改进,不需要增加额外的测量设备和测量步骤。以下参照附图对本发明作进一步详细描述。
参照图1,本发明使用的三维扫描系统包括:待扫描的目标物体、投影仪P和摄像机V。投影仪P投射编码模板到目标物体上,摄像机V同步拍摄编码模板投影到目标物体后发生的对应变形图像。
参照图2,本发明的实现步骤如下:
步骤1、设计周期互质混合编码的三维扫描系统所需投影的三幅模板图像T1、T2、T3,具体步骤如下:
1a)设定三幅结合相移条纹和方波条纹的周期互质混合编码图像模板中正弦条纹的直流分量为A,正弦条纹幅度为B;
1b)设定三幅模板图像T1、T2、T3的初始相移δi分别为δ1=0,δ2=2π/3,δ3=-2π/3;
1c)分别定义三幅模板图像T1、T2、T3中像素点(x,y)处的灰度值:
其中,表示模板图像中像素点(x,y)处方波条纹的灰度值,C为模板图像中方波条纹幅度,mod(a,b)是指a除以b得到的余数;Imix(x,y,δi)为第i幅模板中(x,y)处像素点的灰度值,i=1,2,3;Tps和Tsq分别为三幅模板图像T1、T2、T3中正弦条纹周期和方波条纹周期,它们是互为质数的整数。
按上述步骤设计得到的模板图像T1、T2、T3,如图3所示,其中图3(a)为第一幅模板图像T1、图3(b)为第二幅模板图像T2、图3(c)为第三幅模板图像T3。
步骤2、将投影仪P与摄像机V水平放置,并使两者光轴平行,用投影仪P将所述的三幅模板图像T1、T2、T3依次投射到目标物体上。
步骤3、通过摄像机V同步拍摄三幅模板投影到目标物体后发生的对应变形图像并传回至计算机;
步骤4、计算机对摄像机捕获到的对应变形图像进行解码,得到这些变形图像中正弦条纹上像素点(x,y)处截断相位φps(x,y)和方波条纹上像素点(x,y)处截断相位φsq(x,y),具体步骤如下:
4a)将三幅模板图像T1、T2、T3中像素点(x,y)处正弦条纹的灰度值表示为:
其中,Ips(x,y,δi)为第i幅模板中像素点(x,y)处正弦条纹的灰度值,三幅模板图像T1、T2、T3的初始相移δi,i=1,2,3,取δ1=0,δ2=2π/3,δ3=-2π/3,A和B分别为模板图像中正弦条纹的直流分量和振幅,Tps是模板图像中正弦条纹的周期,单位为像素;
4b)将摄像机捕获到的对应变形图像上像素点(x,y)处正弦条纹灰度值表示为:
其中,和分别为背景光强和正弦条纹调制深度,Φps(x,y)为三幅模板图像T1、T2、T3中正弦条纹被目标物体调制后的相位,即正弦条纹截断相位φps(x,y)的延拓相位;
对进行正弦函数和差化积展开,得到如下式子:
式3)可以看成是的Fourier级数展开,其中为直流分量,和为一阶谐波分量,根据三角函数的正交性有:
从式4)可以得到:
4c)将三幅模板图像T1、T2、T3中像素点(x,y)处方波条纹的灰度值表示为:
其中,C为模板图像中方波条纹幅度,Tsq为方波条纹周期,mod(x,Tsq)表示x对方波条纹周期Tsq求余,即x除以Tsq的余数;
由于模板图像T1、T2、T3中方波条纹是均匀变化的,所以方波条纹的延拓相位Φsq(x,y)是随着x变化的线性函数,它们的关系如下:
其中,φsq(x,y)为模板图像T1、T2、T3中像素点(x,y)处方波条纹的截断相位,n为模板图像中方波条纹的延拓数;
4d)根据步骤4a)和步骤4c),得到三幅模板图像T1、T2、T3中像素点(x,y)处的灰度值:
其中,Imix(x,y,δi)为第i幅模板图像中像素点(x,y)处的灰度值,δi为第i幅模板图像的初始相移,i=1,2,3,取δ1=0,δ2=2π/3,δ3=-2π/3,Tps和Tsq分别为正弦条纹周期和方波条纹周期,它们是互为质数的整数;
4e)记摄像机捕获到的对应变形图像上像素点(x,y)处的灰度值为:
其中,i=1,2,3,和分别为背景光强和正弦条纹调制深度,Φps(x,y)为三幅模板图像T1、T2、T3中正弦条纹被目标物体调制后的相位,即正弦条纹截断相位φps(x,y)的延拓相位,是由这些模板图像中像素点(x,y)处方波条纹的灰度值Isq(x,y)在目标物体表面发生反射的结果,对于确定的(x,y),其为常数;
4f)计算机对摄像机捕获到的对应变形图像进行解码,得到这些变形图像中像素点(x,y)处正弦条纹截断相位φps(x,y)和方波条纹截断相位φps(x,y):
4f1)根据式5),将正弦条纹延拓相位Φps(x,y)表示为:
分析式10)可以知道,对于确定的像素点(x,y)来说为常数,如果对每个都加上一个常数等式的值不会改变,故将式10)变为:
通过式11),将延拓相位Φps(x,y)表示为:
其中,φps(x,y)为正弦条纹截断相位,n为正弦条纹周期延拓数;
4f2)对式9)两边进行求和,得到下式:
整理式13)得到模板图像中方波条纹被目标物体调制后的灰度值
其中为背景光强,对于一个确定的像素点(x,y)而言,其是一个常数。
根据方波条纹在幅度上加上或者减去一个常数并不改变其形状的特性,计算方波条纹截断相位φsq(x,y):
其中,为二维Gabor滤波器求相位操作,即将二维Gabor函数与进行卷积得到卷积结果G(x,y),为第i幅变形图像中(x,y)处像素点的灰度值,i=1,2,3,imag(·)为取复数的虚部操作,real(·)为取复数的实部操作。
至此,通过式12)和式15)可以分别计算正弦条纹截断相位φps(x,y)和方波条纹截断相位φsq(x,y)。
步骤4得到的截断相位如图4所示,图4(a)中表示解码出的正弦条纹截断相位φps(x,y),图4(b)表示解码出的方波条纹截断相位φsq(x,y)。
步骤5、根据正弦条纹截断相位φps和方波条纹截断相位φsq,求解得到目标物体在世界坐标系下的三维信息。
5a)计算正弦条纹相位的中间值Δxps:
其中,Tps为模板图像中正弦条纹周期,[·]表示按四舍五入方式取整;
5b)计算方波条纹相位的中间值Δxsq:
其中,Tsq为模板图像中方波条纹周期,[·]表示按四舍五入方式取整;
5c)分别计算正弦条纹周期Tps的方波条纹周期数论倒数和方波条纹周期Tsq的正弦条纹周期数论倒数
其中,Tps和Tsq分别为三幅模板T1、T2、T3中正弦条纹周期和方波条纹周期,它们是互为质数的整数,mod表示取余;
5d)计算变形图像中像素点(x,y)对应在模板图像中的位置(x0,y0):
其中,mod(a,b)表示a除以b的余数,a代表b代表Tps×Tsq。
根据编码阶段是沿着x方向进行编码的,则变形图像中像素点(x,y)与对应在模板图像中的位置(x0,y0)在y方向上保持不变,即y0=y;
5e)利用三角测距原理计算得到目标物体在世界坐标系下的三维信息:
其中,Xw,Yw,Zw分别为目标物体在世界坐标系下沿x轴、y轴、z轴的三维坐标值,世界坐标系的原点选在摄像机的光心,x轴沿摄像机成像平面的水平方向,y轴沿摄像机成像平面的垂直方向,z轴与摄像机成像平面垂直,x、y分别为变形图像中像素点(x,y)的行、列坐标,x0为变形图像中像素点(x,y)对应在模板图像中的位置(x0,y0)的行坐标,f为摄像机的焦距,b为摄像机V光心和投影仪P光心的水平距离,M为变形图像的总行数,N为变形条纹图像总列数,MT为模板图像T1、T2、T3的总行数,hcp为摄像机单个像素所代表的实际高度,其值为θV为摄像机V的垂直视场角,wcp为摄像机单个像素所代表的实际宽度,hpp为投影仪模版中单个像素所代表的实际高度。
本发明的效果通过如下仿真进一步说明:
利用上述步骤得到的目标物体的三维信息,对目标物体进行仿真重构,其效果如图5所示,从图5中可以看到,本发明重构出的目标物体不含噪声点,而且所需时间仅为传统相移法的几百分之一,因此本发明提出的基于周期互质混合编码的目标物体三维信息获取方法具有精度高、分辨率高、速度快的优点。
以上描述仅是本发明的一个具体事例,并不构成对本发明的任何限制。显然对于本领域的专业人员来说,在了解了本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变,但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。
Claims (4)
1.一种基于周期互质混合编码的目标物体三维信息获取方法,包括如下步骤:
(1)设计周期互质混合编码的三维扫描系统所需投影的三幅模板图像T1、T2、T3:
1a)设定三幅周期互质混合编码模板图像中正弦条纹直流分量为A,正弦条纹幅度为B;
1b)设定三幅模板图像T1、T2、T3的初始相移δi分别为δ1=0,δ2=2π/3,δ3=-2π/3;
1c)分别定义三幅模板中任意像素点(x,y)处的灰度值:
其中,表示模板图像中像素点(x,y)处方波条纹的灰度值,C为模板图像中方波条纹幅度,mod(a,b)是指a除以b得到的余数;Imix(x,y,δi)为第i幅模板图像中像素点(x,y)处的灰度值,i=1,2,3;Tps和Tsq分别为三幅模板图像T1、T2、T3中正弦条纹周期和方波条纹周期,它们是互为质数的整数;
(2)将投影仪P与摄像机V水平放置,并使两者光轴平行,用投影仪P将三幅模板图像T1、T2、T3依次投射到目标物体上;
(3)通过摄像机V同步拍摄三幅模板图像投影到目标物体后发生的对应变形图像并传回至计算机;
(4)计算机对摄像机捕获到的对应变形图像进行解码,得到这些变形图像中正弦条纹截断相位φps和方波条纹截断相位φsq;
(5)根据上步所解码出的正弦条纹截断相位φps和方波条纹截断相位φsq,进 行基于中国剩余定理的匹配,求解出变形图像中像素点(x,y)对应在模板图像中的位置(x0,y0),再利用三角测距原理求解得到目标物体在世界坐标系下的三维信息:
其中,Xw,Yw,Zw分别为目标物体在世界坐标系下沿x轴、y轴、z轴的三维坐标值,世界坐标系的原点选在摄像机的光心,x轴沿摄像机成像平面的水平方向,y轴沿摄像机成像平面的垂直方向,z轴与摄像机成像平面垂直,x、y分别为变形图像中像素点(x,y)的行、列坐标,x0为变形图像中像素点(x,y)对应在模板图像中的位置(x0,y0)的行坐标,f为摄像机的焦距,b为摄像机V光心和投影仪P光心的水平距离,M为变形图像的总行数,N为变形条纹图像总列数,MT为模板图像T1、T2、T3的总行数,hcp为摄像机单个像素所代表的实际高度,其值为θV为摄像机V的垂直视场角,wcp为摄像机单个像素所代表的实际宽度,hpp为投影仪模版中单个像素所代表的实际高度。
2.根据权利要求1所述的基于周期互质混合编码的目标物体三维信息获取方法,其中所述步骤(3)中的变形图像其像素点(x,y)处的灰度值为:
其中,和分别为背景光强和正弦条纹调制深度,Φps(x,y)为三幅模板图像T1、T2、T3中正弦条纹被目标物体调制后的相位,即正弦条纹截断相位φps(x,y)的延拓相位;是由这些模板图像中像素点(x,y)处方波条纹的灰度值Isq(x,y)在目标物体表面发生反射的结果,对于确定的(x,y)来说为常数。
3.根据权利要求1所述的基于周期互质混合编码的目标物体三维信息获取方法,其中步骤(4)通过如下公式计算得到:
计算正弦条纹的截断相位φps(x,y):
计算方波条纹的截断相位φsq(x,y):
其中,为二维Gabor滤波器求相位操作,即将二维Gabor函数与进行卷积得到卷积结果G(x,y),为第i幅变形图像中(x,y)处像素点的灰度值,δi为这三幅模板图像的初始相移,取δ1=0,δ2=2π/3,δ3=-2π/3,i=1,2,3,imag(·)为取复数的虚部操作,real(·)为取复数的实部操作。
4.根据权利要求1所述的基于周期互质混合编码的目标物体三维信息获取方法,其中步骤(5)所述的求解出变形图像中像素点(x,y)对应在模板图像中的位置(x0,y0),按如下步骤进行:
4a)计算正弦条纹相位的中间值Δxps:
其中,φps为正弦条纹截断相位,Tps为模板图像中正弦条纹周期,[·]表示按四舍五入方式取整;
4b)计算方波条纹相位的中间值Δxsq:
其中,φsq为方波条纹截断相位,Tsq为模板图像中方波条纹周期,[·]表示按四舍五入方式取整;
4c)分别计算正弦条纹周期Tps的方波条纹周期数论倒数和方波条纹周期Tsq的正弦条纹周期数论倒数
其中,Tps和Tsq分别为三幅模板图像T1、T2、T3中正弦条纹周期和方波条纹周期,它们是互为质数的整数,mod表示取余;
4d)计算每个像素点(x,y)对应在模板图像中的位置(x0,y0):
其中,mod(a,b)表示a除以b的余数,a代表b代表Tps×Tsq;
根据编码阶段是沿着x方向进行编码的,则变形图像中像素点(x,y)与对应在模板图像中的位置(x0,y0)在y方向上保持不变,即y0=y。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410178890.5A CN103940370B (zh) | 2014-04-29 | 2014-04-29 | 基于周期互质混合编码的目标物体三维信息获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410178890.5A CN103940370B (zh) | 2014-04-29 | 2014-04-29 | 基于周期互质混合编码的目标物体三维信息获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103940370A CN103940370A (zh) | 2014-07-23 |
CN103940370B true CN103940370B (zh) | 2016-08-17 |
Family
ID=51188141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410178890.5A Expired - Fee Related CN103940370B (zh) | 2014-04-29 | 2014-04-29 | 基于周期互质混合编码的目标物体三维信息获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103940370B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105303572B (zh) * | 2015-10-26 | 2018-05-25 | 西安电子科技大学 | 基于主被动结合的深度信息获取方法 |
CN105844633B (zh) * | 2016-03-21 | 2019-03-26 | 西安电子科技大学 | 基于De序列和相位编码的单帧结构光深度获取方法 |
CN105844623A (zh) * | 2016-03-21 | 2016-08-10 | 西安电子科技大学 | 基于De序列混合编码的目标物体深度信息获取方法 |
CN107024988A (zh) * | 2017-03-20 | 2017-08-08 | 宇龙计算机通信科技(深圳)有限公司 | 一种基于用户动作执行操作的方法及装置 |
CN107990846B (zh) * | 2017-11-03 | 2020-01-31 | 西安电子科技大学 | 基于单帧结构光的主被动结合深度信息获取方法 |
CN108534714B (zh) * | 2018-03-09 | 2019-09-24 | 南昌航空大学 | 基于正弦和二进制条纹投影的快速三维测量方法 |
CN110500957B (zh) * | 2019-09-10 | 2021-09-14 | 中国科学院苏州纳米技术与纳米仿生研究所 | 一种主动三维成像方法、装置、设备以及存储介质 |
CN113686263B (zh) * | 2021-07-23 | 2022-05-03 | 安徽农业大学 | 一种基于相移条纹投影的刚性运动物体三维测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813461A (zh) * | 2010-04-07 | 2010-08-25 | 河北工业大学 | 基于复合彩色条纹投影的绝对相位测量方法 |
CN102519390A (zh) * | 2011-12-21 | 2012-06-27 | 哈尔滨理工大学 | 三编码周期灰度梯形相移结构光三维信息获取方法 |
CN103453852A (zh) * | 2013-09-08 | 2013-12-18 | 西安电子科技大学 | 三维扫描系统中的快速相位展开方法 |
CN103697815A (zh) * | 2014-01-15 | 2014-04-02 | 西安电子科技大学 | 基于相位编码的混频结构光三维信息获取方法 |
-
2014
- 2014-04-29 CN CN201410178890.5A patent/CN103940370B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813461A (zh) * | 2010-04-07 | 2010-08-25 | 河北工业大学 | 基于复合彩色条纹投影的绝对相位测量方法 |
CN102519390A (zh) * | 2011-12-21 | 2012-06-27 | 哈尔滨理工大学 | 三编码周期灰度梯形相移结构光三维信息获取方法 |
CN103453852A (zh) * | 2013-09-08 | 2013-12-18 | 西安电子科技大学 | 三维扫描系统中的快速相位展开方法 |
CN103697815A (zh) * | 2014-01-15 | 2014-04-02 | 西安电子科技大学 | 基于相位编码的混频结构光三维信息获取方法 |
Non-Patent Citations (4)
Title |
---|
Flexible 3D shape measurement using projector defocusing: extended measurement range;Song Zhang;《OPTICS LETTERS》;20100319;第35卷(第7期);第934-936页 * |
Period Coded Phase Shifting Strategy for Real-time 3-D Structured Light Illumination;Yongchang Wang等;《IEEE TRANSACTIONS ON IMAGE PROCESSING》;20111019;第20卷(第11期);第3001-3013页 * |
Three-dimensional vision based on a combination of gray-code and phase-shift light projection:analysis and compensation of the systematic errors;Giovanna Sansoni等;《APPLIED OPTICS》;19991101;第38卷(第31期);第6565-6573页 * |
Unambiguous 3D measurement from speckle-embedded fringe;Yueyi Zhang等;《APPLIED OPTICS》;20131108;第52卷(第32期);第7797-7805页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103940370A (zh) | 2014-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103940370B (zh) | 基于周期互质混合编码的目标物体三维信息获取方法 | |
CN105303572B (zh) | 基于主被动结合的深度信息获取方法 | |
CN105203044B (zh) | 以计算激光散斑为纹理的立体视觉三维测量方法及系统 | |
CN105844633B (zh) | 基于De序列和相位编码的单帧结构光深度获取方法 | |
CN106931910A (zh) | 一种基于多模态复合编码和极线约束的高效三维图像获取方法 | |
CN103697815A (zh) | 基于相位编码的混频结构光三维信息获取方法 | |
CN103453852B (zh) | 三维扫描系统中的快速相位展开方法 | |
CN107990846B (zh) | 基于单帧结构光的主被动结合深度信息获取方法 | |
CN108596008B (zh) | 针对三维人脸测量的面部抖动补偿方法 | |
CN104680534B (zh) | 基于单帧复合模板的物体深度信息获取方法 | |
CN106643562A (zh) | 基于时域空域混合编码的结构光条纹投射方法 | |
CN106500626A (zh) | 一种手机立体成像方法及三维成像手机 | |
CN116718133A (zh) | 近距离单点结构光三维测量方法 | |
Sui et al. | Active stereo 3-D surface reconstruction using multistep matching | |
CN108895979B (zh) | 线段编码的结构光深度获取方法 | |
Yang et al. | Single-shot dense depth sensing with frequency-division multiplexing fringe projection | |
Fernandez et al. | Absolute phase mapping for one-shot dense pattern projection | |
CN115290004B (zh) | 一种基于压缩感知与hsi的水下并行单像素成像方法 | |
Cheng et al. | 3D object scanning system by coded structured light | |
Fu et al. | Single-shot colored speckle pattern for high accuracy depth sensing | |
CN105844623A (zh) | 基于De序列混合编码的目标物体深度信息获取方法 | |
Zhang et al. | Evaluation of coded structured light methods using ground truth | |
Dizeu et al. | Frequency shift triangulation: a robust fringe projection technique for 3D shape acquisition in the presence of strong interreflections | |
Fernandez et al. | A novel structured light method for one-shot dense reconstruction | |
Wijenayake et al. | An error correcting 3D scanning technique using dual pseudorandom arrays |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 |