用于线激光位置标定的标定板及线激光相机测量系统标定方法
【专利摘要】本发明提供一种用于线激光位置标定的标定板及线激光相机测量系统标定方法，该用于线激光位置标定的标定板包括第一标定区域和第二标定区域，所述第一标定区域的高度大于第二标定区域的高度，该线激光相机测量系统标定方法包括以下步骤：S1、通过相机获取所述标定板的图像，在所述第二标定区域选取较准点，通过校准点求得相机的畸变系数、相机坐标系到世界坐标系的旋转平移矩阵及相机焦距；S2、打开激光器，激光器在所述第一标定区域和第二标定区域各投影出一条激光线，计算每条所述激光线上的任意两个点的坐标，通过所述坐标求得对应激光线的方程；S3、通过所述激光线的方程计算激光平面的方程以及激光平面与标定板平面的夹角。
【专利说明】
用于线激光位置标定的标定板及线激光相机测量系统标定 方法
技术领域
[0001] 本发明涉及三维测量技术领域，具体涉及一种用于线激光位置标定的标定板及线 激光相机测量系统标定方法。
【背景技术】
[0002] 三维重建是视觉领域中非常重要的一部分，摄像机加线激光的表面三维重建就是 三维重建众多方法中的一种，其结构相对双目重建复杂一点，但是其精度比较高，常用于工 业产品中各种缺陷检测或高精度表面重建环境。
[0003] 在使用单目相机加线激光进行三维重建过程中，相机内外参数的标定(畸变系数、 焦距、旋转平移矩阵等的求取）和线激光位置(光平面在相机坐标系或世界坐标系中的方 程)的标定极其重要，这些参数直接影响到重建出来的效果好与否，精度高与低，线激光位 置标定的方法主要分为三维标定和二维标定。在三维祀标标定中，需要制作高精度的三维 靶标，立体靶标的大小需要与测量视场基本相当，存在加工难度大、费用高等缺点，未能在 生产实际中广泛应用。在二维平面标定中，二维平面标定板的制作虽然简单，但其缺点在于 需要两次摆放标定板，两次求取旋转平移矩阵，标定过程过于繁琐，计算量大。
【发明内容】
[0004] 本发明的目的是提供一种线激光相机测量系统标定方法，以解决现有夹具功能单 一的问题。
[0005] 为解决上述技术问题，本发明提供一种用于线激光位置标定的标定板，包括第一 标定区域和第二标定区域，所述第一标定区域的高度大于第二标定区域的高度。
[0006] 提供一种用于线激光位置标定的标定板，包括第一标定区域和第二标定区域，所 述第一标定区域的高度大于第二标定区域的高度。
[0007] 提供一种采用上述标定板的线激光相机测量系统标定方法，包括以下步骤：
[0008] Sl、通过相机获取所述标定板的图像；
[0009] S2、在所述第二标定区域选取较准点，通过校准点求得相机的畸变系数、相机坐标 系到世界坐标系的旋转平移矩阵及相机焦距；
[0010] S3、打开激光器，激光器在所述第一标定区域和第二标定区域各投影出一条激光 线；
[0011] S4、计算所述每条激光线上的任意两个点的坐标；
[0012] S5、通过所述坐标求得对应激光线的方程；
[0013] S6、通过所述激光线的方程计算激光平面的方程以及激光平面与标定板平面的夹 角。
[0014] 本发明的有益效果为:综合了三维和二维标定的特点，提出了一种基于梯形标定 板的标定方法，可用于工业现场标定，降低了标定过程和运算的复杂度，提高了速度、稳定 性及实用性。
【附图说明】
[0015] 此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，在这 些附图中使用相同的参考标号来表示相同或相似的部分，本申请的示意性实施例及其说明 用于解释本申请，并不构成对本申请的不当限定。在附图中：
[0016] 图1示出了本申请标定板的示意图；
[0017]图2示出了本申请的流程示意图；
[0018] 图3示出了本申请相机模型的透视投影原理图；
[0019] 图4示出了本申请世界坐标系到相机坐标系的旋转平移示意图；
[0020] 图5示出了本申请提取较准点的示意图；
[0021 ]图6示出了本申请线激光标定时的透视原理图。
【具体实施方式】
[0022] 为了使本发明的设计方法、技巧更加清楚明白，以下将结合附图以及实施例对本 发明作详细的介绍。这里所描述的具体是实施方案仅用于解释本发明，并不用于限定本发 明。
[0023] 根据本申请的一个实施例，如图1所示，该用于线激光位置标定的标定板包括第一 标定区域7和第二标定区域8,所述第一标定区域7的高度大于第二标定区域8的高度，其高 度差为h;该用于线激光位置标定的第一标定区域7和第二标定区域8上均分布有大小相同 的棋盘格，且棋盘格的边长已知。
[0024] 根据本申请的一个实施例，如图2所示，该线激光相机测量系统标定方法，包括以 下步骤：
[0025] 第一步，畸变系数k的求取：
[0026] 将标定板放置在相机视野中，通过相机获取一张清晰的标定板图像，选取标定板 上第二标定区域8内同一直线上的四点作为较准点来标定相机的畸变系数k，如图5所示，选 取了对角线上的A，B，C，D四点来标定相机的畸变系数k。
[0027]如图3所示的透视投影原理图，其中0wXwywZw为世界坐标系，也就是现实中的坐标 系，在标定过程中，将标定板第二标定区域8的一个顶点作为世界坐标系的原点，如图5所 示，Xw轴和yw轴分别与棋盘格边缘重合，Z w轴垂直与标定板面。Oxyz为相机坐标系，其中Oz轴 为相机光轴。VXY为图像物理坐标系，坐标系建立在CCD图像阵列上，以mm为单位，其中为 光轴与图像的交点。〇〃uv为图像像素坐标系，也是建立在CCD图像阵列上的，以像素为单位， 其中〇〃为图像左上角起始点。相机坐标系、图像物理坐标系以及图像像素坐标系的X轴和Y 轴分别平行。在图像像素坐标系中的坐标为(11〇，仰），00 /的长度等于相机的焦距€。设一点 P，其世界坐标系为P (Xw，yx，Zw)相机坐标系为P (X，y，z)，P点在图像上的理想投影点所对应 的物理坐标应为Pu (Xu，Yu )，但由于镜头存在畸变，所以其实际点坐标应该为畸变后的一点 Pd(Xd，Yd)〇
[0028]图像像素坐标系和图像物理坐标系存在如下关系：
[0029]
(I)
[0030] 其中，Dx和Dy分别是图像中每个像素的实际长和宽的尺寸，（UQ，Vo)为图像中心坐 标，即光轴与图像的交点，也是图像物理坐标系的坐标原点。
[0031] 由于相机镜头存在一定的畸变，设畸变系数为k，可得到等式：
[0032]
(2)
[0033] 其中，#二_ + 1|，Χ，Υ为理想条件下的图像物理坐标，即不考虑畸变情况下的坐 标;Xd，Yd为实际获取的图像物理坐标，即畸变后的坐标。
[0034]如图5所示，由于A，B，C，D四点都在标定板上，且棋盘格大小边长都是已知的，所以 A，B，C，D四点在世界坐标系中的点坐标是已知的，且其Zw坐标都为0。设A，B，C，D在图像中的 理想投影点（不考虑畸变彡为六:^^心^^根据透视投影原理和交比不变性可得到如下等 式：
[0035]
[0036]等式（3)是在理想情况下建立的方程，而A，B，C，D在图像中的实际投影点并不是 △1，81，&，01，而是其畸变后的点六<1，8 (1心，0(1。由于4,8，(：，0四点的坐标都是已知的，且从等式 (2)可知道A1，B1，C1，图像物理坐标可以用Ad，Bd，Cd，Dd的图像物理坐标以及参数k的来表 示，Ad，Bd，Cd，Dd的图像物理坐标可从其图像像素坐标变换中求取出来，而Ad，Bd，Cd，Dd图像像 素坐标可从图像中获得，所以可以综合等式(3)、等式(2)以及等式(1)求得畸变系数k。
[0037] 第二步，相机外部参数以及相机焦距f的求取：
[0038] 其中相机外部参数是世界坐标系到相机坐标系的旋转矩阵R以及平移矩阵T。
[0039] 如图3所示，根据投影变换原理，可得到如下方程：
[0040]
[0041]
[0042]
[0043]其中P为一比例因子，f为相机焦距，X，Y为图像物理坐标系中点的坐标，x，y，z为相 机坐标系中点的坐标。
[0044] 两个坐标系之间可以通过一个坐标系平移和旋转达到和另一个坐标系相重合。如 图4所示，为世界坐标系到相机坐标系的旋转平移示意图。
[0045]
[0046] 等式(6)中R为3X3正交单位矩阵，代表世界坐标系与相机坐标系之间的旋转角度 关系，T为3 Xl的平移向量，其中
[0047] T=(tx，ty，tz)T (7)
[0048]
[0049] 在等式(8)中，α为绕Xw轴的旋转角，β为绕yw轴的旋转角，γ为绕Zw轴的旋转角。各 个角度旋转的正方向定义为从坐标系原点沿各坐标轴正方向观察时的逆时针旋转方向。如 图4所不，将世界坐标系沿Xw轴平移t x，沿轴平移ty，沿Zw轴平移tz，然后绕Xw轴旋转α，绕yw 轴旋转队绕2|轴旋转γ，即可与相机坐标系重合。
[0050] 八
[0051]
[0052] 将等式(9)和等式(7)代入等式(6)可以得到：
[0053]
[0054] 将等式（10)和（1)代入(5)，可以得到：
[0055；
[0056] 在等式(11)中，ri~r9，tx，ty，t z，f为未知参数，也就是需要标定的部分参数。
[0057] 根据等式(10)和等式(4)可以获得以下等式：
[0058]
[0059] 由于世界坐标系是建立在标定板上的，而A，B，C，D四个点都在标定板上，所以其~ =0,继而等式（12)可改为：
[0060]
[0061 ] 任意选取标定板上六个已知较准点代入等式（13)中可以求出ri，Γ2，Γ4，Γ5，tx，ty， 由于旋转矩阵为正交单位阵，根据单位正交矩阵的特性，可以通过ri，Γ2，Γ4，Γ5求得Γ3，Γ6， Π,YSjT^o
[0062] 从等式(10)和等式(4)还可以获得以下等式：
[0063]
[0064] 结合等式（14)和等式（l)，f，tz可以很容易的被求解出来。到此，相机的参数标定 基本完成。
[0065]第三步，激光线投影方程的计算：
[0066]标定板和相机保持刚才的相对位置不动，打开激光器，激光器将在标定板上投影 出两条激光线，如图6所示，图中标注1为相机平面，2表示激光平面，3为激光平面与标定板 第一标定区域7的交线L2,4为标定板第一标定区域7,5为激光平面与标定板第二标定区域8 的交线L 1,6是标定板的第二标定区域8，h为第一标定区域7和第二标定区域8已知的高度 差，Θ为激光平面与标定板平面的夹角。激光平面交标定板平面于L1，1^两条直线，其中LHi 于标定板第二标定区域8，L2位于标定板第一标定区域7,同时L1,1^均位于相机的某一个平 面上。
[0067] SL1上的一点丹，可知道它在图像像素坐标系中的坐标为P1 (m，V1)，可以通过等式 (11)求得其在世界坐标系中的坐标为P1(PnPpO)，同理，可求得1^上另一点&的世界坐标Q1 (Qx，Qy，0 )，通过P1，Q1的点坐标，可求得激光线1^在世界坐标系中的方程。同理，可以在激光 线Ud:取不同的两点P 2，Q2,由于标定板高度差为已知h，再结合等式（11)可求得P2，Q2两点在 世界坐标系中的坐标P 2(Px，Py，h)和Q2(Qx，Q y，h)，继而求得激光线1^2在世界坐标系中的方 程。
[0068]第四步，激光平面方程及其与标定板平面夹角的计算：
[0069]由于L1,1^都属于同一个激光平面，且L1,1^相互平行，所以可通过L 1,1^直线的方 程，求得激光平面的方程。由于世界坐标系是建立在标定板上的，所以标定板平面为〇wXwyw 平面，标定板平面方程和激光平面方程均已知，通过数学知识很容易求得两已知平面的夹 角θ，即激光平面与标定板平面之间的夹角θ。
[0070]结合相机参数f，k以及旋转平移矩阵和激光平面的方程，可求得激光平面上任意 一点在世界坐标系中的坐标，从而可求得激光平面与被测物体表面交线上点的在世界坐标 系中的三维坐标，通过被测物体和激光线结构光的相对运动来扫描物体轮廓，从而重建出 被测物体的三维表面轮廓。
【主权项】
1. 一种用于线激光位置标定的标定板，其特征在于，包括第一标定区域和第二标定区 域，所述第一标定区域的高度大于第二标定区域的高度。2. 根据权利要求1所述的用于线激光位置标定的标定板，其特征在于，所述第一标定区 域和第二标定区域上均分布有大小相同的棋盘格。3. -种线激光相机测量系统标定方法，其特征在于，包括以下步骤： 51、 通过相机获取所述标定板的图像； 52、 在所述第二标定区域选取较准点，通过校准点求得相机的畸变系数、相机坐标系到 世界坐标系的旋转平移矩阵及相机焦距； 53、 打开激光器，激光器在所述第一标定区域和第二标定区域各投影出一条激光线； 54、 计算所述每条激光线上的任意两个点的坐标； 55、 通过所述坐标求得对应激光线的方程； 56、 通过所述激光线的方程计算激光平面的方程以及激光平面与标定板平面的夹角。4. 根据权利要求3所述的线激光相机测量系统标定方法，其特征在于，计算所述畸变系 数的公式为：其中，q2 = Xd2+Yd2，X、Y为不考虑畸变情况下的图像物理坐标，Xd、Yd为实际获取的图像物 理坐标，即畸变后的坐标，k为畸变系数。
【文档编号】G06T7/00GK106056620SQ201610420487
【公开日】2016年10月26日
【申请日】2016年6月13日
【发明人】刘桂华, 张华 , 游斌相, 龙惠民
【申请人】西南科技大学