CN101881605B - 基于相位编码技术的光学三维测量方法 - Google Patents

基于相位编码技术的光学三维测量方法 Download PDF

Info

Publication number
CN101881605B
CN101881605B CN201010190225XA CN201010190225A CN101881605B CN 101881605 B CN101881605 B CN 101881605B CN 201010190225X A CN201010190225X A CN 201010190225XA CN 201010190225 A CN201010190225 A CN 201010190225A CN 101881605 B CN101881605 B CN 101881605B
Authority
CN
China
Prior art keywords
grating
phase
coding
optical
encoding
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
CN201010190225XA
Other languages
English (en)
Other versions
CN101881605A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201010190225XA priority Critical patent/CN101881605B/zh
Publication of CN101881605A publication Critical patent/CN101881605A/zh
Application granted granted Critical
Publication of CN101881605B publication Critical patent/CN101881605B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

一种基于相位编码技术的光学三维测量方法,属三维测量技术领域。由以下步骤组成:(1)计算机生成光栅编码条纹图像;(2)用数字投影设备投射光栅编码条纹图像;(3)采集光栅条纹图像并利用相移算法进行编码求解和三维测量;其特征在于:步骤(1)中生成的所述光栅编码条纹图像的每个光栅周期内像素个数等于所采用的相移算法步数;且所述光栅编码条纹图像的每个光栅周期内像素灰度需满足正弦或余弦曲线分布,且像素灰度值为相应正弦或余弦函数的极值点。本方法相对传统正弦相移编码方法,对数字投影设备的伽玛非线性影响具有更强的抗干扰性,编码具有更高稳定性、可靠性和正确性,在物体表面光学三维轮廓精确测量方面有着重要作用。

Description

基于相位编码技术的光学三维测量方法
所属技术领域
一种基于相位编码技术的光学三维测量方法,涉及计算机视觉测量领域,属于三维测量方法和仪器技术领域,尤其涉及光学三维测量系统的结构光编码技术。
背景技术
基于光学的非接触式三维物体形貌测量技术得到了长足的发展,该技术已经广泛应用于产品设计与制造、质量检测与控制、机器人视觉等领域。另外,随着数字投影技术的进步,基于数字投影设备的结构光三维测量技术得到快速的扩展,在整个光学三维测量技术领域中占据越来越重要的地位。该方法:首先利用计算机生成光栅编码条纹图像;再使用数字投影设备投射光栅编码条纹图像;最后采集光栅条纹图像并进行编码求解和三维测量。
基于正弦的相移算法具有较高的测量精度,使得该种方法在光学轮廓测量系统中得到广泛应用,一般正弦条纹光强函数用数学表达式可表示为:
Figure BSA00000146239300011
其中Ii(x,y)为相机采集的第i幅光栅条纹图像中像素坐标为(x,y)的光强值,I′(x,y)为背景光强,I″(x,y)为条纹幅值,
Figure BSA00000146239300012
为待测相位值,δi为第i次相移量。经过多步相移后待测相位
Figure BSA00000146239300013
计算值为:
Figure BSA00000146239300014
其中ai和bi取决于相移的次数,求解得
Figure BSA00000146239300015
在主值[-π,+π]内,呈现出有2π锯齿反复跳动的包裹相位,经过相位展开算法求解得整张图像内相位的唯一编码。
设相机采集到的理想的正弦条纹灰度表示为
Figure BSA00000146239300021
设经数字化投影设备投射出,相机实际采集到的条纹灰度表示为
Figure BSA00000146239300022
由于一般的商用数字投影设备的投影引擎具有伽玛非线性特性,使得采集到的条纹灰度产生非正弦性误差,不能得到完全理想的正弦灰度分布,使用这种非理想的光栅灰度条纹进行计算时,会产生相位编码误差,进而导致三维测量精度大大降低,其中,相位编码误差以非正弦性波动的形式表现出来。
针对这种误差相继出现了多种解决方法,主要可以分为两类:
一类是主动式的方法(盖绍彦,达飞鹏.基于数字投影仪的光栅相位自校正方法,自动化学报,34(11):1363-1368(2008).),在计算出的相位中选取局部的误差区域,通过逆向补偿的方式计算出需要投射的条纹灰度,通过改变投射的条纹灰度,从而使采集的到的条纹灰度呈理想的正弦分布,但这类方法理论上可行,实际应用中由于每个像素具有不同的伽马非线性反射,使得该方法很难得到很好的结果。
第二类方法是被动补偿的方法(H.Guo,H.He,and M.Chen.Gamma correction fordigital fringe projection profilometry,Appl.Opt.43,2906-14(2004).),这类方法认为相位编码误差具有相同的周期性,首先提取一个周期或几个周期内的相位误差,建立对应相位的误差补偿相位对照表,实现相位误差标定,然后对求解的相位编码根据误差对照表进行逆向补偿,实现相位的正确编码,这类方法的前提条件是相位编码误差具有周期相同性,通过对一个或多个周期内的相位误差进行统计,建立误差对照表进行对应补偿,然而,实现应用中这种完全相同的周期性现象基本上不存在,因此,通过这类方法并不能完全消除误差,补偿后残余的相位误差对于高精度的三维测量仍然是不能忽略的。因此,免疫于相位误差编码方法的研究是高精度三维测量首要解决的任务。
发明内容
针对现有方法的不足,本发明提出了一种免疫于相移误差的,高精度的结构光编码方法。
一种基于相位编码技术的光学三维测量方法,由以下步骤组成:(1)、计算机生成光栅编码条纹图像;(2)、使用数字投影设备投射光栅编码条纹图像;(3)、采集光栅条纹图像并利用相移算法进行编码求解和三维测量;其特征在于:(4)、所述的光栅编码条纹图像的每个光栅周期内像素个数等于所采用的相移算法步数;(5)、所述的光栅编码条纹图像的每个光栅周期内灰度值需满足正弦曲线分布,并且像素灰度值为正弦函数的极值点。
本发明的编码方法中,(1)每个光栅编码周期内像素个数与相移步数相等,屏蔽了相位误差产生的像素,增强了对相位误差的抗干扰能力。以四步相移算法为例,相位误差在一个光栅周期内呈现四个周期性的正弦波动,每个正弦波动占距光栅周期的四分之一,这些相位误差是由一个光栅周期的四分之一像素计算得到,为了将这些误差减少到最小,可以将一个光栅周期内四分之一的像素个数减少到1,则一个光栅周期内的四个正弦波动间距减少到零,波动误差呈一个点误差;(2)光栅周期内的像素点灰度为正弦函数的些极值点,使用这些点作为光栅内像素的灰度值进行相位计算,满足了反三角函数的收敛特性,大大减少了非线性灰度反映造成的相位求解误差。根据相位计算的公式(2),相位的求解值是反正切函数的值域,根据三角函数的性质,反正切函数值在±π/2处急剧收敛,定义域的剧烈变化不会大幅度的影响值域的值,因此,相位计算时,采用正弦函数的极值点计算,使得值域的值相对稳定不受定义域剧烈变化的大幅度影响,这也正是相位误差正弦性波动的过零点产生的原因;(3)每个光栅编码周期内的灰度值数量少,而且具有较强的对比性,减少了相机采集的灰度误差.本发明的光栅编码条纹灰度值个数较少,以四步相移为例,仅有255,127.5和0三个灰度值,一个周期内灰度值对比明显,提高了相机采集图像灰度的准确性。
本发明方法根据相移算法中反三角函数的收敛特性,选取三角函数中的极值点值作为编码的灰度值,并将易产生相位误差的像素点进行压缩,从原理上将易产生相位误差的像素点屏蔽,提高了相位误差的抗干扰性,保证了相位计算的正确性和稳定性。相比于传统的相移编码方法具有更高的稳定性、可靠性,提高了编码精度。
附图说明
图1为以四步相移算法为例本发明的光栅编码条纹,灰度图表示。
图2为本发明的光栅编码条纹灰度值曲线图。
图3为相机采集到的本发明的光栅编码灰度图。
图4为使用相移编码计算的相位灰度图。
图5为传统的光栅编码条纹灰度曲线图。
图6为本发明与传统光栅编码误差对比图
具体实施方案
下面对本发明方法结合附图做进一步详细说明。本发明首次使用免疫于相移误差的光栅编码,对多幅编码图像使用传统的相移算法进行求解,既保证了相移算法的高精度,又实现相位误差的消除。
(1)、计算机生成光栅编码条纹图像:通过程序编制生成一系列编码条纹图像,编码图像幅数与使用的相移算法(三步相移算法、四步相移算法、五步相移
算法等)的步数n相等,n幅图像内相同像素坐标点的相位值相差
Figure BSA00000146239300041
图像里每个光栅周期内的像素个数与相移步数n也要相等,并且一个周期内的条纹灰度值满足正弦曲线分布,以四步相移算法为例,用灰度图表示的光栅编码条纹可表示为:Ii(x,y)=I′(x,y)+I″(x,y)cos[(x*2π)/4+i*π/2],i=0,1,...,3,其中Ii(x,y)为相机采集的第i幅光栅条纹图像中像素坐标为(x,y)的光强灰度值,I′(x,y)为背景光强,I″(x,y)为条纹幅值,
Figure BSA00000146239300042
为待测相位值,δi=i*π/2为第i次相移量,绘制出沿x方向的周期性的光栅条纹,取四幅图像中的一幅图像,图像中部分像素的灰度图像表示如图1所示,图像中共有三种灰度值,最黑和最白的灰度图分别是0和255,灰色的图像部分灰度是127.5,沿x方向将灰度值进行曲线表示如图2所示,编码值曲线呈三角形,一个周期内最高点,中间点和最低点在x方向间隔是一个像素,各值分别是图像的灰度值255,127.5和0。Ii(x,y)和δi为已知量,I′(x,y),I″(x,y)和
Figure BSA00000146239300043
三个为未知量,采集n编码图像,使用多n步相移算法(n>=3)即可实现相位
Figure BSA00000146239300044
的求解。
(2)、数字投影设备投射。将本发明设计的编码条纹如图1,通过数字式投影设备进行投射出来,数字投影设备包括,LED投影机、DMD投影机、LCD投影机、LCOS投影机等。调整数字投影设备的分辨率,使投射的图像和分辨率匹配,调整焦距,使投射的图像清晰,另外,根据被测物体表面信息,尽量使投射出的编码条纹具有明显的对比性。
(3)、图像采集和相位解码。根据投射的编码条纹是灰度信息或彩色信息,使用对应的数字式摄像机进行图像采集,首先根据投射的分辨率选择对应匹配的摄像机,然后调整相机的焦距、光圈、曝光时间等参数,保证采集的编码图像清晰,对采集到的多幅编码图像根据公式(2)进行相位编码计算,以四步相移算法为例,对采集到的四幅图像分别为:
Ij(x,y)=I′(x,y)+I″(x,y)cos[(x*2π)/4+j*π/2],j=0,1,...,3,
其中一幅图像显示如图3所示,使用四步相移算法(3)解相位,将相位进行[0255]归一化进行灰度显示如图4所示。
Figure BSA00000146239300051
本发明相比于传统的相位编码方法具有以下优点:(1)每个光栅编码周期内像素个数与相移步数相等,屏蔽了相位误差产生的像素,增强了对相位误差的抗干扰能力;(2)光栅周期内的像素点灰度为正弦函数的些极值点,使用这些点作为光栅内像素的灰度值进行相位计算,满足了反三角函数的收敛特性,大大减少了非线性灰度反映造成的相位求解误差;(3)每个光栅编码周期内的灰度值数量少,而且具有较强的对比性,减少了相机采集的灰度误差。为了验证本发明方法的优越性,与传统的相移编码方法进行了对比实验,图5是传统的光栅编码的灰度曲线图,按以上步骤分别计算相位编码并和理想的相位进行比较,统计其相位误差对比图如图6所示,图中红色线为传统编码的相位误差,蓝色线为本发明的方法计算的相位误差,可以看出使用本发明的方法能得到很好的相位编码,传统的相位编码方法会有正弦性的误差波动,进而会导致三维测量的误差。对于相位编后的包裹相位使用相位展开算法可以实现周期性截断的移除,获得三维测量所用的绝对相位编码,这样就实现了物体表面三维轮廓测量的结构光编码。

Claims (3)

1.一种基于相位编码技术的光学三维测量方法,由以下步骤组成:
(1)、计算机生成光栅编码条纹图像;
(2)、使用数字投影设备投射光栅编码条纹图像;
(3)、采集光栅条纹图像并利用相移算法进行编码求解和三维测量;
其特征在于:
步骤(1)中生成的所述的光栅编码条纹图像的每个光栅周期内像素个数等于所采用的相移算法步数;且所述的光栅编码条纹图像的每个光栅周期内像素灰度需满足正弦或余弦曲线分布,并且像素灰度值为相应正弦或余弦函数的极值点。
2.根据权利要求1所述的一种基于相位编码技术的光学三维测量方法,其特征在于:所述的光栅编码条纹图像为灰度图像。
3.根据权利要求1所述的一种基于相位编码技术的光学三维测量方法,其特征在于:所述的光栅编码条纹图像为彩色图像。
CN201010190225XA 2010-06-02 2010-06-02 基于相位编码技术的光学三维测量方法 Active CN101881605B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010190225XA CN101881605B (zh) 2010-06-02 2010-06-02 基于相位编码技术的光学三维测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010190225XA CN101881605B (zh) 2010-06-02 2010-06-02 基于相位编码技术的光学三维测量方法

Publications (2)

Publication Number Publication Date
CN101881605A CN101881605A (zh) 2010-11-10
CN101881605B true CN101881605B (zh) 2011-11-30

Family

ID=43053676

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010190225XA Active CN101881605B (zh) 2010-06-02 2010-06-02 基于相位编码技术的光学三维测量方法

Country Status (1)

Country Link
CN (1) CN101881605B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014000170A1 (en) * 2012-06-27 2014-01-03 Honeywell International Inc. Encoded information reading terminal with micro-projector
CN103791856A (zh) * 2014-01-21 2014-05-14 南京理工大学 基于四幅光栅条纹图像的相位求解与去包裹方法

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072700B (zh) * 2010-11-18 2013-03-20 华中科技大学 一种基于投影莫尔原理的共面度测量系统
CN102243103A (zh) * 2011-04-13 2011-11-16 西安交通大学 一种物体颜色及三维轮廓的快速测量方法
JP6099115B2 (ja) * 2011-10-26 2017-03-22 学校法人福岡工業大学 三次元表面検査装置および三次元表面検査方法
CN102519396B (zh) * 2011-12-21 2014-11-05 哈尔滨理工大学 三个灰度对称线性编码周期的采样点三维信息获取方法
CN103292734B (zh) * 2013-05-27 2015-10-28 华中科技大学 相位测量系统中伽玛值标定方法
CN103383360B (zh) * 2013-07-29 2016-01-13 重庆理工大学 一种薄带连铸坯表面缺陷正弦光栅相移检测装置及检测方法
CN103759673B (zh) * 2014-01-21 2016-07-06 南京理工大学 基于双频三灰阶正弦光栅条纹投影的时间相位去包裹方法
CN104567730B (zh) * 2015-01-15 2017-11-07 四川大学 一种时空二元编码产生正弦结构光场的方法
CN104697468B (zh) * 2015-03-30 2017-11-21 南昌航空大学 基于倾斜相位编码条纹的光学三维测量方法
CN106017339B (zh) * 2016-06-06 2019-04-23 河北工业大学 一种在非完全约束系统中投影非均匀条纹的三维测量方法
CN106091988B (zh) * 2016-06-17 2018-09-11 太原科技大学 一种高温锻件三维尺寸测量精度的补偿方法
CN106352817A (zh) * 2016-08-31 2017-01-25 李志伟 一种非接触四轮定位仪及其定位方法
CN106441161B (zh) * 2016-11-01 2018-12-11 哈尔滨工程大学 一种基于周期编码的快速相位解缠方法
CN106840039B (zh) * 2017-03-10 2019-03-26 合肥工业大学 一种基于复合光栅投影的快速三维测量方法
CN107339954B (zh) * 2017-05-23 2019-09-13 南昌航空大学 基于周期不同步正弦条纹加相位编码条纹的三维测量方法
CN107421467B (zh) * 2017-07-17 2019-09-20 苏州优纳科技有限公司 结构光三维测量方法及系统
CN108253907B (zh) * 2018-02-01 2020-07-21 深圳市易尚展示股份有限公司 基于希尔伯特变换相位误差校正的三维测量方法和装置
CN109341589B (zh) * 2018-10-17 2020-08-04 深圳市华汉伟业科技有限公司 一种光栅图像投影方法、三维重建方法及三维重建系统
CN110132430B (zh) * 2019-03-29 2021-02-23 黑龙江科技大学 相移法两级编码高精度绝对相位获取方法
CN110455220B (zh) * 2019-08-02 2021-07-06 西安理工大学 一种蝶式定日镜曲面度的测量装置及其测量方法
CN111174730B (zh) * 2020-01-07 2021-07-16 南昌航空大学 一种基于相位编码的快速相位解包裹方法
CN112184788B (zh) * 2020-09-16 2023-11-07 西安邮电大学 一种四步相移的主值相位提取方法
CN113639644A (zh) * 2021-08-12 2021-11-12 武汉维斯克科技有限公司 用于条纹投影轮廓测量的域映射简单伽马计算方法
CN113532330B (zh) * 2021-08-28 2022-10-04 哈尔滨理工大学 一种相位格雷码三维测量方法
CN114279356B (zh) * 2021-12-09 2023-07-25 南京信息工程大学 一种用于三维测量的灰度条纹图案设计方法
CN114754701B (zh) * 2022-02-09 2023-04-25 湖北大学 一种基于Goldstein枝切法的相位展开方法
CN114838674A (zh) * 2022-07-06 2022-08-02 三亚学院 一种用于结构光三维面形测量重建的相位校正方法
CN116295073B (zh) * 2023-02-10 2024-02-09 南京航空航天大学 大型航空复合材料成型模具变形测量装置、方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10321883A1 (de) * 2003-05-07 2004-12-09 Universität Stuttgart Verfahren und Sensor zur hochgenauen optischen Abtastung
CN101089547A (zh) * 2007-07-11 2007-12-19 华中科技大学 一种基于彩色结构光的二维三频解相测量方法
CN101245994A (zh) * 2008-03-17 2008-08-20 南京航空航天大学 物体表面三维轮廓结构光测量系统的标定方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100839167B1 (ko) * 2007-09-18 2008-06-17 주식회사 엔씨비네트웍스 모아레 무늬 발생기를 적용한 위상천이 영사식 3차원형상측정장치 및 그 방법

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10321883A1 (de) * 2003-05-07 2004-12-09 Universität Stuttgart Verfahren und Sensor zur hochgenauen optischen Abtastung
CN101089547A (zh) * 2007-07-11 2007-12-19 华中科技大学 一种基于彩色结构光的二维三频解相测量方法
CN101245994A (zh) * 2008-03-17 2008-08-20 南京航空航天大学 物体表面三维轮廓结构光测量系统的标定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JP特开2009-75104A 2009.04.09
孙军华等.基于相移的彩色结构光编码三维扫描技术.《光学技术》.2008,第34卷(第1期),第122-125页. *
崔海华等.A Three-step system calibration procedure with error compensation for 3D shape measurement.《CHINESE OPTICS LETTERS》.2010,第8卷(第1期),第33-37页. *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014000170A1 (en) * 2012-06-27 2014-01-03 Honeywell International Inc. Encoded information reading terminal with micro-projector
US9773145B2 (en) 2012-06-27 2017-09-26 Honeywell International, Inc. Encoded information reading terminal with micro-projector
CN103791856A (zh) * 2014-01-21 2014-05-14 南京理工大学 基于四幅光栅条纹图像的相位求解与去包裹方法

Also Published As

Publication number Publication date
CN101881605A (zh) 2010-11-10

Similar Documents

Publication Publication Date Title
CN101881605B (zh) 基于相位编码技术的光学三维测量方法
CN107607060B (zh) 一种应用于光栅三维投影测量中的相位误差补偿方法
CN106840036B (zh) 一种适用于快速三维形貌测量的二元结构光优化方法
CN101813461B (zh) 基于复合彩色条纹投影的绝对相位测量方法
CN102322823B (zh) 基于相位级次自编码的光学三维测量方法
CN108955571B (zh) 双频外差与相移编码相结合的三维测量方法
CN104197861B (zh) 基于结构光灰度向量的三维数字成像方法
CN101236066B (zh) 投影光栅的自校正方法
US20220107173A1 (en) Phase-shifting phase measurement error correction method based on pixel tracing of object raster images
CN104567730B (zh) 一种时空二元编码产生正弦结构光场的方法
CN105300317B (zh) 基于正弦和三角波条纹投影的三维测量方法
CN101986098A (zh) 基于三色光栅投影的傅里叶变换三维测量法
CN109141291A (zh) 一种快速相位解包裹算法
CN103292733B (zh) 一种基于相移和三视张量的对应点查找方法
CN108195313A (zh) 一种基于光强响应函数的高动态范围三维测量方法
CN110174079B (zh) 一种基于四步相移编码型面结构光的三维重建方法
CN102261896A (zh) 一种基于相位测量的物体三维形貌测量方法及系统
CN108596008B (zh) 针对三维人脸测量的面部抖动补偿方法
CN110207620A (zh) 一种通过不同频率确定数字光栅投影结构光级数的三维重建方法
CN104390608A (zh) 基于投影栅相位法的结构光三维形状构建方法
CN109523627A (zh) 一种基于泰勒指数表达式型面结构光的三维重建方法
CN109631798A (zh) 一种基于π相移方法的三维面形垂直测量方法
CN102494637A (zh) 三个编码周期的三基色梯形相移三维信息获取方法
CN102853783A (zh) 一种高精度多波长三维测量方法
CN108061517B (zh) 基于莫尔序列光栅的面结构光解相方法

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