CN107133989B - 一种三维扫描系统参数标定方法 - Google Patents
一种三维扫描系统参数标定方法 Download PDFInfo
- Publication number
- CN107133989B CN107133989B CN201710445013.3A CN201710445013A CN107133989B CN 107133989 B CN107133989 B CN 107133989B CN 201710445013 A CN201710445013 A CN 201710445013A CN 107133989 B CN107133989 B CN 107133989B
- Authority
- CN
- China
- Prior art keywords
- camera
- calibration
- dimensional scanning
- scanning system
- parameters
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
- G06T7/85—Stereo camera calibration
Abstract
本发明公开了一种三维扫描系统参数标定方法,包括:从由三维扫描系统拍摄获得的标定物不同角度的多组标定物图像中,获取标定点数据;基于获得的标定点数据对三维扫描系统的各个相机的参数分别进行初步标定,得到初步标定结果;基于初步标定结果,获得三维扫描系统中第一相机与第二相机之间进行位置坐标变换的变换关系;基于第一相机和第二相机的初步标定结果以及变换关系,对第一相机、第二相机的参数进行迭代运算,获得三维扫描系统的最终参数结果。本发明三维扫描系统参数标定方法,在对各单个相机进行初步标定后,从系统角度出发对标定参数进行全局优化,使获得的参数标定结果针对于系统精度更高。
Description
技术领域
本发明涉及计算机视觉技术领域,特别是涉及一种三维扫描系统参数标定方法。
背景技术
长久以来,对真实物体几何尺寸的精确测量主要采用依赖于游标卡尺、千分尺、角度尺等测量工具的接触式人工测量方式,而这种测量方式无法测量形状不规则的物体表面,并且对于文物古迹等要求非接触的测量对象也存在技术瓶颈。
随着计算机视觉技术的发展以及新技术新需求的驱动,三维扫描系统成为非接触式测量的重要手段,被广泛应用于工业设计领域的辅助制造与检验、医疗领域的手术定位与康复、游戏娱乐领域的实景建模与仿真以及考古领域的遗址保护与复原等。
三维扫描系统工作过程主要包括数据采集和三维模型重构两部分,三维扫描仪由结构光光源和双目相机组成,在数据采集过程中,当扫描对象将结构光光线反射进互成一定夹角的双目相机中时,便可得到扫描对象的一组立体图像。随着三维扫描仪位置的不断变化,立体图像的数量不断增多,这些立体图像分别记录了扫描对象不同位置、不同角度的图像信息,进一步利用依次获得的多组立体图像,依据三角测量原理,构建出三维模型。
三维扫描系统在进行三维模型重构时,涉及的首要关键问题就是三维扫描系统的标定问题,即如何获得三维扫描系统中相机的内部参数和外部参数,以实现基于相机内部、外部参数的三角测量原理,进而完成三维模型重构。
现有技术中,采用的众多标定方法只针对三维扫描系统中的单个相机自身进行,比如利用张正友标定法实现相机内、外参数的标定,但获得的标定参数相对于系统整体来说精度不高。
发明内容
本发明的目的是提供一种三维扫描系统参数标定方法,在对各单个相机进行初步标定后,从系统角度出发对标定参数进行全局优化,使获得的参数标定结果针对于系统精度更高。
为实现上述目的,本发明提供如下技术方案:
一种三维扫描系统参数标定方法,包括:
从由三维扫描系统拍摄获得的标定物不同角度的多组标定物图像中,获取标定点数据;
基于获得的所述标定点数据对所述三维扫描系统的各个相机的参数分别进行初步标定,得到初步标定结果;
基于所述初步标定结果,获得所述三维扫描系统中第一相机与第二相机之间进行位置坐标变换的变换关系;
基于所述第一相机和所述第二相机的初步标定结果以及所述变换关系,对所述第一相机、所述第二相机的参数进行迭代运算,获得所述三维扫描系统的最终参数结果。
可选地,所述获取标定点数据包括:
以标定物建立世界坐标系,以所述标定物上设置的标定点的空间坐标作为所述标定点的世界坐标;
从所述标定物图像中识别出所述标定点,并提取所述标定点的像素坐标;
以所述标定点的世界坐标和像素坐标的组合,作为所述标定点的标定点数据。
可选地,基于获得的所述标定点数据对所述三维扫描系统单个相机的参数进行初步标定,得到初步标定结果包括:
迭代初始化,输入相应相机的标定点数据、设定的所述相机的初始内部参数以及根据所述标定点数据和所述初始内部参数计算得到的所述相机的外部参数;
重投影运算,对应所述相机拍摄的各标定物图像,将标定点的空间坐标重投影到相应的像素坐标系,得到各标定物图像的重投影像素坐标矩阵,并计算得到所述重投影像素坐标矩阵对内部参数各项分别求导获得的雅克比矩阵Jfm、Jcm、Jkm和对外部参数各项分别求导获得的雅克比矩阵Jom、Jtm,以及将输入的像素坐标与所述重投影像素坐标做差运算得到的像素偏差矩阵,记为单列矩阵Dm,其中m表示第m幅标定物图像,m∈[1,p];
建立海森矩阵,构建矩阵Am=[Jfm Jcm Jkm]T和Bm=[Jom Jtm]T,并构建海森矩阵H以及用于迭代运算的雅克比偏差单列矩阵V,如下:
进行迭代运算获得初步标定结果,具体包括:计算所述相机的标定参数偏差矩阵,计算表达式为:P=[1-(1-α)(i+1)]H-1V,其中,i表示迭代次数,α表示平滑系数,α大于0小于1;β表示迭代精度系数,迭代精度系数β的表达式为:
其中,下标new表示计算得到的新参数,标定参数偏差矩阵P中前8项为内部参数偏差,按顺序依次包括等效焦距fu、fv,镜头成像中心像素坐标cu、cv,畸变系数k1、k2、k3、k4,之后诸项为多组外部参数;
根据更新得到的内部参数以及初始化时输入的各标定物图像的标定点数据,更新所述相机对应各标定物图像的外部参数,进行迭代运算直到满足预设条件,则结束迭代运算过程,获得所述相机的初步标定结果。
可选地,所述基于所述初步标定结果,获得所述三维扫描系统中第一相机与第二相机之间进行位置坐标变换的变换关系包括:
基于所述初步标定结果,获得所述三维扫描系统中第一相机与第二相机之间进行坐标变换的旋转向量O和平移向量T;
具体包括:根据以下计算式计算得到p组结果,对p组结果求平均值得到旋转矩阵R和平移向量T,根据旋转矩阵R得到旋转向量O;
其中,R1和t1表示第一相机的外部参数,R2和t2表示第二相机的外部参数。
可选地,所述基于所述第一相机和所述第二相机的初步标定结果以及所述变换关系,对所述第一相机、所述第二相机的参数进行迭代运算,获得所述三维扫描系统的最终参数结果,包括:
标定点匹配,对所述第一相机和所述第二相机同一时刻获取的同一幅标定物图像中获取的标定点进行匹配;
迭代初始化,输入经初步标定获得的所述第一相机、所述第二相机的内部参数,所述第一相机的外部参数,旋转向量和平移向量,以及匹配后的标定点的世界坐标和像素坐标;
重投影运算,将第一相机与第二相机的匹配标定点,从世界坐标系经坐标变换和增加畸变分别投影到第一相机和第二相机的像素坐标系,重投影运算中所用的所述第二相机外部参数由所述第一相机的外部参数经所述旋转向量和所述平移向量计算得到;重投影运算得到各组标定物图像对应所述第一相机、所述第二相机的重投影像素坐标矩阵以及雅克比矩阵,所述雅克比矩阵包括所述第一相机的重投影像素坐标矩阵对所述第一相机的内部参数、外部参数求导得到的雅克比矩阵Jflm、Jclm、Jklm、Jolm、Jtlm,所述第二相机的重投影像素坐标矩阵对所述第二相机的内部参数求导得到的雅克比矩阵Jfrm、Jcrm、Jkrm,所述第二相机的重投影像素坐标矩阵对所述第一相机外部参数求导得到的雅克比矩阵Jorm、Jtrm,以及所述第二相机的重投影像素坐标矩阵对所述旋转向量、平移向量求导得到的雅克比矩阵Jodm、Jtdm,将所述第一相机、所述第二相机匹配后的标定点像素坐标与所述第一相机、所述第二相机重投影像素坐标做差运算,并写成单列矩阵Em,其中m表示第m幅标定物图像,m∈[1,p];
建立标定参数全局优化偏差运算矩阵J以及像素偏差单列矩阵W,如下:
进行标定参数全局优化运算,计算标定参数全局优化偏差矩阵,计算表达式为:Q=(JTJ)-1JTW,其中,Q表示标定参数偏差矩阵;精度系数γ的计算表达式为:
其中,O表示旋转向量,T表示平移向量,下标new表示计算得到的新的数据;
根据初始化时输入的所述第一相机、所述第二相机的内部参数、所述第一相机的外部参数、所述旋转向量和平移向量,以及所述标定参数全局优化偏差矩阵Q,更新所述第一相机、所述第二相机的内部参数、所述第一相机的外部参数、所述旋转向量和所述平移向量,进行迭代运算直到满足预设条件,则结束迭代运算过程,获得所述三维扫描系统的最终参数。
可选地,还包括:获取所述三维扫描系统的结构光平面的平面方程,对所述三维扫描系统的结构光进行标定。
可选地,具体包括:
由三维扫描系统拍摄多幅经不同位置平面截取的结构光光线图像;
以获得的各相机的标定结果,对结构光光线图像进行畸变校正和极线校正,提取图像中结构光光线成像点的像素坐标;
进行空间位置重构,计算所述第一相机拍摄的图像中结构光光线成像点对应的实际物点的空间坐标(X1,Y1,Z1),计算表达式为:
其中:RT1为旋转平移矩阵,u1、v1表示结构光光线在所述第一相机中成像点像素坐标,u2、v2表示结构光光线在所述第二相机中成像点像素坐标,fu1、fv1表示所述第一相机的等效焦距,cu1、cv1表示所述第一相机的镜头成像中心像素坐标,fu2表示所述第二相机的u轴等效焦距,cu2表示所述第二相机图像的镜头成像中心u轴像素坐标;
利用在所述第一相机坐标系下重构出的成像点空间坐标(X1,Y1,Z1)拟合平面方程,所述平面方程的形式为Z=a+bX+cY,其中a,b,c为平面方程系数。
可选地,所述标定物包括平面标定图案,标定图案包括圆点和围绕圆点的带有编码信息的间断圆环,图案中圆点的圆心坐标作为标定点的世界坐标。
可选地,所述从所述标定物图像中识别出所述标定点,并提取所述标定点的像素坐标包括:
对标定物图像进行边缘检测,得到二值化的轮廓边缘图像;
对轮廓边缘图像进行椭圆拟合,获得轮廓椭圆的中心像素坐标;
对椭圆进行圆度判断和倾斜度判断,初步提取出有效的标定圆点像素坐标;
对标定圆点周围的编码圆环进行解码,得到标定圆点编号,确定标定圆点在标定板的位置,以与相应的世界坐标匹配;
对初步提取的标定圆点进行随机抽样,进行一致性判断,去除提取出的错误的标定圆点,得到最终有效的标定圆点及其像素坐标。
由上述技术方案可知,本发明所提供的三维扫描系统参数标定方法,从由三维扫描系统拍摄获得的标定物不同角度的多组标定物图像中,获取标定点数据,基于获得的所述标定点数据对所述三维扫描系统的各个相机的参数分别进行初步标定,得到初步标定结果,然后基于所述初步标定结果,获得所述三维扫描系统中第一相机与第二相机之间进行位置坐标变换的变换关系,基于第一相机和第二相机的初步标定结果以及所述变换关系,对第一相机、第二相机的参数进行迭代运算,获得三维扫描系统的最终参数结果。
本发明三维扫描系统参数标定方法,在对各单个相机进行初步标定后,从系统角度出发对标定参数进行全局优化,使获得的参数标定结果针对于系统精度更高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种三维扫描系统参数标定方法的流程图;
图2为本发明实施例中获取标定点数据采用的一种标定物示意图;
图3为本发明实施例中获取标定点数据时建立的坐标系;
图4为本发明又一实施例提供的一种三维扫描系统参数标定方法的流程图。
具体实施方式
为了使本技术领域的人员更好地理解本发明中的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
本发明实施例提供一种三维扫描系统参数标定方法,包括:
从由三维扫描系统拍摄获得的标定物不同角度的多组标定物图像中,获取标定点数据;
基于获得的所述标定点数据对所述三维扫描系统的各个相机的参数分别进行初步标定,得到初步标定结果;
基于所述初步标定结果,获得所述三维扫描系统中第一相机与第二相机之间进行位置坐标变换的变换关系;
基于所述第一相机和所述第二相机的初步标定结果以及所述变换关系,对所述第一相机、所述第二相机的参数进行迭代运算,获得所述三维扫描系统的最终参数结果。
可以看出,本实施例三维扫描系统参数标定方法,从由三维扫描系统拍摄获得的标定物不同角度的多组标定物图像中,获取标定点数据,基于获得的所述标定点数据对所述三维扫描系统的各个相机的参数分别进行初步标定,得到初步标定结果,然后基于所述初步标定结果,获得所述三维扫描系统中第一相机与第二相机之间进行位置坐标变换的变换关系,基于第一相机和第二相机的初步标定结果以及所述变换关系,对第一相机、第二相机的参数进行迭代运算,获得三维扫描系统的最终参数结果。
本实施例三维扫描系统参数标定方法,在对各单个相机进行初步标定后,从系统角度出发对标定参数进行全局优化,使获得的参数标定结果针对于系统精度更高。
下面结合具体实施方式对本实施例三维扫描系统参数标定方法进行详细说明。
请参考图1,本实施例提供的三维扫描系统参数标定方法,包括步骤:
S10:从由三维扫描系统拍摄获得的标定物不同角度的多组标定物图像中,获取标定点数据。
由三维扫描系统以不同角度对标定物拍摄图像,三维扫描系统的各个相机均拍摄图像,获得不同角度的多组标定物图像。
从标定物不同角度的多组标定物图像中分别获取标定点数据,其中,在标定物上设置有空间坐标已知的标定点,标定点数据包括标定点在图像中的像素坐标和在空间中的世界坐标(即空间坐标)。
具体的,获取标定点数据包括以下步骤:
S100:以标定物建立世界坐标系,以所述标定物上设置的标定点的空间坐标作为所述标定点的世界坐标。
S101:从所述标定物图像中识别出所述标定点,提取所述标定点的像素坐标。
S102:以所述标定点的世界坐标和像素坐标的组合,作为所述标定点的标定点数据。
例如,请参考图2和图3,图2为本实施例中获取标定点数据采用的一种标定物示意图,图3为本实施例中获取标定点数据时建立的坐标系,在图3中Ow-XwYwZw表示世界坐标系,O-XcYcZc表示相机坐标系,OI-xy表示图像坐标系,Oo-uv表示像素坐标系。若标定物选取如图2所示的平面标定图案,标定图案包括圆点和围绕圆点的带有编码信息的间断圆环,则标定点数据的获取在于提取标定圆点的圆心坐标,其世界坐标在Zw方向始终为0,Xw和Yw方向的圆心坐标已知。
其像素坐标可以通过以下步骤获取:首先,对标定物图像进行边缘检测,得到二值化的轮廓边缘图像;其次,对轮廓边缘图像进行椭圆拟合,获得轮廓椭圆的中心像素坐标;然后,对椭圆进行圆度判断和倾斜度判断,初步提取出有效的标定圆点像素坐标;接下来,对标定圆点周围的编码圆环进行解码,得到标定圆点编号,从而确定标定圆点在标定板的位置,与相应的世界坐标匹配;最后,进一步对初步提取的标定圆点进行随机抽样,进行一致性判断,去除提取出的错误的标定圆点,得到最终有效的标定圆点及其像素坐标。
S11:基于获得的所述标定点数据对所述三维扫描系统的各个相机的参数分别进行初步标定,得到初步标定结果。
对三维扫描系统中单个相机进行参数标定包括计算得到单个相机的内部参数,以及单个相机拍摄的每幅标定物图像对应的外部参数,其中,内部参数包括等效焦距fu、fv,镜头成像中心像素坐标cu、cv,以及畸变系数k1、k2、k3、k4。外部参数包括单个相机以不同角度拍摄标定物时,获得的每幅标定物图像对应的世界坐标系到相机坐标系的旋转向量o和平移向量t。
本实施例方法中,优选的,对单个相机标定的计算过程采用对相机相应的内部参数、外部参数迭代运算逼近真值的方法。
具体的,对单个相机的内部参数、外部参数迭代运算的具体过程为:
S110:迭代初始化,输入相应相机的标定点数据、设定的所述相机的初始内部参数以及根据所述标定点数据和所述初始内部参数计算得到的所述相机的外部参数。
输入的标定点数据包括标定点的世界坐标和像素坐标。输入的相机初始内部参数由人为设定。
S111:重投影运算,对应所述相机拍摄的各标定物图像,将标定点的空间坐标重投影到相应的像素坐标系,得到各标定物图像的重投影像素坐标矩阵,并计算得到所述重投影像素坐标矩阵对内部参数各项分别求导获得的雅克比矩阵Jfm、Jcm、Jkm和对外部参数各项分别求导获得的雅克比矩阵Jom、Jtm,以及将输入的像素坐标与所述重投影像素坐标做差运算得到的像素偏差矩阵,记为单列矩阵Dm,其中m表示第m幅标定物图像,m∈[1,p]。
其中,将标定物图像中标定点的空间坐标,由世界坐标系经内部参数、外部参数作用,重投影到像素坐标系,得到各标定物图像的重投影像素坐标矩阵。
S112:建立海森矩阵,构建矩阵Am=[Jfm Jcm Jkm]T和Bm=[Jom Jtm]T,并构建海森矩阵H以及用于迭代运算的雅克比偏差单列矩阵V:
所述海森矩阵H为迭代运算中用以计算标定参数偏差的运算矩阵,对于获取的p幅标定物图像的标定点数据,经步骤S111获得了p组雅克比矩阵和像素偏差矩阵,基于此来构建海森矩阵H。
S113:进行迭代运算获得初步标定结果,具体包括:计算所述相机的标定参数偏差矩阵,计算表达式为:P=[1-(1-α)(i+1)]H-1V,其中,i表示迭代次数,α表示平滑系数,α大于0小于1;β表示迭代精度系数,迭代精度系数β的表达式为:
其中,下标new表示计算得到的新参数,标定参数偏差矩阵P中前8项为内部参数偏差,按顺序依次包括等效焦距fu、fv,镜头成像中心像素坐标cu、cv,畸变像素k1、k2、k3、k4,之后诸项为多组外部参数偏差;
根据更新得到的内部参数以及初始化时输入的各标定物图像的标定点数据,更新所述相机对应各标定物图像的外部参数,进行迭代运算直到迭代精度系数达到预设值,或者迭代次数达到阈值极限,则结束迭代运算过程,获得所述相机的初步标定结果。
这样完成对三维扫描系统中单个相机的初步标定。
S12:基于所述初步标定结果,获得所述三维扫描系统中第一相机与第二相机之间进行位置坐标变换的变换关系。
第一相机与第二相机之间进行坐标变换是指第一相机与第二相机的坐标系之间的变换,包括进行坐标变换的旋转向量O和平移向量T。
具体过程包括:根据以下计算式计算得到p组结果,对p组结果求平均值得到旋转矩阵R和平移向量T,根据旋转矩阵R得到旋转向量O;
其中,R1和t1表示第一相机的外部参数,R2和t2表示第二相机的外部参数。
根据上步骤计算得到的第一相机、第二相机的p组标定物图像对应的外部参数进行计算,根据上述计算式计算得到p组结果。其中,旋转矩阵R经罗德里格斯变换得到相应旋转向量O,
S13:基于所述第一相机和所述第二相机的初步标定结果以及所述变换关系,对所述第一相机、所述第二相机的参数进行迭代运算,获得所述三维扫描系统的最终参数结果。
本方法中,以第一相机的内部参数、外部参数为基准,通过计算获得的第一相机与第二相机之间的旋转向量和平移向量,计算出第二相机的外部参数,使进行迭代运算时第二相机关联于第一相机,这样通过迭代运算实现标定参数全局优化。
具体包括以下过程:
S130:标定点匹配,对所述第一相机和所述第二相机同一时刻获取的同一幅标定物图像中获取的标定点进行匹配。
这样保证第一相机、第二相机从各组标定物图像中提取到相同的标定点。
S131:迭代初始化,输入经初步标定获得的所述第一相机、所述第二相机的内部参数,所述第一相机的外部参数,所述旋转向量和平移向量,以及匹配后的标定点的世界坐标和像素坐标。
S132:重投影运算,将第一相机与第二相机的匹配标定点,从世界坐标系经坐标变换和增加畸变分别投影到第一相机和第二相机的像素坐标系,重投影运算中所用的第二相机外部参数由所述第一相机的外部参数经所述旋转向量和所述平移向量计算得到。
重投影运算得到各组标定物图像对应所述第一相机、所述第二相机的重投影像素坐标矩阵以及雅克比矩阵,所述雅克比矩阵包括所述第一相机的重投影像素坐标矩阵对所述第一相机的内部参数、外部参数求导得到的雅克比矩阵Jflm、Jclm、Jklm、Jolm、Jtlm,所述第二相机的重投影像素坐标矩阵对所述第二相机的内部参数求导得到的雅克比矩阵Jfrm、Jcrm、Jkrm,所述第二相机的重投影像素坐标矩阵对所述第一相机外部参数求导得到的雅克比矩阵Jorm、Jtrm,以及所述第二相机的重投影像素坐标矩阵对所述旋转向量、平移向量求导得到的雅克比矩阵Jodm、Jtdm,将所述第一相机、所述第二相机匹配后的标定点像素坐标与所述第一相机、所述第二相机重投影像素坐标做差运算,并写成单列矩阵Em,其中m表示第m幅标定物图像,m∈[1,p]。
S133:建立标定参数全局优化偏差运算矩阵J以及像素偏差单列矩阵W,如下:
S134:进行标定参数全局优化运算,计算标定参数全局优化偏差矩阵,计算表达式为:Q=(JTJ)-1JTW;γ表示迭代精度系数,迭代精度系数γ的表达式为:
其中,O表示旋转向量,T表示平移向量,下标new表示计算得到的新数据;
根据初始化时输入的所述第一相机、所述第二相机的内部参数、所述第一相机的外部参数、所述旋转向量和平移向量,以及所述标定参数全局优化偏差矩阵Q,更新所述第一相机、所述第二相机的内部参数、所述第一相机的外部参数、所述旋转向量和平移向量,进行迭代运算直到满足预设条件,则结束迭代运算过程,获得所述三维扫描系统的最优参数。
这样,获得全局优化后高精度的第一相机、第二相机的内部、外部参数。
其中,进行迭代运算直到迭代精度系数达到预设值,或者迭代次数达到阈值,则结束迭代运算过程。
进一步优选的,在上述实施例内容的基础上,本实施例三维扫描系统参数标定方法,还包括步骤S14:获取所述三维扫描系统的结构光平面的平面方程,对所述三维扫描系统的结构光进行标定。
具体的,对三维扫描系统结构光进行标定的过程如下:
S140:由三维扫描系统拍摄多幅经不同位置平面截取的结构光光线图像。
S141:以获得的各相机的标定结果,对结构光光线图像进行畸变校正和极线校正,提取图像中结构光光线成像点的像素坐标。
S142:进行空间位置重构,计算所述第一相机拍摄的图像中结构光光线成像点对应的实际物点的空间坐标(X1,Y1,Z1),计算表达式为:
其中:RT1为旋转平移矩阵,u1、v1表示结构光光线在所述第一相机中成像点像素坐标,u2、v2表示结构光光线在所述第二相机中成像点像素坐标,fu1、fv1表示所述第一相机的等效焦距,cu1、cv1表示所述第一相机的镜头成像中心像素坐标,fu2表示所述第二相机的u轴等效焦距,cu2表示所述第二相机图像的镜头成像中心u轴像素坐标。
S143:利用在所述第一相机坐标系下重构出的成像点空间坐标(X1,Y1,Z1)拟合平面方程,所述平面方程的形式为Z=a+bX+cY,其中a,b,c为平面方程系数。
将重构出的在第一相机坐标系下结构光的空间坐标(X1,Y1,Z1)全部带入平面方程,构成结构光平面的非齐次线性方程组,方程组中方程个数远大于未知数个数,故对所述非齐次线性方程组求近似解,则可得到平面方程系数,即拟合出结构光平面的平面方程。
本实施例三维扫描系统参数标定方法,在对各单个相机进行初步标定后,从系统角度出发对标定参数进行全局优化,使获得的参数标定结果针对于系统精度更高。并进一步对三维扫描系统中的结构光进行标定,实现三维扫描系统的一次性全参数高精度标定。
本发明三维扫描系统参数标定方法,一次性获取多幅不同角度拍摄的标定物图像,以及多幅经不同位置平面截取的结构光光线图像,可完成对三维扫描系统中的全部关键参数的整体标定和优化工作。与传统方法相比,本标定方法针对于三维扫描系统,将各单个相机作为统一系统上的整体考虑,在获取相机之间的旋转平移向量后,对相机的标定参数进行全局优化,使相机参数针对于三维扫描系统具有更高精度。同时,本发明中还进行了对三维扫描系统中结构光的标定,使得本发明一次性的完成了三维扫描系统的全参数标定工作,具有高效、精准的优点。
以上对本发明所提供的一种三维扫描系统参数标定方法进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。
Claims (8)
1.一种三维扫描系统参数标定方法,其特征在于,包括:
从由三维扫描系统拍摄获得的标定物不同角度的多组标定物图像中,获取标定点数据;
基于获得的所述标定点数据对所述三维扫描系统的各个相机的参数分别进行初步标定,得到初步标定结果;
基于所述初步标定结果,获得所述三维扫描系统中第一相机与第二相机之间进行位置坐标变换的变换关系;
基于所述第一相机和所述第二相机的初步标定结果以及所述变换关系,对所述第一相机、所述第二相机的参数进行迭代运算,获得所述三维扫描系统的最终参数结果;
所述基于所述第一相机和所述第二相机的初步标定结果以及所述变换关系,对所述第一相机、所述第二相机的参数进行迭代运算,获得所述三维扫描系统的最终参数结果,包括:
标定点匹配,对所述第一相机和所述第二相机同一时刻获取的同一幅标定物图像中获取的标定点进行匹配;
迭代初始化,输入经初步标定获得的所述第一相机、所述第二相机的内部参数,所述第一相机的外部参数,旋转向量和平移向量,以及匹配后的标定点的世界坐标和像素坐标;
重投影运算,将第一相机与第二相机的匹配标定点,从世界坐标系经坐标变换和增加畸变分别投影到第一相机和第二相机的像素坐标系,重投影运算中所用的所述第二相机的外部参数由所述第一相机的外部参数经所述旋转向量和所述平移向量计算得到;重投影运算得到各组标定物图像对应所述第一相机、所述第二相机的重投影像素坐标矩阵以及雅克比矩阵,所述雅克比矩阵包括所述第一相机的重投影像素坐标矩阵对所述第一相机的内部参数、外部参数求导得到的雅克比矩阵Jflm、Jclm、Jklm、Jolm、Jtlm,所述第二相机的重投影像素坐标矩阵对所述第二相机的内部参数求导得到的雅克比矩阵Jfrm、Jcrm、Jkrm,所述第二相机的重投影像素坐标矩阵对所述第一相机外部参数求导得到的雅克比矩阵Jorm、Jtrm,以及所述第二相机的重投影像素坐标矩阵对所述旋转向量、平移向量求导得到的雅克比矩阵Jodm、Jtdm,将所述第一相机、所述第二相机匹配后的标定点像素坐标与所述第一相机、所述第二相机重投影像素坐标做差运算,并写成单列矩阵Em,其中m表示第m幅标定物图像,m∈[1,p];
建立标定参数全局优化偏差运算矩阵J以及像素偏差单列矩阵W,如下:
进行标定参数全局优化运算,计算标定参数全局优化偏差矩阵,计算表达式为:Q=(JTJ)-1JTW,其中,Q表示标定参数偏差矩阵;精度系数γ的计算表达式为:
其中,O表示旋转向量,T表示平移向量,下标new表示计算得到的新的数据;
根据初始化时输入的所述第一相机、所述第二相机的内部参数、所述第一相机的外部参数、所述旋转向量和所述平移向量以及所述标定参数全局优化偏差矩阵Q,更新所述第一相机、所述第二相机的内部参数、所述第一相机的外部参数,以及所述旋转向量和所述平移向量,进行迭代运算直到满足预设条件,则结束迭代运算过程,获得所述三维扫描系统的最终参数。
2.根据权利要求1所述的三维扫描系统参数标定方法,其特征在于,所述获取标定点数据包括:
以标定物建立世界坐标系,以所述标定物上设置的标定点的空间坐标作为所述标定点的世界坐标;
从所述标定物图像中识别出所述标定点,并提取所述标定点的像素坐标;
以所述标定点的世界坐标和像素坐标的组合,作为所述标定点的标定点数据。
3.根据权利要求1所述的三维扫描系统参数标定方法,其特征在于,基于获得的所述标定点数据对所述三维扫描系统单个相机的参数进行初步标定,得到初步标定结果包括:
迭代初始化,输入相应相机的标定点数据、设定的所述相机的初始内部参数以及根据所述标定点数据和所述初始内部参数计算得到的所述相机的外部参数;
重投影运算,对应所述相机拍摄的各标定物图像,将标定点的空间坐标重投影到相应的像素坐标系,得到各标定物图像的重投影像素坐标矩阵,并计算得到所述重投影像素坐标矩阵对内部参数各项分别求导获得的雅克比矩阵Jfm、Jcm、Jkm和对外部参数各项分别求导获得的雅克比矩阵Jom、Jtm,以及将输入的像素坐标与所述重投影像素坐标做差运算得到的像素偏差矩阵,记为单列矩阵Dm,其中m表示第m幅标定物图像,m∈[1,p];
构建矩阵Am=[Jfm Jcm Jkm]T和Bm=[Jom Jtm]T,并构建海森矩阵H以及用于迭代运算的雅克比偏差单列矩阵V,如下:
进行迭代运算获得初步标定结果,具体包括:计算所述相机的标定参数偏差矩阵,计算表达式为:P=[1-(1-α)(i+1)]H-1V,其中,i表示迭代次数,α表示平滑系数,α大于0小于1;β表示迭代精度系数,迭代精度系数β的表达式为:
其中,下标new表示计算得到的新参数,标定参数偏差矩阵P中前8项为内部参数偏差,按顺序依次包括等效焦距fu、fv,镜头成像中心像素坐标cu、cv,畸变系数k1、k2、k3、k4,之后诸项为多组外部参数;
根据更新得到的内部参数以及初始化时输入的各标定物图像的标定点数据,更新所述相机对应各标定物图像的外部参数,进行迭代运算直到满足预设条件,则结束迭代运算过程,获得所述相机的初步标定结果。
5.根据权利要求1-4任一项所述的三维扫描系统参数标定方法,其特征在于,还包括:获取所述三维扫描系统的结构光平面的平面方程,对所述三维扫描系统的结构光进行标定。
6.根据权利要求5所述的三维扫描系统参数标定方法,其特征在于,具体包括:
由三维扫描系统拍摄多幅经不同位置平面截取的结构光光线图像;
以获得的各相机的标定结果,对结构光光线图像进行畸变校正和极线校正,提取图像中结构光光线成像点的像素坐标;
进行空间位置重构,计算所述第一相机拍摄的图像中结构光光线成像点对应的实际物点的空间坐标(X1,Y1,Z1),计算表达式为:
其中:RT1为旋转平移矩阵,u1、v1表示结构光光线在所述第一相机中成像点像素坐标,u2、v2表示结构光光线在所述第二相机中成像点像素坐标,fu1、fv1表示所述第一相机的等效焦距,cu1、cv1表示所述第一相机的镜头成像中心像素坐标,fu2表示所述第二相机的u轴等效焦距,cu2表示所述第二相机图像的镜头成像中心u轴像素坐标;
利用在所述第一相机坐标系下重构出的成像点空间坐标(X1,Y1,Z1)拟合平面方程,所述平面方程的形式为Z=a+bX+cY,其中a,b,c为平面方程系数。
7.根据权利要求2所述的三维扫描系统参数标定方法,其特征在于,所述标定物包括平面标定图案,标定图案包括圆点和围绕圆点的带有编码信息的间断圆环,图案中圆点的圆心坐标作为标定点的世界坐标。
8.根据权利要求7所述的三维扫描系统参数标定方法,其特征在于,所述从所述标定物图像中识别出所述标定点,并提取所述标定点的像素坐标包括:
对标定物图像进行边缘检测,得到二值化的轮廓边缘图像;
对轮廓边缘图像进行椭圆拟合,获得轮廓椭圆的中心像素坐标;
对椭圆进行圆度判断和倾斜度判断,初步提取出有效的标定圆点像素坐标;
对标定圆点周围的编码圆环进行解码,得到标定圆点编号,确定标定圆点在标定板的位置,以与相应的世界坐标匹配;
对初步提取的标定圆点进行随机抽样,进行一致性判断,去除提取出的错误的标定圆点,得到最终有效的标定圆点及其像素坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710445013.3A CN107133989B (zh) | 2017-06-12 | 2017-06-12 | 一种三维扫描系统参数标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710445013.3A CN107133989B (zh) | 2017-06-12 | 2017-06-12 | 一种三维扫描系统参数标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107133989A CN107133989A (zh) | 2017-09-05 |
CN107133989B true CN107133989B (zh) | 2020-11-06 |
Family
ID=59734330
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710445013.3A Active CN107133989B (zh) | 2017-06-12 | 2017-06-12 | 一种三维扫描系统参数标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107133989B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6879168B2 (ja) * | 2017-11-01 | 2021-06-02 | オムロン株式会社 | 3次元測定装置、3次元測定方法及びプログラム |
CN109872368B (zh) * | 2017-12-04 | 2023-05-16 | 广东虚拟现实科技有限公司 | 图像处理方法、装置及测试系统 |
CN108805976B (zh) * | 2018-05-31 | 2022-05-13 | 武汉中观自动化科技有限公司 | 一种三维扫描系统及方法 |
CN108921902B (zh) * | 2018-06-06 | 2020-05-08 | 贝壳技术有限公司 | 一种修正结构光相机偏差的方法及装置 |
CN109556578A (zh) * | 2018-12-06 | 2019-04-02 | 成都天睿特科技有限公司 | 一种无人机旋扫式测量拍摄方法 |
TWI680436B (zh) | 2018-12-07 | 2019-12-21 | 財團法人工業技術研究院 | 深度相機校正裝置及其方法 |
CN109489553B (zh) * | 2018-12-27 | 2020-10-16 | 中国科学院长春光学精密机械与物理研究所 | 一种空间标志点库的生成方法、装置、设备及存储介质 |
CN111736235A (zh) | 2019-03-25 | 2020-10-02 | 同方威视技术股份有限公司 | Ct设备的几何参数标定件及标定方法 |
TWI720447B (zh) * | 2019-03-28 | 2021-03-01 | 財團法人工業技術研究院 | 影像定位方法及其系統 |
CN110135376A (zh) * | 2019-05-21 | 2019-08-16 | 北京百度网讯科技有限公司 | 确定图像传感器的坐标系转换参数的方法、设备和介质 |
CN110580724B (zh) * | 2019-08-28 | 2022-02-25 | 贝壳技术有限公司 | 一种对双目相机组进行标定的方法、装置和存储介质 |
CN111521132B (zh) * | 2020-04-14 | 2021-04-30 | 长春理工大学 | 一种新型自标定激光扫描投影方法 |
CN111973204B (zh) * | 2020-08-04 | 2022-09-06 | 上海涛影医疗科技有限公司 | 一种纳入重力的新型双平板x光机的校准方法 |
CN112419421B (zh) * | 2020-09-25 | 2022-11-04 | 河海大学 | 一种相机标定方法和装置 |
CN112258593B (zh) * | 2020-10-23 | 2024-02-02 | 明峰医疗系统股份有限公司 | 单目相机下ct或pet-ct智能定位扫描方法 |
CN112489109B (zh) * | 2020-11-19 | 2022-10-21 | 广州视源电子科技股份有限公司 | 一种三维成像系统方法、装置及三维成像系统 |
CN112634377A (zh) * | 2020-12-28 | 2021-04-09 | 深圳市杉川机器人有限公司 | 扫地机器人的相机标定方法、终端和计算机可读存储介质 |
CN113269836A (zh) * | 2021-03-30 | 2021-08-17 | 深圳市世宗自动化设备有限公司 | 3d相机标定方法、装置、计算机设备及其存储介质 |
CN113298883A (zh) * | 2021-06-08 | 2021-08-24 | 清德智体(北京)科技有限公司 | 用于对多个相机进行标定的方法、电子设备和存储介质 |
CN114037768A (zh) * | 2021-12-29 | 2022-02-11 | 天远三维(天津)科技有限公司 | 一种多套跟踪式扫描仪联合标定的方法和装置 |
CN116499364B (zh) * | 2023-06-30 | 2023-09-12 | 济南作为科技有限公司 | 盘煤仪三维激光点云调畸变方法及系统 |
CN116852382A (zh) * | 2023-09-04 | 2023-10-10 | 青岛理工大学 | 一种轴孔装配机器人末端姿态快速调整系统和方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5173796A (en) * | 1991-05-20 | 1992-12-22 | Palm Steven G | Three dimensional scanning system |
ES2327212T3 (es) * | 2000-08-25 | 2009-10-27 | 3Shape A/S | Procedimiento y aparato para un escaneado optico tridimensional de superficies interiores. |
CN100583151C (zh) * | 2006-09-22 | 2010-01-20 | 东南大学 | 三维扫描系统中双摄像机标定方法 |
CN101763643A (zh) * | 2010-01-07 | 2010-06-30 | 浙江大学 | 一种结构光三维扫描仪系统自动标定方法 |
CN102663767B (zh) * | 2012-05-08 | 2014-08-06 | 北京信息科技大学 | 视觉测量系统的相机参数标定优化方法 |
CN104299261B (zh) * | 2014-09-10 | 2017-01-25 | 深圳大学 | 人体三维成像方法及系统 |
CN105551039B (zh) * | 2015-12-14 | 2017-12-08 | 深圳先进技术研究院 | 结构光三维扫描系统的标定方法及装置 |
CN106780628A (zh) * | 2016-12-24 | 2017-05-31 | 大连日佳电子有限公司 | 基于混合畸变模型的高精度摄像机标定方法 |
-
2017
- 2017-06-12 CN CN201710445013.3A patent/CN107133989B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107133989A (zh) | 2017-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107133989B (zh) | 一种三维扫描系统参数标定方法 | |
CN110276808B (zh) | 一种单相机结合二维码测量玻璃板不平度的方法 | |
CN110288642B (zh) | 基于相机阵列的三维物体快速重建方法 | |
Furukawa et al. | Accurate camera calibration from multi-view stereo and bundle adjustment | |
Yang et al. | Robust and accurate surface measurement using structured light | |
CN113409382B (zh) | 车辆损伤区域的测量方法和装置 | |
CN110044374B (zh) | 一种基于图像特征的单目视觉测量里程的方法及里程计 | |
CN109580649B (zh) | 一种工程结构表面裂缝识别与投影修正方法及系统 | |
CN107588721A (zh) | 一种基于双目视觉的零件多尺寸的测量方法及系统 | |
CN113129430B (zh) | 基于双目结构光的水下三维重建方法 | |
CN109579695B (zh) | 一种基于异构立体视觉的零件测量方法 | |
CN111784778A (zh) | 基于线性求解非线性优化的双目相机外参标定方法和系统 | |
CN111189416B (zh) | 基于特征相位约束的结构光360°三维面形测量方法 | |
CN116309829B (zh) | 一种基于多目视觉的长方体扫描体群解码和位姿测量方法 | |
CN111524195A (zh) | 一种掘进机截割头定位中的相机标定方法 | |
CN113505626A (zh) | 一种快速三维指纹采集方法与系统 | |
CN113129384A (zh) | 基于一维编码靶标的双目视觉系统灵活标定方法 | |
CN114170284B (zh) | 基于主动标志点投射辅助的多视图点云配准方法 | |
CN106500625A (zh) | 一种远心立体视觉测量装置及其应用于物体三维形貌微米级精度测量的方法 | |
CN113963067B (zh) | 一种采用小靶标对大视场视觉传感器进行标定的标定方法 | |
CN111028280A (zh) | 井字结构光相机系统及进行目标有尺度三维重建的方法 | |
CN110838146A (zh) | 一种共面交比约束的同名点匹配方法、系统、装置及介质 | |
Wijenayake et al. | Dual pseudorandom array technique for error correction and hole filling of color structured-light three-dimensional scanning | |
CN107610215B (zh) | 一种高精度多角度口腔三维数字成像模型构建方法 | |
CN110992259B (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 |