1. Introduction
In recent years, efforts have been made to address the traffic congestion caused by urbanization. Major cities in China are accelerating the planning and construction of fast and high-capacity transportation modes such as subways and light rail. By the end of 2022, a total of 308 urban rail transit operation lines will be opened in 55 cities in mainland China, with a total length of 10,287.45 km [
1]. Of these, subway operating lines account for 77.84%, with a total length of 8008.17 km. Other standard urban rail transit operating lines account for 22.16%, with a total length of 2279.28 km. In 2022, the length of newly added operating routes was 1080.63 km, and the cumulative passenger volume for the year reached 19.302 billion people. The total length of the construction plan currently being implemented is 6675.57 km. By the end of the 14th Five Year Plan period, the scale of urban rail transit operation lines will be close to 13,000 km, with over 60 operating cities. The operation scale of urban rail transit continues to expand, and its role as a backbone in public transportation is becoming more apparent. As the role of subway and rail transit in Chinese cities has become increasingly important, its safety has become extremely important. The subway project has huge investment, complex construction, and a long construction period and is considered a century-old project. Therefore, during construction, completion, and operation and maintenance, the quality requirements are very strict. With the increase in the number of subways, the workload of subway surveys is enormous, and the time windows for subway surveys during operation are relatively short. The traditional survey methods include using the total station to detect ellipticity and using a ruler to detect platform dislocation. The total station is an improved version of the theodolite with auto-leveling and electronic rangefinder. However, these methods do not have automation functions and need manual surveys. In the process of measurement, the instrument needs to be moved and leveled frequently, which leads to low efficiency. Therefore, it is necessary to conduct research on the efficient, accurate, and automated daily maintenance and survey of subway tunnels to meet the requirements of subway construction and operation processes.
Currently, mobile 3D laser scanners are commonly used tools for tunnel surveys. These can quickly obtain the surface information of the scanning area and generate point cloud data. Point cloud data are a collection of discrete points that represent the surface information of objects in three-dimensional space. They include three-dimensional coordinates and laser reflection intensity for each point. In point cloud data, the X-axis of the coordinate system represents the horizontal direction of the tunnel cross-section, the Z-axis represents the vertical direction, and the Y-axis represents the direction in which the car moves forward (i.e., the mileage direction). The point cloud reflection intensity refers to the reflected energy intensity of each point recorded by a 3D laser scanner during the scanning process. When the laser beam is irradiated on the target surface, some of the light is reflected back by the target object. The laser scanner will measure and record the energy intensity of this reflected light, and the reflection intensity of different objects varies. The reflection intensity value is also affected by factors such as distance and equipment parameters. In the tunnel, the influence of atmospheric conditions on reflection intensity is temporarily not considered.
The information of daily maintenance and surveys of subway can be quickly extracted from point cloud data, which can effectively carry out subway surveys every day. To meet the accuracy requirements of various surveys, point cloud data need to be denoised. In most daily maintenance and inspection projects, such as the ovality and step of tunnels, only point cloud information of pipe segments is required. However, the tunnel point cloud data obtained via scanning include point cloud data on segments, pipelines, tracks, supports, and various ancillary equipment. At this time, for survey work, the point cloud information of pipelines, tracks, supports, and various ancillary equipment will be regarded as noise. Noise can have a significant impact on the analysis results, so it is necessary to remove the noise.
In addition to three-dimensional and intensity information, there is also more basic information in point cloud data, such as color information, point cloud density, and normal vectors. Currently, utilizing three-dimensional information for denoising is the most common approach, but usually, further processing is carried out on three-dimensional information. Zhu [
2] obtained curvature features from the three-dimensional information of point clouds and then performed denoising on the point clouds. Some authors consider point cloud density. Zhao [
3] proposed a point cloud denoising algorithm based on curvature and point cloud density. For points used for surface reconstruction, Zama [
4] proposed using a denoising method based on point cloud density. For ring forgings containing noise, Zhang [
5] extracted profile lines based on point cloud normal vectors. Liu [
6] proposed using point cloud normal vectors and curvatures to smooth out noise in hydraulic torque converters. These denoising methods all process the three-dimensional information of point clouds, requiring a large amount of computation for each point, resulting in a longer processing time.
At present, there are many denoising methods using 3D information, such as filter-based denoising, optimization-based denoising, and deep learning-based denoising. Wei [
7] uses the dual normal filtering algorithm to denoise irregular surfaces. Mattei [
8], Dai [
9], and Zheng [
10] all specialize in processing sharp features of point clouds, using MRPCA denoising, position filtering denoising, and guided normal filtering denoising, respectively. Liang [
11] uses an adaptive optimal neighborhood algorithm for denoising scattering point cloud models. Gao [
12] transformed the point cloud denoising problem into an optimization problem of the prior regularization of graphic signal smoothness. Sun [
13] uses L0 minimization for denoising, but the boundaries cannot be well handled. Some authors use denoising networks for denoising. Li [
14] targets unordered noise, while Yan [
15] targets noise in the engine model. The above denoising methods can meet the needs in their respective fields, but most methods only use three-dimensional coordinate information for denoising, while other basic information of point clouds is rarely used.
The core method of Irfan [
16] and Cao [
17] denoises through color information, but Irfan further combines geometric information to construct k-NN graphs for denoising. However, there are some limitations in denoising through color information, and it is difficult to remove the noise of the same color. Choi [
18] and Clode [
19] denoise by setting a threshold range for the reflection intensity value, but the calculation of the threshold range is still entirely based on the reflection intensity information. Most of these methods only use one piece of information of the point cloud for denoising. This paper also uses reflection intensity to denoise, so it is necessary to study the problems existing in intensity denoising. Although only using the threshold value of reflection intensity can remove some noises, there are still many noises with similar reflection intensity values that cannot be removed.
Faced with this situation, we propose a denoising method for tunnel point clouds, which includes reflection intensity threshold denoising, joint shape and reflection intensity denoising, and shape denoising. Among them, the joint shape and reflection intensity denoising method comprehensively considers the correlation between the tunnel geometry and reflection intensity of the point cloud. This method is also suitable for other situations where the strength value changes uniformly, such as the measurement of a circular shield working shaft. This proposed method aims to provide high-quality point clouds and provide a reliable basis for the subsequent analysis of and application to tunnels.
Before denoising the point cloud data of shield tunnels in this paper, preprocessing is required. The first step in preprocessing is to segment each segment ring from the point cloud data. A segment ring is a complete lining ring composed of multiple segments. This approach can break down a large computational task into multiple small tasks, thereby improving computational efficiency. Next, unfold the point cloud data of each segmented segment, and partition the point cloud data into multiple parts along the unfolding direction. Once the pre-processing of the point cloud data is complete, the denoising process begins. This paper has three steps for point cloud denoising. The first step is the threshold denoising method of reflection intensity. This step involves setting an appropriate reflection intensity threshold range to remove point cloud data with intensity values exceeding that threshold range. The next step is to denoise the shape and reflection intensity jointly. In this method, the intensity value of the unfolded point cloud data is processed by the optimized quartile range method and the Z score method. The last step is the geometric shape denoising method. This step involves setting a distance threshold based on the design radius of the segment and the fitted elliptical curve and removing point cloud data outside the distance threshold. By using the above steps to denoise the point cloud data of shield tunnels, the quality and accuracy of the data can be improved, providing a reliable foundation for subsequent analysis and application.
The organizational structure of this paper is as follows. In the
Section 1, the background and research questions are introduced. The
Section 2 provides a detailed introduction to the preprocessing and denoising steps of point cloud data. The feasibility of the algorithm was verified through practical engineering examples, which are outlined in the
Section 3. The
Section 4 verifies the denoising results of the algorithm. The
Section 5 provides some conclusions and future plans.
2. Denoising Methods
This section introduces each step of the denoising method in detail. First, before denoising, some processing of point clouds is beneficial to efficient denoising. The first step is to segment each segment ring from the point cloud data. The second step is to unfold the point cloud of the segment. Next, denoise the point cloud. The third step is to process point clouds by the reflection intensity threshold method. The fourth step is to combine the denoising methods of shape and reflection intensity to process the point cloud. The fifth step is the geometric shape denoising method to process the point cloud. The denoising flowchart is shown in
Figure 1.
2.1. Shield Tunnel Point Cloud Segmentation
The complete point cloud data of a shield tunnel usually comprise several hundred gigabytes, which is a large amount of data to process. Therefore, before denoising the point cloud, it is necessary to segment the point cloud data, which helps to reduce the complexity of data, improve processing efficiency, and better manage and utilize large-scale point cloud datasets. By segmenting point clouds, a series of data subsets can be obtained, making further processing, analysis, and application more effective. In shield tunnels, there are assembly gaps between adjacent ring segments. These gaps will cause changes in the scanning distance, leading to changes in the intensity value of point cloud data, as shown in
Figure 2. This paper utilizes the variation information of point cloud intensity values to accurately segment the point cloud of pipe segments. In order to avoid the impact of equipment such as pipelines, tracks, and signal boxes in the tunnel on point cloud segmentation, the point cloud data at the top of the tunnel are selected for analysis.
However, using point cloud data from the entire top of the tunnel for calculation can result in significant computational complexity. Therefore, some cloud points (6 degrees) are selected from the top of the tunnel to replace the point clouds at the top of the whole tunnel. Point cloud segmentation is calculated by the selected part of the point cloud, as shown in
Figure 3a. It can reduce computational complexity and improve work efficiency while ensuring accuracy.
Point cloud data are grouped in groups of 40 mm on the mileage (
Y axis), and the average value of the intensity values of all points in each group is calculated as the intensity value of this group of points, as shown in
Figure 3b. Calculate the difference between the intensity value of each group and its five groups of adjacent intensity values before and after. By comparing the differences, the intensity value that is much smaller than the adjacent intensity value is calculated. The grouping position of these smaller strength values is the assembling position between each ring segment. To verify the accuracy of the position, the distance between these smaller intensity values is calculated. If the separation distance is the width of segment, then the positions are accurate. The point cloud is segmented into single-ring segments at the assembly position, as shown in
Figure 3c.
2.2. Unfolding of Point Cloud of Shield Tunnel Segment
After dividing the point cloud data, in order to perform denoising processing, it is necessary to unfold the point cloud data. Because the mobile 3D laser scanner travels on the track of the tunnel, the center of the scanner is horizontally located on the axis of the tunnel. However, in the vertical direction, due to the height limitation of the scanner, the center of the scanner is not located on the axis of the tunnel, as shown in
Figure 4. Therefore, the
Y-axis in the coordinate system of the point cloud data does not coincide with the axis of the tunnel. Before unfolding, it is necessary to process the Z-values of all point clouds so that the
Y-axis in the coordinate system of the point cloud data coincides with the axis of the tunnel. In addition, ensure that the coordinate origin is located on the axis of the tunnel.
After processing the point cloud coordinates, project the point cloud onto the XOZ plane. The coordinates of any point P in the projected point cloud are (x, z). The radian is used as the unit of point cloud unfolding, which is an angular unit and dimensionless [
20]. During the unfolding process, the radian is based on the negative direction of the
Z-axis, and counterclockwise rotation is considered positive, namely,
, as shown in
Figure 5. In the figure, A is −π/2, which is approximately −1.57, and B is 3 π/2, which is approximately 4.71. The formula for calculating radians is shown in Equation (1). During the unfolding process of the tunnel point cloud, other values in the point cloud data remain unchanged. Due to the obstruction of the track, the point cloud of the bottom segment cannot be scanned, so the bottom point cloud is considered noise. Therefore, point clouds within the range of 0.5 rad to 8.5 rad were selected for denoising.
2.3. Reflection Intensity Threshold Denoising
After preprocessing the point cloud data, the point cloud is first denoised using the reflection intensity threshold method. The shield tunnel contains objects such as segments, pipelines, tracks, supports, and other equipment. They all have different materials, so their surface reflectivity and reflection intensity will also vary. By setting a threshold range for reflection intensity, only points within the threshold range are retained. The points that exceed the threshold are taken as noise and removed. When determining the threshold range of the point cloud reflection intensity, it is necessary to consider the distribution of the point cloud intensity and the reflection intensity of the segments and adjust it based on specific application scenarios and data characteristics. The threshold range can be found directly in the point cloud, through the different intensity values of points at different positions in the segment, or the distribution of intensity values in the point cloud is counted, and the threshold range can be selected according to the concentration range of intensity values.
However, relying solely on reflection intensity threshold denoising methods cannot remove all noise points. In the point cloud of a shield tunnel, noise may have reflection intensity values similar to the target point. So, the reflection intensity threshold denoising can only remove noise with obvious differences and cannot remove noise with similar intensity values. Therefore, when applying the threshold denoising method for reflection intensity values, it is necessary to comprehensively consider other denoising methods and data characteristics in order to achieve better denoising results.
2.4. Intensity and Geometry Combined Denoising
The joint denoising of the shape and reflection intensity is also a method of denoising through point cloud intensity values, which combines the geometric shape of point clouds and reflection intensity values. Due to the distance effect, objects of the same material may produce different reflection intensities, while objects of different materials may produce the same reflection intensity. This seriously affects the effectiveness of using the reflection intensity threshold for denoising. To address this issue, this method not only takes into account the variation in strength values of each point itself but also takes into account the geometric shape characteristics of the tunnel point cloud. For the geometric shape characteristics of tunnel point clouds, this paper considers the variation of strength values between each point and its adjacent points. When it is necessary to extract the same material point cloud from a point cloud, a situation may occur where the points are geometrically distant from each other, but the intensity values of the points are the same (the material and distance have a comprehensive impact on the intensity values).
For mobile 3D laser scanners in shield tunneling subway tunnels, the instrument center point during the scanning process is not at the center of the cross-section design circle. Therefore, when the laser head of the scanner rotates once, the distance from the laser to the segment surface at different angles changes uniformly. If the strength value of the same material is only affected by the change in scanning distance, and this change is uniform, the intensity value will also vary uniformly. Therefore, after unfolding the point cloud, the point cloud intensity value of the segment should exhibit a uniformly changing curve. When an abnormal intensity value occurs in a uniform variation of intensity values, the point corresponding to the abnormal intensity value is noise. This abnormal intensity value may be caused by two factors: first, different target materials, which lead to different reflectivities, and second, different distances. Based on this situation, this paper proposes a joint denoising approach of shape and reflection intensity. First, the point cloud is processed according to the method in
Section 2.2. After that, the point cloud data are divided into 200 groups of points according to a group of 0.04 rad in radian. Finally, the intensity values of each group of points are processed by IQR and Z-Score.
2.4.1. Interquartile Range
The IQR (interquartile range) is a robust method used to label outliers [
21]. The IQR method for detecting outliers was developed by John Tukey, who was a pioneer in exploratory data analysis. Arrange the dataset from small to large, and partition it into four equal parts; the values at the three dividing points are quartiles, namely, Q1 (the lower quartile), Q2 (the median), and Q3 (the upper quartile). The box-whisker plot uses quartiles to plot the shape of data, as shown in
Figure 6. The box represents Q1 and Q3, and the lines inside the box represent Q2, which is the median. The interquartile range (hence the name of the outlier detection method) is the distance between Q1 and Q3 (box edges), Q3 − Q1. Tukey considered data located outside Q3 + 1.5 × IQR(Upper 1.5 × IQR whisker) or Q1 − 1.5 × IQR(Lower 1.5 × IQR whisker) to be outliers. The whisker line extends to the last data point within the boundary.
In order to meet the needs of joint denoising of shape and reflection intensity, the box-whisker plot is improved, as shown in
Figure 7. First, set the size range for the box-whisker plot boundary. The next step is to detect the boundary value of the box-whisker plot. Finally, smooth and continuous processing is carried out between adjacent box-whisker plots. The detailed steps are as follows:
Calculate the box-whisker plot of each group of points: Calculate the interquartile range of each group of points, and calculate the box-whisker plot of each group of points according to the quartile, so that each group of points has a box-whisker plot.
Set the size range of the box-whisker plot boundary: Calculate the difference between Q3 + 1.5 × IQR and Q1 − 1.5 × IQR and set the boundary size range. In each group of points, the rotation of the scanner is small, so the laser distance changes little, and thus, the intensity value changes little. Therefore, the segment strength values in each group of points are concentrated near one value. By setting the maximum limit range and the minimum limit range, it can prevent the inclusion of abnormal values and the removal of target values.
Boundary value detection: Calculate the distance from the edge point to the nearest boundary, and if it is less than a certain value, take the edge point as the boundary value. This ensures that multiple filters do not delete all values, retaining more valid points and point cloud data accuracy.
Check the position of the box-whisker plot: It is necessary to check the position of each box-whisker plot to ensure that each box-whisker plot and its adjacent box-whisker plot meet the requirements of smoothness and continuity. When smoothing is detected, the slope of Q2 between adjacent boxes is calculated, and the slope threshold is set. If the slope exceeds the threshold, it needs to be smoothed. In this paper, the interpolation method is adopted to deal with it, and Q2 of the box-whisker plot exceeding the threshold is set as the interpolation of Q2 of the previous box-whisker plot and the next box-whisker plot. The boundary of the box-whisker plot changes with Q2. The reason why the positions of some adjacent box-whisker plot are quite different is that there are almost no segment points in this group, which are basically noise caused by obstructions.
2.4.2. Z-Score
The Z-score describes a data point in relation to the mean and standard deviation of a set of points [
22]. It maps the data to a distribution with a mean of 0 and a standard deviation of 1. The purpose of using Z-scores is to eliminate the influence of the data position and scale, allowing for direct comparison across different datasets. The logic behind the Z-score outlier detection method is that once the data are centered and rescaled, any values that deviate significantly from zero should be considered outliers. The Z-score method relies on the mean and standard deviation of a set of data to measure its central tendency and dispersion. However, since the mean and standard deviation are highly influenced by outliers, they are not robust measures. The presence of outliers causing skewness is one of the main reasons for identifying and removing outliers from a dataset. Another drawback of the Z-score method is its poor performance on small datasets. In fact, if the number of items in the dataset is less than 12, the Z-score method fails to detect outliers. This led to the development of the modified Z-score method, which overcomes these limitations, shown in Equation (2). This method is particularly suitable for skewed data or data with a non-normal distribution and a small number of observations. The MAD (median absolute deviation) is a robust measure of central tendency and dispersion used in the modified Z-score method. One advantage of the modified Z-score method is that it uses the median and MAD instead of the mean and standard deviation. The median and MAD are robust measures that provide a more robust assessment of central tendency and dispersion.
2.5. Geometric Denoising
Geometric shape denoising is a method based on the geometric features of shield tunnel point clouds, which denoises by setting different ranges of circles and ellipses [
23]. This method identifies and removes noise by using the distance between the point in the point cloud and the center of the tunnel. This method uses circles and ellipses as distance thresholds. Identify points within the distance threshold as valid points, and identify points outside the distance threshold as noise and remove them. By adjusting the size of circles and ellipses, the range of denoising can be controlled. First, the ellipse is fitted once to determine the center point of the tunnel segment. Afterward, by increasing and decreasing a certain value of the design radius, a circular range is formed, and points outside the circular range are treated as noise and removed. Next, fit the ellipse again to obtain the center point of the segment and the elliptical equation. By increasing and decreasing the values of the major and minor axes of the ellipse, two ellipses are obtained, forming a circular range. Similarly, treat points outside the torus as noise and remove them.
The advantages of geometric shape denoising methods are that they are simple, intuitive, and easy to implement. However, this method also has some limitations and precautions. The selection of the distance threshold needs to be adjusted according to the characteristics of specific data and application needs. If the threshold selection is too small, the target point may be erroneously removed. If the threshold selection is too large, it may not be possible to effectively remove noise points. Therefore, in this method, geometric shape denoising is used as a supplement to the intensity value denoising method. By comprehensively applying multiple denoising methods to process point clouds, better denoising effects and data quality can be achieved.
3. Experimental Results
In order to verify the practicality of the algorithm, this paper takes a shield tunnel in a certain section as an example to carry out the experiment. The tunnel had a total length of 575 m and was constructed using the shield tunneling method. It had an inner diameter of 5.4 m, an outer diameter of 6 m, and a segment width of 1.5 m. We used a mobile 3D laser scanner to obtain point cloud data of shield tunnel segments. We rendered the point cloud based on its intensity values, as shown in
Figure 8a,b. By analyzing the point cloud strength values of shield tunnel segments, we found that the strength values were mainly concentrated within the range of [100, 160], as shown in
Figure 8c.
Because the segmentation with complete point cloud data requires a huge amount of data, it will lead to extremely low computational efficiency. Point cloud segmentation helps to reduce the complexity of data and improve the processing efficiency. So, according to the method in
Section 2.1, some point cloud data above the top of the tunnel are selected for segmentation, and the selected data are shown in
Figure 9.
We projected the selected point cloud data onto the XOY plane, as shown in
Figure 10a, and drew a mileage intensity map of the point cloud, as shown in
Figure 10b. The horizontal axes in
Figure 10a,b are the same and can be observed correspondingly. Observation shows that with every 1.5 m change in mileage, the point cloud intensity value suddenly decreased. These reduced strength values corresponded to the assembly positions between adjacent ring segments. These assembly gaps can cause changes in the scanning distance, resulting in changes in the intensity value of the point cloud data obtained through scanning.
According to
Section 2.1, the intensity values were partitioned into 200 groups by mileage, and the average value of each intensity value was calculated. We drew an average intensity map and highlighted the data with significant differences between the average intensity and the surrounding area in red, as shown in
Figure 11. The mileage of the red difference data in
Figure 11 is compared with
Figure 10, where the red difference data are located at the junction of segments. The red difference data interval was about 1.5 m, which was consistent with the actual width of the shield tunnel segments.
The shield tunnel point cloud was cut out according to the mileage of the red difference data, as shown in
Figure 12a,b. The intensity distribution of a ring point cloud was analyzed, as shown in
Figure 12c.
According to the method given in
Section 2.2, we unfolded the point cloud data of the first ring shield tunnel. We rendered the point cloud based on intensity values, as shown in
Figure 13. From
Figure 13, it can be clearly seen that the strength of the pipe segment gradually changed with the curvature. The change in the intensity value was due to the scanner not being located at the center of the tunnel section, resulting in changes in the distance as the laser reached the surface of the pipe from different angles.
After unfolding the point cloud data into a plane, the intensity value of the point cloud could be added as an additional dimension, as shown in
Figure 14. Then, we drew a single-loop radian–intensity diagram, as shown in
Figure 15. The intensity value of the segment point cloud was only affected by the distance of uniform variation, so the intensity value showed a trend of uniform variation. In
Figure 15, the intensity values are a uniformly changing curve, but there are many abrupt abnormal intensity values in the curve. The intensity values of these mutations are the noise to be removed.
First, according to the method in
Section 2.3, noise was removed by the intensity threshold method. As shown in
Figure 12c and
Figure 13, we set the threshold range of point cloud intensity values to [100, 160] and removed points outside the threshold range.
Next, the intensity values of the single-ring point cloud were partitioned into 200 groups according to radians. Then, according to the method in
Section 2.4, denoising was carried out through the combination of shape and reflection intensity. For each segment point cloud intensity value, we denoised using the optimized IQR and Z scores. The denoising method removed the point cloud of intensity values that varied dramatically in a uniformly changing curve. It could remove points that exhibited drastic changes in a uniformly changing intensity curve. The radian-intensity graph after denoising is shown in
Figure 16a.
Figure 16b shows the radian-intensity plot of the removed noise.
We compared the tunnel point cloud cross-sectional images before and after denoising. An image from before denoising is shown in
Figure 17a, and an image from after denoising is shown in
Figure 17b. From the point cloud cross-sectional view, it can be seen that many noise points such as pipelines and supports in the tunnel before denoising were mostly removed after denoising.
The reflection intensity map of the point cloud after denoising is shown in
Figure 18a, and the reflection intensity map of the noisy point cloud is shown in
Figure 18b. The unfolded point cloud after denoising is shown in
Figure 19a. We expanded the noisy point cloud, as shown in
Figure 19b. It can be clearly seen from
Figure 18 and
Figure 19 that the denoising algorithm achieved the almost complete removal of noise in pipelines, brackets, screw holes, and pipe assembly areas.
The distribution map of the point cloud intensity after denoising is shown in
Figure 20a. The point cloud intensity distribution of noise is shown in
Figure 20b. From
Figure 20b, it can be seen that the point cloud within the intensity threshold contained a large amount of noise, which could not be denoised solely by relying on the intensity threshold.
According to the method in
Section 2.5, we applied geometric shape denoising to the point cloud. First, we performed elliptical fitting on the denoised segment point cloud data to obtain the center position of the segment [
24,
25]. We increased or decreased the design radius of the shield tunnel segment by a certain value at the center point of the segment, forming a circular range of inner and outer circle. We considered the points outside the ring as noise, while retaining the points inside the ring, as shown in
Figure 21.
After denoising, when zooming in on the point cloud data, protrusions and unevenness could still be observed, as shown in
Figure 22. This was because in the denoising processing, a threshold based on the tunnel design radius was used. However, due to the deformation of the segments, some noise that should have been removed may still have remained in the data.
Next, we fit the ellipse to the point cloud data of the tunnel segment and obtained the center point and ellipse equation of the tunnel segment. We increased and decreased the major and minor axes of the ellipse by a certain value to obtain two ellipses, one large and one small, forming a circular range. We considered the points outside the ring as noise, while retaining the points inside the ring, as shown in
Figure 23. After denoising, we zoomed in on the details of the point cloud, as shown in
Figure 24. From the figure, it can be seen that there were no protrusions or other noise points in the segment point cloud data.
4. Discussion
In order to verify the reliability of the denoising method in the paper, this section will verify it in three ways. First, the most intuitive method was to visualize the point cloud and observe the visual effect of the point cloud data before and after denoising, as shown in
Figure 25. For example, the visual effect of denoising could be evaluated by observing the smoothness, detail retention, shape changes, and other aspects of the point cloud.
The second method of verification was to calculate the noise removal rate for evaluation. The noise removal rate quantifies the algorithm’s ability to remove noise by calculating the proportion or percentage of noise successfully removed by the denoising algorithm in the total noise. The noise removed by manual denoising is regarded as all the noise in the point cloud. Every point in the tunnel point cloud was numbered. According to the corresponding number of each point, the noise removed by this algorithm accounted for 99.77% of all noise. It was proved that the denoising effect of this algorithm was excellent. In the correspondence of each point number, it was found that the points removed by the algorithm contained some non-noise points. This paper verifies whether this will affect the denoising effect by cross-section shape preservation.
In the third method, this paper evaluates the degree to which the algorithm maintained the tunnel cross-section shape by cross-section shape preservation. This could be achieved by calculating the shape deviation between the denoised point cloud and the manually denoised point cloud. Ellipse fitting was carried out on the denoised two groups of point clouds, and the major axis and minor axis of the ellipse were obtained. The quantitative results of shape deviation were calculated by the major axis and minor axis of two fitted ellipses. The maximum relative deviation (MRD) is defined as the maximum relative difference between this method and manually denoised point cloud fitted elliptical data. The average deviation (AD) is defined as the average deviation between this method and manually denoised point cloud fitted elliptical data. The mean squared error (MSE) is defined as the mean squared error between this method and manually denoised point cloud fitted elliptical data, as shown in Equation (3). The comparison results of two fitted ellipse data are shown in
Table 1. The results indicated that the shape deviation was small and the cross-section shape preservation was good. At the same time, it was proved that some non-noise points removed by the algorithm had no influence on the denoising effect.
In the equations,
represents the ellipse fitting data obtained from manual denoising,
represents the ellipse fitting data obtained from filtering denoising.
The algorithm took 2 min to denoise the single-ring point cloud (the algorithm is operated on i5 processor by Matlab2022b), while the manual denoising took 15 min. In this paper, the experimental tunnel was 575 m long, including 383 ring segments. Theoretically, if the algorithm processed 16 single-ring point clouds in parallel, it would take 47 min to process the point cloud data of this tunnel. According to the computer performance, the algorithm could achieve higher efficiency by processing more single-ring point clouds in parallel. This paper proves the reliability of the algorithm results from three methods and verifies that the algorithm can be used as the basis for realizing efficient, comprehensive, and automatic tunnel measurement.