-
车载移动测量系统是一种多传感器集成系统,各传感器均有独立坐标系,系统进行测量工作时,系统需根据数学模型将激光扫描仪获取的角度和距离信息转换成三维空间坐标,再经过多次坐标转换,获得具有地理参考信息的激光点云数据。转换过程中涉及到的坐标系有:激光扫描仪坐标系(Scanner’s Own Coordinate System,SOCS)、载体坐标系(Body Frame System,BFS)、当地水平坐标系(Local-Level Frame System,LLS)、地心地固坐标系(Earth Centered Earth Fixed Frame System,ECEF),各坐标系间转换关系如图1所示。
车载移动测量系统通过激光扫描仪获取原始扫描数据后,经过图1所示的一系列坐标转换获得具有地理参考信息的点云数据,定位模型如公式(1)所示:
$$ {{{P}}_{ECEF}} = {{{T}}^{ECEF}} + \left[ {\left[ {{{{P}}_{SOCS}} {{R}}_{SOCS}^{BFS} + {{T}}_{SOCS}^{BFS}} \right] {{R}}_{BFS}^{LLS}} \right]{{R}}_{LLS}^{ECEF} $$ (1) 式中:
${{{P}}_{ECEF}} = {[ \!\!{\begin{array}{*{20}{c}} X&Y&Z \end{array}} \!\!]_{ECEF}}$ 为激光脚点在地心地固坐标系下的坐标;${{{T}}^{ECEF}} = [\!\! {\begin{array}{*{20}{c}} {{X_0}}&{{Y_0}}&{{Z_0}} \end{array}}\!\!]$ 为当地水平坐标系原点在地心地固坐标系下平移向量;${{{P}}_{SOCS}}\left( {\rho ,\theta ,\phi } \right)$ 为激光脚点在激光扫描坐标系下坐标,$\rho $ 为目标点到扫描中心斜距,$\theta $ 为激光扫描仪测量目标点时的竖直角,$\phi $ 为激光扫描仪测量目标点时的水平角;${{R}}_{SOCS}^{BFS}\left( {\alpha ,\beta ,\gamma } \right)$ 为扫描仪与载体坐标系间的旋转矩阵,$\alpha $ ,$\beta $ ,$\gamma $ 是安置参数角元素;${{T}}_{SOCS}^{BFS} = [\!\!\! {\begin{array}{*{20}{c}} {\Delta X}\!\!&\!\!{\Delta Y}\!\!&\!\!{\Delta Z} \end{array}} \!\!\!]$ 为扫描仪与载体坐标系间的平移向量;${{R}}_{BFS}^{LLS}\left( {r,p,h} \right)$ 为载体坐标系与当地水平坐标系间的旋转矩阵,$r$ ,$p$ ,$h$ 分别为组合导航系统所测的翻滚角、俯仰角、航向角;${{R}}_{LLS}^{ECEF}\left( {B,L} \right)$ 为当地水平坐标系与地心地固坐标系间的旋转矩阵,$B$ ,$L$ 分别为当地水平坐标系原点${O_L}$ 在地心地固坐标系下的大地纬度和大地经度。 -
车载移动测量系统集成后,受加工和安装精度的限制,安置误差难以避免且无法通过直接观测手段获得误差大小。根据扫描点定位模型可以看出,系统最终精度受多种因素影响,此外还受GNSS天线相位中心与IMU中心杆臂值量取误差影响、控制点坐标误差影响,基站仪器高量取误差影响,以上误差在自检校方法中无法准确体现,且容易引入人为误差。为更好地实现车载移动测量安置误差流程化检校以及对系统精度进行完善的评估,文中依托高精度三维检校场利用不同朝向参考面特征实现车载移动测量系统安置误差检校。
-
高精度三维检校场建立在山东科技大学校园内,包含教学楼、实验楼、标志性建筑物、学生公寓和部分地面等。检校场附近为通畅、宽阔的直行道路且两侧树木矮小,采集区域GNSS信号良好、可完整扫描参考面。参考面模型平整且边缘易识别,包含不同方向、不同远近的平面,分布如图2所示。
检校场测量遵循“先控制后碎部”的原则,在高精度GNSS控制网基础上,使用Leica TM30高精度全站仪采用免棱镜模式进行参考面测量,每个参考面上均匀分布15~30个测量点,车载移动测量系统安置误差精度及可靠性很大程度上取决于检校场模型精度,文中对所有参考平面进行平面拟合并分析误差,所得拟合误差如图3所示。校园内93个朝向不同的拟合参考面误差主要集中于1~8 mm,中误差为3.7 mm,所测参考面能够满足绝大多数车载移动测量系统检校需求。
-
由系统定位公式可知,最终的点云数据涉及多个传感器提供的数据,但各传感器均存在误差。为提高检校模型稳定性,获得最佳参数,文中选用附有参数的条件平差模型,将每个传感器的观测值误差一并考虑到解算模型中。
观测值:
$$ {{L}} = \left[ {\begin{array}{*{20}{c}} \rho &\theta &\phi &r&p&h&{{X_0}}&{{Y_0}}&{{Z_0}} \end{array}} \right] $$ (2) 安置参数:
$$ {{X}} = \left[ {\begin{array}{*{20}{c}} \alpha &\beta &\gamma &{\Delta X}&{\Delta Y}&{\Delta Z} \end{array}} \right] $$ (3) 参考平面方程如公式(4)所示:
$$ Ax + By + Cz - D = 0 $$ (4) 式中:A、B、C、D为参考面方程参数。
通过车载移动测量系统对参考平面进行数据采集,每个参考面上的扫描点云均应满足对应参考面方程,这里以点到平面距离平方和最小作为约束条件,基于平面的数学约束模型如公式(5)所示:
$$ \begin{split} {f_p} = & {f_p}\left( {{\hat{L}},{\hat{X}}} \right) = {f_p}\left( {{{L}} + {{V}},{{{X}}^{\rm{0}}} + {\hat{x}}} \right) = \\ & \frac{{{{\left( {{A_p}{x_i} + {B_p}{y_i} + {C_p}{z_i} - {D_p}} \right)}^2}}}{{A_p^2 + B_p^2 + C_p^2}} \end{split} $$ (5) 式中:
${{{X}}^{\rm{0}}}$ 表示检校参数初始值;${{V}}$ 为观测值改正数;${\hat{X}}$ 为检校参数初始值的改正数;${x_i}$ 、${y_i}$ 、${z_i}$ 为激光脚点i在地心地固坐标系下坐标; Ap、Bp、Cp、Cp为i扫描点所在平面对应的约束参考面P的平面方程参数。 -
安置参数中包含角元素和线元素,在车载移动测量系统中,安置角误差为主要影响因素。初始参数获取过程中,线元素精度要高于角元素,因传感器精度不同,线元素误差基本与设计值一致,通常在毫米级。另外,考虑到角元素与线元素对点云精度影响具有相关性,所以采用分步平差方法,首先解算影响较大的安置角误差,在解算结果的基础上再以平移量作为参数,将剩余误差再分配进行平移参数求解,最终实现安置误差高精度检校。
车载移动测量系统作为一种多传感器集成设备,各传感器本身也存在一定误差。为简化模型,文中认为不同传感器测得观测值之间不具有相关性[25],权阵可表示为对角阵,如公式(6)所示:
$$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\; {{P}} = {{diag}}\left( {\begin{array}{*{20}{c}} {{{P}}_p^1}&{{{P}}_p^2}& \cdots &{{{P}}_p^{i - 1}}&{{{P}}_p^i} \end{array}} \right) \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\; {{P}}_p^i = {{diag}}\left( {\begin{array}{*{20}{c}} {{{P}}_i^1}&{{{P}}_i^2}& \cdots &{{{P}}_i^{j - 1}}&{{{P}}_i^j} \end{array}} \right) \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\; {{P}}_i^j = {{diag}}\left( {\begin{array}{*{20}{c}} {{{P}}_j^1}&{{{P}}_j^2}& \cdots &{{{P}}_j^{k - 1}}&{{{P}}_j^k} \end{array}} \right) \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\; {\begin{aligned} {{P}}_j^k = & {{diag}} \left( {\frac{1}{{\sigma _{{\rho _k}}^2}}} \quad {\frac{1}{{\sigma _{{\theta _k}}^2}}} \quad {\frac{1}{{\sigma _{{\phi _k}}^2}}} \quad {\frac{1}{{\sigma _{{r _k}}^2}}} \quad {\frac{1}{{\sigma _{{p _k}}^2}}} \right. \\ & \left. {\frac{1}{{\sigma _{{h _k}}^2}}} \quad {\frac{1}{{\sigma _{{X_0}_k}^2}}} \quad {\frac{1}{{\sigma _{{Y_0}_k}^2}}} \quad {\frac{1}{{\sigma _{{Z_0}_k}^2}}} \right) \end{aligned} } \end{array} $$ (6) 式中:i表示检校实验中第i个平面集;j表示平面集i中第j个平面约束模型;k表示第平面约束模型j上第k个激光脚点。参数解算过程涉及的
${{A}}$ 、${{B}}$ 、${{W}}$ 矩阵的形式与${{P}}$ 类似,${{A}}$ 为对角矩阵,${{B}}$ 、${{W}}$ 为列矩阵,为简化形式,下文只写出第k个激光脚点的方程。首先将安置角作为参数,将函数模型线性化,得安置角改正数条件方程如公式(7)所示:
$$ \mathop {{A}}\limits_{n,m} \mathop {{V}}\limits_{m} + \mathop {{{{B}}}}\limits_{n,{r}} \mathop {{{{\hat{X}}}}}\limits_{{r},} + \mathop {{W}}\limits_{n} =\mathop {{0}}\limits_{n} $$ (7) 系数矩阵
${{A}}$ :$$ \begin{aligned} {{A}}_j^k = & \left[ {\frac{{\partial f_p^k}}{{\partial \rho }}} \quad {\frac{{\partial f_p^k}}{{\partial \theta }}} \quad {\frac{{\partial f_p^k}}{{\partial \phi }}} \quad {\frac{{\partial f_p^k}}{{\partial r }}} \right. \\ &\left. {\frac{{\partial f_p^k}}{{\partial p }}} \quad {\frac{{\partial f_p^k}}{{\partial h }}} \quad {\frac{{\partial f_p^k}}{{\partial {X_0}}}} \quad {\frac{{\partial f_p^k}}{{\partial {Y_0}}}} \quad {\frac{{\partial f_p^k}}{{\partial {Z_0}}}} \right] \end{aligned} $$ (8) 系数矩阵
${{B}}$ :$$ {{B}}_j^k = \left[ {\begin{array}{*{20}{c}} {\dfrac{{\partial f_p^k}}{{\partial \alpha }}}&{\dfrac{{\partial f_p^k}}{{\partial \beta }}}&{\dfrac{{\partial f_p^k}}{{\partial \gamma }}} \end{array}} \right] $$ (9) 闭合差向量
${{W}}$ :$$ W_j^k = \frac{{{{\left( {a_p^j{x_k} + b_p^j{y_k} + c_p^j{z_k} - d_p^j} \right)}^2}}}{{{{\left( {a_p^j} \right)}^2} + {{\left( {b_p^j} \right)}^2} + {{\left( {c_p^j} \right)}^2}}} $$ (10) 采用最小二乘平差法求解,观测值、安置角改正数:
$$ {\hat{X}} = - {\left( {{{{B}}^{\rm{T}}}{{\left( {{{A}}{{{P}}^{ - 1}}{{{A}}^{\rm{T}}}} \right)}^{ - 1}}{{B}}} \right)^{ - 1}}{{{B}}^{\rm{T}}}{\left( {{{A}}{{{P}}^{ - 1}}{{{A}}^{\rm{T}}}} \right)^{ - 1}}{{W}} $$ (11) $$ {{V}} = - {{{P}}^{ - 1}}{{{A}}^{\rm{T}}}{\left( {{{A}}{{{P}}^{ - 1}}{{{A}}^{\rm{T}}}} \right)^{ - 1}}\left( {{{W}} + {{B\hat X}}} \right) $$ (12) 将改正数代入,迭代求解出角元素平差值:
$$ {\hat{L}} = {{L}} + {{V}} $$ (13) $$ {\hat{X}} = {{{X}}^0} + {\hat{X}} $$ (14) 当安置角求解完成后,将其作为已知量对平移参数进行求解。求解过程与安置角求解过程类似,最后迭代求出平移参数。
-
为了评估安置参数检校结果,使用车载移动测量系统采集了两组实验数据:第一组为检校场数据,用于验证检校方法和评估内符合精度;第二组为外场数据,区域内包括利用RTK、全站仪等手段采集的检核点以及少量靶标球,用于评估检校后移动测量系统的外符合精度。实验所用的移动测量系统为青岛秀山移动测量有限公司自主研发的VSurs-Q型轻便型移动测量系统,如图4所示。
-
为了最大化体现出安置误差影响,通常采集路线均以往返扫描作为一组并尽可能获取完整的参考面扫描点云数据。
在外业数据采集过程中,需进行必要的机动以保证组合导航精度。采集完成后,需要对组合导航处理结果进行精度检核,采用Waypoint Inertial Explorer 8.80软件的紧组合算法对基站和流动站数据进行处理。一般要求后处理结果要基本达到组合导航系统标称精度,避免引入其他误差造成检校结果较差。若POS处理精度较差,则重新安排检校实验。
-
外场数据采集区域如图5所示,采集区域往返路线长度约15 km,主要分为建筑物采集区和道路采集区两部分,区域内检核点主要采用全站仪结合CORS手段获取。建筑采集区包括广场、教学楼,包含大量的规则标志物,部分特征如图5(b)、(d)所示,因主要在广场区域采集,故GNSS观测条件较好;道路采集区包含道路标线、交通标志牌、防撞桶、栅栏等特征地物,考虑到存在特征点提取误差,数据采集时还布设了高精度标靶球作为检核之一,标靶球点位采用静态观测获得,部分特征如图5(a)、(c)所示,道路两侧分布有较高大的树木且部分区域建有高楼,比较接近一般工程作业环境。移动测量系统扫描点密度随着距离的增加而变小,故检核点主要分布在系60 m扫描半径范围内以保证准确提取特征点。
使用移动测量系统对外场采集区域进行往返测量,采集路线如图5中黄线所示。当采集完成后进行POS解算,后处理组合导航位置及姿态精度整体良好,高程精度优于0.055 m,平面精度优于0.025 m,航向角精度优于0.28 ',横滚角、俯仰角精度优于0.14 '。为方便比对,测得点云数据转换至投影坐标系下与检核点进行外符合精度评估。
-
车载移动测量数据采集完成,POS后处理精度满足要求后,开始进行安置参数检校处理,参数解算处理流程如图6所示,安置参数初始值可以通过系统设计模型中获得或者利用高精度工业测量系统进行标定。
首先,使用初始参数和组合导航处理结果对扫描仪原始数据进行融合变换,生成扫描点云数据和观测值索引文件;然后,利用RANSAC算法从每趟路线扫描的点云数据中提取出参考面对应的扫描面并进行标记;根据提取的扫描点云面的索引信息从观测值索引文件中检索出观测值信息,生成观测值信息文件;通过编号将每个扫描点云面与参考面进行匹配,获得约束方程;至此完成数据准备工作,开始进行参数解算,第一步构建安置角参数求解平差模型,利用最小二乘迭代求解安置角参数,第二步是将解算出来的安置角参数作为已知值,构建安置平移参数求解平差模型,迭代求解安置平移参数,最终输出安置参数解算结果并进行精度评定。计算过程与最终计算结果如表1所示。
表 1 安置参数解算结果
Table 1. Solution result of boresight parameters
ΔX/m ΔY/m ΔZ /m α/rad β/rad γ/rad Initial value 0.067 0.306 0.210 1.577 79 –1.041 09 0.022 07 Angle parameter Iteration1 0.067 0.306 0.210 1.581 73 –1.051 20 0.015 75 Iteration2 0.067 0.306 0.210 1.582 49 –1.049 65 0.016 67 Iteration3 0.067 0.306 0.210 1.582 40 –1.049 60 0.016 57 Iteration4 0.067 0.306 0.210 1.582 40 –1.049 61 0.016 56 Translation parameter Iteration1 0.069 0.306 0.210 1.582 40 –1.049 61 0.016 56 Iteration2 0.070 0.307 0.208 1.582 40 –1.049 61 0.016 56 Final results 0.070 0.307 0.208 1.582 40 –1.049 61 0.016 56 将计算结果与初始安置参数比对,检校前后α角相差–0.26°、β角相差0.48°、γ角相差0.31°,若不进行高精度检校,车载移动测量系统获取数据的精度会受到极大影响,将无法满足高精度测量任务。造成这种误差的主要原因是标定的传感器轴系与实际轴系存在偏差。这种偏差对角度影响较大,对平移量较小,通过平移参数改正量即可得到验证,检校后的平移参数改正量为毫米级,基本符合实际加工安装情况。从迭代过程看出,分步平差法在一定程度上分离了安置参数中的角元素和线元素,可避免解算过程中二者耦合带来其他不确定因素影响。
-
文中设计了一种内符合精度评估方法,通过将检校后的扫描点投影到匹配的参考特征面上,计算出扫描点和投影点之间的坐标和距离偏差来评估检校结果,如图7所示。图中所有数据均在地心地固坐标系下,C为对应参考面,红色箭头为平面法向量
${{n}}$ ,${P_i}$ 为激光脚点坐标,${P_s}$ 为${P_i}$ 对应的投影点坐标。理想情况下,参考面的点云数据应与参考面C重合,但由于存在多种随机误差和残留系统误差影响,检校后的激光扫描点仍未必完全在参考面C上。误差的计算方式为${\rm{d}}x = {X_i} - {X_S}$ ,${\rm{d}}y = {Y_i} - {Y_S}$ ,${\rm{d}}z = {Z_i} - {Z_S}$ ,${\rm{d}}s = \cos \left\langle {{{n}},\overrightarrow {{P_i}{P_s}} } \right\rangle$ $\sqrt {\rm{d}}{x^2} + {\rm{d}}{y^2} + {\rm{d}}{z^2}$ 。内符合精度评估结果如图8所示,参与计算的扫描点数共239 409个。图8(a)、(b)分别是检校前、检校后四个评估残差对应的频率分布直方图。由图8(a)知,安置误差检校前,距离偏差
${\rm{d}}s$ 分布于±0.170 m区间内,中误差为0.050 m,误差较为分散,内符合精度较差。由图8(b)知,安置误差检校后,距离偏差${\rm{d}}s$ 分布于±0.035 m区间内,中误差为0.007 m。误差分布集中,几乎接近正态分布,由此可以看出因安置参数误差存在而导致的系统差被基本消除,检校后的数据内符合精度较高。 -
在检校过程中,为了较好地消除安置误差,在作业场地选择、检校路线设计方面均采用较高标准,最大限度地保证检校过程中不引入其他误差。然而,在实际作业过程中,环境较为复杂,系统实际精度一般低于检校精度。为验证测量精度的可靠性,对系统外符合精度进行评估,选择泰山广场和前湾港路作为精度评估实验区域,通过比对同名特征点坐标差值的方式进行外符合精度评估。数据采集前,已预先通过CORS和全站仪获取部分特征点坐标。数据采集完成后,解算成投影坐标系下点云,在点云中拾取同名特征点进行比对,方向位偏差分布如图9所示。
此次参与外符合精度评估的球形靶标、建筑物角点、交通标志线角点共99个特征点,其中球形靶标27个,建筑物角点46个,交通标志线角点26个。按公式(15)计算,外符合精度为0.024 m。
$$ \sigma {\rm{ = }}\sqrt {\frac{{\sum\limits_{i = 1}^m {\left[ {{\rm{d}}x_i^2 + {\rm{d}}y_i^2 + {\rm{d}}z_i^2} \right]} }}{m}} $$ (15) 前湾港路为双向六车道,车载移动测量系统沿道路双向各扫描一次,往返扫描点云数据见图10,蓝色为自西向东扫描数据,红色为自东向西扫描数据,从水平和垂直方向分别选取特征地物往返扫描结果,图(a)为花坛中的树木,图(b)为道路中央栅栏,图(c)为交通标志牌,图(d)为道路横向剖面图。图中往返两次点云数据的重合度良好。由此可看出,利用文中方法检校后,车载移动测量系统的安置误差基本消除。
Boresight parameters calibration method of VMLS system based on reference planar features constraint
-
摘要: 车载移动测量系统是一种多传感器高度集成的测量设备,系统精度不仅取决于集成的传感器精度,还受激光扫描仪与组合导航系统之间安置参数检校的准确度影响。考虑到安置参数检校方法的便捷、有效性以及系统最终精度评估,提出一种基于参考面特征约束的车载移动测量系统安置参数检校方法。该方法根据包含系统安置参数的激光扫描点定位方程,利用参考面上的激光扫描点到参考面方程距离偏差最小作为约束条件,同时考虑到安置参数旋转量与偏移量间存在相关性,采用分步解算方法将旋转和平移量进行分开求解。最后,通过采集检校场和外场数据进行系统内符合和外符合精度评估。实验结果表明:该方法能够有效的消除安置误差影响,检校后内符合精度为0.007 m,外符合精度为0.024 m。Abstract: The vehicle mobile laser scanning (VMLS) system is a highly integrated multi-sensor measurement system. The accuracy of the VMLS system depends on the accuracy of not only the integrated sensor, but also the calibration of boresight parameters between the laser scanner and the integrated navigation system. Considering the convenience and effectiveness of the calibration method of the boresight parameters and the accuracy evaluation of the VMLS system, a calibrating method of the boresight parameters of the VMLS system based on the reference planar features constraint was proposed. The proposed method was based on the direct georeferencing of lidar measurements including the boresight parameters. The equation used the minimum distance deviation from laser footpoint to reference plane as the constraint. In addition, considering the correlation between rotation and offset amount of the boresight parameters, the stepwise solving method was composed to separate boresight angle and offset parameters. Finally, through the datum acquired from calibration field and check fields, the internal and external accuracy of VMLS system were evaluated. The experimental results show that the proposed method can effectively eliminate the influence of boresight errors. After calibration, the internal accuracy is 0.007 m and the external accuracy is 0.024 m.
-
表 1 安置参数解算结果
Table 1. Solution result of boresight parameters
ΔX/m ΔY/m ΔZ /m α/rad β/rad γ/rad Initial value 0.067 0.306 0.210 1.577 79 –1.041 09 0.022 07 Angle parameter Iteration1 0.067 0.306 0.210 1.581 73 –1.051 20 0.015 75 Iteration2 0.067 0.306 0.210 1.582 49 –1.049 65 0.016 67 Iteration3 0.067 0.306 0.210 1.582 40 –1.049 60 0.016 57 Iteration4 0.067 0.306 0.210 1.582 40 –1.049 61 0.016 56 Translation parameter Iteration1 0.069 0.306 0.210 1.582 40 –1.049 61 0.016 56 Iteration2 0.070 0.307 0.208 1.582 40 –1.049 61 0.016 56 Final results 0.070 0.307 0.208 1.582 40 –1.049 61 0.016 56 -
[1] Yang Bisheng, Liang Fuxun, Huang Ronggang. Progress, challenges and perspectives of 3D LiDAR point cloud processing [J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1509−1516. (in Chinese) [2] Wang Yu, Cheng Liang, Chen Yanming, et al. Building point detection from vehicle-borne LiDAR data based on voxel group and horizontal hollow analysis [J]. Remote Sensing, 2016, 8(5): 419. doi: 10.3390/rs8050419 [3] Yao Lianbi, Wang Zhifei, Sun Haili. Design and implementation of vehicle laser scanner's external parameter calibration [J]. Journal of Tongji University (Natural Science), 2016, 44(1): 161−166. (in Chinese) [4] Glennie C, Lichti D D. Static calibration and analysis of the velodyne HDL-64E S2 for high accuracy mobile scanning [J]. Remote Sensing, 2010, 2(6): 1610−1624. doi: 10.3390/rs2061610 [5] Habib A, Kersting A, Ruifang Z, et al. LiDAR strip adjustment using conjugate linear features in overlapping strips[D]. Beijing:Proceedings of International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, 2008: 385-390. [6] Habib A, Bang K, Kersting A. Impact of lidar system calibration on the relative and absolute accuracy of derived point cloud[C]//The Euro COW 2010, 2010, 2: 283. [7] Maas H G. Least-squares matching with airborne laserscanning data in a TIN structure [J]. International Archives of Photogrammetry and Remote Sensing, 2000, 33(3): 548−555. [8] Lee J, Yu K, Kim Y, et al. Adjustment of discrepancies between LIDAR data strips using linear features [J]. IEEE Geoscience and Remote Sensing Letters, 2007, 4(3): 475−479. doi: 10.1109/LGRS.2007.898079 [9] Zhang Ka, Sheng Yehua, Ye Chun, et al. Absolute calibration and precision analysis for vehicle-borne 3D data acquiring system integrated with GPS, INS and CCD-camera [J]. Geomatics and Information Science of Wuhan University, 2008, 33(1): 55−59. (in Chinese) [10] Nie Qian, Chen Weimin, Chen Changjun. Research on absolute calibration for vehicle-borne 3D laser imaging system [J]. Bulletin of Surveying and Mapping, 2013, 11: 74−77. (in Chinese) [11] Wang Li, Li Guangyun, Yang Yongxin, et al. One-site calibration of mounted parameter method for 3D mobile laser scanning [J]. Infrared and Laser Engineering, 2016, 45(11): 1106005. (in Chinese) doi: 10.3788/IRLA201645.1106005 [12] Kang Yongwei, Zhong Ruofei, Wu Yu. Research of calibrating vehicle laser scanner′s external parameters [J]. Infrared and Laser Engineering, 2008, 37(S): 249−253. (in Chinese) [13] Ye Zetian, Yang Yong, Zhao Wenji, et al. Dynamic calibration of exterior orientations for vehicle GPS/ IMU/ LS laser imaging system [J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(3): 345−350, 372. (in Chinese) [14] Pusztai Z, Eichhardt I, Hajder L. Accurate calibration of multi-lidar-multi-camera systems [J]. Sensors, 2018, 18(7): 2139. doi: 10.3390/s18072139 [15] Jian Z, Ravani B. Boresight calibration of mobile laser scanner using an external fixture [J]. Journal of Applied Geodesy, 2016, 10(3): 0011. [16] Zhang Jing, Jiang Wanshou, Jiang San. Automated airborne lidar system calibration using virtual tie point model [J]. Acta Geodaetica et Cartographica Sinica, 2013, 42(3): 389−396. (in Chinese) [17] Li Z, Tan J, Liu H. Rigorous boresight self-calibration of mobile and UAV lidar scanning systems by strip adjustment [J]. Remote Sensing, 2019, 11(4): 442. doi: 10.3390/rs11040442 [18] Chan T O, Lichti D D, Glennie C L. Multi-feature based boresight self-calibration of a terrestrial mobile mapping system [J]. Isprs Journal of Photogrammetry and Remote Sensing, 2013, 82: 112−124. doi: 10.1016/j.isprsjprs.2013.04.005 [19] Heinz E, Eling C, Wieland M, et al. Analysis of different reference plane setups for the calibration of a mobile laser scanning system [J]. Proceedings of the Internationalen Ingenieurvermessungskurs, 2017: 25−29. [20] Tian Xiangrui, Xu Lijun, Xu Teng, et al. Calibration of installation angles for mobile lidar scanner system [J]. Infrared and Laser Engineering, 2014, 43(10): 3292−3297. (in Chinese) [21] Lu X, Feng C, Ma Y, et al. Calibration method of rotation and displacement systematic errors for ship-borne mobile surveying systems [J]. Survey Review, 2017, 51(364): 78−86. [22] Le Scouarnec R, Touzé T, Lacambre J B, et al. A new reliable boresight calibration method for mobile laser scanning applications [J]. ISPRS - International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, 2014, XL-3/W1: 67−72. doi: 10.5194/isprsarchives-XL-3-W1-67-2014 [23] Rieger P, Studnicka N, Pfennigbauer M, et al. Boresight alignment method for mobile laser scanning systems [J]. Journal of Applied Geodesy, 2010, 4(1): 13−21. [24] Zhang Haixiao, Zhong Ruofei, Sun Haili. External parameter calibration method of vehicle laser scanning system based on planar features [J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(12): 1640−1649. (in Chinese) [25] Skaloud J, Lichti D. Rigorous approach to bore-sight self-calibration in airborne laser scanning [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2006, 61(1): 47−59. doi: 10.1016/j.isprsjprs.2006.07.003