CN110207621B - 一种基于变相移的四步相移的面结构光解相方法 - Google Patents

一种基于变相移的四步相移的面结构光解相方法 Download PDF

Info

Publication number
CN110207621B
CN110207621B CN201910482281.1A CN201910482281A CN110207621B CN 110207621 B CN110207621 B CN 110207621B CN 201910482281 A CN201910482281 A CN 201910482281A CN 110207621 B CN110207621 B CN 110207621B
Authority
CN
China
Prior art keywords
phase
phase shift
image
acquired image
variable
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
CN201910482281.1A
Other languages
English (en)
Other versions
CN110207621A (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN201910482281.1A priority Critical patent/CN110207621B/zh
Publication of CN110207621A publication Critical patent/CN110207621A/zh
Application granted granted Critical
Publication of CN110207621B publication Critical patent/CN110207621B/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
    • G01B11/2518Projection by scanning of the object
    • G01B11/2527Projection by scanning of the object with phase change by in-plane movement of the patern

Abstract

本发明公开了一种基于变相移的四步相移的面结构光解相方法,通过生成一组包含变相移的四步相移信息的条纹光栅图像,并通过相对应的相位解算算法对被采集图像进行相位信息的重新解算而得到相位信息。该基于变相移的四步相移的面结构光解相方法仅使用四幅竖直条纹光栅图像,不需要更多的图像或者额外的预知信息,在保证结构光三维测量中本身的精度和优点下,提高了测量速度。同时,本发明方法在抗噪性上也有较好的效果,算法具有较好的鲁棒性。

Description

一种基于变相移的四步相移的面结构光解相方法
技术领域
本发明涉及一种面结构光解相方法,具体涉及一种基于变相移的四步相移的面结构光解相方法。
背景技术
随着科技的进步,高精度的三维测量技术需求越来越大,根据测量时是否与被测量物体接触,可将三维测量技术分为接触式和非接触式。非接触式三维测量技术同样具备高精度的优点,同时还具备传统接触式测量技术所不具备的优点:如无破坏性、高效率、工作距离大等。近年来非接触式三维测量技术在文物保护、逆向工程、虚拟现实等领域也得到越来越多的应用。
在非接触式三维测量中,基于相位结构光的主动式三维测量方法以高精度、高速度、易于实现等优点逐渐成为三维测量领域的研究热门。
基于面结构光的三维测量原理:在由投影仪和相机组成的结构光测量系统中,通过投影仪将具有已知相位信息的面结构光投射到被测物体表面,再由相机同步采集经过被测物体表面高度调制后的条纹光栅图像,通过对采集的条纹光栅图像的解相来得到相应的相位分布,再根据系统标定中得到的相位和三维坐标的关系重建出物体的三维形貌。
存在的问题:在面结构光三维测量的实际应用中,从采集的条纹光栅图像中解出准确的相位信息对三维测量和重建精度有直接的影响。传统的四步相移多频投影测量方式,在计算相位级整数部分的时候,由于随机噪声、非线性误差和过亮过暗等因素的影响,在求取相位级数的时候,由于误差的不断传递放大,导致所求得的相位级数整数部分存在跳跃性误差,导致最终求出相位中出现2pi的整数倍跳跃误差,造成相位缺失,在实际测量中,这种跳跃误差往往较多且不易消除。
发明内容
有鉴于此,本发明的目的在于提供一种基于变相移的四步相移的面结构光解相方法,该方法的主要内容是通过生成一组包含变相移的四步相移信息的条纹光栅图像,并通过相对应的相位解算算法对被采集图像进行相位信息的重新解算而得到相位信息。该基于变相移的四步相移的面结构光解相方法仅使用四幅图像,不需要更多的图像或者额外的预知信息。因此,本发明方法可以灵活的应用于高速三维测量。
本发明的一种基于变相移的四步相移的面结构光解相方法,包括如下步骤:
步骤1:确定用于解包裹相位的节距p;
步骤2:根据步骤1确定的节距生成一组变相移的四步相移条纹光栅图像,分别记为投影图像1、投影图像2、投影图像3和投影图像4;
步骤3:将投影图像1、投影图像2、投影图像3和投影图像4分别通过投影仪投射到物体表面,并通过相机采集图像,获得采集图像1、采集图像2、采集图像3和采集图像4;
步骤4:根据步骤3中获得采集图像1、采集图像2、采集图像3和采集图像4,计算对应采集图像的各点相移α(i,j);
步骤5:根据步骤4中获得相移α(i,j),计算采集图像的条纹顺序k;
步骤6:根据步骤4中获得相移α(i,j),计算采集图像的包裹相位φ'(i,j);
步骤7:根据步骤5中获得条纹顺序k和步骤6中获得包裹相位φ'(i,j),计算采集图像的绝对相位φ(i,j)。
进一步,其中步骤2中,所述包含有变相移的四步相移的条纹光栅图像通过下列(1)至(4)式求得:
I1(i,j)=A+Bcos[φ(i,j)+α(i,j)] (1)
I2(i,j)=A-Bcos[φ(i,j)+α(i,j)] (2)
I3(i,j)=A-Bsin[φ(i,j)+α(i,j)] (3)
I4(i,j)=A+Bsin[φ(i,j)-α(i,j)] (4)
其中,In(i,j)表示该组中第n幅图中(i,j)点处的灰度值大小,n=1,2,3,4;A代表图像的平均灰度值,B代表图像的灰度调制;φ(i,j)为该组图像中对应的各点绝对相位,α(i,j)为该组图像中对应的各点相移。
进一步,其中步骤4中,采集图像中对应的各点相移通过下列(5)式求得:
Figure GDA0002803574290000021
其中,
Figure GDA0002803574290000022
进一步,其中步骤5中,采集图像中对应的条纹顺序k通过下列(6)式求得:
Figure GDA0002803574290000023
其中,round表示四舍五入取整运算。
进一步,其中步骤6中,包裹相位φ'(i,j)通过下列(7)式求得:
φ'(i,j)=arctan2[(I3-A),(I2-A)]-α(i,j) (7)
其中,
Figure GDA0002803574290000031
进一步,其中步骤7中,绝对相位φ(i,j)通过下列(8)式求得:
φ(i,j)=k×2π+φ'(i,j) (8)
本发明提供的一种基于变相移的四步相移的面结构光解相方法有益效果在于:该方法将生成的图像表示为具有量化相移的特殊正弦形式,仅需要四幅图像同时提取包裹相位和条纹顺序。此外,对于数字条纹投影系统而言,该方法不需要任何预先知道的信息,最终逐像素的获得绝对相位图,可提高整个快速测量过程的操作效率。该解相方法操作简单,能够提高三维快速测量的速度,适合大规模的推广和应用。
附图说明
下面结合附图和实施例对本发明作进一步描述:
图1是本发明的测量系统的示意图。
图2是本发明的该基于变相移的四步相移的面结构光解相方法的流程示意图。
图3是本发明的该基于变相移的四步相移的面结构光解相方法的变相移四步相移条纹光栅图像生成示意图。
具体实施方式
如图1至图3所示是依本发明的发明精神提供的一种基于变相移的四步相移的面结构光解相方法,其通过一个结构光测量系统来实现,该测量系统包括一个工业相机1、一个投影仪2以及一个控制终端3,该工业相机1和该投影仪2布置在被测物体4前方,并且该工业相机1和该投影仪2均连接于该控制终端3。该控制终端3能够控制该投影仪2向该被测量物体投射具有相位信息的条纹光栅,该工业相机1能够实时采集经过该被测量物体表面调制的条纹图像,并将其发送到该终端3以进行分析和后续的操作,例如进行该方法的操作。
结合图1所示,本发明公开一种基于变相移的四步相移的面结构光解相方法,具体实现步骤如下:
步骤1:如图1所示的面结构光测量系统,主要由该工业相机1、该投影仪2、该控制终端3组成,其中该工业相机1和该投影仪2布置在被测物体4前方,其布置位置应保证该工业相机1拍摄范围覆盖该投影仪2的投影范围。在该优选的实施例中选取新颖的变相移的四步相移解相方法,确定展开频率数M=1。根据相关研究,选用于解包裹相位的节距p,则该优选的实施例需采集一组该节距下的图像,该组图像共有四幅竖直条纹光栅图像。设定每幅待投射图像的宽度和高度分别为W、H。
步骤2:依据本发明的发明精神提供的一种基于变相移的四步相移的面结构光解相方法,该优选实施例设定生成的结构光条纹光栅表达式中待定参数:
A、B
Figure GDA0002803574290000041
Figure GDA0002803574290000042
其中,p是条纹节距,代表着在一个周期内的像素总数。N代表一幅完整图像中条纹总的个数。int()表示对括号内的数进行四舍五入取整操作。
则对应于p节距的该组竖直条纹光栅(i,j)处光强的表达式通过下列(1)至(4)式求得:
Figure GDA0002803574290000043
Figure GDA0002803574290000044
Figure GDA0002803574290000045
Figure GDA0002803574290000046
根据上述该步骤1中预设p,可生成一组变相移的四步相移竖直条纹光栅图像。对应于该优选的实施例生成的结构光条纹光栅,展示在附图3中。
步骤3:根据本发明的发明精神提供的一种基于变相移的四步相移的结构光解相方法,由附图1中该投影仪2将步骤2中生成的四张结构光条纹光栅图像投射至被测物体4,通过该工业相机1同步采集经过被测物体表面调制后相位变形的条纹图像,并将采集图像存储至该控制终端3。
步骤4:依据本发明的发明精神提供的一种基于变相移的四步相移的面结构光解相方法,计算步骤3中被采集图像的相位。具体过程为:选取图像中有效信息区域,解算节距p的四幅图像的水平包裹相位。设该优选实施例所得各组采集图像的灰度值为:
对应p的采集图像光栅灰度值:
Figure GDA0002803574290000047
其中,v表示竖直条纹光栅。根据上述假设所采集图像的水平相移分布通过下列(5)式求得:
Figure GDA0002803574290000051
其中,
Figure GDA0002803574290000052
步骤5:基于步骤4中获得相移α(i,j),条纹顺序k通过下列(6)式求得:
Figure GDA0002803574290000053
其中,round表示四舍五入取整运算。
步骤6:基于步骤4中获得相移α(i,j),包裹相位φ'(i,j)通过下列(7)式求得:
φ'(i,j)=arctan2[(I3 v-A),(I2 v-A)]-α(i,j) (7)
其中,
Figure GDA0002803574290000054
步骤7:基于步骤5中获得条纹顺序k和步骤6中获得包裹相位φ'(i,j),绝对相位φ(i,j)通过下列(8)式求得:
φ(i,j)=k×2π+φ'(i,j) (8)
总的来说,本发明在面结构光三维测量系统中,提供一种基于变相移的四步相移的面结构光解相方法,该方法使用一种图像频率,生成一组包含变相移的四步相移编码信息的条纹光栅图像,精确测量三维信息。该基于变相移的四步相移的面结构光解相方法仅使用四幅图像,不需要更多的图像或者额外的预知信息,测量速度更快。由于没有数字条纹投影系统的预设几何约束,该方法的测量范围更加广泛。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (5)

1.一种基于变相移的四步相移的面结构光解相方法,其特征在于,包括如下步骤:
步骤1:确定用于解包裹相位的节距p;
步骤2:根据步骤1确定的节距生成一组变相移的四步相移条纹光栅图像,分别记为投影图像1、投影图像2、投影图像3和投影图像4;所述包含有变相移的四步相移的条纹光栅图像通过下列(1)至(4)式求得:
I1(i,j)=A+B cos[φ(i,j)+α(i,j)] (1)
I2(i,j)=A-B cos[φ(i,j)+α(i,j)] (2)
I3(i,j)=A-B sin[φ(i,j)+α(i,j)] (3)
I4(i,j)=A+B sin[φ(i,j)-α(i,j)] (4)
其中,In(i,j)表示该组中第n幅图中(i,j)点处的灰度值大小,n=1,2,3,4;A代表图像的平均灰度值,B代表图像的灰度调制;φ(i,j)为该组图像中对应的各点绝对相位,α(i,j)为该组图像中对应的各点相移;
步骤3:将投影图像1、投影图像2、投影图像3和投影图像4分别通过投影仪投射到物体表面,并通过相机采集图像,获得采集图像1、采集图像2、采集图像3和采集图像4;
步骤4:根据步骤3中获得采集图像1、采集图像2、采集图像3和采集图像4,计算对应采集图像的各点相移α(i,j);
步骤5:根据步骤4中获得各点相移α(i,j),计算采集图像的条纹顺序k;
步骤6:根据步骤4中获得各点相移α(i,j),计算采集图像的包裹相位φ'(i,j);
步骤7:根据步骤5中获得条纹顺序k和步骤6中获得包裹相位φ'(i,j),计算采集图像的绝对相位φ(i,j)。
2.根据权利要求1所述的一种基于变相移的四步相移的面结构光解相方法,其特征在于:其中步骤4中,采集图像中对应的各点相移通过下列(5)式求得:
Figure FDA0002803574280000011
其中,
Figure FDA0002803574280000012
3.根据权利要求1所述的一种基于变相移的四步相移的面结构光解相方法,其特征在于:其中步骤5中,采集图像中对应的条纹顺序k通过下列(6)式求得:
Figure FDA0002803574280000021
其中,round表示四舍五入取整运算,f为图像的条纹周期数。
4.根据权利要求1所述的一种基于变相移的四步相移的面结构光解相方法,其特征在于:其中步骤6中,包裹相位φ'(i,j)通过下列(7)式求得:
φ'(i,j)=arctan2[(I3-A),(I2-A)]-α(i,j) (7)
其中,
Figure FDA0002803574280000022
5.根据权利要求1所述的一种基于变相移的四步相移的面结构光解相方法,其特征在于:其中步骤7中,绝对相位φ(i,j)通过下列(8)式求得:
φ(i,j)=k×2π+φ'(i,j) (8)
其中,k为条纹顺序,φ'(i,j)为包裹相位。
CN201910482281.1A 2019-06-04 2019-06-04 一种基于变相移的四步相移的面结构光解相方法 Active CN110207621B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910482281.1A CN110207621B (zh) 2019-06-04 2019-06-04 一种基于变相移的四步相移的面结构光解相方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910482281.1A CN110207621B (zh) 2019-06-04 2019-06-04 一种基于变相移的四步相移的面结构光解相方法

Publications (2)

Publication Number Publication Date
CN110207621A CN110207621A (zh) 2019-09-06
CN110207621B true CN110207621B (zh) 2021-03-26

Family

ID=67790742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910482281.1A Active CN110207621B (zh) 2019-06-04 2019-06-04 一种基于变相移的四步相移的面结构光解相方法

Country Status (1)

Country Link
CN (1) CN110207621B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6100984A (en) * 1999-06-11 2000-08-08 Chen; Fang Surface measurement system with a laser light generator
JP4830871B2 (ja) * 2007-01-26 2011-12-07 パナソニック電工株式会社 3次元形状計測装置及び3次元形状計測方法
JP2010039364A (ja) * 2008-08-07 2010-02-18 Sony Corp 画像表示装置および画像表示方法、画像処理装置および画像処理方法、並びにプログラム
CN105698708B (zh) * 2016-01-29 2018-06-01 四川大学 一种三维视觉重建方法
CN108061517B (zh) * 2017-12-21 2020-08-11 重庆大学 基于莫尔序列光栅的面结构光解相方法
CN108955571B (zh) * 2018-05-16 2019-07-09 南昌航空大学 双频外差与相移编码相结合的三维测量方法
CN109523627B (zh) * 2018-10-24 2023-07-18 重庆大学 一种基于泰勒指数表达式型面结构光的三维重建方法

Also Published As

Publication number Publication date
CN110207621A (zh) 2019-09-06

Similar Documents

Publication Publication Date Title
CN108534714B (zh) 基于正弦和二进制条纹投影的快速三维测量方法
CN109489585B (zh) 基于改进多频条纹结构光的三维测量方法
CN108955574B (zh) 一种三维测量方法及系统
CN101813461B (zh) 基于复合彩色条纹投影的绝对相位测量方法
CN108955571B (zh) 双频外差与相移编码相结合的三维测量方法
CN109523627B (zh) 一种基于泰勒指数表达式型面结构光的三维重建方法
Song et al. A new phase unwrapping algorithm based on three wavelength phase shift profilometry method
CN109307483B (zh) 一种基于结构光系统几何约束的相位展开方法
CN105300317B (zh) 基于正弦和三角波条纹投影的三维测量方法
CN110174079B (zh) 一种基于四步相移编码型面结构光的三维重建方法
CN111174730B (zh) 一种基于相位编码的快速相位解包裹方法
CN110207620B (zh) 一种通过不同频率确定数字光栅投影结构光级数的三维重建方法
CN107036556A (zh) 基于分段量化相位编码的结构光三维测量方法
CN107339954A (zh) 基于周期不同步正弦条纹加相位编码条纹的三维测量方法
CN108061517B (zh) 基于莫尔序列光栅的面结构光解相方法
Yu et al. Unequal-period combination approach of gray code and phase-shifting for 3-D visual measurement
CN110849290A (zh) 基于形态学操作的分段量化编码强度的三维测量方法
CN112880589A (zh) 基于双频相位编码的光学三维测量方法
CN110428459A (zh) 一种基于数字顺序编码的相位解包裹的方法
CN106840039B (zh) 一种基于复合光栅投影的快速三维测量方法
CN110207621B (zh) 一种基于变相移的四步相移的面结构光解相方法
CN114252026B (zh) 调制三维编码于周期边缘的三维测量方法及系统
CN110440714A (zh) 一种基于多频和二进制条纹的相位解包裹方法
CN114234852B (zh) 基于最佳映射点集匹配的多目结构光三维测量方法及系统
CN114152203B (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