CN110702025B - 一种光栅式双目立体视觉三维测量系统及方法 - Google Patents

一种光栅式双目立体视觉三维测量系统及方法 Download PDF

Info

Publication number
CN110702025B
CN110702025B CN201910460042.6A CN201910460042A CN110702025B CN 110702025 B CN110702025 B CN 110702025B CN 201910460042 A CN201910460042 A CN 201910460042A CN 110702025 B CN110702025 B CN 110702025B
Authority
CN
China
Prior art keywords
light
point
dimensional
sensor
measurement model
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
Application number
CN201910460042.6A
Other languages
English (en)
Other versions
CN110702025A (zh
Inventor
孙军华
谢艳霞
何荣芳
高国防
程晓琦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Institute of Chemical Material of CAEP
Original Assignee
Beihang University
Institute of Chemical Material of CAEP
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University, Institute of Chemical Material of CAEP filed Critical Beihang University
Priority to CN201910460042.6A priority Critical patent/CN110702025B/zh
Publication of CN110702025A publication Critical patent/CN110702025A/zh
Application granted granted Critical
Publication of CN110702025B publication Critical patent/CN110702025B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2433Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring outlines by shadow casting

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明涉及一种光栅式双目立体视觉三维测量系统及方法,包括:左右工业相机构建双目立体视觉传感器,左右工业相机又分别与多线激光器组成左右结构光传感器;多线激光器在物体表面投射多条光条,左右工业相机采集物体左右图像;基于Hessian矩阵法,提取左右图像的光条中心点,构建左右特征点集;利用近邻反投影光条中心线约束和极线约束,获得左右特征点集中的匹配点对;基于双目立体视觉测量模型重构匹配点对,左右结构光视觉测量模型重构左右特征点集中未匹配上的光条中心点;融合双目立体视觉测量模型和左右结构光视觉测量模型的三维重构结果,从而获得物体完整的三维形貌。本发明能够更加快速、准确地获得左右图像的匹配点对,并且在遮挡、反光、阴影干扰的情况下,依然能够获得物体表面完整的三维形貌信息。

Description

一种光栅式双目立体视觉三维测量系统及方法
技术领域
本发明涉及三维测量技术领域,特别涉及一种光栅式双目立体视觉三维测量系统及方法。
背景技术
光栅式双目立体视觉三维测量技术能够实时、高精度地获得物体表面的三维形貌信息,被广泛应用于物体表面重构、缺陷检测、逆向工程等领域。
物体表面三维形貌检测技术主要包括两种:接触式检测和非接触式检测。三坐标测量机作为一种常用的接触式检测方法,具有操作复杂、效率低下、容易损伤被测件表面等缺点。而非接触式检测方法主要包括几何光线法、干涉发等。然而这些方法仅适用于光学元件的面型检测。
随着图像处理和计算机技术的快速发展,视觉测量在三维重构领域取得了广泛的应用。郑顺义等人在专利号为201710044195.3的发明专利“一种实时生成模型的手持三维扫描仪”中设计了一种实时生成模型手持三维扫描仪,该手持式三维扫描仪由两个对称放置的相机和激光器构成。该手持式三维扫描仪首先采集图像,然后根据前方交会原理重构物体表面点的三维坐标,并进行物体表面三维点拼接,从而能够获得物体表面的三维数据。熊才华等人在专利号201710141383.8的发明专利“一种基于双目立体视觉的物体三维尺寸简便测量方法”中提出了一种基于双目立体视觉的物体三维尺寸简便测量方法,该方法自动识别拍摄角度,并基于立体视觉技术快速获取物体的三维尺寸。黄剑艇等人在专利号201810572129.8的发明专利“基于红外传感器、智能扫描系统实现三维扫描的方法”中提出了一种结构简单且可扫描出物体表面质感的方法,该方法首先通过红外传感器在物体周围定义一个三维扫描框,然后获得扫描框内的空间结构和颜色信息,并选择相应的模型造型和材质,通过云端处理,最终获得物体三维模型。虽然上述方法能快速获得光滑表面或者复杂程度低的物体的较高三维特征,但是对于一些结构复杂、表面高度起伏变化剧烈的物体三维测量精度低下。雷彦章等人在论文《一种单双目视觉系统结合的三维测量方法》中提出了一种单双目结合的三维测量方法,该方法将两个相机与投影仪分别组成两个基于投影栅相位法的单目,利用单目的测量结果填补双目测量结果中的空洞。该方法虽然能够获得较高精度的物体三维表面轮廓,但是基于投影栅相位法的解相过程复杂、获得匹配点的计算量大。
对于光栅式双目立体视觉三维测量系统而言,其难点在于如何快速获得唯一的匹配点对。谢光辉等人在论文《一种光栅式双目立体视觉传感器光条匹配方法》中提出了一种光条匹配方法。该方法将将光栅式双目立体视觉传感器分解为两个独立的结构光传感器,并利用两个光栅式结构光传感器分别测量不同光条中心点的三维坐标,构成多个光条中心点集,光条中心点集间欧氏距离最小且对应结构光平面相同的作为匹配光条。然而,该方法计算量大,处理速度慢,且存在大量的计算冗余,无法满足实际的测量需求。魏振忠等人在论文《光栅式双目立体视觉传感器光条快速匹配算法》中提出了一种基于空间三维搜索的快速光条匹配方法。该方法将光栅式双目立体视觉传感器当作两个光栅结构光传感器,并基于可见光条组合模式分析,获得参与光条匹配运算的光条及模型,计算各光条中心点集合间的平均距离,距离最小值对应的光条即为匹配光条。虽然该方法减少了计算量,提高了计算效率,但是该方法的可见光条组合模式分析的有效性是基于光条连续且完全可见的前提,若在遮挡干扰或物体表面高度差变化明显的情况下,该方法将会造成大量的误匹配。
综上所述,现有的基于相移光栅的三维测量方法解相过程复杂,计算量大。而现有的光栅式双目立体视觉三维测量方法无法适用于表面高度起伏变化明显物体的测量。此外,在遮挡、反光、阴影干扰的情况下,三维重构结果存在较多空洞,且大量的误匹配导致三维重构精度无法满足工程实践需求。
发明内容
本发明技术解决问题:克服现有技术的不足,提供一种光栅式双目立体视觉三维测量系统及方法,利用近邻反投影光条中心线约束和极线约束进行对应点匹配,融合双目立体视觉测量模型和左右结构光视觉测量模型的三维重构结果,不仅将匹配过程由二维匹配优化到点对点匹配,提高了对应点匹配的效率和准确性,同时在遮挡、反光、阴影干扰的情况下,依然能够获得物体表面完整的三维形貌信息。
为达到上述目的,本发明的技术方案是这样实现的:
本发明的一种光栅式双目立体视觉三维测量系统,包括双目立体视觉传感器和左右结构光传感器、基于Hessian矩阵的光条中心点提取算法模块、基于近邻反投影光条中心点约束和极线约束的特征点匹配算法模块,双目立体视觉传感器和左右结构光传感器的三维重构结果融合模块;其中:
双目立体视觉传感器和左右结构光传感器,包括:左工业相机、右工业相机和一个多线激光器,三者按照左工业相机、多线激光器和右工业相机的顺序摆放,使得多线激光器投射的多根光条均位于左右工业相机的视场中间;其中,左右工业相机组成一个双目立体视觉传感器,而左右工业相机分别与多线激光器构成左右结构光传感器;多线激光器向物体表面投射多条光条,左右工业相机分别采集物体的左右图像;
基于Hessian矩阵的光条中心点提取算法模块:提取左右图像中光条中心点坐标,构建左、右特征点集;
基于近邻反投影光条中心线约束和极线约束的光条中心点匹配算法模块,遍历左特征点集,对于每一个光条中心点,计算对应右图像的近邻反投影光条中心线和极线,将近邻反投影光条中心线与极线的交点作为侯匹配点,若匹配成功,构成匹配点对,并标记右特征点集中相应的两个光条中心点;否则,则标记左特征点集中的该光条中心点;
双目立体视觉传感器和左右结构光传感器的三维重构结果融合模块:基于双目立体视觉测量模型重构匹配点对;采用左结构光视觉测量模型重构左特征点集中被标记的光条中心点;遍历右特征点集,将右特征点集中未被标记的光条中心点带入右结构光视觉测量模型中,获得对应的空间点的三维坐标;并将左右结构光视觉测量模型的重构结果转换到双目立体视觉传感器坐标系下。
所述双目立体视觉传感器和左右结构光传感器的实现如下:
(1)分别对双目立体视觉传感器和左右结构光传感器进行标定,获得左右工业相机的内参、双目立体视觉传感器的结构参数及左右结构光传感器的光平面方程,其中,左工业相机内参Kl,右工业相机内参Kr分别为:
Figure GDA0002681061380000031
其中,αl和αr分别为左工业相机、右工业相机u轴的有效焦距,βl和βr分别为左工业相机、右工业相机v轴的有效焦距,(u0l,v0l)和(u0r,v0r)分别为左工业相机、右工业相机的光学中心,γl和γr分别为左工业相机、右工业相机u轴与v轴不垂直因子;双目立体视觉传感器的结构参数,即旋转矩阵R,平移矢量T为:
Figure GDA0002681061380000032
T=[tx,ty,tz]T
其中,rs(s=1,...,9)为旋转矩阵R的第s个元素,tx,ty和tz分别为平移矢量T的三个分量;左结构传感器光平面方程:
aix+biy+ciz+di=0 i=0,1,...,m
其中,(ai,bi,ci)为第i个左结构传感器光平面的法向量,di为第i个左结构传感器光平面的截距,(x,y,z)为左结构传感器光平面上的三维点,m为多线激光器的投射的光条根数;
右结构传感器光平面方程:
a′ix′+b′iy′+c′iz′+d′i=0 i=0,1,...,m
其中,(a′i,b′i,c′i)为第i个右结构传感器光平面的法向量,d′i为第i个右结构传感器光平面的截距,(x′,y′,z′)为右结构传感器光平面上的三维点;
其中,m为多线激光器投射的光条根数;
(2)多线激光器向物体表面投射多根光条,左右工业相机从不同角度同时采集物体左右图像,左右图像中均包含携带物体表面高度信息的光条。
所述基于Hessian矩阵的光条中心点提取算法模块中,通过求取图像
Hessian矩阵的最大特征向量,获得左右图像中光条的中心点坐标,从而构建左右特征点集。
所述基于近邻反投影光条中心点约束和极线约束的光条中心点匹配算法模块具体实现过程如下:
(1)遍历左特征点集中的光条中心点(ulj,vlj),基于双目立体视觉测量模型,计算光条中心点(ulj,vlj)对应右图像上的极线Ej,具体计算过程如下:
Figure GDA0002681061380000041
其中,F为3×3的基本矩阵,
Figure GDA0002681061380000042
为光条中心点的齐次坐标形式;
(2)基于左结构光视觉测量模型,获得m个三维点坐标Pi=(xi,yi,zi)T,具体计算过程如下:
Figure GDA0002681061380000043
其中,
Figure GDA0002681061380000044
为光条中心点的齐次坐标形式,Kl为左工业相机内参,ρ为比例因子,(ai,bi,ci)为第i个左结构传感器光平面的法向量,di为第i个左结构传感器光平面的截距;
(3)将m个三维点坐标Pi=(xi,yi,zi)T反投影到右图像上,获得在右图像上的m个反投影点坐标(uri,vri);
(4)对于每个反投影点(uri,vri),搜索右特征点集中与该反投影点欧氏距离最小的两个光条中心点,获得对应的近邻反投影中心线li,同时标记右特征点集中的两个光条中心点;
(5)计算极线Ej与m根近邻反投影中心线li的交点(u′ri,v′ri),将得到的m个交点作为光条中心点(ulj,vlj)在右图像上的侯匹配点,构成m对侯匹配点对;
(6)将m个侯匹配点对代入双目立体视觉测量模型中,实现三维重构,获得m个对应的三维坐标P′i=(x′i,y′i,z′i)T
(7)将m个交点(u′ri,v′ri)代入右结构光视觉测量模型中,得到m个三维坐标P″i=(x″i,y″i,z″i)T,具体计算过程如下:
Figure GDA0002681061380000051
其中,
Figure GDA0002681061380000052
为交点的齐次坐标形式,Kr为右工业相机内参,ρ为比例因子,(a′i,b′i,c′i)为第i个右结构传感器光平面的法向量,di′为第i个右结构传感器光平面的截距;
(8)将三维点Pi=(xi,yi,zi)T,P′i=(x′i,y′i,z′i)T,P″i=(x″i,y″i,z″i)T一一对应,构成m个三维点对{Pi,P′i,P″i},计算每个三维点对中两两三维点之间的欧氏距离之和,若最小的欧氏距离之和小于给定的阈值,则光条中心点(ulj,vlj)与对应的交点(u′ri,v′ri)为一对匹配点对;否则,则表明光条中心点(ulj,vlj)在右图像上无匹配点,并在左特征点集中标记该光条中心点(ulj,vlj)。
所述双目立体视觉传感器和左右结构光传感器的三维重构结果融合模块具体实现过程如下:
(1)基于双目立体视觉测量模型重构匹配点对所对应空间点的三维坐标Pb
(2)采用左结构光视觉测量模型重构左特征点集中所有被标记的光条中心点Pl
(3)将右特征点集中所有未被标记的光条中心点带入右结构光视觉测量模型中,获得对应的空间点的三维坐标Pr
(4)由于基于双目立体视觉测量模型的重构结果Pb与基于左结构光视觉测量模型的重构结果Pl均位于双目立体视觉传感器坐标系下,而基于右结构光视觉测量模型的重构结果Pr位于右工业相机坐标系下,需要将基于右结构光视觉测量模型的重构结果转到双目立体视觉传感器坐标系下,得到融合后的三维重构结果FinalCloud,其具体实现过程如下:
Figure GDA0002681061380000061
其中,R和T分别为双目立体视觉传感器的旋转矩阵和平移矢量,Pr′为位于双目立体视觉传感器坐标系下的右结构光视觉测量模型的重构结果。
本发明的一种光栅式双目立体视觉三维测量方法,步骤如下:
(1)多线激光器向物体表面投射多条光条,左右工业相机分别采集物体的左右图像;
(2)通过求取左右Hessian矩阵的最大特征向量,获得左右图像中光条的中心点坐标,构建左、右特征点集;
(3)为了获得准确唯一的匹配,提高三维重构精度,在现有极线约束的基础上增加近邻反投影光条中心点约束,将匹配过程从二维匹配优化为点到点匹配过程,遍历左特征点集中的光条中心点(ulj,vlj),基于双目立体视觉测量模型,计算光条中心点(ulj,vlj)对应右图像上的极线Ej,具体计算过程如下:
Figure GDA0002681061380000062
其中,F为3×3的基本矩阵,
Figure GDA0002681061380000063
为光条中心点的齐次坐标形式;
(4)基于左结构光视觉测量模型,获得m个三维点坐标Pi=(xi,yi,zi)T,具体计算过程如下:
Figure GDA0002681061380000064
其中,
Figure GDA0002681061380000065
为光条中心点的齐次坐标形式,Kl为左工业相机内参,ρ为比例因子,(ai,bi,ci)为第i个左结构传感器光平面的法向量,di为第i个左结构传感器光平面的截距;
(5)将m个三维点坐标Pi=(xi,yi,zi)T反投影到右图像上,获得在右图像上的m个反投影点坐标(uri,vri);
(6)对于每个反投影点(uri,vri),搜索右特征点集中与该反投影点欧氏距离最小的两个光条中心点,获得对应的近邻反投影中心线li,同时标记右特征点集中的两个光条中心点;
(7)计算极线Ej与m根近邻反投影中心线li的交点(u′ri,v′ri),将得到的m个交点作为光条中心点(ulj,vlj)在右图像上的侯匹配点,构成m对侯匹配点对;
(8)将m个侯匹配点对代入双目立体视觉测量模型中,实现三维重构,获得m个对应的三维坐标P′i=(x′i,y′i,z′i)T
(9)将m个交点(u′ri,v′ri)代入右结构光视觉测量模型中,得到m个三维坐标P″i=(x″i,y″i,z″i)T,具体计算过程如下:
Figure GDA0002681061380000071
其中,
Figure GDA0002681061380000072
为交点的齐次坐标形式,Kr为右工业相机内参,ρ为比例因子,(a′i,b′i,c′i)为第i个右结构传感器光平面的法向量,d′i为第i个右结构传感器光平面的截距;
(10)将三维点Pi=(xi,yi,zi)T,P′i=(x′i,y′i,z′i)T,P″i=(x″i,y″i,z″i)T一一对应,构成m个三维点对{Pi,P′i,P″i},计算每个三维点对中两两三维点之间的欧氏距离之和,若最小的欧氏距离之和小于给定的阈值,则光条中心点(ulj,vlj)与对应的交点(u′ri,v′ri)为一对匹配点对,将基于双目立体视觉测量模型重构的三维坐标P′i=(x′i,y′i,z′i)T作为其对应空间点的三维坐标;否则,表明由于遮挡、反光、阴影干扰导致光条中心点(ulj,vlj)在右图像上无匹配点,将基于左结构光视觉测量模型重构的三维点Pi=(xi,yi,zi)T作为其对应空间点的三维坐标;
(11)当左特征点集中所有光条中心点遍历完之后,右特征点集中未被标记的光条中心点带入右结构光视觉测量模型中,获得对应的空间点的三维坐标M;
(12)在遮挡、反光、阴影干扰的情况下,由于双目立体视觉传感器视场的局限性,基于双目立体视觉测量模型的重构结果会存在缺失,因此需要将双目立体视觉传感器与左右结构光传感器的三维重构结果融合;基于双目立体视觉测量模型的重构结果P′与基于左结构光视觉测量模型的重构结果P均位于双目立体视觉传感器坐标系下,而基于右结构光视觉测量模型的重构结果M位于右工业相机坐标系下;需要将基于右结构光视觉测量模型的重构结果转到双目立体视觉传感器坐标系下,从而得到融合后的完整的三维重构结果FinalCloud,其具体实现过程如下:
Figure GDA0002681061380000073
其中,R和T分别为双目立体视觉传感器的旋转矩阵和平移矢量,M′为位于双目立体视觉传感器坐标系下的右结构光视觉测量模型的重构结果。
本发明与现有技术相比的优点在于:
(1)现有的对应点匹配算法,计算量大,存在大量冗余,甚至在遮挡干扰或物体表面高度起伏变化明显的情况下会造成大量的误匹配。而本发明在现有极线约束的基础上增加近邻反投影光条中心点约束,将匹配过程从二维匹配优化为点到点匹配过程,大大提高了对应点匹配的效率,在遮挡、物体表面高度起伏变化明显的情况下依旧能获得准确唯一的匹配。
(2)在遮挡、反光、阴影干扰的情况下,由于双目立体视觉传感器视场的局限性,基于双目立体视觉测量模型的重构结果会存在缺失。现有的三维测量方法主要采用多次扫描填补空洞的方式,耗时且存在大量的冗余数据。本发明在不增加硬件的前提下,将光栅式双目立体视觉传感器等效为一个双目立体视觉传感器和左右结构光传感器,并将基于双目立体视觉测量模型的重构结果与左右结构光视觉测量模型的重构结果融合,从而填补了基于双目立体视觉测量模型的重构结果中的空洞,无需多次扫描,依然能获得物体表面完整的三维形貌信息。
附图说明
图1为本发明一种光栅式双目立体视觉三维测量系统及方法的总体实现流程图;
图2为光栅式双目立体视觉三维测量系统;其中,①为左工业相机,②为右工业相机,③为激光器,④为控制模块,⑤为补光光源,⑥为按键,⑦为待测物体;
图3为基于近邻反投影光条中心点约束和极线约束的光条中心点匹配算法示意图;
图4为金属圆柱筒外表面的局部三维重构结果。
具体实施方式
以下结合附图及具体实施实例对本发明再作进一步详细的说明。
如图1所示,本发明的具体实现包括以下步骤:
步骤1:光栅式双目立体视觉三维测量系统如图2所示,主要包括左工业相机1、右工业相机2、7线激光器3、控制模块4、补光光源5及按键6。左右工业相机按照从左到右的顺序交叉摆放,7激光器3安装在左右3工业相机中间,使得7线激光器3投射的7根光条均位于左右工业相机的视场中间。7线激光器向待测物体7表面投射7根光条,左右工业相机分别采集待测物体7的左右图像。
将左右工业相机构建一个双目立体视觉传感器,左右工业相机又分别与7线激光器组成左右结构光传感器,并分别对双目立体视觉传感器和左右结构光传感器进行标定,获得左右工业相机的内参、双目立体视觉传感器的结构参数及左右结构光传感器的光平面方程,其中,左工业相机内参Kl,右工业相机内参Kr分别为:
Figure GDA0002681061380000091
双目立体视觉传感器的结构参数(旋转矩阵R,平移矢量T)为:
Figure GDA0002681061380000092
T=[-273.26083,-34.86622,115.48466]T
左结构传感器光平面方程:
Figure GDA0002681061380000093
右结构传感器光平面方程:
Figure GDA0002681061380000094
步骤2:7线激光器在金属圆柱筒表面投射7根光条,光条受金属圆柱筒表面高度调制,携带了金属圆柱筒表面的高度信息,左右工业相机采集金属圆柱筒左右图像。
步骤3:基于Hessian矩阵的光条中心点提取算法中,通过求取图像Hessian矩阵的最大特征向量,获得左右图像中光条的中心点坐标,从而构建左右特征点集。具体过程在文章“C.Steger,“Analytical and empirical performance evaluation of subpixel lineand edge detection”,Empir.Eval.Methods Comput.Vis.,pp.188-210,Apr.1998.”中有详细描述。
步骤4:如图3所示,遍历左特征点集中的光条中心点(ulj,vlj),基于双目立体视觉测量模型,计算其对应右图像上的极线Ej;基于左结构光视觉测量模型,获得7个三维点坐标Pi=(xi,yi,zi)T
Figure GDA0002681061380000101
其中,
Figure GDA0002681061380000102
为光条中心点的齐次坐标形式。
将7个三维点坐标Pi=(xi,yi,zi)T反投影到右图像上,获得在右图像上的7个反投影点坐标(uri,vri);若反投影点坐标(uri,vri)超出右图像的范围,则剔除。
对于每个反投影点(uri,vri),根据k近邻算法(k-nearest neighbors,KNN)搜索右特征点集中与该反投影点欧氏距离最小的两个光条中心点,从而获得对应的7根近邻反投影中心线li,同时标记右特征点集中的两个光条中心点。其中,k近邻算法在文章“CoverT.,Hart,P.,“Nearest neighbor pattern classification”,IEEE Trans.Inf.Theor.,pp.21-27,1968.”中有详细介绍。
计算极线Ej与7根近邻反投影中心线li的交点(u′ri,v′ri),将得到的7个交点作为光条中心点(ulj,vlj)在右图像上的侯匹配点,构成7对侯匹配点对;
将7个侯匹配点对代入双目立体视觉测量模型中,实现三维重构,获得7个对应的三维坐标P′i=(x′i,y′i,z′i)T
将7个交点(u′ri,v′ri)代入右结构光视觉测量模型中,得到7个三维坐标P″i=(x″i,y″i,z″i)T,具体计算过程如下:
Figure GDA0002681061380000111
其中,
Figure GDA0002681061380000112
为交点的齐次坐标形式;
将三维点Pi=(xi,yi,zi)T,P′i=(x′i,y′i,z′i)T,P″i=(x″i,y″i,z″i)T一一对应,构成m个三维点对{Pi,P′i,P″i},计算每个三维点对中两两三维点之间的欧氏距离之和,若最小的欧氏距离之和小于给定的阈值4.5,则光条中心点(ulj,vlj)与对应的交点(u′ri,v′ri)为一对匹配点对,将基于双目立体视觉测量模型重构的三维坐标P′i=(x′i,y′i,z′i)T作为其对应空间点的三维坐标;否则,表明光条中心点(ulj,vlj)在右图像上无匹配点,将基于左结构光视觉测量模型重构的三维点Pi=(xi,yi,zi)T作为其对应空间点的三维坐标;
当左特征点集中所有光条中心点遍历完之后,将右特征点集中未被标记的光条中心点带入右结构光视觉测量模型中,获得对应的空间点的三维坐标M。
步骤5:将基于右结构光视觉测量模型的重构结果转到双目立体视觉传感器坐标系下,从而得到融合后的三维重构结果FinalCloud。金属圆柱筒外表面的局部三维重构结果如图4所示。
Figure GDA0002681061380000113
以上所述,仅为本发明的较佳实施例而已,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,所有这些改进和变换都应属于本发明的保护范围。

Claims (6)

1.一种光栅式双目立体视觉三维测量系统,其特征在于:包括双目立体视觉传感器和左右结构光传感器、基于Hessian矩阵的光条中心点提取算法模块、基于近邻反投影光条中心点约束和极线约束的特征点匹配算法模块,双目立体视觉传感器和左右结构光传感器的三维重构结果融合模块;其中:
双目立体视觉传感器和左右结构光传感器,包括:左工业相机、右工业相机和一个多线激光器,三者按照左工业相机、多线激光器和右工业相机的顺序摆放,使得多线激光器投射的多根光条均位于左右工业相机的视场中间;其中,左右工业相机组成一个双目立体视觉传感器,而左右工业相机分别与多线激光器构成左右结构光传感器;多线激光器向物体表面投射多条光条,左右工业相机分别采集物体的左右图像;
基于Hessian矩阵的光条中心点提取算法模块:提取左右图像中光条中心点坐标,构建左、右特征点集;
基于近邻反投影光条中心线约束和极线约束的光条中心点匹配算法模块,遍历左特征点集,对于每一个光条中心点,计算对应右图像的近邻反投影光条中心线和极线,将近邻反投影光条中心线与极线的交点作为侯匹配点,若匹配成功,构成匹配点对,并标记右特征点集中相应的两个光条中心点;否则,则标记左特征点集中的该光条中心点;
双目立体视觉传感器和左右结构光传感器的三维重构结果融合模块:基于双目立体视觉测量模型重构匹配点对;采用左结构光视觉测量模型重构左特征点集中被标记的光条中心点;遍历右特征点集,将右特征点集中未被标记的光条中心点带入右结构光视觉测量模型中,获得对应的空间点的三维坐标;并将左右结构光视觉测量模型的重构结果转换到双目立体视觉传感器坐标系下。
2.根据权利要求1所述的一种光栅式双目立体视觉三维测量系统,其特征在于:所述双目立体视觉传感器和左右结构光传感器的实现如下:
(1)分别对双目立体视觉传感器和左右结构光传感器进行标定,获得左右工业相机的内参、双目立体视觉传感器的结构参数及左右结构光传感器的光平面方程,其中,左工业相机内参Kl,右工业相机内参Kr分别为:
Figure FDA0002681061370000011
其中,αl和αr分别为左工业相机、右工业相机u轴的有效焦距,βl和βr分别为左工业相机、右工业相机v轴的有效焦距,(u0l,v0l)和(u0r,v0r)分别为左工业相机、右工业相机的光学中心,γl和γr分别为左工业相机、右工业相机u轴与v轴不垂直因子;双目立体视觉传感器的结构参数,即旋转矩阵R,平移矢量T为:
Figure FDA0002681061370000021
T=[tx,ty,tz]T
其中,rs(s=1,...,9)为旋转矩阵R的第s个元素,tx,ty和tz分别为平移矢量T的三个分量;左结构传感器光平面方程:
aix+biy+ciz+di=0 i=0,1,...,m
其中,(ai,bi,ci)为第i个左结构传感器光平面的法向量,di为第i个左结构传感器光平面的截距,(x,y,z)为左结构传感器光平面上的三维点,m为多线激光器的投射的光条根数;
右结构传感器光平面方程:
a′ix′+b′iy′+c′iz′+d′i=0 i=0,1,...,m
其中,(a′i,b′i,c′i)为第i个右结构传感器光平面的法向量,d′i为第i个右结构传感器光平面的截距,(x′,y′,z′)为右结构传感器光平面上的三维点;
(2)多线激光器向物体表面投射多根光条,左右工业相机从不同角度同时采集物体左右图像,左右图像中均包含携带物体表面高度信息的光条。
3.根据权利要求1所述的一种光栅式双目立体视觉三维测量系统,其特征在于:所述基于Hessian矩阵的光条中心点提取算法模块中,通过求取图像Hessian矩阵的最大特征向量,获得左右图像中光条的中心点坐标,从而构建左右特征点集。
4.根据权利要求1所述的一种光栅式双目立体视觉三维测量系统,其特征在于:所述基于近邻反投影光条中心点约束和极线约束的光条中心点匹配算法模块具体实现过程如下:
(1)遍历左特征点集中的光条中心点(ulj,vlj),基于双目立体视觉测量模型,计算光条中心点(ulj,vlj)对应右图像上的极线Ej,具体计算过程如下:
Figure FDA0002681061370000022
其中,F为3×3的基本矩阵,
Figure FDA0002681061370000031
为光条中心点的齐次坐标形式;
(2)基于左结构光视觉测量模型,获得m个三维点坐标Pi=(xi,yi,zi)T,具体计算过程如下:
Figure FDA0002681061370000032
其中,
Figure FDA0002681061370000033
为光条中心点的齐次坐标形式,Kl为左工业相机内参,ρ为比例因子,(ai,bi,ci)为第i个左结构传感器光平面的法向量,di为第i个左结构传感器光平面的截距;
(3)将m个三维点坐标Pi=(xi,yi,zi)T反投影到右图像上,获得在右图像上的m个反投影点坐标(uri,vri);
(4)对于每个反投影点(uri,vri),搜索右特征点集中与该反投影点欧氏距离最小的两个光条中心点,获得对应的近邻反投影中心线li,同时标记右特征点集中的两个光条中心点;
(5)计算极线Ej与m根近邻反投影中心线li的交点(u′ri,v′ri),将得到的m个交点作为光条中心点(ulj,vlj)在右图像上的侯匹配点,构成m对侯匹配点对;
(6)将m个侯匹配点对代入双目立体视觉测量模型中,实现三维重构,获得m个对应的三维坐标P′i=(x′i,y′i,z′i)T
(7)将m个交点(u′ri,v′ri)代入右结构光视觉测量模型中,得到m个三维坐标Pi″=(x″i,y″i,z″i)T,具体计算过程如下:
Figure FDA0002681061370000034
其中,
Figure FDA0002681061370000035
为交点的齐次坐标形式,Kr为右工业相机内参,(a′i,b′i,c′i)为第i个右结构传感器光平面的法向量,d′i为第i个右结构传感器光平面的截距;
(8)将三维点Pi=(xi,yi,zi)T,Pi′=(x′i,y′i,z′i)T,Pi″=(x″i,y″i,z″i)T一一对应,构成m个三维点对{Pi,Pi′,Pi″},计算每个三维点对中两两三维点之间的欧氏距离之和,若最小的欧氏距离之和小于给定的阈值,则光条中心点(ulj,vlj)与对应的交点(u′ri,v′ri)为一对匹配点对;否则,则表明光条中心点(ulj,vlj)在右图像上无匹配点,并在左特征点集中标记该光条中心点(ulj,vlj)。
5.根据权利要求1所述的一种光栅式双目立体视觉三维测量系统,其特征在于:所述双目立体视觉传感器和左右结构光传感器的三维重构结果融合模块具体实现过程如下:
(1)基于双目立体视觉测量模型重构匹配点对所对应空间点的三维坐标Pb
(2)采用左结构光视觉测量模型重构左特征点集中所有被标记的光条中心点Pl
(3)将右特征点集中所有未被标记的光条中心点带入右结构光视觉测量模型中,获得对应的空间点的三维坐标Pr
(4)由于基于双目立体视觉测量模型的重构结果Pb与基于左结构光视觉测量模型的重构结果Pl均位于双目立体视觉传感器坐标系下,而基于右结构光视觉测量模型的重构结果Pr位于右工业相机坐标系下,需要将基于右结构光视觉测量模型的重构结果转到双目立体视觉传感器坐标系下,得到融合后的三维重构结果FinalCloud,其具体实现过程如下:
Figure FDA0002681061370000041
其中,R和T分别为双目立体视觉传感器的旋转矩阵和平移矢量,Pr′为位于双目立体视觉传感器坐标系下的右结构光视觉测量模型的重构结果。
6.一种光栅式双目立体视觉三维测量方法,其特征在于,步骤如下:
(1)多线激光器向物体表面投射多条光条,左右工业相机分别采集物体的左右图像;
(2)通过求取左右Hessian矩阵的最大特征向量,获得左右图像中光条的中心点坐标,构建左、右特征点集;
(3)为了获得准确唯一的匹配,提高三维重构精度,在现有极线约束的基础上增加近邻反投影光条中心点约束,将匹配过程从二维匹配优化为点到点匹配过程,遍历左特征点集中的光条中心点(ulj,vlj),基于双目立体视觉测量模型,计算光条中心点(ulj,vlj)对应右图像上的极线Ej,具体计算过程如下:
Figure FDA0002681061370000042
其中,F为3×3的基本矩阵,
Figure FDA0002681061370000043
为光条中心点的齐次坐标形式;
(4)基于左结构光视觉测量模型,获得m个三维点坐标Pi=(xi,yi,zi)T,具体计算过程如下:
Figure FDA0002681061370000051
其中,
Figure FDA0002681061370000052
为光条中心点的齐次坐标形式,Kl为左工业相机内参,ρ为比例因子,(ai,bi,ci)为第i个左结构传感器光平面的法向量,di为第i个左结构传感器光平面的截距;
(5)将m个三维点坐标Pi=(xi,yi,zi)T反投影到右图像上,获得在右图像上的m个反投影点坐标(uri,vri);
(6)对于每个反投影点(uri,vri),搜索右特征点集中与该反投影点欧氏距离最小的两个光条中心点,获得对应的近邻反投影中心线li,同时标记右特征点集中的两个光条中心点;
(7)计算极线Ej与m根近邻反投影中心线li的交点(u′ri,v′ri),将得到的m个交点作为光条中心点(ulj,vlj)在右图像上的侯匹配点,构成m对侯匹配点对;
(8)将m个侯匹配点对代入双目立体视觉测量模型中,实现三维重构,获得m个对应的三维坐标Pi′=(x′i,y′i,z′i)T
(9)将m个交点(u′ri,v′ri)代入右结构光视觉测量模型中,得到m个三维坐标Pi″=(x″i,y″i,z″i)T,具体计算过程如下:
Figure FDA0002681061370000053
其中,
Figure FDA0002681061370000054
为交点的齐次坐标形式,Kr为右工业相机内参,ρ为比例因子,(a′i,b′i,c′i)为第i个右结构传感器光平面的法向量,d′i为第i个右结构传感器光平面的截距;
(10)将三维点Pi=(xi,yi,zi)T,Pi′=(x′i,y′i,z′i)T,Pi″=(x″i,y″i,z″i)T一一对应,构成m个三维点对{Pi,Pi′,Pi″},计算每个三维点对中两两三维点之间的欧氏距离之和,若最小的欧氏距离之和小于给定的阈值,则光条中心点(ulj,vlj)与对应的交点(u′ri,v′ri)为一对匹配点对,将基于双目立体视觉测量模型重构的三维坐标Pi′=(x′i,y′i,z′i)T作为其对应空间点的三维坐标;否则,表明由于遮挡、反光、阴影干扰导致光条中心点(ulj,vlj)在右图像上无匹配点,将基于左结构光视觉测量模型重构的三维点Pi=(xi,yi,zi)T作为其对应空间点的三维坐标;
(11)当左特征点集中所有光条中心点遍历完之后,右特征点集中未被标记的光条中心点带入右结构光视觉测量模型中,获得对应的空间点的三维坐标M;
(12)在遮挡、反光、阴影干扰的情况下,由于双目立体视觉传感器视场的局限性,基于双目立体视觉测量模型的重构结果会存在缺失,因此需要将双目立体视觉传感器与左右结构光传感器的三维重构结果融合;基于双目立体视觉测量模型的重构结果P′与基于左结构光视觉测量模型的重构结果P均位于双目立体视觉传感器坐标系下,而基于右结构光视觉测量模型的重构结果M位于右工业相机坐标系下;需要将基于右结构光视觉测量模型的重构结果转到双目立体视觉传感器坐标系下,从而得到融合后的完整的三维重构结果FinalCloud,其具体实现过程如下:
Figure FDA0002681061370000061
其中,R和T分别为双目立体视觉传感器的旋转矩阵和平移矢量,M′为位于双目立体视觉传感器坐标系下的右结构光视觉测量模型的重构结果。
CN201910460042.6A 2019-05-30 2019-05-30 一种光栅式双目立体视觉三维测量系统及方法 Active CN110702025B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910460042.6A CN110702025B (zh) 2019-05-30 2019-05-30 一种光栅式双目立体视觉三维测量系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910460042.6A CN110702025B (zh) 2019-05-30 2019-05-30 一种光栅式双目立体视觉三维测量系统及方法

Publications (2)

Publication Number Publication Date
CN110702025A CN110702025A (zh) 2020-01-17
CN110702025B true CN110702025B (zh) 2021-03-19

Family

ID=69193046

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910460042.6A Active CN110702025B (zh) 2019-05-30 2019-05-30 一种光栅式双目立体视觉三维测量系统及方法

Country Status (1)

Country Link
CN (1) CN110702025B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111854642B (zh) * 2020-07-23 2021-08-10 浙江汉振智能技术有限公司 基于随机点阵的多线激光三维成像方法及系统
CN112304245A (zh) * 2020-11-22 2021-02-02 西北工业大学 一种框肋零件外形偏差测量方法
CN113008158B (zh) * 2021-03-25 2023-02-24 烟台大学 多线激光轮胎花纹深度测量方法
CN113983933B (zh) * 2021-11-11 2022-04-19 易思维(杭州)科技有限公司 一种多线激光传感器的标定方法
CN114098985A (zh) * 2021-11-29 2022-03-01 北京柏惠维康科技有限公司 患者与患者医学影像的空间匹配方法、装置、设备及介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005509879A (ja) * 2001-11-21 2005-04-14 マプヴィジョン オーユー エルテーデー 三次元計測における対応ポイントを決定する方法
CN102062588A (zh) * 2009-11-11 2011-05-18 中国科学院沈阳自动化研究所 一种计算机双目视觉义齿扫描装置及其三维重建方法
CN105091782A (zh) * 2015-05-29 2015-11-25 南京邮电大学 一种基于双目视觉的多线激光器光平面标定方法
CN105157602A (zh) * 2015-07-13 2015-12-16 西北农林科技大学 基于机器视觉的远程三维扫描系统及方法
CN106802138A (zh) * 2017-02-24 2017-06-06 杭州先临三维科技股份有限公司 一种三维扫描系统及其扫描方法
CN107167118A (zh) * 2017-06-05 2017-09-15 南京航空航天大学 一种基于非编码平行多线的稳定实时激光测量方法
CN107907048A (zh) * 2017-06-30 2018-04-13 长沙湘计海盾科技有限公司 一种基于线结构光扫描的双目立体视觉三维测量方法
CN108151671A (zh) * 2016-12-05 2018-06-12 杭州先临三维科技股份有限公司 一种三维数字成像传感器、三维扫描系统及其扫描方法
CN108267097A (zh) * 2017-07-17 2018-07-10 杭州先临三维科技股份有限公司 基于双目三维扫描系统的三维重构方法和装置
CN108375346A (zh) * 2018-02-11 2018-08-07 西安知象光电科技有限公司 一种基于微振镜分光式投影装置的多线激光三维扫描方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005509879A (ja) * 2001-11-21 2005-04-14 マプヴィジョン オーユー エルテーデー 三次元計測における対応ポイントを決定する方法
CN102062588A (zh) * 2009-11-11 2011-05-18 中国科学院沈阳自动化研究所 一种计算机双目视觉义齿扫描装置及其三维重建方法
CN105091782A (zh) * 2015-05-29 2015-11-25 南京邮电大学 一种基于双目视觉的多线激光器光平面标定方法
CN105157602A (zh) * 2015-07-13 2015-12-16 西北农林科技大学 基于机器视觉的远程三维扫描系统及方法
CN108151671A (zh) * 2016-12-05 2018-06-12 杭州先临三维科技股份有限公司 一种三维数字成像传感器、三维扫描系统及其扫描方法
CN106802138A (zh) * 2017-02-24 2017-06-06 杭州先临三维科技股份有限公司 一种三维扫描系统及其扫描方法
CN107167118A (zh) * 2017-06-05 2017-09-15 南京航空航天大学 一种基于非编码平行多线的稳定实时激光测量方法
CN107907048A (zh) * 2017-06-30 2018-04-13 长沙湘计海盾科技有限公司 一种基于线结构光扫描的双目立体视觉三维测量方法
CN108267097A (zh) * 2017-07-17 2018-07-10 杭州先临三维科技股份有限公司 基于双目三维扫描系统的三维重构方法和装置
CN108375346A (zh) * 2018-02-11 2018-08-07 西安知象光电科技有限公司 一种基于微振镜分光式投影装置的多线激光三维扫描方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
光栅式双目立体视觉传感器的立体匹配方法;魏振忠等;《红外与激光工程》;20100425;第39卷(第2期);第330-334页 *
基于同心圆合成图像匹配的双目视觉标定;侯俊捷等;《光学学报》;20120310;第32卷(第3期);第0315003-1-0315003-5页 *

Also Published As

Publication number Publication date
CN110702025A (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CN110702025B (zh) 一种光栅式双目立体视觉三维测量系统及方法
CN110288642B (zh) 基于相机阵列的三维物体快速重建方法
CN110230998B (zh) 基于线激光和双目相机的快速精密三维测量方法和装置
CN100430690C (zh) 利用单数码相机自由拍摄进行物体三维测量的方法
CN110009667B (zh) 基于罗德里格斯变换的多视点云全局配准方法
CN115345822A (zh) 一种面向航空复杂零件的面结构光自动化三维检测方法
CN105574812B (zh) 多角度三维数据配准方法及装置
CN114998499A (zh) 一种基于线激光振镜扫描的双目三维重建方法及系统
CN104835158A (zh) 基于格雷码结构光与极线约束的三维点云获取方法
CN105241397A (zh) 基于结构光的实时测量拼接方法及其设备
CN103366360A (zh) 信息处理设备和信息处理方法
CN104408762A (zh) 利用单目和二维平台获取物体图像信息及三维模型的方法
KR20160121509A (ko) 2개의 카메라로부터의 곡선의 세트의 구조형 광 매칭
Fernandez et al. Planar-based camera-projector calibration
He et al. 3D Surface reconstruction of transparent objects using laser scanning with LTFtF method
CN110940295A (zh) 基于激光散斑极限约束投影的高反射物体测量方法及系统
CA3233222A1 (en) Method, apparatus and device for photogrammetry, and storage medium
CN112184793B (zh) 深度数据的处理方法、装置及可读存储介质
CN116295113A (zh) 一种融合条纹投影的偏振三维成像方法
CN116242277A (zh) 基于全场三维视觉的电源机箱结构件尺寸自动测量方法
CN103884294A (zh) 一种宽视场的红外光三维形貌测量的方法及其装置
CN118334239A (zh) 基于条纹投影测量的管道三维重建方法和设备
Yamauchi et al. Calibration of a structured light system by observing planar object from unknown viewpoints
Grifoni et al. 3D multi-modal point clouds data fusion for metrological analysis and restoration assessment of a panel painting
CN117450955A (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