CN106461380A - 一种基于自适应条纹的投影仪镜头畸变校正方法及其系统 - Google Patents
一种基于自适应条纹的投影仪镜头畸变校正方法及其系统 Download PDFInfo
- Publication number
- CN106461380A CN106461380A CN201680000846.2A CN201680000846A CN106461380A CN 106461380 A CN106461380 A CN 106461380A CN 201680000846 A CN201680000846 A CN 201680000846A CN 106461380 A CN106461380 A CN 106461380A
- Authority
- CN
- China
- Prior art keywords
- phase
- plane
- striped
- projector
- projection
- 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 title claims abstract description 39
- 238000012937 correction Methods 0.000 title claims abstract description 30
- 230000003044 adaptive effect Effects 0.000 title claims abstract 6
- 238000005259 measurement Methods 0.000 claims abstract description 49
- 238000013507 mapping Methods 0.000 claims abstract description 41
- 230000010363 phase shift Effects 0.000 claims abstract description 19
- 230000003287 optical effect Effects 0.000 claims abstract description 16
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 230000006978 adaptation Effects 0.000 claims description 49
- 238000009826 distribution Methods 0.000 claims description 17
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000009897 systematic effect Effects 0.000 claims description 7
- 238000002360 preparation method Methods 0.000 abstract description 2
- 238000002372 labelling Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 7
- 239000000919 ceramic Substances 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 229910052571 earthenware Inorganic materials 0.000 description 4
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 238000010587 phase diagram Methods 0.000 description 3
- 210000001215 vagina Anatomy 0.000 description 3
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 229910052573 porcelain Inorganic materials 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- 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
- G01B11/2504—Calibration devices
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供一种基于自适应条纹的投影仪镜头畸变校正方法,其中,所述方法包括:利用条纹投影测量系统获取标准平板的绝对相位;根据所述绝对相位,计算由于投影仪的镜头畸变引起的相位误差;定义垂直于投影仪镜头光轴的一平面G,根据所述相位误差计算平面G与投影平面的映射关系;根据所述映射关系,生成变频、相移自适应条纹用于投影,使得投影仪投出的条纹在平面G上为标准的正弦条纹。本发明还提供一种基于自适应条纹的投影仪镜头畸变校正系统。本发明提供的技术方案能够改善测量精度、降低了计算复杂度,能够缩短前期准备的时间,能够适用于任意设置的条纹投射测量系统。
Description
技术领域
本发明涉及光学三维测量领域,尤其涉及一种基于自适应条纹的投影仪镜头畸变校正方法及其系统。
背景技术
条纹投影轮廓术,作为光学三维测量技术的分支,由于具有非接触性、全场性、快速性、高分辨率等特点,使其在工业制造、医学诊断、文化遗产、影视娱乐等领域得到了日益广泛的应用。它的工作原理是:利用投影仪向待测物体表面投射一系列的载频条纹图(如通常为正弦条纹图),待测物体的深度信息会对条纹图调制而使之发生变形,然后利用条纹分析技术对相机采集的变形条纹图进行相位解调即可获得包含待测物体深度信息的相位图,最后根据标定好的相位深度映射关系即可由相位图求得待测物体的深度分布,实现对待测物体表面形貌的三维测量。
然而,由于待测物体的三维形貌是根据相位计算获得,相位图的质量是决定三维测量结果精度的关键因素之一。而在实际测量系统中,由于制造和装配误差,投影仪和相机的镜头畸变误差是不可避免。在这种情况下,投影仪的镜头畸变将会使得投射的正弦条纹图变形——即在垂直于投影仪光轴的平面上会得到变周期的条纹图。与定周期的条纹图(如正弦条纹图)相比较,如果将这种条纹图投射到待测物体表面,解调的相位中将会引入由于投影仪的镜头畸变引起的变形相位,导致最后的三维测量结果产生误差。
目前,针对投影仪的镜头畸变带来测量误差,现有方法大多先标定投影仪,然后根据标定的镜头畸变系数修改待投影的条纹图,使经过投影仪投射后的条纹图分布是标准的正弦条纹。然而,由于投影仪不能直接采集图像,精确地标定投影仪仍然存在一定难度。此外,现有标定方法通常假定投影平面的重投影误差满足正态分布,采用极大似然估计法标定投影仪的内外部参数,然而,由于相机像面与投影平面之间映射关系是非线性,在投影平面上无法获得具有统计特性的重投影误差,在这种情况,采用极大似然估计法标定投影仪的内外参数存在一定的误差。
因此,如何实现高精度、易操作的投影仪镜头畸变校正成为条纹投影轮廓术中的一项技术难题,亟待解决。
发明内容
有鉴于此,本发明的目的在于提供一种基于自适应条纹的投影仪镜头畸变校正方法及其系统,旨在解决现有技术中投影仪镜头畸变校正的精度较低且不易操作的问题。
本发明提出一种基于自适应条纹的投影仪镜头畸变校正方法,主要包括:
步骤一、利用条纹投影测量系统获取标准平板的绝对相位;
步骤二、根据所述绝对相位,计算由于投影仪的镜头畸变引起的相位误差;
步骤三、定义垂直于投影仪镜头光轴的一平面G,根据所述相位误差计算平面G与投影平面的映射关系;
步骤四、根据所述映射关系,生成变频、相移自适应条纹用于投影,使得投影仪投出的条纹在平面G上为标准的正弦条纹。
另一方面,本发明还提供一种基于自适应条纹的投影仪镜头畸变校正系统,所述系统包括:
相位获取模块,用于利用条纹投影测量系统获取标准平板的绝对相位;
误差计算模块,用于根据所述绝对相位,计算由于投影仪的镜头畸变引起的相位误差;
映射计算模块,用于定义垂直于投影仪镜头光轴的一平面G,根据所述相位误差计算平面G与投影平面的映射关系;
条纹生成模块,用于根据所述映射关系,生成变频、相移自适应条纹用于投影,使得投影仪投出的条纹在平面G上为标准的正弦条纹。
本发明提供的技术方案与现有技术相比,其显著优点包括:(1)本发明无需直接标定投影仪,避免了标定误差对测量结果的影响,能够改善测量精度、降低了计算复杂度;(2)本发明无需多次改变平板的空间方位和姿态,只需将标准平板任意摆放在测量空间范围内一次,即可生成自适应条纹用于投影,能够缩短前期准备的时间;(3)本发明无需计算测量系统的参数,如投影仪和相机光轴之间的夹角等,能够适用于任意设置的条纹投射测量系统。
附图说明
图1为本发明一实施方式中基于自适应条纹的投影仪镜头畸变校正方法流程图;
图2为本发明一实施方式中基于自适应条纹的投影仪镜头畸变校正系统10的内部结构示意图;
图3为本发明一实施方式中利用自适应条纹测量标准平板的结果示意图;
图4为本发明一实施方式中利用自适应条纹测量标准球的测量结果示意图;
图5为本发明一实施方式中投影仪和相机标定的重投影误差分布示意图;
图6为本发明一实施方式中利用相位辅助法标定投影仪后对陶瓷平面的测量结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明采用如下技术方案:首先利用条纹投射测量系统测量标准平板(例如平面度小于0.01mm)的绝对相位;其次根据测量的绝对相位,计算由投影仪的镜头畸变引起的相位误差;接着定义垂直于投影仪镜头光轴的某一平面G,根据相位误差建立平面G与投影平面之间的映射关系;最后利用所建立的映射关系,生成自适应条纹用于投影,即可抵消投影仪镜头畸变对条纹图的影响,在平面G上获得标准的正弦条纹。这种方法无需标定投影仪和测量系统参数,如投影仪的镜头畸变系数、投影仪和相机光轴之间的夹角等,能够降低计算复杂度、改善测量结果的精度,适用于任意设置的条纹投射测量系统。
以下将对本发明所提供的一种基于自适应条纹的投影仪镜头畸变校正方法进行详细说明。
请参阅图1,为本发明一实施方式中基于自适应条纹的投影仪镜头畸变校正方法流程图。
在步骤S1中,步骤一、利用条纹投影测量系统获取标准平板的绝对相位。
在本实施方式中,所述步骤一具体包括:
利用投影仪向标准平板分别投射一套竖直和水平的变频、相移正弦条纹图,然后采用条纹分析技术分析相机记录的正弦条纹图,并获得标准平板的绝对相位(Φm,h,Φm,v)。
在本实施方式中,将标准平板任意放置在条纹投射系统的测量空间内,利用投影仪向平板分别投射一套变频、相移的水平和竖直正弦条纹图。假设利用相机记录的第k帧条纹图强度可表示为其中,I′(i,j)表示平均强度,I″(i,j)表示调制度,Φ(i,j)表示待求的相位,δk=2πN/k表示相邻两帧条纹图之间的相移量,N表示相移次数。根据相移算法,包裹相位Φw可由下式计算求得:
其中,利用相位解包裹算法对Φw展开,可获得待测平板的两个绝对相位Φm,h和Φm,v。
在步骤S2中,步骤二、根据所述绝对相位,计算由于投影仪的镜头畸变引起的相位误差。
在本实施方式中,所述步骤二具体包括:
利用公式拟合测量的绝对相位获得载波相位分布,其中,Φc(i,j)表示在像素点(i,j)处拟合的载波相位,r,s,t,u,v表示与系统参数相关的系数;
然后将测量的绝对相位(Φm,h,Φm,v)减去拟合的载波相位Φc(i,j),即可获得由于投影仪的镜头畸变引起的相位误差(ΔΦa,h,ΔΦa,v)。
在本实施方式中,假设忽略投影仪和相机的镜头畸变影响,则平板的绝对相位可用如下公式表示Φc(i,j)表示在像素点(i,j)处拟合的载波相位,r,s,t,u,v表示与系统参数相关的系数。然而,对于实际的条纹投射系统,投影仪和相机的镜头畸变对测量结果的影响不可忽略,标准平板绝对相位的测量结果包括三部分:由于条纹投射测量系统设置引起的载波相位,以及分别由投影仪和相机的镜头畸变误差引起的相位误差。由于相机标定方法非常成熟,且精度高。因此,根据相机标定结果对采集的条纹图预处理——即去畸变校正,即可消除相机镜头畸变对测量相位的影响。投影仪的镜头畸变量通常都很小,很难直接通过裸眼观察发现投射的条纹图变形。在这种情况下,利用公式拟合测量的绝对相位,可获得载波相位分布。最后将测量的绝对相位减去拟合的载波相位,可获得由于投影仪的镜头畸变引起的相位误差ΔΦa,其中,ΔΦa=Φm-Φc。显然,这样容易获得在投射水平和竖直相移条纹后,由于投影仪的镜头畸变引起的相位误差(ΔΦa,h,ΔΦa,v)。
在步骤S3中,步骤三、定义垂直于投影仪镜头光轴的一平面G,根据所述相位误差计算平面G与投影平面的映射关系。
在本实施方式中,所述步骤三具体包括:
利用步骤一中测量的绝对相位(Φm,h,Φm,v),建立投影平面与相机像面之间等相位点对的映射关系为其中,(s1,t1)为投影平面上的像素点,(s1,t1)像素点的相位与相机像面上的像素点(i,j)的相位相等;
根据步骤二中得到的相位误差(ΔΦa,h,ΔΦa,v)建立三维数据点集(s1,t1,ΔΦa,h(i,j))和(s1,t1,ΔΦa,v(i,j));
利用泽尼克多项式(Zernike Polynomials)拟合上述两个三维数据点集(s1,t1,ΔΦa,h(i,j))和(s1,t1,ΔΦa,v(i,j)),根据拟合的系数计算投影平面上整数像素点(s,t)的相位误差(ΔΦa,h(s,t),ΔΦa,v(s,t));
通过公式建立平面G和投影平面之间的映射关系,其中,(s2,t2)为平面G在投影平面上的映射点。
在步骤S4中,根据所述映射关系,生成变频、相移自适应条纹用于投影,使得投影仪投出的条纹在平面G上为标准的正弦条纹。
在本实施方式中,所述步骤四具体包括:
定义期望在平面G上获得的条纹图的相位分布;
根据步骤三中建立的映射关系以及平面G上的相位分布,计算投影平面上整数像素点的相位
利用公式生成自适应条纹图,其中,表示自适应条纹图在像素点(s,t)处的强度,Imax表示自适应条纹图的最大强度。
本发明提供的一种基于自适应条纹的投影仪镜头畸变校正方法,首先获取标准平板的绝对相位,然后利用绝对相位计算由于投影仪的镜头畸变引起的相位误差,接着定义垂直于投影仪光轴的某一平面G,利用相位误差建立平面G与投影平面之间的映射关系,最后根据建立的映射关系,生成自适应条纹图用于投影,即可投影仪镜头畸变对条纹图的影响,在平面G上获得标准的正弦条纹图。这种方法无需标定投影仪和测量系统参数,如投影仪的镜头畸变系数、投影仪和相机光轴之间的夹角等,能够降低计算复杂度、改善测量结果的精度,适用于任意设置的条纹投射测量系统。
以下将对本发明所提供的一种基于自适应条纹的投影仪镜头畸变校正系统进行详细说明。
请参阅图2,所示为本发明一实施方式中基于自适应条纹的投影仪镜头畸变校正系统10的结构示意图。
在本实施方式中,基于自适应条纹的投影仪镜头畸变校正系统10,主要包括相位获取模块11、误差计算模块12、映射计算模块13以及条纹生成模块14。
相位获取模块11,用于利用条纹投影测量系统获取标准平板的绝对相位。
在本实施方式中,所述相位获取模块11具体用于:
利用投影仪向标准平板分别投射一套竖直和水平的变频、相移正弦条纹图,然后采用条纹分析技术分析相机记录的正弦条纹图,并获得标准平板的绝对相位(Φm,h,Φm,v)。
在本实施方式中,将标准平板任意放置在条纹投射系统的测量空间内,利用投影仪向平板分别投射一套变频、相移的水平和竖直正弦条纹图。假设利用相机记录的第k帧条纹图强度可表示为其中,I′(i,j)表示平均强度,I″(i,j)表示调制度,Φ(i,j)表示待求的相位,δk=2πN/k表示相邻两帧条纹图之间的相移量,N表示相移次数。根据相移算法,包裹相位Φw可由下式计算求得:
其中,利用相位解包裹算法对Φw展开,可获得待测平板的两个绝对相位Φm,h和Φm,v。
误差计算模块12,用于根据所述绝对相位,计算由于投影仪的镜头畸变引起的相位误差。
在本实施方式中,所述误差计算模块12具体用于:
利用公式拟合测量的绝对相位获得载波相位分布,其中,Φc(i,j)表示在像素点(i,j)处拟合的载波相位,r,s,t,u,v表示与系统参数相关的系数;
然后将测量的绝对相位(Φm,h,Φm,v)减去拟合的载波相位Φc(i,j),即可获得由于投影仪的镜头畸变引起的相位误差(ΔΦa,h,ΔΦa,v)。
在本实施方式中,假设忽略投影仪和相机的镜头畸变影响,则平板的绝对相位可用如下公式表示Φc(i,j)表示在像素点(i,j)处拟合的载波相位,r,s,t,u,v表示与系统参数相关的系数。然而,对于实际的条纹投射系统,投影仪和相机的镜头畸变对测量结果的影响不可忽略,标准平板绝对相位的测量结果包括三部分:由于条纹投射测量系统设置引起的载波相位,以及分别由投影仪和相机的镜头畸变误差引起的相位误差。由于相机标定方法非常成熟,且精度高。因此,根据相机标定结果对采集的条纹图预处理——即去畸变校正,即可消除相机镜头畸变对测量相位的影响。投影仪的镜头畸变量通常都很小,很难直接通过裸眼观察发现投射的条纹图变形。在这种情况下,利用公式拟合测量的绝对相位,可获得载波相位分布。最后将测量的绝对相位减去拟合的载波相位,可获得由于投影仪的镜头畸变引起的相位误差ΔΦa,其中,ΔΦa=Φm-Φc。显然,这样容易获得在投射水平和竖直相移条纹后,由于投影仪的镜头畸变引起的相位误差(ΔΦa,h,ΔΦa,v)。
映射计算模块13,用于定义垂直于投影仪镜头光轴的一平面G,根据所述相位误差计算平面G与投影平面的映射关系。
在本实施方式中,所述映射计算模块13具体用于:
利用所述绝对相位(Φm,h,Φm,v),建立投影平面与相机像面之间等相位点对的映射关系为其中,(s1,t1)为投影平面上的像素点,(s1,t1)像素点的相位与相机像面上的像素点(i,j)的相位相等;
根据得到的相位误差(ΔΦa,h,ΔΦa,v)建立三维数据点集(s1,t1,ΔΦa,h(i,j))和(s1,t1,ΔΦa,v(i,j));
利用泽尼克多项式拟合上述两个三维数据点集(s1,t1,ΔΦa,h(i,j))和(s1,t1,ΔΦa,v(i,j)),根据拟合的系数计算投影平面上整数像素点(s,t)的相位误差(ΔΦa,h(s,t),ΔΦa,v(s,t));
通过公式建立平面G和投影平面之间的映射关系,其中,(s2,t2)为平面G在投影平面上的映射点。
条纹生成模块14,用于根据所述映射关系,生成变频、相移自适应条纹用于投影,使得投影仪投出的条纹在平面G上为标准的正弦条纹。
在本实施方式中,所述条纹生成模块14具体用于:
定义期望在平面G上获得的条纹图的相位分布;
根据建立的映射关系以及平面G上的相位分布,计算投影平面上整数像素点的相位
利用公式生成自适应条纹图,其中,表示自适应条纹图在像素点(s,t)处的强度,Imax表示自适应条纹图的最大强度。
本发明提供的一种基于自适应条纹的投影仪镜头畸变校正系统10,无需标定投影仪和测量系统参数,如投影仪的镜头畸变系数、投影仪和相机光轴之间的夹角等,能够降低计算复杂度、改善测量结果的精度,适用于任意设置的条纹投射测量系统。
为了测试基于自适应条纹的投影仪镜头畸变误差校正方法的有效性,对一陶瓷板(平面度小于0.01mm)测量。首先利用本发明所提方法生成一套变频(条纹宽度分别为13,14,15个像素)、相移(四步相移)自适应条纹,然后将其依次投射到标准平板上,并通过相机记录条纹图,采用条纹分析技术处理条纹图像。图3(a)为实际测量相位减去载波相位后的残差分布,图3(b)为重建的平板三维数据。为了计算重建结果的精度,采用最小二乘拟合算法拟合平板的测量结果,图3(c)为测量三维数据与拟合结果的偏差分布。同时,为了与常规条纹投射测量方法比较,利用计算机生成一套标准的相移正弦条纹图,并投射到陶瓷平板上测量。图3(d)为投射标准正弦相移条纹图后恢复的相位误差。从图中可知,投影仪的镜头畸变会导致投射的条纹图变形,使得相位恢复结果中存在误差,影响三维重建结果的精度,而采用本发明所提方法能够有效抑制投影仪镜头畸变对条纹图的影响,获得准确的三维重建结果。
此外,为了检验本发明方法的重复性和精度,将陶瓷板任意摆放在11处不同的位置,在每个位置处,利用本发明所提方法向陶瓷板投射自适应条纹图,并重建三维测量结果,最后利用最小二乘拟合算法计算测量数据与拟合结果之间的偏差。表1统计了11处测量结果的偏差的PV(Peak to valley)值和RMS(Root-mean-square)值。
表1.陶瓷平面11次测量的统计结果
从表中可知,本发明所提方法能够在系统测量空间内获得稳定的测量结果。同时,为了检验本发明所提方法的精度,将陶瓷平板安装在高精度的移动导轨(分辨率为1.25μm),通过该导轨将陶瓷板分别移动在不同的位置。在每个位置处,采用本发明所提方法向陶瓷板投射自适应条纹图,并重建三维测量结果。图4为陶瓷平板中心区域某一截线轮廓的测量结果。从图中可知,利用本发明所提方法的重建误差为±0.1mm。
最后,为了与现有投影仪标定方法比较,采用相位辅助法标定投影仪。图5为投影仪和相机标定结果的重投影误差,从图中可知,投影仪标定结果的重投影误差比相机标定结果的重投影误差大,这一结果证明投影仪标定方法的精度明显低于相机的标定精度。利用上述标定结果对同一陶瓷平板测量,并采用文献中所述方法重建三维数据,图6为重建三维数据与最小二乘拟合结果的偏差分布,从图中可知重建结果误差的PV和RMS值分别为0.378mm和0.035mm。此外,图中边缘区域的重建误差明显大于中心区域,这主要是因为投影仪的标定结果存在一定的误差,假设实际像点与像面中心点的距离为d,当采用非线性畸变模型计算与实际像点对应的理想像点时,计算误差随着d的增大而增加,最终导致边缘重建结果的精度明显低于中心区域。上述测量结果间接证明了本发明所提方法能够有效抑制投影仪的镜头畸变误差,提高三维重建结果的精度。
值得注意的是,上述实施例中,所包括的各个单元只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。
另外,本领域普通技术人员可以理解实现上述各实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,相应的程序可以存储于一计算机可读取存储介质中,所述的存储介质,如ROM/RAM、磁盘或光盘等。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于自适应条纹的投影仪镜头畸变校正方法,其特征在于,所述方法包括:
步骤一、利用条纹投影测量系统获取标准平板的绝对相位;
步骤二、根据所述绝对相位,计算由于投影仪的镜头畸变引起的相位误差;
步骤三、定义垂直于投影仪镜头光轴的一平面G,根据所述相位误差计算平面G与投影平面的映射关系;
步骤四、根据所述映射关系,生成变频、相移自适应条纹用于投影,使得投影仪投出的条纹在平面G上为标准的正弦条纹。
2.如权利要求1所述的基于自适应条纹的投影仪镜头畸变校正方法,其特征在于,所述步骤一具体包括:
利用投影仪向标准平板分别投射一套竖直和水平的变频、相移正弦条纹图,然后采用条纹分析技术分析相机记录的正弦条纹图,并获得标准平板的绝对相位(Φm,h,Φm,v)。
3.如权利要求2所述的基于自适应条纹的投影仪镜头畸变校正方法,其特征在于,所述步骤二具体包括:
利用公式拟合测量的绝对相位获得载波相位分布,其中,Φc(i,j)表示在像素点(i,j)处拟合的载波相位,r,s,t,u,v表示与系统参数相关的系数;
然后将测量的绝对相位(Φm,h,Φm,v)减去拟合的载波相位Φc(i,j),即可获得由于投影仪的镜头畸变引起的相位误差(ΔΦa,h,ΔΦa,v)。
4.如权利要求3所述的基于自适应条纹的投影仪镜头畸变校正方法,其特征在于,所述步骤三具体包括:
利用步骤一中测量的绝对相位(Φm,h,Φm,v),建立投影平面与相机像面之间等相位点对的映射关系为其中,(s1,t1)为投影平面上的像素点,(s1,t1)像素点的相位与相机像面上的像素点(i,j)的相位相等;
根据步骤二中得到的相位误差(ΔΦa,h,ΔΦa,v)建立三维数据点集(s1,t1,ΔΦa,h(i,j))和(s1,t1,ΔΦa,v(i,j));
利用泽尼克多项式拟合上述两个三维数据点集(s1,t1,ΔΦa,h(i,j))和(s1,t1,ΔΦa,v(i,j)),根据拟合的系数计算投影平面上整数像素点(s,t)的相位误差(ΔΦa,h(s,t),ΔΦa,v(s,t));
通过公式建立平面G和投影平面之间的映射关系,其中,(s2,t2)为平面G在投影平面上的映射点。
5.如权利要求4所述的基于自适应条纹的投影仪镜头畸变校正方法,其特征在于,所述步骤四具体包括:
定义期望在平面G上获得的条纹图的相位分布;
根据步骤三中建立的映射关系以及平面G上的相位分布,计算投影平面上整数像素点的相位
利用公式生成自适应条纹图,其中,表示自适应条纹图在像素点(s,t)处的强度,Imax表示自适应条纹图的最大强度。
6.一种基于自适应条纹的投影仪镜头畸变校正系统,其特征在于,所述系统包括:
相位获取模块,用于利用条纹投影测量系统获取标准平板的绝对相位;
误差计算模块,用于根据所述绝对相位,计算由于投影仪的镜头畸变引起的相位误差;
映射计算模块,用于定义垂直于投影仪镜头光轴的一平面G,根据所述相位误差计算平面G与投影平面的映射关系;
条纹生成模块,用于根据所述映射关系,生成变频、相移自适应条纹用于投影,使得投影仪投出的条纹在平面G上为标准的正弦条纹。
7.如权利要求6所述的基于自适应条纹的投影仪镜头畸变校正系统,其特征在于,所述相位获取模块具体用于:
利用投影仪向标准平板分别投射一套竖直和水平的变频、相移正弦条纹图,然后采用条纹分析技术分析相机记录的正弦条纹图,并获得标准平板的绝对相位(Φm,h,Φm,v)。
8.如权利要求7所述的基于自适应条纹的投影仪镜头畸变校正系统,其特征在于,所述误差计算模块具体用于:
利用公式拟合测量的绝对相位获得载波相位分布,其中,Φc(i,j)表示在像素点(i,j)处拟合的载波相位,r,s,t,u,v表示与系统参数相关的系数;
然后将测量的绝对相位(Φm,h,Φm,v)减去拟合的载波相位Φc(i,j),即可获得由于投影仪的镜头畸变引起的相位误差(ΔΦa,h,ΔΦa,v)。
9.如权利要求8所述的基于自适应条纹的投影仪镜头畸变校正系统,其特征在于,所述映射计算模块具体用于:
利用所述绝对相位(Φm,h,Φm,v),建立投影平面与相机像面之间等相位点对的映射关系为其中,(s1,t1)为投影平面上的像素点,(s1,t1)像素点的相位与相机像面上的像素点(i,j)的相位相等;
根据得到的相位误差(ΔΦa,h,ΔΦa,v)建立三维数据点集(s1,t1,ΔΦa,h(i,j))和(s1,t1,ΔΦa,v(i,j));
利用泽尼克多项式拟合上述两个三维数据点集(s1,t1,ΔΦa,h(i,j))和(s1,t1,ΔΦa,v(i,j)),根据拟合的系数计算投影平面上整数像素点(s,t)的相位误差(ΔΦa,h(s,t),ΔΦa,v(s,t));
通过公式建立平面G和投影平面之间的映射关系,其中,(s2,t2)为平面G在投影平面上的映射点。
10.如权利要求9所述的基于自适应条纹的投影仪镜头畸变校正系统,其特征在于,所述条纹生成模块具体用于:
定义期望在平面G上获得的条纹图的相位分布;
根据建立的映射关系以及平面G上的相位分布,计算投影平面上整数像素点的相位
利用公式生成自适应条纹图,其中,表示自适应条纹图在像素点(s,t)处的强度,Imax表示自适应条纹图的最大强度。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2016/097664 WO2018040017A1 (zh) | 2016-08-31 | 2016-08-31 | 一种基于自适应条纹的投影仪镜头畸变校正方法及其系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106461380A true CN106461380A (zh) | 2017-02-22 |
CN106461380B CN106461380B (zh) | 2019-04-05 |
Family
ID=58215970
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201680000846.2A Active CN106461380B (zh) | 2016-08-31 | 2016-08-31 | 一种基于自适应条纹的投影仪镜头畸变校正方法及其系统 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN106461380B (zh) |
WO (1) | WO2018040017A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112050751A (zh) * | 2020-07-17 | 2020-12-08 | 深圳大学 | 一种投影仪标定方法、智能终端及存储介质 |
CN112184788A (zh) * | 2020-09-16 | 2021-01-05 | 西安邮电大学 | 一种四步相移的主值相位提取方法 |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109798841B (zh) * | 2019-03-05 | 2023-09-22 | 盎锐(上海)信息科技有限公司 | 相机与投影仪的标定系统及方法 |
CN111156928B (zh) * | 2020-02-07 | 2023-04-18 | 武汉玄景科技有限公司 | 一种基于dlp投影的光栅三维扫描仪莫尔条纹消除方法 |
CN111476788B (zh) * | 2020-04-27 | 2023-08-25 | 武汉精立电子技术有限公司 | 一种显示屏夹层缺陷检测方法及系统 |
CN111784842B (zh) * | 2020-06-29 | 2024-04-12 | 北京百度网讯科技有限公司 | 三维重建方法、装置、设备及可读储存介质 |
CN111768341B (zh) * | 2020-07-01 | 2023-09-19 | 成都国科微电子有限公司 | 一种图像处理方法及图像处理装置 |
CN112070842B (zh) * | 2020-07-28 | 2023-03-21 | 安徽农业大学 | 一种基于正交编码条纹的多摄像机全局标定方法 |
CN112014408B (zh) * | 2020-09-04 | 2024-04-12 | 东莞市盟拓智能科技有限公司 | 一种基于结构光原理对pcb板进行重建的检测方法 |
CN112581605B (zh) * | 2020-12-24 | 2024-05-14 | 西安中科光电精密工程有限公司 | 一种结构光三维重建校正方法及装置 |
CN113191963B (zh) * | 2021-04-02 | 2022-08-05 | 华中科技大学 | 一种无需附加操作的投影仪残余畸变全场标定方法及装置 |
CN113358062B (zh) * | 2021-05-31 | 2022-08-09 | 湖北工业大学 | 三维重建相位误差补偿方法 |
CN115063341B (zh) * | 2022-04-02 | 2024-04-26 | 安徽农业大学 | 一种基于相移条纹的苹果果梗和花萼检测方法 |
CN115115788B (zh) * | 2022-08-12 | 2023-11-03 | 梅卡曼德(北京)机器人科技有限公司 | 三维重建方法、装置及电子设备、存储介质 |
CN115393507B (zh) * | 2022-08-17 | 2023-12-26 | 梅卡曼德(北京)机器人科技有限公司 | 三维重建方法、装置及电子设备、存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1242906A (zh) * | 1997-08-27 | 2000-01-26 | Ldt激光展示技术公司 | 视频图象中的几何图象误差的补偿方法及实施此方法的装置 |
US20050024598A1 (en) * | 2003-07-30 | 2005-02-03 | Casio Computer Co., Ltd. | Projector and method of projecting projection image |
CN101419708A (zh) * | 2008-12-08 | 2009-04-29 | 北京航空航天大学 | 一种基于一维靶标的结构光参数标定方法 |
CN101949693A (zh) * | 2010-08-03 | 2011-01-19 | 河北工业大学 | 一种三维成像系统的标定方法 |
US20110013247A1 (en) * | 2008-03-26 | 2011-01-20 | Hoya Corporation | Two-dimensional scanning projector |
CN102508578A (zh) * | 2011-10-09 | 2012-06-20 | 清华大学深圳研究生院 | 投影定位装置及方法、交互系统和交互方法 |
CN105091750A (zh) * | 2015-07-30 | 2015-11-25 | 河北工业大学 | 一种基于双四步相移的投影仪标定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9459094B2 (en) * | 2013-09-12 | 2016-10-04 | Hong Kong Applied Science and Technology Research Institute Company Limited | Color-encoded fringe pattern for three-dimensional shape measurement |
-
2016
- 2016-08-31 CN CN201680000846.2A patent/CN106461380B/zh active Active
- 2016-08-31 WO PCT/CN2016/097664 patent/WO2018040017A1/zh active Application Filing
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1242906A (zh) * | 1997-08-27 | 2000-01-26 | Ldt激光展示技术公司 | 视频图象中的几何图象误差的补偿方法及实施此方法的装置 |
US20050024598A1 (en) * | 2003-07-30 | 2005-02-03 | Casio Computer Co., Ltd. | Projector and method of projecting projection image |
US20110013247A1 (en) * | 2008-03-26 | 2011-01-20 | Hoya Corporation | Two-dimensional scanning projector |
CN101419708A (zh) * | 2008-12-08 | 2009-04-29 | 北京航空航天大学 | 一种基于一维靶标的结构光参数标定方法 |
CN101949693A (zh) * | 2010-08-03 | 2011-01-19 | 河北工业大学 | 一种三维成像系统的标定方法 |
CN102508578A (zh) * | 2011-10-09 | 2012-06-20 | 清华大学深圳研究生院 | 投影定位装置及方法、交互系统和交互方法 |
CN105091750A (zh) * | 2015-07-30 | 2015-11-25 | 河北工业大学 | 一种基于双四步相移的投影仪标定方法 |
Non-Patent Citations (1)
Title |
---|
崔艳军 等: "《条纹投影三维测量的Gamma畸变校正方法》", 《光学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112050751A (zh) * | 2020-07-17 | 2020-12-08 | 深圳大学 | 一种投影仪标定方法、智能终端及存储介质 |
CN112050751B (zh) * | 2020-07-17 | 2022-07-22 | 深圳大学 | 一种投影仪标定方法、智能终端及存储介质 |
CN112184788A (zh) * | 2020-09-16 | 2021-01-05 | 西安邮电大学 | 一种四步相移的主值相位提取方法 |
CN112184788B (zh) * | 2020-09-16 | 2023-11-07 | 西安邮电大学 | 一种四步相移的主值相位提取方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2018040017A1 (zh) | 2018-03-08 |
CN106461380B (zh) | 2019-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106461380B (zh) | 一种基于自适应条纹的投影仪镜头畸变校正方法及其系统 | |
CN107607060B (zh) | 一种应用于光栅三维投影测量中的相位误差补偿方法 | |
CN105783775B (zh) | 一种镜面及类镜面物体表面形貌测量装置与方法 | |
US20220221270A1 (en) | A calibration method for fringe projection systems based on plane mirrors | |
Zhang | Review of single-shot 3D shape measurement by phase calculation-based fringe projection techniques | |
CN104330052B (zh) | 外差式三频不等步相移解相位方法 | |
Karpinsky et al. | High-resolution, real-time 3D imaging with fringe analysis | |
Zappa et al. | Static and dynamic features of Fourier transform profilometry: A review | |
CN109059806B (zh) | 一种基于红外条纹的镜面物体三维面形测量装置及方法 | |
CN109141291A (zh) | 一种快速相位解包裹算法 | |
CN103994732B (zh) | 一种基于条纹投影的三维测量方法 | |
Mao et al. | A multi-frequency inverse-phase error compensation method for projector nonlinear in 3D shape measurement | |
CN103791856A (zh) | 基于四幅光栅条纹图像的相位求解与去包裹方法 | |
CN108955571A (zh) | 双频外差与相移编码相结合的三维测量方法 | |
Zhang et al. | Full-field phase error analysis and compensation for nonsinusoidal waveforms in phase shifting profilometry with projector defocusing | |
Xiaoling et al. | Calibration of a fringe projection profilometry system using virtual phase calibrating model planes | |
CN112097670B (zh) | 高反光物体三维面型测量方法和测量设备 | |
Yu et al. | High sensitivity fringe projection profilometry combining optimal fringe frequency and optimal fringe direction | |
Wang et al. | Motion-induced error reduction for phase-shifting profilometry with phase probability equalization | |
Li et al. | Complex object 3D measurement based on phase-shifting and a neural network | |
Li et al. | An improved 2+ 1 phase-shifting algorithm | |
Zhu et al. | Alternate iterative least-squares algorithm based on nonuniform phase shifting for suppressing nonlinearity errors in fringe projection profilometry | |
Wang et al. | Nonlinear high-order harmonics correction for phase measuring profilometry | |
KR101001894B1 (ko) | 컬러 투사 모아레 기법을 이용한 3차원 형상 측정 장치 및 방법 | |
Fu et al. | Least-squares calibration method for fringe projection profilometry with some practical considerations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |