一种相机与多线激光的联合标定方法

1.本发明属于图像处理与多传感器标定技术领域，特别涉及一种相机与多线激光的联合标定方法。


背景技术：

2.在机器视觉中，相机与线激光的联合测量一直是一种解决方法，为了让结果测量结果更准确，通常需要对相机与激光线进行联合标定。
3.但是目前的各种标定方式，标定简单的精度不够，精度足够的又不够简单方便。于是本发明提出了一种全新的相机与多线激光的联合标定方法，兼顾了准确性与简单易实施特点。


技术实现要素：

4.为了克服以上技术问题，本发明的目的在于提供一种相机与多线激光的联合标定方法，可以快速且准确的标定出相机与多线激光的空间位置，从而使后续的测量结果更准确，具有操作方便，标定速度快，标定结果准确的特点。
5.为了实现上述目的，本发明采用的技术方案是：
6.一种相机与多线激光的联合标定方法，包括以下步骤；
7.步骤s1，相机的标定，图像的采集，加载与去畸变预处理；
8.步骤s2，通过标定板上的特征点与相机的内外参，拟合相机坐标系下的标定板平面；
9.步骤s3，提取亚像素级别的激光点，并对这些点进行聚类或者分类，确定该点归于哪条激光线。并拟合出各条激光线在像素坐标系下的方程ax+by+c＝0；
10.步骤s4，计算激光线与标定板上特征点在水平或者垂直方向上连线的交点；
11.步骤s5，根据s2与s5得到的已知数据，计算s4中的2d交点在相机坐标系下的3d坐标；
12.步骤s6，拟合各条2d激光线在相机坐标系下所属的3d激光平面；
13.步骤s7，通过激光平面与相机x轴延长线的交点，确定等效的激光发射点。
14.所述步骤s1具体为，首先对相机进行单独标定，获得相机的内参，对每一张图像进行rectify矫正，消除因为镜头带来的畸变，得到相机坐标系下的真实坐标；保证激光线完全打在标定板平面上；保证标定板中的所有目标点都显示在图像熵；同一位置,在激光器打开与激光器关闭情况下，分别采集两组图像；对采集到的图像进行矫正处理，去除畸变；
15.得到的矫正后的所有的激光线关闭情况下的各组i类图像；根据已知标定板中所有特征点的3d坐标，与得到的相机的内参，求解出各组图像中的相机的外参(相机的旋转r与平移t)，通过计算，获取相机坐标系下，各组i类图像中的标定板上的各个特征点的3d空间坐标；
16.公式如下：pc＝tw
‑
>c*pw，其中，pc为需要计算得到的相机坐标系下的3d坐标，tw
‑
>c为计算的到的各组照片的相机外参，pw为已知的当前组照片的标定板系统坐标系下的3d坐标。
17.所述步骤s2通过步骤s1得到各组i类图像中的的特征点，在相机坐标系下3d坐标，进行空间平面拟合，剔除误差较大的3d点后，最终得到各组照片中的标定板所属的平面坐标，公式为ax+by+cz+d＝0，其中x,y,z为变量，a,b,c,d为计算得到的常量。
18.所述步骤s3，获取步骤s1中得到的矫正后的所有的激光线打开情况下的各组ii类图像；
19.使用灰度重心法，提取得到各组ii类图像中的激光线中的亚像素中心点；
20.对得到的所有亚像素激光点进行聚类，得到图像中的每条激光线归属与哪条激光线；
21.各个激光线投射在标定板上都是一条直线，所以可以根据得到的各个激光线所属的2d激光点坐标，拟合出各组ii类图像中所属的各条激光线的2d坐标，坐标公式为ax+by+c＝0。其中x,y为变量，a,b,c为计算得到的常量，直线的斜率为
‑
a/b。
22.所述步骤s4，根据s3得到的斜率，计算各组i类图像中，各条激光线与标定板横/竖方向上特征点的连线与延长线上的交点；
23.其中，如果斜率>1，说明该组激光线更接近垂直，则计算各行特征点的直线公式ax+by+c＝0；
24.如果斜率<1，说明该组激光线更接近垂直，则计算各列特征点的直线公式ax+by+c＝0；
25.根据斜率计算各条激光线与标定板横/竖方向上特征点的连线与延长线上的交点；
26.其中，如果斜率>1，说明该组激光线更接近垂直，则计算与水平方向连线与延长线上的交点；
27.如果斜率<1，说明该组激光线更接近垂直，则计算与竖直方向连线与延长线上的交点。
28.所述步骤s5根据拟合的标定板平面，与得到的各组ii类图像中的各条激光线的交点，与张正友的小孔成像模型公式，即可得到最终结果，得到相机坐标系下的各条激光线的空间方程。
29.所述步骤s7，根据激光平面方程，与相机坐标的x轴坐标的交点，就是虚拟的激光发射点，激光入射角度可以根据需要进行定义；
30.循环执行步骤s5和步骤s6，遍历各个激光线，即可分别得到每条激光线相对于相机位置的激光平面方程。
31.本发明的有益效果。
32.本发明可以在一次性的拍摄多组照片后，直接并准确标定出各条激光线的空间平面坐标，对使用相机与多线激光器进行经过的空间测量有重大意义。
附图说明：
33.图1为本发明流程示意图。
具体实施方式
34.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
35.[1]特质标定板的描述如下：
[0036]
标定板更大，中心部分是传统的棋盘格标定板，但拥有更大的边缘空白区域。中心是棋盘格，边缘是白色图涂料。
[0037]
[2]步骤s1，首先对相机进行单独标定，获得相机的内参。
[0038]
[3]对每一张图像进行rectify矫正，消除因为镜头带来的畸变，得到相机坐标系下的真实坐标。
[0039]
[4]采集图像：
[0040]
保证激光线完全打在标定板平面上。
[0041]
保证标定板中的所有目标点都显示在图像熵。
[0042]
同一位置,在激光器打开与激光器关闭情况下，分别采集两组图像。
[0043]
[5]矫正：
[0044]
对采集到的图像进行矫正处理，去除畸变。
[0045]
[6]获取步骤[5]中得到的矫正后的所有的激光线关闭情况下的各组i类图像。
[0046]
此时，根据已知标定板中所有特征点的3d坐标，与步骤[2]得到的相机的内参，求解出各组图像中的相机的外参(相机的旋转r与平移t)。
[0047]
[7]通过计算，获取相机坐标系下，各组i类图像中的标定板上的各个特征点的3d空间坐标。
[0048]
公式如下：pc＝tw
‑
>c*pw，其中，pc为需要计算得到的相机坐标系下的3d坐标，tw
‑
>c为步骤[6]计算的到的各组照片的相机外参，pw为已知的当前组照片的标定板系统坐标系下的3d坐标。
[0049]
[8]步骤s2，通过步骤[7]得到各组i类图像中的的特征点，在相机坐标系下3d坐标，进行空间平面拟合，剔除误差较大的3d点后，最终得到各组照片中的标定板所属的平面坐标，公式为ax+by+cz+d＝0，其中x,y,z为变量，a,b,c,d为计算得到的常量。
[0050]
[9]步骤s3，获取步骤[5]中得到的矫正后的所有的激光线打开情况下的各组ii类图像。
[0051]
[10]使用灰度重心法，提取得到各组ii类图像中的激光线中的亚像素中心点。
[0052]
[11]对步骤[10]中得到的所有亚像素激光点进行聚类，得到图像中的每条激光线归属与哪条激光线。
[0053]
[12]各个激光线，投射在标定板上都是一条直线。所以可以根据步骤[11]得到的各个激光线所属的2d激光点坐标，拟合出各组ii类图像中所属的各条激光线的2d坐标，坐标公式为ax+by+c＝0。其中x,y为变量，a,b,c为计算得到的常量，直线的斜率为
‑
a/b。
[0054]
[13]步骤s4，根据步骤[12]得到的斜率，计算各组i类图像中，各条激光线与标定板横/竖方向上特征点的连线与延长线上的交点。
[0055]
其中，如果斜率>1，说明该组激光线更接近垂直，则计算各行特征点的直线公式ax
+by+c＝0；
[0056]
如果斜率<1，说明该组激光线更接近垂直，则计算各列特征点的直线公式ax+by+c＝0；
[0057]
[14]根据步骤[12]得到的斜率，计算各条激光线与标定板横/竖方向上特征点的连线与延长线上的交点。
[0058]
其中，如果斜率>1，说明该组激光线更接近垂直，则计算与水平方向连线与延长线上的交点；
[0059]
如果斜率<1，说明该组激光线更接近垂直，则计算与竖直方向连线与延长线上的交点；
[0060]
[15]步骤s5，根据步骤[8]拟合的标定板平面，与步骤[14]得到的各组ii类图像中的各条激光线的交点，与张正友的小孔成像模型公式，即可按照如下公式，即可得到最终结果，得到相机坐标系下的各条激光线的空间方程。
[0061]
[16]推倒公式如下图：
[0062]
[0063][0064]
其中，zc是z方向上的深度d；
[0065]
(xc,yc,zc)是相机坐标系下的3d空间坐标，也就似乎最终的结果；cameramatrix为3x4的相机内参矩阵；
[0066]
u,v为步骤[14]得到的特征点与各个激光线在像素坐标系下的交点，为输入的已知变量；
[0067]
fx,fy是相机的内参；
[0068]
u0,v0是相机坐标系下的中心点；
[0069]
a,b,c,d是步骤[8]拟合的标定板平面的计算参数。
[0070]
[17]通过最小二乘法，求解如下矩阵，解出来各个交点(u,v)所对应的像素坐标系下的3d坐标(xc,yc,zc)。
[0071][0072][0073]
[18]步骤s6，遍历当前组ii类图像的所有交点，循环统计出3d激光点在当前组ii类图像所对应的标定平面下的，相机坐标系下的3d空间坐标(xc,yc,zc)。
[0074]
[19]遍历各组ii类图像的所有交点，循环统计出3d激光点在各个ii类图像所对应的各个标定平面下的，相机坐标系下的3d空间坐标(xc,yc,zc)。
[0075]
对本次遍历的结果进行3d平面拟合，即可得到该激光平面相对于相机的空间坐标，即所要求的结果。
[0076]
[20]步骤s7，根据步骤[19]得到的激光平面方程，与相机坐标的x轴坐标的交点，就是虚拟的激光发射点，激光入射角度可以根据需要进行定义。
[0077]
[21]循环执行步骤[16]到步骤[20]，遍历各个激光线，即可分别得到每条激光线相对于相机位置的激光平面方程。