CN113532330B - 一种相位格雷码三维测量方法 - Google Patents
一种相位格雷码三维测量方法 Download PDFInfo
- Publication number
- CN113532330B CN113532330B CN202110999148.0A CN202110999148A CN113532330B CN 113532330 B CN113532330 B CN 113532330B CN 202110999148 A CN202110999148 A CN 202110999148A CN 113532330 B CN113532330 B CN 113532330B
- Authority
- CN
- China
- Prior art keywords
- phase
- phase shift
- gray code
- cosine
- code
- 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é
Abstract
一种相位格雷码三维测量方法属于结构光三维测量技术领域。该方法包括以下步骤:首先,生成相位格雷码的四步余弦相移条纹图案;然后,依次投射各位余弦相移条纹图案并同步采集余弦相移条纹图像;第三,根据四步相移法获取余弦相移条纹图像的相位码字;第四,将相位码字转换为相位格雷码字;第五,将相位格雷码字转换为二进制码字;第六,将二进制码转换为十进制码;最后,根据三角法使用十进制码计算被测表面三维坐标、形成三维图像。本发明有益效果在于:相比格雷码三维测量方法,本发明方法具有更强的抗干扰能力,有效地减少了强度干扰带来的解码误差及其导致的测量误差,明显提升了三维成像效果。
Description
技术领域
一种相位格雷码三维测量方法属于结构光三维测量技术领域。
背景技术
图像因信息量大而得到广泛应用。然而,传统的二维图像因损失了深度信息而难以反映出真实的三维客观世界。因此,基于各种原理的三维测量技术应运而生。其中,属于光学非接触三维测量技术的结构光三维测量方法最具实用性和发展前景。
结构光三维测量方法向被测景物表面投射编码光图案,然后采集经被测景物表面调制后的编码光图像,再根据三角法通过图像处理和运算得到被测表面的三维坐标。格雷码其码字由0和1组成而编码简单;该编解码不仅在最低位出错仅导致±1的误差,而且在任意位其跳变点左右不超过最低位码字宽度的范围内出错也仅导致±1的误差,从而保证了错误最小化;另外,它采用时间域编码而不存在译码误差累计问题。因此,格雷码最适合用于解决高度剧烈变化或不连续表面的三维测量难题。
但是,实际中格雷码图像因测量系统硬件噪声、环境干扰以及被测表面物理和几何特性的影响必然存在畸变,其跳变边缘即亮(+1)暗(-1)的边界在物理上也必然存在过渡区,这些干扰因素均给格雷码图像带来强度误差,结果导致格雷码误差和三维测量误差。因此,提升格雷码的抗干扰能力有益于减少格雷码误差及其导致的测量误差、提升三维成像质量。
发明内容
为了提升格雷码的抗干扰能力,本发明公开了一种相位格雷码三维测量方法,相比于现有格雷码三维测量方法,本发明方法将强度转换为相位后再通过相位确定码字,有效地减少了强度干扰带来的解码误差及其导致的测量误差,明显提升了三维成像效果。
本发明的目的是这样实现的:
一种相位格雷码三维测量方法,包括以下步骤:
步骤a、生成n位相位格雷码中第i位的四步余弦相移条纹图案:
其中,i=1,2,…,n;(x,y)为图案像素点坐标;Ii1(x,y)为第一步余弦相移条纹图案像素强度,Ii2(x,y)为第二步余弦相移条纹图案像素强度,Ii3(x,y)为第三步余弦相移条纹图案像素强度,Ii4(x,y)为第四步余弦相移条纹图案像素强度;Ai(x,y)为图案像素平均强度;Bi(x,y)为图案像素调制强度;为第i位的相位码字,且:
其中,Gi(x,y)∈{0,1};
其中,round()为四舍五入取整函数;
步骤f、将二进制码转换为十进制码Dc:
步骤g、根据三角法使用十进制码计算被测表面三维坐标、形成三维图像。
有益效果:相比格雷码三维测量方法,本发明方法抗干扰能力更强,有效地减少了强度干扰带来的解码误差及其导致的测量误差,明显提升了三维成像效果。
附图说明
图1为本发明方法流程图。
图2为第i位相位格雷码的第一步余弦相移条纹图案中一行像素的强度曲线。
图3为第i位相位格雷码的第二步余弦相移条纹图案中一行像素的强度曲线。
图4为第i位相位格雷码的第三步余弦相移条纹图案中一行像素的强度曲线。
图5为第i位相位格雷码的第四步余弦相移条纹图案中一行像素的强度曲线。
图6为第i位相位格雷码中一行相位码字曲线。
图7为第i位相位格雷码中一行相位格雷码字曲线。
图8为实验一中被测平面。
图9为实验一中干扰为0%时本发明方法得到的三维测量结果(不出现干扰导致的测量误差)。
图10为实验一中干扰为25%时本发明方法得到的三维测量结果(不出现干扰导致的测量误差)。
图11为实验一中干扰为49%时本发明方法得到的三维测量结果(从此开始出现干扰导致的测量误差)。
图12为实验一中干扰为0%时格雷码方法得到的三维测量结果(不出现干扰导致的测量误差)。
图13为实验一中干扰为25%时格雷码方法得到的三维测量结果(从此开始出现干扰导致的测量误差)。
图14为实验一中干扰为49%时格雷码方法得到的三维测量结果(出现了更多的干扰导致的测量误差)。
图15为实验二中被测浮雕圆盘。
图16为实验二中干扰为0%时本发明方法得到的三维测量结果(不出现干扰导致的测量误差)。
图17为实验二中干扰为9%时本发明方法得到的三维测量结果(不出现干扰导致的测量误差)。
图18为实验二中干扰为18%时本发明方法得到的三维测量结果(从此开始出现干扰导致的测量误差)。
图19为实验二中干扰为0%时格雷码方法得到的三维测量结果(不出现干扰导致的测量误差)。
图20为实验二中干扰为9%时格雷码方法得到的三维测量结果(从此开始出现干扰导致的测量误差)。
图21为实验二中干扰为18%时格雷码方法得到的三维测量结果(出现了更多的干扰导致的测量误差)。
具体实施方式
下面结合附图对本发明具体实施方式作进一步详细描述。
本发明一种相位格雷码三维测量方法,其具体流程如图1所示,包括以下步骤:
步骤a、生成n位相位格雷码中第i位的四步余弦相移条纹图案:
其中,i=1,2,…,n;(x,y)为图案像素点坐标;Ii1(x,y)为第一步余弦相移条纹图案像素强度,Ii2(x,y)为第二步余弦相移条纹图案像素强度,Ii3(x,y)为第三步余弦相移条纹图案像素强度,Ii4(x,y)为第四步余弦相移条纹图案像素强度;Ai(x,y)为图案像素平均强度;Bi(x,y)为图案像素调制强度;为第i位的相位码字;
第i位的四步余弦相移条纹图案中,各行像素的强度曲线相同,其中一行像素的强度曲线如图2至图5所示;
其中,Ti为第i位余弦相移条纹的周期,且Ti+1=Ti/2;floor()为向下取整函数;
其中,Gi(x,y)∈{0,1};各行的第i位相位格雷码字相同,其中一行相位格雷码字曲线如图7所示;
其中,round()为四舍五入取整函数;
步骤e、将相位格雷码字转换为二进制码字Bi:
B1=G1
步骤f、将二进制码转换为十进制码Dc:
步骤g、根据三角法使用十进制码计算被测表面三维坐标、形成三维图像。
需要说明的是,本发明技术方案所对应的技术领域为结构光三维测量技术领域,就本领域技术人员而言,本发明方法中每一个步骤的具体参数和采用的仪器设备,本领域技术人员能够根据专业知识进行选择和应用,在本发明中没有必要进行具体的数字说明和仪器设备描述;本发明最后一个步骤中根据三角法使用十进制码计算被测表面三维坐标、形成三维图像,属于本领域中众所周知的成熟技术,本领域技术人员完全能够自主实现,本发明已经做到充分公开。
下面分析说明本发明相位格雷码三维测量方法与格雷码三维测量方法相比,原理上具有抗干扰能力强的技术优势。
实际测量中,格雷码三维测量系统存在硬件噪声和环境干扰,被测表面物理和几何特性存在不一致,格雷码图案跳变边缘在物理上存在过渡区,这些干扰因素均给格雷码图像带来强度误差。
现有格雷码三维测量方法因采用强度确定码字而对强度误差敏感,因此容易产生格雷码误差及其导致的测量误差。
本发明所提出的方法先将强度转换为相位后再通过相位确定码字,因为相位抗干扰能力在原理上优于强度,所以本发明对强度误差不敏感,能有效减少强度干扰带来的解码误差及其导致的测量误差。
下面通过数值模拟仿真实验来证明,本发明方法与格雷码方法相比,具有抗干扰能力强的技术优势。
仿真实验采用数值模拟方式进行,实验重复次数为107,实验过程如下:
采用本文发明方法生成标准的一位相位格雷码四步相移条纹图像,然后在图像上施加均匀分布随机干扰ΔI(x)∈[-γ(2Bi),γ(2Bi)],其中误差系数γ表征了像素干扰与其调制强度的比例关系,接着通过解码依次得到干扰后图像的相位码和相位格雷码,再通过与施加干扰前的标准条纹图像的相位格雷码对比来判断是否出错、得到出错概率。
采用格雷码方法生成标准的一位格雷码条纹图像,然后施加同样的干扰,接着通过解码得到干扰后图像的格雷码,再通过与施加干扰前的标准条纹图像的格雷码对比来判断是否出错、得到出错概率。
实验结果参见表1。
表1数值模拟仿真实验中采用本发明方法和格雷码方法的编解码出错概率
由实验结果可知,编码图像强度误差超过±25.9%时格雷码开始出错而超过±50%时相位格雷码才开始出错,而且同一干扰下本发明方法的出错概率明显低于格雷码方法,因此本发明方法的抗干扰能力远高于格雷码三维测量方法。
下面通过仿真测量实验来证明,本发明相位格雷码三维测量方法与格雷码三维测量方法相比,具有抗干扰能力强的技术优势。
所采用的仿真测量实验系统在3DMAX仿真环境中构建,它由分辨率为1024×768像素的投影仪和分辨率为2048×1536像素的相机所组成。
使用仿真测量实验系统实现了本发明方法和格雷码三维测量方法,其中相位格雷码和格雷码两者码字宽度均为2个像素、位数均为9位。
针对平面和浮雕圆盘进行仿真测量,而且在本发明方法中的相移条纹图像和格雷码方法中的编码图像施加同一均匀分布随机干扰ΔI(x)∈{-γ(2Bi),γ(2Bi)}。
需要说明的是,格雷码三维测量方法为结构光三维测量技术领域中的常用方法,本领域技术人员对其充分知晓、能够根据专业知识自主实现,在本发明中没有必要进行具体描述。
仿真测量实验一、平面测量实验
被测平面和采用本发明方法的仿真测量实验结果如图8至图11所示。
采用格雷码方法的仿真测量实验结果如图12至图14所示。
通过比较实验结果可知,本发明方法在干扰达到或超过±49%时出现干扰导致的测量误差,而格雷码方法在干扰达到或超过±25%时就出现干扰导致的测量误差,这验证了采用本发明方法具有更强的抗干扰能力、可得到更准确的测量结果。
仿真测量实验二、浮雕圆盘测量实验
被测浮雕圆盘和采用本发明方法的仿真测量实验结果如图15至图18。
采用格雷码的仿真测量实验结果如图19至图21所示。
通过比较实验结果可知,本发明方法在干扰达到或超过±18%时出现干扰导致的测量误差,而格雷码方法在干扰达到或超过±9%时就出现干扰导致的测量误差,这进一步验证了采用本发明方法具有更强的抗干扰能力、可得到更准确的测量结果。
Claims (1)
1.一种相位格雷码三维测量方法,其特征在于:包括以下步骤:
步骤a、生成n位相位格雷码中第i位的四步余弦相移条纹图案:
其中,i=1,2,…,n;(x,y)为图案像素点坐标;Ii1(x,y)为第一步余弦相移条纹图案像素强度,Ii2(x,y)为第二步余弦相移条纹图案像素强度,Ii3(x,y)为第三步余弦相移条纹图案像素强度,Ii4(x,y)为第四步余弦相移条纹图案像素强度;Ai(x,y)为图案像素平均强度;Bi(x,y)为图案像素调制强度;为第i位的相位码字,且:
其中,Gi(x,y)∈{0,1};
其中,round()为四舍五入取整函数;
步骤f、将二进制码转换为十进制码Dc:
步骤g、根据三角法使用十进制码计算被测表面三维坐标、形成三维图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110999148.0A CN113532330B (zh) | 2021-08-28 | 2021-08-28 | 一种相位格雷码三维测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110999148.0A CN113532330B (zh) | 2021-08-28 | 2021-08-28 | 一种相位格雷码三维测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113532330A CN113532330A (zh) | 2021-10-22 |
CN113532330B true CN113532330B (zh) | 2022-10-04 |
Family
ID=78092191
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110999148.0A Active CN113532330B (zh) | 2021-08-28 | 2021-08-28 | 一种相位格雷码三维测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113532330B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114234849B (zh) * | 2021-12-20 | 2022-07-15 | 广东工业大学 | 一种调制格雷码信息于周期边缘的三维测量方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303229A (zh) * | 2007-05-09 | 2008-11-12 | 哈尔滨理工大学 | 基于边缘格雷码和线移的结构光3d测量技术 |
CN101881605A (zh) * | 2010-06-02 | 2010-11-10 | 南京航空航天大学 | 基于相位编码技术的光学三维测量方法 |
CN102519390A (zh) * | 2011-12-21 | 2012-06-27 | 哈尔滨理工大学 | 三编码周期灰度梯形相移结构光三维信息获取方法 |
CN104075669A (zh) * | 2014-04-14 | 2014-10-01 | 哈尔滨理工大学 | 一种用于复杂表面编码光测量的容错编解码方法 |
CN105890546A (zh) * | 2016-04-22 | 2016-08-24 | 无锡信捷电气股份有限公司 | 基于正交格雷码和线移相结合的结构光三维测量方法 |
CN108332685A (zh) * | 2018-02-05 | 2018-07-27 | 哈尔滨理工大学 | 一种编码结构光三维测量方法 |
WO2018171385A1 (zh) * | 2017-03-24 | 2018-09-27 | 南京理工大学 | 一种基于改进的傅立叶变换轮廓技术的超快三维形貌测量方法及其系统 |
CN108955574A (zh) * | 2018-07-09 | 2018-12-07 | 广东工业大学 | 一种三维测量方法及系统 |
CN109186476A (zh) * | 2018-10-26 | 2019-01-11 | 广东工业大学 | 一种彩色结构光三维测量方法、装置、设备及存储介质 |
CN111207694A (zh) * | 2020-01-13 | 2020-05-29 | 南昌航空大学 | 一种双步相移法结合相位编码的三维测量方法 |
CN111707192A (zh) * | 2020-07-08 | 2020-09-25 | 中国科学院长春光学精密机械与物理研究所 | 正弦相移非对称结合格雷码的结构光编解码方法及装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6373963B1 (en) * | 1998-02-05 | 2002-04-16 | Textile/Clothing Technology Corporation | Systems, methods and computer program for measuring the surface contour of an object |
US20070115484A1 (en) * | 2005-10-24 | 2007-05-24 | Peisen Huang | 3d shape measurement system and method including fast three-step phase shifting, error compensation and calibration |
US20110080471A1 (en) * | 2009-10-06 | 2011-04-07 | Iowa State University Research Foundation, Inc. | Hybrid method for 3D shape measurement |
CN112146596B (zh) * | 2020-08-31 | 2022-01-28 | 南昌航空大学 | 一种最优量化相位编码三维测量方法 |
CN112880589B (zh) * | 2021-01-18 | 2022-04-01 | 南昌航空大学 | 基于双频相位编码的光学三维测量方法 |
-
2021
- 2021-08-28 CN CN202110999148.0A patent/CN113532330B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303229A (zh) * | 2007-05-09 | 2008-11-12 | 哈尔滨理工大学 | 基于边缘格雷码和线移的结构光3d测量技术 |
CN101881605A (zh) * | 2010-06-02 | 2010-11-10 | 南京航空航天大学 | 基于相位编码技术的光学三维测量方法 |
CN102519390A (zh) * | 2011-12-21 | 2012-06-27 | 哈尔滨理工大学 | 三编码周期灰度梯形相移结构光三维信息获取方法 |
CN104075669A (zh) * | 2014-04-14 | 2014-10-01 | 哈尔滨理工大学 | 一种用于复杂表面编码光测量的容错编解码方法 |
CN105890546A (zh) * | 2016-04-22 | 2016-08-24 | 无锡信捷电气股份有限公司 | 基于正交格雷码和线移相结合的结构光三维测量方法 |
WO2018171385A1 (zh) * | 2017-03-24 | 2018-09-27 | 南京理工大学 | 一种基于改进的傅立叶变换轮廓技术的超快三维形貌测量方法及其系统 |
CN108332685A (zh) * | 2018-02-05 | 2018-07-27 | 哈尔滨理工大学 | 一种编码结构光三维测量方法 |
CN108955574A (zh) * | 2018-07-09 | 2018-12-07 | 广东工业大学 | 一种三维测量方法及系统 |
CN109186476A (zh) * | 2018-10-26 | 2019-01-11 | 广东工业大学 | 一种彩色结构光三维测量方法、装置、设备及存储介质 |
CN111207694A (zh) * | 2020-01-13 | 2020-05-29 | 南昌航空大学 | 一种双步相移法结合相位编码的三维测量方法 |
CN111707192A (zh) * | 2020-07-08 | 2020-09-25 | 中国科学院长春光学精密机械与物理研究所 | 正弦相移非对称结合格雷码的结构光编解码方法及装置 |
Non-Patent Citations (6)
Title |
---|
3-D Shape measurement based on the unequal-period combination approach of phase coding and sinusoidal fringe;Fu Ling等;《Journal of Modern Optics》;20191122;第66卷(第21期);第2083-2092页 * |
Phase-shifting profilometry combined with Gray-code patterns projection: unwrapping error removal by an adaptive median filter;Zheng DongLiang等;《Optics express》;20171231;第25卷(第5期);第4700-4713页 * |
Unequal-period combination approach of gray code and phase-shifting for 3-D visual measurement;Yu Shuang等;《Optics Communications》;20190901;第374卷;第97-106页 * |
基于格雷码和多步相移法的双目立体视觉三维测量技术研究;王兵等;《计算机测量与控制》;20180531;第26卷(第05期);第25-28页 * |
基于格雷码和相移算法的动态三维测量技术;浦婷婷;《光学学报》;20210315;第41卷(第05期);第33-43页 * |
彩色结构光三维测量拼接与点云目标识别方法研究;周浩源;《中国优秀硕士学位论文全文数据库(电子期刊)》;20210630;l138-383 * |
Also Published As
Publication number | Publication date |
---|---|
CN113532330A (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110645919B (zh) | 一种基于空域二值编码的结构光三维测量方法 | |
CN109974626B (zh) | 一种基于相移量编码条纹级次的结构光三维测量方法 | |
CN101326549B (zh) | 检测数字图像中的条纹的方法 | |
CN110230997B (zh) | 一种基于改进单调法的阴影区相位噪声校正方法 | |
CN112146596B (zh) | 一种最优量化相位编码三维测量方法 | |
CN110849290B (zh) | 基于形态学操作的分段量化编码强度的三维测量方法 | |
CN112880589B (zh) | 基于双频相位编码的光学三维测量方法 | |
CN113532330B (zh) | 一种相位格雷码三维测量方法 | |
CN111307066B (zh) | 一种分区间处理的相位展开方法 | |
CN114941999A (zh) | 一种用于结构光投影的二进制编码条纹设计方法 | |
CN113514009B (zh) | 一种移位阶梯相位码与相移条纹非对称组合三维测量方法 | |
Cai et al. | An improved gray-level coding method for absolute phase measurement based on half-period correction | |
Chen et al. | 3D measurement method based on S-shaped segmental phase encoding | |
Gaertner et al. | New highly efficient binary codes for structured light methods | |
CN116608794B (zh) | 一种抗纹理3d结构光成像方法、系统、装置及存储介质 | |
CN112802084A (zh) | 基于深度学习的三维形貌测量方法、系统和存储介质 | |
CN115830154B (zh) | 一种基于二倍角相位编码的解包裹方法 | |
CN113345039A (zh) | 三维重建量化结构光相位图像编码方法 | |
CN116182745A (zh) | 一种基于标记点消隐与提取的三维形貌、形变测量方法 | |
CN116774428A (zh) | 一种基于非对称误差扩散核函数的结构光编码方法 | |
Zheng et al. | A universal self-correcting approach for abnormal jump errors in absolute phase retrieval | |
CN113091647B (zh) | 一种用于高反光表面工件的结构光三维扫描方法 | |
CN114279356B (zh) | 一种用于三维测量的灰度条纹图案设计方法 | |
CN114998409A (zh) | 一种自适应结构光测量方法、装置、电子设备及介质 | |
CN111340957B (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 |