多视点云数据同步配准新方法

朱瑞芳，方 勇

(上海大学通信与信息工程学院，上海 200071)

0 引言

目前，三维激光扫描技术的测量能力、速度、自动化程度、数据处理效率以及整体经济效益均明显高于其他测量技术。由于被测目标物体表面的互相遮挡以及测量设备视域的局限，要得到物体表面完整的几何信息，需从多个角度对物体进行测量。多视扫描导致所得点云数据没有相同的参照系，配准的目的正是将其转换到统一的坐标系下。所以说点云数据的配准是真实世界物体三维模型重建过程中非常重要的环节，数据匹配的精度直接制约着激光扫描技术在工程中的推广应用。

对于三维数据的配准方案，许多学者做了大量的研究工作。最具代表性的是Besl等提出的迭代最近点算法［1］(iterative closest point，ICP)。该算法利用牛顿迭代原理寻找2组点云对应的最近点对，再利用点对的欧氏距离作为目标函数迭代求得三维的刚体变换。ICP算法对参与配准的点云有较严格的限制，如2组点云之间的初始位置不能相差太大，要求其中一个点集是另一个点集的子集等。袁建英［2］、戴静兰［3］等对 ICP 算法做了改进，虽然配准精度很高，但海量点云数据的迭代计算速度较慢［4］。

针对以上局限性，本文提出基于空间相似变换的多视点三维激光扫描点云数据的同步配准方法。通过引入空间相似变换原理，运用罗德里格旋转矩阵实现点云数据的同步配准，利用Real works Survey软件研究同步配准方案的准确性和可靠性。

1 空间相似变换原理

坐标系之间的转换关键就是要求得到高精度的布尔莎7参数［5］，包括3个平移参数、3个旋转参数和1个比例参数。设有2个空间直角坐标系分别为O－XYZ和O'－X'Y'Z'(图1)，其坐标系原点不一致，存在3个平移参数ΔX，ΔY，ΔZ;它们间的坐标轴也相互不平行，存在3 个旋转参数 εX，εY，εZ。

图1 不同空间直角坐标系Fig．1 Different spatial rectangular coordinate system s

显然，2个坐标系需通过坐标轴的平移和旋转，变换至同一坐标系下，坐标系的转换关系如下:

在3D空间下，设2个不同坐标系间的对应点坐标分别为X，Y，Z和X'，Y'，Z'，则

式中:ΔX，ΔY和ΔZ为平移参数;λ为比例参数;R为正交旋转矩阵，即

式中:r11=cosεYcosεZ;r12=－cosεYsinεZ;r13=－sinεY;r21=cosεXsinεZ－sinεXsinεYcosεZ;r22=cosεXcosεZ+sinεXsinεYsinεZ;r23=－sinεXcosεY;r31=sinεXsinεZ+cosεXsinεYcosεZ;r32=sinεXcosεZ－cosεXsinεYsinεZ;r33=cosεXsinεY;εX，εY，εZ分别为绕X，Y，Z轴旋转的旋转参数。

若已知式(1)7个参数，利用式(1)就可方便地进行2个坐标系下的坐标转换。多数情况下，7个参数的值，至少需要7个独立方程根据式(1)反算得到。因此，首先对式(1)进行泰勒级数展开，使其线性化［6－8］，化简后得到观测方程

式中:δx，δy，δz为原点重心化坐标;λ1，λ2，λ3为线性化后比例参数系数;ann(1≤n≤13)为旋转角度参数。将式(3)写成间接平差形式，按平差解算得

式中:V为公共点对配准误差;A为系数阵;x=(ATA)－1(ATL);L为由常数项组成的矩阵。

在空间相似变换计算过程中，采用线性化公式需迭代过程。在式(1)中有7个未知数，至少需3对及以上不在同一直线上的控制点(2个坐标系中的公共点)采用间接平差法才可求出坐标转换参数。

2 多视点同步配准实现

多视点三维激光扫描数据同步配准流程如图2所示。

图2 多视点三维激光扫描数据同步配准流程Fig．2 Flow chart of multi－view 3D point cloud data synchronization registration

将空间相似变换结果作为初始参数值进行计算，根据各测站间的联系，对变换的初始参数进行联合平差，使各测站间的配准更精确，减少误差积累。

3 计算实例及分析

数据来源于佛像1∶10前期制作模型(图3)，扫描仪器为Trimble GXTM型地面三维激光扫描仪。

图3 原始佛像Fig．3 Original buddha

扫描时采用标定物来进行后期的数据配准工作，外置标定物为扫描仪系统自带的标靶。为了确定标靶数据在采集过程中空间三维坐标值保持不变，首先在佛像周围合适的位置建立一个坐标控制网，然后将标靶放至控制点中。本次数据采集在扫描物体周围设立了4个控制点，并记录其坐标值。

图4是佛像的点云数据，所扫描的4测站点云以不同颜色显示。

图4 佛像的点云数据Fig．4 Point cloud data of buddha

如图5所示，在控制网上共设立4个标靶，依次命名为A，B，C，D，拟合出的标靶中心坐标值为测站配准的原始数据，通过这些数据将12个测站的点云数据旋转变换至统一坐标系下，达到点云数据的同步配准目的。提取4个测站(前Q，后下HX，左上ZS，右上YS)坐标作为本次实验扫描的原始数据，用以检验配准方法的实用性和可靠性。

图5 4测站扫描示意图Fig．5 Schematic diagram of scanning by means of four survey stations

3．1 利用Realworks Survey软件进行配准

Realworks Survey为Trimble GXTM型地面三维激光扫描仪系统自带的后处理软件，该软件的配准模式中支持3种配准方式:基于点云的配准、基于目标的配准和通过大地参考的配准［9］。基于点云的配准是指在2幅扫描图像中，手动找出3个同名点，软件将自动进行搜索计算，最终实现2幅点云的配准;基于目标的配准是指在2幅或者多幅点云图像中，任意2幅点云有2个或2个以上同名标靶，将标靶作为目标进行标靶匹配，软件将自动将多幅点云配准到一起;大地参考的配准是先将标靶的大地坐标输入到软件中，然后通过标靶将点云数据配准到大地坐标系中(图6)。

图6 Realworks Survey配准示意图Fig．6 Schematic diagram of registration in Realworks Survey

本次实验选取基于目标的配准工具方案(图6)，以测站Q坐标数据作为本次实验的基准测站，将测站YS，HX和ZS坐标数据作为配准测站，自动配准至基准测站坐标系下，配准结果如图7所示。

图7 Realworks Survey点云配准效果Fig．7 Registration effect of point cloud by Realworks Survey

经Realworks Survey软件配准后，自动生成配准报告，根据测站和目标的配准报告，从中读出配准过程的基本信息及配准精度指标(表1)。

表1 Realworks Survey软件配准报告Tab．1 Registrated report of the soft wear Realwork Survey

3．2 本文方法配准

首先对实验数据进行预处理，然后提取标靶中心三维坐标值，最后对目标点云数据进行配准。2测站之间变换至少需要3对公共点作为配准的起始数据，才能求解出变换参数。本次实验数据中均采用4个标靶作为配准的公共点对，这样在2测站之间的配准过程就存在一个平差过程。针对本次实验数据，设计了基于4组公共点对的平差程序，并且根据参数初始值的求取原理，进一步编写程序，大角度空间直角坐标系转换参数求取界面如图8所示。

图8 大角度空间直角坐标系转换参数求取软件界面Fig．8 Interface of transform parameters for wide－angle spatial rectangular coordination system

其中第1，2，3，4点的公共点坐标分别对应为A，B，C，D标靶拟合中心点坐标值;k为比例参数;ΔX，ΔY，ΔZ为平移参数;an，bn，cn为罗德里格矩阵组成的旋转矩阵R的参数。根据自编程序可求出各测站旋转变换至基准参考系下(本实验以Q测站坐标系作为基准坐标系)的变换参数(表2)。

根据表2的旋转矩阵R，可得出罗德里格矩阵的旋转参数a，b，c(表3)。

表2 转换参数Tab．2 Transform parameters (m)

表3 罗德里格矩阵参数Tab．3 Rodriguez matrix parameters (m)

根据罗德里格旋转参数、平移参数和比例参数值，经过变换可求出初始变换后的标靶坐标中心三维坐标值(表4)，本文方法配准结果如图9所示。

图9 本文方法点云配准效果Fig．9 Registration effect of point cloud by the proposed method in this paper

表4 配准后标靶中心坐标值Tab．4 Coordinate values of target center points after registration (mm)

3．3 精度对比及分析

运用归一化残差平方和、最大残差、残差均方差及单位权中误差作为精度评定的标准，比较2种配准方案的精度及可靠性(表5)。

表5 同步配准前后残差值Tab．5 Residual error before and after synchronous registration

通过基于Realworks Survy软件配准和本文方法配准这2种方案的精度指标对比，可以得到:与软件配准相比，本文方法在某些方面还有待提高，但也相差无几，这2种方案的精度大体相当，说明本文方法具有可行性;在实际应用中，本文方法的收敛速度快，有较好的应用效果。

4 结论

1)基于空间相似变换的多视点云同步配准算法具有理论严密性，而且证明了三维激光扫描多视点云数据同步配准方案的可行性和实用性，应用中有良好的配准精度，在实例中评价指标与软件配准精度相当。Realworks Survey是优秀的商业软件，编写严密，经过了较长时间的实践和改进，可以作为其他算法的评判标准。本文算法满足应用精度，是一个很好的尝试。

2)本文方法计算时没有考虑标靶中心点拟合误差的影响，在通过Matlab软件计算误差方程过程中没有足够的迭代次数导致精度稍逊，是本文需要改进之处，另外如何将标靶的拟合误差影响加入到同步配准模型中，也是本文下一步需做的工作。

［1］ Besl P J，Mckay N D．A method for registration of 3－D shapes［J］．IEEE Transactions on Pattern Analysis and Machine Intelligence，1992，14(2):239－256．

［2］ 袁建英，刘先勇，刘 伟，等．改进ICP算法实现多视点云精确配准研究［J］．传感器与微系统，2008，27(5):27－30．Yuan JY，Liu X Y，Liu W，et al．Research on precise registration of multi－views based on improved ICP algorithm［J］．Transducer and Microsystem Technology，2008，27(5):27－30．

［3］ 戴静兰，陈志扬，叶修梓．ICP算法在点云配准中的应用［J］．中国图形图象学报，2007，12(3):517－521．Dai JL，Chen Z Y，Ye X Z．The application of ICP algorithm in point cloud alignment［J］．Journal of Image and Graphics，2007，12(3):517－521．

［4］ 徐源强，高井祥，张 丽，等．地面三维激光扫描的点云配准误差研究［J］．大地测量与地球动力学，2011，31(2):129－132．Xu Y Q，Gao JX，Zhang L，etal．Research on point cloud registration error of terrestrial laser scanning［J］．Journal of Geodesy and Geodynamics，2011，31(2):129－132．

［5］ 楚振庭，赵 健．国土资源工作中测绘数据坐标转换研究［J］．现代测绘，2010，23(1):20－23．Chu Z T，Zhao J．Research on conversion of survey data coordinates in land resource work［J］．Modern Surveying and Mapping，2010，23(1):20－23．

［6］ 周拥军，寇建新．空间相似变换的直接解法及比较研究［J］．工程勘察，2010(11):62－66，74．Zhou Y J，Kou JX．Comparative study on closed－form solutions for spatial similar transformation［J］．Geotechnical Investigation and Survey，2010(11):62－66，74．

［7］ 张 东，黄 腾，陈建华，等．基于罗德里格矩阵的三维激光扫描点云配准算法［J］．测绘科学，2012，37(1):156－157，173．Zhang D，Huang T，Chen J H，et al．An algorithm for terrestrial laser scanning point clouds registration based on rodriguez matrix［J］．Science of Surveying and Mapping，2012，37(1):156－157，173．

［8］ 张 卡，张道俊，盛业华，等．三维坐标转换的两种方法及其比较研究［J］．数学的实践与认识，2008，38(23):121－128．Zhang K，Zhang D J，Sheng Y H，etal．Research on two methods of three dimensional coordinate transformation and their comparison［J］．Mathematics in Practice and Theory，2008，38(23):121－128．

［9］ Trimble RealWorks Survey Technical Notes［OL］．Trimble navigation limited，2009．http://trl.trimble.com/docushare/dsweb/Get/Document－429608/022504－107B_PPP_RealWorks_TN_1109_sec．pdf#xml=http://tknbe．trimble．com:8080/pdf?id=9f6156c3ac＆pg=0＆iphl=based:target:tool:－:registrations:tooling:registration:tools．