CN113028989B - 物体的三维信息获取方法及装置 - Google Patents

物体的三维信息获取方法及装置 Download PDF

Info

Publication number
CN113028989B
CN113028989B CN202110244233.6A CN202110244233A CN113028989B CN 113028989 B CN113028989 B CN 113028989B CN 202110244233 A CN202110244233 A CN 202110244233A CN 113028989 B CN113028989 B CN 113028989B
Authority
CN
China
Prior art keywords
phase
gray code
pixel point
fringe pattern
pattern
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
CN202110244233.6A
Other languages
English (en)
Other versions
CN113028989A (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.)
Suzhou Tztek Precision Co ltd
Tztek Technology Co Ltd
Original Assignee
Suzhou Tztek Precision Co ltd
Tztek Technology Co Ltd
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 Suzhou Tztek Precision Co ltd, Tztek Technology Co Ltd filed Critical Suzhou Tztek Precision Co ltd
Priority to CN202110244233.6A priority Critical patent/CN113028989B/zh
Publication of CN113028989A publication Critical patent/CN113028989A/zh
Application granted granted Critical
Publication of CN113028989B publication Critical patent/CN113028989B/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
    • 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/254Projection of a pattern, viewing through a pattern, e.g. moiré

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供一种物体的三维信息获取方法及装置,其中,三维信息获取方法包括:步骤S1,对待测量物体投射相移条纹图和格雷码条纹图;步骤S2,捕捉制后的所述相移条纹图和所述格雷码条纹图;步骤S3,解码得到相位图和格雷码码值图;步骤S4,将所述相位图以及所述格雷码码值图中的像素点分类为第一像素点或第二像素点;步骤S5,确定第一像素点和第二像素点的相位阶数,从而得到相位阶数图,根据所述相位阶数图解析所述相位图,从而得到解包相位图;步骤S6,通过所述解包相位图,获取待测量物体的三维信息。根据本发明的三维信息获取方法,能够准确地获取物体的三维信息,成本低且稳定性高。

Description

物体的三维信息获取方法及装置
技术领域
本发明涉及三维传感领域,具体涉及一种物体的三维信息获取方法及装置。
背景技术
光学三维传感技术是通过特定的光学结构,通过主动发射具有某种特性的光信号,被接收装置捕捉后,对信号进行分析并获取其三维信息的对欧几里得世界感知的重要手段。主要的代表技术有基于光的干涉的干涉仪,基于三角测量原理的,基于飞行时间的TOF等。
干涉仪由于其高昂的研制成本和苛刻的测试环境,纳米级的精度,一般用于实验室,作为精密计量的装置。飞行时间法,主要分为dTOF和iTOF。其中dTOF的核心芯片为SPAD(Single Photon Avalanche Diode单光子雪崩二极管),直接计算光子的到达时间来确定距离的远近。而iTOF主要是利用信号的调制解调,其信号载体为时间本身,通过解调来间接获得飞行时间,异曲同工解决了对于距离的感知。TOF具有感知距离远,结构简单的特点,但其精度较低,一般在mm级,无法用于精密的工业测量领域。
基于三角测量原理的技术,一般分为点,线和面结构光。点激光顾名思义,是通过打出去一个点,对该点返回的信号被相机CCD/CMOS捕捉后,通过投射,相机和空间点组成的三角形进行几何求解。线激光则是通过投射打出一个扇面,该扇面在物体上体现为一条被物体表面调制的轮廓,轮廓线同样被相机捕捉,通过提取激光线中心像素坐标,同样构成相机,投影,空间点的三角关系,得到物体表面该点的空间坐标。而面结构光技术,则是打出去一个面,相对于点和线激光,其具有一次成像范围广,速度快,数据稠密等特点。按照行业分,消费级领域一般使用由VCSEL为代表的散斑结构光,和工业领域以美国德州仪器(TexasInstrument,简称TI)公司的数字投影技术。
但是,目前获取的三维信息中的边缘信息并不准确,导致三维信息精度不高,不能满足客户高精度需求。
发明内容
有鉴于此,本发明提供一种物体的三维信息获取方法及装置,能够准确地获取物体的三维信息,成本低且稳定性高。
为解决上述技术问题,一方面,本发明提供物体的三维信息获取方法,包括如下步骤:
步骤S1,对待测量物体投射相移条纹图和格雷码条纹图,以由待测量物体对所述相移条纹图和所述格雷码条纹图进行调制;
步骤S2,捕捉经所述待测量物体调制后的所述相移条纹图和所述格雷码条纹图;
步骤S3,对捕捉到的经调制后的所述相移条纹图进行解码,得到相位图,对捕捉到的经调制后的所述格雷码条纹图进行解码,得到格雷码码值图,其中,所述相位图中的像素点的位置与所述格雷码码值图中的像素点的位置一一对应;
步骤S4,将所述相位图以及所述格雷码码值图中的像素点分类为第一像素点或第二像素点,所述第一像素点位于所述格雷码码值图的非边缘区域、所述格雷码码值图的跳边沿区域或所述相位图的跳边沿区域,所述第二像素点为所述格雷码码值图或所述相位图中除了所述第一像素点外的其他像素点;
步骤S5,根据所述第一像素点在所述格雷码码值图中的数值确定第一像素点的相位阶数,通过空域的信息对所述第二像素点进行矫正得到第二像素点的相位阶数,根据所述第一像素点和所述第二像素点的相位阶数得到相位阶数图,根据所述相位阶数图解析所述相位图,从而得到解包相位图;
步骤S6,通过所述解包相位图,获取待测量物体的三维信息。
进一步地,所述步骤S5中:所述第二像素点的相位阶数为所述像素第二像素点的前一个像素点的相位阶数,
即所述第一像素点的相位阶数k(x,y)=G(x,y),所述第二像素点的相位阶数k(x,y)=k(x-1,y),
其中,G(x,y)为所述格雷码码值图坐标为(x,y)的像素点的格雷码码值,k(x,y)为坐标为(x,y)的像素点的相位阶数,k(x-1,y)为坐标为(x-1,y)的像素点的相位阶数。
进一步地,所述步骤S3中:
根据|G(x,y)-k(x-1,y)|>1,判断出此像素点位于所述格雷码码值图的跳边沿区域,即为第一像素点;
根据
Figure BDA0002963483960000031
判断出此像素点位于所述相位图的跳边沿区域,即为第一像素点,
其中,
Figure BDA0002963483960000032
为所述相位图中坐标为(x,y)的像素点的相位值,
Figure BDA0002963483960000033
为所述相位图中坐标为(x-1,y)的像素点的相位值,T为相应噪声阈值。
进一步地,T=π。
进一步地,所述步骤S5中:
解包相位图中像素点的相位值
Figure BDA0002963483960000034
进一步地,所述步骤S1中,通过数字光处理方式进行投射所述相移条纹图和所述格雷码条纹图。
进一步地,所述相移条纹图由四步相移法生成,所述格雷码条纹图与所述相移条纹图相对应。
进一步地,所述相移条纹图为正弦相移条纹图。
进一步地,所述步骤S2中对所述相移条纹图进行解码,得到所述相位图包括:
对所述相移条纹图通过反正切函数进行解码,得到所述相位图。
第二方面,本发明提供物体的三维信息获取装置,包括:
投射模块,所述投射模块用于对待测量物体投射相移条纹图和格雷码条纹图,以由待测量物体对所述相移条纹图和所述格雷码条纹图进行调制;
捕捉模块,所述捕捉模块用于捕捉经所述待测量物体调制后的所述相移条纹图和所述格雷码条纹图;
处理模块,所述处理模块用于对捕捉到的经调制后的所述相移条纹图进行解码,得到相位图,对捕捉到的经调制后的所述格雷码条纹图进行解码,得到格雷码码值图,其中,所述相位图中的像素点的位置与所述格雷码码值图中的像素点的位置一一对应;
分类模块,所述分类模块用于将所述相位图以及所述格雷码码值图中的像素点分类为第一像素点或第二像素点,所述第一像素点位于所述格雷码码值图的非边缘区域、格雷码码值图的跳边沿区域或所述相位图的跳边沿区域,所述第二像素点为所述格雷码码值图或所述相位图中除了所述第一像素点外的其他像素点;
解析模块,所述解析模块用于根据所述第一像素点在所述格雷码码值图中的数值确定第一像素点的相位阶数,通过空域的信息对所述第二像素点进行矫正得到第二像素点的相位阶数,根据所述第一像素点和所述第二像素点的相位阶数得到相位阶数图,根据所述相位阶数图解析所述相位图,从而得到解包相位图;
获取模块,所述获取模块用于通过所述解包相位图,获取待测量物体的三维信息。
本发明的上述技术方案至少具有如下有益效果之一:
根据本发明的物体的三维信息获取方法,通过将像素点分为可信度较高的第一像素点和可信度较低的第二像素点,并通过空域信息对第二像素点进行矫正,能够准确地获取解包相位图,从而准确地获取待测量物体的三维信息。
附图说明
图1为根据本发明一实施例的物体的三维信息获取方法的流程图;
图2为根据本发明实施例的物体的三维信息获取装置的示意图;
图3为相位图和格雷码码值图边缘示意图;
图4为截止频率等于信号主频率情况下的方波信号对于投影镜头低通下的数值模拟图;
图5为截止频率等于主频率一半情况下的方波信号对于投影镜头低通下的数值模拟图;
图6为不同频率下的误码率的示意图;
图7为基于相位图边缘的格雷码边缘模糊矫正示意图;
图8为跳边沿导致相位边缘无法被识别到的示意图;
图9为跳边沿的矫正示意图;
图10为方波滤波前后边缘示意图;
图11矫正后的相位图、格雷码码值图及解包相位图。
具体实施方式
下面将结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
下面,结合图1说明根据本发明实施例的物体的三维信息获取方法。
如图1所示,本发明实施例的物体的三维信息获取方法,包括:
步骤S1,对待测量物体投射相移条纹图和格雷码条纹图,以由待测量物体对所述相移条纹图和所述格雷码条纹图进行调制。
相移条纹图和格雷码条纹图结合投射至待测量物体,具有较短的投射时间,高信号鲁棒性,是给物体表面添加纹理信息较为理想的手段。
待测量物体对相移条纹图和格雷码图进行调制,调制后的相移条纹图和格雷码图能够体现待测量物体的三维形貌。
投射相移条纹图和格雷码条纹图通常通过投影镜头向待测量物体投射,经本发明人分析,格雷码条纹图的边缘通过投影镜头,存在边缘模糊和不对称问题,无法保证格雷码条纹图和相位图的边缘完全对齐,从而导致后续的相位解析错误。投影镜头具有低通滤波特性,相移条纹图的信号是单频信号,投影镜头的低通滤波特性造成的失真并不明显,但是格雷码条纹图的信号本质是跳边沿的方波信号,方波信号便于存在除主频外的奇次谐波,而投影镜头的低通滤波特性,导致方波信号的边缘不仅存在模糊,更显示出模糊的不对称性。例如,如图3所示,投射的相位条纹图和格雷码条纹图经过解码后的图3中(a)图的相位图边缘与图3中(b)图的格雷码码值图边缘并不对齐。
步骤S2,捕捉经所述待测量物体调制后的所述相移条纹图和所述格雷码条纹图。
例如,通过相机等捕捉经所述待测量物体调制后的所述相移条纹图和所述格雷码条纹图。由此,能够接收到反应待测量物体的三维形貌的相移条纹图和格雷码条纹图,便于后续的解码。
步骤S3,对捕捉到的经调制后的所述相移条纹图进行解码,得到相位图,对捕捉到的经调制后的所述格雷码条纹图进行解码,得到格雷码码值图,其中,所述相位图中的像素点的位置与所述格雷码码值图中的像素点的位置一一对应。以上的解码过程为现有技术,在此不再详述。相位图中的像素点的位置与格雷码码值图中的像素点的位置一一对应,便于后续的矫正。
步骤S4,将所述相位图以及所述格雷码码值图中的像素点分类为第一像素点或第二像素点,所述第一像素点位于所述格雷码码值图的非边缘区域、所述格雷码码值图的跳边沿区域或所述相位图的跳边沿区域,所述第二像素点为所述格雷码码值图或所述相位图中除了所述第一像素点外的其他像素点。
将像素点分类为第一像素点和第二像素点,第一像素点即可信度较高的区域的像素点,第二像素点即可信度较低的区域的像素点,便于针对性地对容易解码错误的区域进行矫正,避免对所有像素点进行矫正而造成无法矫正正确的情况。
步骤S5,根据所述第一像素点在所述格雷码码值图中的数值确定第一像素点的相位阶数,通过空域的信息对所述第二像素点进行矫正得到第二像素点的相位阶数,根据所述第一像素点和所述第二像素点的相位阶数得到相位阶数图,根据所述相位阶数图解析所述相位图,从而得到解包相位图。
其中,空域信息为通过某种可信度较高的种子点,通过规划路径解码,即通过一些可信度高的像素点去矫正可信度低的像素点。除了空域信息矫正,还可以利用时域的信息进行矫正,即投射一系列相位阶数具有对应关系的格雷码条纹,但是通过下文的验证,通过时域的信息无法解决格雷码条纹图边缘存在模糊和不对称的情况。
也就是说,第一像素点的相位阶数为格雷码码值图中的数值,第二像素点由于格雷码码值图中的数值并不准确,通过空域信息进行矫正,从而得到第二像素点准确的相位阶数值。通过准确的第一像素点和第二像素点的相位阶数值确定相位阶数图,通过相位阶数图解析相位图,从而得到准确的解包相位图。
步骤S6,通过所述解包相位图,获取待测量物体的三维信息。通过解包相位图获取待测量物体的三维信息为现有技术,在此不再详述。
以上的物体的三维信息获取方法,通过对将像素点分为可信度较高的第一像素点和可信度较低的第二像素点,并通过空域信息对第二像素点进行矫正。由此,能够准确地获取解包相位图,从而准确地获取待测量物体的三维信息。
根据本发明一些实施例,所述步骤S5中:所述第二像素点的相位阶数为所述像素第二像素点的前一个像素点的相位阶数,即所述第一像素点的相位阶数k(x,y)=G(x,y),所述第二像素点的相位阶数k(x,y)=k(x-1,y),
其中,G(x,y)为所述格雷码码值图坐标为(x,y)的像素点的格雷码码值,k(x,y)为坐标为(x,y)的像素点的相位阶数,k(x-1,y)为坐标为(x-1,y)的像素点的相位阶数。
也就是说,通过空域信息对第二像素点进行矫正是通过前一个像素点的相位阶数作为此第二像素点的相位阶数。由此,能够准确获取第二像素点的相位阶数。
进一步地,所述步骤S3中:根据|G(x,y)-k(x-1,y)|>1,判断出此像素点位于所述格雷码码值图的跳边沿区域,即为第一像素点;根据
Figure BDA0002963483960000071
Figure BDA0002963483960000072
判断出此像素点位于所述相位图的跳边沿区域,即为第一像素点,
其中,
Figure BDA0002963483960000073
为所述相位图中坐标为(x,y)的像素点的相位值,
Figure BDA0002963483960000074
为所述相位图中坐标为(x-1,y)的像素点的相位值,T为相位噪声阈值。
由此,能够简便且准确地获知此像素点为第一像素点,便于将像素点准确地分类为第一像素点或第二像素点。
优选地,T=π。相位噪声阈值T(对相位噪声具有的容忍性),可以根据三维测量的精度进行选择,T小于2π,例如为1.5π,0.8π等。
进一步地,所述步骤S5中:解包相位图中像素点的相位值
Figure BDA0002963483960000075
Figure BDA0002963483960000076
也就是说,利用相位阶数图中像素点的相位阶数值k(x,y)解析相位图中的像素点的相位值
Figure BDA0002963483960000077
从而得到解包相位图中的像素点的相位值Φ(x,y)。由此,能够准确地解析出解包相位图。
根据本发明一些实施例,所述步骤S1中,通过数字光处理(DLP)方式进行投射所述相移条纹图和所述格雷码条纹图。
以DMD(数字微镜器件)为代表的DLP数字投影方式具有一次成像范围广、速度快、稳定性高及数据稠密等特点,能够以较好的方式投射相移条纹图和格雷码条纹图。
进一步地,所述相移条纹图由四步相移法生成,所述格雷码条纹图与所述相移条纹图相对应。
四歩相移法具有一定滤波作用,对相位噪声具有一定的抑制作用,格雷码条纹图与相移条纹图相对应,便于解码后的相位图和格雷码码值图中的像素点一一对应。
进一步地,所述相移条纹图为正弦相移条纹图。
进一步地,所述步骤S2中对所述相移条纹图进行解码,得到所述相位图包括:对所述相移条纹图通过反正切函数进行解码,得到所述相位图。
由此,能够准确地获取相位图,便于准确地获取待测量物体的三维信息。
下面,结合图2,说明根据本发明实施例的物体的三维信息获取装置1000。
如图2所示,本发明实施例的物体的三维信息获取装置1000包括:投射模块1001、捕捉模块1002、处理模块1003、分类模块1004、解析模块1005及获取模块1006。
投射模块1001用于对待测量物体投射相移条纹图和格雷码条纹图,以由待测量物体对所述相移条纹图和所述格雷码条纹图进行调制。投影模块可以包括投影机等。
捕捉模块1002用于捕捉经所述待测量物体调制后的所述相移条纹图和所述格雷码条纹图。捕捉模块可以包括相机等。
处理模块1003用于对捕捉到的经调制后的所述相移条纹图进行解码,得到相位图,对捕捉到的经调制后的所述格雷码条纹图进行解码,得到格雷码码值图,其中,所述相位图中的像素点的位置与所述格雷码码值图中的像素点的位置一一对应。
分类模块1004用于将所述相位图以及所述格雷码码值图中的像素点分类为第一像素点或第二像素点,所述第一像素点位于所述格雷码码值图的非边缘区域、格雷码码值图的跳边沿区域或所述相位图的跳边沿区域,所述第二像素点为所述格雷码码值图或所述相位图中除了所述第一像素点外的其他像素点。
解析模块1005用于根据所述第一像素点在所述格雷码码值图中的数值确定第一像素点的相位阶数,通过空域的信息对所述第二像素点进行矫正得到第二像素点的相位阶数,根据所述第一像素点和所述第二像素点的相位阶数得到相位阶数图,根据所述相位阶数图解析所述相位图,从而得到解包相位图。
获取模块1006用于通过所述解包相位图,获取待测量物体的三维信息。
由此,能够准确获取待测量物体的三维信息。
以下分析通过时域的信息为什么无法解决格雷码条纹图边缘模糊和不对称的问题。
假设投影镜头的传递函数呈现高斯滤波器特性,那么在一维情况下,不妨假设
Figure BDA0002963483960000091
为零偏置高斯分布,且其最高点为h(0)=1,如果令-3db为截止频率,且频率分辨率为fres,那么其截止频率为
Figure BDA0002963483960000092
如果信号频率可给出初始值例如fs,故可根据信号频率,来选择合适低通特性的投影镜头,做定制化设计。如下所示
Figure BDA0002963483960000093
故本文根据式4-3来生成不同截止频率下的高斯滤波器。
如图4所示,当截止频率等于信号主频率情况下的方波信号对投影镜头的低通下的数值模拟。图4中(a)图为原始方波信号,图4中(b)图为方波信号的频率和低通滤波器,图4中(c)图为经过低通滤波后的方波信号,图4中(d)图为二值化后的方波。从图4中可以看出,当当截止频率等于信号主频率的时候,即使产生模糊,变成第二行左所示图像,原始信号也能被较好还原。其主要原因在于,该情况下,方波的边缘仅仅是变圆滑,但依旧是关于直流分量对称的,而该直流分量是用来做二值化的阈值。
如图5所示,当截止频率等于信号主频率一半情况下的方波信号对投影镜头的低通下的数值模拟。图5中(a)图为原始方波信号,图5中(b)图为方波信号的频率和低通滤波器,图5中(c)图为经过低通滤波后的方波信号,图5中(d)图为二值化后的方波。从图5中可以看出,截至频率低于主频,信号本身存在较大的跳动,如图5中(c)图所示,其并非关于直流分量对称,故使用直流分量作为阈值后的方波,如图5中(d)图所示,和原始信号相差甚远,无法还原原始信号。
如图6所示,不同频率的误码率示意图。图6中(a)图为差分信号经过同一个高斯低通滤波后,且截至频率为信号主频的一半,图6中(b)图为单方向格雷码下的误码率,图6中(c)图为差分格雷码下的误码率,图6中(d)图为投影镜头滤波下的差分信号示意图,其中,单方向指只投放0°信号,差分指相差180°的两路信号。从图6中(a)图和(b)图可以看出,即使主频被截取,其差分信号依旧是关于直流分量对称,从图6中(c)图和(d)图可以看出,即使使用了差分信号,其错码率依旧没有改善,原因示意在图6中(d)图,差分信号并没有改变其误差带。
以上模拟结果表明,在低通滤波下,格雷码边缘存在模糊和不对称性,无法使用时域的信息(时域上通过增加差分信号的方式)去矫正,故有本发明利用空域的信息来矫正。
以下具体说明本发明的矫正方法。
本发明主要利用相位图和格雷码码值图。
相位图:相位图上每个像素的相位是通过四象限反三角函数arctan2计算得到,其存在2π的断差,体现在相位图上即为边缘信息。故对于相位图而言,其边缘信息非常明确。
格雷码解码图:由于格雷码本身是一种二进制码,那么其在非边缘地区的鲁棒性是非常强的。原因在于,表征格雷码的是最亮和最暗,故用合适的阈值(一般会投射一张全白和全黑的图,其均值作为每个点的阈值),如果物体表面反射率不是很低的情况下,是非常容易得到像素点对应的格雷码值是0还是1。但是,在边缘地区,由于投影镜头的低通特性,或者说是调制传递过程中,会滤除一部分的高频成分,导致其边缘变的圆滑,且存在基于其直流分量的不对称性。导致其边缘无法准确得到是0还是1。
那么本文提到的解决方案,刚好是取长补短,有机结合两者的优势和劣势,来矫正格雷码边缘的模糊不对称性问题。
首先,定义某点的差分信号为后差分
Figure BDA0002963483960000111
假设格雷码在该点的码值对应二进制为G(x,y),那么做以下操作
Figure BDA0002963483960000112
其中k(x,y)为相位的阶数,T表示阈值(小于2π),为了对相位噪声具有一定的容忍性,可考虑取较大的数值,例如π。
从图7中(a)图看出,当格雷码边缘出现在非相位边缘时,会被逐个纠正,直到遇到相位边缘。由于格雷码在非边缘的高鲁棒性,故其很少出现00101000111等码值连续跳动的情况;可认为格雷码边缘结果是比较准确的。
该方法的优点:对相位图误差不敏感。
1)如果相位错码出现在格雷码稳定区域,但我们只关心格雷码变换区域,故不在考虑内。
2)边缘正确。如图7中(a)图所示,解码完全正确。
3)边缘错误。如图7中(b)图所示,如果相位由于噪声,导致其边缘提取错误,例如提前一个像素,那么会导致格雷码边缘提前一个像素被锁定,但不会影响后续的码值。
从图8中可以看出,在跳边沿(Jump edge)处,很可能通过差分式无法识别,因为相位相差在阈值范围内。导致出现图8所示情况,后续码值连续解错。
那么为了解决这个问题。可通过加强限制条件为
Figure BDA0002963483960000113
由于Jump edge,本身相位就不准,所以该点最后并不会被用到。但为了保证k(x,y)的连续性,保证后续像素的正确性,故不作处理,最终如图9中(a)图和(b)图所示。
其次,在像素分辨率极高的情况下,很少会遇到这种边界。再者,跳变后,落在边缘的可能性极小,故保持格雷码值不变,即大概率保证了其解码的正确性,较为合理。
如图10所示,格雷码边缘落后相位边缘下的方波滤波前后边缘示意图。在特定的系统下,判断该现象是否有可能发生,继续对方波进行模拟。
对不同的截止频率进行模拟,如图10中(a)图所示,曲线表示截止频率从0.4倍的信号频率到0.8倍的信号频率下的方波滤波后的形态。图10中(b)图和(c)图分别为图10中(a)图的下降沿和上升沿的局部放大图。
fc/fs=0.4的线表示二值的阈值。故可看到,无论是上升沿还是下降沿,曲线穿过阈值的时刻先于方波边缘,即格雷码边缘永远先于相位图边缘。
所以,不需要讨论当相位图边缘先于格雷码边缘的情况,在物理上不存在这种情况。
以下验证矫正的效果。该实现拍摄对象是拍摄了一张白纸,投射一组正弦相移图和一组对应的格雷码条纹图,通过解码得到一张相位图和一张格雷码码值图。从图11中(a)图的相位图、图11中的(b)图的格雷码码值图及图11中(c)图的解包相位图可以看出,格雷码码值图的边缘和相位图的边缘严格对齐,解包相位图的相位连续平滑,无毛刺。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种物体的三维信息获取方法,其特征在于,包括如下步骤:
步骤S1,对待测量物体投射相移条纹图和格雷码条纹图,以由待测量物体对所述相移条纹图和所述格雷码条纹图进行调制;
步骤S2,捕捉经所述待测量物体调制后的所述相移条纹图和所述格雷码条纹图;
步骤S3,对捕捉到的经调制后的所述相移条纹图进行解码,得到相位图,对捕捉到的经调制后的所述格雷码条纹图进行解码,得到格雷码码值图,其中,所述相位图中的像素点的位置与所述格雷码码值图中的像素点的位置一一对应;
步骤S4,将所述相位图以及所述格雷码码值图中的像素点分类为第一像素点或第二像素点,所述第一像素点位于所述格雷码码值图的非边缘区域、所述格雷码码值图的跳边沿区域或所述相位图的跳边沿区域,所述第二像素点为所述格雷码码值图或所述相位图中除了所述第一像素点外的其他像素点;
步骤S5,根据所述第一像素点在所述格雷码码值图中的数值确定第一像素点的相位阶数,通过空域的信息对所述第二像素点进行矫正得到第二像素点的相位阶数,根据所述第一像素点和所述第二像素点的相位阶数得到相位阶数图,根据所述相位阶数图解析所述相位图,从而得到解包相位图;
步骤S6,通过所述解包相位图,获取待测量物体的三维信息。
2.根据权利要求1所述的物体的三维信息获取方法,其特征在于,所述步骤S5中:所述第二像素点的相位阶数为所述第二像素点的前一个像素点的相位阶数,
即所述第一像素点的相位阶数k(x,y)=G(x,y),所述第二像素点的相位阶数k(x,y)=k(x-1,y),
其中,G(x,y)为所述格雷码码值图坐标为(x,y)的像素点的格雷码码值,k(x,y)为坐标为(x,y)的像素点的相位阶数,k(x-1,y)为坐标为(x-1,y)的像素点的相位阶数。
3.根据权利要求2所述的物体的三维信息获取方法,其特征在于,所述步骤S3中:
根据|G(x,y)-k(x-1,y)|>1,判断出此像素点位于所述格雷码码值图的跳边沿区域,即为第一像素点;
根据
Figure FDA0003707417380000021
判断出此像素点位于所述相位图的跳边沿区域,即为第一像素点,
其中,
Figure FDA0003707417380000022
为所述相位图中坐标为(x,y)的像素点的相位值,
Figure FDA0003707417380000023
为所述相位图中坐标为(x-1,y)的像素点的相位值,T为相应噪声阈值。
4.根据权利要求3所述的物体的三维信息获取方法,其特征在于,T=π。
5.根据权利要求4所述的物体的三维信息获取方法,其特征在于,所述步骤S5中:
解包相位图中像素点的相位值
Figure FDA0003707417380000024
6.根据权利要求1所述的物体的三维信息获取方法,其特征在于,所述步骤S1中,通过数字光处理方式进行投射所述相移条纹图和所述格雷码条纹图。
7.根据权利要求6所述的物体的三维信息获取方法,其特征在于,
所述相移条纹图由四步相移法生成,所述格雷码条纹图与所述相移条纹图相对应。
8.根据权利要求7所述的物体的三维信息获取方法,其特征在于,所述相移条纹图为正弦相移条纹图。
9.根据权利要求8所述的物体的三维信息获取方法,其特征在于,所述步骤S2中对所述相移条纹图进行解码,得到所述相位图包括:
对所述相移条纹图通过反正切函数进行解码,得到所述相位图。
10.一种物体的三维信息获取装置,其特征在于,包括:
投射模块,所述投射模块用于对待测量物体投射相移条纹图和格雷码条纹图,以由待测量物体对所述相移条纹图和所述格雷码条纹图进行调制;
捕捉模块,所述捕捉模块用于捕捉经所述待测量物体调制后的所述相移条纹图和所述格雷码条纹图;
处理模块,所述处理模块用于对捕捉到的经调制后的所述相移条纹图进行解码,得到相位图,对捕捉到的经调制后的所述格雷码条纹图进行解码,得到格雷码码值图,其中,所述相位图中的像素点的位置与所述格雷码码值图中的像素点的位置一一对应;
分类模块,所述分类模块用于将所述相位图以及所述格雷码码值图中的像素点分类为第一像素点或第二像素点,所述第一像素点位于所述格雷码码值图的非边缘区域、格雷码码值图的跳边沿区域或所述相位图的跳边沿区域,所述第二像素点为所述格雷码码值图或所述相位图中除了所述第一像素点外的其他像素点;
解析模块,所述解析模块用于根据所述第一像素点在所述格雷码码值图中的数值确定第一像素点的相位阶数,通过空域的信息对所述第二像素点进行矫正得到第二像素点的相位阶数,根据所述第一像素点和所述第二像素点的相位阶数得到相位阶数图,根据所述相位阶数图解析所述相位图,从而得到解包相位图;
获取模块,所述获取模块用于通过所述解包相位图,获取待测量物体的三维信息。
CN202110244233.6A 2021-03-05 2021-03-05 物体的三维信息获取方法及装置 Active CN113028989B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110244233.6A CN113028989B (zh) 2021-03-05 2021-03-05 物体的三维信息获取方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110244233.6A CN113028989B (zh) 2021-03-05 2021-03-05 物体的三维信息获取方法及装置

Publications (2)

Publication Number Publication Date
CN113028989A CN113028989A (zh) 2021-06-25
CN113028989B true CN113028989B (zh) 2022-08-16

Family

ID=76468004

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110244233.6A Active CN113028989B (zh) 2021-03-05 2021-03-05 物体的三维信息获取方法及装置

Country Status (1)

Country Link
CN (1) CN113028989B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114719781B (zh) 2022-06-08 2022-09-16 广东工业大学 一种三维测量方法及相关装置
CN115830154B (zh) * 2023-02-22 2023-05-02 南京信息工程大学 一种基于二倍角相位编码的解包裹方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100443854C (zh) * 2006-09-15 2008-12-17 东南大学 三维扫描系统中基于格雷码的相位展开方法
US20110080471A1 (en) * 2009-10-06 2011-04-07 Iowa State University Research Foundation, Inc. Hybrid method for 3D shape measurement
CN102175179A (zh) * 2011-02-23 2011-09-07 东南大学 一种人体表面轮廓三维重建的方法与装置
CN105890546A (zh) * 2016-04-22 2016-08-24 无锡信捷电气股份有限公司 基于正交格雷码和线移相结合的结构光三维测量方法
EP3516328B1 (en) * 2016-09-21 2023-05-03 Philip M. Johnson Non-contact coordinate measuring machine using hybrid cyclic binary code structured light
CA2945256C (en) * 2016-10-13 2023-09-05 Lmi Technologies Inc. Fringe projection for in-line inspection
CN107607060B (zh) * 2017-08-24 2019-09-10 东南大学 一种应用于光栅三维投影测量中的相位误差补偿方法
CN108955573B (zh) * 2018-06-05 2020-03-06 天津大学 一种编码结构光的阶次无损校正方法
JP2020041822A (ja) * 2018-09-06 2020-03-19 株式会社Yoods 3次元計測システム、3次元計測カメラ、3次元計測方法及びプログラム
CN109186476B (zh) * 2018-10-26 2021-01-26 广东工业大学 一种彩色结构光三维测量方法、装置、设备及存储介质
CN109740402B (zh) * 2019-01-02 2022-04-22 苏州天准科技股份有限公司 基于频谱信息的二维码方向与位置获取方法
JP2020139869A (ja) * 2019-02-28 2020-09-03 キヤノン株式会社 計測装置、算出方法、システム及びプログラム
CN110068273B (zh) * 2019-05-21 2020-08-14 苏州天准软件有限公司 一种基于3d模型的激光融合标定方法
US10659056B1 (en) * 2019-06-13 2020-05-19 Omnivision Technologies, Inc. Gray code counting signal distribution system
CN111724362A (zh) * 2020-06-12 2020-09-29 七海测量技术(深圳)有限公司 一种基于视觉条纹技术的高反面外观检测方法

Also Published As

Publication number Publication date
CN113028989A (zh) 2021-06-25

Similar Documents

Publication Publication Date Title
JP6162225B2 (ja) 能動奥行き検知のためのアフィン不変空間マスクの送信
CN113028989B (zh) 物体的三维信息获取方法及装置
US9557167B2 (en) Three-dimensional-shape measurement apparatus, three-dimensional-shape measurement method, and non-transitory computer-readable storage medium
US12013228B2 (en) Phase-shifting phase measurement error correction method based on pixel tracing of object raster images
US10670829B2 (en) Imaging device
US9948920B2 (en) Systems and methods for error correction in structured light
US20170329012A1 (en) Optoelectronic modules for distance measurements and/or multi-dimensional imaging
EP3271893A1 (en) Systems and methods for enhanced depth map retrieval for moving objects using active sensing technology
CN103069250A (zh) 三维测量设备、三维测量方法和计算机程序
CN107408306B (zh) 一种产生对象的深度图信息的方法、设备及可读介质
CN111308482B (zh) 基于编码调制图像的滤波连续波飞行时间测量
US20160284102A1 (en) Distance measurement apparatus, distance measurement method, and storage medium
US20220066004A1 (en) Distortion determination apparatus and method of determining a distortion
CN112082508A (zh) 反光物体的三维测量方法、装置及终端设备
CN112824934B (zh) 基于调制光场的tof多径干扰去除方法、系统、设备及介质
WO2022066583A1 (en) Decoding an image for active depth sensing to account for optical distortions
Shi et al. Depth sensing with coding-free pattern based on topological constraint
CN116520348A (zh) 基于调制光场的深度成像系统、方法、设备及介质
KR20000053779A (ko) 2차원 격자무늬를 이용한 3차원 형상측정시스템
CN113327317B (zh) 三维点云图获取方法、装置、电子设备和存储介质
Gupta Gray code composite pattern structured light illumination
He Temporal phase unwrapping towards high-accuracy and high-speed fringe projection profilometry
FI98959C (fi) Menetelmä ja laite pinnan korkeustiedon määrittämiseksi pinnan tasomaisuuden mittauksessa

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