基于同名特征地物的移动测量系统激光扫描仪标定方法
【技术领域】
[0001] 本发明涉及一种基于同名特征地物的车载移动测量系统激光扫描仪标定方法，属 于测绘信息领域。
【背景技术】
[0002] 移动测量技术是当今测绘界最为前沿的技术之一，移动测量系统成为数字城市、 传统测绘等领域获取数据的新手段，具有十分广泛的应用。
[0003] 移动测量系统采集的数据主要包括POS数据，激光点云数据和影像数据。通过激 光点云数据可以直接获得地物坐标及特征信息，进而进行DSM生成、三维建模等后续生产。 但获得高精度的点云数据的前提是激光扫描仪已经进行了精确标定。
[0004] 传统的激光扫描仪标定都需要特定的标定场地，在实际应用中有很大限制，降低 了数据采集的效率，不利于移动测量系统的推广和应用。另外，从标定场到作业区域的运输 过程，也有可能使激光扫描仪的实际参数与标定结果产生一定的误差，从而降低了激光点 云的精度。
【发明内容】
[0005] 本发明提供了一种基于同名特征地物的移动测量系统激光扫描仪的标定方法，解 决了上述【背景技术】中的不足，大大提高了移动测量系统数据采集的效率，有效的推动了移 动测量系统的应用，具有非常巨大的经济和社会价值。
[0006] 实现本发明上述目的所采用的技术方案为：
[0007] -种基于同名特征地物的移动测量系统激光扫描仪的标定方法，其特征在于包括 以下步骤：
[0008] (1)选取点云中同名点、线、面特征地物对应的激光点；
[0009] (2)获得这些激光点对应的源数据；
[0010] (3)获取激光扫描仪待标定参数的初始值；
[0011] (4)根据不同类型的特征地物，选取不同的标定方程，采用最小二乘平差算法，计 算激光扫描仪的最优参数，完成标定。
[0012] 步骤（1)中同名特征地物对应的激光点获取方法为人工选取；
[0013] 步骤（2)中获取这些激光点对应的源数据的方法为：根据激光点的三维坐标， 提取激光点对应的源数据，包括激光点发射瞬间POS的位置和姿态（posX、posY、posZ、 Heading、Pitch、Roll)以及激光点在激光扫描仪坐标系下的坐标（scanX、scanY、scanZ)。
[0014] 步骤（3)中的扫描仪待标定参数初始值一般为扫描仪设计和安装的理论值。
[0015] 步骤（4)中，根据同名特征地物具有相同特性，列出标定方程，建立参数解算模 型。
[0016] 假设激光头扫描坐标系在载体坐标系中的位置为T (XQ，YQ，Z。)τ，姿态为（Heading， Pitch，Roll)，其旋转矩阵为R。同名点、线、面标定方程如下：
[0017] a)同名点标定方程
[0018] 同名点和在大地坐标系中坐标应相同，即
[0019] Rpl (R*Prawl+T) +Tpl= R p2 (R*Praw2+T) +Tp2 (I. I)
[0020] 其中：
[0021] Rpl、Rp2:为两次采集时，载体坐标系在大地坐标系中的姿态角所构成的
[0022] 旋转矩阵；
[0023] Tpl、Tp2为两次采集时，载体坐标系在大地坐标系中的位置；
[0024] Prawl、Pmw2:为两次采集时，同名点在激光扫描仪坐标系下的坐标；
[0025] 误差方程为
[0027] 根据I. 1~1. 2可列误差方程，约束条件为：每对同名点坐标差最小。η对同名点， 总共列3*η个方程，未知数R、T共6个，故至少两对控制点可解算未知数。
[0028] b)同名线标定方程
[0029] 设同名线有η条，第i条线的方程为
[0031] 每个点的世界坐标系坐标为
[0032] P (Xp，Yp，Zp) = Rp (R*Praw+T) +Tp (1.4)
[0033] 该点在对应直线上的垂足为Qfei+mjtq，bi+njtq，Ci+pjtq)。
[0034] Q点在过P点且与对应直线垂直的平面上，可求出tq:
[0038] 根据1. 4~1. 6可列误差方程，约束条件为：每个点距离其对应的直线距离最小。 η条直线，每个线特征取m(m> 1)个激光点观测数据，总共列n*m个方程，未知数为
[0039]
[0040] 共6+6*n个，故n*m > 6*n+6时，方可计算出未知数。
[0041] c)同名面标定方程
[0042] 设同名面有η个，第i个面的方程为
[0043] aiX+biY+CiZ+di= 0
[0044] 其中 i = 1，2, · · ·，η。向量
[0046] 为面的单位法向量，即
[0048] 每个点的世界坐标系坐标为
[0049] Pw (X，Y，Z) = Rp (R*Praw+T) +Tp (1.8)
[0050] 该点与平面的距离为
[0051] V = aiX+biY+CiZ+di (1. 9)
[0052] 根据1. 7~1. 9可列误差方程，约束条件为：每个点距离其对应的平面距离最小。 η个平面，每个平面取m(m > 2)个激光点观测数据，总共列n*m个方程，未知数为
[0053]
[0054] 共6+4*n个，故n*m > 4*n+6时，方可计算出未知数。
[0055] 根据以上方程，采用最小二乘算法，求解出激光扫描仪待标定的最优参数。
[0056] 本发明所提供的方法，利用作业区域现有的点、线、面特征，直接根据点与点、线与 线、面与面间的关系，建立模型，解算出激光扫描仪待标定的参数。相对于传统的采用标定 场对激光扫描仪进行标定的方法，本发明提供的方法不需要专门的标定场，有特征点面地 物的地方皆可标定，节省了专门去标定场标定的时间，克服了需要标定却没有标定场的难 题，提高了移动测量系统作业效率，大力推广了移动测量系统的应用，具有很高的实际应用 价值。
【附图说明】
[0057] 图1 :为基于同名点特征地物的车载移动测量系统激光扫描仪标定流程图。
[0058] 图2 :为基于同名线特征地物的车载移动测量系统激光扫描仪标定流程图。
[0059] 图3 :为基于同名面特征地物的车载移动测量系统激光扫描仪标定流程图。
【具体实施方式】
[0060] 下面结合具体实施例对本发明做详细具体的说明。
[0061] 1.基于同名点特征地物的车载移动测量系统激光扫描仪标定方法
[0062] 本实施例中的同名点特征为道