CN118031904A - Expressway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud - Google Patents

Expressway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud Download PDF

Info

Publication number
CN118031904A
CN118031904A CN202410439311.1A CN202410439311A CN118031904A CN 118031904 A CN118031904 A CN 118031904A CN 202410439311 A CN202410439311 A CN 202410439311A CN 118031904 A CN118031904 A CN 118031904A
Authority
CN
China
Prior art keywords
point cloud
cloud data
laser point
data
dimensional
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202410439311.1A
Other languages
Chinese (zh)
Other versions
CN118031904B (en
Inventor
南轲
叶朋飞
李升甫
贾洋
达乾龙
易守东
蒲慧龙
汪致恒
杨天宇
孙晓鹏
罗文韬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Department of Transportation Highway Planning Prospecting and Design Research Institute
Original Assignee
Sichuan Department of Transportation Highway Planning Prospecting and Design Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sichuan Department of Transportation Highway Planning Prospecting and Design Research Institute filed Critical Sichuan Department of Transportation Highway Planning Prospecting and Design Research Institute
Priority to CN202410439311.1A priority Critical patent/CN118031904B/en
Publication of CN118031904A publication Critical patent/CN118031904A/en
Application granted granted Critical
Publication of CN118031904B publication Critical patent/CN118031904B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C7/00Tracing profiles
    • G01C7/06Tracing profiles of cavities, e.g. tunnels
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/26Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/30Noise filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/36Applying a local operator, i.e. means to operate on image points situated in the vicinity of a given point; Non-linear local filtering operations, e.g. median filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/50Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/761Proximity, similarity or dissimilarity measures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/764Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/50Context or environment of the image
    • G06V20/56Context or environment of the image exterior to a vehicle by using sensors mounted on the vehicle
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/64Three-dimensional objects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Health & Medical Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Nonlinear Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention belongs to the technical field of tunnel clearance measurement, and particularly relates to a highway tunnel clearance measurement method and device based on a vehicle-mounted laser point cloud. According to the method, a vehicle-mounted three-dimensional laser scanning system is used for acquiring three-dimensional laser point cloud data in an expressway tunnel, point cloud registration is carried out on the three-dimensional laser point cloud data based on actual measurement coordinates of a layout target, filtering is carried out on the three-dimensional laser point cloud data after registration, surface layer point cloud data are extracted, a ground irregular triangular network and a lining irregular triangular network are constructed by utilizing the surface layer point cloud data, and altitude values of the lining irregular triangular network and the ground irregular triangular network are differed to obtain clearance information of the tunnel. Compared with a manual measurement method, the method has higher measurement efficiency, higher integrity of the measured tunnel clearance information and convenience for repeated measurement.

Description

基于车载激光点云的高速公路隧道净空测量方法及装置Highway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud

技术领域Technical Field

本发明属于隧道净空测量技术领域,特别涉及基于车载激光点云的高速公路隧道净空测量方法及装置。The invention belongs to the technical field of tunnel clearance measurement, and in particular relates to a highway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud.

背景技术Background technique

隧道净空是指隧道内轮廓线所包围的空间,包括公路隧道建筑限界、通风及其它功能所需的断面积。断面形状和大小应根据结构设计力求得到最经济值。净空所包括的其它断面中,有通风机或通风管道、照明灯具及其它设备、监控设备和运行管理设备、电缆沟或电缆桥架、防灾设备等断面,以及富裕量和施工允许误差等。Tunnel clearance refers to the space enclosed by the tunnel contour, including the cross-sectional area required for the construction limit, ventilation and other functions of the highway tunnel. The cross-sectional shape and size should strive to obtain the most economical value according to the structural design. Other cross-sectional areas included in the clearance include ventilators or ventilation ducts, lighting fixtures and other equipment, monitoring equipment and operation management equipment, cable trenches or cable trays, disaster prevention equipment, etc., as well as margins and construction tolerances.

高速公路隧道作为交通基础设施中特别重要的一部分,对我国国民经济提升以及社会快速发展具有非常积极的推动作用。在高速公路不断新建的同时,过去三十年已经建设完成的高速公路隧道越来越多,长期的超负荷运营和养护维护不到位等情况造成许多在意隧道在运营期间都出现了不同程度的病害,如衬砌结构强度不足、开裂,路面开裂、错台、隆起,路面仰拱及仰拱填充厚度不足等等。在对病害处置的过程中通常采用衬砌工字钢钢架加固、重铺路面等等处治措施,这类处治方法都会对隧道净空造成侵占,减小行车道可用净空,对车辆行驶造成安全威胁。与此同时,由于隧道的变形,在病害处置的过程中可能会对道路标线的位置进行调整,新设计标线位置的隧道净空也是关注重点。因此,快速、便捷、准确的提取运营高速公路隧道净空是隧道病害处置和评估隧道运营安全性的重要环节。As a particularly important part of transportation infrastructure, highway tunnels have played a very positive role in promoting the improvement of my country's national economy and rapid social development. While new highways are being built, more and more highway tunnels have been built in the past three decades. Long-term overload operation and inadequate maintenance have caused many tunnels to suffer from varying degrees of defects during operation, such as insufficient lining structure strength and cracking, pavement cracking, misalignment, bulges, insufficient pavement inverts and invert filling thickness, etc. In the process of dealing with the defects, lining I-beam steel frame reinforcement, resurfacing of the road surface and other treatment measures are usually adopted. These treatment methods will encroach on the tunnel clearance, reduce the available clearance of the lane, and pose a safety threat to vehicle driving. At the same time, due to the deformation of the tunnel, the position of the road marking may be adjusted during the process of disease treatment, and the tunnel clearance at the newly designed marking position is also a focus of attention. Therefore, the rapid, convenient and accurate extraction of the clearance of operating highway tunnels is an important part of tunnel disease treatment and evaluation of tunnel operation safety.

随着激光雷达技术的快速发展以及在公路工程中的普及应用,在生产过程中,高速公路隧道净空测量主要分为两种方式,一种是通过人工方式上路实地测量,采用测距仪、全站仪等测量设备,这种方法工作效率低,断道时间长,对高速路网通行能力影响大,只能对特定位置进行测量,隧道净空信息不完整;另一种是通过三维激光扫描技术,目前这种方多采用站式激光雷达扫描设备,作业时间长,且主要服务于施工阶段的隧道净空收敛监测,作业过程繁琐,无法适应运营阶段隧道标线反复调整的净空测量需求。With the rapid development of LiDAR technology and its popular application in highway engineering, in the production process, highway tunnel clearance measurement is mainly divided into two methods. One is through manual on-the-road field measurement, using rangefinders, total stations and other measuring equipment. This method has low work efficiency, long road interruption time, and great impact on the traffic capacity of the highway network. It can only measure specific locations, and the tunnel clearance information is incomplete. The other is through three-dimensional laser scanning technology. At present, this method mostly uses station-type LiDAR scanning equipment, which has a long operation time and mainly serves the tunnel clearance convergence monitoring in the construction phase. The operation process is cumbersome and cannot meet the clearance measurement needs of repeated adjustment of tunnel markings in the operation phase.

发明内容Summary of the invention

本发明的目的在于克服现有高速公路隧道净空测量方法存在的测量效率低、所测隧道净空信息完整度不高和难以重复测量的问题,提供基于车载激光点云的高速公路隧道净空测量方法及装置。The purpose of the present invention is to overcome the problems of low measurement efficiency, low integrity of measured tunnel clearance information and difficulty in repeated measurement in existing highway tunnel clearance measurement methods, and to provide a highway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud.

为了实现上述发明目的,本发明提供了以下技术方案:In order to achieve the above-mentioned object of the invention, the present invention provides the following technical solutions:

基于车载激光点云的高速公路隧道净空测量方法,该方法包括如下步骤:A highway tunnel clearance measurement method based on vehicle-mounted laser point cloud, the method comprises the following steps:

S1:在待测高速公路隧道内布设若干标靶,测量得到各标靶中心的实测坐标,并使用车载激光扫描系统获取多批次的所述待测高速公路隧道的三维激光点云数据,所述三维激光点云数据包含所述标靶的点云数据;S1: placing a number of targets in the highway tunnel to be measured, measuring the measured coordinates of the center of each target, and using a vehicle-mounted laser scanning system to obtain multiple batches of three-dimensional laser point cloud data of the highway tunnel to be measured, wherein the three-dimensional laser point cloud data includes the point cloud data of the targets;

S2:去除所述三维激光点云数据中的异常点和噪声点后,基于所述各标靶中心的实测坐标对所述三维激光点云数据进行点云配准;S2: after removing abnormal points and noise points in the three-dimensional laser point cloud data, performing point cloud registration on the three-dimensional laser point cloud data based on the measured coordinates of the centers of the targets;

S3:对点云配准后的三维激光点云数据进行分割和滤波,并计算生成地面不规则三角网和衬砌不规则三角网;S3: segment and filter the 3D laser point cloud data after point cloud registration, and calculate and generate the ground irregular triangulated network and lining irregular triangulated network;

S4:将待测高速公路隧道的标线线位的连续线位数据或设计线位的连续线位数据转化为点数据,所述点数据所在位置的衬砌不规则三角网的高程值与地面不规则三角网的高程值作差,得到所述点数据所在位置的净空值。S4: Convert the continuous line position data of the marking line position or the continuous line position data of the design line position of the highway tunnel to be tested into point data, and subtract the elevation value of the lining irregular triangulated network at the location of the point data from the elevation value of the ground irregular triangulated network to obtain the clearance value at the location of the point data.

优选地,所述点云配准包括对不同批次的三维激光点云数据进行点云粗配准,和对点云粗配准后的三维激光点云数据进行点云精配准。Preferably, the point cloud registration includes performing rough point cloud registration on different batches of three-dimensional laser point cloud data, and performing fine point cloud registration on the three-dimensional laser point cloud data after the rough point cloud registration.

优选地,所述点云粗配准包括:Preferably, the point cloud coarse registration comprises:

S21:在当前批次的三维激光点云数据中,人工标记若干个覆盖单个标靶的点云数据的矩形区域,构建所述矩形区域对应的三维统计直方图,并对所述三维统计直方图进行平均化处理,得到基准三维统计直方图;S21: manually marking a number of rectangular areas of point cloud data covering a single target in the current batch of three-dimensional laser point cloud data, constructing a three-dimensional statistical histogram corresponding to the rectangular area, and averaging the three-dimensional statistical histogram to obtain a reference three-dimensional statistical histogram;

S22:按预设步距将未标记区域的三维激光点云数据,标记为与步骤S21中矩形区域尺寸相同的若干待匹配的矩形区域,构建所述待匹配的矩形区域对应的三维统计直方图,并与基准三维统计直方图进行相似度匹配,获取所有覆盖单个标靶的点云数据的矩形区域,以得到所有标靶的中心点坐标;其中,标靶的中心点坐标为所述覆盖单个标靶的点云数据的矩形区域内的三维激光点云数据的几何质心;S22: marking the three-dimensional laser point cloud data of the unmarked area into a plurality of rectangular areas to be matched with the same size as the rectangular area in step S21 according to a preset step size, constructing a three-dimensional statistical histogram corresponding to the rectangular area to be matched, and performing similarity matching with the reference three-dimensional statistical histogram, obtaining all rectangular areas of point cloud data covering a single target, so as to obtain the center point coordinates of all targets; wherein the center point coordinates of the target are the geometric centroid of the three-dimensional laser point cloud data within the rectangular area of point cloud data covering a single target;

S23:将各个所述标靶的中心点坐标纠正至与该标靶中心的实测坐标相同,得到该批次点云粗配准后的三维激光点云数据。S23: Correcting the coordinates of the center point of each target to be the same as the measured coordinates of the center of the target, and obtaining the three-dimensional laser point cloud data after the rough registration of the batch of point clouds.

优选地,所述三维统计直方图的构建方法包括:将每个矩形区域内的三维激光点云数据投影至该矩形区域内的三维激光点云数据的法向量方向的平面上,形成点云强度图像,生成含有灰度和梯度信息的三维统计直方图。Preferably, the method for constructing the three-dimensional statistical histogram includes: projecting the three-dimensional laser point cloud data within each rectangular area onto a plane in the direction of the normal vector of the three-dimensional laser point cloud data within the rectangular area to form a point cloud intensity image, and generating a three-dimensional statistical histogram containing grayscale and gradient information.

优选地,所述点云精配准包括:计算所有批次点云粗配准后的三维激光点云数据的中误差,选取中误差最小时对应批次的三维激光点云数据作为基准点云数据,采用ICP算法将其余批次的三维激光点云数据配准至所述基准点云数据处,得到点云配准后的三维激光点云数据。Preferably, the point cloud fine registration includes: calculating the mean square error of the three-dimensional laser point cloud data after the rough registration of all batches of point clouds, selecting the three-dimensional laser point cloud data of the corresponding batch when the mean square error is the smallest as the reference point cloud data, and using the ICP algorithm to register the remaining batches of three-dimensional laser point cloud data to the reference point cloud data to obtain the three-dimensional laser point cloud data after the point cloud registration.

优选地,所述S3包括:将所述点云配准后的三维激光点云数据分割为下部地面点云数据和上部衬砌点云数据;分别提取所述下部地面点云数据的表层点云数据和上部衬砌点云数据的表层点云数据,根据所述表层点云数据分别计算生成地面不规则三角网和衬砌不规则三角网。Preferably, S3 includes: dividing the three-dimensional laser point cloud data after the point cloud registration into lower ground point cloud data and upper lining point cloud data; respectively extracting the surface point cloud data of the lower ground point cloud data and the surface point cloud data of the upper lining point cloud data, and respectively calculating and generating a ground irregular triangulated network and a lining irregular triangulated network according to the surface point cloud data.

优选地,所述下部地面点云数据的表层点云数据使用渐进三角网加密滤波算法进行提取。Preferably, the surface point cloud data of the lower ground point cloud data is extracted using a progressive triangulation encryption filtering algorithm.

优选地,所述上部衬砌点云数据的表层点云数据提取方法,包括:从上部衬砌点云数据中,提取低于预设高程阈值的点云数据作为低点点云簇,再使用渐进三角网加密滤波算法提取上部衬砌点云数据的底部表层点云数据,将所述低点点云簇与所述底部表层点云数据合并,得到上部衬砌点云数据的表层点云数据。Preferably, the method for extracting surface point cloud data of the upper lining point cloud data comprises: extracting point cloud data below a preset elevation threshold from the upper lining point cloud data as a low-point point cloud cluster, and then using a progressive triangulation encryption filtering algorithm to extract the bottom surface point cloud data of the upper lining point cloud data, and merging the low-point point cloud cluster with the bottom surface point cloud data to obtain the surface point cloud data of the upper lining point cloud data.

优选地,所述地面不规则三角网和衬砌不规则三角网采用Delaunay三角网生长算法生成。Preferably, the ground irregular triangulated network and the lining irregular triangulated network are generated by using a Delaunay triangulated network growth algorithm.

基于车载激光点云的高速公路隧道净空测量装置,包括至少一个处理器,以及与所述至少一个处理器通信连接的存储器;所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行上述基于车载激光点云的高速公路隧道净空测量方法。A highway tunnel clearance measurement device based on vehicle-mounted laser point cloud includes at least one processor and a memory communicatively connected to the at least one processor; the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the above-mentioned highway tunnel clearance measurement method based on vehicle-mounted laser point cloud.

与现有技术相比,本发明的有益效果:Compared with the prior art, the present invention has the following beneficial effects:

本发明使用三维激光扫描系统获取高速公路隧道内部的三维激光点云数据,并基于布设标靶的实测坐标对三维激光点云数据进行点云配准,并对配准后的三维激光点云数据进行滤波,提取表层点云数据,并利用表层点云数据构建地面不规则三角网和衬砌不规则三角网,对衬砌不规则三角网和地面不规则三角网的高程值作差即可得到隧道的净空信息。相较于人工测量的方法,本发明方法通过一次测量工作便可获取全隧道区域各个位置的净高信息,所测的隧道净空信息完整度更高;在线位调整变更的情况下,只需重新离散线位并根据不规则三角网的高程值重新赋值作差便可得到变更后的位置的净空值,测量效率更高且便于重复测量。The present invention uses a three-dimensional laser scanning system to obtain three-dimensional laser point cloud data inside a highway tunnel, and performs point cloud registration on the three-dimensional laser point cloud data based on the measured coordinates of the deployed targets, and filters the registered three-dimensional laser point cloud data, extracts surface point cloud data, and uses the surface point cloud data to construct a ground irregular triangulated network and a lining irregular triangulated network, and the clearance information of the tunnel can be obtained by subtracting the elevation values of the lining irregular triangulated network and the ground irregular triangulated network. Compared with the manual measurement method, the method of the present invention can obtain the clear height information of each position in the entire tunnel area through one measurement, and the measured tunnel clearance information is more complete; in the case of line position adjustment and change, it is only necessary to re-discretize the line position and re-assign the difference according to the elevation value of the irregular triangulated network to obtain the clearance value of the changed position, which has higher measurement efficiency and is convenient for repeated measurement.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为本发明实施例1的基于车载激光点云的高速公路隧道净空测量方法的流程图;FIG1 is a flow chart of a highway tunnel clearance measurement method based on vehicle-mounted laser point cloud according to Embodiment 1 of the present invention;

图2为本发明实施例1的构建的三维统计直方图的示意图;FIG2 is a schematic diagram of a three-dimensional statistical histogram constructed in Example 1 of the present invention;

图3为本发明实施例1的完全封闭的高速公路单向隧道在横向截面方向上点云粗配准后的三维激光点云数据;FIG3 is a three-dimensional laser point cloud data of a completely enclosed one-way highway tunnel in Example 1 of the present invention after rough registration of point clouds in the transverse cross-sectional direction;

图4为本发明实施例1的完全封闭的高速公路单向隧道在横向截面方向上点云精配准后的三维激光点云数据;FIG4 is a three-dimensional laser point cloud data of a completely enclosed one-way highway tunnel in Example 1 of the present invention after precise registration of point clouds in the transverse cross-sectional direction;

图5为本发明实施例1的对完全封闭的高速公路单向隧道构建的地面不规则三角网;FIG5 is a ground irregular triangulated network constructed for a completely enclosed one-way tunnel of a highway according to Example 1 of the present invention;

图6为本发明实施例1的完全封闭的高速公路单向隧道构建的衬砌不规则三角网。FIG. 6 is an irregular triangulated network of the lining constructed for a completely enclosed one-way highway tunnel according to Example 1 of the present invention.

具体实施方式Detailed ways

下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。The present invention is further described in detail below in conjunction with test examples and specific implementation methods. However, this should not be understood as the scope of the above subject matter of the present invention being limited to the following embodiments, and all technologies realized based on the content of the present invention belong to the scope of the present invention.

实施例1Example 1

基于车载激光点云的高速公路隧道净空测量方法,其流程图如图1所示,包括以下步骤:The highway tunnel clearance measurement method based on vehicle-mounted laser point cloud is shown in the flowchart in Figure 1 and includes the following steps:

S1:在待测高速公路隧道内布设若干标靶,测量得到各标靶中心的实测坐标,并使用车载激光扫描系统获取多批次的待测高速公路隧道的三维激光点云数据,所述三维激光点云数据包含布设标靶的点云数据。S1: several targets are arranged in the highway tunnel to be tested, the measured coordinates of the center of each target are measured, and a vehicle-mounted laser scanning system is used to obtain multiple batches of three-dimensional laser point cloud data of the highway tunnel to be tested, wherein the three-dimensional laser point cloud data includes point cloud data of the arranged targets.

在高速公路隧道内,车载三维激光扫描系统难以获取GNSS(Global NavigationSatellite System,全球导航卫星系统)数据,因此对所采集的三维激光点云数据常常依赖标靶来进行点云配准;通常情况下标靶布设的密度与隧道的长度和线型复杂度成正比。本实施例根据待测高速公路隧道的长度、宽度和线型,设计适当的采集车速和标靶布设方案,根据实际测量经验,通常每隔10~20米布设一个标靶;而为了获取更密集的点云数据,采集车速通常与隧道宽度成反比。布设好标靶以后,使用全站仪测量所布设标靶中心的空间三维坐标作为实测坐标,并通过车载三维激光扫描系统扫描待测隧道以采集包含布设标靶的点云数据的三维激光点云数据,本实施例所采用的车载三维激光扫描系统型号为AS900-HL。In highway tunnels, it is difficult for the vehicle-mounted 3D laser scanning system to obtain GNSS (Global Navigation Satellite System) data, so the collected 3D laser point cloud data often relies on targets for point cloud registration; usually, the density of target layout is proportional to the length and linear complexity of the tunnel. This embodiment designs an appropriate collection speed and target layout plan based on the length, width and linearity of the highway tunnel to be measured. According to actual measurement experience, a target is usually laid out every 10 to 20 meters; and in order to obtain more dense point cloud data, the collection speed is usually inversely proportional to the width of the tunnel. After the target is laid out, the total station is used to measure the spatial three-dimensional coordinates of the center of the laid target as the measured coordinates, and the vehicle-mounted 3D laser scanning system is used to scan the tunnel to be measured to collect 3D laser point cloud data containing the point cloud data of the laid target. The model of the vehicle-mounted 3D laser scanning system used in this embodiment is AS900-HL.

对于高速公路隧道的车载三维激光点云数据采集工作,需要根据隧道环境和运营状态的不同制定差异化的采集策略:待测隧道若为完全封闭的高速公路单向隧道,则在快慢车道往返采集三维激光点云数据;若为完全封闭的高速公路双向隧道,则在双向的快车道往返采集三维激光点云数据;若为半封闭的高速公路单向隧道,则在封闭车道往返采集三维激光点云数据;若为半封闭的高速公路双向隧道,则在封闭车道各向往返采集三维激光点云数据。For the vehicle-mounted 3D laser point cloud data collection work of highway tunnels, it is necessary to formulate differentiated collection strategies according to the different tunnel environments and operating conditions: if the tunnel to be tested is a completely closed one-way tunnel on a highway, the 3D laser point cloud data is collected back and forth in the fast and slow lanes; if it is a completely closed two-way tunnel on a highway, the 3D laser point cloud data is collected back and forth in the fast lanes in both directions; if it is a semi-closed one-way tunnel on a highway, the 3D laser point cloud data is collected back and forth in the closed lanes; if it is a semi-closed two-way tunnel on a highway, the 3D laser point cloud data is collected back and forth in the closed lanes in all directions.

使用车载激光扫描系统单程、连续地扫描待测高速公路隧道全程可以得到单批次的三维激光点云数据;为了保证三维激光点云数据的覆盖完整性,至少需要对待测高速公路隧道进行两次单程连续的扫描,以采集多批次的三维激光点云数据。A single batch of 3D laser point cloud data can be obtained by using a vehicle-mounted laser scanning system to continuously scan the entire highway tunnel to be tested in a single pass. In order to ensure the coverage integrity of the 3D laser point cloud data, at least two single-pass continuous scans are required for the highway tunnel to be tested to collect multiple batches of 3D laser point cloud data.

S2:去除步骤S1中获取的三维激光点云数据中的异常点和噪声点后,基于各标靶中心的实测坐标对三维激光点云数据进行点云配准。S2: After removing abnormal points and noise points in the three-dimensional laser point cloud data obtained in step S1, the three-dimensional laser point cloud data is aligned based on the measured coordinates of the center of each target.

采集得到的三维激光点云数据中可能存在与周围点具有显著不同特征或属性的异常点,和由干扰因素引起的噪声点,为了提高点云数据的质量和准确性,本实施例采用临近点搜索法,即在搜索半径内的所有点云数据中,将高程差值大于某一阈值的点认定为异常点和噪点,进行去除,对去除噪声点和异常点后的三维激光点云数据再进行点云配准。The collected 3D laser point cloud data may contain abnormal points with significantly different characteristics or attributes from the surrounding points, and noise points caused by interference factors. In order to improve the quality and accuracy of the point cloud data, this embodiment adopts a neighboring point search method, that is, among all the point cloud data within the search radius, points with elevation differences greater than a certain threshold are identified as abnormal points and noise points, which are removed, and the 3D laser point cloud data after removing the noise points and abnormal points is then subjected to point cloud registration.

本实施例点云配准工作分为点云粗配准和点云精配准,先根据使用全站仪测得的隧道坐标系下的实测标靶中心,对所有批次的三维激光点云数据进行点云粗配准,在配准点云的同时,为点云精配准提供基准点云数据;然后再根据基准点云数据对其余批次粗配准后的三维激光点云数据使用ICP算法进行点云精配准,进一步减小不同批次三维激光点云数据间的配准误差,提高点云数据质量。The point cloud registration work in this embodiment is divided into point cloud coarse registration and point cloud fine registration. First, according to the measured target center in the tunnel coordinate system measured by the total station, all batches of 3D laser point cloud data are coarsely registered. While registering the point cloud, reference point cloud data is provided for point cloud fine registration. Then, according to the reference point cloud data, the ICP algorithm is used to perform point cloud fine registration on the remaining batches of coarsely registered 3D laser point cloud data, so as to further reduce the registration error between different batches of 3D laser point cloud data and improve the quality of point cloud data.

其中,点云粗配准按以下步骤实施:Among them, the point cloud coarse registration is implemented in the following steps:

S21:在当前批次的三维激光点云数据中,人工标记若干个覆盖单个标靶的点云数据的矩形区域,构建所述矩形区域对应的三维统计直方图,并对所述三维统计直方图进行平均化处理,得到基准三维统计直方图。S21: In the current batch of three-dimensional laser point cloud data, manually mark several rectangular areas of point cloud data covering a single target, construct a three-dimensional statistical histogram corresponding to the rectangular area, and average the three-dimensional statistical histogram to obtain a benchmark three-dimensional statistical histogram.

本实施例先将人工标记的每个覆盖单个标靶的点云数据的矩形区域内的三维激光点云数据,投影至该矩形区域内的三维激光点云数据的法向量方向的平面上,形成点云强度图像,然后对点云强度图像中的图像灰度和梯度信息进行统计,生成三维统计直方图,再将这些矩形区域生成的三维统计直方图进行平均化处理,以消除偶然差异,得到基准三维统计直方图;含有灰度和梯度信息的三维统计直方图更能反映该矩形区域内点云数据的特征,提高相似度匹配的准确性,本实施例构建的三维统计直方图的示意图如图2所示,X轴灰度表示对点云强度图像中标记的矩形区域的图像像元灰度值的统计,灰度值分为[0,32)、[32,64)、[64,96)、[96,128)、[128,160)、[160,192)、[192,224)、[224-255]共8个数值区间。Y轴梯度表示点云强度图像中标记的外接矩形的图像像元梯度方向的统计,该梯度方向是某一像元梯度方向与图像中心点方向的顺时针差值,将梯度方向差值分为[0,45)、[45,90)、[90,135)、[135,180)、[180,225)、[225,270)、[270,315)、[315,360)共8个数值区间。Z轴数量表示点云强度图像中对应的某一灰度区间和梯度区间的像元个数。In this embodiment, the three-dimensional laser point cloud data within each rectangular area of the point cloud data covering a single target is first projected onto a plane in the normal vector direction of the three-dimensional laser point cloud data within the rectangular area to form a point cloud intensity image, and then the image grayscale and gradient information in the point cloud intensity image are statistically analyzed to generate a three-dimensional statistical histogram. The three-dimensional statistical histograms generated by these rectangular areas are then averaged to eliminate accidental differences and obtain a reference three-dimensional statistical histogram. The three-dimensional statistical histogram containing grayscale and gradient information can better reflect the characteristics of the point cloud data within the rectangular area and improve the accuracy of similarity matching. A schematic diagram of the three-dimensional statistical histogram constructed in this embodiment is shown in FIG2 , where the grayscale on the X-axis represents the statistics of the grayscale values of the image pixels in the rectangular area marked in the point cloud intensity image, and the grayscale values are divided into 8 numerical intervals: [0,32), [32,64), [64,96), [96,128), [128,160), [160,192), [192,224), and [224-255]. The Y-axis gradient represents the statistics of the image pixel gradient direction of the circumscribed rectangle marked in the point cloud intensity image. The gradient direction is the clockwise difference between the gradient direction of a certain pixel and the direction of the image center point. The gradient direction difference is divided into 8 value intervals: [0,45), [45,90), [90,135), [135,180), [180,225), [225,270), [270,315), [315,360). The Z-axis quantity represents the number of pixels in a corresponding grayscale interval and gradient interval in the point cloud intensity image.

S22:按预设步距将未标记区域的三维激光点云数据,标记为与步骤S21中矩形区域尺寸相同的若干待匹配的矩形区域,构建待匹配的矩形区域内三维激光点云数据的三维统计直方图,并按照预设的相似度阈值逐一与步骤S21中构建的基准三维统计直方图进行相似度匹配,本实施例所设步距为0.05米,相似度阈值为0.9,该相似度阈值以上的矩形区域即为覆盖单个标靶的点云数据的矩形区域,这些矩形区域内的三维激光点云数据的几何质心,即为本实施例中标靶的中心点坐标。S22: Mark the three-dimensional laser point cloud data of the unmarked area as a number of rectangular areas to be matched with the same size as the rectangular area in step S21 according to a preset step size, construct a three-dimensional statistical histogram of the three-dimensional laser point cloud data in the rectangular area to be matched, and perform similarity matching with the reference three-dimensional statistical histogram constructed in step S21 one by one according to a preset similarity threshold. The step size set in this embodiment is 0.05 meters, and the similarity threshold is 0.9. The rectangular area above the similarity threshold is the rectangular area of the point cloud data covering a single target, and the geometric centroid of the three-dimensional laser point cloud data in these rectangular areas is the center point coordinate of the target in this embodiment.

S23:将S22中得到的标靶的中心点坐标纠正至与该标靶中心的实测坐标相同,得到该批次点云粗配准后的三维激光点云数据。本实施例中,针对完全封闭的高速公路单向隧道,在横向截面方向上,点云粗配准后的三维激光点云数据如图3所示。S23: Correct the coordinates of the center point of the target obtained in S22 to be the same as the measured coordinates of the center of the target, and obtain the three-dimensional laser point cloud data after the rough registration of the batch of point clouds. In this embodiment, for a completely closed one-way tunnel of a highway, in the transverse cross-sectional direction, the three-dimensional laser point cloud data after the rough registration of the point cloud is shown in FIG3 .

点云粗配准后,三维激光点云数据中标靶的中心点坐标与该标靶在隧道空间坐标系中的实测坐标仍存在一定误差,且不同标靶的误差也不同,本实施例通过计算所有批次粗配准后的三维激光点云数据的中误差,选取中误差最小时对应批次的三维激光点云数据作为基准点云数据,采用ICP算法将其余批次的三维激光点云数据配准至基准点云数据处,完成点云精配准。After the point cloud is roughly aligned, there is still a certain error between the coordinates of the center point of the target in the 3D laser point cloud data and the measured coordinates of the target in the tunnel space coordinate system, and the errors of different targets are also different. In this embodiment, the mean error of all batches of 3D laser point cloud data after rough alignment is calculated, and the 3D laser point cloud data of the corresponding batch with the smallest mean error is selected as the reference point cloud data. The ICP algorithm is used to align the remaining batches of 3D laser point cloud data to the reference point cloud data to complete the point cloud precise alignment.

中误差计算公式如下:The calculation formula of mean error is as follows:

;

其中,RMSE为点云数据的中误差,n为标靶数量,xi为三维激光点云数据中标靶的中心点坐标与实测的标靶中心点坐标之差,为该批次三维激光点云数据中所有标靶的中心点坐标与实测中心点坐标之差的平均值。Among them, RMSE is the mean error of the point cloud data, n is the number of targets, xi is the difference between the center point coordinates of the target in the 3D laser point cloud data and the measured center point coordinates of the target, It is the average value of the difference between the center point coordinates of all targets in this batch of 3D laser point cloud data and the measured center point coordinates.

ICP算法(Iterative Closest Point,最近点迭代算法),其核心思想是通过迭代的方式,不断优化两个点云之间的中误差,从而找到最佳的点云配准结果。The core idea of the ICP algorithm (Iterative Closest Point) is to continuously optimize the mean error between two point clouds through iteration to find the best point cloud registration result.

具体地,本实施例中首先在待配准点云数据中选择初始点集Pi,并在基准点云数据中搜索离Pi最近的点,形成点集Qi,使||Qi-Pi||最小,其中||·||表示两点集的欧几里得距离;计算点集Pi变换到点集Qi所需的旋转矩阵R和平移矩阵t,利用旋转矩阵R和平移矩阵t对待纠正点云数据进行旋转平移变换,得到纠正后的点云Pi';再次计算Pi'与Qi的中误差,若中误差≤设定阈值,则精配准完成,否则,将Pi'作为初始点集,重新进行上述迭代计算,直至中误差≤设定阈值或超过最大迭代次数;本实施例中设定阈值为1.5cm,针对完全封闭的高速公路单向隧道,在横向截面方向上,点云精配准后的三维激光点云数据如图4所示。Specifically, in this embodiment, an initial point set Pi is first selected from the point cloud data to be registered, and the point closest to Pi is searched in the reference point cloud data to form a point set Qi , so that || Qi - Pi || is minimized, where ||·|| represents the Euclidean distance between the two point sets; the rotation matrix R and translation matrix t required to transform the point set Pi to the point set Qi are calculated, and the rotation matrix R and translation matrix t are used to perform rotation and translation transformation on the point cloud data to be corrected to obtain the corrected point cloud Pi '; the mean square error between Pi ' and Qi is calculated again, and if the mean square error is ≤ the set threshold, the precise registration is completed, otherwise, Pi ' is used as the initial point set, and the above iterative calculation is performed again until the mean square error is ≤ the set threshold or the maximum number of iterations is exceeded; in this embodiment, the threshold is set to 1.5 cm, and for a completely closed one-way tunnel of a highway, the three-dimensional laser point cloud data after precise registration of the point cloud in the transverse section direction is shown in Figure 4.

S3:对点云配准后的三维激光点云数据进行分割和滤波,并生成地面不规则三角网和衬砌不规则三角网。S3: Segment and filter the 3D laser point cloud data after point cloud registration, and generate ground irregular triangulated network and lining irregular triangulated network.

通常在1m至1.5m处将点云配准后的三维激光点云数据分割为下部地面点云数据和上部衬砌点云数据;为了提高净空信息测量的准确性,需要对下部地面点云数据和上部衬砌点云数据进行滤波,以提取下部地面点云数据的表层点云数据和上部衬砌点云数据的表层点云数据。Usually, the 3D laser point cloud data after point cloud registration is divided into lower ground point cloud data and upper lining point cloud data at 1m to 1.5m; in order to improve the accuracy of clearance information measurement, the lower ground point cloud data and the upper lining point cloud data need to be filtered to extract the surface point cloud data of the lower ground point cloud data and the surface point cloud data of the upper lining point cloud data.

具体地,本实施例基于渐进三角网加密滤波算法提取下部地面点云数据的表层点云数据和上部衬砌点云数据的表层点云数据。该算法基本思想是从初始的地面点开始,逐步将离散的点云数据组成不规则三角网。本实施例首先将整体点云数据分块,以每个块中的高程值最低的点作为初始的地面种子点,并以此为基础构建初始的不规则三角网。接着以这个初始的不规则三角网为基础,对剩余的全部点云数据进行迭代加密。即首先找到待判断的点在xOy平面上的投影坐标所对应的三角面,计算该三角面的坡度,如果坡度大于预设的地形倾角阈值,则使用该点相对于该三角形内的最高点的镜像点来判断该点的类别;如果三角面的坡度小于预设的地形倾角阈值,则计算待判断点到该三角面的距离,以及该待判断点与距其最近的三角面的顶点的连线与该三角面之间的夹角角度,如果距离与夹角均小于预设的距离阈值和角度阈值,则该点视为表层点云数据,否则为非表层点云数据。在构建新的不规则三角网之后重复上述过程,直至没有新的待判断点加入或达到最大的迭代次数为止。如果某个点被判别为表层点云数据,则将其加入不规则三角网;如果被判别为非表层点云数据,则将其滤除。Specifically, this embodiment extracts the surface point cloud data of the lower ground point cloud data and the surface point cloud data of the upper lining point cloud data based on the progressive triangulated network encryption filtering algorithm. The basic idea of the algorithm is to start from the initial ground point and gradually form the discrete point cloud data into an irregular triangulated network. This embodiment first divides the overall point cloud data into blocks, takes the point with the lowest elevation value in each block as the initial ground seed point, and constructs the initial irregular triangulated network based on this. Then, based on this initial irregular triangulated network, all the remaining point cloud data are iteratively encrypted. That is, first find the triangular face corresponding to the projection coordinates of the point to be judged on the xOy plane, calculate the slope of the triangular face, if the slope is greater than the preset terrain inclination threshold, use the mirror point of the point relative to the highest point in the triangle to determine the category of the point; if the slope of the triangular face is less than the preset terrain inclination threshold, calculate the distance from the point to be judged to the triangular face, and the angle between the line connecting the point to be judged and the vertex of the nearest triangular face and the triangular face. If the distance and the angle are both less than the preset distance threshold and angle threshold, the point is regarded as surface point cloud data, otherwise it is non-surface point cloud data. After constructing a new irregular triangulated network, repeat the above process until no new points to be judged are added or the maximum number of iterations is reached. If a point is judged as surface point cloud data, it is added to the irregular triangulated network; if it is judged as non-surface point cloud data, it is filtered out.

对于下部地面点云数据,由于隧道地面可能存在标志标牌、运营车辆、杂物等非表层点云数据,因此需要使用渐进三角网加密滤波算法对其进行滤波分类;首先将下部地面点云数据分块,选取每一块中高程值最低的点作为种子点建立初始稀疏三角网,设定距离阈值与角度阈值,逐渐加密三角网,滤除非表层点云数据,提取出下部地面点云数据的表层点云数据。For the lower ground point cloud data, since there may be non-surface point cloud data such as signs, operating vehicles, and debris on the tunnel ground, it is necessary to use a progressive triangulation encryption filtering algorithm to filter and classify it; first, the lower ground point cloud data is divided into blocks, and the point with the lowest elevation value in each block is selected as the seed point to establish an initial sparse triangulation network, set the distance threshold and angle threshold, gradually encrypt the triangulation network, filter out non-surface point cloud data, and extract the surface point cloud data of the lower ground point cloud data.

对于上部衬砌点云数据,由于隧道顶部有通风、照明设施,若隧道衬砌已开展病害处置,还可能有工字钢钢架等装置,因此需要先根据待测隧道的顶部设施的情况设定合理的高程阈值,从上部衬砌点云数据中,提取低于预设高程阈值的点云数据作为低点点云簇,再使用渐进三角网加密滤波算法提取上部衬砌点云数据的底部表层点云数据:首先将上部衬砌点云数据分块,选取每一块中高程值最低的点作为种子点建立初始稀疏三角网,设定距离阈值与角度阈值,逐渐加密三角网,滤除非表层点云数据,提取出上部衬砌点云数据的底部表层点云数据。最后将提取得到的低点点云簇与底部表层点云数据合并,得到上部衬砌点云数据的表层点云数据。For the upper lining point cloud data, since there are ventilation and lighting facilities on the top of the tunnel, if the tunnel lining has been treated for defects, there may also be I-beam steel frames and other devices. Therefore, it is necessary to first set a reasonable elevation threshold according to the conditions of the top facilities of the tunnel to be tested, and extract the point cloud data below the preset elevation threshold from the upper lining point cloud data as the low-point point cloud cluster, and then use the progressive triangulation encryption filtering algorithm to extract the bottom surface point cloud data of the upper lining point cloud data: first, the upper lining point cloud data is divided into blocks, and the point with the lowest elevation value in each block is selected as the seed point to establish the initial sparse triangulation network, set the distance threshold and angle threshold, gradually encrypt the triangulation network, filter out the non-surface point cloud data, and extract the bottom surface point cloud data of the upper lining point cloud data. Finally, the extracted low-point point cloud cluster is merged with the bottom surface point cloud data to obtain the surface point cloud data of the upper lining point cloud data.

根据下部地面点云数据的表层点云数据和上部衬砌点云数据的表层点云数据即可分别构建地面不规则三角网和衬砌不规则三角网。According to the surface point cloud data of the lower ground point cloud data and the surface point cloud data of the upper lining point cloud data, the ground irregular triangulated network and the lining irregular triangulated network can be constructed respectively.

本实施例使用Delaunay三角网生长算法分别构建地面不规则三角网和衬砌不规则三角网,首先在提取出的表层点云数据中任取一点作为起点,连接与该点最近的一个点作为初始基线;在初始基线的右侧,查找与基线距离最短的第三个点,该点应满足自身与基线两个端点形成的三角形的外接圆内不包含其它点;以初始基线的两个端点和找到的第三点作为顶点生成一个Delaunay三角形;以新生成的三角网的两边作为新的基线继续生成三角形,直到表层点云数据中所有的点都被处理完毕,由这些生成的三角形构成了不规则三角网,本实施例中,针对完全封闭的高速公路单向隧道,所生成的地面不规则三角网和衬砌不规则三角网分别如图5和图6所示。This embodiment uses the Delaunay triangulation growth algorithm to construct a ground irregular triangulation network and a lining irregular triangulation network respectively. First, any point in the extracted surface point cloud data is selected as the starting point, and the point closest to the point is connected as the initial baseline; on the right side of the initial baseline, the third point with the shortest distance to the baseline is found, and the point should satisfy that the circumscribed circle of the triangle formed by itself and the two endpoints of the baseline does not contain other points; a Delaunay triangle is generated with the two endpoints of the initial baseline and the third point found as vertices; triangles are continuously generated with the two sides of the newly generated triangulation network as new baselines until all points in the surface point cloud data are processed, and an irregular triangulation network is formed by these generated triangles. In this embodiment, for a completely closed one-way tunnel of a highway, the generated ground irregular triangulation network and lining irregular triangulation network are shown in Figures 5 and 6 respectively.

S4:将待测高速公路隧道的标线线位的连续线位数据或设计线位的连续线位数据转化为点数据,所述点数据所在位置的衬砌不规则三角网的高程值与地面不规则三角网的高程值作差,得到所述点数据所在位置的净空值。S4: Convert the continuous line position data of the marking line position or the continuous line position data of the design line position of the highway tunnel to be tested into point data, and subtract the elevation value of the lining irregular triangulated network at the location of the point data from the elevation value of the ground irregular triangulated network to obtain the clearance value at the location of the point data.

通过表层点云数据生成的衬砌不规则三角网和地面不规则三角网的顶点是自带高程值的,因此在测量标线线位或设计线位上某一点的净空值时,只需要在不规则三角网中找到包含该点的三角形,基于该点所对应的三角形的顶点的高程值,使用线性插值法计算得到该点在衬砌不规则三角网和地面不规则三角网中的高程值,该点对应的衬砌不规则三角网的高程值减去该点对应的地面不规则三角网的高程值,即为该点所在位置的净空值,将待测高速公路隧道的标线线位的连续线位数据或设计线位的连续线位数据转化为点数据,再通过本发明方法测量所有标线线位或设计线位上的离散的点数据的净空值,就能得到待测高速公路隧道的完整净空信息。The vertices of the lining irregular triangulated network and the ground irregular triangulated network generated by the surface point cloud data have their own elevation values. Therefore, when measuring the clearance value of a certain point on the marking line position or the design line position, it is only necessary to find the triangle containing the point in the irregular triangulated network. Based on the elevation values of the vertices of the triangle corresponding to the point, the elevation value of the point in the lining irregular triangulated network and the ground irregular triangulated network is calculated using the linear interpolation method. The elevation value of the lining irregular triangulated network corresponding to the point minus the elevation value of the ground irregular triangulated network corresponding to the point is the clearance value of the location of the point. The continuous line position data of the marking line position or the continuous line position data of the design line position of the highway tunnel to be measured is converted into point data, and then the clearance values of the discrete point data on all marking line positions or design line positions are measured by the method of the present invention, so that the complete clearance information of the highway tunnel to be measured can be obtained.

实施例2Example 2

本实施例提供一种基于车载激光点云的高速公路隧道净空测量装置,包括至少一个处理器,以及与所述至少一个处理器通信连接的存储器;所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行实施例1中的基于车载激光点云的高速公路隧道净空测量方法。This embodiment provides a highway tunnel clearance measurement device based on vehicle-mounted laser point cloud, including at least one processor and a memory communicatively connected to the at least one processor; the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the highway tunnel clearance measurement method based on vehicle-mounted laser point cloud in Example 1.

以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims (10)

1.基于车载激光点云的高速公路隧道净空测量方法,其特征在于,包括如下步骤:1. A highway tunnel clearance measurement method based on vehicle-mounted laser point cloud, characterized in that it includes the following steps: S1:在待测高速公路隧道内布设若干标靶,测量得到各标靶中心的实测坐标,并使用车载激光扫描系统获取多批次的所述待测高速公路隧道的三维激光点云数据,所述三维激光点云数据包含所述标靶的点云数据;S1: placing a number of targets in the highway tunnel to be measured, measuring the measured coordinates of the center of each target, and using a vehicle-mounted laser scanning system to obtain multiple batches of three-dimensional laser point cloud data of the highway tunnel to be measured, wherein the three-dimensional laser point cloud data includes the point cloud data of the targets; S2:去除所述三维激光点云数据中的异常点和噪声点后,基于所述各标靶中心的实测坐标对所述三维激光点云数据进行点云配准;S2: after removing abnormal points and noise points in the three-dimensional laser point cloud data, performing point cloud registration on the three-dimensional laser point cloud data based on the measured coordinates of the centers of the targets; S3:对点云配准后的三维激光点云数据进行分割和滤波,并生成地面不规则三角网和衬砌不规则三角网;S3: segment and filter the 3D laser point cloud data after point cloud registration, and generate ground irregular triangulated network and lining irregular triangulated network; S4:将待测高速公路隧道的标线线位的连续线位数据或设计线位的连续线位数据转化为点数据,所述点数据所在位置的衬砌不规则三角网的高程值与地面不规则三角网的高程值作差,得到所述点数据所在位置的净空值。S4: Convert the continuous line position data of the marking line position or the continuous line position data of the design line position of the highway tunnel to be tested into point data, and subtract the elevation value of the lining irregular triangulated network at the location of the point data from the elevation value of the ground irregular triangulated network to obtain the clearance value at the location of the point data. 2.如权利要求1所述的基于车载激光点云的高速公路隧道净空测量方法,其特征在于,所述点云配准包括对不同批次的三维激光点云数据进行点云粗配准,和对点云粗配准后的三维激光点云数据进行点云精配准。2. The highway tunnel clearance measurement method based on vehicle-mounted laser point cloud as described in claim 1 is characterized in that the point cloud alignment includes coarse alignment of different batches of three-dimensional laser point cloud data and fine alignment of the three-dimensional laser point cloud data after the coarse alignment of the point clouds. 3.如权利要求2所述的基于车载激光点云的高速公路隧道净空测量方法,其特征在于,所述点云粗配准包括:3. The highway tunnel clearance measurement method based on vehicle-mounted laser point cloud according to claim 2, characterized in that the point cloud coarse registration comprises: S21:在当前批次的三维激光点云数据中,人工标记若干个覆盖单个标靶的点云数据的矩形区域,构建所述矩形区域对应的三维统计直方图,并对所述三维统计直方图进行平均化处理,得到基准三维统计直方图;S21: manually marking a number of rectangular areas of point cloud data covering a single target in the current batch of three-dimensional laser point cloud data, constructing a three-dimensional statistical histogram corresponding to the rectangular area, and averaging the three-dimensional statistical histogram to obtain a reference three-dimensional statistical histogram; S22:按预设步距将未标记区域的三维激光点云数据,标记为与步骤S21中矩形区域尺寸相同的若干待匹配的矩形区域,构建所述待匹配的矩形区域对应的三维统计直方图,并与基准三维统计直方图进行相似度匹配,获取所有覆盖单个标靶的点云数据的矩形区域,以得到所有标靶的中心点坐标;其中,标靶的中心点坐标为所述覆盖单个标靶的点云数据的矩形区域内的三维激光点云数据的几何质心;S22: marking the three-dimensional laser point cloud data of the unmarked area into a plurality of rectangular areas to be matched with the same size as the rectangular area in step S21 according to a preset step length, constructing a three-dimensional statistical histogram corresponding to the rectangular area to be matched, and performing similarity matching with the reference three-dimensional statistical histogram, obtaining all rectangular areas of point cloud data covering a single target, so as to obtain the center point coordinates of all targets; wherein the center point coordinates of the target are the geometric centroid of the three-dimensional laser point cloud data within the rectangular area of point cloud data covering a single target; S23:将各个所述标靶的中心点坐标纠正至与该标靶中心的实测坐标相同,得到该批次点云粗配准后的三维激光点云数据。S23: Correcting the coordinates of the center point of each target to be the same as the measured coordinates of the center of the target, and obtaining the three-dimensional laser point cloud data after the rough registration of the batch of point clouds. 4.如权利要求3所述的基于车载激光点云的高速公路隧道净空测量方法,其特征在于,所述三维统计直方图的构建方法包括:将每个矩形区域内的三维激光点云数据投影至该矩形区域内的三维激光点云数据的法向量方向的平面上,形成点云强度图像,生成含有灰度和梯度信息的三维统计直方图。4. The highway tunnel clearance measurement method based on vehicle-mounted laser point cloud as described in claim 3 is characterized in that the method for constructing the three-dimensional statistical histogram includes: projecting the three-dimensional laser point cloud data in each rectangular area onto a plane in the direction of the normal vector of the three-dimensional laser point cloud data in the rectangular area to form a point cloud intensity image, and generating a three-dimensional statistical histogram containing grayscale and gradient information. 5.如权利要求2所述的基于车载激光点云的高速公路隧道净空测量方法,其特征在于,所述点云精配准包括:计算所有批次点云粗配准后的三维激光点云数据的中误差,选取中误差最小时对应批次的三维激光点云数据作为基准点云数据,采用ICP算法将其余批次的三维激光点云数据配准至所述基准点云数据处,得到点云配准后的三维激光点云数据。5. The highway tunnel clearance measurement method based on vehicle-mounted laser point cloud as described in claim 2 is characterized in that the point cloud fine registration includes: calculating the mean square error of the three-dimensional laser point cloud data after the rough registration of all batches of point clouds, selecting the three-dimensional laser point cloud data of the corresponding batch when the mean square error is the smallest as the reference point cloud data, and using the ICP algorithm to register the remaining batches of three-dimensional laser point cloud data to the reference point cloud data to obtain the three-dimensional laser point cloud data after the point cloud registration. 6.如权利要求1-5任一所述的基于车载激光点云的高速公路隧道净空测量方法,其特征在于,所述S3包括:将所述点云配准后的三维激光点云数据分割为下部地面点云数据和上部衬砌点云数据;分别提取所述下部地面点云数据的表层点云数据和上部衬砌点云数据的表层点云数据,根据所述表层点云数据分别计算生成地面不规则三角网和衬砌不规则三角网。6. The method for measuring clearance of a highway tunnel based on a vehicle-mounted laser point cloud as described in any one of claims 1 to 5, characterized in that S3 comprises: dividing the three-dimensional laser point cloud data after the point cloud registration into lower ground point cloud data and upper lining point cloud data; respectively extracting the surface point cloud data of the lower ground point cloud data and the surface point cloud data of the upper lining point cloud data, and respectively calculating and generating a ground irregular triangulated network and a lining irregular triangulated network according to the surface point cloud data. 7.如权利要求6所述的基于车载激光点云的高速公路隧道净空测量方法,其特征在于,所述下部地面点云数据的表层点云数据使用渐进三角网加密滤波算法进行提取。7. The highway tunnel clearance measurement method based on vehicle-mounted laser point cloud as described in claim 6 is characterized in that the surface point cloud data of the lower ground point cloud data is extracted using a progressive triangulation encryption filtering algorithm. 8.如权利要求6所述的基于车载激光点云的高速公路隧道净空测量方法,其特征在于,所述上部衬砌点云数据的表层点云数据提取方法,包括:从上部衬砌点云数据中,提取低于预设高程阈值的点云数据作为低点点云簇,再使用渐进三角网加密滤波算法提取上部衬砌点云数据的底部表层点云数据,将所述低点点云簇与所述底部表层点云数据合并,得到上部衬砌点云数据的表层点云数据。8. The method for measuring clearance of a highway tunnel based on a vehicle-mounted laser point cloud as claimed in claim 6 is characterized in that the method for extracting surface point cloud data of the upper lining point cloud data comprises: extracting point cloud data below a preset elevation threshold from the upper lining point cloud data as a low-point point cloud cluster, then using a progressive triangulation encryption filtering algorithm to extract bottom surface point cloud data of the upper lining point cloud data, merging the low-point point cloud cluster with the bottom surface point cloud data to obtain surface point cloud data of the upper lining point cloud data. 9.如权利要求6所述的基于车载激光点云的高速公路隧道净空测量方法,其特征在于,所述地面不规则三角网和衬砌不规则三角网采用Delaunay三角网生长算法生成。9. The highway tunnel clearance measurement method based on vehicle-mounted laser point cloud according to claim 6, characterized in that the ground irregular triangulation network and lining irregular triangulation network are generated using a Delaunay triangulation growth algorithm. 10.基于车载激光点云的高速公路隧道净空测量装置,其特征在于,包括至少一个处理器,以及与所述至少一个处理器通信连接的存储器;所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1至9中任一项所述的基于车载激光点云的高速公路隧道净空测量方法。10. A highway tunnel clearance measurement device based on vehicle-mounted laser point cloud, characterized in that it includes at least one processor and a memory communicatively connected to the at least one processor; the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the highway tunnel clearance measurement method based on vehicle-mounted laser point cloud described in any one of claims 1 to 9.
CN202410439311.1A 2024-04-12 2024-04-12 Expressway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud Active CN118031904B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410439311.1A CN118031904B (en) 2024-04-12 2024-04-12 Expressway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410439311.1A CN118031904B (en) 2024-04-12 2024-04-12 Expressway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud

Publications (2)

Publication Number Publication Date
CN118031904A true CN118031904A (en) 2024-05-14
CN118031904B CN118031904B (en) 2024-06-25

Family

ID=90997219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410439311.1A Active CN118031904B (en) 2024-04-12 2024-04-12 Expressway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud

Country Status (1)

Country Link
CN (1) CN118031904B (en)

Citations (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102564393A (en) * 2011-12-28 2012-07-11 北京工业大学 Method for monitoring and measuring full section of tunnel through three-dimensional laser
CN102980531A (en) * 2012-12-07 2013-03-20 中国铁道科学研究院 Volume measurement method and device based on three-dimensional laser scanning
CN104075691A (en) * 2014-07-09 2014-10-01 广州市城市规划勘测设计研究院 Method for quickly measuring topography by using ground laser scanner based on CORS (Continuous Operational Reference System) and ICP (Iterative Closest Point) algorithms
CN105574929A (en) * 2015-12-15 2016-05-11 电子科技大学 Single vegetation three-dimensional modeling method based on ground LiDAR point cloud data
CN106401643A (en) * 2016-08-31 2017-02-15 铁道第三勘察设计院集团有限公司 Tunnel back-break detection method based on laser-point cloud
CN108389250A (en) * 2018-03-08 2018-08-10 武汉大学 The method for quickly generating building cross-section diagram based on point cloud data
CN108917712A (en) * 2018-07-10 2018-11-30 湖南城市学院 A kind of Tunnel automation monitoring system and method based on three-dimensional laser scanning technique
CN109459439A (en) * 2018-12-06 2019-03-12 东南大学 A kind of Tunnel Lining Cracks detection method based on mobile three-dimensional laser scanning technique
CN109631786A (en) * 2018-12-14 2019-04-16 青岛理工大学 Three-dimensional laser scanning underground engineering similar material simulation test surface layer deformation method
CN110244321A (en) * 2019-04-22 2019-09-17 武汉理工大学 A detection method of road passable area based on 3D lidar
US20190318536A1 (en) * 2016-06-20 2019-10-17 Ocean Maps GmbH Method for Generating 3D Data Relating to an Object
CN110390687A (en) * 2019-07-29 2019-10-29 四川大学 A method for measuring river erosion and deposition based on 3D laser scanning
US20200019815A1 (en) * 2018-07-16 2020-01-16 Here Global B.V. Clustering for k-anonymity in location trajectory data
WO2020114466A1 (en) * 2018-12-05 2020-06-11 中国铁建重工集团股份有限公司 Tunnel point cloud data analysis method and system
CN211783409U (en) * 2020-03-26 2020-10-27 中交公路规划设计院有限公司 Automatic measuring device for settlement of vault of tunnel
US20200364929A1 (en) * 2019-05-13 2020-11-19 Wuhan University Multi-story indoor structured three-dimensional modeling method and system
CN112254637A (en) * 2020-10-13 2021-01-22 成都天佑智隧科技有限公司 Tunnel excavation surface scanning device and detection method based on various fusion data
CN113012205A (en) * 2020-11-17 2021-06-22 浙江华云电力工程设计咨询有限公司 Three-dimensional reconstruction method based on multi-source data fusion
CN113280798A (en) * 2021-07-20 2021-08-20 四川省公路规划勘察设计研究院有限公司 Geometric correction method for vehicle-mounted scanning point cloud under tunnel GNSS rejection environment
CN114549879A (en) * 2022-04-25 2022-05-27 四川省公路规划勘察设计研究院有限公司 A method for target recognition and center point extraction of tunnel vehicle scanning point cloud
CN114638954A (en) * 2022-02-22 2022-06-17 深圳元戎启行科技有限公司 Point cloud segmentation model training method, point cloud data segmentation method and related device
CN115203778A (en) * 2022-05-17 2022-10-18 中铁二十局集团第三工程有限公司 Tunnel overbreak and underexcavation detection method and device, terminal equipment and storage medium
CN115343299A (en) * 2022-10-18 2022-11-15 山东大学 Lightweight highway tunnel integrated detection system and method
CN115657049A (en) * 2022-10-18 2023-01-31 山东大学 Tunnel vehicle-mounted laser radar positioning and deviation rectifying method and system
CN116659460A (en) * 2023-05-06 2023-08-29 中交第二公路勘察设计研究院有限公司 Rapid generation method for laser point cloud slice of road cross section

Patent Citations (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102564393A (en) * 2011-12-28 2012-07-11 北京工业大学 Method for monitoring and measuring full section of tunnel through three-dimensional laser
CN102980531A (en) * 2012-12-07 2013-03-20 中国铁道科学研究院 Volume measurement method and device based on three-dimensional laser scanning
CN104075691A (en) * 2014-07-09 2014-10-01 广州市城市规划勘测设计研究院 Method for quickly measuring topography by using ground laser scanner based on CORS (Continuous Operational Reference System) and ICP (Iterative Closest Point) algorithms
CN105574929A (en) * 2015-12-15 2016-05-11 电子科技大学 Single vegetation three-dimensional modeling method based on ground LiDAR point cloud data
US20190318536A1 (en) * 2016-06-20 2019-10-17 Ocean Maps GmbH Method for Generating 3D Data Relating to an Object
CN106401643A (en) * 2016-08-31 2017-02-15 铁道第三勘察设计院集团有限公司 Tunnel back-break detection method based on laser-point cloud
CN108389250A (en) * 2018-03-08 2018-08-10 武汉大学 The method for quickly generating building cross-section diagram based on point cloud data
CN108917712A (en) * 2018-07-10 2018-11-30 湖南城市学院 A kind of Tunnel automation monitoring system and method based on three-dimensional laser scanning technique
US20200019815A1 (en) * 2018-07-16 2020-01-16 Here Global B.V. Clustering for k-anonymity in location trajectory data
WO2020114466A1 (en) * 2018-12-05 2020-06-11 中国铁建重工集团股份有限公司 Tunnel point cloud data analysis method and system
CN109459439A (en) * 2018-12-06 2019-03-12 东南大学 A kind of Tunnel Lining Cracks detection method based on mobile three-dimensional laser scanning technique
CN109631786A (en) * 2018-12-14 2019-04-16 青岛理工大学 Three-dimensional laser scanning underground engineering similar material simulation test surface layer deformation method
CN110244321A (en) * 2019-04-22 2019-09-17 武汉理工大学 A detection method of road passable area based on 3D lidar
US20200364929A1 (en) * 2019-05-13 2020-11-19 Wuhan University Multi-story indoor structured three-dimensional modeling method and system
CN110390687A (en) * 2019-07-29 2019-10-29 四川大学 A method for measuring river erosion and deposition based on 3D laser scanning
CN211783409U (en) * 2020-03-26 2020-10-27 中交公路规划设计院有限公司 Automatic measuring device for settlement of vault of tunnel
CN112254637A (en) * 2020-10-13 2021-01-22 成都天佑智隧科技有限公司 Tunnel excavation surface scanning device and detection method based on various fusion data
CN113012205A (en) * 2020-11-17 2021-06-22 浙江华云电力工程设计咨询有限公司 Three-dimensional reconstruction method based on multi-source data fusion
CN113280798A (en) * 2021-07-20 2021-08-20 四川省公路规划勘察设计研究院有限公司 Geometric correction method for vehicle-mounted scanning point cloud under tunnel GNSS rejection environment
CN114638954A (en) * 2022-02-22 2022-06-17 深圳元戎启行科技有限公司 Point cloud segmentation model training method, point cloud data segmentation method and related device
CN114549879A (en) * 2022-04-25 2022-05-27 四川省公路规划勘察设计研究院有限公司 A method for target recognition and center point extraction of tunnel vehicle scanning point cloud
CN115203778A (en) * 2022-05-17 2022-10-18 中铁二十局集团第三工程有限公司 Tunnel overbreak and underexcavation detection method and device, terminal equipment and storage medium
CN115343299A (en) * 2022-10-18 2022-11-15 山东大学 Lightweight highway tunnel integrated detection system and method
CN115657049A (en) * 2022-10-18 2023-01-31 山东大学 Tunnel vehicle-mounted laser radar positioning and deviation rectifying method and system
CN116659460A (en) * 2023-05-06 2023-08-29 中交第二公路勘察设计研究院有限公司 Rapid generation method for laser point cloud slice of road cross section

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
WANG, KY 等: "Adaptively unsupervised seepage detection in tunnels from 3D point clouds", 《STRUCTURE AND INFRASTRUCTURE ENGINEERING》, no. 11, 18 November 2022 (2022-11-18), pages 15732479 *
张力学: "三维激光扫描技术在地铁隧道净空检测中的应用研究", 《市政技术》, vol. 39, no. 6, 15 June 2021 (2021-06-15), pages 94 - 99 *
张春森 等: "多视几何无人机影像堆体体积量算", 《西安科技大学学报》, vol. 39, no. 1, 15 January 2019 (2019-01-15), pages 124 - 129 *

Also Published As

Publication number Publication date
CN118031904B (en) 2024-06-25

Similar Documents

Publication Publication Date Title
CN107092020B (en) Road roughness monitoring method based on UAV LiDAR and high-resolution images
CN104075691B (en) A method for measuring topography with a terrestrial laser scanner
WO2021232463A1 (en) Multi-source mobile measurement point cloud data air-ground integrated fusion method and storage medium
CN109100719B (en) Terrain map joint mapping method based on satellite-borne SAR (synthetic aperture radar) image and optical image
CN109900338B (en) Method and device for measuring volume of pavement pit
CN105844629A (en) Automatic segmentation method for point cloud of facade of large scene city building
CN108845569A (en) Generate semi-automatic cloud method of the horizontal bend lane of three-dimensional high-definition mileage chart
CN104123730A (en) Method and system for remote-sensing image and laser point cloud registration based on road features
CN116793305A (en) A tunnel detection system and method based on three-dimensional laser scanning
CN112070756B (en) A Pavement Stereo-Disease Measurement Method Based on UAV Oblique Photography
CN114549879A (en) A method for target recognition and center point extraction of tunnel vehicle scanning point cloud
CN114283070B (en) Method for manufacturing terrain section by fusing unmanned aerial vehicle image and laser point cloud
CN105956542B (en) High-resolution remote sensing image road extraction method based on statistical matching of structural wire harnesses
CN114119863A (en) Method for automatically extracting street tree target and forest attribute thereof based on vehicle-mounted laser radar data
Zhou et al. Highway alignments extraction and 3D modeling from airborne laser scanning point clouds
CN114821522A (en) Urban road cross slope and super height value calculation method based on vehicle-mounted laser point cloud data
CN104729529B (en) The method and system that map surveying systematic error judges
JP4619504B2 (en) 3D digital map generator
CN118274784A (en) Road surface non-uniform settlement detection method and device based on high-density gray point cloud
CN111783192B (en) Complex terrain field flat earthwork calculation method based on oblique photography live-action model
Rufei et al. Research on a pavement pothole extraction method based on vehicle-borne continuous laser scanning point cloud
CN118395577B (en) A variable milling parameter design method and device for asphalt pavement
Gao et al. Automatic geo-referencing mobile laser scanning data to UAV images
CN118031904B (en) Expressway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud
CN116740307A (en) Smart city three-dimensional model construction method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant