Application of the Terrestrial Laser Scanning in Slope Deformation Monitoring: Taking a Highway Slope as an Example

: Slope deformation monitoring is the prerequisite for disaster risk assessment and engineering control. Terrestrial laser scanning (TLS) is highly applicable to this field. Coarse registration method of point cloud based on scale-invariant feature transform (SIFT) feature points and ﬁne registration method based on the k-dimensional tree (K-D tree) improved iterative closest point (ICP) algorithm were proposed. The results show that they were superior to other algorithms (such as speeded-up robust features (SURF) feature points, Harris feature points, and Levenberg-Marquardt (LM) improved ICP algorithm) when taking the Stanford Bunny as an example, and had high applicability in coarse and ﬁne registration. In order to integrate the advantages of point measurement and surface measurement, an improved point cloud comparison method was proposed and the optimal model parameters were determined through model tests. A case study was conducted on the left side of the K146 + 150 point at S236 Boshan section, Shandong Province, and research results show that from 14 August 2018 and 9 November 2019, the overall deformation of the slope was small with a maximum value of 0.183 m, and the slope will continue to maintain a stable state without special inducing factors such as earthquake, heavy rainfall and artiﬁcial excavation.


Introduction
The original stress state of natural or artificial slopes will change in the process of geological evolution and engineering construction, resulting in stress redistribution and stress concentration [1]. In order to adapt to the new stress state, the rock and soil of slope will undergo different forms and scales of deformation (such as unloading rebound and creep) [2][3][4]. Revealing the probability and extent of disastrous deformation are the prerequisites for slope stability analysis [5]. Traditional slope deformation monitoring methods include manual survey, global navigation satellite system (GNSS) measurement, total station measurement and close-range photogrammetry [6]. Among them, manual survey is unsuitable for high and steep slopes, and the results depend on human interpretation of the terrain, which is highly subjective [7]; GNSS measurements and total station measurements are single-point measurement methods, which have difficulties in determining the deformation of the areas without monitoring points, moreover, once the monitoring points are destroyed, the continuity of the data will be affected [8,9]. In close-range photogrammetry, each image covers only a small part of the scene and the data processing involves cumbersome processes such as orienteering and inlay, resulting in low work efficiency [10,11]. Slope deformation monitoring based on terrestrial laser scanning (TLS) can provide discrete three-dimensional data of the slope surface, thereby avoiding the locality and one-sidedness of stress-strain analysis based on data from single monitoring point [12]. These advantages ensure that TLS-based slope deformation monitoring has broad application prospects [13,14].
Affected by factors such as the performance of TLS facility and range of slope, the whole information of some slopes needs to be acquired through multi-station scanning, and the point cloud data of each station need to be registered to the same coordinate system [15,16]. Point cloud registration includes feature-based registration and featureless registration [17], of which feature extraction is the core of the feature-based registration algorithm, and representative studies include: Zhu et al. [18] proposed an image registration algorithm named belief propagation (BP)-scale-invariant feature transform (SIFT), where the keypoint matching of SIFT descriptors was formulated as a global optimization problem and provided a suboptimum solution using BP; Chen and Shang [19] presented an improved SIFT algorithm on characteristic statistical distributions and consistency constraints since the feature descriptor section of conventional SIFT does not make full use of the local feature information; Wang et al. [20] proposed an algorithm for automatic extraction of diagonal marks by combining template matching and gradient peaks, which improved the robustness of corner point detection; Lu et al. [21] proposed a super 4-points congruent sets (4PCS) algorithm combining the scale invariant features in light of the low efficiency of the traditional initial registration method; and Ran and Xu [22] proposed a point cloud registration method based on SIFT and geometry feature in order to improve the registration accuracy. The object was extracted based on threshold segmentation, the initial corresponding points were acquired based on the SIFT operator, and furthermore, the curvature was performed to remove the wrong points.
Besl and McKay [23] proposed the iterative closest point (ICP) algorithm, which by far is the most widely used featureless point cloud registration algorithm. It calculates the rotation and translation transformation matrix of the point cloud to be registered by creating a corresponding point pair with the smallest Euclidean distance between the reference point cloud and point cloud to be registered [24]. The ICP algorithm needs to fulfill two basic conditions: (1) There is a common area between the reference point cloud and point cloud to be registered; and (2) The initial relative deviation between the reference point cloud and point cloud to be registered is not large [25,26]. Some scholars have improved the ICP algorithm and representative studies include Wei et al. [27], who proposed an ICP algorithm based on point cloud homography and explained the method of establishing homography point pairs; Li and Song [28], who proposed a vector information registration method based on facet triangles in the stereolithography (STL) format file and introduced a dynamic adjustment factor to speed up the iterative convergence; Cheng et al. [29] proposed an improved ICP algorithm, which used a point-to-point correspondence instead of point-to-surface correspondence to reduce the calculation cost; Wu et al. [30] proposed a novel robust scale ICP algorithm by introducing the maximum correntropy criterion (MCC) as the similarity measure when the point sets had a large number of outliers and noises; Combès and Prima [31] presented an efficient expectation-maximisation (EM)-ICP algorithm for non-linear pairwise registration of large 3D point sets, compared to other methods using the same "EM-ICP" framework, four key modifications leading to an efficient algorithm: (1) truncation of the cost function; (2) symmetrization of the point-to-point correspondences; (3) specification of priors on these correspondences using differential geometry; and (4) efficient encoding of deformations using the reproducing kernel Hilbert space (RKHS) theory and the Fourier analysis.
Slope deformation measurement based on TLS determines the spatial position of the comparison point relative to the reference point by processing the reference point cloud (point cloud of the first phase) and the comparison point cloud (point cloud of the second phase) [32]. It includes four types of methods: point cloud superposition method, point cloud comparison method, center of gravity method, and fitting method [33]. Among them, the point cloud superposition method and point cloud comparison method belong to surface measurement as their monitoring points do not necessarily correspond to one-to-one, and the center of gravity method and fitting method can achieve one-to-one correspondence of monitoring points, equivalent to conventional feature point measurement [34]. Representative studies include Schürch et al. [35], who realized landslide stability monitoring by calculating the landslide volume in different periods based on the point cloud and digital elevation model (DEM); Zhu et al. [36] proposed a method for surface deformation monitoring of high-risk slopes based on TLS that extracted the normal vector of the reference points and calculated the center of gravity of the comparison point cloud, which were used as the basis for comparison of surface deformation; and Hu et al. [37] took the front edge of Jinpingzi landslide II region as the research object, arranged target bodies in "surface monitoring areas" and performed deformation analysis on the target bodies.
Slope deformation monitoring based on TLS can still be improved in the following aspects: (1) Regarding point cloud registration, existing feature extraction algorithms have disadvantages such as low extraction efficiency, require multiple parameters, and have poor real-time performance for large scenes; existing improved ICP algorithms have disadvantages such as high time and space complexity and are prone to matching errors and slow iterative convergence; (2) Regarding slope deformation measurement, point measurement methods have high accuracies, but only the key parts can be monitored and the deformations in areas without monitoring points cannot be revealed, whereas surface measurement methods can describe the deformation of each part of the slope surface, but the accuracies are low. Thus, a slope deformation measurement method integrating the advantages of point measurement and surface measurement does not yet exist.
In light of this, the slope on the left side of the K146 + 150 point at S236 Boshan section, Shandong Province was taken as the research object and the point clouds were collected by TLS, coarse registration of point cloud based on SIFT features and fine registration based on k-dimensional tree (K-D tree) improved ICP algorithm were conducted respectively, and the slope deformations were calculated by the improved point cloud comparison method with the aim to provide a direct basis for disaster risk assessment and engineering control.

Methodology
The goal of point cloud coarse registration is to reduce the rotation and translation dislocation between the reference point cloud and point cloud to be registered, and provide a good initial value for fine registration [17,31]. The goal of fine registration is to obtain the optimal registration parameters and minimize the errors [38]. In this paper, the corresponding point pairs were obtained by extracting the SIFT features, the error point pairs were eliminated based on the random sample consensus (RANSAC) algorithm, the initial matrix was determined by utilizing the singular value decomposition (SVD) algorithm to complete the coarse registration. In order to improve the overall point cloud data registration speed, K-D tree was used to search for the nearest point in the ICP algorithm to complete the fine registration, as shown in Figure 1. necessarily correspond to one-to-one, and the center of gravity method and fitting method can achieve one-to-one correspondence of monitoring points, equivalent to conventional feature point measurement [34]. Representative studies include Schürch et al. [35], who realized landslide stability monitoring by calculating the landslide volume in different periods based on the point cloud and digital elevation model (DEM); Zhu et al. [36] proposed a method for surface deformation monitoring of high-risk slopes based on TLS that extracted the normal vector of the reference points and calculated the center of gravity of the comparison point cloud, which were used as the basis for comparison of surface deformation; and Hu et al. [37] took the front edge of Jinpingzi landslide II region as the research object, arranged target bodies in "surface monitoring areas" and performed deformation analysis on the target bodies. Slope deformation monitoring based on TLS can still be improved in the following aspects: (1) Regarding point cloud registration, existing feature extraction algorithms have disadvantages such as low extraction efficiency, require multiple parameters, and have poor real-time performance for large scenes; existing improved ICP algorithms have disadvantages such as high time and space complexity and are prone to matching errors and slow iterative convergence; (2) Regarding slope deformation measurement, point measurement methods have high accuracies, but only the key parts can be monitored and the deformations in areas without monitoring points cannot be revealed, whereas surface measurement methods can describe the deformation of each part of the slope surface, but the accuracies are low. Thus, a slope deformation measurement method integrating the advantages of point measurement and surface measurement does not yet exist.
In light of this, the slope on the left side of the K146 + 150 point at S236 Boshan section, Shandong Province was taken as the research object and the point clouds were collected by TLS, coarse registration of point cloud based on SIFT features and fine registration based on k-dimensional tree (K-D tree) improved ICP algorithm were conducted respectively, and the slope deformations were calculated by the improved point cloud comparison method with the aim to provide a direct basis for disaster risk assessment and engineering control.

Methodology
The goal of point cloud coarse registration is to reduce the rotation and translation dislocation between the reference point cloud and point cloud to be registered, and provide a good initial value for fine registration [17,31]. The goal of fine registration is to obtain the optimal registration parameters and minimize the errors [38]. In this paper, the corresponding point pairs were obtained by extracting the SIFT features, the error point pairs were eliminated based on the random sample consensus (RANSAC) algorithm, the initial matrix was determined by utilizing the singular value decomposition (SVD) algorithm to complete the coarse registration. In order to improve the overall point cloud data registration speed, K-D tree was used to search for the nearest point in the ICP algorithm to complete the fine registration, as shown in Figure 1.

SIFT Feature Extraction
SIFT feature extraction includes spatial feature point extraction, feature point detection, determination of the feature point direction, and determination of the feature point descriptor [39]. The main processes include: Assume that the coordinate of a point in the point cloud is P(x, y, z) and the spatial scale of the point cloud is L(x, y, z, σ). The relationship between the point coordinate and the spatial scale is as shown in Equation (1).
where σ is the spatial scale factor of the point cloud; G(x, y, z, σ) is the three-dimensional Gaussian kernel function of the point cloud; and the calculation method is as shown in Equation (2).
When constructing a point cloud Gaussian pyramid, the Different-of-Gaussian (DoG) operators at the k scale are defined as shown in Equation (3).
where s ∈ [0, 1]. Using the Lowe operator to linearly fit the Gaussian difference-scale space function and eliminate extreme points with strong edge effects, the remaining extreme points are the feature points. Equation (4) is used to calculate the amplitude m(x, y, z), azimuth θ(x, y, z), and angle ϕ(x, y, z) from the feature point ( Based on the azimuth and angle of the neighborhood of each feature point, define the amplitude as the weight and smooth the histogram, and define the peak direction of the histogram as the feature point direction. Define a four-dimensional vector (σ, m, θ, ϕ) for any feature point, divide its neighborhood into sub-domains of n*n*n, and count the frequency of occurrence of all neighborhood points in the sub-domain in each direction. The frequency of occurrence is the feature point descriptor.

K-D Tree Improved ICP Algorithm
In the K-D tree improved ICP algorithm, the following methods are used to search the nearest point [40][41][42]: (1) Compare the node to be checked with the value of the determined splitting dimension. If the value of the determined splitting dimension is larger, it enters the left sub-tree, and if it is smaller, it enters the right sub-tree. In this manner, loop to the leaf nodes of the binary tree and find the nearest point in the same subspace of the node to be checked along this search path; (2) Perform a "backtrack" operation for each point. If the other subspaces of the nodes along the search path have closer points, skip to the subspace nodes to search the nearest point; and (3) Repeat the above two steps until the search path is empty.

Verification of the Methodology
In order to verify the point cloud data registration method, this study took the Stanford Bunny as the research object to perform the coarse registration based on SIFT feature points, SURF feature points, and Harris feature points, and compared the quantity of extracted feature points, accuracy Appl. Sci. 2020, 10, 2808 5 of 14 rates, and elapsed times [19]. Furthermore, the elapsed times and mean square errors (RMSEs) of the K-D tree improved ICP algorithm and Levenberg-Marquardt (LM) improved ICP algorithm were verified respectively in the fine registration [30].  According to Figure 2, the quantity, accuracy rate, and elapsed time of each kind of feature point were calculated respectively, as shown in Table 1. As can be seen from Table 1, since the Stanford Bunny is a kind of standard model data and contains no noise [43], both the reference point cloud and point cloud to be registered can find strict points with the same name, therefore, the registration accuracies of the three algorithms were relatively high. Comparably, the extraction accuracy rate of the SIFT feature points reached 100% and took the least time, and the overall performance was better than the SURF feature points and Harris feature points. Therefore, it is reasonable and feasible to perform point cloud coarse registration based on SIFT feature points.

Fine Registration
After applying SIFT feature points to perform coarse registration, fine registration utilizing the K-D tree improved ICP algorithm and LM improved ICP algorithm were performed respectively, of which, the LM improved ICP algorithm can solve the problem of slow convergence of the classic ICP algorithm when the Gauss-Newton approximation error is large [44]. The registration results of the two algorithms are shown in Figure 3. According to Figure 2, the quantity, accuracy rate, and elapsed time of each kind of feature point were calculated respectively, as shown in Table 1. As can be seen from Table 1, since the Stanford Bunny is a kind of standard model data and contains no noise [43], both the reference point cloud and point cloud to be registered can find strict points with the same name, therefore, the registration accuracies of the three algorithms were relatively high. Comparably, the extraction accuracy rate of the SIFT feature points reached 100% and took the least time, and the overall performance was better than the SURF feature points and Harris feature points. Therefore, it is reasonable and feasible to perform point cloud coarse registration based on SIFT feature points.

Fine Registration
After applying SIFT feature points to perform coarse registration, fine registration utilizing the K-D tree improved ICP algorithm and LM improved ICP algorithm were performed respectively, of which, the LM improved ICP algorithm can solve the problem of slow convergence of the classic ICP algorithm when the Gauss-Newton approximation error is large [44]. The registration results of the two algorithms are shown in Figure 3. Appl. Sci. 2020, 10, x FOR PEER REVIEW 6 of 14 The advantages and disadvantages of the two algorithms were evaluated by the elapsed time and RMSE; the results are as shown in Table 2. As can be seen from Table 2, the K-D tree improved ICP algorithm was better than the LM improved ICP algorithm. Therefore, it is reasonable and feasible to use the K-D tree improved ICP algorithm to perform fine registration.

Methodology
The integration of point measurement and surface measurement can not only accurately measure the deformation of a single monitoring point, but also obtain the overall deformation of the slope surface [32][33][34]45]. In order to integrate the advantages of point measurement and surface measurement, an improved point cloud comparison method was proposed. This method used the normal vector of the reference point cloud as the axis to construct a prism, and determined the position of the comparison point cloud by calculating the center of gravity in the prism and extracted the three-dimensional deformation information along the normal direction.
The main ideas are as follows [33,46,47]: (1) For any reference point i(x, y, z), construct a covariance matrix based on the neighborhood point set in the sphere with the reference point as the center and R as the radius, and find the eigenvector corresponding to the smallest eigenvalue, that is, the normal vector N of the reference point; (2) Construct a prism with the normal vector N as the axis, d as the side length, and H as the height. The prism contains a reference point set and a comparison point set, of which the reference point i is located in the reference point set; (  The advantages and disadvantages of the two algorithms were evaluated by the elapsed time and RMSE; the results are as shown in Table 2. As can be seen from Table 2, the K-D tree improved ICP algorithm was better than the LM improved ICP algorithm. Therefore, it is reasonable and feasible to use the K-D tree improved ICP algorithm to perform fine registration.

Methodology
The integration of point measurement and surface measurement can not only accurately measure the deformation of a single monitoring point, but also obtain the overall deformation of the slope surface [32][33][34]45]. In order to integrate the advantages of point measurement and surface measurement, an improved point cloud comparison method was proposed. This method used the normal vector of the reference point cloud as the axis to construct a prism, and determined the position of the comparison point cloud by calculating the center of gravity in the prism and extracted the three-dimensional deformation information along the normal direction.
The main ideas are as follows [33,46,47]: (1) For any reference point i(x, y, z), construct a covariance matrix based on the neighborhood point set in the sphere with the reference point as the center and R as the radius, and find the eigenvector corresponding to the smallest eigenvalue, that is, the normal vector N of the reference point; (2) Construct a prism with the normal vector N as the axis, d as the side length, and H as the height. The prism contains a reference point set and a comparison point set, of which the reference point i is located in the reference point set; (3) Project the comparison point set onto the normal vector N, and calculate the barycentric coordinates of the projected points. The barycentric point is the point of the same name corresponding to the reference point i, and the distance of the point pair of the same name is calculated as the deformation amount of the reference point; and (4) Obtain the model parameters with the highest accuracy of slope deformation measurement through model tests. Conduct study on the research object to obtain the amount of slope deformation during the monitoring period based on the optimal parameters. The calculation principle is as shown in Figure 4.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 14 deformation during the monitoring period based on the optimal parameters. The calculation principle is as shown in Figure 4. It is important to note that differences between the reference point cloud and comparison point cloud may be caused by non-slope deformations such as the growth of vegetation, so it is very important to filter the influence of the vegetation on the deformation measurement. An interesting ground filtering method was developed by Li et al. [48], who proposed a novel filter of LiDAR point clouds based on geodesic transformations of mathematical morphology, where the algorithm made it unnecessary to select different window sizes or determine the maximum window size, which could enhance the robustness and automation for unknown environments. In this paper, this algorithm was adopted in ground filtering before the slope deformation measurement.

Model Tests
In order to assess the applicability of the improved point cloud comparison method and analyze the calculation errors, this study conducted model tests to determine the model parameters with the highest measurement accuracy [49,50].
(1) Test model The main body of the test model was a rectangular gray iron plate of 1.2 m × 1.8 m. The inclination of the iron plate was set to 45° to simulate the slope. A black cone magnet was placed on the iron plate to simulate the slope deformation by moving the position of the magnet. The displacement of the magnet was controlled by a Vernier caliper, as shown in Figure 5. It is important to note that differences between the reference point cloud and comparison point cloud may be caused by non-slope deformations such as the growth of vegetation, so it is very important to filter the influence of the vegetation on the deformation measurement. An interesting ground filtering method was developed by Li et al. [48], who proposed a novel filter of LiDAR point clouds based on geodesic transformations of mathematical morphology, where the algorithm made it unnecessary to select different window sizes or determine the maximum window size, which could enhance the robustness and automation for unknown environments. In this paper, this algorithm was adopted in ground filtering before the slope deformation measurement.

Model Tests
In order to assess the applicability of the improved point cloud comparison method and analyze the calculation errors, this study conducted model tests to determine the model parameters with the highest measurement accuracy [49,50].
(1) Test model The main body of the test model was a rectangular gray iron plate of 1.2 m × 1.8 m. The inclination of the iron plate was set to 45 • to simulate the slope. A black cone magnet was placed on the iron plate to simulate the slope deformation by moving the position of the magnet. The displacement of the magnet was controlled by a Vernier caliper, as shown in Figure 5 (2) Contents of tests When conducting the slope deformation measurement, the parameters that affect the accuracy mainly include the interval of scanning points t, radius of the sphere used for normal estimation R, side length of the prism d, and height of the prism H [51,52]. The commonly used values of these parameters are shown in Table 3.

Data Analysis
Taking one parameter as the variable and the other three as fixed parameters to analyze the calculation results [32], there were four cases in total: (1) Taking the interval of scanning points t as the variable Taking the average of the 27 sets of calculation results corresponding to a certain t value, and obtaining the relationship between t and n: t and n have a positive correlation, that is, the smaller t is, the smaller n is, but t < 0.05 m will significantly increase the time consumed for scanning and deformation calculation, reducing the applicability of the model, as shown in Figure 6a. (2) Contents of tests When conducting the slope deformation measurement, the parameters that affect the accuracy mainly include the interval of scanning points t, radius of the sphere used for normal estimation R, side length of the prism d, and height of the prism H [51,52]. The commonly used values of these parameters are shown in Table 3. Table 3. Commonly used values of the parameters.

Model Parameter Common Values
Interval  [33]. Thirty additional calculations were added in the subsequent data analysis stage. A total of 111 calculations were performed to obtain 111 sets of deformations l and errors n, where error n was calculated as shown in Equation (5).

Data Analysis
Taking one parameter as the variable and the other three as fixed parameters to analyze the calculation results [32], there were four cases in total: (1) Taking the interval of scanning points t as the variable Taking the average of the 27 sets of calculation results corresponding to a certain t value, and obtaining the relationship between t and n: t and n have a positive correlation, that is, the smaller t is, the smaller n is, but t < 0.05 m will significantly increase the time consumed for scanning and deformation calculation, reducing the applicability of the model, as shown in Figure 6a (2) Taking the radius of the sphere used for normal estimation R as the variable Similarly, the relationship between R and n is as follows: R and n are inversely correlated, that is, the smaller R is, the larger n is, but R > 0.20 m will greatly increase the time consumed for deformation calculation, as shown in Figure 6b.  Figure 6d.
To conclude, the optimal model parameters with the highest accuracy for slope deformation measurement are as follows: interval of scanning points t = 0.05 m, radius of the sphere used for normal estimation R = 0.20 m, side length of the prism d = 0.20 m, and height of the prism H = 0.50 m.

Research Object
The research object was the left side of the K146 + 150 point at S236 Boshan section, Shandong Province. The gradient was about 80°, the aspect was about 110°, the relative height difference was (2) Taking the radius of the sphere used for normal estimation R as the variable Similarly, the relationship between R and n is as follows: R and n are inversely correlated, that is, the smaller R is, the larger n is, but R > 0.20 m will greatly increase the time consumed for deformation calculation, as shown in Figure 6b.

Research Object
The research object was the left side of the K146 + 150 point at S236 Boshan section, Shandong Province. The gradient was about 80 • , the aspect was about 110 • , the relative height difference was about 30 m, and the width was about 94.5 m. Cambrian limestone was exposed on the surface and there was almost no vegetation coverage. The location of the research object is shown in Figure 7.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 14 about 30 m, and the width was about 94.5 m. Cambrian limestone was exposed on the surface and there was almost no vegetation coverage. The location of the research object is shown in Figure 7. From 4-6 July 2018, Boshan district had been hit by severe rainstorms, rockfall disasters had occurred on the research object, and a number of tension cracks at the top of the slope appeared to penetrate. Once a large-scale instability occurs, it may block S236 and threaten the lives and property of nearby residents and pedestrians.
Two-phase scanning of the research object was performed with a Leica HDS-8800 3D laser scanner in 14 August 2018 and 9 November 2019, respectively. In order to reduce errors and improve accuracy, the pouring position of the concrete pile on the open ground opposite the research object was determined by using total station measurements. Data were collected once a week for a total of four times. The data analysis results verified that the pouring position was fixed. A screw was used to fix the 3d laser scanner on the pile and the coordinate system used was the coordinate system of the 3d laser scanner itself [37]. The scanning process included two steps: coarse scanning and fine scanning [53]. The coarse scanning and fine scanning both used two stations to obtain the whole slope information. Of which, the purpose of coarse scanning was to estimate the scanning time and determine whether the whole slope could be identified [54]. The interval of scanning point was set to 0.20 m for a total of 53,385 points. The purpose of fine scanning was to provide the point cloud data required for deformation measurement [55]. The interval of scanning points was 0.05 m, as determined in Section 3 in this paper, and a total of 851,467 points were scanned.

Calculation Results
Point cloud data registration was carried out using the method proposed in Section 2 of this paper. The total time taken was 34.7 s, and the registration results are shown in Figure 8.  From 4-6 July 2018, Boshan district had been hit by severe rainstorms, rockfall disasters had occurred on the research object, and a number of tension cracks at the top of the slope appeared to penetrate. Once a large-scale instability occurs, it may block S236 and threaten the lives and property of nearby residents and pedestrians.
Two-phase scanning of the research object was performed with a Leica HDS-8800 3D laser scanner in 14 August 2018 and 9 November 2019, respectively. In order to reduce errors and improve accuracy, the pouring position of the concrete pile on the open ground opposite the research object was determined by using total station measurements. Data were collected once a week for a total of four times. The data analysis results verified that the pouring position was fixed. A screw was used to fix the 3d laser scanner on the pile and the coordinate system used was the coordinate system of the 3d laser scanner itself [37]. The scanning process included two steps: coarse scanning and fine scanning [53]. The coarse scanning and fine scanning both used two stations to obtain the whole slope information. Of which, the purpose of coarse scanning was to estimate the scanning time and determine whether the whole slope could be identified [54]. The interval of scanning point was set to 0.20 m for a total of 53,385 points. The purpose of fine scanning was to provide the point cloud data required for deformation measurement [55]. The interval of scanning points was 0.05 m, as determined in Section 3 in this paper, and a total of 851,467 points were scanned.

Calculation Results
Point cloud data registration was carried out using the method proposed in Section 2 of this paper. The total time taken was 34.7 s, and the registration results are shown in Figure 8.
Taking R = 0.20 m, d = 0.20 m, and H = 0.50 m as the model parameters, the slope deformation calculations were carried out using the method proposed in Section 3 of this paper, and the calculation results are shown in Figure 9, where positive values represent upward deformations in the normal direction, and negative values represent downward deformations in the normal direction [45,47]. scanned.

Calculation Results
Point cloud data registration was carried out using the method proposed in Section 2 of this paper. The total time taken was 34.7 s, and the registration results are shown in Figure 8.  Taking R = 0.20 m, d = 0.20 m, and H = 0.50 m as the model parameters, the slope deformation calculations were carried out using the method proposed in Section 3 of this paper, and the calculation results are shown in Figure 9, where positive values represent upward deformations in the normal direction, and negative values represent downward deformations in the normal direction [45,47]. It can be seen from Figure 9 that from 14 August 2018 and 9 November 2019, a total of four rockfall disasters occurred at the wide cracks on the surface of the slope; the overall deformation of the slope was small with a maximum value of 0.183 m and occurred at the middle of the wide tension crack in the upper part of the slope. According to Li et al. [13] and Li et al. [56], the slope will continue to maintain a stable state without special inducing factors such as earthquake, heavy rainfall, and artificial excavation.

Conclusions
The coarse registration method of point cloud based on SIFT feature points and the fine registration method based on the K-D tree improved ICP algorithm were proposed. The results showed that they were superior to other algorithms (such as SURF feature points, Harris feature points, and LM improved ICP algorithm) when taking the Stanford Bunny as an example, and have high applicability in coarse and fine registration. In order to integrate the advantages of point measurement and surface measurement, an improved point cloud comparison method was proposed and the optimal model parameters were determined through model tests. A case study was conducted on the left side of the K146 + 150 point at S236 Boshan section, Shandong Province, where research results showed that from 14 August 2018 and 9 November 2019, the overall deformation of the slope was small with a maximum value of 0.183 m, and the slope will continue to maintain a stable state without special inducing factors such as earthquake, heavy rainfall, and artificial excavation.
Author Contributions: C.Y. and H.L. initiated and developed the concepts related to this research work; C.Y. and Y.L. proposed the point cloud registration method; C.Y. and Z.H. proposed the improved point cloud comparison method; C.Y. wrote the paper draft under Li's guidance. All authors discussed the results and commented on the manuscript. It can be seen from Figure 9 that from 14 August 2018 and 9 November 2019, a total of four rockfall disasters occurred at the wide cracks on the surface of the slope; the overall deformation of the slope was small with a maximum value of 0.183 m and occurred at the middle of the wide tension crack in the upper part of the slope. According to Li et al. [13] and Li et al. [56], the slope will continue to maintain a stable state without special inducing factors such as earthquake, heavy rainfall, and artificial excavation.

Conclusions
The coarse registration method of point cloud based on SIFT feature points and the fine registration method based on the K-D tree improved ICP algorithm were proposed. The results showed that they were superior to other algorithms (such as SURF feature points, Harris feature points, and LM improved ICP algorithm) when taking the Stanford Bunny as an example, and have high applicability in coarse and fine registration. In order to integrate the advantages of point measurement and surface measurement, an improved point cloud comparison method was proposed and the optimal model parameters were determined through model tests. A case study was conducted on the left side of the K146 + 150 point at S236 Boshan section, Shandong Province, where research results showed that from 14 August 2018 and 9 November 2019, the overall deformation of the slope was small with a maximum value of 0.183 m, and the slope will continue to maintain a stable state without special inducing factors such as earthquake, heavy rainfall, and artificial excavation.
Author Contributions: C.Y. and H.L. initiated and developed the concepts related to this research work; C.Y. and Y.L. proposed the point cloud registration method; C.Y. and Z.H. proposed the improved point cloud comparison method; C.Y. wrote the paper draft under Li's guidance. All authors discussed the results and commented on the manuscript. All authors have read and agreed to the published version of the manuscript.