CN101576379B - 基于二维彩色标靶的主动投影三维测量系统快速标定方法 - Google Patents
基于二维彩色标靶的主动投影三维测量系统快速标定方法 Download PDFInfo
- Publication number
- CN101576379B CN101576379B CN2009100592516A CN200910059251A CN101576379B CN 101576379 B CN101576379 B CN 101576379B CN 2009100592516 A CN2009100592516 A CN 2009100592516A CN 200910059251 A CN200910059251 A CN 200910059251A CN 101576379 B CN101576379 B CN 101576379B
- Authority
- CN
- China
- Prior art keywords
- image
- color
- projector
- target
- camera
- 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
Abstract
一种用于主动投影照明三维形貌测量系统中的投影仪和摄像机同时快速标定方法,利用合理的颜色设计,在彩色图像的RGB模型基础上,利用投影仪将一幅具有彩色特征点的图像投影到一个具有彩色特征点的二维平面标靶上,彩色摄像机一次同时获取标靶上两个彩色图案经颜色空间混色后的图像,由图像分离得到标靶上制作的特征点图像信息和投影的特征点图像信息,然后分别用分离出的标靶特征图像标定摄像机,用分离出的投影特征点图像标定投影仪,最后利用双目立体视觉原理标定出摄像机和投影仪的相对位置关系,从而完成整个主动照明三维形貌测量系统的标定。由于在每一个标定位置只需要拍摄一幅图像,该方法具有测量数据少,标定速度快的优点。
Description
技术领域
本发明涉及三维面形测量系统中的系统标定领域,尤其涉及基于主动投影照明三维测量系统的摄像机和投影仪同时标定。
背景技术
摄像机标定是三维重建中必不可少的一个步骤,标定结果的好坏直接影响着三维测量的精度和三维重建结果的好坏,因此研究摄像机的标定方法具有重要的理论研究意义和实际应用价值。
摄像机模型解决的是二维图像信息和三维空间坐标转换的问题,一般来说有线性模型和非线性模型两种。线性摄像机模型(针孔相机成像模型)是应用最广泛的一种,很多视觉测量技术都是建立在针孔成像模型的基础之上。非线性模型也称为畸变模型,是在线性模型的基础上考虑成像镜头的光学畸变后建立的模型。
(1)线性模型摄像机的标定
摄像机定标,也就是求解其成像投影变换矩阵M,一般都需要在摄像机前放一个特制的标定参照物,摄像机获取该物体的图像,并由此计算摄像机的投影变换矩阵M。标定参照物的每一个特征点相对于世界坐标系GCS的位置在制作时应精确测定,世界坐标系GCS可选为参照物的物体坐标系。在得到这些已知点在图像上的投影位置后,可以计算出摄像机的投影变换矩阵M,之后无须计算出摄像机的全部内外部参数,采用隐式校准的办法可以直接将像面坐标映射成空间坐标。
其中,(Xwi,Ywi,Zwi)为空间第i个点的坐标,(ui,vi,1)为第i点的图像像素坐标,mij为投影矩阵M的第i行第i列元素。公式1包含了三个方程:
将公式2中的第一式除第三式,第二式除第三式分别消去Zci后,可得如下两个关于mij的线性方程:
公式3表示,若标定块上有n个已知点,并已知它们的空间坐标(Xwi,,Ywi,Zwi)(i=1,2,...,n)与它们的图像点坐标(ui,vi)(i=1,2,...,n),则就可获得2n个关于M矩阵元素的线性方程,其矩阵形式的方程为:
由4式可见,M矩阵乘以任意不为零的常数并不影响(Xw,,Yw , Zw)与(u,v)的关系,因此,在上式中可以指定m34=1,从而得到关于M矩阵其它元素的2n个线性方程,将4式简写成:
Km=U 5
其中,K为上式左边2n×11矩阵,m为未知的11维向量,U为右边的2n维向量,K,U为已知向量。当2n>11时,可用最小二乘法求出上述线性方程的解:
m=(KTK)-1KTU 6
m向量与m34=1构成了所求解的M矩阵。由上可见,由空间6个以上已知点与它们的图像点坐标,可以求出M矩阵。在一般的定标工作中,使标定块上有数十个已知点,使方程的个数大大超过未知数的个数,从而用最小二乘法求解以降低误差造成的影响。
(2)张正友的平面标定法
线性模型不能准确的描述成像几何关系,在精度要求比较高的场合需要考虑非线性的影响,即需要采用非线性模型。非线性模型考虑了镜头的光学畸变对成像的影响,故使精度有较大的提高。非线性模型需要解非线性方程,从而增大了解的难度和解的不确定度,从而也限制了非线性模型的应用。
基于上面的考虑,Tsai和Roger提出了著名的两步法,该法先采用解析方法直接计算线性部分参数,然后以这些线性参数作为非线性优化的初值,对其余参数进行迭代优化,故称为两步法。该法一方面克服了传统的线性模型的不足,考虑了径向畸变,提高了标定精度,另一方面又通过解析法得到初值,从而减少了优化的次数,提高了运算速度。
1998年张正友针对径向畸变,提出了一种可以利用多幅平面模板标定摄像机所有内外参数的方法。该方法需要摄像机从不同的角度拍摄平面模板的多幅图像(至少三幅),不过摄像机与平面模板间可以自由地运动,而且运动参数无需知道。由于平面模板上每个特征点与其图像上相应的像点之间存在一个对应关系,这个关系可用一个映射矩阵来表示,所以对于每幅图像,就都可以确定一个映射矩阵,这就为内部参数的求解提供了两个约束条件。该算法也是基于 两步法的思想,即先有一个线性解法求出部分参数的初始值,然后考虑径向畸变(一阶和二阶)并基于极大似然准则对线性结果进行非线性优化最后利用计算好的内部参数和平面模板映射矩阵求出外部参数。
在主动投影照明式光学三维面形测量系统中常用投影仪进行光栅图像的投影,其具有相移准确、精度高等优点。这样的主动投影照明式双目立体视觉三维测量系统,投影仪的标定就成为一个新的问题。
文献“Novel method for structured light system calibration”(Song Zhang,PeisenS.Huang,Optical Engineering,2006)提出一种投影仪标定方法,其将投影仪看成是摄像机的倒置,这样投影仪标定模型就和摄像机相同。其利用投影仪投射出的正弦条纹的相位信息,按照张正友的平面方法来实现系统的标定。该方法需要在两个垂直方向上投影正弦条纹且使用条纹相移术,故在每一个平面标定位置至少需要六幅条纹图。此方法需要进行相位展开运算,数据处理较繁琐。
发明内容
本发明提供一种基于二维彩色标靶的主动投影照明光学三维面形测量系统快速标定方法。本发明适合现场快速标定,可应用于基于结构光主动投影的三维测量系统中。
本发明采用如下技术方案:
标定靶上的图案由三基色中任意两种颜色组成,其一为背景色,另一为特征点的颜色。使用打印、印刷、光电制作或者机械制作方式获得特征图案,然后将该图像附着在平面上即制作成标靶。
根据标靶上的颜色组成,由颜料和光的混色原理,投影图像中白色背景上的 特征点颜色应该为三基色中标靶图案不包含的那种颜色的补色。
将投影图像投影到标靶上,投影图像的光颜色和标靶图案本身的颜料颜色按照光和颜料的反射原理形成一个图案,由RGB模型的图像记录原理,利用事先的颜色设计,标靶上物理存在的特征点信息会包含在组成标靶颜色的两种基色相对应的RGB颜色通道中,投影图像在标靶上的变形图像则包含在另外的一个基色对应的颜色通道中,从而达到了空间复用的目的。
在对多个位置拍摄多幅标靶上的混色图像(每个位置只拍摄一幅图像)后,分离出与组成标靶图案的两基色相对应的两颜色通道信息,可知此即为标定板上固有图像的灰度图,利用此两图像组中的任一组可以提取出特征点的信息,从而标定出摄像机的内外参数。剩下的那个通道的图像信息即是投影图像在标靶上的变形图像,利用此图像组,由标定的摄像机内外参数,计算出拍摄的变形投影图像像素坐标系中一点(up,vp)在世界坐标系中的位置,可由下式求得:
式中Mc为摄像机的投影变换矩阵,(uc,vc)为点(up,vp)在摄像机拍摄图像对应特征点的像素坐标。在得到投影图像在标靶上的变形图像角点的世界坐标位置后,可由8式求解出投影仪的投影变换矩阵Mp,从而完成投影仪的标定。
最后,在完成投影仪的标定后,即可利用双目立体标定方法标定出系统相关参数,从而完成系统的标定。
实验证明,本方法具有标定速度快,所需数据少的优点,可以应用于主动 光三维测量系统的标定中。
附图说明
图1测量系统中各坐标关系的示意图
图2采用红蓝颜色搭配、特征点为棋盘格的二维彩色平面标靶图
图3采用白、品红颜色搭配、特征点为棋盘格的投影图像
图4拍摄得到的空间复用的混色图像
图5分离的红通道图像
图6分离的蓝通道图像
图7分离的绿通道图像
图8采用上述颜色搭配和特征点的系统标定流程图
图9三维点数据在拟合平面上的分布
图10测量用的半圆柱物体
图11测量物体的一幅变形条纹图
图12彩色标靶法标定后重建的物体三维面形
图13本专利方法和传统的相位方法恢复面形沿X轴的切线比较
具体实施方式
下面结合一个具体实例进行说明,本实例并不包括本专利的所有内容。
在本实例中,制作在平面标靶上的图像是红蓝相间的棋盘格,如图2所示,投影到标靶上的图像是白和品红相间的棋盘格图像,如图3所示,彩色摄像机获取投影图像到达标靶后的混色图像,如图4所示,利用两个棋盘格合理的颜色设计,可以在红蓝通道分离出如图5和6所示的、对应物理存在于标靶上的 棋盘格图像,在绿通道分离出如图7所示的、投影到标靶上的棋盘格图像;然后由标靶上物理存在的棋盘格图像标定摄像机,由投影到标靶上的棋盘格图像标定投影仪;最后利用双目立体视觉原理标定出摄像机和投影仪的相对位置关系,从而完成系统的标定。
下面结合图8,具体说明标定步骤:
第一步:标定数据的获取
投影仪投射出白和品红相间的棋盘格图像,红蓝相间的棋盘格特征标靶分别放置在测量范围内若干不同位置,用摄像机拍摄记录不同位置处的混色图像,每个位置只需要拍摄一幅图像。设利用摄像机获得的任一幅彩色图像为P,将P的RGB三通道信息分离,为三幅灰度图像PR,PG,PB,按照RGB彩色模型的原理并根据事先的颜色搭配设计,PR应为标定靶标中的红颜色信息,如图5所示,同理PB应为标定靶标中的蓝色信息如图6所示,PG信息对应于投影仪投影白品红相间的图像到标定靶标表面后形成的图像,也成为黑白相间的棋盘格图案,如图7所示。
第二步:摄像机的标定
PR、PB中的标定靶标就成为黑白相间的棋盘格图案,对应于靶标表面物理存在的棋盘格图像,由角点提取算法提取出棋盘格图案上的特征点求出投影变换矩阵,完成摄像机的内外参数标定。
第三步:借助摄像机实现投影仪的标定
对PG提取角点的像素坐标(uc,vc),由第二步所求的摄像机内外参数可以计算出对应像素点的世界坐标系位置(Xw,Yw,Zw),接着提取相对应点在投影前图像上的像素坐标位置(up,vp)(即投影仪系统的像素坐标系位置),由(Xw,Yw,Zw)和(up,vp)求出投影变换矩阵,就可以标定出投影仪的内外参数。
第四步:由(uc,vc)、(Xw,Yw,Zw)、(up,vp)按双目立体原理标定出摄像机和投影仪的相对位置关系,从而完成系统的标定。
测量中使用的摄像机为佳能E503单反相机,使用的相机分辨率为1728×1152pixel。投影仪为三洋液晶投影仪,分辨率为1024×768pixel。标定出的摄像机和投影仪的内外参数如表1、表2和公式9、公式10所示。
表1标定出的摄像机内参数
表2标定出的投影仪内参数
系统旋转矩阵如下:
在整个系统标定结束后,为了验证标定方法的正确性,由投影仪投射一幅有计算机生成的有若干标记点的图像到白色平面上,再由摄像机拍摄获取图像。分别提取拍摄图像和投影图像相应标记点的像素坐标位置,由于系统已经标定,利用三角法计算得出这些点的三维坐标,理论上其应该位于同一个平面内。利用这些点由最小二乘法拟合出一个平面,然后,分别计算这些点到这个平面的距 离求出最大距离值。实验中提取并计算了156个标记点,点到拟合平面的最大距离为0.8039mm。图9为用文中方法标定系统后156个点在拟合平面两侧的分布情况。
同时利用标定后的系统测量了一个半圆柱物体的三维面形。传统的相位标定方法重建的半圆柱物体面形对应的棱和直径测量值为65.8798mm和65.4848mm,,本专利方法标定后重建物体面形的棱和直径测量值为65.4785mm和65.1873mm。图13是用本专利方法和相位方法标定的同一系统后,测量同一物体的重建面形在沿X轴方向的切线形状比较。
从测量结果可以看出分别使用本文方法和传统方法标定的同一系统测量的同一物体面形数据符合的较好。与现有方法相比,本发明的优点在于标定过程简单快捷,标定中每个平面位置只需要一幅图像,所需要处理的数据少。
Claims (4)
1.一种用于主动投影照明三维形貌测量系统中投影仪和摄像机同时快速标定的方法,其特征在于:在彩色图像的RGB模型基础上,利用合理的颜色设计,由测量系统中的投影仪将一幅具有彩色特征点的投影图像投影到一个用彩色特征点制作的二维平面标靶上,标靶上的图案称为标靶图案,由三基色中任意两种颜色组成,其一为背景色,另一为特征点的颜色,投影图像中白色背景上的特征点颜色应该为三基色中标靶图案不包含的那种颜色的补色,彩色摄像机一次同时获取平面标靶表面两个特征颜色图案的空间混色图像,再按照RGB颜色模型分离此图像,分离出与组成标靶图案的两基色相对应的两颜色通道信息,从而标定出摄像机的内外参数,剩下的那个通道的图像信息即是投影图像在标靶上的变形图像,用提取出来的投影图像标定投影仪,最后利用双目立体视觉原理标定出摄像机和投影仪的相对位置关系,从而完成整个主动投影照明三维形貌测量系统的标定,利用颜料和光的混色原理,使RGB模型的三通道中分别包含标靶图案信息和投影图像的彩色特征点颜色信息,从而达到空间复用的目的。
2.根据权利要求1所述的一种用于主动投影照明三维形貌测量系统中投影仪和摄像机同时快速标定的方法,其特征在于通过颜色设计,利用RGB模型,将标定两个设备所需图像信息混色进行空间复用,从而在每一个标定平面位置只需要拍摄一幅图像,即可完成投影仪的内部参数标定和投影仪与摄像机的系统标定,与传统的基于相位方法的标定方法相比,具有测量数据少,标定速度快的优点。
3.根据权利要求1所述的一种用于主动投影照明三维形貌测量系统中投影仪和摄像机同时快速标定的方法,其特征在于采用具有彩色特征点的投影图像和标靶图案,若标靶图案采用红蓝颜色搭配图像,则投影图像为白和品红颜色搭配图像;若标靶图案采用红绿颜色搭配图像,则投影图像为白黄颜色搭配图像;若标靶图案采用绿蓝颜色搭配图像,则投影图像为白青颜色搭配图像。
4.根据权利要求1所述的一种用于主动投影照明三维形貌测量系统中投影仪和摄像机同时快速标定的方法,其特征在于在图像分离时,若所采用的标靶图案本身为红蓝图像、投影图像本身为白品红图像,则拍摄的彩色图像中红、蓝通道分离的图像用来标定摄像机,绿通道图像可以用来标定投影仪;若所采用的标靶图案本身为红绿图像、投影图像本身为白黄图像,则拍摄的彩色图像中红、绿通道分离的图像用来标定摄像机,蓝通道图像可以用来标定投影仪;若所采用的标靶图案本身为绿蓝图像、投影图像本身为白青图像,则拍摄的彩色图像中绿、蓝通道分离的图像用来标定摄像机,红通道图像可以用来标定投影仪。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100592516A CN101576379B (zh) | 2009-05-12 | 2009-05-12 | 基于二维彩色标靶的主动投影三维测量系统快速标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100592516A CN101576379B (zh) | 2009-05-12 | 2009-05-12 | 基于二维彩色标靶的主动投影三维测量系统快速标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101576379A CN101576379A (zh) | 2009-11-11 |
CN101576379B true CN101576379B (zh) | 2011-04-27 |
Family
ID=41271374
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100592516A Expired - Fee Related CN101576379B (zh) | 2009-05-12 | 2009-05-12 | 基于二维彩色标靶的主动投影三维测量系统快速标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101576379B (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102201126B (zh) * | 2010-03-24 | 2013-02-13 | 联想(北京)有限公司 | 一种图像处理方法、系统及终端 |
US8736674B2 (en) | 2010-09-23 | 2014-05-27 | Dolby Laboratories Licensing Corporation | Method and system for 3D display calibration with feedback determined by a camera device |
CN102325411A (zh) * | 2011-09-01 | 2012-01-18 | 中国科学院深圳先进技术研究院 | 智能台灯及其进行互动的方法 |
CN102521822B (zh) * | 2011-10-25 | 2013-11-06 | 南京大学 | 一种用于机器视觉自动标定的主动发光式标靶及其标定方法 |
CN103559710B (zh) * | 2013-11-05 | 2016-06-29 | 重庆安钻理科技股份有限公司 | 一种用于三维重建系统的标定方法 |
EP3141954B1 (en) * | 2014-05-06 | 2024-04-24 | Ningbo Sunny Opotech Co., Ltd. | Light-deflection three-dimensional imaging device and projection device, and application thereof |
CN110716381B (zh) * | 2014-05-06 | 2021-10-15 | 宁波舜宇光电信息有限公司 | 转光三维成像装置和投射装置及其应用 |
CN105987805B (zh) * | 2015-02-02 | 2019-06-21 | 宁波舜宇光电信息有限公司 | 一种投影结构光系统的检测方法 |
CN104551411B (zh) * | 2014-11-18 | 2016-06-01 | 南京大学 | 一种双目立体视觉引导下的激光振镜加工系统的标定方法 |
CN106840251B (zh) * | 2015-12-07 | 2020-04-14 | 中国电力科学研究院 | 一种用于低压电流互感器外观检测的三维扫描系统 |
CN106296705A (zh) * | 2016-08-17 | 2017-01-04 | 上海交通大学 | 面向无人机成像色度的快速自动标定装置与方法 |
CN106651961B (zh) * | 2016-12-09 | 2019-10-11 | 中山大学 | 一种基于彩色立体标定物的无人机标定方法及系统 |
CN107082049A (zh) * | 2017-01-16 | 2017-08-22 | 王忠亮 | 车辆车窗控制装置 |
CN106657974B (zh) * | 2017-02-27 | 2024-02-09 | 北京图森智途科技有限公司 | 双目摄像机的控制方法及装置、双目摄像机 |
CN107063129B (zh) * | 2017-05-25 | 2019-06-07 | 西安知象光电科技有限公司 | 一种阵列式并行激光投影三维扫描方法 |
CN107478172B (zh) * | 2017-06-20 | 2019-08-16 | 南京航空航天大学 | 基于双目视觉的激光三维曲线轮廓定位投影方法 |
CN110298878B (zh) * | 2018-03-21 | 2021-03-12 | 北京猎户星空科技有限公司 | 一种目标物体三维位姿的确定方法、装置及电子设备 |
CN110823515B (zh) * | 2018-08-14 | 2022-02-01 | 宁波舜宇光电信息有限公司 | 一种结构光投射模组多工位检测装置及其检测方法 |
CN109360234B (zh) * | 2018-08-29 | 2020-08-21 | 清华大学 | 基于总体不确定度的三维场景重建方法及系统 |
CN109373912B (zh) * | 2018-12-21 | 2021-03-02 | 福州大学 | 一种基于双目视觉的非接触式六自由度位移测量方法 |
CN109900223B (zh) * | 2019-04-18 | 2021-10-08 | 盎锐(上海)信息科技有限公司 | 用于投影光栅建模的成像方法及装置 |
CN110360950B (zh) * | 2019-06-04 | 2021-07-06 | 重庆大学 | 五步相序黑白数字光栅编码方法 |
CN112179292B (zh) * | 2020-11-20 | 2022-07-08 | 苏州睿牛机器人技术有限公司 | 一种基于投影仪的线结构光视觉传感器标定方法 |
-
2009
- 2009-05-12 CN CN2009100592516A patent/CN101576379B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101576379A (zh) | 2009-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101576379B (zh) | 基于二维彩色标靶的主动投影三维测量系统快速标定方法 | |
CN202074952U (zh) | 基于单相机-单投影仪的三维形貌和彩色纹理获取系统 | |
CN107063129B (zh) | 一种阵列式并行激光投影三维扫描方法 | |
CN107003109B (zh) | 校准装置、校准方法、光学装置、摄影装置、投影装置、测量系统以及测量方法 | |
CN101308012B (zh) | 双单目白光三维测量系统标定方法 | |
CN100573040C (zh) | 物体表面三维轮廓结构光测量系统的标定方法 | |
US6125197A (en) | Method and apparatus for the processing of stereoscopic electronic images into three-dimensional computer models of real-life objects | |
CN104596439A (zh) | 一种基于相位信息辅助的散斑匹配三维测量方法 | |
CN110207614B (zh) | 一种基于双远心相机匹配的高分辨高精度测量系统及方法 | |
CN102445165B (zh) | 基于单幅彩色编码光栅的立体视觉测量方法 | |
CN109186491A (zh) | 基于单应性矩阵的平行多线激光测量系统及测量方法 | |
CN106091983B (zh) | 包含扫描方向信息的线结构光视觉测量系统完整标定方法 | |
CN109767476A (zh) | 一种自动对焦双目摄像头标定及深度计算方法 | |
US6195455B1 (en) | Imaging device orientation information through analysis of test images | |
CN104376558A (zh) | 一种基于长方体的Kinect深度相机的内参标定方法 | |
WO2000000926A1 (en) | Method and apparatus for capturing stereoscopic images using image sensors | |
CN102376089A (zh) | 一种标靶校正方法及系统 | |
CN106767528A (zh) | 一种基于彩色圆环标定板的光栅三维测量系统的标定方法 | |
CN103968782A (zh) | 一种基于彩色正弦结构光编码的实时三维测量方法 | |
CN104778694A (zh) | 一种面向多投影拼接显示的参数化自动几何校正方法 | |
CN103530880A (zh) | 基于投影高斯网格图案的摄像机标定方法 | |
CN104596444B (zh) | 一种基于编码图案投影的三维照相系统及方法 | |
CN112489109B (zh) | 一种三维成像系统方法、装置及三维成像系统 | |
CN104501741B (zh) | 一种用于三维面形测量的正交光栅相移方法 | |
CN104634323B (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 | ||
C14 | Grant of patent or utility model | ||
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: 20110427 Termination date: 20180512 |