CN109506590A - 一种边界跃变相位误差快速定位方法 - Google Patents
一种边界跃变相位误差快速定位方法 Download PDFInfo
- Publication number
- CN109506590A CN109506590A CN201811617247.2A CN201811617247A CN109506590A CN 109506590 A CN109506590 A CN 109506590A CN 201811617247 A CN201811617247 A CN 201811617247A CN 109506590 A CN109506590 A CN 109506590A
- Authority
- CN
- China
- Prior art keywords
- phase
- frequency
- error
- formula
- striped
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及结构光三维测量技术领域,具体为一种边界跃变相位误差快速定位方法,包括如下步骤:S1、采用不同频率正弦条纹图进行投影测量,其中频率指单幅条纹图中条纹数目,各频率之间满足一定关系,确保通过计算最终能获得单频相位图;S2、通过相位展开技术获取各频率展开相位图,其中包含单频相位图;S3、利用各个频率展开相位,计算时间相位梯度;S4、评估时间相位梯度与单周期相位的残差,对残差进行阈值化处理,定位跃变区域相位误差。本发明提供了一种简单便捷适用性更广的边界跃变相位误差快速定位方法。
Description
技术领域
本发明涉及结构光三维测量技术领域,具体为一种边界跃变相位误差快速定位方法。
背景技术
结构光3D面型检测技术,因其具有快速,非接触,高精度等优点,被认为最具前景的光学三维面形测量技术之一,该技术利用投影设备将条纹投影到待测物体表面,被测物体形貌的变化导致条纹发生扭曲,形变,对这种条纹形变进行分析来获取高度信息。在对不连续区域进行测量时,常常采用多频条纹投影,利用时间相位展开算法有效避免错误相位的传播,然而错误相位依然存在与相位图中。不连续区域的高度跃变往往容易造成较大相位误差,高度跃变区域相位误差主要源于跃变边缘被相机单个像元积分曝光,导致条纹图在边缘区域出现失真,无法准确反映相位影响。为了便于后期对3D数据进一步处理,如相位误差评估、相位修复、3D数据处理等工作,所有对跃变相位误差区域的定位显得非常必要。
传统方法利用多频展开相位对每个像素点在时间序列上进行线性拟合,通过比较拟合曲线的标准差是否超过阈值,判断该点是否为跃变边缘相位误差点,由于线性拟合方法往往需要多组数据,也需要多组不同频率条纹,即便是所需改进后的逆指数时间相位展开方法一般投影条纹也不会少于6种频率。但在实际测量中,为提高测量速度,往往只采用双频或三频投影,因此单个像素上进行线性拟合的数据点仅仅2个或3个,此时就失去了拟合的意义,因此传统跃变相位误差定位方法在实际使用中受到限制;为此,我们提出一种边界跃变相位误差快速定位方法。
发明内容
本发明的目的在于提供一种边界跃变相位误差快速定位方法,以解决上述背景技术中提出的现有结构光测量系统中对物体形貌跃变边界相位误差定位方法的不足的问题。
为实现上述目的,本发明提供如下技术方案:
一种边界跃变相位误差快速定位方法,包括如下步骤:
S1、采用不同频率正弦条纹图进行投影测量,其中频率指单幅条纹图中条纹数目,各频率之间满足一定关系,确保通过计算最终能获得单频相位图;
S11、利用计算机编码生成空间频率为f1、f2、f3...fn-1、fn的正弦条纹图,n为大于1的整数,各频率选取需确保通过计算最终能获得单频相位图;
S2、通过相位展开技术获取各频率展开相位图,其中包含单频相位图;
S21、将编码生成的频率为f1、f2、f3...fn-1、fn的正弦条纹图投影至待测物体表面;
S22、利用相机捕获不同频率变形条纹图;
S23、利用时间相位展开技术获得各频率展开相位为:并通过计算得到单频相位为
S3、利用各个频率展开相位,计算时间相位梯度;
S31、利用展开相位计算时间相位梯度具体计算方式如下所示:
S4、评估时间相位梯度与单周期相位的残差,对残差进行阈值化处理,定位跃变区域相位误差。
优选的,所述步骤S4包括如下步骤:
S41、取为频率ft的条纹在(x,y)位置的展开相位,为单频相位在(x,y)位置的相位值,则在理想情况下有如下关系:
而实际测量过程中,在高度跃变区域由于条纹遮挡以及相机像元的积分曝光过程,使得跃变边界条纹信息失真,导致边界区域相位值存在较大误差,且条纹信息的失真程度与投影条纹频率相关,不同频率条纹失真情况不同。单频条纹由于条纹灰度变化平缓,相位值动态范围小,因此跃变误差影响极小。由于跃变误差与频率相关,取频率为ft的相位图在(x,y)位置的跃变误差为Dt(x,y),增加具体相位解调过程相位偏移量St,该偏移量与为整体偏移,与像素位置(x,y)无关,进一步式(2)可变为如下:
将式(3)代入式(1)有:
如式(4)所示,其中
Et(x,y)项与跃变区域所造成的相位误差相关,Mt与整体偏移量相关与位置(x,y)无关;
S42、利用式(1)计算获得的时间相位梯度与单频相位做差可得到,残差项:
其中偏移量Mt为常数,考虑跃变区域相对于整个测量区域而言,为小范围区域,Mt可视Rt(x,y)的均值,因此,根据式(1)、(5)可以得到跃变区域误差相关项Et(x,y)的表达式:
S43、式(6)给出了跃变区域误差相关项Et(x,y)的表达式,Et(x,y)与频率ft+1和ft相关,当投影频率为f1、f2、f3...fn-1、fn,n种频率时,为更加合理的利用所有频率,对所有Et(x,y)求均值得到E(x,y)如下所示:
式(7)中,abs{*}为求绝对值函数,对所得跃变区域误差相关项E(x,y)设置阈值T,当E(x,y)在(μ-T,μ+T)之外时,则判定为跃变相位误差区域。
与现有技术相比,本发明的有益效果是:本发明提供了一种简单便捷适用性更广的边界跃变相位误差快速定位方法;在对跃变区域相位误差进行定位时,无需进行复杂的线性拟合,仅仅需要计算不同频率时间方向的相位梯度并与单频率相位进行比较即可。该发明计算过程简单,没有复杂的拟合算法仅仅涉及矩阵间的加、减法;对数据要求量少,当投影频率数为2时,通过优化误差判定阈值,同样能够实现跃变区域相位误差定位效果。
附图说明
图1为本发明原理流程框图;
图2为本发明实施例中H(x,y)仿真面型的三维图和俯视图;
图3为本发明实施例中I1、I2、I3三种不同频率条纹经过相同高度调整后的条纹图;
图4为本发明实施例中对应展开相位;
图5为本发明实施例中R1(x,y)、R2(x,y)对应相位梯度残差图;
图6为本发明实施例中E(x,y)跃变相位误差定位区域。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1-6,本发明提供一种技术方案:
一种边界跃变相位误差快速定位方法,包括如下步骤:
S1、采用不同频率正弦条纹图进行投影测量,其中频率指单幅条纹图中条纹数目,各频率之间满足一定关系,确保通过计算最终能获得单频相位图;
S11、利用计算机编码生成空间频率为f1、f2、f3...fn-1、fn的正弦条纹图,n为大于1的整数,各频率选取需确保通过计算最终能获得单频相位图;
S2、通过相位展开技术获取各频率展开相位图,其中包含单频相位图;
S21、将编码生成的频率为f1、f2、f3...fn-1、fn的正弦条纹图投影至待测物体表面;
S22、利用相机捕获不同频率变形条纹图;
S23、利用时间相位展开技术获得各频率展开相位为:并通过计算得到单频相位为
S3、利用各个频率展开相位,计算时间相位梯度;
S31、利用展开相位计算时间相位梯度具体计算方式如下所示:
S4、评估时间相位梯度与单周期相位的残差,对残差进行阈值化处理,定位跃变区域相位误差;
S41、取为频率ft的条纹在(x,y)位置的展开相位,为单频相位在(x,y)位置的相位值,则在理想情况下有如下关系:
而实际测量过程中,在高度跃变区域由于条纹遮挡以及相机像元的积分曝光过程,使得跃变边界条纹信息失真,导致边界区域相位值存在较大误差,且条纹信息的失真程度与投影条纹频率相关,不同频率条纹失真情况不同。单频条纹由于条纹灰度变化平缓,相位值动态范围小,因此跃变误差影响极小。由于跃变误差与频率相关,取频率为ft的相位图在(x,y)位置的跃变误差为Dt(x,y),增加具体相位解调过程相位偏移量St,该偏移量与为整体偏移,与像素位置(x,y)无关,进一步式(2)可变为如下:
将式(3)代入式(1)有:
如式(4)所示,其中
Et(x,y)项与跃变区域所造成的相位误差相关,Mt与整体偏移量相关与位置(x,y)无关;
S42、利用式(1)计算获得的时间相位梯度与单频相位做差可得到,残差项:
其中偏移量Mt为常数,考虑跃变区域相对于整个测量区域而言,为小范围区域,Mt可视Rt(x,y)的均值,因此,根据式(1)、(5)可以得到跃变区域误差相关项Et(x,y)的表达式:
S43、式(6)给出了跃变区域误差相关项Et(x,y)的表达式,Et(x,y)与频率ft+1和ft相关,当投影频率为f1、f2、f3...fn-1、fn,n种频率时,为更加合理的利用所有频率,对所有Et(x,y)求均值得到E(x,y)如下所示:
式(7)中,abs{*}为求绝对值函数,对所得跃变区域误差相关项E(x,y)设置阈值T,当E(x,y)在(μ-T,μ+T)之外时,则判定为跃变相位误差区域。
实施例
如图2-6,如上的边界跃变相位误差快速定位方法的一种实施例:
一、如图2所示,计算机仿真生成最高高度为10cm的斜面体H,左图为仿真面型3D图,右图为仿真面型俯视图;
二、计算机仿真,采用三步相移法计算机生成空间频率为f1、f2、f3的条纹图,其中f1=29、f2=34、f3=40。将仿真体H对条纹图进行调制后获得变形条纹图,如图3所示,I1、I2、I3分别为f1、f2、f3对于的变形条纹图。为仿真实际相机拍摄过程的曝光积分效果,对I1、I2、I3进行窗口为3*3大小的高斯滤波。另外考虑随机噪声,在动态范围为0~255的变形条纹图上加标准差为1的高斯噪声;
三、利用外差法和多步相移技术解调变形条纹图,获得各个频率展开相位以及单频相位如图4所示;
四、利用以及单频相位结合公式(6),计算E1(x,y)、E2(x,y)如图5所示;
五、利用E1(x,y)、E2(x,y),结合公式(7),求得E(x,y),计算E(x,y)的均方根δ,取阈值T=3*δ,得到阈值处理后的E(x,y),即为跃变边界相位误差区域,如图6所示,白色区域为跃变边界相位误差区域。
本发明提供了一种简单便捷适用性更广的边界跃变相位误差快速定位方法;在对跃变区域相位误差进行定位时,无需进行复杂的线性拟合,仅仅需要计算不同频率时间方向的相位梯度并与单频率相位进行比较即可。该发明计算过程简单,没有复杂的最小二乘拟合算法,仅仅涉及矩阵间的加、减法;也无需额外采集数据,对数据要求量少,适用大部分时间相位展开条纹投影测量技术。当投影频率数为2时,通过优化误差判定阈值,同样能够实现跃变区域相位误差定位效果。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (2)
1.一种边界跃变相位误差快速定位方法,其特征在于,包括如下步骤:
S1、采用不同频率正弦条纹图进行投影测量,其中频率指单幅条纹图中条纹数目,各频率之间满足一定关系,确保通过计算最终能获得单频相位图;
S11、利用计算机编码生成空间频率为f1、f2、f3...fn-1、fn的正弦条纹图,n为大于1的整数,各频率选取需确保通过计算最终能获得单频相位图;
S2、通过相位展开技术获取各频率展开相位图,其中包含单频相位图;
S21、将编码生成的频率为f1、f2、f3...fn-1、fn的正弦条纹图投影至待测物体表面;
S22、利用相机捕获不同频率变形条纹图;
S23、利用时间相位展开技术获得各频率展开相位为:并通过计算得到单频相位为
S3、利用各个频率展开相位,计算时间相位梯度;
S31、利用展开相位计算时间相位梯度具体计算方式如下所示:
S4、评估时间相位梯度与单周期相位的残差,对残差进行阈值化处理,定位跃变区域相位误差。
2.根据权利要求1所述的一种边界跃变相位误差快速定位方法,其特征在于,所述步骤S4包括如下步骤:
S41、取为频率ft的条纹在(x,y)位置的展开相位,为单频相位在(x,y)位置的相位值,则在理想情况下有如下关系:
而实际测量过程中,在高度跃变区域由于条纹遮挡以及相机像元的积分曝光过程,使得跃变边界条纹信息失真,导致边界区域相位值存在较大误差,且条纹信息的失真程度与投影条纹频率相关,不同频率条纹失真情况不同。单频条纹由于条纹灰度变化平缓,相位值动态范围小,因此跃变误差影响极小。由于跃变误差与频率相关,取频率为ft的相位图在(x,y)位置的跃变误差为Dt(x,y),增加具体相位解调过程相位偏移量St,该偏移量与为整体偏移,与像素位置(x,y)无关,进一步式(2)可变为如下:
将式(3)代入式(1)有:
如式(4)所示,其中
Et(x,y)项与跃变区域所造成的相位误差相关,Mt与整体偏移量相关与位置(x,y)无关;
S42、利用式(1)计算获得的时间相位梯度与单频相位做差可得到,残差项:
其中偏移量Mt为常数,考虑跃变区域相对于整个测量区域而言,为小范围区域,Mt可视Rt(x,y)的均值,因此,根据式(1)、(5)可以得到跃变区域误差相关项Et(x,y)的表达式:
S43、式(6)给出了跃变区域误差相关项Et(x,y)的表达式,Et(x,y)与频率ft+1和ft相关,当投影频率为f1、f2、f3...fn-1、fn,n种频率时,为更加合理的利用所有频率,对所有Et(x,y)求均值得到E(x,y)如下所示:
式(7)中,abs{*}为求绝对值函数,对所得跃变区域误差相关项E(x,y)设置阈值T,当E(x,y)在(μ-T,μ+T)之外时,则判定为跃变相位误差区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811617247.2A CN109506590B (zh) | 2018-12-28 | 2018-12-28 | 一种边界跃变相位误差快速定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811617247.2A CN109506590B (zh) | 2018-12-28 | 2018-12-28 | 一种边界跃变相位误差快速定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109506590A true CN109506590A (zh) | 2019-03-22 |
CN109506590B CN109506590B (zh) | 2020-10-27 |
Family
ID=65755593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811617247.2A Active CN109506590B (zh) | 2018-12-28 | 2018-12-28 | 一种边界跃变相位误差快速定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109506590B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110793463A (zh) * | 2019-09-25 | 2020-02-14 | 西安交通大学 | 一种基于相位分布的解包裹相位误差检测与校正方法 |
CN113393481A (zh) * | 2021-06-10 | 2021-09-14 | 湖南大学 | 基于边缘检测的快速相位展开方法、装置、设备及介质 |
CN115002344A (zh) * | 2022-05-11 | 2022-09-02 | 昆山丘钛微电子科技股份有限公司 | 误差曲面拟合方法、装置、设备及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1537534A2 (en) * | 2002-09-12 | 2005-06-08 | Nline Corporation | System and method for acquiring and processing complex images |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1695166A (zh) * | 2002-09-12 | 2005-11-09 | 恩莱因公司 | 用于获取和处理复合图像的系统及方法 |
JP4756024B2 (ja) * | 2007-11-12 | 2011-08-24 | アンリツ株式会社 | 三次元形状測定装置 |
CN100541148C (zh) * | 2007-11-21 | 2009-09-16 | 北京理工大学 | 一种光纤傅立叶变换白光干涉测量法 |
CN101650164B (zh) * | 2009-09-07 | 2011-02-09 | 东南大学 | 三维扫描系统中基于错位条纹的相位展开方法 |
CN102620685B (zh) * | 2012-03-23 | 2014-11-26 | 东南大学 | 一种基于史托克维尔变换的改进窗口傅里叶三维测量法 |
CN102607466B (zh) * | 2012-03-29 | 2014-10-01 | 天津大学 | 高反射自由曲面光栅投影快速非接触测量方法及装置 |
CN103487441B (zh) * | 2013-09-24 | 2015-09-30 | 电子科技大学 | 一种用于硅晶片缺陷检测和面形测量的方法 |
CN103727895B (zh) * | 2014-01-16 | 2017-02-22 | 电子科技大学 | 基于单帧彩色复合光栅条纹反射的镜面三维面形测量方法 |
CN103940371A (zh) * | 2014-05-12 | 2014-07-23 | 电子科技大学 | 一种用于跃变物体的高精度三维面形测量的方法 |
CN104279980A (zh) * | 2014-10-20 | 2015-01-14 | 电子科技大学 | 一种基于智能拍照手机的镜面三维面形测量系统 |
CN104482877B (zh) * | 2014-12-03 | 2017-02-01 | 深圳大学 | 动态物体三维成像中的运动补偿方法与系统 |
TW201800719A (zh) * | 2016-06-23 | 2018-01-01 | 林嘉宇 | 物體辨識系統及辨識物體的方法 |
US10731972B2 (en) * | 2016-11-18 | 2020-08-04 | Arizona Board Of Regents On Behalf Of The University Of Arizona | Instantaneous phase mapping deflectometry |
-
2018
- 2018-12-28 CN CN201811617247.2A patent/CN109506590B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1537534A2 (en) * | 2002-09-12 | 2005-06-08 | Nline Corporation | System and method for acquiring and processing complex images |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110793463A (zh) * | 2019-09-25 | 2020-02-14 | 西安交通大学 | 一种基于相位分布的解包裹相位误差检测与校正方法 |
CN113393481A (zh) * | 2021-06-10 | 2021-09-14 | 湖南大学 | 基于边缘检测的快速相位展开方法、装置、设备及介质 |
CN115002344A (zh) * | 2022-05-11 | 2022-09-02 | 昆山丘钛微电子科技股份有限公司 | 误差曲面拟合方法、装置、设备及介质 |
CN115002344B (zh) * | 2022-05-11 | 2024-04-19 | 昆山丘钛微电子科技股份有限公司 | 误差曲面拟合方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109506590B (zh) | 2020-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109658515B (zh) | 点云网格化方法、装置、设备及计算机存储介质 | |
KR101475382B1 (ko) | 광학적 3차원 측량의 자기 적응 윈도우 푸리에 위상추출방법 | |
Sagawa et al. | Dense one-shot 3D reconstruction by detecting continuous regions with parallel line projection | |
Wieneke | Improvements for volume self-calibration | |
CN106032976A (zh) | 基于波长选择的三条纹投影相位展开方法 | |
CN109506590A (zh) | 一种边界跃变相位误差快速定位方法 | |
Garcia et al. | Consistent stereo-assisted absolute phase unwrapping methods for structured light systems | |
CN110163817B (zh) | 一种基于全卷积神经网络的相位主值提取方法 | |
EP2881702A1 (en) | Image processing system, and image processing method | |
CN109141291A (zh) | 一种快速相位解包裹算法 | |
CN103822598B (zh) | 地基sar在时间去相关严重区域的形变监测方法 | |
CN105547190B (zh) | 基于双角度单频率条纹投影的三维形貌测量方法及装置 | |
CN111563952B (zh) | 基于相位信息与空间纹理特征实现立体匹配方法及系统 | |
CN109631797B (zh) | 一种基于相移技术的三维重构无效区域快速定位方法 | |
CN110268223A (zh) | 三维形状计测装置、三维形状计测方法以及程序 | |
CN111947599B (zh) | 基于学习的条纹相位恢复和散斑相关的三维测量方法 | |
CN105066904B (zh) | 基于相位梯度阈值的流水线产品三维面型检测方法 | |
JP2024507089A (ja) | 画像のコレスポンデンス分析装置およびその分析方法 | |
CN104019767B (zh) | 一种基于灰度直方图的自适应高动态范围三维测量方法 | |
CN110223384A (zh) | 一种白光干涉三维形貌重建方法、装置、系统及存储介质 | |
CN111998799A (zh) | 一种多频投影三维面形测量的运动区域检测方法及系统 | |
Qi et al. | Invalid-point removal based on epipolar constraint in the structured-light method | |
CN105588518B (zh) | 基于双角度多频率条纹投影的三维形貌获取方法及装置 | |
Wang et al. | A phase unwrapping method based on sparse depth for depth perception | |
CN205718873U (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 |