Ground-based LiDAR measurement system integrating BeiDou/GNSS positioning and 5G communication

融合北斗/GNSS定位和5G通讯的地基激光雷达测量系统 Ground-based LiDAR measurement system integrating BeiDou/GNSS positioning and 5G communication 1 2 3 4 first-author 郭 明（1981-），男，湖北黄陂人，博士，教授，博士生导师，2004 年于江汉大学获得学士学位，2008 年于北京建筑工程学院获得硕士学位，2011 年于武汉大学获得博士学位，主要从事激光雷达技术应用及三维地理信息系统方面的研究，研究方向为文化遗产数字化保护、激光雷达测绘技术、移动测量技术、建筑健康监测等。E-mail： guoming@bucea.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436305&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436302&type= 22.01332855 32.00399780 郭 明 （1981-），男，湖北黄陂人，博士，教授，博士生导师，2004 年于江汉大学获得学士学位，2008 年于北京建筑工程学院获得硕士学位，2011 年于武汉大学获得博士学位，主要从事激光雷达技术应用及三维地理信息系统方面的研究，研究方向为文化遗产数字化保护、激光雷达测绘技术、移动测量技术、建筑健康监测等。E-mail： guoming@bucea.edu.cn guoming@bucea.edu.cn 1 5 corresp E-mail： 1033547194@qq.com E-mail： 1033547194@qq.com 郭可才（1983-），男，湖北武汉人，高级工程师，2006年毕业于鄂州大学，现为北京申信达成科技有限公司总经理，主要从事移动LiDAR扫描、工程测量和自动化技术方面的研究。E-mail： 1033547194@qq.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436308&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436312&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436309&type= 22.01333237 32.00400543 郭可才 （1983-），男，湖北武汉人，高级工程师，2006年毕业于鄂州大学，现为北京申信达成科技有限公司总经理，主要从事移动LiDAR扫描、工程测量和自动化技术方面的研究。E-mail： 1033547194@qq.com 1033547194@qq.com 1 1 北京建筑大学 测绘与城市空间信息学院，北京 102616 北京建筑大学 测绘与城市空间信息学院，北京 102616 School of Geomatics and Urban Spatial Informatics， Beijing University of Civil Engineering and Architecture， Beijing 102616， China School of Geomatics and Urban Spatial Informatics， Beijing University of Civil Engineering and Architecture， Beijing 102616， China 代表性建筑与古建筑数据库教育部工程研究中心，北京 100044 代表性建筑与古建筑数据库教育部工程研究中心，北京 100044 Engineering Research Centre of Representative Building and Architectural Heritage Database， Ministry of Education， Beijing 100044， China Engineering Research Centre of Representative Building and Architectural Heritage Database， Ministry of Education， Beijing 100044， China 现代城市测绘国家测绘地理信息局重点实验室，北京 100044 现代城市测绘国家测绘地理信息局重点实验室，北京 100044 Key Laboratory of Modern Urban Surveying and Mapping， National Administration of Surveying Mapping and Geoinformation， Beijing 100044， China Key Laboratory of Modern Urban Surveying and Mapping， National Administration of Surveying Mapping and Geoinformation， Beijing 100044， China 建筑遗产精细重构与健康监测北京市重点实验室 北京建筑大学，北京 100044 建筑遗产精细重构与健康监测北京市重点实验室 北京建筑大学，北京 100044 Beijing Key Laboratory for Architectural Heritage Fine Reconstruction & Health Monitoring， Beijing University of Civil Engineering and Architecture， Beijing 100044， China Beijing Key Laboratory for Architectural Heritage Fine Reconstruction & Health Monitoring， Beijing University of Civil Engineering and Architecture， Beijing 100044， China 北京申信达成科技有限公司，北京 102444 北京申信达成科技有限公司，北京 102444 Beijing Shenxindacheng Technology Co.， Beijing 102444， China Beijing Shenxindacheng Technology Co.， Beijing 102444， China 针对海量激光点云异地实时扫描处理困难和现场绝对坐标系的配准问题，结合北斗/GNSS定位系统和5G通讯技术研制一套低成本地基激光雷达测量系统。该测量系统由激光雷达、高精度电机、北斗/GNSS接收模块和5G模组集成，激光雷达采集点云信息，高精度电机获取角度信息，北斗/GNSS接收模块接收时间并通过自主研发的时间同步模块，与高精度电机与激光雷达进行时间标记，获取带有时间标签的点云和角度文件，并通过5G通讯技术实时传输。终端通过自主研发的数据预处理软件，通过线性插值算法进行多帧点云时空配准，室外根据北斗/GNSS坐标进行多站激光点云粗配准，室内无北斗/GNSS环境下通过单站特征点进行多站点云粗配准，再利用临近迭代算法进行精配准以完成整体配准，并在自主研发的实时点云管理与可视化系统上进行可视化。实验表明该测量系统使用5G通讯技术传输速率50 Mbit/s，可以实现异地实时扫描传输，配准后的点云误差在3 mm以下，可为数字孪生、物质文化遗产监测、特大异形建筑施工运维分析等实时应用领域提供数字基础设施。 In response to the difficulties of processing laser point cloud off-site real-time scanning and on-site absolute coordinate system alignment， a low-cost ground-based LiDAR measurement system is developed by combining the BeiDou/GNSS positioning system and 5G communication technology. The developed measurement system is integrated by LiDAR， a high-precision motor， BeiDou/GNSS receiver module， and a 5G module. LiDAR collects point cloud information， the high-precision motor acquires angle information， BeiDou/GNSS receiver module receives the time， the high-precision motor acquires time stamps using a self-developed time synchronization module， and LiDAR obtains the point cloud and angle file with time tag before transmitting in real-time via 5G communication technology. The terminal performs spatial and temporal alignment of multi-frame point clouds by linear interpolation algorithm through self-developed data pre-processing software， coarse alignment of the multi-station laser point clouds based on BeiDou/GNSS coordinates outdoors， coarse alignment of the multi-station point clouds by single-station feature points in an indoor environment without BeiDou/GNSS， fine alignment using the proximity iteration algorithm to complete overall alignment， and visualization on the self-developed real-time point cloud management and visualization system. The experimental results show that the measurement system achieves real-time off-site scanning and transmission at a transmission rate of 50 Mbit/s by using 5G communication technology. The error of the aligned point cloud is less than 3 mm， which can provide digital infrastructure for real-time applications such as digital twin， monitoring of material and cultural heritage， analysis of construction and operation， and maintenance of large buildings. 北斗 5G 地基激光雷达 点云配准 可视化 BeiDou 5G ground-based LIDAR point cloud registration 3D visualization 1 引 言 1 以激光点云为代表的三维地理信息数据已经成为国民经济建设的基石，行业应用领域对三维地理信息的需求量也日益剧增 ［ 1 1 - 2 2 1-2 ］ 。泛在测绘技术获取的点云数据是数字基础设施的核心，激光雷达测量技术是泛在测绘技术的重要组成部分，已经成为空间数据获取的重要手段 ［ 3 3 - 5 5 3-5 ］ ，但激光扫描获取海量点云数据处理困难，难以为数字孪生 ［ 6 6 ］ 、物质文化遗产监测 ［ 7 7 - 8 8 7-8 ］ 、特大异形建筑施工运维分析 ［ 9 9 ］ 等需要远程自动采集、存储和传输数据的实时应用领域提供数字基础设施。第五代移动通信技术（5G）时代的到来，以毫米波通信为代表的技术引起广泛关注，5G具有高速率、低时延、大容量等特征，采用5G通信技术更能保证数据的可靠性和及时性 ［ 10 10 - 11 11 10-11 ］ 。同时，我国自主研发的北斗导航系统已全面开通服务，可为全球用户提供高精度定位、导航和授时等服务 ［ 12 12 - 15 15 12-15 ］ 。为实现海量激光点云异地实时传输，结合北斗/GNSS定位和5G技术研制一款地基激光雷达测量系统。 对于北斗卫星导航系统和5G通讯技术的研究，目前大部分专家学者都将其应用于城市交通、无人驾驶等方面，对于北斗卫星导航系统和5G通讯技术在地基激光雷达方面的应用还比较稀缺。市场上现有地基激光雷达，如FARO公司的FARO Focus S350激光雷达，集GPS、指南针、高度传感器及双轴补偿器扫为一体，扫描距离达到350 m，扫描速度高达每秒976 000点。RIEGL公司推出的VZ-400i地面激光雷达，内置惯导，提供了3G/4G通讯接口、Wi-Fi和以太网交互硬件，可通过WiFi和4G将激光点云数据实时导入云端，配备处理平台可以同时进行数据在线拼接。该公司推出的VZ系列可进行超长测距的VZ-6000三维激光扫描仪，提供超过6 000 m的超长距离测量能力，甚至可以在沙尘、雾天、雨天、雪天等能见度较低的情况下使用并进行多重目标回波的识别，适用于雪地和冰川测量。天宝公司研制的TrimbleTX8三维激光扫描仪测距范围约为120 m，凭借LightningTM技术，数据获取速度高达每秒一百万个精确激光点，可以获得具有良好完整性的数据结果，内置集成HDR相机可实现快速图像获取用于生成彩色扫描，大大提高外业生产率。国内如数字绿土公司推出的LiPod大空间激光雷达集激光扫描、影像采集、GNSS ［ 16 16 ］ 坐标获取为一体，无需标靶特征，短时间内即可快速采集数据，测距精度为±3 cm，并通过WIFI热点的方式进行设备连接。北科天绘的U-Arm系列地面三维激光激光雷达，内置GNSS板卡，数据采集频率最高可达600 KHz。以上市场上主流公司研制推出的激光雷达虽大都内置GNSS板卡，但没有一款结合5G通讯技术，更不具备海量激光点云实时处理系统，无法进行海量数据的实时传输处理 ［ 17 17 - 22 22 17-22 ］ 。为了填补这一市场空白，解决海量激光点云实时处理困难问题，本文融合北斗/GNSS定位和5G技术研制了一款地基激光雷达测量系统。该激光雷达测量采用独特的非重复扫描方式，能获得更高的点云密度和探测距离，可探测远至260 m外的物体，数据采集频率可达100 KHz，内置北斗/GNSS接收模块和5G模组，可通过5G实时云传输，并利用北斗/GNSS接收模块获取时间进行时间同步，获取坐标信息并结合线性插值算法、罗德里格矩阵及ICP算法进行点云实时配准获取完整点云，并采用分层调度策略进行可视化。 本文分析了海量激光点云实时处理需求，提出结合5G通讯技术和北斗/GNSS定位模块研制一款地基激光雷达测量系统。利用激光雷达模块获取激光点云，高精度电机获取角度信息，北斗/GNSS获取时间并进行时间同步，通过5G技术实现海量激光点云异地实时传输，终端通过线性插值算法进行多帧点云配准，根据北斗/GNSS提供坐标中心进行多站激光点云粗配准，室内无GNSS环境则根据室内特征进行多站点云数据粗配准，并采用ICP算法进行精配准以完成整体配准，通过LOD调度策略实现海量激光点云高效可视化。实验表明，该测量系统可通过5G通讯技术异地实时发送海量激光点云数据，并配备数据预处理程序，填补了利用北斗/GNSS和5G技术研制激光雷达这一领域的空白，具备相当市场价值。 2 地基激光雷达测量系统集成 1 2.1　系统总体设计 2 本文采用多线式三维激光扫描仪作为激光雷达传感器，使用国产高精度电机作为激光雷达的旋转平台，配置指北针和双水准气泡并结合北斗卫星导航系统及5G通讯模组研制地基激光雷达测量系统。指北针和双水准气泡保证扫描初始指北，激光雷达扫描获取被测场景点云信息，由于激光雷达采集到的所有的帧数据都是同一个角度，为得到完整的点云数据配备外部高精度电机，以获取角度信息。通过北斗/GNSS接收模块获取时间信息，嵌入式微控制器将时间信息解码成时间戳，并利用自主研发的时间同步模块将时间信息同步给激光雷达和高精度电机。通过5G模组，将实时扫描数据发送给远程终端，终端接收数据后通过自主研发的数据预处理程序进行处理得到完整点云数据。对处理后的激光点云数据进行数据分割及格式转换，得到具有多细节层次的3DTiles瓦片结构，并结合LOD调度策略实现激光点云数据高效可视化。地基激光雷达测量系统整体设计如 图1 图1 所示。 2.2　系统硬件集成 2 如 图2 图2 所示，测量系统整体以三脚架为支撑，并配置指北针和双水准气泡，保证每站扫描初始方向一致，矫正电机零位，电机零位和激光雷达初始位置一致，扫描阶段激光雷达旋转中心和电机旋转轴一致，保证了旋转过程中相对位置一致。系统主控芯片采用STM32F103CBT6，通过自主编写程序激光雷达启动，同时通过嵌入式微控制器操控高精度电机进行旋转工作，北斗/GNSS接收机接收GPRMC时间信号和PPS秒脉冲数据，通过RS232总线将时间信号传输给嵌入式微控制器，同时将PPS信号传输并转换成嵌入式微控制器能够识别和采集的TTL串口信号。嵌入式微控制器将时间信号解码转化成UNIX时间戳，最终通过并通过485接口传输给高精度电机、激光雷达。计算机获取带有时间标签文件，利用5G通讯技术将采集数据实时发送到到远程终端，终端通过自主研发的配准程序进行数据处理后并进行可视化展示。 3 激光雷达测量系统时空配准 1 3.1　多模块时间同步 2 为了获取完整点云，以时间标为索引，通过对标时间标签，使点云帧数据旋转对应的角度，进而实现点云的配准。激光雷达扫描得到的点云可记录海量的空间三维点坐标，但点云并不带有时间标签，需要接入外部时间，而且本测量系统采用的激光雷达需要输入GNRMC/GPRMC格式的时间数据，并通过pps秒脉冲信号的方式进行时间同步；高精度电机也同样不配备时间设备，同样需要外部时间信息，因此可以通过接入北斗模块的方式进行时间同步，北斗/GNSS接收机将产生的PPS信号输入到激光雷达的PPS信号输入接口，时间字符串通过RS232串口发送到激光雷达的时间数据接口，频率与PPS信号相同。激光雷达接收到一对PPS和时间字符串之后，内部高精度定时器与PPS对准，将对准的某一PPS对应的时间字符串作为激光雷达的时间基准，实现激光雷达时间同步，同时嵌入式微控制器将时间信号解码转化成时间戳通过485接口传输高精度电机，实现电机时间同步。 室内工作室通常是处于无GNSS的环境，如何在无GNSS下对点云进行配准，以获取全部点云也是需要解决的问题。为解决该问题，利用嵌入式开发，给嵌入式微控制器编写了一个可以产生虚拟GNRMC/GPRMC时间和PPS秒脉冲信号的功能。嵌入式微控制器产生PPS秒脉冲通过RS232转485的方式将PPS信号输入激光雷达，产生的虚拟GNRMC/GPRMC时间通过SDK协议传输给激光雷达，使激光雷达的点云数据带有时间标签；嵌入式微控制器还接收高精度电机返回的角度数据同时打标虚拟GNRMC/GPRMC时间并传输给计算机，使得计算机得到带有时间标签的角度数据。 3.2　点云数据整体配准 2 3 . 2 . 1　多帧点云数据配准 3 激光雷达扫描获取得到单帧点云文件，由于激光雷达不配备旋转电机，因此采集到的所有的帧数据都会是同一个角度，若想得到完整的点云数据则需要对采集到的点云数据匹配角度信息。本文采用时空同步方法进行多帧点云配准，安装阶段电机零位和激光雷达初始位置保持一致，实验前矫正电机零位，使电机零位固定，扫描阶段，激光雷达旋转中心和电机旋转轴一致，保证了旋转过程中相对位置一致。并采用时间同步思路，根据时间信息将外部电机获取角度信息匹配给点云数据。主要方法是遍历文件寻找两个文件对应时间，通过对标时间标签使相同时间的每帧点云数据匹配相应角度数据，并采用线性插值算法使得无法对标时间的点云数据获取对应角度得到完整点云文件。 多帧点云配准对标时间获取空间变换参数进而完成配准，初始点 X （ x ， y ， z ，1）与配准点 X '（ x ' ， y '， z '，1）关系如 式（1） 式（1） 所示： <math display="block" id="M1"><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mtable columnalign="left"><mtr><mtd><msup><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msup></mtd></mtr><mtr><mtd><msup><mrow><mi mathvariant="bold-italic">y</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msup></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mtable columnalign="left"><mtr><mtd><msup><mrow><mi mathvariant="bold-italic">z</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msup></mtd></mtr><mtr><mtd><mn mathvariant="normal">1</mn></mtd></mtr></mtable></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi mathvariant="bold">A</mi></mtd><mtd><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi mathvariant="bold">A</mi></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi mathvariant="bold">A</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi mathvariant="bold">A</mi></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">1</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">1</mn></mtd></mtr></mtable></mrow></mfenced><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mtable columnalign="left"><mtr><mtd><mi mathvariant="bold-italic">x</mi></mtd></mtr><mtr><mtd><mi mathvariant="bold-italic">y</mi></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mtable columnalign="left"><mtr><mtd><mi mathvariant="bold-italic">z</mi></mtd></mtr><mtr><mtd><mn mathvariant="normal">1</mn></mtd></mtr></mtable></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436331&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436330&type= 54.01733780 18.71133423 ， <math display="block" id="M2"><mi mathvariant="bold-italic">A</mi><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">A</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mfrac><mrow><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">A</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">A</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfrac><mfenced separators="|"><mrow><mi mathvariant="bold-italic">T</mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436335&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436334&type= 48.34466553 8.97466660 ， 式中： A 为角度， T 为时间。最终配准效果如 图3 图3 所示。 3 . 2 . 2　多站点云数据整体配准 3 由于被测场景复杂、仪器设备影响等原因，通常需要多个视角扫描，才能获取完整扫描场景点云，但多次扫描的点云可能不在同一坐标系下，存在着空间旋转错位关系。因此需对不同位置、不同视角扫描得到的点云通过坐标变换进行配准，将这些点云合并到同一个坐标系下，形成完整点云数据。室外处于GNSS环境下可根据北斗/GNSS测站信息作为坐标中心点进行粗配准，再使用ICP算法精配准，室内往往处于无北斗/GNSS环境下，但室内特征明显，可通过特征进行粗配准再使用ICP算法精配准。室外点云配准流程如 图4 图4 所示。 点云数据粗配准是将两站点云大致配准到同一坐标系下，为精配准提供良好的初始值。地基激光雷达测量系统配备指北针和双水准气泡，设定每站扫描启动阶段回到指北针初始位置，双水准气泡在近距离范围内水平与垂直定向，保证每站都以指北针方向为初始方向，结合北斗/GNSS模块提供的北斗/GNSS测站信息作为坐标中心点，在模块安装阶段将其换算到激光雷达模块中心，从而保证多站点云坐标能够转换到同一坐标系下，完成多站激光点云粗配准。 室内特征配准根据SIFT特征点提取算法获取单站点云内的关键点，以关键点坐标欧式距离表征的三维位置相似度为基础，来估算它们的对应关系，确定同名特征点，一般由3个或3个以上不同的同名特征点计算出空间转换参数进而完成配准，相邻两测站中的同名特征点分别为 <math id="M3"><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mfenced separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">y</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">z</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436357&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436355&type= 25.39999962 4.82600021 与 <math id="M4"><mi mathvariant="bold-italic">X</mi><mo>=</mo><mfenced separators="|"><mrow><mi mathvariant="bold-italic">x</mi><mo>,</mo><mi mathvariant="bold-italic">y</mi><mo>,</mo><mi mathvariant="bold-italic">z</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436363&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436360&type= 20.48933220 4.82600021 ，两点间偏移量为 <math id="M5"><mi mathvariant="normal">Δ</mi><mi mathvariant="bold-italic">X</mi><mo stretchy="false">(</mo><mi mathvariant="normal">Δ</mi><mi mathvariant="bold-italic">x</mi><mo>,</mo><mi mathvariant="normal">Δ</mi><mi mathvariant="bold-italic">y</mi><mo>,</mo><mi mathvariant="normal">Δ</mi><mi mathvariant="bold-italic">z</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436368&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436365&type= 25.90800095 3.97933316 ，两点之间有如 式（3） 式（3） 所示的关系： <math display="block" id="M6"><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mfenced separators="|"><mrow><mi mathvariant="bold-italic">λ</mi><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">X</mi><mo>+</mo><mi mathvariant="normal">Δ</mi><mi mathvariant="bold-italic">X</mi></mrow></mfenced><mo>=</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436373&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436371&type= 33.52799988 4.06400013 0， 其中： λ 是尺度参数，点云变换中尺度不变即 λ =1； <math id="M7"><mi mathvariant="normal">Δ</mi><mi mathvariant="bold-italic">X</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436379&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436376&type= 5.41866684 2.79399991 为偏移量，重心坐标分别为 X m 和 X n ，重心化坐标为： <math display="block" id="M8"><msub><mrow><mover accent="true"><mi mathvariant="bold-italic">X</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>,</mo><mover accent="true"><mi mathvariant="bold-italic">X</mi><mo>¯</mo></mover><mo>=</mo><mi mathvariant="bold-italic">X</mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436384&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436380&type= 42.75666428 4.65666676 . 坐标重心化后，同名点间的坐标转化如 式（5） 式（5） ： <math display="block" id="M9"><msub><mrow><mover accent="true"><mi mathvariant="bold-italic">X</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">R</mi><mover accent="true"><mi mathvariant="bold-italic">X</mi><mo>¯</mo></mover><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436388&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436385&type= 20.91266632 4.23333359 . 线性化后得到旋转参数的误差方程： <math display="block" id="M10"><mi mathvariant="bold-italic">V</mi><mo>=</mo><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">t</mi><mo>-</mo><mi mathvariant="bold-italic">L</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436392&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436390&type= 18.88066673 2.96333337 ， 其中： <math id="M11"><mi mathvariant="bold-italic">A</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mo>-</mo><msub><mrow><mover><mrow><mi mathvariant="bold-italic">z</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">z</mi></mtd><mtd><mo>-</mo><msub><mrow><mover><mrow><mi mathvariant="bold-italic">y</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">y</mi></mtd></mtr><mtr><mtd><mo>-</mo><msub><mrow><mover><mrow><mi mathvariant="bold-italic">z</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">z</mi></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><msub><mrow><mover><mrow><mi mathvariant="bold-italic">x</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mi mathvariant="bold-italic">x</mi></mtd></mtr><mtr><mtd><msub><mrow><mover><mrow><mi mathvariant="bold-italic">y</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mi mathvariant="bold-italic">y</mi></mtd><mtd><msub><mrow><mover><mrow><mi mathvariant="bold-italic">x</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mi mathvariant="bold-italic">x</mi></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436397&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436395&type= 58.75866318 18.79599953 ， <math id="M12"><mi mathvariant="bold-italic">t</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mi mathvariant="bold-italic">a</mi></mtd></mtr><mtr><mtd><mi mathvariant="bold-italic">b</mi></mtd></mtr><mtr><mtd><mi mathvariant="bold-italic">c</mi></mtd></mtr></mtable></mrow></mfenced><mo>,</mo><mi mathvariant="bold-italic">L</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mover><mrow><mi mathvariant="bold-italic">x</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mover><mrow><mi mathvariant="bold-italic">x</mi></mrow><mo>¯</mo></mover></mtd></mtr><mtr><mtd><msub><mrow><mover><mrow><mi mathvariant="bold-italic">y</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">y</mi></mtd></mtr><mtr><mtd><msub><mrow><mover><mrow><mi mathvariant="bold-italic">z</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">z</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436402&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436400&type= 34.03599930 18.03400040 ，求解罗德里格参数构造矩阵 R 后，平移参数的误差方程为： <math display="block" id="M13"><mi mathvariant="bold-italic">V</mi><mo>=</mo><mi mathvariant="normal">Δ</mi><mi mathvariant="bold-italic">X</mi><mo>-</mo><mfenced separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">X</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436407&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436404&type= 35.47533035 3.97933316 . 粗配准后点云大致重合但精度不高，利用ICP算法再次配准提高配准精度。给定相邻两站点云 <math id="M14"><mi mathvariant="bold-italic">P</mi><mo>=</mo><msubsup><mrow><mfenced open="{" close="}" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub></mrow></mfenced></mrow><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow/><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>P</mi></mrow></msub></mrow></msub></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436411&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436409&type= 18.45733452 5.16466665 和 <math id="M15"><mi mathvariant="bold-italic">X</mi><mo>=</mo><msubsup><mrow><mfenced open="{" close="}" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub></mrow></mfenced></mrow><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow/><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>X</mi></mrow></msub></mrow></msub></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436416&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436414&type= 19.13466644 5.16466665 ，在 X 点集内搜索距离 P 中的每一个点 <math id="M16"><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436420&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436418&type= 3.21733332 3.72533321 的最近点 <math id="M17"><msub><mrow><mi mathvariant="bold-italic">Y</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436423&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436421&type= 3.55599999 3.72533321 ，建立点对的映射关系， <math id="M18"><mi mathvariant="bold-italic">P</mi><mo>=</mo><msubsup><mrow><mfenced open="{" close="}" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub></mrow></mfenced></mrow><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow/><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>P</mi></mrow></msub></mrow></msub></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436429&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436426&type= 18.45733452 5.16466665 的对应点集即为 <math id="M19"><mi mathvariant="bold-italic">Y</mi><mo>=</mo><msubsup><mrow><mfenced open="{" close="}" separators="|"><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub></mrow></mfenced></mrow><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow/><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>Y</mi></mrow></msub></mrow></msub></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436433&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436431&type= 18.79599953 5.16466665 ，即 Y = C （ P ， X ）其中 C 为对应点的坐标变换，不停迭代直到满足精度为止。粗配准后的点云误差在7 mm以下，经过ICP算法再次配准的误差在3 mm以下。本节以校园内教学楼为实验区域，使用地基激光雷达测量系统，对教学楼一面进行布站扫描，并使用自主研发数据预处理程序进行处理得到完整点云，实现了自动化扫描与配准，有效节省了数据处理时间，配准后的数据未出现分层情况，精度满足各实时应用领域的需求。室外点云整体配准结果如 图5 图5 所示。 4 实时点云管理与可视化系统 1 4.1　激光点云数据5G实时传输 2 由于Wi-Fi网络的抗干扰性差，4G网络的传输速率和时延的限制，造成海量激光点云实时传输困难，而5G网络的高带宽、低延时、大容量等特性正可以解决这些技术瓶颈。地基激光雷达测量系统工作获取海量激光点云点云文件和角度文件，将采集的全部数据打包压缩，传输给5G模块，当5G通讯模块接受到数据后，利用云端服务器，通过内网穿透映射传送给远程终端。通过测试证明使用5G模组传输速率可达50 Mbit/s，能够实现海量激光点云实时传输，具体传输流程如 图6 图6 所示。 4.2　终端激光点云数据加载与可视化 2 本文采用自主开发的实时点云管理与可视化系统，主要由点云数据传输模块、点云数据预处理模块、GIS分析模块、功能管理模块、三维场景模块、三维比对模块组成，各个模块间各自分工又相互协作，可对地基激光测量系统获取的激光点云数据进一步管理和应用。系统架构如 图7 图7 所示。 系统主界面如 图8 图8 所示，用户可进行三维数据分析、数据加载、实时监测以及三维场景可视化等。 为了实现海量激光点云高效可视化对点云构建具有多细节层次的3DTiles瓦片结构，通过八叉树分割算法将场景文件分割成多细节层次散列文件，将散列文件转化为pnts瓦片数据格式，并将瓦片数据的地理空间位置、旋转矩阵、几何误差等信息按照八叉树空间数据结构生成索引文件，最终生成3DTiles瓦片结构。激光点云构建3D Tiles的流程如 图9 图9 所示。 同时在可视化过程中引入LOD调度策略，使用屏幕空间误差 <math id="M20"><mi>S</mi><mi>S</mi><mi>E</mi><mo>=</mo><mi>g</mi><mi>e</mi><mi>o</mi><mi>m</mi><mi>e</mi><mi>t</mi><mi>r</mi><mi>i</mi><mi>c</mi><mi>E</mi><mi>r</mi><mi>r</mi><mi>o</mi><mi>r</mi><mo>×</mo><mi>h</mi><mi>e</mi><mi>i</mi><mi>g</mi><mi>h</mi><mi>t</mi><mo>/</mo><mi>d</mi><mi>i</mi><mi>s</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi>c</mi><mi>e</mi><mo>×</mo><mi>s</mi><mi>s</mi><mi>e</mi><mi>D</mi><mi>e</mi><mi>n</mi><mi>o</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>a</mi><mi>t</mi><mi>o</mi><mi>r</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436480&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40436478&type= 94.57266998 3.89466691 来决定瓦片的细节层次，去加载调度不同级别的瓦片数据。数据可视化效果如 图10 图10 所示。 5 结 论 1 本文将5G通讯技术与北斗/GNSS定位模块相结合研制了一款价格低廉、携带方便的地基激光雷达测量系统，并配备自主研制的时间同步模块、点云预处理程序及实时点云管理与可视化系统，实现海量激光点云从采集、处理到应用一体化。实验结果表明，5G通信技术可以有效解决海量激光点云远程实时传输问题，配备自主研制点云预处理程序可供用户快速获取完整点云数据，可为数字孪生、物质文化遗产监测、特大异形建筑施工运维分析等实时应用领域提供数字基础设施，具有广泛的应用价值和较大的经济效益。 参考文献 刘大同 ， 郭凯 ， 王本宽 ， 等 . 数字孪生技术综述与展望 ［J］. 仪器仪表学报 ， 2018 ， 39 （ 11 ）： 1 - 10 . 数字孪生技术综述与展望 LIU D T ， GUO K ， WANG B K ， et al . Summary and perspective survey on digital twin technology ［J］. Chinese Journal of Scientific Instrument ， 2018 ， 39 （ 11 ）： 1 - 10 . （in Chinese） Summary and perspective survey on digital twin technology 伍锡如 ， 薛其威 . 基于激光雷达的无人驾驶系统三维车辆检测 ［J］. 光学 精密工程 ， 2022 ， 30 （ 4 ）： 489 - 497 . doi: 10.37188/OPE.20223004.0489 http://dx.doi.org/10.37188/OPE.20223004.0489 基于激光雷达的无人驾驶系统三维车辆检测 WU X R ， XUE Q W . 3D vehicle detection for unmanned driving systerm based on lidar ［J］. Opt. Precision Eng. ， 2022 ， 30 （ 4 ）： 489 - 497 . （in Chinese） . doi: 10.37188/OPE.20223004.0489 http://dx.doi.org/10.37188/OPE.20223004.0489 3D vehicle detection for unmanned driving systerm based on lidar GUO M ， SUN M X ， ZHOU T F ， et al . Novel trajectory optimization algorithm of vehicle-borne LiDAR mobile measurement system ［J］. Sensors and Materials ， 2020 ， 32 （ 11 ）： 3935 . doi: 10.18494/sam.2020.3052 http://dx.doi.org/10.18494/sam.2020.3052 Novel trajectory optimization algorithm of vehicle-borne LiDAR mobile measurement system 乔俊峰 ， 周沅桢 ， 王永 ， 等 . 三维激光扫描测体积技术及其应用进展 ［J］. 激光与红外 ， 2021 ， 51 （ 9 ）： 1115 - 1122 . doi: 10.3969/j.issn.1001-5078.2021.09.001 http://dx.doi.org/10.3969/j.issn.1001-5078.2021.09.001 三维激光扫描测体积技术及其应用进展 QIAO J F ， ZHOU Y ZH ， WANG Y ， et al . Advances in 3D laser scanning volume measurement technology and application ［J］. Laser & Infrared ， 2021 ， 51 （ 9 ）： 1115 - 1122 . （in Chinese） . doi: 10.3969/j.issn.1001-5078.2021.09.001 http://dx.doi.org/10.3969/j.issn.1001-5078.2021.09.001 Advances in 3D laser scanning volume measurement technology and application GUO M ， ZHOU Y Q ， ZHAO J H ， et al . Urban Geospatial Information Acquisition Mobile Mapping System based on close-range photogrammetry and IGS site calibration ［J］. Geo-Spatial Information Science ， 2021 ， 24 （ 4 ）： 558 - 579 . doi: 10.1080/10095020.2021.1924084 http://dx.doi.org/10.1080/10095020.2021.1924084 Urban Geospatial Information Acquisition Mobile Mapping System based on close-range photogrammetry and IGS site calibration 张旭辉 ， 吕欣媛 ， 王甜 ， 等 . 数字孪生驱动的掘进机器人决策控制系统研究 ［J］. 煤炭科学技术 ， 2022 ， 50 （ 7 ）： 36 - 49 . 数字孪生驱动的掘进机器人决策控制系统研究 ZHANG X H ， LYU X Y ， WANG T ， et al . Research on decision control system of tunneling robot driven by digital twin ［J］. Coal Science and Technology ， 2022 ， 50 （ 7 ）： 36 - 49 . （in Chinese） Research on decision control system of tunneling robot driven by digital twin MING ， GUO . Normal cloud model theory-based comprehensive fuzzy assessment of wooden pagoda safety ［J］. Journal of Cultural Heritage ， 2022 ， 55 ： 1 - 10 . doi: 10.1016/j.culher.2022.02.001 http://dx.doi.org/10.1016/j.culher.2022.02.001 Normal cloud model theory-based comprehensive fuzzy assessment of wooden pagoda safety 郭明 ， 闫冰男 ， 周腾飞 ， 等 . 激光雷达技术在应县木塔形变分析中的应用 ［J］. 建筑科学与工程学报 ， 2020 ， 37 （ 2 ）： 109 - 117 . 激光雷达技术在应县木塔形变分析中的应用 GUO M ， YAN B N ， ZHOU T F ， et al . Application of LiDAR technology in deformation analysis of Yingxian wooden pagoda ［J］. Journal of Architecture and Civil Engineering ， 2020 ， 37 （ 2 ）： 109 - 117 . （in Chinese） Application of LiDAR technology in deformation analysis of Yingxian wooden pagoda 郭明 ， 孙梦溪 ， 黄明 ， 等 . 融合激光雷达与无人机的特大钢结构高精度测量 ［J］. 光学 精密工程 ， 2021 ， 29 （ 5 ）： 989 - 998 . doi: 10.37188/OPE.20212905.0989 http://dx.doi.org/10.37188/OPE.20212905.0989 融合激光雷达与无人机的特大钢结构高精度测量 GUO M ， SUN M X ， HUANG M ， et al . High-precision measurement of steel structure based on LiDAR and UAV ［J］. Opt. Precision Eng. ， 2021 ， 29 （ 5 ）： 989 - 998 . （in Chinese） . doi: 10.37188/OPE.20212905.0989 http://dx.doi.org/10.37188/OPE.20212905.0989 High-precision measurement of steel structure based on LiDAR and UAV NAVARRO-ORTIZ J ， ROMERO-DIAZ P ， SENDRA S ， et al . A survey on 5G usage scenarios and traffic models ［J］. IEEE Communications Surveys & Tutorials ， 2020 ， 22 （ 2 ）： 905 - 929 . doi: 10.1109/comst.2020.2971781 http://dx.doi.org/10.1109/comst.2020.2971781 A survey on 5G usage scenarios and traffic models DANGI R ， LALWANI P ， CHOUDHARY G ， et al . Study and investigation on 5G technology： a systematic review ［J］. Sensors （Basel， Switzerland） ， 2021 ， 22 （ 1 ）： 26 . doi: 10.3390/s22010026 http://dx.doi.org/10.3390/s22010026 Study and investigation on 5G technology： a systematic review 谢军 ， 庄建楼 ， 康成斌 . 基于北斗系统的物联网技术与应用 ［J］. 南京航空航天大学学报 ， 2021 ， 53 （ 3 ）： 329 - 337 . 基于北斗系统的物联网技术与应用 XIE J ， ZHUANG J L ， KANG CH B . Internet of Things technology and application based on Beidou system ［J］. Journal of Nanjing University of Aeronautics & Astronautics ， 2021 ， 53 （ 3 ）： 329 - 337 . （in Chinese） Internet of Things technology and application based on Beidou system 高志钰 ， 郭进义 ， 刘杰 . 北斗在地壳形变监测中的应用进展 ［J］. 测绘通报 ， 2022 （ 3 ）： 32 - 35 . doi: 10.3969/j.issn.0494-0911.2022.3.chtb202203007 http://dx.doi.org/10.3969/j.issn.0494-0911.2022.3.chtb202203007 北斗在地壳形变监测中的应用进展 GAO ZH Y ， GUO J Y ， LIU J . Application progress of BDS in monitoring crustal deformation ［J］. Bulletin of Surveying and Mapping ， 2022 （ 3 ）： 32 - 35 . （in Chinese） . doi: 10.3969/j.issn.0494-0911.2022.3.chtb202203007 http://dx.doi.org/10.3969/j.issn.0494-0911.2022.3.chtb202203007 Application progress of BDS in monitoring crustal deformation XIE J ， KANG C B . Engineering innovation and the development of the BDS-3 navigation constellation ［J］. Engineering ， 2021 ， 7 （ 5 ）： 37 - 49 . doi: 10.1016/j.eng.2021.04.002 http://dx.doi.org/10.1016/j.eng.2021.04.002 Engineering innovation and the development of the BDS-3 navigation constellation WEISSPETER . Welcome to the global navigation multi-constellation ［J］. Engineering ， 2021 ， 7 （ 4 ）： 421 - 423 . doi: 10.1016/j.eng.2021.03.006 http://dx.doi.org/10.1016/j.eng.2021.03.006 Welcome to the global navigation multi-constellation 郭明 ， 周玉泉 ， 陈才 ， 等 . 北斗导航授时的移动激光雷达测量系统时间同步装置设计 ［J］. 红外与激光工程 ， 2020 ， 49 （ S2 ）： 33 - 42 . 北斗导航授时的移动激光雷达测量系统时间同步装置设计 GUO M ， ZHOU Y Q ， CHEN C ， et al . Design of time synchronization device for mobile LiDAR measurement system based on BeiDou navigation timing ［J］. Infrared and Laser Engineering ， 2020 ， 49 （ S2 ）： 33 - 42 . （in Chinese） Design of time synchronization device for mobile LiDAR measurement system based on BeiDou navigation timing LI L ， CAO X Y ， SUN J . Three-dimensional point cloud registration based on normal vector angle ［J］. Journal of the Indian Society of Remote Sensing ， 2019 ， 47 （ 4 ）： 585 - 593 . doi: 10.1007/s12524-018-0918-4 http://dx.doi.org/10.1007/s12524-018-0918-4 Three-dimensional point cloud registration based on normal vector angle ZHANG J Y ， YAO Y X ， DENG B L . Fast and robust iterative closest point ［J］. IEEE Transactions on Pattern Analysis and Machine Intelligence ， 2022 ， 44 （ 7 ）： 3450 - 3466 . Fast and robust iterative closest point 郭清达 ， 全燕鸣 ， 姜长城 ， 等 . 应用摄像机位姿估计的点云初始配准 ［J］. 光学 精密工程 ， 2017 ， 25 （ 6 ）： 1635 - 1644 . doi: 10.3788/ope.20172506.1635 http://dx.doi.org/10.3788/ope.20172506.1635 应用摄像机位姿估计的点云初始配准 GUO Q D ， QUAN Y M ， JIANG CH CH ， et al . Initial registration of point clouds using camera pose estimation ［J］. Opt. Precision Eng. ， 2017 ， 25 （ 6 ）： 1635 - 1644 . （in Chinese） . doi: 10.3788/ope.20172506.1635 http://dx.doi.org/10.3788/ope.20172506.1635 Initial registration of point clouds using camera pose estimation GUO M ， YAN B N ， WANG G L ， et al . Iterative automatic global registration algorithm for multi-view point cloud of underground tunnel space ［J］. Measurement and Control ， 2021 ， 54 （ 3/4 ）： 385 - 395 . doi: 10.1177/00202940211003935 http://dx.doi.org/10.1177/00202940211003935 Iterative automatic global registration algorithm for multi-view point cloud of underground tunnel space 曾如铁 ， 徐金鸿 . 基于PCL的点云数据粗配准算法研究 ［J］. 测绘科学技术 ， 2019 （ 1 ）： 18 - 23 . doi: 10.12677/gst.2019.71004 http://dx.doi.org/10.12677/gst.2019.71004 基于PCL的点云数据粗配准算法研究 ZENG R T ， XU J H . PCL-based coarse registration algorithm for point cloud data ［J］. Geomatics Science and Technology ， 2019 （ 1 ）： 18 - 23 . （in Chinese） . doi: 10.12677/gst.2019.71004 http://dx.doi.org/10.12677/gst.2019.71004 PCL-based coarse registration algorithm for point cloud data 王欣 ， 张明明 ， 于晓 ， 等 . 应用改进迭代最近点方法的点云数据配准 ［J］. 光学 精密工程 ， 2012 ， 20 （ 9 ）： 2068 - 2077 . doi: 10.3788/ope.20122009.2068 http://dx.doi.org/10.3788/ope.20122009.2068 应用改进迭代最近点方法的点云数据配准 WANG X ， ZHANG M M ， YU X ， et al . Point cloud registration based on improved iterative closest point method ［J］. Opt. Precision Eng. ， 2012 ， 20 （ 9 ）： 2068 - 2077 . （in Chinese） . doi: 10.3788/ope.20122009.2068 http://dx.doi.org/10.3788/ope.20122009.2068 Point cloud registration based on improved iterative closest point method 10.37188/OPE.20233104.0450 TN967.1 A 1004-924X（2023）04-0450-09 郭明,齐慧慧,郭可才等.融合北斗/GNSS定位和5G通讯的地基激光雷达测量系统[J].光学精密工程,2023,31(04):450-458. 郭明,齐慧慧,郭可才等.融合北斗/GNSS定位和5G通讯的地基激光雷达测量系统[J].光学精密工程,2023,31(04):450-458. GUO Ming,QI Huihui,GUO Kecai,et al.Ground-based LiDAR measurement system integrating BeiDou/GNSS positioning and 5G communication[J].Optics and Precision Engineering,2023,31(04):450-458. GUO Ming,QI Huihui,GUO Kecai,et al.Ground-based LiDAR measurement system integrating BeiDou/GNSS positioning and 5G communication[J].Optics and Precision Engineering,2023,31(04):450-458. 国家自然科学基金资助项目（No. 41971350）；北京建筑大学金字塔人才培养工程建大杰青培养计划资助项目（No. JDJQ20200307）；北京市未来城市设计高精尖创新中心资助项目（No. UDC2019031724） 2022-09-26 2022-10-26 2023-02-25 2023-02-28 光学精密工程 4 31