CN108692677B - 基于反射光场偏振与空时编码的表面三维重建方法和系统 - Google Patents
基于反射光场偏振与空时编码的表面三维重建方法和系统 Download PDFInfo
- Publication number
- CN108692677B CN108692677B CN201810321331.3A CN201810321331A CN108692677B CN 108692677 B CN108692677 B CN 108692677B CN 201810321331 A CN201810321331 A CN 201810321331A CN 108692677 B CN108692677 B CN 108692677B
- Authority
- CN
- China
- Prior art keywords
- shape
- light field
- normal vector
- dimensional
- data
- 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.)
- Expired - Fee Related
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于反射光场偏振与空时编码的表面三维重建方法和系统,能够有效降低累计误差,实现复杂表面的高分辨高精度三维重建,且配置灵活,成本较低。该方法包括:对待测物体进行可调控精度的空时编码三维重建以获取第一三维面形数据Dshape,并计算其法向量场Nshape;采用圆偏振光照明待测物体表面,计算反射光场法向量场Npolar;通过Nshape纠正反射光场法向量分量入射方位角φ,得到第一优化的反射光场法向量场Nc polar;通过Dshape纠正反射光场法向量相关的折射率nm,得到第二优化的反射光场法向量场Nrpolar;基于圆偏振度和Dshape联合引导梯度积分得到第二三维面形数据Dpolar;将Dshape和Dpolar进行数据融合,得到待测物体表面的第三三维面形数据Dhh。
Description
技术领域
本发明涉及光学三维测量与检测技术领域,尤其涉及一种基于反射光场偏振与空时编码的表面三维重建方法和系统。
背景技术
精密测量与检测技术一直是设备制造、加工的重要工具,贯穿到产品的整个生产加工、质量控制、零件装配、缺陷检测等各个环节,是高精度几何零部件和复杂曲面加工的依据和保证,也是工业制造技术发展的基础和前提。三维测量技术作为现代非接触精密测量与检测技术不可缺少的重要手段之一,被广泛应用于汽车制造、航空航天、逆向工程等重要领域,在推动经济发展和社会进步方面起着举足轻重的作用。随着第四次工业革命的到来,如何在复杂环境下,以较低的成本实现复杂表面高分辨高精度三维测量成为新一代精密测量与检测技术发展的重要方向。
被动三维测量基于自然光照明的二维图像分析亮度、阴影、纹理、像流、几何特征等各种三维线索重建待测物体三维面形或者深度数据,主要应用于精度要求不高的场合。其中,光度立体视觉、阴影三维恢复这类方法在可控光照和满足某些假设的条件下,通过解算表面法向量场来重建面形,具有很高的局部分辨率,可以检测到亚毫米甚至更高的几何特征,但获得是相对三维面形数据,在重建模型不合理或者解算方法适应性欠佳等情况下,容易产生形貌失真。
基于结构光照明的主动三维测量技术通过投射激光或者其它形式的点、线、面的结构光,分析受待测物体表面调制的光场图像的相位、强度、调制度或者对比度等信息进行三维重建。与被动三维测量技术相比,主动三维测量具有很高的测量精度和全场形貌还原能力,然而其局部测量分辨率有限,一般只能达到1mm左右,对于具有丰富信息的更高频细节难以探测。采用干涉、全息等原理的测量技术主要面向微纳结构的高分辨三维重建,尽管分辨率可以达到纳米级,但其测量范围小,不适用于宏观尺度的复杂表面三维测量。
目前,从单一信息源实现宏观尺度的高分辨高精度三维重建很难突破其技术瓶颈。将光度立体、阴影形状恢复这类技术与深度数据结合以实现复杂表面高分辨高精度无形貌扭曲的三维重建是一种行之有效的解决途径。然而,光度立体视觉要求估计光源的信息、假设朗伯照明模型等苛刻条件,导致用于实际的复杂环境的精密测量与检测当中仍然面累计误差大、成本高等临诸多问题。
发明内容
本发明的目的之一至少在于,针对如何克服上述现有技术存在的问题,提供一种基于反射光场偏振与空时编码的表面三维重建方法和系统,能够有效降低累计误差,实现复杂表面的高分辨高精度三维重建,且配置灵活,成本较低。
为了实现上述目的,本发明采用的技术方案包括以下各方面。
一种基于反射光场偏振与空时编码的表面三维重建方法,其包括:
对待测物体进行可调控精度的空时编码三维重建以获取第一三维面形数据Dshape,并计算其法向量场Nshape;采用圆偏振光照明待测物体表面,计算反射光场法向量场Npolar;通过第一三维面形数据法向量场Nshape纠正反射光场法向量分量入射方位角得到第一优化的反射光场法向量场Nc polar;通过第一三维面形数据Dshape纠正反射光场法向量相关的折射率nm,得到第二优化的反射光场法向量场Nr polar;基于圆偏振度和Dshape联合引导梯度积分得到第二三维面形数据Dpolar;基于标准球拟合建立异源数据的定量关系,并将Dshape和Dpolar进行数据融合,得到待测物体表面的第三三维面形数据Dhh。
一种基于反射光场偏振与空时编码的表面三维重建系统,其包括:一个空时编码光场发生装置、一个圆偏振光发生装置、两台工业相机,以及计算装置;其中,每台工业相机包括偏振滤光片和相位延迟片;所述计算装置用于执行前述任一方法。
综上所述,由于采用了上述技术方案,本发明至少具有以下有益效果:
通过分析待测物体表面对入射光场偏振特性的作用机制,采用圆偏振光场照明待测表面,很好地消除了法向量场入射角的不确定性;通过空时编码光场调制的三维面形重建技术,获得可调控精度的三维面形数据,满足不同高频细节待测物体表面的三维重建需求,为消除偏振法向量场的不确定性和深度优化提供依据;采用空时编码三维面形数据和偏振度可靠性排序联合引导的梯度积分算法,有效降低了传统梯度积分方法引入的累计误差,从而实现复杂表面高分辨(例如分辨率达到0.1mm)高精度(例如精度达到0.01mm)的三维重建;而且,系统配置灵活性强,无需严格同步,硬件成本较低。
附图说明
图1是根据本发明实施例的基于反射光场偏振与空时编码的表面三维重建方法的流程图。
图2是根据本发明实施例获取第一三维面形数据的流程图。
具体实施方式
下面结合附图及实施例,对本发明进行进一步详细说明,以使本发明的目的、技术方案及优点更加清楚明白。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
根据本发明实施例的基于反射光场偏振与空时编码的表面三维重建系统包括:一个空时编码光场发生装置、一个圆偏振光发生装置、两台黑白工业相机(包括偏振滤光片和相位延迟片),以及计算装置。其中,空时编码光场发生装置的产生的投影与两台相机无需严格同步。该方法通过圆偏振光场照明待测物体表面,以反射光场的偏振特性为信息载体提取法向量场,以可调控精度的空时编码三维重建结果为依据纠正反射光场法向量场,并建立最小生成树路径引导反射光场梯度积分实现高分辨三维重建,最后将二者的三维重建结果进行数据融合,实现高分辨高精度复杂表面的三维重建。
图1示出了根据本发明实施例的基于反射光场偏振与空时编码的表面三维重建方法的流程图。该实施例的方法包括以下步骤:
步骤101:对待测物体进行可调控精度的空时编码三维重建以获取第一三维面形数据Dshape,并计算其法向量场Nshape。
步骤102:采用圆偏振光照明待测物体表面,计算反射光场法向量场Npolar
具体地,可以利用斯托克斯(Stokes)矢量表征偏振光的光强和偏振态。当一束用斯托克斯向量表示的入射光Sin=[S0 S1 S2 S3]T,其与x轴成夹角α,经过与快轴和x轴平行的相位延迟片后相位延迟为其光强分布表示为:
I(α,φ)=0.5[S0+S1cos(2α)cosφ+S2sin(2α)cosφ-S3sin(2α)sinφ]
其中,S0表示入射光的总光强;S1表示偏振方向为0度和90度方向光强差;S2为45度和135度方向的光强差;S3表示左圆偏光和右旋圆偏光的光强差。
则入射光场的线偏振度DOLP(0~1表示非线性偏振~完全线性偏振),圆偏振度DOCP以及偏振角AOP(表示偏振椭圆主轴与参考轴之间的夹角)可以表示为:
假设照明光场从空气(折射率nair≈1)中以入射角θin照明折射率为nm的待测物体表面,折射角为θre,则根据折射定律有:
nairsinθin=nmsinθre
根据菲涅尔定律,入射光场经过待测物体表面的调制,其偏振态将会发生改变,入射光场偏振态的变化采用穆勒矩阵Mreflect(θin,θre)表征:
为了计算反射光场的Stokes向量,引入一个关于θ的旋转矩阵Mrotate(θ)将入射光场的Stokes向量变换到待测物体表面切平面坐标系:
联立穆勒矩阵Mreflect(θin,θre)和旋转矩阵Mrotate(θ),建立入射光场和反射光场的偏振特性关系:
Sout=[Mrotate(-θ)Mreflect(θin,θre)Mrotate(θ)]Sin
某一像素点p的归一化法向量为:
进一步通过计算包括每一像素点的归一化法向量的集合,获得反射光场法向量场Npolar。
示例性地,建立优化目标函数:
通过操作符B[·]对已经基于圆偏振光场纠正的反射光场法向量Npolar进行优化,同时引入平滑控制因子χ,然后基于优化的结果对仅有不确定性法向量场Npolar进行纠正,获取第一优化的反射光场法向量场Nc polar:
Nc polar=B(Npolar)
步骤104:通过第一三维面形数据Dshape纠正反射光场法向量相关的折射率nm,得到第二优化的反射光场法向量场Nr polar
上述步骤103中反射光场法向量场经过不确定性纠正后Nc ploar仅仅消除了符号不确定性,因此可以进一步通过第一三维面形数据Dshape求解的法向量场Mask⊙Nshape,建立一个以旋转矩阵R为优化参数的目标函数,深度优化Nc ploar:
Nc polar=R(Nc polar)
基于上述第一优化的反射光场法向量场Nc polar并结合圆偏振光场照明模型求解入射角θin,最后对偏振法向量场进行更新:Nr polar=Nc polar。进而得到第二优化的反射光场法向量场Nr polar,作为最终的偏振法向量场。
步骤105:基于圆偏振度和Dshape联合引导梯度积分得到第二三维面形数据Dpolar
示例性地,根据三维重建过程,基于Dshape相关系数δcorr和圆偏振度DOCP联合归一化可靠性排序,引导定义一个2D积分路径的最小生成树图形Ggraph,其节点对应重建表面的每个像素点,其边缘对应梯度,每个节点均可以沿着某积分路径到达,然后根据Ggraph进行梯度积分得到Dpolar,其具有很高的局部分辨率,其单位为像素。
步骤106:基于标准球拟合建立异源数据的定量关系,并将Dshape和Dpolar进行数据融合,得到待测物体表面的第三三维面形数据Dhh。
具体地,采用一个经检验的已知面形数据分布的标准球作为测量对象,分别用上述步骤测得其面形数据和通过选取测量可靠的区域的数据建立两种异源数据的比例因子s,将Dshape和Dpolar统一到同一尺度的坐标系下:
Maskshape{0,1},Maskpolar∈{0,1}
去掉冗余的数据并进行重新采样,从而得到待测物体表面的高分辨高精度三维重建数据Dhh。
图2示出了根据本发明实施例进行可调控精度的空时编码三维重建以获取第一三维面形数据Dshape的流程图。其包括如下步骤:
步骤201:获取多对空域时域不相关且受待测物体表面形状调制的立体图像序列。具体地,可以基于双目立体视觉原理,通过两台工业相机对待测物体表面进行拍摄,获取多对空域时域不相关的立体图像序列,这些图像受待测物体表面形状调制而变形。例如,可以通过调整立体图像序列的数量,以及系统三维重建参数(例如,延迟角度、偏振角度、偏振图像数量等)来适应不同复杂面形和不同几何特征物体高分辨高精度三维重建。其中,可以利用反射楔形镜调制固定模式的二元编码图形以在测量空间产生空时无关连续扫描光场。
步骤202:对立体图像序列进行二值化处理,示例性地,可以将每个序列匹配窗口内所有像素灰度值的平均值作为阈值进行二值化计算,从而得到二值立体图像序列对。
步骤203:对获取的二值立体图像序列对进行异或处理,完成对应点立体匹配,得到反映待测物体表面面形的视差图;
步骤204:根据系统的标定参数,对获得的视差图进行三维重建,得到第一三维面形数据Dshape,其表示了待测物体绝对整体形貌数据,单位为毫米。
在各种实施例中,可以采用最近邻搜索算法计算第一三维面形数据Dshape的法向量场Nshape。具体地,可以根据受待测物体调制的空时序列图像和立体匹配的双重选择标准,选取第一三维面形数据Dshape的高可靠性区域Mask;通过最近邻kNN搜索算法找出Mask区域三维点pi=[pix,piy,piz]T邻域的点集P={p1,p2,p3,...,pn},pi∈R3,对应矩阵P=[p1,p2,…,pn]T。对点pi,从其k邻域Qi={qi1,qi2,qi3,...,qik}(qik∈P且qik≠pi)估算该点的法向量ni=[nix,niy,niz]T。Qi对应的矩阵形式为Qi=[qi1,qi2,…,qik]T,
以ni为优化变量,建立一个代价函数:
以上所述,仅为本发明具体实施方式的详细说明,而非对本发明的限制。相关技术领域的技术人员在不脱离本发明的原则和范围的情况下,做出的各种替换、变型以及改进均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于反射光场偏振与空时编码的表面三维重建方法,其特征在于,所述方法包括:
对待测物体进行可调控精度的空时编码三维重建以获取第一三维面形数据Dshape,并计算其法向量场Nshape;采用圆偏振光照明待测物体表面,计算反射光场法向量场Npolar;通过第一三维面形数据法向量场Nshape纠正反射光场法向量分量入射方位角得到第一优化的反射光场法向量场Ncpolar;通过第一三维面形数据Dshape纠正反射光场法向量相关的折射率nm,得到第二优化的反射光场法向量场Nrpolar;基于圆偏振度和Dshape联合引导梯度积分得到第二三维面形数据Dpolar;基于标准球拟合建立异源数据的定量关系,并将Dshape和Dpolar进行数据融合,得到待测物体表面的第三三维面形数据Dhh。
2.根据权利要求1所述的方法,其特征在于,所述获取第一三维面形数据Dshape包括:
获取多对空域时域不相关且受待测物体表面形状调制的立体图像序列;对立体图像序列进行二值化处理,得到二值立体图像序列对;对获取的二值立体图像序列对进行异或处理,完成对应点立体匹配,得到反映待测物体表面面形的视差图;根据系统的标定参数,对获得的视差图进行三维重建,得到第一三维面形数据Dshape。
3.根据权利要求2所述的方法,其特征在于,计算第一三维面形数据Dshape的法向量场Nshape包括:
根据受待测物体调制的空时序列图像和立体匹配的双重选择标准,选取第一三维面形数据Dshape的高可靠性区域Mask;通过最近邻kNN搜索算法找出Mask区域三维点pi=[pix,piy,piz]T邻域的点集P={p1,p2,p3,...,pn},pi∈R3,对应矩阵P=[p1,p2,…,pn]T;对点pi,从其k邻域Qi={qi1,qi2,qi3,...,qik}(qik∈P且qik≠pi)估算该点的法向量ni=[nix,niy,niz]T;Qi对应的矩阵形式为Qi=[qi1,qi2,…,qik]T,
以ni为优化变量,建立一个代价函数:
4.根据权利要求3所述的方法,其特征在于,所述方法进一步包括:
I(α,φ)=0.5[S0+S1 cos(2α)cosφ+S2 sin(2α)cosφ-S3 sin(2α)sinφ]
其中,S0表示入射光的总光强;S1表示偏振方向为0度和90度方向光强差;S2为45度和135度方向的光强差;S3表示左圆偏光和右旋圆偏光的光强差。
6.根据权利要求5所述的方法,其特征在于,所述方法包括:
基于Dshape相关系数δcorr和圆偏振度DOCP联合归一化可靠性排序,引导定义一个2D积分路径的最小生成树图形Ggraph,其节点对应重建表面的每个像素点,其边缘对应梯度,每个节点均可以沿着某积分路径到达,然后根据Ggraph进行梯度积分得到Dpolar。
8.根据权利要求1至7任一项所述的方法,其特征在于,所述方法包括:根据待测物体表面几何特征和面形复杂度调整立体图像序列的数量以及三维重建参数。
9.一种基于反射光场偏振与空时编码的表面三维重建系统,其特征在于,所述系统包括:一个空时编码光场发生装置、一个圆偏振光发生装置、两台工业相机,以及计算装置;其中,每台工业相机包括偏振滤光片和相位延迟片;
所述计算装置用于执行前述任一项权利要求所述的方法。
10.根据权利要求9所述的系统,其特征在于,所述系统包括反射楔形镜,用于调制固定模式的二元编码图形以在测量空间产生空时无关连续扫描光场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810321331.3A CN108692677B (zh) | 2018-04-11 | 2018-04-11 | 基于反射光场偏振与空时编码的表面三维重建方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810321331.3A CN108692677B (zh) | 2018-04-11 | 2018-04-11 | 基于反射光场偏振与空时编码的表面三维重建方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108692677A CN108692677A (zh) | 2018-10-23 |
CN108692677B true CN108692677B (zh) | 2020-03-13 |
Family
ID=63845636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810321331.3A Expired - Fee Related CN108692677B (zh) | 2018-04-11 | 2018-04-11 | 基于反射光场偏振与空时编码的表面三维重建方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108692677B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109357632B (zh) * | 2018-12-26 | 2020-09-15 | 河北工业大学 | 一种镜面物体三维形貌测量方法及装置 |
CN110021067B (zh) * | 2019-03-22 | 2022-12-06 | 嘉兴超维信息技术有限公司 | 一种基于镜面反射梯度偏振光构建三维人脸法线的方法 |
CN109884665A (zh) * | 2019-03-27 | 2019-06-14 | 北京理工大学 | 偏振成像入射角歧义解消除方法及应用 |
CN113029040B (zh) * | 2021-02-26 | 2022-09-27 | 中国科学院深圳先进技术研究院 | 一种偏振相位偏折测量方法和装置 |
CN113108720B (zh) * | 2021-04-06 | 2023-07-21 | 北京工业大学 | 一种基于线偏振光和条纹反射的表面三维重建方法 |
CN113706693B (zh) * | 2021-08-27 | 2023-11-17 | 武汉大学 | 一种微光条件下偏振三维重建方法 |
CN117523112B (zh) * | 2024-01-05 | 2024-05-17 | 深圳市宗匠科技有限公司 | 三维模型建立方法及其系统、控制设备及存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3553652B2 (ja) * | 1994-08-19 | 2004-08-11 | 茂樹 小林 | 形状計測装置、検査装置、及び製品の製造方法 |
CN105571518B (zh) * | 2016-01-19 | 2017-10-17 | 大连理工大学 | 基于折射图像偏差的三维信息视觉测量方法 |
CN105806257B (zh) * | 2016-03-12 | 2019-02-22 | 上海大学 | 一种高反射物体表面光场偏折术测量系统和方法 |
CN106683188A (zh) * | 2016-11-17 | 2017-05-17 | 长春理工大学 | 透明目标的双表面三维重构方法、装置和系统 |
CN107657604A (zh) * | 2017-09-06 | 2018-02-02 | 西安交通大学 | 一种基于近场非标准光源的磨痕三维形貌原位获取方法 |
-
2018
- 2018-04-11 CN CN201810321331.3A patent/CN108692677B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN108692677A (zh) | 2018-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108692677B (zh) | 基于反射光场偏振与空时编码的表面三维重建方法和系统 | |
WO2021027719A1 (zh) | 一种基于反射镜的条纹投影系统标定方法 | |
CN110487216B (zh) | 一种基于卷积神经网络的条纹投影三维扫描方法 | |
Yamazaki et al. | Simultaneous self-calibration of a projector and a camera using structured light | |
CN112053432B (zh) | 一种基于结构光与偏振的双目视觉三维重建方法 | |
US8836766B1 (en) | Method and system for alignment of a pattern on a spatial coded slide image | |
CN111028295A (zh) | 一种基于编码结构光和双目的3d成像方法 | |
CN113108720B (zh) | 一种基于线偏振光和条纹反射的表面三维重建方法 | |
US20210254968A1 (en) | Method and System for Automatic Focusing for High-Resolution Structured Light 3D Imaging | |
EP3631757B1 (en) | Deriving topology information of a scene | |
CN116295113A (zh) | 一种融合条纹投影的偏振三维成像方法 | |
Liu et al. | Accurate measurement method for tube’s endpoints based on machine vision | |
Wang et al. | Pixel-wise phase unwrapping with adaptive reference phase estimation for 3-D shape measurement | |
Wang et al. | Light-field-assisted phase unwrapping of fringe projection profilometry | |
CN115661226B (zh) | 镜面物体的三维测量方法、计算机可读存储介质 | |
JP2004537732A (ja) | 干渉縞を投影し絶対位相写像を評価することによる三次元画像法 | |
CN117053716A (zh) | 一种圆形孔径干涉图轮廓参数自动检测方法 | |
CN116824047A (zh) | 一种基于条纹图像增强的结构光三维重建方法 | |
Wong et al. | 3D model reconstruction by constrained bundle adjustment | |
CN114076579A (zh) | 一种基于偏振成像的三维粗糙度检测装置及方法 | |
Okazaki et al. | A projector-camera system for high-quality synthesis of virtual reflectance on real object surfaces | |
Xu et al. | Global calibration method of a camera using the constraint of line features and 3D world points | |
CN112325799A (zh) | 一种基于近红外光投影的高精度三维人脸测量方法 | |
Li et al. | Combining encoded structured light and photometric stereo for underwater 3D reconstruction | |
CN116486120B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200313 |