CN111207694B - 一种双步相移法结合相位编码的三维测量方法 - Google Patents
一种双步相移法结合相位编码的三维测量方法 Download PDFInfo
- Publication number
- CN111207694B CN111207694B CN202010032926.4A CN202010032926A CN111207694B CN 111207694 B CN111207694 B CN 111207694B CN 202010032926 A CN202010032926 A CN 202010032926A CN 111207694 B CN111207694 B CN 111207694B
- Authority
- CN
- China
- Prior art keywords
- phase
- phi
- fringe
- phase shift
- wrapped
- 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
Links
Images
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
- G01B11/25—Measuring 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/254—Projection of a pattern, viewing through a pattern, e.g. moiré
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种双步相移法结合相位编码的三维测量方法,由双步相移原理、相位融合方法、相位编码原理三大关键部分组成。本发明的优点是:与传统的相移算法相比,双步相移算法能够更大地减少测量误差,且能较好用于实际测量;直接融合两幅包裹相位图的相位信息,只需对一幅包裹相位图进行相位展开,一定程度上减少了出错率;与基于强度编码的方法相比,相位编码的方法对表面对比度、环境光和相机噪声不敏感,具有较强的鲁棒性。
Description
技术领域
本发明涉及一种关于相位编码条纹的三维测量方法,属于光电检测领域,尤其涉及一种双步相移法结合相位编码的三维测量方法。
背景技术
条纹投影轮廓术因其高速、高精度、非接触、全场测量、快速信息获取等优点,在三维测量中有重要意义,已广泛应用于工业制造、文物保护、医疗等诸多领域。
通过对国内外研究现状及发展动向分析研究,传统的三维测量技术己经发展较为成熟,但多为静态测量,无法实现运动物体的三维测量和形貌重构。近年来高速、实时和高精度的三维测量在物体形变分析、工业自动检测、医学诊断、人脸识别等方面都有着广泛的应用,并且随着采集设备、投影设备、高速处理器的性能提升,高速、实时和高分辨率的三维测量方法成为研究的热点。因此,提高三维测量的精度成为一个亟待解决的问题。。
本发明提出一种双步相移法结合相位编码的三维测量方法,此方法需要投影6幅正弦条纹图获得包裹相位,虽然投影的图数增多了,但其精度得到了提升,且误差率大大减小,适用于那些对于精度和准确度要求高的场合。其次,本发明采用的直接融合两幅包裹相位图的相位信息融合方法比传统的相位融合方法更能准确的展开包裹相位,能有效的减少包裹相位图中的噪声或者不连续现象,使相位展开变得容易。
本发明所述的一种双步相移法结合相位编码的三维测量方法正是在这一背景下提出的。
发明内容
本发明的目的在于提出一种双步相移法结合相位编码的三维测量方法,此方法较传统的正弦条纹加相位编码的三维测量方法,在测量精度上有了明显提高,适用于那些对于精度和准确度要求高的场合。
一种双步相移法结合相位编码的三维测量方法,包括以下步骤:
步骤一、利用计算机生成三幅初始相位为0的正弦条纹图、三幅初始相位为π/3的正弦条纹、三幅阶梯相位编码条纹图;
步骤二、将步骤一中所有条纹图通过投影仪投射依次到参考平面和物体上,再通过相机采集投影到参考平面和物体上的六幅正弦条纹图和三幅阶梯相位编码条纹图;
步骤三、利用六幅正弦条纹图得到物体的包裹相位,利用阶梯相位编码条纹图得到物体的条纹级次,进而求得参考平面和物体的连续相位;
步骤四、通过获得的连续相位并通过相位-高度公式求得物体表面每一点的高度信息。
所述步骤一的具体实施方式:
A、编写正弦条纹:
实际测量时先采用三步相移算法进行测量,则正弦条纹图可以表示为:
Ii(x,y)=Ia(x,y)+Ib(x,y)cos[φ(x,y)+2πn/3] (1)
其中,Ia(x,y)背景光强,Ib(x,y)是调制幅度,φ(x,y)是理想相位信息;
三步相移算法的相位误差函数可近似表示为:
Δφ(x,y)≈c1 sin[3φ(x,y)] (2)
其中c1为一个常数;
Δφ'(x,y)≈-c2 sin[3φ(x,y)+3×π/3]
=c2 sin[3φ(x,y)]=-Δφ(x,y) (6)
则此时实际包裹相位可以由两种形式表示:
由投影的两组正弦条纹图可以求得两组实际包裹相位信息φ′1和φ′2,可以削弱或去除三步相移算法对应的相位误差,以得到理想的正确相位φ(x,y);
B、进行相位融合:
直接融合两幅包裹相位图的相位信息融合,可以削弱或去除三步相移算法对应的相位误差,以得到理想的正确相位φ(x,y),根据两组包裹相位信息φ′1和φ′2之间的关系,直接融合两组包裹相位信息,方法如下式:
C、编写相位编码条纹:
相位编码方法的码字需预先设计好,首先设置初始相位,如下式:
其中,floor[·]表示不大于x的最大整数;mod(x,y)是x除以y所得的余数;S为投影仪的横向分辨率;P为条纹间距,即每个周期的像素数;N为设定码字的个数;
将相位归一化到[-π,π],相位归一化公式为:
将归一化相位嵌入相位编码条纹中,把每个像素的相位值代入正弦函数,以获得每个点的条纹强度值,条纹强度分布如下式所示:
An(x,y)=A(x,y)+B(x,y)cos[φs(x,y)+2π(k-1)/3] (12)
其中,A(x,y)为平均亮度,B(x,y)是调制亮度,φs(x,y)为待求相位。生成A1(x,y)、A2(x,y)、A3(x,y)三幅阶梯相位编码条纹图。
所述步骤三的具体实施方式为:
通过反正切三步相移算法计算阶梯相位:
根据A结果,由三步相移法计算得到包裹相位φ′1和φ′2:
根据B结果,融合得到φ(x,y),
根据公式(12)和(14),获得阶梯相位φs(x,y),确定条纹级次k,
k(x,y)=Round[N(φ(x,y)+π)/2π] (15)
其中,Round[]为确定最接近的整数;
根据下式将包裹相位φ(x,y)转换为绝对相位Φ(x,y):
Φ(x,y)=φ(x,y)+2π×k(x,y) (16)
进行相位解包裹,分别得到参考面的连续相位值α和带物体的连续相位值β。
所述步骤四的具体实施方式为:
通过公式:
得到物体的连续相位差,从而利用相位-高度公式:
其中d为投影仪与相机之间的距离,l0为投影仪与相机到参考平面的距离为,f0为参考平面上的正弦条纹频率,最后得到物体表面每一点的高度信息。
本发明的优点是:
(1)与传统的相移算法相比,双步相移算法能够更大地减少测量误差,且能较好用于实际测量;
(2)直接融合两幅包裹相位图的相位信息,只需对一幅包裹相位图进行相位展开,一定程度上减少了出错率;
(3)与基于强度编码的方法相比,相位编码的方法对表面对比度、环境光和相机噪声不敏感,具有较强的鲁棒性。
附图说明
图1为本发明的三维测量系统示意图。
图2为本发明的两组正弦条纹图。
图3为本发明的经过相位融合后的包裹相位图。
图4为本发明的三幅相位条纹图。
图5为本发明的包裹相位和条纹级次。
具体实施方式
以下结合附图说明对本发明的实施例作进一步详细描述,但本实施例并不用于限制本发明,凡是采用本发明的相似结构及其相似变化,均应列入本发明的保护范围。
本发明的一种双步相移法结合相位编码的三维测量方法的三维测量系统如图1所示,包括DLP投影仪1、CCD相机2、计算机3、测量支架4、参考平面5和待测物体6;DLP投影仪1和CCD相机2放在测量支架4上;DLP投影仪1、CCD相机2分别通过数据线连接计算机3;待测物体6放在参考平面5上;计算机3内包含图像采集卡、投影软件、测量软件。DLP投影仪1光轴和CCD相机2光轴相交于O点。DLP投影仪1和CCD相机2为同一高度,它们之间的距离为d,它们到参考平面的距离为l0。
一种双步相移法结合相位编码的三维测量方法,包括以下步骤:
(一)编写正弦条纹
实际测量时先采用三步相移算法进行测量,则正弦条纹图可以表示为:
Ii(x,y)=Ia(x,y)+Ib(x,y)cos[φ(x,y)+2πn/3] (1)
其中,Ia(x,y)背景光强,Ib(x,y)是调制幅度,φ(x,y)是理想相位信息。
三步相移算法的相位误差函数可近似表示为:
Δφ(x,y)≈c1 sin[3φ(x,y)] (2)
其中c1为一个常数。
Δφ'(x,y)≈-c2 sin[3φ(x,y)+3×π/3]
=c2 sin[3φ(x,y)]=-Δφ(x,y) (6)
则此时实际包裹相位可以由两种形式表示:
由投影的两组正弦条纹图可以求得两组实际包裹相位信息φ′1和φ′2,可以削弱或去除三步相移算法对应的相位误差,以得到理想的正确相位φ(x,y)。需投影的两组三幅正弦条纹图如图2所示。
(二)进行相位融合
直接融合两幅包裹相位图的相位信息融合,可以削弱或去除三步相移算法对应的相位误差,以得到理想的正确相位φ(x,y)。根据两组包裹相位信息φ′1和φ′2之间的关系,直接融合两组包裹相位信息,方法如下式:
(三)编写相位编码条纹
相位编码方法的码字需预先设计好,首先设置初始相位,如下式:
其中,floor[·]表示不大于x的最大整数;mod(x,y)是x除以y所得的余数;S为DLP投影仪1的横向分辨率;P为条纹间距,即每个周期的像素数;N为设定码字的个数。
将相位归一化到[-π,π],相位归一化公式为:
将归一化相位嵌入相位编码条纹中,把每个像素的相位值代入正弦函数,以获得每个点的条纹强度值,条纹强度分布如下式所示:
An(x,y)=A(x,y)+B(x,y)cos[φs(x,y)+2π(k-1)/3] (12)
其中,A(x,y)为平均亮度,B(x,y)是调制亮度,φs(x,y)为待求相位。生成A1(x,y)、A2(x,y)、A3(x,y)三幅阶梯相位编码条纹图,如图4所示。
通过反正切三步相移算法计算阶梯相位:
(四)采集图片
用投影仪依次投射上述生成的正弦条纹和相位相位编码条纹于参考平面5和待测物体6上,用CCD相机2采集这三组条纹图,将结果送入计算机3中。
(五)三维重建
根据(一)结果,由三步相移法计算得到包裹相位φ′1和φ′2:
根据(二)结果,融合得到φ(x,y)。
根据公式(12)和(14),获得阶梯相位φs(x,y),确定条纹级次k,如图3所示:
k(x,y)=Round[N(φ(x,y)+π)/2π] (15)
其中,Round[]为确定最接近的整数。
根据下式将包裹相位φ(x,y)转换为绝对相位Φ(x,y):
Φ(x,y)=φ(x,y)+2π×k(x,y) (16)
进行相位解包裹,分别得到参考平面5的连续相位值α和待测物体6的连续相位值β。
通过公式:
得到待测物体6的连续相位差,从而利用相位-高度公式:
其中DLP投影仪1与CCD相机2之间的距离为d,DLP投影仪1与CCD相机2到参考平面的距离为l0,f0为参考平面5上的正弦条纹频率,最后得到待测物体6表面每一点的高度信息。
不局限于此,任何不经过创造性劳动想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书所限定的保护范围为准。
Claims (1)
1.一种双步相移法结合相位编码的三维测量方法,其特征在于,包括以下步骤:
步骤一、利用计算机生成三幅初始相位为0的正弦条纹图、三幅初始相位为π/3的正弦条纹图、三幅阶梯相位编码条纹图;
步骤二、将步骤一中所有条纹图通过投影仪依次投射到参考平面和物体上,再通过相机采集投影到参考平面和物体上的六幅正弦条纹图和三幅阶梯相位编码条纹图;
步骤三、利用六幅正弦条纹图得到物体的包裹相位,利用阶梯相位编码条纹图得到物体的条纹级次,进而求得参考平面和物体的连续相位;
步骤四、通过获得的连续相位并通过相位-高度公式求得物体表面每一点的高度信息;
步骤一的具体实施方式如下:
A、编写正弦条纹:
实际测量时先采用三步相移算法进行测量,则正弦条纹图可以表示为:
Ii(x,y)=Ia(x,y)+Ib(x,y)cos[φ(x,y)+2πn/3] (1)
其中,Ia(x,y)背景光强,Ib(x,y)是调制幅度,φ(x,y)是理想相位信息;
三步相移算法的相位误差函数可近似表示为:
Δφ(x,y)≈c1sin[3φ(x,y)] (2)
其中c1为一个常数;
其中c2为一个常数;
φ1'(x,y)=φ(x,y)+Δφ(x,y) (5)
Δφ'(x,y)≈-c2sin[3φ(x,y)+3×π/3]
=c2sin[3φ(x,y)]=-Δφ(x,y) (6)
则此时实际包裹相位可以由两种形式表示:
由投影的两组正弦条纹图可以求得两组实际包裹相位φ′1(x,y)和φ'2(x,y),可以削弱或去除三步相移算法对应的相位误差,以得到理想的正确相位φ(x,y);
B、进行相位融合:
直接融合两幅包裹相位图的相位,可以削弱或去除三步相移算法对应的相位误差,以得到理想的正确相位φ(x,y),根据两组包裹相位φ′1(x,y)和φ'2(x,y)之间的关系,直接融合两组包裹相位,方法如下式:
C、编写相位编码条纹:
相位编码方法的码字需预先设计好,首先设置初始相位,如下式:
其中,floor[]表示不大于x的最大整数;mod(x,y)是x除以y所得的余数;S为投影仪的横向分辨率;P为条纹间距,即每个周期的像素数;N为设定码字的个数;
将相位归一化到[-π,π],相位归一化公式为:
将归一化相位嵌入相位编码条纹中,把每个像素的相位值代入正弦函数,以获得每个点的条纹强度值,条纹强度分布如下式所示:
An(x,y)=A(x,y)+B(x,y)cos[φs(x,y)+2π(n-1)/3] (12)
其中,A(x,y)为平均亮度,B(x,y)是调制亮度,φs(x,y)为待求相位,生成A1(x,y)、A2(x,y)、A3(x,y)三幅阶梯相位编码条纹图。
步骤三的具体实施方式如下:
通过反正切三步相移算法计算阶梯相位:
根据A结果,由三步相移法计算得到包裹相位φ′1(x,y)和φ'2(x,y):
根据B结果,融合得到φ(x,y),
根据公式(12)和(14),获得阶梯相位φs(x,y),确定条纹级次k,
k(x,y)=Round[N(φ(x,y)+π)/2π] (15)
其中,Round[]为确定最接近的整数;
根据下式将包裹相位φ(x,y)转换为绝对相位Φ(x,y):
Φ(x,y)=φ(x,y)+2π×k(x,y) (16)
进行相位解包裹,分别得到参考平面的连续相位值α和物体的连续相位值β;
步骤四的具体实施方式如下:
通过公式:
得到物体的连续相位差,从而利用相位-高度公式:
其中d为投影仪与相机之间的距离,l0为投影仪与相机到参考平面的距离,f0为参考平面上的正弦条纹频率,最后得到物体表面每一点的高度信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010032926.4A CN111207694B (zh) | 2020-01-13 | 2020-01-13 | 一种双步相移法结合相位编码的三维测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010032926.4A CN111207694B (zh) | 2020-01-13 | 2020-01-13 | 一种双步相移法结合相位编码的三维测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111207694A CN111207694A (zh) | 2020-05-29 |
CN111207694B true CN111207694B (zh) | 2021-06-04 |
Family
ID=70784303
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010032926.4A Active CN111207694B (zh) | 2020-01-13 | 2020-01-13 | 一种双步相移法结合相位编码的三维测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111207694B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112082510A (zh) * | 2020-09-03 | 2020-12-15 | 南昌航空大学 | 一种降低凹面镜测量中中间轴部分面形恢复错误的方法 |
CN113029039B (zh) * | 2021-01-27 | 2022-04-12 | 四川大学 | 一种基于彩色编码的三维测量绝对相位矫正方法 |
CN113390365B (zh) * | 2021-07-09 | 2022-12-06 | 西安工业大学 | 一种基于随机三步移相干涉的面形重构方法 |
CN113532330B (zh) * | 2021-08-28 | 2022-10-04 | 哈尔滨理工大学 | 一种相位格雷码三维测量方法 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4830871B2 (ja) * | 2007-01-26 | 2011-12-07 | パナソニック電工株式会社 | 3次元形状計測装置及び3次元形状計測方法 |
JP5391053B2 (ja) * | 2009-12-24 | 2014-01-15 | ローランドディー.ジー.株式会社 | 3次元形状計測方法および3次元形状計測装置 |
CN102322823B (zh) * | 2011-09-13 | 2013-01-23 | 四川大学 | 基于相位级次自编码的光学三维测量方法 |
US10812694B2 (en) * | 2013-08-21 | 2020-10-20 | Faro Technologies, Inc. | Real-time inspection guidance of triangulation scanner |
CN103575237B (zh) * | 2013-11-12 | 2016-05-18 | 南昌航空大学 | 一种相移编码方法 |
CN104697468B (zh) * | 2015-03-30 | 2017-11-21 | 南昌航空大学 | 基于倾斜相位编码条纹的光学三维测量方法 |
CN106017339B (zh) * | 2016-06-06 | 2019-04-23 | 河北工业大学 | 一种在非完全约束系统中投影非均匀条纹的三维测量方法 |
CN106840036B (zh) * | 2016-12-30 | 2019-10-01 | 浙江四点灵机器人股份有限公司 | 一种适用于快速三维形貌测量的二元结构光优化方法 |
CN107036556B (zh) * | 2017-04-28 | 2019-12-13 | 南昌航空大学 | 基于分段量化相位编码的结构光三维测量方法 |
CN107339954B (zh) * | 2017-05-23 | 2019-09-13 | 南昌航空大学 | 基于周期不同步正弦条纹加相位编码条纹的三维测量方法 |
CN107367244B (zh) * | 2017-07-21 | 2019-06-25 | 南京理工大学 | 一种基于时间相位解缠的设计最优条纹序列方法 |
JP7090446B2 (ja) * | 2018-03-22 | 2022-06-24 | 株式会社キーエンス | 画像処理装置 |
CN109186476B (zh) * | 2018-10-26 | 2021-01-26 | 广东工业大学 | 一种彩色结构光三维测量方法、装置、设备及存储介质 |
CN109579742B (zh) * | 2018-11-26 | 2020-05-29 | 南昌航空大学 | 一种s型分段相位编码结构光三维测量方法 |
-
2020
- 2020-01-13 CN CN202010032926.4A patent/CN111207694B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111207694A (zh) | 2020-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111207694B (zh) | 一种双步相移法结合相位编码的三维测量方法 | |
CN107036556B (zh) | 基于分段量化相位编码的结构光三维测量方法 | |
CN112880589B (zh) | 基于双频相位编码的光学三维测量方法 | |
CN108955574B (zh) | 一种三维测量方法及系统 | |
CN110645919B (zh) | 一种基于空域二值编码的结构光三维测量方法 | |
CN102322823B (zh) | 基于相位级次自编码的光学三维测量方法 | |
CN110849290B (zh) | 基于形态学操作的分段量化编码强度的三维测量方法 | |
CN111174730B (zh) | 一种基于相位编码的快速相位解包裹方法 | |
CN111207692B (zh) | 一种改进的分段阶梯相位编码三维测量方法 | |
CN108955571B (zh) | 双频外差与相移编码相结合的三维测量方法 | |
CN107339954B (zh) | 基于周期不同步正弦条纹加相位编码条纹的三维测量方法 | |
CN109579742B (zh) | 一种s型分段相位编码结构光三维测量方法 | |
CN110174079B (zh) | 一种基于四步相移编码型面结构光的三维重建方法 | |
CN113514009B (zh) | 一种移位阶梯相位码与相移条纹非对称组合三维测量方法 | |
An et al. | The absolute phase retrieval based on the rotation of phase-shifting sequence | |
Wang et al. | A 3D shape measurement method based on novel segmented quantization phase coding | |
An et al. | Temporal phase unwrapping based on unequal phase-shifting code | |
CN112504165A (zh) | 一种基于双边滤波优化的复合立体相位展开方法 | |
Zheng et al. | Joint coding strategy of the phase domain and intensity domain for absolute phase retrieval | |
CN112146596A (zh) | 一种最优量化相位编码三维测量方法 | |
CN114152203B (zh) | 基于相位编码结构光的轴承内外径尺寸测量方法 | |
Zhang et al. | Three-dimensional shape measurement based on speckle-embedded fringe patterns and wrapped phase-to-height lookup table | |
CN108534715A (zh) | 一种应用于不同波长三维测量的展开相位解相方法 | |
WO2023236725A1 (zh) | 一种三维测量方法、设备和存储介质 | |
CN112945089A (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 |