CN105551039B - 结构光三维扫描系统的标定方法及装置 - Google Patents

结构光三维扫描系统的标定方法及装置 Download PDF

Info

Publication number
CN105551039B
CN105551039B CN201510925237.5A CN201510925237A CN105551039B CN 105551039 B CN105551039 B CN 105551039B CN 201510925237 A CN201510925237 A CN 201510925237A CN 105551039 B CN105551039 B CN 105551039B
Authority
CN
China
Prior art keywords
index
mark
dimensional
plane
scanning system
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
CN201510925237.5A
Other languages
English (en)
Other versions
CN105551039A (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.)
Song Zhan
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
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 Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201510925237.5A priority Critical patent/CN105551039B/zh
Priority to PCT/CN2015/098936 priority patent/WO2017101150A1/zh
Publication of CN105551039A publication Critical patent/CN105551039A/zh
Application granted granted Critical
Publication of CN105551039B publication Critical patent/CN105551039B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis

Abstract

本发明公开了一种结构光三维扫描系统的标定方法及装置,该方法包括:根据初始标定参数,对标定平面扫描,得到标定平面的三维点云数据;标定平面上设置有多个均匀分布的标记点;根据第一指标、第二指标和第三指标与标定参数的关系,找到最优标定参数;第一指标为:根据标记点对应的三维点云数据,计算得到的每相临两个标记点之间的距离与实际距离差异的平均值;第二指标为:根据标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与实际夹角角度差异的平均值;第三指标为:根据标定平面的三维点云数据,计算得到的标定平面的平面度与实际平面度的差异。上述技术方案提高了标定结果的准确性。

Description

结构光三维扫描系统的标定方法及装置
技术领域
本发明涉及三维扫描系统标定技术领域,特别涉及一种结构光三维扫描系统的标定方法及装置。
背景技术
目前,基于投影仪和相机架构的结构光三维扫描系统是最为广泛的非接触式三维测量手段,该技术基于立体视觉中的三角测量原理,通过投影仪投射特定编码光学图案,由相机获取投射图像,通过图像解码方法获取投射特征,进而建立起投影仪与相机的匹配关系,实现三维重建过程。
基于投影仪和相机架构的结构光三维扫描系统进行三维重建,涉及到的首个关键问题就是立体视觉系统的标定问题,即如何获取投影仪、相机的内部参数及外部参数,进而,由该参数建立三角测量函数,实现三维坐标的计算。目前常用的标定策略大概有两种:第一种是利用三维标定物进行标定,该标定物具有精确的已知三维信息,但该方法的缺点是三维标定物需要特殊制作,且精度要求很高,因此应用较少;第二种方法是用标定棋盘格,以棋盘格的角点为标定特征,该方法包括:首先,对相机进行标定,然后,根据相机标定结果,由投影仪投射出若干特征点到棋盘格平面上,计算出投射的特征点的三维信息,再对投影仪进行标定,因为这种方法比较易于实现,通过打印的棋盘格贴在平面物体上即可操作,因而为业内广泛使用。
基于棋盘格的结构光三维扫描系统标定方法的基本策略是:使用精确打印的棋盘格图像,将其粘贴在一个标准平面物体表面,首先,将棋盘格放在某个位置,由相机拍摄棋盘格图像,通过角点检测算法提取出角点信息,进而由投影仪投射出一个棋盘格图像,再通过角点检测算法检测出投射的棋盘格的角点信息;然后改变棋盘格标定物的姿态或距离,如此重复多次,获取足够多的标定图像信息(常用需15个不同姿态或位置的图像对)。在标定过程中,首先基于打印棋盘格图像对相机完成标定过程,之后就可以估算出棋盘格的空间平面三维信息,进而通过对投射棋盘格角点在空间平面上的三维信息进行计算,利用投射棋盘格在投影仪平面上的已知角点信息,即可完成对投影仪的标定,以及相机与投影仪的外部参数的估计。
中国专利申请号为CN201410164584的申请即公开了一种高精度的投影仪-摄像机标定系统及标定方法,该发明主要的方法是:用摄像机标定方法标定摄像机,得到摄像机内参数;向标定板投影纯白图案并与标定板的图案叠加,捕获标定区域图像;利用摄像机内参数对图像畸变校正后提取标定区域图像中的角点坐标;根据角点的对应关系估计摄像机像平面和标定板平面之间的单应性矩阵;依序向标定板投影不同特定棋盘图案并分别与标定板的图案叠加,分别捕获标定区域图像;对标定区域图像差分、滤波处理后提取标定板平面上的角点坐标;取角点坐标平均值后应用单应性矩阵映射角点至标定板平面;按角点的获取情况重复上述步骤,利用摄像机标定方法标定投影仪。这种标定方法已为业内广泛使用和研究,在此不再赘述。
然而,基于棋盘格的结构光三维扫描系统标定方法存在的主要问题包括:1)棋盘格需要精确的打印和制作,保证精确的尺寸,且棋盘格角点检测过程本身存在的检测误差也会对最终的标定结果产生影响;2)粘贴棋盘格的平面需要较高的平面度,因为这种标定方法是建立在平面假设基础之上的;3)棋盘格所包含的角点数量是有限的,其覆盖范围很难囊括整个投影范围和相机拍摄范围,没有角点存在的图像区域,其畸变参数等是难以估计的;4)现有标定方法是以角点在平面上的重映射误差最小为最终优化准则的,在实际测量中,这种误差准则与三角测量的结果存在不一致性,即使重投影误差很小也难以保证最终测量结果达到最优误差;5)标定过程和操作过程存在的人为操作误差导致结果的不确定性增加。
综上所述,现有标定方法存在的问题主要包括:1)标定过程有诸多些问题需要注意,如标定板的制作,标定过程摆放的位置,标定图像的数量等等,如果经验不足,很容易得到较差的标定参数;2)标定参数优化是基于检测到角点到标定平面的映射2D图像坐标误差最小为准则的,而实际测量中是以三维坐标距离等三维尺度为标准的,因此标定优化误差很难反映出实际的标定参数优劣。因此,利用现有标定方法对结构光三维扫描系统标定,标定结果不准确。
发明内容
本发明实施例提供了一种结构光三维扫描系统的标定方法,用以提高标定结果的准确性,该方法包括:
根据初始标定参数,对标定平面扫描,得到标定平面的三维点云数据;所述标定平面上设置有多个均匀分布的标记点;
根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数;所述第一指标为:根据所述标记点对应的三维点云数据,计算得到的每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;所述第二指标为:根据所述标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;所述第三指标为:根据标定平面的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异。
本发明实施例还提供了一种结构光三维扫描系统的标定装置,用以提高标定结果的准确性,该装置包括:
三维点云数据获取模块,用于根据初始标定参数,对标定平面扫描,得到标定平面的三维点云数据;所述标定平面上设置有多个均匀分布的标记点;
最优标定参数计算模块,用于根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数;所述第一指标为:根据所述标记点对应的三维点云数据,计算得到的每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;所述第二指标为:根据所述标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;所述第三指标为:根据标定平面的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异。
与现有技术相比较,本发明实施例提供的技术方案至少具有以下有益技术效果:
首先,本发明技术方案,不需要像现有技术中获取足够多的标定图像信息,仅需获得初始的标定参数,在该初始的标定参数的基础上,根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数即可,这样本发明技术方案不仅实施高效简单,而且标定结果准确可靠;
其次,通过建立一个标定平面,该标定平面上设置有多个均匀分布的标记点,采用初始标定参数,对该标定平面进行一次性完整扫描,得到标定平面上的三维点云数据;然后,在三维点云数据中找到每个标记点对应的三维点云数据;接着,根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数,第一指标和第二指标是根据标记点对应的三维点云数据得到的,第三指标也是根据标定平面的三维点云数据进行计算得到的,也就是说,寻找最优标定参数的过程是在三维空间中进行的,这样可以提高标定结果的准确性;另外,寻找结构光三维扫描系统的最优标定参数考虑了如下三个指标:第一指标:每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;第二指标:每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;以及第三指标:标定平面的平面度与所述标定平面的实际平面度的差异,这样保证了获取最优的系统标定参数。
通过上述可知,本发明实施例提供的技术方案,提高了标定结果的准确性,进而提高了结构光三维扫描系统的三维测量精度和可靠性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明实施例中结构光三维扫描系统的标定方法的流程示意图;
图2是本发明实施例中使用的标定平面的示意图;
图3是本发明实施例中结构光三维扫描系统的标定装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
在实际工作中,发明人发现现有结构光三维扫描系统的标定方法存在的问题主要有两点:1)标定过程有诸多些问题需要注意,如标定板的制作,标定过程摆放的位置,标定图像的数量等等,如果经验不足,很容易得到较差的标定参数;2)现有方法的标定参数优化是基于检测到角点到标定平面的映射2D图像坐标误差最小为准则的,而实际测量中是以三维坐标距离等三维尺度为标准的,因此标定优化误差很难反映出实际的标定参数优劣。
由于发明人发现了上述技术问题,本发明提出了一种结构光三维扫描系统的标定方法,所使用的优化目标为一标准平面(如玻璃白板),平面上预先标记出若干个标记点,标记点之间的尺寸为精确打印所得,通过对该标定板的一次性3D扫描,通过初始标定参数,获取其三维点云数据,通过第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数,整个过程简单易行,得到的最优标定参数能保证整个扫描系统在测量精度上有较大的提升,结果稳定可靠,对初始标定参数的精度要求不高,可以广泛应用于现有的结构光系统标定过程,用于提高标定结果的准确性和稳定性。下面进行详细说明。
图1是本发明实施例中结构光三维扫描系统的标定方法的流程示意图;如图1所示,该方法包括如下步骤:
步骤101:根据初始标定参数,对标定平面扫描,得到标定平面的三维点云数据;所述标定平面上设置有多个均匀分布的标记点;
步骤102:根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数;所述第一指标为:根据所述标记点对应的三维点云数据,计算得到的每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;所述第二指标为:根据所述标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;所述第三指标为:根据标定平面的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异。
具体实施时,在上述步骤101中,首先要获取结构光三维扫描系统的初始标定参数,其可以通过如下步骤实现:
参数定义如下:
c-相机,p-投影仪;
m—图像坐标,M-三维坐标,(相对于相机为Mc,相对于投影仪为Mp);
k-径向畸变,R,T—相机投影仪之间旋转,平移参数。
用Mc(p)=[Xc(p),Yc(p),Zc(p)]T表示物体上某个点的三维空间坐标,其中下标c(p),c代表的是相对于摄像机(相机)坐标系,p代表的是相对于投影仪坐标系。
那么,我们根据投射投影模型,可以将其图像坐标表示为
在考虑了径向畸变和切向误差之后有,对图像坐标进行反畸变处理,
其中,
而式子(2)中,代表的是切向畸变误差。
那么,最后可以得到:
其中,K是相机/投影仪的内部参数矩阵,即:
fc(p)x,fc(p)y,表示x-y方向像素尺度因子,r表示像素平面X-Y轴夹角畸变因子,cc(p)x,cc(p)y,表示图像中心点坐标。
上文中说过我们还标定摄像机和投影仪的相对位置,那么有:
上述公式描述了基本的结构光系统中,相机-投影仪的模型和参数,在本发明中,我们首先利用传统的摄像机标定方法进行相机标定,之后通过利用摄像机去捕捉投影仪投射棋盘图案到平面上的图案。建立起投影仪坐标和摄像机坐标系之间的一个单应矩阵,将投影仪建立与相机一样的模型,我们可以得到相机-投影仪的内部、外部初始标定参数。
在得到结构光系统的标定结果之后,对需找到的对应点利用三角原理可以得到其三维空间点的位置是:
再利用式子(5)我们可以计算出其X,Y坐标。
基于上述方法,我们可以得到相机-投影仪的内部及外部参数,及实现每个图像点的三维坐标计算,在此基础上,我们完成一次对标定平面(图2)的扫描过程,得到高密度的点云数据和图2中标记点的中心点(本发明提到的计算三个指标时,都是根据该标记点的中心点的坐标来计算的)图像坐标及三维坐标,在此基础上,我们可以寻找最优标定参数,例如:标定参数的优化,即在获取初始标定参数之后,就可以寻找最优标定参数,例如:进行标定参数的优化,其大致过程可以如下:
采用初始标定参数,完成对标定平面的一次完整扫描,此时,可以获取该标定平面的三维点云数据;从扫描图像中检测出图2中所有标记点的图像坐标,在三维扫描数据中找到标记点的图像坐标对应的三维坐标;建立全局优化函数,该函数包括相机和投影机的内部参数以及全部外部参数,其优化指标包括三个:平面度:基于所有重建点的平面拟合误差,即根据标定平面上三维点云数据计算出的平面度与所述标定平面的实际平面度的差异;尺寸精度:检测出的标记点的尺寸误差,即每两个标记点之间的距离与实际距离的差异;角度误差:检测出的标记点的夹角误差,即每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与实际夹角角度的差异。最后,以初始标定参数为全局优化函数初始值,给定相关标定参数的上下限范围阈值,执行全局优化函数,直至所有误差指标其中之一或任意组合最小时最小,即可获得最优标定参数。
下面详细介绍上述优化标定参数的步骤:
具体实施时,本发明实施例中,根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数,可以通过上述提到的建立全局优化函数的方法得到,例如:
建立所述第一指标、第二指标和第三指标与标定参数的全局优化函数;
以所述初始标定参数为所述全局优化函数的初始值,迭代找到最优标定参数,对于每个迭代周期均执行以下操作:
根据所述标记点对应的三维点云数据,计算所述第一指标和第二指标;根据标定平面的三维点云数据,计算所述第三指标;直到找到第一指标、第二指标和第三指标中一个或多个最小时对应的标定参数,作为结构光三维扫描系统的最优标定参数。
本发明实施例提供的技术方案,不需要像现有技术中获取足够多的标定图像信息,仅需获得初始的标定参数,在该初始的标定参数的基础上,根据建立的全局优化函数,进行进一步的标定参数优化,找到三个指标的其中之一或任意组合最小时对应的标定参数,作为最优标定参数,这样本发明技术方案不仅实施高效简单,而且标定结果准确可靠。
当然,除了通过建立第一指标、第二指标和第三指标与标定参数的全局优化函数,来迭代找到最优标定参数的方式之外,还可以有多种方式,例如,可以通过建立一个第一指标、第二指标和第三指标与标定参数的表格等等方式,来找到最优标定参数。
具体实施时,上述“以所述初始标定参数为所述全局优化函数的初始值,迭代找到最优标定参数”可以包括:
以所述初始标定参数为全局优化函数的初始值,设定标定参数的范围阈值,迭代找到最优标定参数。
在进行优化的过程中,在初始标定参数基础上,设定的优化标定参数范围的上下限的优点是:可以有效避免优化进入局部最小,同时提高运行效率。
具体实施时,在上述步骤101中,标定平面的构建可以有多种方式,例如图2所示,标定平面可以包括:平面物体以及粘贴在平面物体上的、满足结构光三维扫描系统扫描范围的白纸;所述白纸上设置有多个均匀分布的标记点。
具体实施时,制作过程可以为:根据所需标定的结构光系统的扫描范围设计一张白色的网格点图像,将其打印后粘贴与平面物体上,其主要目的是在白色表面上标记出若干精确尺寸参考点。当然,标定平面也可以是特意制作的、设置有多个均匀分布的标记点的玻璃白板。
具体实施时,所述标记点的形状可以为圆形。标记点的形状可以为圆形的优点是:在计算上述提到的三个指标时,以计算两个标记点之间的距离为例,需要计算两个标记点的中心点之间的距离,同理,在计算其它两个指标时,也是以中心点为标准进行计算的,圆形的中心点就是圆心,三维空间里面是球心,方便寻找中心点。实际实施中,也可选择其它形状的标记点,只要方便精确计算三个指标即可。
本发明实施例提到的标定参数可以包括:结构光三维扫描系统中相机和摄影机的内部参数和全部外部参数。本发明实施例中的全局优化函数包含了所有的相机、投影仪的内部参数,以及12个外部参数(旋转矩阵9个,平移向量3个),其优化目标函数最小化准则包含三个指标,下面对这三个指标进行详细介绍:
a)平面度误差指标E1(即为:根据标定平面上所有点的三维点云数据,计算出的标定平面的平面度与所述标定平面的实际平面度的差异):
平面度具体的计算过程可以为:由于重建的点云数量巨大,可达几百万,为了提高优化函数的运行速度,具体实施时,可以根据从标定平面上均匀采样的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异,即可以为:对点云采用随机采样或者均匀采样,得到较为稀疏的点云数据,对这些点云数据进行最小二乘平面拟合,得到标定平面的平面度误差E1;
b)尺寸误差指标E2(即为:根据每个标记点对应的三维点云数据,计算出的每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值):
具体的计算过程可以为:由于标定平面上的标记点的尺寸为精确已知(在制作标定平面时设计的尺寸数据),因此,我们计算所有标记点之间的距离误差,取平均值,计算的距离值与真实距离值的误差E2;
c)角度误差E3(即为:每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值):
具体实施时,可以采用角度(根据三维点云数据,计算得到的)的余弦值与标准的角度(实际的夹角角度,如90度)进行比较后的差异,作为角度误差值:
E3=∑|cosθcal-cosθreal|,其中,θcal为根据三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度,θreal为所述两条直线之间的实际夹角角度(即,标定平面精确设计打印时的标准尺寸)。
基于上述三个误差标准,我们就可以构造出全局优化函数的目标函数如下:
Fobj=∑(Ed+α·Ep+β·Ea); (7)
其中,α,β满足:Ed≈α·Ep≈β·Ea;Ed为第一指标(尺寸误差):每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;Ea为第二指标(角度误差):根据所述标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;Ep为第三指标(平面度误差):根据标定平面的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异。α和β是经验参数,用来将三个误差调整到一个数量级,比如尺寸误差是0.1mm左右的话,平面度误差常在1mm左右,角度误差在0.1度左右的话,则可设置为α=0.1和β=1。
最后,我们可以建立全局优化函数如下:
minFobj(parasSLS)
s.t. (8)
其中,是在初始标定参数基础上,设定的标定参数的上限,是在初始标定参数基础上,设定的标定参数的下限,这样可以有效避免优化进入局部最小,同时提高运行效率。
具体实施时,在对上述全局优化函数(8)求解时:全局优化函数(8)属于典型的非线性、多目标优化问题,我们可以借助现有的优化工具(例如:MATLAB工具)对其进行优化,并最终获取能够满足三个优化准则(三个指标)对应的最优标定参数。在所建立的全局优化函数中,我们综合考虑了三个指标,具体实施,亦可根据实际需要,选择单个或任意个指标进行优化处理。
通过上述可知,本发明实施提供的技术方案,不需要像现有技术中获取足够多的标定图像信息,仅需获得初始的标定参数,在该初始的标定参数的基础上,根据建立的全局优化函数,进行进一步的标定参数优化,找到三个指标中一个或多个最小时对应的标定参数,作为最优标定参数,这样本发明技术方案不仅实施高效简单,而且标定结果准确可靠;
本发明实施提供的技术方案,通过建立一个标定平面,该标定平面上设置有多个均匀分布的标记点,采用初始标定参数,对该标定平面进行一次性完整扫描,得到标定平面上所有点的三维点云数据;然后,在三维点云数据中找到每个标记点对应的三维点云数据;接着,建立全局优化函数,以初始标定参数为全局优化函数的初始值,迭代计算找到三个指标的一个或多个最小时对应的标定参数,作为最优标定参数。计算三个指标是根据标定平面上所有点的三维点云数据和每个标记点对应的三维点云数据进行计算的,也就是说,优化标定参数的过程是在三维空间中进行的,这样可以提高标定结果的准确性;另外,全局优化函数的目标函数包括以下三个关键几何属性指标的其中之一或任意组合:标定平面的平面度,每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异、每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异,这样保证了获取最优的系统标定参数。
通过上述可知,本发明实施例提供的技术方案,提高了标定结果的准确性,进而提高了结构光三维扫描系统的三维测量精度和可靠性。
基于同一发明构思,本发明实施例还提供了一种结构光三维扫描系统的标定装置,如下面的实施例。由于结构光三维扫描系统的标定装置解决问题的原理与结构光三维扫描系统的标定方法相似,因此结构光三维扫描系统的标定装置的实施可以参见结构光三维扫描系统的标定方法的实施,重复之处不再赘述。以下所使用的,术语“单元”、“装置”或者“模块”可以实现预定功能的软件和/或硬件的组合。以下实施例所描述的装置可以以软件、硬件,或者软件和硬件的组合的实现,也是可能并被构想的。
图3是本发明实施例中结构光三维扫描系统的标定装置的结构示意图,如图3所示,该装置包括:
三维点云数据获取模块10,用于根据初始标定参数,对标定平面扫描,得到标定平面的三维点云数据;所述标定平面上设置有多个均匀分布的标记点;
最优标定参数计算模块20,用于根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数;所述第一指标为:根据所述标记点对应的三维点云数据,计算得到的每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;所述第二指标为:根据所述标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;所述第三指标为:根据标定平面的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异。
在一个实施例中,所述最优标定参数计算模块具体用于:
建立所述第一指标、第二指标和第三指标与标定参数的全局优化函数;
以所述初始标定参数为所述全局优化函数的初始值,迭代找到最优标定参数,对于每个迭代周期均执行以下操作:
根据所述标记点对应的三维点云数据,计算所述第一指标和第二指标;根据标定平面的三维点云数据,计算所述第三指标;直到找到第一指标、第二指标和第三指标中一个或多个最小时对应的标定参数,作为结构光三维扫描系统的最优标定参数。
经过发明人大量的实验认证,本发明实施例效果很理想,和设计的预期一致。本发明实施例提供的技术方案的有益技术效果为:
1)使用了一个简单的带有若干尺寸标记点的标定平面作为结构光三维扫描系统标定参数优化的实验目标,实施简单,结果可靠。
2)提出了一种全局参数优化算法,基于给定的标定平面和标记点参数,建立了包含平面度、尺寸、角度三个目标函数的全局优化函数,通过一次平面扫描数据,即可完成对全部参数的优化过程。
3)与现有的标定参数优化方面不同的是,本发明的优化方法是在实测的三维空间中进行的,而且其优化参数包含了所有的系统参数,确保能够获取最优的系统标定参数,将系统测量误差降到最低。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块、装置或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种结构光三维扫描系统的标定方法,其特征在于,包括:
根据初始标定参数,对标定平面扫描,得到标定平面的三维点云数据;所述标定平面上设置有多个均匀分布的标记点;
根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数;所述第一指标为:根据所述标记点对应的三维点云数据,计算得到的每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;所述第二指标为:根据所述标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;所述第三指标为:根据标定平面的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异;
根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数,包括:
建立所述第一指标、第二指标和第三指标与标定参数的全局优化函数;
以所述初始标定参数为所述全局优化函数的初始值,迭代找到最优标定参数,对于每个迭代周期均执行以下操作:
根据所述标记点对应的三维点云数据,计算所述第一指标和第二指标;根据标定平面的三维点云数据,计算所述第三指标;直到找到第一指标、第二指标和第三指标中一个或多个最小时对应的标定参数,作为结构光三维扫描系统的最优标定参数。
2.如权利要求1所述的结构光三维扫描系统的标定方法,其特征在于,以所述初始标定参数为所述全局优化函数的初始值,迭代找到最优标定参数,包括:
以所述初始标定参数为全局优化函数的初始值,设定标定参数的范围阈值,迭代找到最优标定参数。
3.如权利要求1至2任一权利要求所述的结构光三维扫描系统的标定方法,其特征在于,所述标记点的形状为圆形。
4.如权利要求1至2任一权利要求所述的结构光三维扫描系统的标定方法,其特征在于,所述标定平面包括:平面物体以及粘贴在平面物体上的、满足结构光三维扫描系统扫描范围的白纸;所述白纸上设置有多个均匀分布的标记点。
5.如权利要求1至2任一权利要求所述的结构光三维扫描系统的标定方法,其特征在于,所述标定参数包括:结构光三维扫描系统中相机和摄影机的内部参数和外部参数。
6.如权利要求1至2任一权利要求所述的结构光三维扫描系统的标定方法,其特征在于,所述第三指标具体为:根据从标定平面上均匀采样的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异。
7.如权利要求1至2任一权利要求所述的结构光三维扫描系统的标定方法,其特征在于,所述全局优化函数为:
minFobj(parasSLS)
使得:
其中,Fobj为全局优化函数的目标函数:Fobj=∑(Ed+α·Ep+β·Ea);α,β满足:Ed≈α·Ep≈β·Ea是标定参数的上限,是标定参数的下限;Ed为第一指标:每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;Ea为第二指标:根据所述标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;Ep为第三指标:根据标定平面的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异;α和β是经验参数,用来调整Ed、Ea和Ep的数量级。
8.一种结构光三维扫描系统的标定装置,其特征在于,包括:
三维点云数据获取模块,用于根据初始标定参数,对标定平面扫描,得到标定平面的三维点云数据;所述标定平面上设置有多个均匀分布的标记点;
最优标定参数计算模块,用于根据第一指标、第二指标和第三指标与标定参数的关系,找到结构光三维扫描系统的最优标定参数;所述第一指标为:根据所述标记点对应的三维点云数据,计算得到的每相临两个标记点之间的距离与所述两个标记点之间的实际距离差异的平均值;所述第二指标为:根据所述标记点对应的三维点云数据,计算得到的每个标记点和它相临的两个标记点连接形成的两条直线之间的夹角角度与所述两条直线之间的实际夹角角度差异的平均值;所述第三指标为:根据标定平面的三维点云数据,计算得到的标定平面的平面度与所述标定平面的实际平面度的差异;
所述最优标定参数计算模块具体用于:
建立所述第一指标、第二指标和第三指标与标定参数的全局优化函数;
以所述初始标定参数为所述全局优化函数的初始值,迭代找到最优标定参数,对于每个迭代周期均执行以下操作:
根据所述标记点对应的三维点云数据,计算所述第一指标和第二指标;根据标定平面的三维点云数据,计算所述第三指标;直到找到第一指标、第二指标和第三指标中一个或多个最小时对应的标定参数,作为结构光三维扫描系统的最优标定参数。
CN201510925237.5A 2015-12-14 2015-12-14 结构光三维扫描系统的标定方法及装置 Active CN105551039B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201510925237.5A CN105551039B (zh) 2015-12-14 2015-12-14 结构光三维扫描系统的标定方法及装置
PCT/CN2015/098936 WO2017101150A1 (zh) 2015-12-14 2015-12-25 结构光三维扫描系统的标定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510925237.5A CN105551039B (zh) 2015-12-14 2015-12-14 结构光三维扫描系统的标定方法及装置

Publications (2)

Publication Number Publication Date
CN105551039A CN105551039A (zh) 2016-05-04
CN105551039B true CN105551039B (zh) 2017-12-08

Family

ID=55830215

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510925237.5A Active CN105551039B (zh) 2015-12-14 2015-12-14 结构光三维扫描系统的标定方法及装置

Country Status (2)

Country Link
CN (1) CN105551039B (zh)
WO (1) WO2017101150A1 (zh)

Families Citing this family (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106918302B (zh) * 2017-03-13 2019-09-27 东华大学 一种三维扫描仪自适应标定方法
CN107133989B (zh) * 2017-06-12 2020-11-06 中国科学院长春光学精密机械与物理研究所 一种三维扫描系统参数标定方法
JP7101616B2 (ja) * 2017-08-25 2022-07-15 クリス・シンライ・リウ 標的要素の位置を特定するためのマシンビジョンシステムおよび方法
CN108089198B (zh) * 2017-12-11 2023-12-19 同方威视技术股份有限公司 三维扫描装置、机器人及数据处理方法
CN109920003B (zh) * 2017-12-12 2023-09-15 广东虚拟现实科技有限公司 相机标定检测方法、装置以及设备
CN110440708A (zh) * 2018-05-04 2019-11-12 苏州玻色智能科技有限公司 一种用于三维白光扫描设备的标准件及其标定方法
CN108921902B (zh) * 2018-06-06 2020-05-08 贝壳技术有限公司 一种修正结构光相机偏差的方法及装置
CN109061609B (zh) * 2018-08-01 2021-09-07 歌尔光学科技有限公司 Tof模组标定装置及方法
JP7139953B2 (ja) * 2019-01-08 2022-09-21 オムロン株式会社 三次元形状測定装置、三次元形状測定方法及びプログラム
CN109977466B (zh) * 2019-02-20 2021-02-02 深圳大学 一种三维扫描视点规划方法、装置及计算机可读存储介质
CN111833392A (zh) * 2019-04-16 2020-10-27 杭州思看科技有限公司 标记点多角度扫描方法、系统及装置
CN110207622B (zh) * 2019-06-30 2021-04-20 武汉尺子科技有限公司 一种三维扫描仪标定方法及系统
CN110598239B (zh) * 2019-07-19 2023-11-17 中国建设基础设施有限公司 基于轨行区点云大数据的应用方法
CN110849389B (zh) * 2019-10-24 2023-08-29 西北大学 一种基于空间球的两个三维点云采集系统外参数标定方法
CN112750165B (zh) * 2019-10-29 2023-12-15 商汤集团有限公司 参数标定方法、智能驾驶方法及其装置、设备和存储介质
WO2021081963A1 (zh) * 2019-10-31 2021-05-06 深圳市大疆创新科技有限公司 标定板、测试角分辨率的方法、装置及计算机存储介质
CN111161339B (zh) * 2019-11-18 2020-11-27 珠海随变科技有限公司 一种间距测量方法、装置、设备以及计算机可读介质
CN110930382A (zh) * 2019-11-19 2020-03-27 广东博智林机器人有限公司 基于标定板特征点提取的点云拼接精度评估方法及系统
TWI722703B (zh) 2019-12-09 2021-03-21 財團法人工業技術研究院 投影設備與投影校正方法
CN111028280B (zh) * 2019-12-09 2022-06-03 西安交通大学 井字结构光相机系统及进行目标有尺度三维重建的方法
CN111179339B (zh) * 2019-12-13 2024-03-08 深圳市瑞立视多媒体科技有限公司 基于三角测量的坐标定位方法、装置、设备及存储介质
GB202305331D0 (en) 2019-12-18 2023-05-24 Motional Ad Llc Camera-to-lidar calibration and validation
CN111445534B (zh) * 2020-04-03 2023-04-21 北京四维智联科技有限公司 一种用于全自动泊车的相机高精度自动标定方法及装置
CN111486802B (zh) * 2020-04-07 2021-04-06 东南大学 基于自适应距离加权的旋转轴标定方法
CN111735487B (zh) * 2020-05-18 2023-01-10 清华大学深圳国际研究生院 传感器、传感器标定方法与设备、存储介质
CN111627113B (zh) * 2020-05-25 2023-09-19 易视智瞳科技(深圳)有限公司 3d线扫建模水平角度的标定方法、装置、设备及存储介质
CN112082512B (zh) * 2020-09-08 2023-04-14 深圳广成创新技术有限公司 一种相位测量偏折术的标定优化方法、装置及计算机设备
CN112270746B (zh) * 2020-11-06 2023-06-16 太原科技大学 基于邻域协方差特征参数阈值的铝合金3d打印点云精简算法
CN112581605A (zh) * 2020-12-24 2021-03-30 西安中科光电精密工程有限公司 一种结构光三维重建校正方法及装置
WO2023277889A1 (en) * 2021-06-29 2023-01-05 Hewlett-Packard Development Company, L.P. Calibration of 3d scanning systems
CN113487680B (zh) * 2021-06-29 2022-07-26 湖南泽塔科技有限公司 一种单幅图像内外参数标定方法及系统
CN113970297A (zh) * 2021-11-19 2022-01-25 江苏智库智能科技有限公司 堆垛机货位列数据的标定测量方法
CN114037768A (zh) * 2021-12-29 2022-02-11 天远三维(天津)科技有限公司 一种多套跟踪式扫描仪联合标定的方法和装置
CN115937069B (zh) * 2022-03-24 2023-09-19 北京小米移动软件有限公司 零件检测方法、装置、电子设备及存储介质
CN116772746A (zh) * 2023-08-17 2023-09-19 湖南视比特机器人有限公司 利用射灯图案检测进行平面度轮廓度测量方法及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101673399A (zh) * 2009-09-29 2010-03-17 浙江工业大学 编码结构光三维视觉系统的标定方法
CN103985157A (zh) * 2014-05-30 2014-08-13 深圳先进技术研究院 一种结构光三维扫描方法及系统

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN100573586C (zh) * 2008-02-21 2009-12-23 南京航空航天大学 一种双目立体测量系统的标定方法
CN100573040C (zh) * 2008-03-17 2009-12-23 南京航空航天大学 物体表面三维轮廓结构光测量系统的标定方法
CN101419061B (zh) * 2008-12-08 2011-06-29 北京航空航天大学 一种镜像式结构光视觉测量系统和测量方法
CN101697233B (zh) * 2009-10-16 2012-06-06 长春理工大学 一种基于结构光的三维物体表面重建方法
CN101763643A (zh) * 2010-01-07 2010-06-30 浙江大学 一种结构光三维扫描仪系统自动标定方法
US10008007B2 (en) * 2012-09-20 2018-06-26 Brown University Method for generating an array of 3-D points
CN103942796B (zh) * 2014-04-23 2017-04-12 清华大学 一种高精度的投影仪‑摄像机标定系统及标定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101673399A (zh) * 2009-09-29 2010-03-17 浙江工业大学 编码结构光三维视觉系统的标定方法
CN103985157A (zh) * 2014-05-30 2014-08-13 深圳先进技术研究院 一种结构光三维扫描方法及系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Improving calibration accuracy of structured light systems using planebased residual error compensation;AD Han et al.,;《Optical Engineering》;20131031;第52卷(第10期);摘要,第1、3-4章 *
Novel calibration method for structured-light system with an out-of-focus projector;Beiwen Li et al.;《Applied Optics》;20140523;第53卷(第16期);第3415-3426页 *
Optimal Image-Based Euclidean Calibration ofStructured Light Systems in General Scenes;Nathaniel Bird et al.;《IEEE TRANSACTIONS ON AUTOMATION SCIENCE AND ENGINEERING》;20111031;第8卷(第4期);第815-823页 *
结构光测量系统误差分析;许丽 等;《光电工程》;20081031;第35卷(第10期);第26-31页 *

Also Published As

Publication number Publication date
CN105551039A (zh) 2016-05-04
WO2017101150A1 (zh) 2017-06-22

Similar Documents

Publication Publication Date Title
CN105551039B (zh) 结构光三维扫描系统的标定方法及装置
CN108151671B (zh) 一种三维数字成像传感器、三维扫描系统及其扫描方法
CN109000582A (zh) 跟踪式三维扫描装置的扫描方法及系统、存储介质、设备
CN104517291B (zh) 基于目标同轴圆特征的位姿测量方法
CN106780619B (zh) 一种基于Kinect深度相机的人体尺寸测量方法
CN104835144B (zh) 利用一个球的球心的像及正交性求解摄像机内参数的方法
CN102364299B (zh) 一种多个结构光投影三维型面测量头的标定技术
CN108288292A (zh) 一种三维重建方法、装置及设备
CN102589530B (zh) 基于二维相机和三维相机融合的非合作目标位姿测量方法
CN107133989A (zh) 一种三维扫描系统参数标定方法
CN103106661B (zh) 空间二条相交直线线性求解抛物折反射摄像机内参数
CN102054276B (zh) 物体三维几何重建相机定标方法和系统
CN109816724A (zh) 基于机器视觉的三维特征提取方法及装置
CN103106662B (zh) 空间一条直线线性求解抛物折反射摄像机内参数
CN106990776A (zh) 机器人归航定位方法与系统
He et al. Eccentricity error identification and compensation for high-accuracy 3D optical measurement
CN109360230A (zh) 一种基于2d相机与3d相机的图像配准方法及系统
CN107886547A (zh) 一种鱼眼相机标定方法及系统
CN107816942A (zh) 一种基于十字结构光视觉系统的平面二维尺寸测量方法
CN109544628A (zh) 一种指针式仪表的准确读数识别系统及方法
CN102930551B (zh) 利用圆心的投影坐标和极线求解摄像机内参数
CN104200476B (zh) 利用双平面镜装置中的圆周运动求解摄像机内参数的方法
CN109272555A (zh) 一种rgb-d相机的外部参数获得及标定方法
CN105574812A (zh) 多角度三维数据配准方法及装置
CN102914295A (zh) 基于计算机视觉立方体标定的三维测量方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20190116

Address after: 516000 West Side of the 4th Floor of CD Building, No. 2 South Road, Huatai Road, Huiao Avenue, Huizhou City, Guangdong Province

Patentee after: HUIZHOU GOVION TECHNOLOGY CO.,LTD.

Address before: 518055 No. 1068, Xue Yuan Avenue, Xili University Town, Nanshan District, Shenzhen, Guangdong

Patentee before: SHENZHEN INSTITUTES OF ADVANCED TECHNOLOGY

CP03 Change of name, title or address

Address after: 215163 rooms 101, 102, 901 and 902, floor 1, building 11, 198 Jialingjiang Road, high tech Zone, Suzhou City, Jiangsu Province

Patentee after: Gaoshi Technology (Suzhou) Co.,Ltd.

Address before: 516000 West Side of the 4th Floor of CD Building, No. 2 South Road, Huatai Road, Huiao Avenue, Huizhou City, Guangdong Province

Patentee before: HUIZHOU GOVION TECHNOLOGY CO.,LTD.

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 215129 Rooms 101, 102, 901, 902, Floor 9, Building 11, No. 198, Jialing River Road, High tech Zone, Suzhou City, Jiangsu Province

Patentee after: Gaoshi Technology (Suzhou) Co.,Ltd.

Address before: 215163 rooms 101, 102, 901 and 902, floor 1, building 11, 198 Jialingjiang Road, high tech Zone, Suzhou City, Jiangsu Province

Patentee before: Gaoshi Technology (Suzhou) Co.,Ltd.

CP03 Change of name, title or address
TR01 Transfer of patent right

Effective date of registration: 20240111

Address after: 518000 No. 1068 Xueyuan Avenue, Nanshan District, Shenzhen, Guangdong Province

Patentee after: Song Zhan

Address before: 215129 Rooms 101, 102, 901, 902, Floor 9, Building 11, No. 198, Jialing River Road, High tech Zone, Suzhou City, Jiangsu Province

Patentee before: Gaoshi Technology (Suzhou) Co.,Ltd.

TR01 Transfer of patent right