CN111473745A - 一种基于多频相移方案的发光表面微观三维测量方法 - Google Patents
一种基于多频相移方案的发光表面微观三维测量方法 Download PDFInfo
- Publication number
- CN111473745A CN111473745A CN202010577091.0A CN202010577091A CN111473745A CN 111473745 A CN111473745 A CN 111473745A CN 202010577091 A CN202010577091 A CN 202010577091A CN 111473745 A CN111473745 A CN 111473745A
- Authority
- CN
- China
- Prior art keywords
- phase
- phase shift
- light
- emitting surface
- frequency
- 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.)
- Withdrawn
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/2545—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 with one projection direction and several detection directions, e.g. stereo
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
本发明公开了一种基于多频相移技术的发光目标微观三维测量方法。本发明在没有饱和的情况下,使用标准相移算法来计算相位值;在饱和区域,使用广义相移算法来计算包裹相位;对于非饱和强度小于3的过饱和区域,利用低频条纹图像中可能提取的相位来填充最终相位图,以提高测量的完整性。经过相位展开和立体匹配后的双视远心测量***,实现高精度的发光表面的三维重建。本发明从不受光强饱和影响的相移条纹图像的子集计算突出显示区域的相位,提出了一种多频相移方案来提高发光表面最终相位图的完整性,并在此基础上结合显微远心立体视觉***,实现了完整、高精度的三维重建。
Description
技术领域
本发明属于光学测量技术领域,具体涉及一种基于多频相移方案的发光表面微观三维测量方法。
背景技术
结构光和三角测量原理在三维光学计量领域得到了广泛的应用。周期性的正弦条纹图案被投射到被检物体上,被检物体的调制使条纹发生畸变。为了定量地计算调制量并重建目标的三维结果,需要精确地检索在条纹图中编码的相位值。目前,两种常用的相位恢复算法是基于傅里叶变换的算法和基于相移的算法。基于傅里叶变换的动态测量算法是动态测量中常用的算法,而基于相移的动态测量算法由于其像素独立的数学运算特性,更适合于高精度测量。最近的工作表明,采用基于相位的立体匹配方法,可以忽略数字投影仪的本质非线性响应函数,因为不同视图中的相位误差是自动消除的。然而,基于相位的立体匹配方法在处理有光泽表面的物体时容易失败。重建模型的整体性受高光区域的影响,因为无法用密集的条纹图像计算这些区域的相位。
有光泽的表面具有很强的反射性,因此由于数码相机的动态范围有限,光线强度不能线性变换。针对这种情况,最先进的技术之一被称为高动态范围三维形状测量,可分为两类:基于设备的技术和基于算法的技术。对于这组基于设备的技术,设备的最佳参数,例如相机的曝光时间或投影仪的曝光时间是被需要的,以帮助在光亮和黑暗区域形成可见条纹。其他基于光学的方法,例如使用偏振器扫描发光物体,也已被研究,在此基础上可以有效地抑制偏振光强度。此外,也有通过修改相机曝光量的混合方法,但也考虑了引入额外设备、改变观看位置或调整投影仪参数以捕获高速动态范围图像的策略。基于最大强度调制,提出了一种采用高速投影仪的快速高速动态范围解决方案,在700Hz下对光强变化的条纹图像进行投影。
然而,对于有光泽的表面,有时仅仅通过减少曝光时间或投影光的强度,可能不容易解决饱和问题。因此,研究人员还开发了基于算法的技术,主要依靠精心设计的算法,在相机或投影仪曝光时间不允许自由变化或没有附加设备的情况下,从原始条纹图像中提取相位值。
然而,在显微成像中,由于显微投影***的景深较短,被黑色条纹照亮的发光部分不再是纯黑色的,而是会受到白色条纹的影响。在这种情况下,当使用更高频率的条纹时,会有更多的饱和区域。
发明内容
本发明提出了一种基于多频相移方案的发光表面微观三维测量方法。
本发明的具体技术方案如下:
一种基于多频相移方案的发光表面微观三维测量方法,其步骤如下:
步骤一.在发光表面未饱和区域,使用标准相移算法来计算相位值,在发光表面饱和区域,使用广义相移算法来计算包裹相位;
步骤二.对于非饱和强度小于3的过饱和区域,利用低频条纹图像中可能提取的相位来填充最终相位图,以提高测量的完整性;
步骤三.经过相位展开和立体匹配后的双视远心测量***,实现高精度的发光表面的三维重建。
优选的,在步骤一中,基于可控相移量,记录到的相移条纹图表示为式(1),
其中,
优选的,步骤二具体如下:
步骤2.1.考虑到不同条纹周期的饱和程度,在相位展开阶段将参考存储的信息,用于计算图像中每个像素的饱和强度;
步骤2.2.应用对应于方程的广义相移算法,通过一般相移来对部分饱和相移条纹图像的相位进行计算;
步骤2.3.自动融合,用于通过多频率高动态范围保证相位解包裹的正确性。
优选的,步骤三具体如下:
步骤3.4.在立体匹配中完成左右一致性检查后,得到匹配的像素对,实现高精度的发光表面的三维重建。
本发明与现有技术相比,其显著优点:显微条纹投影轮廓术是一种强大的三维测量技术,理论测量精度优于一微米,然而密集条纹的散焦和复杂的表面反射特性通常会导致条纹的光强饱和和质量下降,这使得完整的三维重建变得困难。为了解决这一问题,本发明从不受光强饱和影响的相移条纹图像的子集计算突出显示区域的相位,提出了一种多频相移方案来提高发光表面最终相位图的完整性,并在此基础上结合显微远心立体视觉***,实现了完整、高精度的三维重建。
附图说明
图1为本发明实施例的三维测量方法的流程示意图。
图2为本发明实施例中不同频率条纹饱和度的比较;(a)和(b)未经处理的条纹图像;(c)和(d)放大饱和部分的细节;(e)和(f)显示非饱和状态小于3的区域的指数。
图3为本发明实施例中算法1流程图。
图4为本发明实施例中算法2流程图。
图5为本发明实施例中算法3流程图。
图6为本发明实施例处理部分饱和目标时与传统方法的比较。
图7为本发明实施例在机械表的镀镍板上进行的实验示意图。
图8为本发明实施例在机械表的镀镍板上进行实验的数据图。
图9为本发明实施例在机械表的镀镍板上进行实验的另一个数据图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本申请,并不用于限定本申请。
结合附图,本实施例中,基于多频相移方案的发光表面微观三维测量方法,其步骤如下。
步骤一.在发光表面未饱和区域,使用标准相移算法来计算相位值,在发光表面饱和区域,使用广义相移算法来计算包裹相位。
基于可控相移量,记录到的相移条纹图表示为式(1),
其中,
步骤二.对于非饱和强度小于3的过饱和区域,利用低频条纹图像中可能提取的相位来填充最终相位图,以提高测量的完整性。
在传统的多频相移方法中,利用较低频率的条纹图为相位展开提供参考相位图,使用最高频率的条纹图像决定最终的测量精度。实际上,在饱和区,最终的相位值可以用饱和强度较低的低密度条纹图像的相位值来代替。这样,就可以使三维重建尽可能完整。为此,提出了一种基于多频条纹的高速动态范围表面测量方案。通过图像数据预处理、饱和度检测与补偿算法、相位立体匹配等三个步骤,实现了对发光表面的高精度显微三维测量。
第一步是图像数据预处理阶段,包括图像采集、图像校正和根据条纹频率进行分类。条纹图案依次与触发信号一起投影,以实现相机同步。第二步是主要部分,对应于所提出的基于多频条纹的方案,即通过三种算法来计算展开的相位图,分别如附图3、4、5所示。
算法程序见3、4、5所示,算法1用于计算图像中每个像素的饱和强度。考虑到不同条纹周期的饱和程度,在相位展开阶段将参考存储的信息。算法2是通过一般相移来对部分饱和相移条纹图像的相位计算的算法。我们在每个像素处消除无效的强度,并应用对应于方程的广义相移算法。公式用于相位计算。算法3是一种自动融合方法,用于通过多频率高动态范围进行保证相位解包裹的正确性。
步骤三.经过相位展开和立体匹配后的双视远心测量***,实现高精度的发光表面的三维重建。
从数字投影仪依次投影以水平增加的相位图编码的正弦曲线图。利用所提出的基于多频条纹的方法,可以得到两台摄像机的绝对相位值,用于立体匹配。首先对条纹图像进行远心极性校正。在不失一般性的前提下,将左摄像机作为主摄像机。对于左相机上具有相位值的像素,任务是在右图像的第 行中找到对应的像素 。因为条纹方向是垂直的,所以展开的相位值沿水平方向增加。首先获得积分像素 ,该像素在第行中的相位值最接近,相位为。 然后,基于逆线性插值计算亚像素坐标 ,
在立体匹配中完成左右一致性检查后,得到匹配的像素对,实现高精度的发光表面的三维重建。
图6、图7、图8及图9是本实施例的实验图表。具体的,图6为实施例处理部分饱和目标时与传统方法的比较,其中图(a)表示金属表带上的加盖徽标的边缘图像;(e)表示印刷电路板的边缘图像。图7为采用本实施例方法将机械表在镀镍板上进行的实验示意图,其中,图(a)为样品图片,图(b)为样品底面的条纹图像,图(c)为样品顶部的条纹图像,图(d)为底部的重建3D模型,图(e)为顶部的重建3D模型。图8、图9为采用本实施例方法将机械表在镀镍板上进行实验的数据图,其中,()-()数据的剖视图,用(d)中标记的线表示,()-()数据的剖视图,用(e)中标记的线表示。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (5)
1.一种基于多频相移方案的发光表面微观三维测量方法,其特征在步骤如下:
步骤一.在发光表面未饱和区域,使用标准相移算法来计算相位值,在发光表面饱和区域,使用广义相移算法来计算包裹相位;
步骤二.对于非饱和强度小于3的过饱和区域,利用低频条纹图像中可能提取的相位来填充最终相位图,以提高测量的完整性;
步骤三.经过相位展开和立体匹配后的双视远心测量***,实现高精度的发光表面的三维重建。
4.根据权利要求1所述的基于多频相移方案的发光表面微观三维测量方法,,其特征在于,步骤二具体如下:
步骤2.1.考虑到不同条纹周期的饱和程度,在相位展开阶段将参考存储的信息,用于计算图像中每个像素的饱和强度;
步骤2.2.应用对应于方程的广义相移算法,通过一般相移来对部分饱和相移条纹图像的相位进行计算;
步骤2.3.自动融合,用于通过多频率高动态范围保证相位解包裹的正确性。
5.根据权利要求1所述的基于多频相移方案的发光表面微观三维测量方法,,其特征在于,步骤三具体如下:
步骤3.4.在立体匹配中完成左右一致性检查后,得到匹配的像素对,实现高精度的发光表面的三维重建。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010577091.0A CN111473745A (zh) | 2020-06-23 | 2020-06-23 | 一种基于多频相移方案的发光表面微观三维测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010577091.0A CN111473745A (zh) | 2020-06-23 | 2020-06-23 | 一种基于多频相移方案的发光表面微观三维测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111473745A true CN111473745A (zh) | 2020-07-31 |
Family
ID=71763943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010577091.0A Withdrawn CN111473745A (zh) | 2020-06-23 | 2020-06-23 | 一种基于多频相移方案的发光表面微观三维测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111473745A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111998799A (zh) * | 2020-08-20 | 2020-11-27 | 四川大学 | 一种多频投影三维面形测量的运动区域检测方法及*** |
CN112747693A (zh) * | 2020-12-25 | 2021-05-04 | 南京理工大学智能计算成像研究院有限公司 | 针对高反光物体的一种基于彩色图像检测的三维测量方法 |
CN114326352A (zh) * | 2021-12-31 | 2022-04-12 | 南京理工大学智能计算成像研究院有限公司 | 一种基于数字全息的实时细胞三维分析方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101173668B1 (ko) * | 2009-10-27 | 2012-08-20 | 서울대학교산학협력단 | 다중 공간 주파수를 이용한 3차원 물체의 깊이 측정 방법 및 그 장치 |
CN106091981A (zh) * | 2016-05-27 | 2016-11-09 | 西安交通大学 | 用于高动态范围物体的区域投影光学三维轮廓测量方法 |
CN109974624A (zh) * | 2019-03-22 | 2019-07-05 | 华中科技大学 | 一种基于多频相移的减少投影图像数量的方法 |
CN110207614A (zh) * | 2019-05-28 | 2019-09-06 | 南京理工大学 | 一种基于双远心相机匹配的高分辨高精度测量***及方法 |
CN110702034A (zh) * | 2019-10-25 | 2020-01-17 | 湖北工业大学 | 高反光表面三维面形测量方法、服务器及*** |
-
2020
- 2020-06-23 CN CN202010577091.0A patent/CN111473745A/zh not_active Withdrawn
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101173668B1 (ko) * | 2009-10-27 | 2012-08-20 | 서울대학교산학협력단 | 다중 공간 주파수를 이용한 3차원 물체의 깊이 측정 방법 및 그 장치 |
CN106091981A (zh) * | 2016-05-27 | 2016-11-09 | 西安交通大学 | 用于高动态范围物体的区域投影光学三维轮廓测量方法 |
CN109974624A (zh) * | 2019-03-22 | 2019-07-05 | 华中科技大学 | 一种基于多频相移的减少投影图像数量的方法 |
CN110207614A (zh) * | 2019-05-28 | 2019-09-06 | 南京理工大学 | 一种基于双远心相机匹配的高分辨高精度测量***及方法 |
CN110702034A (zh) * | 2019-10-25 | 2020-01-17 | 湖北工业大学 | 高反光表面三维面形测量方法、服务器及*** |
Non-Patent Citations (1)
Title |
---|
Y. HU, Q. CHEN AND Y. LIANG ET AL: "Microscopic 3D measurement of shiny surfaces based on a multi-frequency phase-shifting scheme", 《OPTICS AND LASERS IN ENGINEERING》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111998799A (zh) * | 2020-08-20 | 2020-11-27 | 四川大学 | 一种多频投影三维面形测量的运动区域检测方法及*** |
CN112747693A (zh) * | 2020-12-25 | 2021-05-04 | 南京理工大学智能计算成像研究院有限公司 | 针对高反光物体的一种基于彩色图像检测的三维测量方法 |
CN112747693B (zh) * | 2020-12-25 | 2022-05-31 | 南京理工大学智能计算成像研究院有限公司 | 针对高反光物体的一种基于彩色图像检测的三维测量方法 |
CN114326352A (zh) * | 2021-12-31 | 2022-04-12 | 南京理工大学智能计算成像研究院有限公司 | 一种基于数字全息的实时细胞三维分析方法 |
CN114326352B (zh) * | 2021-12-31 | 2024-06-04 | 南京理工大学智能计算成像研究院有限公司 | 一种基于数字全息的实时细胞三维分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11808564B2 (en) | Calibration method for fringe projection systems based on plane mirrors | |
Feng et al. | General solution for high dynamic range three-dimensional shape measurement using the fringe projection technique | |
Liu et al. | Real-time 3D surface-shape measurement using background-modulated modified Fourier transform profilometry with geometry-constraint | |
Feng et al. | Fast three-dimensional measurements for dynamic scenes with shiny surfaces | |
CN111473745A (zh) | 一种基于多频相移方案的发光表面微观三维测量方法 | |
US9322643B2 (en) | Apparatus and method for 3D surface measurement | |
CN106204732B (zh) | 动态曝光的三维重建方法及*** | |
Dai et al. | A dual-frequency fringe projection three-dimensional shape measurement system using a DLP 3D projector | |
Hu et al. | Microscopic 3D measurement of shiny surfaces based on a multi-frequency phase-shifting scheme | |
CN108195313A (zh) | 一种基于光强响应函数的高动态范围三维测量方法 | |
CN111238403A (zh) | 一种基于光场子孔径条纹图像的三维重建方法及装置 | |
CN109307483A (zh) | 一种基于结构光***几何约束的相位展开方法 | |
CN111288925A (zh) | 一种基于数字对焦结构照明光场的三维重建方法和装置 | |
CN110692084B (zh) | 用于导出场景的拓扑信息的装置和机器可读存储介质 | |
Guo et al. | 3-D shape measurement by use of a modified Fourier transform method | |
CN110595388B (zh) | 一种基于双目视觉的高动态实时三维测量方法 | |
Wang et al. | Motion-induced error reduction for phase-shifting profilometry with phase probability equalization | |
Ma et al. | Real-time 3-D shape measurement based on radial spatial carrier phase shifting from circular fringe pattern | |
Suresh et al. | PMENet: phase map enhancement for Fourier transform profilometry using deep learning | |
KR20140021765A (ko) | 컴퓨터 홀로그램 생성을 위한 가상시점 깊이영상 합성 방법 | |
CN116608794B (zh) | 一种抗纹理3d结构光成像方法、***、装置及存储介质 | |
CN116592794A (zh) | 一种基于偏振结构光的快速三维重建方法 | |
CN116518880A (zh) | 一种基于傅里叶轮廓术的正则化解包裹方法及*** | |
Zhang et al. | Accurate measurement of high-reflective surface based on adaptive fringe projection technique | |
CN112325799A (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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200731 |