CN111023995A - 一种基于随机两帧相移条纹图的三维测量方法 - Google Patents

一种基于随机两帧相移条纹图的三维测量方法 Download PDF

Info

Publication number
CN111023995A
CN111023995A CN201911128510.6A CN201911128510A CN111023995A CN 111023995 A CN111023995 A CN 111023995A CN 201911128510 A CN201911128510 A CN 201911128510A CN 111023995 A CN111023995 A CN 111023995A
Authority
CN
China
Prior art keywords
fringe pattern
deformed
correction result
pattern
phase shift
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
Application number
CN201911128510.6A
Other languages
English (en)
Other versions
CN111023995B (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN201911128510.6A priority Critical patent/CN111023995B/zh
Publication of CN111023995A publication Critical patent/CN111023995A/zh
Application granted granted Critical
Publication of CN111023995B publication Critical patent/CN111023995B/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/254Projection of a pattern, viewing through a pattern, e.g. moiré
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/10Image enhancement or restoration using non-spatial domain filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/04Indexing scheme for image data processing or generation, in general involving 3D image data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20056Discrete and fast Fourier transform, [DFT, FFT]

Landscapes

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

Abstract

本发明涉及一种基于随机两帧相移条纹图的三维测量方法,包括:向待测物体依次投影第一条纹图和第二条纹图,得到第一变形条纹图和第二变形条纹图;分别对第一变形条纹图和第二变形条纹图进行背景滤除和调制度幅值矫正,得到第一矫正结果和第二矫正结果;利用利萨如椭圆拟合方法对第一矫正结果和第二矫正结果进行处理,得到待测物体的相位信息;根据相位信息,得到待测物体的三维形貌信息。本发明的方法利萨如椭圆拟合技术提取物体的相位信息应用到条纹投影系统,仅采用两帧具有不均匀背景强度和调制度的条纹图重建物体的三维形貌,计算更简便,重建的三维形貌具有较高的鲁棒性和精确度。

Description

一种基于随机两帧相移条纹图的三维测量方法
技术领域
本发明属于基于条纹投影的物体三维测量技术领域,具体涉及一种基于随机两帧相移条纹图的三维测量方法。
背景技术
在三维物体面形测量领域中,基于条纹投影技术的三维物体测量技术在各个领域已被广泛应用,例如,工业检测,人脸检测,医疗整形等,由于其低成本,高精度,高可靠和非接触性成为三维物体测量技术中最具商业化和最有效的方法之一。该技术是利用投影仪将条纹投影到待测物体表面,被测物体形貌的变化导致条纹发生形变,对这种条纹形变进行分析并提取变形条纹图的相位信息进而获取物体的高度信息。基于条纹投影技术的相位提取技术可以分为多帧和单帧的方法。
多帧方法,例如相移轮廓术,该方法至少需要三帧条纹图提取物体的相位,投影仪投影一系列条纹图至待测物体表面,待测物体形貌的变化导致条纹发生形变,由相机捕获变形的条纹图,最后提取变形条纹图的相位信息并根据相位-高度映射关系式获取物体的三维形貌。但是,由于需要多帧条纹图提取物体的相位,难以实现高速三维测量,而且当测量运动的物体时,所捕获的运动物体条纹图之间的相移量不均匀,会导致相位的测量误差,另外,当物体运动时,帧间的时间间隔不可忽略,因此将会导致无法捕获帧间物体的运动信息,因此不能实现实时三维测量。
单帧方法,例如傅里叶轮廓术,该方法仅需一帧条纹图提取物体的相位信息,这种方法在快速测量时具有优势,但通常对噪音和表面纹理变化非常敏感以及由于频谱丢失导致该方法的重建精度并不高。
因此,提出一种可以精确地提取物体的相位信息以及实现快速三维测量的方法具有重要的应用价值和前景。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于随机两帧相移条纹图的三维测量方法。本发明要解决的技术问题通过以下技术方案实现:
本发明提供了一种基于随机两帧相移条纹图的三维测量方法,其特征在于,包括:
向待测物体依次投影第一条纹图和第二条纹图,得到第一变形条纹图和第二变形条纹图;
分别对所述第一变形条纹图和所述第二变形条纹图进行背景滤除和调制度幅值矫正,得到第一矫正结果和第二矫正结果;
利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到所述待测物体的相位信息;
根据所述相位信息,得到所述待测物体的三维形貌信息。
在本发明的一个实施例中,所述第一条纹图和所述第二条纹图之间的相移量为任意相移量。
在本发明的一个实施例中,所述第一变形条纹图的强度为,
Figure BDA0002277616100000021
所述第二变形条纹图的强度为,
Figure BDA0002277616100000037
其中,(x,y)表示图像的像素点坐标,A1(x,y)表示所述第一变形条纹图的背景强度,B1(x,y)表示所述第一变形条纹图的调制度幅值,A2(x,y)表示所述第二变形条纹图的背景强度,B2(x,y)表示所述第二变形条纹图的调制度幅值,
Figure BDA0002277616100000036
表示待测物体的相位信息,δ表示第一变形条纹图与第二变形条纹图的相移量。
在本发明的一个实施例中,分别对所述第一变形条纹图和所述第二变形条纹图进行背景滤除和调制度幅值矫正,得到第一矫正结果和第二矫正结果,包括:
利用高通滤波器分别将所述第一变形条纹图和所述第二变形条纹图的背景滤除,得到第一滤除背景图和第二滤除背景图,其强度表达式如下,
Figure BDA0002277616100000031
Figure BDA0002277616100000032
根据归一化公式分别对所述第一滤除背景图和所述第二滤除背景图进行归一化处理,得到所述第一矫正结果和所述第二矫正结果,其中,
所述归一化公式如下:
Figure BDA0002277616100000033
其中,unorm表示归一化后的向量,<>表示内积,||||表示求范数;
所述第一矫正结果为,
Figure BDA0002277616100000034
所述第二矫正结果为,
Figure BDA0002277616100000035
其中,k=1,2,…,K表示所述变形条纹图中的像素点。
在本发明的一个实施例中,利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到所述待测物体的相位信息,包括:
利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到利萨如椭圆表达式;
根据所述利萨如椭圆表达式,得到所述第一变形条纹图和所述第二变形条纹图之间的相移量;
根据所述相移量,得到所述测物体的相位信息。
在本发明的一个实施例中,利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到利萨如椭圆表达式,包括:
将所述第一矫正结果和所述第二矫正结果分别相加和相减,简化后得到下式,
Figure BDA0002277616100000041
Figure BDA0002277616100000042
根据上式,得到利萨如椭圆表达式,
Figure BDA0002277616100000043
其中,x0=0,y0=0,(x0,y0)表示利萨如椭圆的圆心,
Figure BDA0002277616100000044
Figure BDA0002277616100000045
在本发明的一个实施例中,根据所述利萨如椭圆表达式,得到所述第一变形条纹图和所述第二变形条纹图之间的相移量,包括:
对所述利萨如椭圆表达式进行变形转换,得到二次方程,
Figure BDA0002277616100000051
其中,
Figure BDA0002277616100000052
采用最小二乘法对所述二次方程进行求解,得到系数a和b,
Figure BDA0002277616100000053
根据所述系数a和b,得到所述第一变形条纹图和所述第二变形条纹图之间的相移量,
Figure BDA0002277616100000054
在本发明的一个实施例中,根据所述相移量,得到所述测物体的相位信息,包括:
根据所述相移量,按照下式计算得到所述测物体的相位信息,
Figure BDA0002277616100000055
在本发明的一个实施例中,根据所述相位信息,得到所述待测物体的三维形貌信息,包括:
根据所述相位信息利用相位-高度映射关系,得到所述待测物体的三维形貌信息,
Figure BDA0002277616100000056
其中,l表示在测量系统中相机到待测物体参考平面的距离,d表示在测量系统中相机与投影仪之间的距离,f表示条纹图的空间频率。
与现有技术相比,本发明的有益效果在于:
1、本发明的基于随机两帧相移条纹图的三维测量方法,将利萨如椭圆拟合技术提取物体的相位信息应用到条纹投影系统,采用两帧具有不均匀背景强度和调制度的条纹图重建物体的三维形貌,与单帧傅里叶方法相比,本发明提高了物体的重建精度,解决了傅里叶轮廓术由于频谱丢失导致的重建精度不高的问题,与多帧相移方法相比,本发明所提方法仅用两帧条纹提取物体的相位信息,减少了条纹帧数,计算更简便。
2、本发明的基于随机两帧相移条纹图的三维测量方法,可以应用到动态三维测量中,同时,能够满足工业在线检测、人脸检测识别、文物保护、医疗整形等领域的需求,降低了三维检测经济成本,而且本发明方法重建的待测物体的三维形貌具有较高的鲁棒性,灵活性和精确度。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂,以下特举较佳实施例,并配合附图,详细说明如下。
附图说明
图1是本发明实施例提供的三维测量系统的结构示意图;
图2是本发明实施例提供的三维测量系统的测量原理图;
图3是本发明实施例提供的一种基于随机两帧相移条纹图的三维测量方法的流程图;
图4-图5是本发明实施例提供的一种参考相位图;
图6-图14是本发明实施例提供的一种仿真结果图。
具体实施方式
为了进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及具体实施方式,对依据本发明提出的一种基于随机两帧相移条纹图的三维测量方法进行详细说明。
有关本发明的前述及其他技术内容、特点及功效,在以下配合附图的具体实施方式详细说明中即可清楚地呈现。通过具体实施方式的说明,可对本发明为达成预定目的所采取的技术手段及功效进行更加深入且具体地了解,然而所附附图仅是提供参考与说明之用,并非用来对本发明的技术方案加以限制。
实施例一
请结合参见图1和图2,图1是本发明实施例提供的三维测量系统的结构示意图,图2是本发明实施例提供的三维测量系统的测量原理图。如图所示,三维测量系统的测量原理为利用投影仪将条纹投影到待测物体表面,待测物体形貌的变化导致条纹发生形变,由CCD相机捕获变形的条纹图,根据变形的条纹图提取物体的相位信息,并根据相位-高度映射关系式获取物体的三维形貌信息。
请参见图3,图3是本发明实施例提供的一种基于随机两帧相移条纹图的三维测量方法的流程图,如图所示,本实施例的基于随机两帧相移条纹图的三维测量方法,包括:
S1:向待测物体依次投影第一条纹图和第二条纹图,得到第一变形条纹图和第二变形条纹图;
在本实施例中,所述第一条纹图和所述第二条纹图之间的相移量为任意相移量,所述任意相移量可以是未知的。利用投影仪将计算机生成的所述第一条纹图和所述第二条纹图依次投影至待测物体表面,待测物体形貌的变化导致条纹发生形变,得到所述第一变形条纹图和所述第二变形条纹图,由CCD相机捕获所述第一变形条纹图和所述第二变形条纹图。
所述第一变形条纹图的强度为,
Figure BDA0002277616100000081
所述第二变形条纹图的强度为,
Figure BDA0002277616100000082
其中,(x,y)表示图像的像素点坐标,A1(x,y)表示所述第一变形条纹图的背景强度,B1(x,y)表示所述第一变形条纹图的调制度幅值,A2(x,y)表示所述第二变形条纹图的背景强度,B2(x,y)表示所述第二变形条纹图的调制度幅值,
Figure BDA0002277616100000083
表示待测物体的相位信息,δ表示第一变形条纹图与第二变形条纹图的相移量。
S2:分别对所述第一变形条纹图和所述第二变形条纹图进行背景滤除和调制度幅值矫正,得到第一矫正结果和第二矫正结果;
具体地,包括:
S21:利用高通滤波器分别将所述第一变形条纹图和所述第二变形条纹图的背景滤除,得到第一滤除背景图和第二滤除背景图,其强度表达式如下,
Figure BDA0002277616100000084
Figure BDA0002277616100000085
S22:根据归一化公式分别对所述第一滤除背景图和所述第二滤除背景图进行归一化处理,得到所述第一矫正结果和所述第二矫正结果,其中,
所述归一化公式如下:
Figure BDA0002277616100000091
其中,unorm表示归一化后的向量,<>表示内积,||||表示求范数;
所述第一矫正结果为,
Figure BDA0002277616100000092
所述第二矫正结果为,
Figure BDA0002277616100000093
其中,k=1,2,…,K表示所述变形条纹图中的像素点。
通过归一化处理,可以消除由于所述第一变形条纹图和所述第二变形条纹图的调制度幅值不相等,即,B1(x,y)≠B2(x,y),对所述待测物体的相位信息提取带来的影响。
S3:利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到所述待测物体的相位信息;
具体地,包括:
S31:利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到利萨如椭圆表达式;
在本实施例中,将所述第一矫正结果和所述第二矫正结果分别相加和相减,省略空间相关性(x,y),简化后得到下式,
Figure BDA0002277616100000094
Figure BDA0002277616100000095
将公式(8)和公式(9)进行变换,得到,
Figure BDA0002277616100000101
Figure BDA0002277616100000102
根据
Figure BDA0002277616100000103
对公式(10)和公式(11)进行变换,得到利萨如椭圆表达式,
Figure BDA0002277616100000104
其中,x0=0,y0=0,(x0,y0)表示利萨如椭圆的圆心,
Figure BDA0002277616100000105
Figure BDA0002277616100000106
S32:根据所述利萨如椭圆表达式,得到所述第一变形条纹图和所述第二变形条纹图之间的相移量;
在本实施例中,对公式(12)进行变形转换,得到二次方程,
Figure BDA0002277616100000107
其中,
Figure BDA0002277616100000108
采用最小二乘法对所述二次方程进行求解,得到系数a和b,
Figure BDA0002277616100000109
根据所述系数a和b,得到所述第一变形条纹图和所述第二变形条纹图之间的相移量,
Figure BDA0002277616100000111
S33:根据所述相移量,得到所述测物体的相位信息。
在本实施例中,根据所述相移量,按照下式计算得到所述测物体的相位信息,
Figure BDA0002277616100000112
S4:根据所述相位信息,得到所述待测物体的三维形貌信息。
具体地,包括:
根据所述相位信息利用相位-高度映射关系,得到所述待测物体的三维形貌信息,
Figure BDA0002277616100000113
其中,l表示在测量系统中相机到待测物体参考平面的距离,d表示在测量系统中相机与投影仪之间的距离,f表示条纹图的空间频率。
本实施例的基于随机两帧相移条纹图的三维测量方法,将利萨如椭圆拟合技术提取物体的相位信息应用到条纹投影系统,仅采用两帧具有不均匀背景强度和调制度的条纹图重建物体的三维形貌,与单帧傅里叶方法相比,本发明提高了物体的重建精度,解决了傅里叶轮廓术由于频谱丢失导致的重建精度不高的问题,与多帧相移方法相比,本发明所提方法仅用两帧条纹提取物体的相位信息,减少了条纹帧数,计算更简便。另外,本实施例的方法,可以应用到动态三维测量中,同时,能够满足工业在线检测、人脸检测识别、文物保护、医疗整形等领域的需求,降低了三维检测经济成本,而且本发明方法重建的待测物体的三维形貌具有较高的鲁棒性,灵活性和精确度。
实施例二
本实施例是对实施例一的基于随机两帧相移条纹图的三维测量方法与四步相移方法、傅里叶轮廓术测量方法进行对比仿真试验,以说明本发明的方法具有较高的鲁棒性和精确度。
本实施例提供的参考相位为
Figure BDA0002277616100000121
其大小为801×801,请参见图4和图5,图4和图5是本发明实施例提供的一种参考相位图,其中图4是参考相位的相位图,图5是参考相位的俯视图。
在本实施例中,为了验证本发明方法的有效性和鲁棒性,将所述第一条纹图与所述第二条纹图的背景强度和调制度幅值设置为不均匀变化,其中,所述第一条纹图的背景强度和调制度幅值分别为A1=0.5和B1=0.25,所述第二条纹图的背景强度和调制度幅值分别为A2=0.4和B2=0.3。
请参见图6-图14,图6-图14是本发明实施例提供的一种仿真结果图图。图6-图8分别是四步相移方法的仿真结果图,图6-图8分别是四步相移方法的仿真结果图,其中,图6是重建相位的相位图,图7是重建相位的俯视图,图8是残差图,残差图表示参考相位和重建相位的差值。图9-图11分别是傅里叶轮廓术的仿真结果图,其中,图9是重建相位的相位图,图10是重建相位的俯视图,图11是残差图。图12-图14分别是本发明方法的仿真结果图,其中,图12是重建相位的相位图,图13是重建相位的俯视图,图14是残差图。通过比较图7、图10、和图13,可以看出图10的傅里叶轮廓术方法由于频谱丢失而出现相位细节模糊的现象。通过计算可以得道,图8和图14的均方根值为4.3784×10-14和4.3714×10-14,而图11的均方根值为0.2508,均方根值越小,表示重建相位图的精确度越高。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (9)

1.一种基于随机两帧相移条纹图的三维测量方法,其特征在于,包括:
向待测物体依次投影第一条纹图和第二条纹图,得到第一变形条纹图和第二变形条纹图;
分别对所述第一变形条纹图和所述第二变形条纹图进行背景滤除和调制度幅值矫正,得到第一矫正结果和第二矫正结果;
利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到所述待测物体的相位信息;
根据所述相位信息,得到所述待测物体的三维形貌信息。
2.根据权利要求1所述的基于随机两帧相移条纹图的三维测量方法,其特征在于,所述第一条纹图和所述第二条纹图之间的相移量为任意相移量。
3.根据权利要求1所述的基于随机两帧相移条纹图的三维测量方法,其特征在于,所述第一变形条纹图的强度为,
Figure FDA0002277616090000011
所述第二变形条纹图的强度为,
Figure FDA0002277616090000012
其中,(x,y)表示图像的像素点坐标,A1(x,y)表示所述第一变形条纹图的背景强度,B1(x,y)表示所述第一变形条纹图的调制度幅值,A2(x,y)表示所述第二变形条纹图的背景强度,B2(x,y)表示所述第二变形条纹图的调制度幅值,
Figure FDA0002277616090000013
表示待测物体的相位信息,δ表示第一变形条纹图与第二变形条纹图的相移量。
4.根据权利要求3所述的基于随机两帧相移条纹图的三维测量方法,其特征在于,分别对所述第一变形条纹图和所述第二变形条纹图进行背景滤除和调制度幅值矫正,得到第一矫正结果和第二矫正结果,包括:
利用高通滤波器分别将所述第一变形条纹图和所述第二变形条纹图的背景滤除,得到第一滤除背景图和第二滤除背景图,其强度表达式如下,
Figure FDA0002277616090000021
Figure FDA0002277616090000022
根据归一化公式分别对所述第一滤除背景图和所述第二滤除背景图进行归一化处理,得到所述第一矫正结果和所述第二矫正结果,其中,
所述归一化公式如下:
Figure FDA0002277616090000023
其中,unorm表示归一化后的向量,<>表示内积,|| ||表示求范数;
所述第一矫正结果为,
Figure FDA0002277616090000024
所述第二矫正结果为,
Figure FDA0002277616090000025
其中,k=1,2,…,K表示所述变形条纹图中的像素点。
5.根据权利要求4所述的基于随机两帧相移条纹图的三维测量方法,其特征在于,利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到所述待测物体的相位信息,包括:
利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到利萨如椭圆表达式;
根据所述利萨如椭圆表达式,得到所述第一变形条纹图和所述第二变形条纹图之间的相移量;
根据所述相移量,得到所述测物体的相位信息。
6.根据权利要求5所述的基于随机两帧相移条纹图的三维测量方法,其特征在于,利用利萨如椭圆拟合方法对所述第一矫正结果和所述第二矫正结果进行处理,得到利萨如椭圆表达式,包括:
将所述第一矫正结果和所述第二矫正结果分别相加和相减,简化后得到下式,
Figure FDA0002277616090000031
Figure FDA0002277616090000032
根据上式,得到利萨如椭圆表达式,
Figure FDA0002277616090000033
其中,x0=0,y0=0,(x0,y0)表示利萨如椭圆的圆心,
Figure FDA0002277616090000034
Figure FDA0002277616090000035
7.根据权利要求6所述的基于随机两帧相移条纹图的三维测量方法,其特征在于,根据所述利萨如椭圆表达式,得到所述第一变形条纹图和所述第二变形条纹图之间的相移量,包括:
对所述利萨如椭圆表达式进行变形转换,得到二次方程,
Figure FDA0002277616090000036
其中,
Figure FDA0002277616090000037
采用最小二乘法对所述二次方程进行求解,得到系数a和b,
Figure FDA0002277616090000041
根据所述系数a和b,得到所述第一变形条纹图和所述第二变形条纹图之间的相移量,
Figure FDA0002277616090000042
8.根据权利要求7所述的基于随机两帧相移条纹图的三维测量方法,其特征在于,根据所述相移量,得到所述测物体的相位信息,包括:
根据所述相移量,按照下式计算得到所述测物体的相位信息,
Figure FDA0002277616090000043
9.根据权利要求8所述的基于随机两帧相移条纹图的三维测量方法,其特征在于,根据所述相位信息,得到所述待测物体的三维形貌信息,包括:
根据所述相位信息利用相位-高度映射关系,得到所述待测物体的三维形貌信息,
Figure FDA0002277616090000044
其中,l表示在测量系统中相机到待测物体参考平面的距离,d表示在测量系统中相机与投影仪之间的距离,f表示条纹图的空间频率。
CN201911128510.6A 2019-11-18 2019-11-18 一种基于随机两帧相移条纹图的三维测量方法 Active CN111023995B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911128510.6A CN111023995B (zh) 2019-11-18 2019-11-18 一种基于随机两帧相移条纹图的三维测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911128510.6A CN111023995B (zh) 2019-11-18 2019-11-18 一种基于随机两帧相移条纹图的三维测量方法

Publications (2)

Publication Number Publication Date
CN111023995A true CN111023995A (zh) 2020-04-17
CN111023995B CN111023995B (zh) 2021-08-06

Family

ID=70200479

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911128510.6A Active CN111023995B (zh) 2019-11-18 2019-11-18 一种基于随机两帧相移条纹图的三维测量方法

Country Status (1)

Country Link
CN (1) CN111023995B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111707216A (zh) * 2020-06-30 2020-09-25 西安工业大学 一种基于随机两步相移的面形检测方法
CN114396886A (zh) * 2021-12-29 2022-04-26 湖北大学 基于空分复用编码的三维测量方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10035879A1 (de) * 2000-07-24 2002-02-21 Stefan Bickert Abbildung und Erfassung eines dreidimensionalen Objekts
CN1784587A (zh) * 2003-03-06 2006-06-07 齐戈股份有限公司 使用扫描干涉测量形成复杂表面结构的轮廓
CN101128835A (zh) * 2004-05-28 2008-02-20 东京毅力科创株式会社 光学计量中的形状粗糙度测量
CN101865674A (zh) * 2009-04-16 2010-10-20 陈亮嘉 相位信息读取方法及其三维形貌测量系统
CN102376089A (zh) * 2010-12-09 2012-03-14 深圳大学 一种标靶校正方法及系统
CN103048789A (zh) * 2013-01-06 2013-04-17 泉州师范学院 一种产生长距离近似无衍射栅型线结构光的光学系统
CN104111039A (zh) * 2014-08-08 2014-10-22 电子科技大学 用于任意摆放条纹投影三维测量系统的标定方法
CN104937367A (zh) * 2013-01-17 2015-09-23 赛博光学公司 用于电路板的三维成像的多照相机传感器
CN105300317A (zh) * 2015-10-28 2016-02-03 南昌航空大学 基于正弦和三角波条纹投影的三维测量方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10035879A1 (de) * 2000-07-24 2002-02-21 Stefan Bickert Abbildung und Erfassung eines dreidimensionalen Objekts
CN1784587A (zh) * 2003-03-06 2006-06-07 齐戈股份有限公司 使用扫描干涉测量形成复杂表面结构的轮廓
CN101128835A (zh) * 2004-05-28 2008-02-20 东京毅力科创株式会社 光学计量中的形状粗糙度测量
CN101865674A (zh) * 2009-04-16 2010-10-20 陈亮嘉 相位信息读取方法及其三维形貌测量系统
CN102376089A (zh) * 2010-12-09 2012-03-14 深圳大学 一种标靶校正方法及系统
CN103048789A (zh) * 2013-01-06 2013-04-17 泉州师范学院 一种产生长距离近似无衍射栅型线结构光的光学系统
CN104937367A (zh) * 2013-01-17 2015-09-23 赛博光学公司 用于电路板的三维成像的多照相机传感器
CN104111039A (zh) * 2014-08-08 2014-10-22 电子科技大学 用于任意摆放条纹投影三维测量系统的标定方法
CN105300317A (zh) * 2015-10-28 2016-02-03 南昌航空大学 基于正弦和三角波条纹投影的三维测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YUN BAI: "Numerical investigation of two-phase flame structures in a simplified coal jet flame", 《FUEL》 *
朱勇建: "条纹反射法测量三维面形中的积分算法", 《光子学报》 *
杜虎兵: "二帧随机相移条纹图相移提取算法研究", 《激光与光电子学进展》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111707216A (zh) * 2020-06-30 2020-09-25 西安工业大学 一种基于随机两步相移的面形检测方法
CN111707216B (zh) * 2020-06-30 2021-07-02 西安工业大学 一种基于随机两步相移的面形检测方法
CN114396886A (zh) * 2021-12-29 2022-04-26 湖北大学 基于空分复用编码的三维测量方法
CN114396886B (zh) * 2021-12-29 2023-05-09 湖北大学 基于空分复用编码的三维测量方法

Also Published As

Publication number Publication date
CN111023995B (zh) 2021-08-06

Similar Documents

Publication Publication Date Title
Gorthi et al. Fringe projection techniques: whither we are?
CN103940371A (zh) 一种用于跃变物体的高精度三维面形测量的方法
CN108759709B (zh) 一种适用于表面形貌检测的白光干涉三维重建方法
CN106032976B (zh) 基于波长选择的三条纹投影相位展开方法
CN109186496B (zh) 一种基于移动最小二乘法的三维面形测量方法
CN111023995B (zh) 一种基于随机两帧相移条纹图的三维测量方法
CN111174730B (zh) 一种基于相位编码的快速相位解包裹方法
TWI489101B (zh) 結合三維及二維形貌之量測方法及裝置
Li et al. Shearlet transform for phase extraction in fringe projection profilometry with edges discontinuity
CN107388963B (zh) 将小波分析和低通滤波结合的数字莫尔条纹相位提取方法
CN109579738B (zh) 一种二值条纹离焦投影系统低通滤波特性测量方法
Li et al. A 3D shape retrieval method for orthogonal fringe projection based on a combination of variational image decomposition and variational mode decomposition
CN106767531A (zh) 一种二维运动物体的三维测量方法及装置
CN114136236A (zh) 一种结合极线约束的相位轮廓自适应投影测量方法
Hou et al. Relative gradients for image lighting correction
CN111023996B (zh) 一种单帧动态三维测量方法
Hu et al. A new method and implementation of blind restoration algorithm for moving fuzzy license plate image based on frequency-domain characteristics
Yagnik et al. 3D shape extraction of human face in presence of facial hair: A profilometric approach
Osten 19. Digital Image Processing for Optical Metrology
CN115524311A (zh) 基于斜条纹图案的表面漫反射率和三维形貌一体成像方法
CN111597880B (zh) 一种基于单幅条纹投影的多视角三维指纹获取方法
CN107741204A (zh) 一种用于动态三维测量的条纹增强方法
Kayaba et al. Non-contact full field vibration measurement based on phase-shifting
Guo 3-D shape measurement based on Fourier transform and phase shifting method
Vargas et al. Background intensity removal in structured light three-dimensional reconstruction

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