地面三维激光扫描仪系统误差模型研究及精度分析
来源：《测绘地理信息》2016年10月
作者：蔡建民、花向红、宣伟、吴帮、杨荣华
分析了外部观测条件对地面三维激光扫描仪测距的影响，以空间相似变换为基础，提出了自检校法的优化模型，并给出了模型解算策略。通过实验对RIEGL VZ-400 地面三维激光扫描仪进行检校和测量数据改正，结果表明：将外部观测条件影响因素纳入系统误差模型中后，可以提高地面三维激光扫描仪的数据精度，从而验证了优化模型的可行性与正确性。
地面三维激光测量技术突破了传统单点测量的方法，通过扫描可获得每秒数万点的高精度点云数据[1]。随着这项技术在科学研究、工业生产等领域的广泛应用，工程项目对其测量精度也提出了较高的要求[2]。地面三维激光扫描仪的系统参数是由厂家在理想状况下测量并给定，受仪器构造及生产过程的不确定因素影响，导致仪器实际指标与标称不一致，给测量带来误差。为了快速获取系统误差的改正参数，并利用这些参数对测量数据进行系统性偏差改正，可对三维激光扫描仪系统误差进行建模和标定分析。张毅等[3]提出的11 参数自检校法未考虑外部观测条件因素对激光测距的影响，本文在地面三维激光扫描仪系统误差影响因素的基础上，对激光测距误差进行深入研究，提出了11参数自检校法的优化模型，利用稳定性好、精度高和计算速度快的罗德里格矩阵求解模型参数[4]，详细地给出了模型解算策略。
1系统误差模型
1.1 自检校法
自检校法是以空间相似变换为基础，由11个待定参数构成地面三维激光扫描系统误差模型[5]，其中6个为外部定向参数，另外5 个为仪器系统误差参数。式(1)为自检校法模型：
式中，k为三维激光测距加常数改正量;r为三维激光测距乘常数改正系数[4]。当k=0时，地面三维激光扫描仪为相位式测距;当r=0时，则为脉冲式测距。R(Φ，Ω，Κ)为旋转矩阵;
该模型对测距S进行了加常数和乘常数改正，同时也考虑了视准轴误差和水平轴误差对测角的影响，通过改正后的三个间接观测值重新计算扫描点坐标。但该模型对距离S进行误差改正时，只考虑了激光测距方式[7]所引起的误差，未考虑测量时外部观测条件对激光测距的影响。
1.2 模型优化
地面三维激光扫描仪的测距误差受激光测距方式和外部观测条件影响，对外部观测条件而言其误差主要是由非线性状态、时间漂移、温度漂移、混合像素等因素引起的[8]，可表示为：
式中，ΔSinst为测距误差;t为温度;T为地面三维激光扫描仪开机时长;y为扫描仪的使用年份;ΔSwalk为时间漂移误差;m0、m1、m2、m3 为尺度因子;b为时间效应比例因子;c和c1 为测距和温度比例因子;最后4项为非线性距离制约误差;其余参数为扫描仪校准系数。
本文将ΔSinst作为一项改正参数，对激光测距误差作进一步改正，式(3)为优化后的自检校法模型：
式中，k=k′+ΔSinst。
测距误差一部分为激光测距方式引起的，另一部分为环境因素的影响。考虑到模型中测距误差ΔSinst和k之间存在关联性，将二者作为整体，即加常数k和环境因素ΔSinst作为一个独立的参数k。
1.3 模型解算
利用Taylor级数将式(3)展开，保留一次项，ΔSinst和加常数合并后的线性展开式中含有11个待定参数。
对待定参数分别求偏导数得：
对旋转矩阵求偏导数得：
将以上所求的偏导数构成系数矩阵B，每个偏导数为3×1的子矩阵，则系数矩阵B从标靶1到标靶n可表示为：
然后计算出误差方程常数项L，每个标靶点的误差方程常数项为3×1的子矩阵，常数项矩阵L从标靶1到标靶n可表示为：
假设实验中共有n对同名标靶点参与自检校模型的解算，则共有3n个误差方程，利用间接平差可得待定参数的改正数为：
式(3)的线性展开式中含有自检校模型的待定参数，因此需确定待定参数的初值并采用迭代法进行计算。对于初值问题，优化后的自检校模型中有6个外部定向参数，假设三维激光扫描仪无系统误差，在此条件下由式(1)对同名点进行空间坐标转换，获取旋转矩阵
式中，x0为待定参数初始值;x1、x2、…、xn为迭代计算时待定参数的改正数。
利用式(9)得到各参数的最终值，结合式(3)计算经系统误差改正后的扫描点坐标，并将改正后的扫描点坐标通过空间相似变换转换至全站仪的绝对坐标系下，同全站仪所测量的同名点标靶坐标进行比较和分析。
2实验设计与结果分析
为了验证优化后的自检校法模型的可靠性与正确性，以RIEGL VZ-400 地面三维激光扫描仪(测距精度±5mm，测角精度±1.8″)为实验对象，对其系统误差进行研究和分析。在实验场地中均匀有层次地布设了10个标靶并对其进行编号。首先使用三维激光扫描仪对实验场地进行扫描，并对标靶点进行精扫描，利用该扫描仪的配套软件RISCANPRO 对数据进行处理，获取标靶区域的点云数据，根据点云数据的强度信息采用阈值分割的方法得到标靶区域，对标靶区域点云数据进行强度加权拟合，获取标靶几何中心在三维激光扫描仪下的坐标。然后采用测量精度优于RIEGL VZ-400 的LeicaTCRP1201全站仪(测距精度±1 mm，测角精度为±1″)对准反射贴片的十字中心，获取同名点标靶在独立坐标系下的几何中心坐标，并以此作为绝对坐标对扫描仪进行系统误差检校。表1为扫描仪和全站仪获取的三维坐标。
表1 标靶在扫描仪和全站仪下的坐标
不考虑三维激光扫描仪内部系统误差时，将表1中的A1-A5作为模型参数计算点，利用其标靶坐标可求出三维激光扫描仪内部坐标系转换至全站仪坐标系的缩放因子、旋转矩阵和平移参数。缩放因子λ 为1.000089，旋转矩阵为：
平移参数为：
通过以上的7 个参数(缩放因子λ 未参与迭代计算)，将三维激光扫描仪的坐标全部转换至全站仪的绝对坐标系下，并与全站仪测量的标靶坐标进行比较，结果见表2，其中A6-A10为模型参数检核点。
表2 扫描仪转换后坐标与全站仪坐标之差比较
考虑地面三维激光扫描仪系统误差后，利用A1-A5号标靶坐标并根据式(8)可计算出各系统误差参数。实际测量时环境因素对激光测距的影响小于激光测距方式的影响，根据实验数据的多次组合，当二者权值为1∶3 时效果较为理想，从而获得ΔSinst的具体数值，结果见表3。
利用表3中的模型参数通过式(3)对三维激光扫描仪测量的A1-A10 标靶坐标进行改正，并利用空间相似变换将改正后的标靶坐标转换至全站仪坐标系下，与全站仪测量的标靶坐标进行比较和分析，结果见表4。
表4 扫描仪经系统误差参数改正转换后坐标与全站仪坐标之差比较
由表2和表4 可以看出：地面三维激光扫描仪经系统误差参数改正转换后坐标与全站仪坐标之差与未经系统误差参数改正转换后坐标与全站仪坐标之差相比，数值上各标靶点的X、Y、Z3个方向明显减少。在不考虑全站仪坐标测量误差时，根据表2的数据计算，对于模型参数计算点，X、Y、Z方向中误差和点位中误差分别为0.0049 m、0.0048 m、0.0050m和0.0084 m;对于模型参数检核点，X、Y、Z方向中误差和点位中误差分别为0.0069 m、0.0072m、0.0073m 和0.0123m。根据表4数据计算，对于模型参数计算点，X、Y、Z方向中误差和点位中误差分别为0.0020m、0.0015m、0.0030m和0.0040m;对于模型参数检核点，X、Y、Z方向中误差和点位中误差分别为0.0035 m、0.0042 m、0.0040m和0.0068m。由此可见，模型参数检核点经系统误差改正后点位精度提高约45%，反映出三维激光扫描仪误差的系统性，也证明了优化后的系统误差模型可以提高扫描仪测量数据精度。
3结束语
本文研究了三维激光扫描仪测距时系统误差的来源，分析了系统误差对激光测距的影响。以地面三维激光扫描仪的观测方程为载体对11 参数自检校法提出了优化模型，引入外部条件综合影响因子ΔSinst，并对激光测距进行改正。详细推导了公式，并通过实验验证了优化模型的正确性与可靠性。从实验结果可以看出，经过系统误差改正后，RIEGLVZ-400地面三维激光扫描仪的数据质量有了很大的提高，从而验证了系统误差模型的可行性与正确性，为促进地面三维激光扫描仪系统误差检校体系的发展和完善提供了参考。
推荐阅读
---------------- 推 广 ---------------