CN106056620B - 线激光相机测量***标定方法 - Google Patents

线激光相机测量***标定方法 Download PDF

Info

Publication number
CN106056620B
CN106056620B CN201610420487.8A CN201610420487A CN106056620B CN 106056620 B CN106056620 B CN 106056620B CN 201610420487 A CN201610420487 A CN 201610420487A CN 106056620 B CN106056620 B CN 106056620B
Authority
CN
China
Prior art keywords
laser
camera
calibration
calibration region
equation
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
Application number
CN201610420487.8A
Other languages
English (en)
Other versions
CN106056620A (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.)
Southwest University of Science and Technology
Original Assignee
Southwest University of Science and Technology
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 Southwest University of Science and Technology filed Critical Southwest University of Science and Technology
Priority to CN201610420487.8A priority Critical patent/CN106056620B/zh
Publication of CN106056620A publication Critical patent/CN106056620A/zh
Application granted granted Critical
Publication of CN106056620B publication Critical patent/CN106056620B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30244—Camera pose

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供一种用于线激光位置标定的标定板及线激光相机测量***标定方法，该用于线激光位置标定的标定板包括第一标定区域和第二标定区域，所述第一标定区域的高度大于第二标定区域的高度，该线激光相机测量***标定方法包括以下步骤：S1、通过相机获取所述标定板的图像，在所述第二标定区域选取较准点，通过校准点求得相机的畸变系数、相机坐标系到世界坐标系的旋转平移矩阵及相机焦距；S2、打开激光器，激光器在所述第一标定区域和第二标定区域各投影出一条激光线，计算每条所述激光线上的任意两个点的坐标，通过所述坐标求得对应激光线的方程；S3、通过所述激光线的方程计算激光平面的方程以及激光平面与标定板平面的夹角。

Description

线激光相机测量***标定方法
技术领域
本发明涉及三维测量技术领域，具体涉及一种用于线激光位置标定的标定板及线激光相机测量***标定方法。
背景技术
三维重建是视觉领域中非常重要的一部分，摄像机加线激光的表面三维重建就是三维重建众多方法中的一种，其结构相对双目重建复杂一点，但是其精度比较高，常用于工业产品中各种缺陷检测或高精度表面重建环境。
在使用单目相机加线激光进行三维重建过程中，相机内外参数的标定(畸变系数、焦距、旋转平移矩阵等的求取)和线激光位置(光平面在相机坐标系或世界坐标系中的方程)的标定极其重要，这些参数直接影响到重建出来的效果好与否，精度高与低，线激光位置标定的方法主要分为三维标定和二维标定。在三维靶标标定中，需要制作高精度的三维靶标，立体靶标的大小需要与测量视场基本相当，存在加工难度大、费用高等缺点，未能在生产实际中广泛应用。在二维平面标定中，二维平面标定板的制作虽然简单，但其缺点在于需要两次摆放标定板，两次求取旋转平移矩阵，标定过程过于繁琐，计算量大。
发明内容
本发明的目的是提供一种线激光相机测量***标定方法，以解决现有夹具功能单一的问题。
为解决上述技术问题，本发明提供一种用于线激光位置标定的标定板，包括第一标定区域和第二标定区域，所述第一标定区域的高度大于第二标定区域的高度。
提供一种用于线激光位置标定的标定板，包括第一标定区域和第二标定区域，所述第一标定区域的高度大于第二标定区域的高度。
提供一种采用上述标定板的线激光相机测量***标定方法，包括以下步骤：
S1、通过相机获取所述标定板的图像；
S2、在所述第二标定区域选取较准点，通过校准点求得相机的畸变系数、相机坐标系到世界坐标系的旋转平移矩阵及相机焦距；
S3、打开激光器，激光器在所述第一标定区域和第二标定区域各投影出一条激光线；
S4、计算所述每条激光线上的任意两个点的坐标；
S5、通过所述坐标求得对应激光线的方程；
S6、通过所述激光线的方程计算激光平面的方程以及激光平面与标定板平面的夹角。
本发明的有益效果为：综合了三维和二维标定的特点，提出了一种基于梯形标定板的标定方法，可用于工业现场标定，降低了标定过程和运算的复杂度，提高了速度、稳定性及实用性。
附图说明
此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，在这些附图中使用相同的参考标号来表示相同或相似的部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1示出了本申请标定板的示意图；
图2示出了本申请的流程示意图；
图3示出了本申请相机模型的透视投影原理图；
图4示出了本申请世界坐标系到相机坐标系的旋转平移示意图；
图5示出了本申请提取较准点的示意图；
图6示出了本申请线激光标定时的透视原理图。
具体实施方式
为了使本发明的设计方法、技巧更加清楚明白，以下将结合附图以及实施例对本发明作详细的介绍。这里所描述的具体是实施方案仅用于解释本发明，并不用于限定本发明。
根据本申请的一个实施例，如图1所示，该用于线激光位置标定的标定板包括第一标定区域7和第二标定区域8，所述第一标定区域7的高度大于第二标定区域8的高度，其高度差为h；该用于线激光位置标定的第一标定区域7和第二标定区域8上均分布有大小相同的棋盘格,且棋盘格的边长已知。
根据本申请的一个实施例，如图2所示，该线激光相机测量***标定方法，包括以下步骤：
第一步，畸变系数k的求取：
将标定板放置在相机视野中,通过相机获取一张清晰的标定板图像,选取标定板上第二标定区域8内同一直线上的四点作为较准点来标定相机的畸变系数k，如图5所示，选取了对角线上的A，B，C，D四点来标定相机的畸变系数k。
如图3所示的透视投影原理图，其中Owxwywzw为世界坐标系，也就是现实中的坐标系，在标定过程中，将标定板第二标定区域8的一个顶点作为世界坐标系的原点，如图5所示，xw轴和yw轴分别与棋盘格边缘重合，zw轴垂直与标定板面。Oxyz为相机坐标系，其中Oz轴为相机光轴。O′XY为图像物理坐标系，坐标系建立在CCD图像阵列上，以mm为单位，其中O′为光轴与图像的交点。O″uv为图像像素坐标系，也是建立在CCD图像阵列上的，以像素为单位，其中O″为图像左上角起始点。相机坐标系、图像物理坐标系以及图像像素坐标系的X轴和Y轴分别平行。O′在图像像素坐标系中的坐标为(u0，v0)，OO′的长度等于相机的焦距f。设一点P，其世界坐标系为P(xw，yx，zw)相机坐标系为P(x，y，z)，P点在图像上的理想投影点所对应的物理坐标应为Pu(Xu，Yu)，但由于镜头存在畸变，所以其实际点坐标应该为畸变后的一点Pd(Xd，Yd)。
图像像素坐标系和图像物理坐标系存在如下关系：
其中，Dx和Dy分别是图像中每个像素的实际长和宽的尺寸，(u0，v0)为图像中心坐标，即光轴与图像的交点，也是图像物理坐标系的坐标原点。
由于相机镜头存在一定的畸变，设畸变系数为k，可得到等式：
其中，，X，Y为理想条件下的图像物理坐标，即不考虑畸变情况下的坐标；Xd，Yd为实际获取的图像物理坐标，即畸变后的坐标。
如图5所示，由于A，B，C，D四点都在标定板上，且棋盘格大小边长都是已知的，所以A，B，C，D四点在世界坐标系中的点坐标是已知的，且其zw坐标都为0。设 A，B，C，D在图像中的理想投影点(不考虑畸变)为A1，B1，C1，D1，根据透视投影原理和交比不变性可得到如下等式：
等式(3)是在理想情况下建立的方程，而A，B，C，D在图像中的实际投影点并不是A1，B1，C1，D1，而是其畸变后的点Ad，Bd，Cd，Dd。由于A，B，C，D四点的坐标都是已知的，且从等式(2)可知道A1，B1，C1，D1的图像物理坐标可以用Ad，Bd，Cd，Dd的图像物理坐标以及参数k的来表示，Ad，Bd，Cd，Dd的图像物理坐标可从其图像像素坐标变换中求取出来，而Ad，Bd，Cd，Dd图像像素坐标可从图像中获得，所以可以综合等式 (3)、等式(2)以及等式(1)求得畸变系数k。
第二步，相机外部参数以及相机焦距f的求取：
其中相机外部参数是世界坐标系到相机坐标系的旋转矩阵R以及平移矩阵 T。
如图3所示，根据投影变换原理，可得到如下方程：
写成矩阵形式为：
其中ρ为一比例因子，f为相机焦距，X，Y为图像物理坐标系中点的坐标，x，y，z 为相机坐标系中点的坐标。
两个坐标系之间可以通过一个坐标系平移和旋转达到和另一个坐标系相重合。如图4所示，为世界坐标系到相机坐标系的旋转平移示意图。
等式(6)中R为3×3正交单位矩阵，代表世界坐标系与相机坐标系之间的旋转角度关系，T为3×1的平移向量，其中
T＝(tx，ty，tz)T (7)
在等式(8)中，α为绕xw轴的旋转角，β为绕yw轴的旋转角，γ为绕zw轴的旋转角。各个角度旋转的正方向定义为从坐标系原点沿各坐标轴正方向观察时的逆时针旋转方向。如图4所示，将世界坐标系沿xw轴平移tx，沿yw轴平移ty，沿zw轴平移tz，然后绕xw轴旋转α，绕yw轴旋转β，绕zw轴旋转γ，即可与相机坐标系重合。
令
将等式(9)和等式(7)代入等式(6)可以得到：
将等式(10)和(1)代入(5)，可以得到：
在等式(11)中，r1～r9，tx，ty，tz，f为未知参数，也就是需要标定的部分参数。
根据等式(10)和等式(4)可以获得以下等式：
由于世界坐标系是建立在标定板上的，而A，B，C，D四个点都在标定板上，所以其zw＝0，继而等式(12)可改为：
任意选取标定板上六个已知较准点代入等式(13)中可以求出 r1，r2，r4，r5，tx，ty，由于旋转矩阵为正交单位阵，根据单位正交矩阵的特性，可以通过r1，r2，r4，r5求得r3，r6，r7，r8，r9。
从等式(10)和等式(4)还可以获得以下等式：
结合等式(14)和等式(1)，f，tz可以很容易的被求解出来。到此，相机的参数标定基本完成。
第三步，激光线投影方程的计算：
标定板和相机保持刚才的相对位置不动，打开激光器，激光器将在标定板上投影出两条激光线，如图6所示，图中标注1为相机平面，2表示激光平面， 3为激光平面与标定板第一标定区域7的交线L2，4为标定板第一标定区域7，5 为激光平面与标定板第二标定区域8的交线L1，6是标定板的第二标定区域8，h 为第一标定区域7和第二标定区域8已知的高度差，θ为激光平面与标定板平面的夹角。激光平面交标定板平面于L1，L2两条直线，其中L1位于标定板第二标定区域8，L2位于标定板第一标定区域7，同时L1，L2均位于相机的某一个平面上。
设L1上的一点P1，可知道它在图像像素坐标系中的坐标为P1(u1，v1)，可以通过等式(11)求得其在世界坐标系中的坐标为P1(Px，Py，0），同理，可求得L1上另一点 Q1的世界坐标Q1(Qx，Qy，0)，通过P1，Q1的点坐标，可求得激光线L1在世界坐标系中的方程。同理，可以在激光线L2上取不同的两点P2，Q2，由于标定板高度差为已知 h，再结合等式(11)可求得P2，Q2两点在世界坐标系中的坐标P2(Px，Py，h)和 Q2(Qx，Qy，h)，继而求得激光线L2在世界坐标系中的方程。
第四步，激光平面方程及其与标定板平面夹角的计算：
由于L1，L2都属于同一个激光平面，且L1，L2相互平行，所以可通过L1，L2直线的方程，求得激光平面的方程。由于世界坐标系是建立在标定板上的，所以标定板平面为Owxwyw平面，标定板平面方程和激光平面方程均已知，通过数学知识很容易求得两已知平面的夹角θ，即激光平面与标定板平面之间的夹角θ。
结合相机参数f，k以及旋转平移矩阵和激光平面的方程，可求得激光平面上任意一点在世界坐标系中的坐标，从而可求得激光平面与被测物体表面交线上点的在世界坐标系中的三维坐标，通过被测物体和激光线结构光的相对运动来扫描物体轮廓，从而重建出被测物体的三维表面轮廓。

Claims (2)

1.一种线激光相机测量***标定方法，用于测量线激光位置标定的标定板包括第一标定区域和第二标定区域，第一标定区域的高度大于第二标定区域的高度，且第一标定区域和第二标定区域上均分布有大小相同的棋盘格；其特征在于，包括以下步骤：
S1、通过相机获取所述标定板的图像；
S2、在所述第二标定区域选取校 准点，通过校准点求得相机的畸变系数、相机坐标系到世界坐标系的旋转平移矩阵及相机焦距；
S3、打开激光器，激光器在所述第一标定区域和第二标定区域各投影出一条激光线；
S4、计算每条所述激光线上的任意两个点的坐标；
S5、通过所述坐标求得对应激光线的方程；
S6、通过所述激光线的方程计算激光平面的方程以及激光平面与标定板平面的夹角。
2.根据权利要求1所述的线激光相机测量***标定方法，其特征在于，计算所述畸变系数的公式为：
其中，q2＝Xd 2+Yd 2，X、Y为不考虑畸变情况下的图像物理坐标，Xd、Yd为实际获取的图像物理坐标，即畸变后的坐标，k为畸变系数。
CN201610420487.8A 2016-06-13 2016-06-13 线激光相机测量***标定方法 Expired - Fee Related CN106056620B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610420487.8A CN106056620B (zh) 2016-06-13 2016-06-13 线激光相机测量***标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610420487.8A CN106056620B (zh) 2016-06-13 2016-06-13 线激光相机测量***标定方法

Publications (2)

Publication Number Publication Date
CN106056620A CN106056620A (zh) 2016-10-26
CN106056620B true CN106056620B (zh) 2018-10-19

Family

ID=57168963

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610420487.8A Expired - Fee Related CN106056620B (zh) 2016-06-13 2016-06-13 线激光相机测量***标定方法

Country Status (1)

Country Link
CN (1) CN106056620B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3615955A4 (en) * 2017-04-28 2020-05-13 SZ DJI Technology Co., Ltd. CALIBRATION OF LASER AND VISION SENSORS
CN107911573B (zh) * 2017-08-31 2019-05-03 大连成者科技有限公司 扫描图像的无标定展平方法
CN109029292A (zh) * 2018-08-21 2018-12-18 孙傲 一种容器内表面三维形貌无损检测装置及检测方法
CN109781164B (zh) * 2018-12-28 2021-02-05 长沙长泰机器人有限公司 一种线激光传感器的静态标定方法
CN112292611A (zh) * 2019-05-24 2021-01-29 深圳市速腾聚创科技有限公司 坐标校正方法、装置、计算设备及计算机存储介质
CN111272102A (zh) * 2020-05-06 2020-06-12 中国空气动力研究与发展中心低速空气动力研究所 一种线激光扫描三维测量标定方法
CN111739104B (zh) * 2020-06-24 2024-05-03 深圳市道通科技股份有限公司 一种激光标定***的标定方法、装置以及激光标定***
CN111623727B (zh) * 2020-07-17 2021-08-06 泉州深索思传感器科技有限公司 一种3d线激光轮廓传感器的快速标定方法
CN111830049B (zh) * 2020-07-24 2023-09-22 深圳市道通科技股份有限公司 一种刹车盘磨损区的识别方法及磨损识别***
CN112179292B (zh) * 2020-11-20 2022-07-08 苏州睿牛机器人技术有限公司 一种基于投影仪的线结构光视觉传感器标定方法
CN112620926B (zh) * 2020-12-18 2022-10-21 深圳辰视智能科技有限公司 一种焊点追踪方法、装置及存储介质
CN113655466B (zh) * 2021-07-16 2023-09-29 国家石油天然气管网集团有限公司 一种基于结构光三角测量的标定***和标定方法
CN114373020B (zh) * 2022-01-07 2023-04-28 天津大学 一种基于联合外参估计的线结构光传感器光平面标定方法
CN114565681B (zh) * 2022-03-01 2022-11-22 禾多科技（北京）有限公司 一种相机标定方法、装置、设备、介质及产品

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3799858B2 (ja) * 1999-02-12 2006-07-19 スズキ株式会社 形状認識装置
CN102628671A (zh) * 2012-03-29 2012-08-08 中国人民解放军第二炮兵工程学院 一种基于单摄像机双色线结构光的三维坐标测量方法
CN204461385U (zh) * 2015-02-11 2015-07-08 苏州天准科技股份有限公司 用于测量传感器间数据进行融合标定的标定板
CN105091782A (zh) * 2015-05-29 2015-11-25 南京邮电大学 一种基于双目视觉的多线激光器光平面标定方法
CN105486289A (zh) * 2016-01-31 2016-04-13 山东科技大学 一种激光摄影测量***及相机标定方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3799858B2 (ja) * 1999-02-12 2006-07-19 スズキ株式会社 形状認識装置
CN102628671A (zh) * 2012-03-29 2012-08-08 中国人民解放军第二炮兵工程学院 一种基于单摄像机双色线结构光的三维坐标测量方法
CN204461385U (zh) * 2015-02-11 2015-07-08 苏州天准科技股份有限公司 用于测量传感器间数据进行融合标定的标定板
CN105091782A (zh) * 2015-05-29 2015-11-25 南京邮电大学 一种基于双目视觉的多线激光器光平面标定方法
CN105486289A (zh) * 2016-01-31 2016-04-13 山东科技大学 一种激光摄影测量***及相机标定方法

Also Published As

Publication number Publication date
CN106056620A (zh) 2016-10-26

Similar Documents

Publication Publication Date Title
CN106056620B (zh) 线激光相机测量***标定方法
CN107063129B (zh) 一种阵列式并行激光投影三维扫描方法
CN111750806B (zh) 一种多视角三维测量***及方法
CN110288642B (zh) 基于相机阵列的三维物体快速重建方法
CN106595528B (zh) 一种基于数字散斑的远心显微双目立体视觉测量方法
CN107170043B (zh) 一种三维重建方法
CN106408556B (zh) 一种基于一般成像模型的微小物体测量***标定方法
CN107121109A (zh) 一种基于前镀膜平面镜的结构光参数标定装置及方法
CN106127745B (zh) 结构光3d视觉***与线阵相机的联合标定方法及装置
CN106091984B (zh) 一种基于线激光的三维点云数据获取方法
CN108122191A (zh) 鱼眼图像拼接成全景图像和全景视频的方法及装置
CN102243432A (zh) 全景立体摄像装置
CN109859272A (zh) 一种自动对焦双目摄像头标定方法及装置
CN110191326A (zh) 一种投影***分辨率扩展方法、装置和投影***
CN103115613A (zh) 一种空间三维定位方法
CN102831601A (zh) 基于联合相似性测度和自适应支持权重的立体匹配方法
Zhang et al. Three-dimensional shape measurement using a structured light system with dual cameras
CN111009030A (zh) 一种多视高分辨率纹理图像与双目三维点云映射方法
CN111091599A (zh) 一种基于球体标定物的多相机-投影仪***标定方法
Zou et al. High-accuracy calibration of line-structured light vision sensors using a plane mirror
CN107560554A (zh) 一种基于旋转透镜的三维信息视觉测量方法
CN106643563A (zh) 一种桌面式大视场三维扫描装置及方法
JP2002516443A (ja) ３次元表示のための方法および装置
CN104200476B (zh) 利用双平面镜装置中的圆周运动求解摄像机内参数的方法
CN108734738A (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
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: 20181019

Termination date: 20190613