CN102589478B - 一种应用于多频率三维测量的全局相位解相方法 - Google Patents
一种应用于多频率三维测量的全局相位解相方法 Download PDFInfo
- Publication number
- CN102589478B CN102589478B CN201210048001.4A CN201210048001A CN102589478B CN 102589478 B CN102589478 B CN 102589478B CN 201210048001 A CN201210048001 A CN 201210048001A CN 102589478 B CN102589478 B CN 102589478B
- Authority
- CN
- China
- Prior art keywords
- phase
- cycle
- formula
- overall
- theta
- 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
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于三维机器视觉领域,涉及一种应用于多频率三维测量的全局相位解相方法。该方法通过投射3种周期分别为P1、P2和P3的正弦或者余弦规律变化的光信息,每种光信息需要经过4-8步相移,通过每个周期的相移图案,分别解算每个周期的相位值、两个周期的合成相位值以及三个周期的最终合成相位值。最后通过合成相位分别计算每个周期的全局相位值,进而通过三维重建获得物体的全貌三维坐标。本发明所设计的全局相位解相方法,优于传统的格雷码相位解相方法,可以有效的解决三维测量中被测物体表面颜色不一致的测量难题,无需喷涂显影剂,更加绿色环保。
Description
技术领域
本发明涉及一种应用于多频率三维测量的全局相位解相方法,更具体的说,本发明涉及一种能够用于高精度三维测量的从多个频率中获取被测物体全局相位的解相方法。
背景技术
光学三维测量方法已广泛应用于工业检测、逆向工程、人体扫描、文物保护、服装鞋帽等多个领域,对自由曲面的检测具有速度快、精度高的优势。按照成像照明方式的不同,光学三维测量技术可分为被动三维测量和主动三维测量两大类。在主动三维测量技术中,结构光三维测量技术发展最为迅速,尤其是相位测量轮廓术(Phase Measuring Profilometry,PMP),也被称为相移测量轮廓术(Phase Shifting Profilometry,PSP),是目前三维测量产品中常用的测量方法。相位测量方法是向被测物体上投射固定周期的按照三角函数(正弦或者余弦)规律变化的光亮度图像,此光亮度图像经过大于3步的均匀相移,最好为4-6步均匀相移,向物体投射4-6次光亮度图像,最终完成一个周期的相位移动。物体上面的每个点,经过相移图像的投射后,在图像中会分别获得几个不同的亮度值。此亮度值经过解相运算,会获得唯一的相位值。由于目前采集到的图像的幅面较大,为了提高相位精度,需要向被测物体投射多个周期的相位图,因此,在一副图像中,相同相位值会出现多次。为了在图像中获得唯一的相位值,格雷码方法是常用的辅助解相方法。目前出现的三维测量产品,普遍采用格雷码加相移的光学投射方法,如德国GOM公司的Atos-I型结构光三维测量***、德国Steinbichler公司的COMMET系列结构光三维测量***、德国Breuckmann公司的optoTOP系列结构光三维测量***、北京天远三维科技有限公司的OKIO-II型三维扫描仪、上海数造科技有限公司的3DSS综合型三维扫描仪、天津世纪动力光电科学仪器有限公司的CPOS三维扫描仪等。由于格雷码的编码方法主要靠图像的二值化来进行编码,因此对于物体表面颜色变化较多的情况,一般需要喷涂显影剂才能实现较好的测量效果。为了解决无法喷涂显影剂的三维测量难题,本发明基于多频率的光学投射方法,设计了一种新的基于多频率的全局相位解相方法。
发明内容
本发明提供一种基于多频率的全局相位解相方法,该方法能够应用于高精度三维测量中,可以弥补格雷码解相方式存在的缺陷。
所述的基于多频率的全局相位解相方法的硬件***包括:
用于投射多种频率光信号的光源投射装置,光源投射装置的分辨率为LR×LC;
用于精度控制、图像采集和数据处理的计算机;
用于采集图像的彩色摄像机,图像分辨率为CR×CC,相机个数为1-2个;
用于放置所述的光源和所述的摄像机的扫描平台;
本发明所设计的全局相位解相方法,对于图像中每一个点(x,y)的全局相位值的计算过程,步骤如下:
步骤1:选取合适的P1、P2和P3的值,P1、P2和P3的取值均在0-LR之间,并且P1和P2的合成周期P12以及P2和P3的合成周期P23,以及P12和P23的合成周期P123,经过公式(1)的运算后,满足P123≥LR;
其中:
M(x,y):为(x,y)点所在的相位周期在全局相位中的编码值,该值的计算公式如公式(3)所示;
全局相位解相方法流程图如图1所示。全局解相后的相位就可以直接用于计算三维坐标。
本发明的有益效果是:通过本发明所介绍的全局相位解相方法,可以解决物体表面颜色不一致的难题,无需喷涂显影剂,即可实现复杂颜色环境下的高精度三维测量,可以避免格雷码全局相位解相方法中存在的缺陷。
附图说明
图1:全局相位解相方法流程图;
图2:相位投射波形及周期示意图;
(a)周期为P1的光源投射波形,P1=16;
(b)周期为P1的光源投射4步相移波形,P1=16;
(c)周期为P1的光亮度函数和相位值得对应关系,P1=16;
图3:相位合成示意图;
(a)周期为P2的光源投射波形,P2=18;
(b)周期为P1与周期为P2的波形图,其中:P1=16,P2=18;
(c)周期为P1和P2的两幅波形进行叠加,可以获得周期为P12的波形,P12=144;
(d)周期为P1、P2和P12的三种周期的相位关系图;
图4:相位合成及流程图。
(a)周期为P2、P3以及P23的波形关系图;
(b)周期为P12、P23以及P123的波形关系图;
(c)相位合成流程图;
图5:全局相位解相示意图
(a)P12周期解相分段图;
(b)P1周期局部相位及全局相位示意图;
(c)P2周期局部相位及全局相位示意图;
具体实施方式
多频率的方法是向被测物体投射至少两种周期的正弦或者余弦函数波,每种周期的函数均经过3-8步相移。通过两种周期的相位函数,可以叠加出一种周期更长的波形。
以一种频率的光信号为例,从光源投射的正弦波形的变化规律如公式(4)所示:
其中:
I(x)为投射光强度;
j: 为周期因子,其值变化为:0~PW
PW:为周期长度;
i: 为步长因子,其值变化为:0~N
N: 为相移的步数
I(x)=sin(θ+δ) 公式(5)
在实际测量中,由于背景光的影响,实际采集到的光亮度Ir的公式如公式(6)所示:
Ir(x)=a+bsin(θ+δ) 公式(6)
以四步相移为例,对于某一个像素点p,通过四次光投射,采集到的图像灰度值Irp1、Irp2、Irp3、Irp4如公式(7)所示。
公式(7)
通过公式(7)可以得到公式(8):
因此相位值θ可通过公式(9)获得:
假设投射的光信息周期为P1,设P1=16,则投射的光信息波形如图1(a)所示。对周期为16的光信息进行四步相移的波形如图1(b)所示。按照公式(6),对每个像素点采集的四幅图像进行解相运算,那么周期为16的光亮度函数和相位值得对应关系如图1(c)所示。从图1(c)可以看出,在144个像素范围内,出现了144/16=9个周期,每个周期均从0-2π变化。对于分辨率为800×600的光源发生装置而言,每次可以投射的周期个数为800/16=50个周期。由于周期个数很多,很难在全场范围内确定唯一相位值。
如果再投射一种周期为P2的波形,则通过P1和P2两幅波形的叠加,可以获得比P1和P2周期都大的波形,将此波形的周期记为P12。则P12的值可以通过公式(10)获得。
设P2=18,则P1和P2的叠加周期为144。周期为18的波形图如图2(a)所示。周期为16与周期为18的波形如图2(b)所示,从图中可以看出,在144的位置,两个周期的波形又重合在一起。通过对周期为P1和P2的两幅波形进行叠加,可以获得周期为144的波形,如图2(c)所示,按照公式(9)的相位计算规律,可以获得不同周期的相位关系图如图2(d)所示。
由于在三维测量中使用的光源投射装置的分辨率为LR×LC,假设LR=800,LC=600,即光源的分辨率为800×600,因此,仅仅通过两个频率的光信号,在全场范围内,会得到800/144个相同的相位周期变化,仍然不能获得在全场唯一的相位变换规律。如果将周期P1和P2的值增大,尽管可以构造在全场范围内相位唯一的正弦函数,但是由于光源的亮度变换范围有限,以及相机的图像采集范围有限等原因,测量精度会下降。为了更好提高测量精度,如果再投射一种周期为P3的波形,则,通过P1与P3,或者P2与P3,都可以再次叠加出一种周期更长的波形,如,可以叠加出P13或者P23。本发明以P23的叠加为例,为了在800个像素内获得唯一的相位值,最佳的取值是使得P12与P23的叠加周期P123经过公式(10)的运算后,P123正好等于800。由于800=25×5×5,很难找到三个值P1、P2、P3,经过两次叠加后获得的P123等于800。为了获得全场范围内唯一的相位值,如果P12与P23的叠加周期大于800,也是可以获得在800范围内唯一的相位值,只不过相位值得变化的最大值不是2π,而是小于2π的数值。假设P3的周期为21,则P2与P3的叠加周期为126。P12与P23的叠加周期为1008,由于1008>800,因此,此种周期的选择是合理的。
综上,P1、P2、P3周期选择的原则是按照公式(11)的运算,最终合成P123的长度需要满足的条件为:P123≥LR。
通过对周期为P2和P3的两幅波形进行叠加,按照公式(9)的相位计算规律,可以得到不同周期的相位关系图如图3(a)所示。按照同样的相位叠加原理,对P12和P23进行叠加,可以获得最终的整体相位P123,P12、P23与P123的关系如图3(b)所示。对于三种频率相位的投射和叠加过程流程图,如图3(c)所示。
周期P1的全局相位值表达式如公式(12)所示。
θG(x,y)=θ1(x,y)+2π×M(x,y) 公式(12)
其中:θG(x,y)为(x,y)点的全局相位值;
θ1(x,y)为(x,y)点根据相移公式计算得到的局部相位值;
M(x,y)为(x,y)点所在的相位周期在全局相位中的编码值。
对于最终合成相位的周期P123的长度中,含有个周期为P12的相位周期。而对于每个P12周期来说,里面含有个P1周期,如图4(a)所示。在三个频率投射过程中,对于每个测量点(x,y)来说,共可以计算得到6个相位值,及 因此,如果以周期为P1的相位进行整体展开的话,那么公式(9)中的M(x,y)可以由公式(13)确定。
公式(13)
以本发明选的P1=16,P2=18,P3=21为例,对周期P1的相位进行全局展开,可以根据公式(12)和公式(13)得到公式(14)。
对周期P2的相位进行全局展开,可以根据公式(12)和公式(13)得到公式(15)。
对周期P3的相位进行全局展开,可以根据公式(12)和公式(13)得到公式(16)。
按照公式(14),对P1周期的相位进行展开的相位图如图4(b)所示。
按照公式(15),对P2周期的相位也进行全局展开,P1展开后相位和P2展开后相位对比图如图4(c)所示。
本发明所设计的全局相位解相方法,对于图像中每一个点(x,y)的全局相位值的计算过程,步骤如下:
步骤1:选取合适的P1、P2和P3的值,P1、P2和P3的取值均在0-LR之间,并且P1和P2的合成周期P12以及P2和P3的合成周期P23,以及P12和P23的合成周期P123,经过下式的运算后,满足P123≥LR;
步骤2:利用光源投射装置,向物体投射周期为P1的相移图像,相移图像的步数应大于3步,最好为4-8步,采集每幅图像,并计算(x,y)在P1周期的相位
其中:
M(x,y):为(x,y)点所在的相位周期在全局相位中的编码值,该值的计算公式如下所示;
全局相位解相方法流程图如图1所示。
本发明与现有的基于格雷码的解相方法的最大区别是:格雷码方法是将图像通过阈值分割,转化为0和1两种信号,无法解决物体表面颜色不一致的情况;而本发明所设计的解相方法,不是仅仅通过0和1两种信号来进行编码,而是通过每个点的相位值来进行编码。因此本发明所设计的方法可以解决物体表面颜色不一致的测量难题,无需喷涂显影剂,即可实现高精度的测量。
本发明与现有的多频率解相方法的区别是:已有的多频率解相方法,仅仅根据最终的合成相位,对被测空间进行相位划分,也仅仅利用最终的合成相位来进行最后的三维计算,由于最终合成相位的计算误差较大,因此,仅仅根据最终合成相位进行三维计算会带来很大的计算误差。而本发明所设计的多频率解相方法,不是仅仅根据最终的合成相位来进行三维计算,而是利用最终合成相位,以及两种频率的合成相位等信息,将三种周期的原始信号进行全局解相。由于三种周期的原始相位计算精度由于合成相位的计算精度,因此本发明所设计的全局相位解相方法,可以最大限度的提高测量精度和可靠性。
综上所述,本发明所述解相方法的优点是:
(1)由于不是靠简单的阈值分割来决定相位编码,因此本发明的全局解相准确度高于格雷码解相方法。
(2)由于不是仅仅根据最终的合成相位来进行解相操作,因此,本发明所设计的全局相位解相方法比传统多频率解相方法具有更高的解相精度。
(3)解决物体表面颜色不一致的测量难题,无需喷涂显影剂等着色材料,测量过程绿色环保,也节约了测量中耗材成本;
以上示意性的对本发明及其实施方式进行了描述,该描述没有局限性,附图中所示的也只是本发明的实施方式之一。所以,如果本领域的普通技术人员受其启示,在不脱离本发明创造宗旨的情况下,采用其它形式的同类部件或其它形式的各部件布局方式,不经创造性的设计出与该技术方案相似的技术方案与实施例,均应属于本发明的保护范围。
Claims (4)
1.一种应用于多频率三维测量的全局相位解相方法,其特征是:对于图像中每一个点(x,y)的全局相位值的计算过程,工作步骤如下:
步骤1:选取P1、P2和P3的值,P1、P2和P3的取值均在0-LR之间,LR由光源投射装置的分辨率LR×LC确定,并且P1和P2的合成周期P12以及P2和P3的合成周期P23,以及P12和P23的合成周期P123,经过公式(1)的运算后,满足P123≥LR;
步骤3:利用光源投射装置,向物体投射周期为P2的相移图像,相移图像的步数应大于3步,采集每幅图像,并计算(x,y)在P2周期的相位
步骤4:利用光源投射装置,向物体投射周期为P3的相移图像,相移图像的步数应大于3步,采集每幅图像,并计算(x,y)在P3周期的相位
其中:
M(x,y):为(x,y)点所在的相位周期在全局相位中的编码值,该值的计算如公式(3)所示;
2.如权利要求1所述的方法,其中步骤2中相移图像的步数为4-8步。
3.如权利要求1所述的方法,其中步骤3中相移图像的步数为4-8步。
4.如权利要求1所述的方法,其中步骤4中相移图像的步数为4-8步。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210048001.4A CN102589478B (zh) | 2012-02-29 | 2012-02-29 | 一种应用于多频率三维测量的全局相位解相方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210048001.4A CN102589478B (zh) | 2012-02-29 | 2012-02-29 | 一种应用于多频率三维测量的全局相位解相方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102589478A CN102589478A (zh) | 2012-07-18 |
CN102589478B true CN102589478B (zh) | 2014-05-21 |
Family
ID=46478470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210048001.4A Expired - Fee Related CN102589478B (zh) | 2012-02-29 | 2012-02-29 | 一种应用于多频率三维测量的全局相位解相方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102589478B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103913128B (zh) * | 2014-03-12 | 2016-09-28 | 哈尔滨工程大学 | 一种基于质量图区域分割的快速相位展开方法 |
CN105806259B (zh) * | 2016-04-29 | 2018-08-10 | 东南大学 | 一种基于二值光栅离焦投影的三维测量方法 |
CN106482665B (zh) * | 2016-09-21 | 2018-05-08 | 大连理工大学 | 一种组合点群高精度三维信息视觉测量方法 |
CN106885533A (zh) * | 2017-03-03 | 2017-06-23 | 哈尔滨理工大学 | 三维傅里叶变换胸腹表面测量方法 |
TWI613617B (zh) * | 2017-03-14 | 2018-02-01 | 立體列印之圖像拼接之方法及其系統 | |
CN107339954B (zh) * | 2017-05-23 | 2019-09-13 | 南昌航空大学 | 基于周期不同步正弦条纹加相位编码条纹的三维测量方法 |
CN108534715A (zh) * | 2018-04-11 | 2018-09-14 | 天津工业大学 | 一种应用于不同波长三维测量的展开相位解相方法 |
CN108613637B (zh) * | 2018-04-13 | 2020-04-07 | 深度创新科技(深圳)有限公司 | 一种基于参考图像的结构光***解相方法及*** |
CN108519064B (zh) * | 2018-04-20 | 2019-12-03 | 天津工业大学 | 一种应用于多频率三维测量的反光抑制方法 |
CN109211140B (zh) * | 2018-08-07 | 2020-12-08 | 广东奥普特科技股份有限公司 | 一种基于数理分析的结构光相位展开方法 |
CN109341536A (zh) * | 2018-09-25 | 2019-02-15 | 深圳市艾视铂智能技术有限公司 | 一种基于双目相机和结构光投影的精密三维视觉测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2175234A1 (en) * | 2008-10-13 | 2010-04-14 | Koh Young Technology Inc. | Method of measuring a three-dimensional shape |
CN102155924A (zh) * | 2010-12-17 | 2011-08-17 | 南京航空航天大学 | 基于绝对相位恢复的四步相移方法 |
CN102243103A (zh) * | 2011-04-13 | 2011-11-16 | 西安交通大学 | 一种物体颜色及三维轮廓的快速测量方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5055191B2 (ja) * | 2008-04-24 | 2012-10-24 | パナソニック株式会社 | 3次元形状計測方法および装置 |
CN101655359A (zh) * | 2009-07-01 | 2010-02-24 | 四川大学 | 用非等步相移正交光栅实现实时三维面形测量的方法 |
-
2012
- 2012-02-29 CN CN201210048001.4A patent/CN102589478B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2175234A1 (en) * | 2008-10-13 | 2010-04-14 | Koh Young Technology Inc. | Method of measuring a three-dimensional shape |
CN102155924A (zh) * | 2010-12-17 | 2011-08-17 | 南京航空航天大学 | 基于绝对相位恢复的四步相移方法 |
CN102243103A (zh) * | 2011-04-13 | 2011-11-16 | 西安交通大学 | 一种物体颜色及三维轮廓的快速测量方法 |
Non-Patent Citations (5)
Title |
---|
JP特开2009-264862A 2009.11.12 |
万文博等.干涉条纹图像处理的相位解包新方法.《应用光学》.2011,第32卷(第1期),第70-74页. |
万文博等.用于光干涉图像处理的相位解包新方法.《万方会议论文数据库》.2010,第1-4页. |
干涉条纹图像处理的相位解包新方法;万文博等;《应用光学》;20110131;第32卷(第1期);第70-74页 * |
用于光干涉图像处理的相位解包新方法;万文博等;《万方会议论文数据库》;20101231;第1-4页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102589478A (zh) | 2012-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102589478B (zh) | 一种应用于多频率三维测量的全局相位解相方法 | |
CN106931910B (zh) | 一种基于多模态复合编码和极线约束的高效三维图像获取方法 | |
CN102853783A (zh) | 一种高精度多波长三维测量方法 | |
CN101881605B (zh) | 基于相位编码技术的光学三维测量方法 | |
CN109489585B (zh) | 基于改进多频条纹结构光的三维测量方法 | |
CN102708566B (zh) | 一种新型的单个摄像机和单个投影光源同步标定方法 | |
CN102607465B (zh) | 一种基于彩色相移条纹二次编码的相位去包裹方法 | |
CN101975558A (zh) | 基于彩色光栅投影的快速三维测量方法 | |
CN102881043A (zh) | 一种应用于明暗反差大场景的高精度三维重建方法 | |
CN102261896A (zh) | 一种基于相位测量的物体三维形貌测量方法及*** | |
CN105300317B (zh) | 基于正弦和三角波条纹投影的三维测量方法 | |
Ding et al. | Absolute phase recovery of three fringe patterns with selected spatial frequencies | |
CN101089547A (zh) | 一种基于彩色结构光的二维三频解相测量方法 | |
CN102589479A (zh) | 三维形貌中心摄动复合光栅投影测量方法及装置 | |
CN110174079B (zh) | 一种基于四步相移编码型面结构光的三维重建方法 | |
CN101109616A (zh) | 一种三频外差相移解相方法 | |
CN101441066A (zh) | 彩色条纹编码的相位去包裹方法 | |
CN108519064B (zh) | 一种应用于多频率三维测量的反光抑制方法 | |
CN103791856A (zh) | 基于四幅光栅条纹图像的相位求解与去包裹方法 | |
CN108955571A (zh) | 双频外差与相移编码相结合的三维测量方法 | |
CN112504165A (zh) | 一种基于双边滤波优化的复合立体相位展开方法 | |
CN108534715A (zh) | 一种应用于不同波长三维测量的展开相位解相方法 | |
CN113465545A (zh) | 一种基于高速led阵列的三维测量***及其测量方法 | |
Yu et al. | Color phase order coding and interleaved phase unwrapping for three-dimensional shape measurement with few projected pattern | |
CN106840039B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140521 Termination date: 20150229 |
|
EXPY | Termination of patent right or utility model |