一种激光雷达标定方法及相关设备与流程

1.本说明书涉及激光雷达领域，更具体地说，本发明涉及一种激光雷达标定方法及相关设备。


背景技术：

2.自动驾驶是汽车产业与人工智能、物联网、高性能计算等新一代信息技术深度融合的产物，是当前全球汽车与交通出行领域智能化和网联化发展的主要方向。环境感知是自动驾驶的关键技术主之一。目前应用于自动驾驶汽车环境感知的车载传感器主要有摄像头、激光雷达、毫米波雷达、超声波雷达、惯性导航等。激光雷达是一种激光测距系统，具有：分辨率高、精度高，抗有源干扰能力强，获取的信息量丰富的特点。随着商用gps和imu的发展，通过激光雷达从移动平台上获得高精度的数据已成为现实并被广泛应用。
3.各车载传感器都以一定的角度和位置安装在车辆上，在多个激光雷达传感器存在的情况下，为了实现激光点云数据的有效融合，需要将各激光雷达之间的相对外参进行标定，将各激光雷达转换到统一的时空坐标系下，以便于后续的数据处理分析。当前，多激光雷达标定技术得到了一定的发展，但也仍然存在一些不足：一些标定方法需要配置复杂的标定环境，并借助定制靶标或全站仪等辅助设备，理论上精度较高，但同时受到操作人员技术熟练程度的影响较大。且此方法耗时耗力，且自动化程度低，难以满足自动驾驶的实际需求。另一些方法中，直接利用点云配准获取坐标转换参数并通过迭代收敛得到最优解的方法，容易陷入局部收敛而导致标定结果不佳。


技术实现要素：

4.在发明内容部分中引入了一系列简化形式的概念，这将在具体实施方式部分中进一步详细说明。本发明的发明内容部分并不意味着要试图限定出所要求保护的技术方案的关键特征和必要技术特征，更不意味着试图确定所要求保护的技术方案的保护范围。
5.第一方面，本发明提出一种激光雷达标定方法，上述方法包括：
6.获取每个激光雷达在目标区域的多组点云数据，其中，上述激光雷达的数量至少为两个，每组上述点云数据包括上述目标区域内的同一段平整路面和至少一个相同的标志性物体；
7.对每组上述点云数据进行预处理以获取待处理点云数据，其中，上述预处理包括阈值筛选和下采样操作；
8.利用地面点云数据进行水平校准以获取第一旋转矩阵和第一平移向量，并利用获取的上述第一旋转矩阵和上述第一平移向量对上述待处理点云数据进行转换，其中，上述待处理点云数据包括地面点云数据和非地面点云数据；
9.根据水平校准后的非地面点云数据基于icp算法获取第二旋转矩阵和第二平移向量；
10.基于根据上述第一旋转矩阵、第二旋转矩阵、第一平移向量和第二平移向量确定
每个激光雷达的目标旋转矩阵和目标平移向量以标定上述激光雷达。
11.可选的，上述对每组上述点云数据进行预处理以获取待处理点云数据，包括：
12.基于预设阈值对每组上述点云数据进行数据裁剪以获取包括上述平整路面和上述标志性物体在内的裁剪点云数据；
13.对上述裁剪点云数据进行下采样以获取上述待处理点云数据。
14.可选的，上述方法还包括：
15.利用上述待处理点云数据基于ransac平面检测方法提取上述地面点云数据和上述非地面点云数据；
16.拟合上述地面点云数据获取拟合平面法向量。
17.可选的，上述利用地面点云数据进行水平校准以获取第一旋转矩阵和第一平移向量，并利用获取的上述第一旋转矩阵和上述第一平移向量对上述待处理点云数据进行转换，包括：
18.根据上述地面点云数据、上述拟合平面法向量获取第一旋转矩阵；
19.求取每组上述点云数据在垂直地面方向的平均值作为第一平移向量；
20.基于上述第一旋转矩阵、上述第一平移向量对上述点云数据进行水平校准。
21.可选的，上述根据水平校准后的非地面点云数据基于icp算法获取第二旋转矩阵和第二平移向量，包括：
22.以每个激光雷达对上述标志性物体的大致估算位置关系作为旋转平移参数初始值；
23.利用上述旋转平移参数初始值、上述非地面点云数据通过icp算法进行迭代计算，将满足预设精度要求的旋转矩阵和平移向量作为上述第二旋转矩阵和上述第二平移向量。
24.可选的，上述利用上述旋转平移参数初始值、上述非地面点云数据通过icp算法进行迭代计算，将满足预设精度要求的旋转矩阵和平移向量作为上述第二旋转矩阵和上述第二平移向量，包括：
25.以一组上述非地面点云数据作为目标点云，其余上述非地面点云数据作为待配准点云，并确定目标对应点；
26.基于上述目标对应点采用最小二乘法迭代求取旋转矩阵和平移向量；
27.在点云均方误差小于预设阈值的情况下，将当前求取的旋转矩阵和平移矩阵确定为上述第二旋转矩阵和上述第二平移向量。
28.可选的，上述基于根据上述第一旋转矩阵、第二旋转矩阵、第一平移向量和第二平移向量确定每个激光雷达的目标旋转矩阵和目标平移向量以标定上述激光雷达，包括：
29.将上述第一旋转矩阵右乘上述第二旋转矩阵的结果作为上述目标旋转矩阵；
30.将上述第一平移向量右乘所述第二旋转矩阵的结果和上述第二平移向量之和确定为上述目标平移向量；
31.基于上述目标旋转矩阵和上述目标平移向量标定上述激光雷达。
32.第二方面，本发明还提出一种激光雷达标定装置，包括：
33.第一获取单元，用于获取每个激光雷达在目标区域的多组点云数据，其中，上述激光雷达的数量至少为两个，每组上述点云数据包括上述目标区域内的同一段平整路面和至少一个相同的标志性物体；
34.第二获取单元，用于对每组上述点云数据进行预处理以获取待处理点云数据，其中，上述预处理包括阈值筛选和下采样操作；
35.第三获取单元，利用地面点云数据进行水平校准以获取第一旋转矩阵和第一平移向量，并利用获取的上述第一旋转矩阵和上述第一平移向量对上述待处理点云数据进行转换，其中，上述待处理点云数据包括地面点云数据和非地面点云数据；
36.第四获取单元，用于根据水平校准后的非地面点云数据基于icp算法获取第二旋转矩阵和第二平移向量；
37.标定单元，用于基于根据上述第一旋转矩阵、第二旋转矩阵、第一平移向量和第二平移向量确定每个激光雷达的目标旋转矩阵和目标平移向量以标定上述激光雷达。
38.第三方面，一种电子设备，包括：存储器、处理器以及存储在上述存储器中并可在上述处理器上运行的计算机程序，上述处理器用于执行存储器中存储的计算机程序时实现如上述的第一方面任一项的激光雷达标定方法的步骤。
39.第四方面，本发明还提出一种计算机可读存储介质，其上存储有计算机程序，上述计算机程序被处理器执行时实现第一方面上述任一项的激光雷达标定方法。
40.综上，本技术实施例的激光雷达标定方法包括：获取每个激光雷达在目标区域的多组点云数据，其中，上述激光雷达的数量至少为两个，每组上述点云数据包括上述目标区域内的同一段平整路面和至少一个相同的标志性物体；对每组上述点云数据进行预处理以获取待处理点云数据，其中，上述预处理包括阈值筛选和下采样操作；利用地面点云数据进行水平校准以获取第一旋转矩阵和第一平移向量，并利用获取的第一旋转矩阵和第一平移向量对待处理点云数据进行转换，其中，上述待处理点云数据包括地面点云数据和非地面点云数据；根据水平校准后的非地面点云数据基于icp算法获取第二旋转矩阵和第二平移向量；基于根据上述第一旋转矩阵、第二旋转矩阵、第一平移向量和第二平移向量确定每个激光雷达的目标旋转矩阵和目标平移向量以标定上述激光雷达。本技术实施例提供的激光雷达标定方法，通过分步标定的形式，通过在平整路面采集到的点云数据，提取出地面点云并作水平校准得到相应的第一旋转矩阵和第一平移向量，再利第一旋转矩阵和第一平移向量对非地面点云进行坐标转换后，基于icp点云配准算法进行迭代求解，得到第二旋转矩阵和第二平移向量，基于两次所得坐标转换参数可最终获取所需目标旋转矩阵和目标平移向量。本技术提出的激光雷达标定方法操作简单，不需要设置复杂的标定环境，不需要其他额外的辅助设备，只需选择一个合适的场景采集数据后，利用本技术提供的快速标定方法进行操作和计算即可得到有效的标定参数。
41.本发明的激光雷达标定方法，本发明的其它优点、目标和特征将部分通过下面的说明体现，部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。
附图说明
42.通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本说明书的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
43.图1为本技术实施例提供的一种激光雷达标定方法流程示意图；
44.图2为本技术实施例提供的一种激光雷达标定的坐标系示意图；
45.图3为本技术实施例提供的一种激光雷达标定模型结构示意图；
46.图4为本技术实施例提供的一种激光雷达标定装置结构示意图；
47.图5为本技术实施例提供的一种激光雷达标定电子设备结构示意图。
具体实施方式
48.本技术实施例提供的激光雷达标定方法，通过分步标定的形式，通过在平整路面采集到的点云数据，提取出地面点云并作水平校准得到相应的第一旋转矩阵和第一平移向量，再利第一旋转矩阵和第一平移向量对非地面点云进行坐标转换后，基于icp点云配准算法进行迭代求解，得到第二旋转矩阵和第二平移向量，基于两次所得坐标转换参数可最终获取所需目标旋转矩阵和目标平移向量。本技术提出的激光雷达标定方法操作简单，不需要设置复杂的标定环境，不需要其他额外的辅助设备，只需选择一个合适的场景采集数据后，利用本技术提供的快速标定方法进行操作和计算即可得到有效的标定参数。
49.本技术的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”、“第四”等(如果存在)是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的实施例能够以除了在这里图示或描述的内容以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。
50.请参阅图1，为本技术实施例提供的一种激光雷达标定方法流程示意图，具体可以包括：
51.s110、获取每个激光雷达在目标区域的多组点云数据，其中，上述激光雷达的数量至少为两个，每组上述点云数据包括上述目标区域内的同一段平整路面和至少一个相同的标志性物体；
52.示例性的，随着车辆智能化程度提升，车辆会安装多个激光雷达传感器，各激光雷达传感器在安装过程中，由于受到实际操作环境的限制、以及操作人员的技术熟练度等的影响，其安装位置和角度与理想状态都会有一定程度的偏离。因此需要在出厂时，或者长时间使用后对传感器进行校准。以图2所示的车辆为例，三个激光雷达分别安装于车身顶部前方中间位置以及左右后视镜附近位置，分别为中部主激光雷达、左翼激光雷达和右翼激光雷达。对于常规的小汽车，由于车身高度有限，三个激光雷达基本都是尽量保持与车身平行；但对于自动驾驶小巴或中巴，由于车辆本身有一定高度，为了能够更有效地进行数据采集，需要将中部激光雷达向地面调整合适的角度(即pitch横滚角)，此时pitch角不再是接近于0的数，具体情况视激光雷达和车辆情况而定。
53.以图2中的车辆为例，t
xoy
定义坐标原点为中间主激光雷达中心点正下方在地面的投影点，x轴为指向中间主激光雷达前方且平行于地面，y轴为指向中间主激光雷达左侧(垂直于x轴)，且平行于地面，z轴为垂直于x轴y轴且符合右手法则。t
xoy
为标定的目标坐标系，m
xoy
为中部主激光雷达原始坐标系，l
xoy
和r
xoy
分别表示左翼和右翼激光雷达原始坐标系。需
要说明的是，本技术实施例只是为了举例说明，具体的激光雷达不限于三个，其安装位置也不做限制。
54.在进行激光雷达标定时，需要在在一个路面平整、且有标志性地物或目标的场景内进行数据采集；将采集车开至目标地点，静止状态下开始采集数据，多个激光雷达获取的数据中同时包含同一块平整地面和一定数量的相同的标志性物体，标志性物体可以为路牌，路旁的建筑物等，中部、左部和右部的激光雷达采集的三组点云数据为可以命名为cloud_middle、cloud_left、cloud_right。
55.s120、对每组上述点云数据进行预处理以获取待处理点云数据，其中，上述预处理包括阈值筛选和下采样操作；
56.示例性的，多激光雷达系统中，可能出现不同类型和规格的激光雷达，其获取的数据在数据形式及点云密度上均可能存在差异，为了提高标定精度并加快处理速度，需首先进行点云子集提取和降采样处理，从而获取待处理点云数据。
57.s130、利用地面点云数据进行水平校准以获取第一旋转矩阵和第一平移向量，并利用获取的第一旋转矩阵和第一平移向量对待处理点云数据进行转换，其中，上述待处理点云数据包括地面点云数据和非地面点云数据；
58.示例性的，对待处理点云数据采用相应的算法分割可以获取地面点云数据和非地面点云数据，对地面点云数据进行水平校准可以获取第一旋转矩阵和第一平移向量，并利用获取的第一旋转矩阵和第一平移向量对待处理点云数据进行转换。在标定的过程中，需要对x、y、z三个方向上的位置和绕三个轴的旋转做出标定，在进行水平校准后，z方向上的平移，绕x轴的旋转roll，绕y轴的旋转pitch的校准已经完成。
59.s140、根据水平校准后的非地面点云数据基于icp算法获取第二旋转矩阵和第二平移向量；
60.示例性的，对于水平校准后的非地面点云数据采用icp算法对x、y方向上的平移和绕z的旋转进行标定。需要说明的是，icp(iterative closest point，迭代最近点)算法是基于数据配准法，利用最近点搜索法通过寻找对应点进行旋转和平移优化。
61.s150、基于根据上述第一旋转矩阵、第二旋转矩阵、第一平移向量和第二平移向量确定每个激光雷达的目标旋转矩阵和目标平移向量以标定上述激光雷达。
62.示例性的，根据水平校准得到的第一旋转矩阵和第一平移向量，以及采用icp算法获取的第二平移向量和第二旋转矩阵可以确定每个激光雷达对应的目标旋转矩阵和目标平移向量，从而完成各个激光雷达的标定。
63.综上，本技术实施例提供的激光雷达标定方法，通过分步标定的形式，通过在平整路面采集到的点云数据，提取出地面点云并作水平校准得到相应的第一旋转矩阵和第一平移向量，再利第一旋转矩阵和第一平移向量对非地面点云进行坐标转换后，基于icp点云配准算法进行迭代求解，得到第二旋转矩阵和第二平移向量，基于两次所得坐标转换参数可最终获取所需目标旋转矩阵和目标平移向量。本技术提出的激光雷达标定方法操作简单，不需要设置复杂的标定环境，不需要其他额外的辅助设备，只需选择一个合适的场景采集数据后，利用本技术提供的快速标定方法进行操作和计算即可得到有效的标定参数。
64.在一些示例中，上述对每组上述点云数据进行预处理以获取待处理点云数据，包括：
65.基于预设阈值对每组上述点云数据进行数据裁剪以获取包括上述平整路面和上述标志性物体在内的裁剪点云数据；
66.对上述裁剪点云数据进行下采样以获取上述待处理点云数据。
67.示例性的，对于激光雷达获取的各组点云数据，选择合适的预设阈值进行数据剪裁，剔除多余的环境信息，提取包含地面点云及地物目标点云在内的点云子集(subcloud_m、subcloud_l、subcloud_r)，同时基于voxelgrid方法对高点云密度点云数据进行下采样，使各组点云数据量在同一数量级，便于后续处理。
68.综上，本技术实施例提供的激光雷达标定方法，通过对不同激光雷达获取的点云数据进行预处理，使其拥有相同的等级的数据量，可以使得不同类型的激光雷达数据有效融合，互相标定。
69.在一些示例中，上述方法还包括：
70.利用上述待处理点云数据基于ransac平面检测方法提取上述地面点云数据和上述非地面点云数据；
71.拟合上述地面点云数据获取拟合平面法向量。
72.示例性的，地面点云分割及拟合平面的法向量计算：对于经过阈值处理和下采样处理后的待处理点云数据，分别基于ransac的平面检测方法提取地面点云(pc_ground_m、pc_ground_l、pc_ground_r)，拟合点云平面并获取法向量(normv_m、normv_l、normv_r)，点云子集的其他点云则为非地面点云(pc_nonground_m、pc_nonground_l、pc_nonground_r)。需要说明的是，ransac算法(随机抽样一致算法)假设数据中包含正确数据和异常数据。同时ransac也假设，给定一组正确的数据，存在可以计算出符合这些数据的模型参数的方法。该算法核心思想就是随机性和假设性，随机性是根据正确数据出现概率去随机选取抽样数据，根据大数定律，随机性模拟可以近似得到正确结果。假设性是假设选取出的抽样数据都是正确数据，然后用这些正确数据通过问题满足的模型，去计算其他点，然后对这次结果进行一个评分，从而提取地面点云、拟合点云平面并获取法向量并提出非地面点云。
73.综上，本技术实施例提供的激光雷达标定方法，通过ransac将地面与非地面点云数据有效提取，利用地面数据进行水平标定，从而获取水平标定后的第一旋转矩阵和第一平移向量。
74.在一些示例中，上述利用地面点云数据进行水平校准以获取第一旋转矩阵和第一平移向量，并利用获取的上述第一旋转矩阵和上述第一平移向量对上述待处理点云数据进行转换，包括：
75.根据上述地面点云数据、上述拟合平面法向量获取第一旋转矩阵；
76.求取每组上述点云数据在垂直地面方向的平均值作为第一平移向量；
77.基于上述第一旋转矩阵、上述第一平移向量对上述点云数据进行水平校准。
78.示例性的，地面点云水平校准的旋转矩阵计算：对于获取的三组地面点云，根据三个拟合平面的法向量与水平面的法向量(0,0,1)计算旋转矩阵(rmat1_m、rmat1_l、rmat1_r)，基于此可实现各组地面点云和非地面点云的水平校准。令上述拟合平面的法向量normvm,normv
l
,normvr分别为：
水平面的法向量为：水平面的法向量为：则依据已知法向量，可分别求出旋转至的旋转矩阵，以为例：
79.旋转到的旋转角为θ，有：
[0080][0081]
旋转到的旋转轴为有：
[0082][0083]
旋转轴方向的单位向量可表示为：令其中其中
[0084]
根据罗德里格旋转公式，旋转到的旋转矩阵为：
[0085][0086]
同理可得左翼右翼激光雷达旋转矩阵rmat1
l
、rmat1r。
[0087]
对地面点云平面获取统一的平移参数：经过水平校准的各组地面点云拟合于不同的水平面，并对各组地面点云z坐标值分别计算z值平均值(arvzm,arvz
l
,arvzr)，依据各平均值获取z方向平移参数，基于此实现各激光雷达地面点云z坐标值统一。此时各平移向量v1_transm、v1_trans
l
、v1_transr为：
[0088]
v1_transm＝(0,0,―arvzm)
[0089]
v1_trans
l
＝(0,0,―arvz
l
)
[0090]
v1_transr＝(0,0,―arvzr)
[0091]
多组激光点云的水平校准：基于前面两步的旋转矩阵和平移参数，分别对地面点云和非地面点云进行坐标转换。(完成水平校准和z坐标值统一的地面点云为ground_trans1_m、ground_trans1_l、ground_trans1_r；完成两步操作的非地面点云为pc_nonground_trans1_m、pc_nonground_trans1_l、pc_nonground_trans1_r。其中：
[0092]
pc_nonground_trans1_m＝rmat1m
·
pc_nonground_m+v1_transm[0093]
pc_nonground_trans1_l＝rmat1l
·
pc_nonground_l+v1_trans
l
[0094]
pc_nonground_trans1_r＝rmat1r
·
pc_nonground_r+v1_transr[0095]
综上，本技术提出的激光雷达标定方法，通过激光点云水平校准可有效保证水平方向的准确性和精确度，在激光点云配准过程中由于多组点云数据不完全相同容易陷入局部收敛而导致结果错误，本技术提出的基于地面点云水平校准的分步标定的方法可有效的解决此问题。
[0096]
在一些示例中，上述根据水平校准后的非地面点云数据基于icp算法获取第二旋转矩阵和第二平移向量，包括：
[0097]
以每个激光雷达对上述标志性物体的大致估算位置关系作为旋转平移参数初始值；
[0098]
利用上述旋转平移参数初始值、上述非地面点云数据通过icp算法进行迭代计算，将满足预设精度要求的旋转矩阵和平移向量作为上述第二旋转矩阵和上述第二平移向量。
[0099]
示例性的，对于经过水平校准后的非地面点云，只考虑点云xoy平面内的转换情况，可通过对点云的可视化，直接目视分析标志性物体得到其水平方向旋转角度及平移参数，作为初值。以中部主激光雷达的点云为目标点云，左翼和右翼激光点云为待配准点云即源点云，利用上一步骤中所得的初值，基于icp算法进行点云配准，并进行迭代计算获得转换参数最优解，包括旋转矩阵rmat2m、rmat2
l
、rmat2r和平移参数v2_transm、v2_trans
l
、v2_transr，其中rmat2m为单位阵，v2_transm＝(0,0,0)
t
。需要说明的是，可以以任意一个激光雷达的点云为目标点云，其余的点云数据为待配准点云。
[0100]
综上，本技术提出的激光雷达标定方法，在进行水平校准后通过icp算法可以将剩余的参数进行标定，从而通过分步的方式完成激光雷达6参数的标定工作。
[0101]
在一些示例中，上述利用上述旋转平移参数初始值、上述非地面点云数据通过icp算法进行迭代计算，将满足预设精度要求的旋转矩阵和平移向量作为上述第二旋转矩阵和上述第二平移向量，包括：
[0102]
以一组上述非地面点云数据作为目标点云，其余上述非地面点云数据作为待配准点云，并确定目标对应点；
[0103]
基于上述目标对应点采用最小二乘法迭代求取旋转矩阵和平移向量；
[0104]
在点云均方误差小于预设阈值的情况下，将当前求取的旋转矩阵和平移矩阵确定为上述第二旋转矩阵和上述第二平移向量。
[0105]
示例性的，通过icp算法进行迭代计算的具体过程可以包括步骤s210至步骤s230：
[0106]
s210、寻找源点云与目标点云对应点。以pc_nonground_trans1_m(用pm表示)为目标点云，将源点云pc_nonground_trans1_l(用p
l
表示)，基于标志性物体的估算位置关系作为旋转平移参数初始值进行坐标转换，将变换后的点云p
l’与目标点云pm进行对比，当点云之间的距离值小于设定阈值即为对应点；
[0107]
s220、坐标转换参数优化。待求参数的自由度为6，结合步骤s210中确定的对应点，采用最小二乘法进行解算，优化坐标转换参数；
[0108]
s230、迭代求解得到坐标转换参数最优解。利用s220中得到的优化参数重新寻找对应点并进行最小二乘求解坐标转换参数，循环迭代步骤s210和s220直至满足点云均方误差小于阈值。目标函数为：
[0109][0110]
其中n
p
为对应点总数，在f(r2,v2)小于阈值时对应的r2
l
和v2
l
，即为求得最优解rmat2
l
、v2_trans
l
。右翼激光雷达坐标转换参数rmat2r、v2_transr同理可得。
[0111]
在一些示例中，上述基于根据上述第一旋转矩阵、第二旋转矩阵、第一平移向量和第二平移向量确定每个激光雷达的目标旋转矩阵和目标平移向量以标定上述激光雷达，包括：
[0112]
将上述第一旋转矩阵右乘上述第二旋转矩阵的结果作为上述目标旋转矩阵；
[0113]
将上述第一平移向量右乘所述第二旋转矩阵的结果和上述第二平移向量之和确定为上述目标平移向量；
[0114]
基于上述目标旋转矩阵和上述目标平移向量标定上述激光雷达。
[0115]
示例性的，将上述第一旋转矩阵rmat1右乘上述第二旋转矩阵rmat2的结果作为上述目标旋转矩阵rmat；将上述第一平移向量v1_trans右乘上述第二旋转矩阵rmat2的结果和上述第二平移向量v2_trans之和确定为上述目标平移向量v_trans，具体如下式所示：
[0116]
rmatm＝rmat2m·
rmat1m[0117]
rmat
l
＝rmat2
l
·
rmat1
l
[0118]
rmatr＝rmat2r·
rmat1r[0119]
v_transm＝rmat2m·
v1_transm+v2_transm[0120]
v_trans
l
＝rmat2
l
·
v1_trans
l
+v2_trans
l
[0121]
v_transr＝rmat2r·
v1_transr+v2_transr[0122]
其中，上式中的下角标m代表中部主激光雷达、l代表左翼激光雷达、r代表右翼激光雷达，同时，可依据所需数据格式要求，通过数据转换输出旋转矩阵、四元数、欧拉角rpy(roll、pitch、yaw)等。
[0123]
在一些示例中，请参阅图3，为本技术实施例提供的一种激光雷达标定模型结构示意图，通过数据采集模块获取多个传感器的多组激光点云数据，将激光点云数据利用数据预处理模块进行阈值筛选和下采样获得待处理点云数据，将待处理点云数据通过基于平面检测的激光点云水平校准模块，基于ransac法将点云分割，获得地面点云和非地面点云，对于地面点云采用基于法向量的地面点云水平校准，并将多组地面点云平面统一化处理，从而获取旋转参数1(即第一旋转矩阵)和平移参数1(即第一平移向量)，并通过基于icp算法的点云配准模块，针对非地面点云利用旋转参数1和平移参数1获取坐标转换后的非地面点云，确定配准的坐标转换初值，对进行坐标转换后的非地面点云基于icp算法进行点云配准从而获取旋转参数2(即第二旋转矩阵)和平移参数2(即第二平移向量)，通过标定参数确定模块，根据第一旋转矩阵和第二旋转矩阵确定目标旋转矩阵，再根据第一平移向量和第二平移向量并结合第二旋转矩阵确定目标平移向量，基于激光雷达标定外参的数据格式进行数据转换，从而完成激光雷达的标定工作。
[0124]
请参阅图4，本技术实施例中激光雷达标定装置的一个实施例，可以包括：
[0125]
第一获取单元21，用于获取每个激光雷达在目标区域的多组点云数据，其中，上述激光雷达的数量至少为两个，每组上述点云数据包括上述目标区域内的同一段平整路面和
至少一个相同的标志性物体；
[0126]
第二获取单元22，用于对每组上述点云数据进行预处理以获取待处理点云数据，其中，上述预处理包括阈值筛选和下采样操作；
[0127]
第三获取单元23，利用地面点云数据进行水平校准以获取第一旋转矩阵和第一平移向量，并利用获取的第一旋转矩阵和第一平移向量对待处理点云数据进行转换，其中，上述待处理点云数据包括地面点云数据和非地面点云数据；
[0128]
第四获取单元24，用于根据水平校准后的非地面点云数据基于icp算法获取第二旋转矩阵和第二平移向量；
[0129]
标定单元25，用于基于根据上述第一旋转矩阵、第二旋转矩阵、第一平移向量和第二平移向量确定每个激光雷达的目标旋转矩阵和目标平移向量以标定上述激光雷达。
[0130]
如图5所示，本技术实施例还提供一种电子设备300，包括存储器310、处理器320及存储在存储器320上并可在处理器上运行的计算机程序311，处理器320执行计算机程序311时实现上述激光雷达标定的任一方法的步骤。
[0131]
由于本实施例所介绍的电子设备为实施本技术实施例中一种激光雷达标定装置所采用的设备，故而基于本技术实施例中所介绍的方法，本领域所属技术人员能够了解本实施例的电子设备的具体实施方式以及其各种变化形式，所以在此对于该电子设备如何实现本技术实施例中的方法不再详细介绍，只要本领域所属技术人员实施本技术实施例中的方法所采用的设备，都属于本技术所欲保护的范围。
[0132]
在具体实施过程中，该计算机程序311被处理器执行时可以实现图1对应的实施例中任一实施方式。
[0133]
需要说明的是，在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详细描述的部分，可以参见其它实施例的相关描述。
[0134]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0135]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式计算机或者其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0136]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0137]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或
其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0138]
本技术实施例还提供了一种计算机程序产品，该计算机程序产品包括计算机软件指令，当计算机软件指令在处理设备上运行时，使得处理设备执行如图1对应实施例中的激光雷达标定的流程。
[0139]
计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行计算机程序指令时，全部或部分地产生按照本技术实施例的流程或功能。计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一计算机可读存储介质传输，例如，计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(digital subscriber line，dsl))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。计算机可读存储介质可以是计算机能够存储的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，dvd)、或者半导体介质(例如固态硬盘(solid state disk，ssd))等。
[0140]
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0141]
在本技术所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
[0142]
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0143]
另外，在本技术各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0144]
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本技术各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、磁碟或者光盘等各种可以存储程序代码的介质。
[0145]
以上，以上实施例仅用以说明本技术的技术方案，而非对其限制；尽管参照前述实施例对本技术进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各
实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本技术各实施例技术方案的精神和范围。