1. Introduction
Terrestrial Laser Scanning (TLS) is a powerful geodetic tool ideal for supporting a wide spectrum of applications in many different environments, e.g., [
1,
2]. Registration is a prerequisite in TLS data processing. In this paper, we present a fine registration approach, where after transforming all the scans into the reference coordinate system, point clouds are divided into small blocks, and then, each block is described by a representative point. Then, a rigid body movement and a similarity transformation are used to detect and correct the errors of representative points between two scans. An Iterative Closest Point (ICP) method is applied to estimate the six-parameter transformation, and a Singular Value Decomposition (SVD) algorithm is used to estimate the seven-parameter transformation. The advantages of the proposed method mean that not only the registration errors from scans to the reference coordinate system are processed, but, also, the errors between the scans are considered.
Previous work about registration methods can be found in the following literature. A Least Squares 3D (LS3D) surface matching [
3] utilizes a generalized Gauss–Markov (GM) model and minimizes the sum of the squares of Euclidean distances between two surfaces. The generalized GM model is composed of a GM estimation model and an extension model, where the unknown parameters are introduced as fictitious observations. Later, an integrated model [
4], which considers simultaneously the match surface geometry and intensity information, was introduced. For the LS3D method, the discrete point cloud is represented as a patch-wise surface function [
5]. The difference with other surface matching methods, such as a the Scale Invariant Feature Transform (SIFT) algorithm [
6], is that the LS3D method worked based on the discrete point cloud, while the SIFT algorithm operated based on images. Although a point cloud from a TLS could be generated as an image, the operated datasets between these two methods are different. The presented approach follows closely on the usage of the point clouds from [
7]. The patches defined from discrete point clouds were used as observations in a generalized GM model. Similarly, processing the points firstly, a fine registration of two range images was introduced by first aligning feature points, followed by feature surface elements, which are sets of points corresponding to the image area determined by the feature [
8].
Recently, Grant
et al. [
9] presented a similar approach to [
8]. Instead of any artificial targets, a point-to-plane approach was launched to register with a general least-squares adjustment method. The point-to-plane correspondence was formed from both scans. A point in a scan and the three nearest scanned points in the other scan constructed a transform point and a hypothesized corresponding planar element. Then, the point was transformed to the planar system. These researches mean that firstly processing the point clouds and then using registration methods may be a reasonable way to get more accurate registration results. An attempt was performed earlier in [
10], where a fine registration was estimated with the help of the ICP method. However, there are still some errors existing in the datasets.
For TLS applied in structural monitoring tasks, a uniform framework needs to be established to describe the object surface and to compare the deformation of an object surface between different epochs. Object surfaces could be recognized in the following two methods: those that segment points based on criteria, like the proximity of points, and those that directly estimate surface parameters [
11]. In many cases, the second type of method is suitable and convenient for regular object shapes, like ellipsoids and cylinders, which could be described by several parameters. Thus, in this paper, a quadratic form estimation [
12,
13] is used to describe the shape of the object surface and to establish a uniform framework (see Section 2.2).
Previous work for TLS applied in structural monitoring tasks was often performed in the following two steps (e.g., [
14]). Point clouds are firstly transformed from multiple scans into a reference coordinate system. Then, the object surface is described by a mathematical model for epoch comparisons. However, the errors among different scanner stations were neglected. In order to adjust these errors, the fine registration is introduced. The basic idea of a transformation is to map a point or a feature in one coordinate system to a point or feature in another coordinate system. Many types of functions could be used to specify the mapping relationships: rigid body movement, similarity transformation, affine transformation, and so on. In the field of TLS, the rigid body movement and the similarity transformation are typically applied for data transformations. The rigid body movement estimates six transformation parameters, which includes three rotations and three translations. If there is a scale difference between the scanner and the objects, the similarity transformation can be used (Vosselman and Maas [
15]). For comparisons in this paper, the rigid body transformation and the similarity transformation are used and compared in the process of fine registration.
For the rigid body transformation, an ICP algorithm proposed by Besl and McKay [
16] is commonly used to estimate the six transformation parameters in TLS data processing. A basic version of the ICP method is based on the search of pairs of nearest points in two datasets and, then, to estimate a rigid body movement. The ICP method runs the whole dataset directly in most of the applications (e.g., [
14]), while the LS3D algorithm divides the whole dataset into patches firstly, then iterates to obtain the transformation parameters. At the same time, the ICP and its variant algorithms minimize the Euclidean distances between two point clouds by least-squares and indirectly estimate the rigid body movement, while the LS3D method directly formulates a GeneralizedGM model.
Under a similar assumption that both observations and a datum are contaminated by errors, a SVD algorithm can be used to compute similarity transformation parameters within an Errors-In-Variables (EIV) model [
17], where the EIV model is a specific type of Gauss–Helmert model. In this case, no iteration and no initial values are necessary to solve the transformation parameters. This saves running time when there are thousands of representative points. For long distances between two scans, the seven-parameter transformation may be more suitable, due to estimating the scale factor. In general cases, the SVD algorithm itself can also be used to estimate the six transformation parameters, and the ICP method itself can also be used to estimate the seven transformation parameters. In this paper, the SVD algorithm is used in the fine registration with seven-parameter transformation; while the ICP method is used in the fine registration with six-parameter transformation. Both the two fine registration types can be employed in the proposed procedure, and their suitability will be compared with experimental data.
In the research demonstrated in this paper, the fine registration in TLS applications is processed by first transforming the point clouds into a reference coordinate system, followed by block-to-point estimation; then, the fine registration is performed based on the estimated representative points. The next section of this paper will present the proposed method. Section 3 displays a real application in a scanning task of a dam surface. The conclusion and future work are presented finally.
2. Proposed Method
In this paper, we present a new methodology that extends the work of Eling [
14], Felus and Burtch [
17] and Besl and McKay [
16]. The result is expected to improve the quality of registration and to decrease or eliminate the effect of random and systematic errors between the scans in TLS applications. By block-to-point estimation, points in one block were estimated as a representative point, then six-parameter transformation is applied for fine registration; for comparisons, fine registration with seven-parameter transformation is also proposed. The proposed method processes the errors between different scanner stations, where these errors are neglected in the work of [
14]. Fine registration is employed based on the estimated representative points, not on the original points. Artificial targets are used to estimate the transformation parameters from scans to a local reference system, while the estimated representative points are used in the process of fine registration.
The proposed method is presented in three steps (
Figure 1): Firstly, a combined model is used to transform all the scans to a reference coordinate system with the help of identical points. Secondly, a quadratic form estimation and a segmentation method are introduced to describe the object surface and then divide the object surface into small blocks based on point clouds from TLS. The point cloud in each block is estimated as a representative point. This step is also named block-to-point estimation. Thirdly, based on these representative points, the six-parameter transformation and the seven-parameter transformation are applied to adjust the errors of representative points between multiple scans with the help of the ICP method and the SVD algorithm, respectively.
Figure 1 demonstrates the workflow of TLS applied in structural monitoring tasks. The registration errors between the network and the scans are processed in step 1. Then, the errors between the scanner stations are adjusted in step 3. It may be more reasonable and comprehensive to adjust the errors between the network and the scans by considering both these two groups of errors.
2.2. Block-to-Point Estimation
The registration of point clouds could be achieved, e.g., by applying the ICP method [
16]. However, the explicit point correspondence is hard to confirm, and this causes the lower precision of the registration. In order to compare the possible deformation of objects from different epochs, a uniform framework is established with the help of the estimated parameters of the quadratic form. A segmentation method is then used to divide the object surface into small blocks. In the research presented in this paper, these two steps, named block-to-point estimation, are used to estimate representative points. These points establish explicit correspondences between the two datasets. After that, instead of the original point clouds, the representative points are used to describe the object surface.
A quadratic form could be used by evaluation of determinants and then testing the form parameters to describe the object surface [
12,
13]. The equation for a quadratic form estimation can be written by:
where
xk is the coordinate vector of a single point and
k = 1, 2...,
n,
n denotes the total number of points scanned by TLS. The parameters,
ai, could also be expressed as:
where
i = 1, 2, ...10,
M is the 3 × 3 symmetric coefficient matrix,
m is the coefficient vector and
α is the scalar.
The parameters from
a1 to
a10 could be estimated by a Gauss–Helmert model with several steps of iteration for convergence. To avoid repeating, details are in [
13]. A determinant method, in which an extended form matrix
M* is constructed, is used to estimate the four motion invariant parameters,
δ, Δ,
B and
J [
12].
The extended form matrix is identified:
The functions of the four parameters,
δ, Δ,
B and
J, could be generated as follows:
The shape of the object surface may be determined by looking for the test tree of automatic form recognition [
12]. The advantage of this method is that only several parameters need to be estimated to describe the shape of the object surface.
Based on the estimated parameters of quadratic form estimation, the segmentation is performed in a spherical coordinate system according to the horizontal angles and the vertical angles. The coordinates,
x(
pu,
qv), of the block centers can be derived by:
where
pu and
qv are defined according to the size of the block, with:
where
m and
n are the total numbers of rows and columns of the whole framework, Δ
p and Δ
q are the intervals of neighboring blocks in the vertical direction and the horizontal direction, respectively. The first block starts from the right bottom of the object surface.
According to the point clouds in a block, the least-squares method is used to estimate the representative point in this block. Because both the block size and curvature of the block surface are small, the point clouds in one block usually are distributed on a plane surface. Thus, the normal vectors,
n, and the distance,
d, of a block are defined to evaluate the plane surface:
where
w is the number of points in a block. The covariance matrix,
Σ̂b̂, includes the stochastic information of the estimated parameters and is estimated with the eigenvalues and eigenvectors of the points in a block:
The center point,
xbc, of the block is introduced to estimate the representative point,
x̂r:
A differential function of
n̂ and
d̂ in
Equation (10) is expressed by a matrix,
D:
By means of the variance propagation algorithm, the covariance matrix,
Σx̂rx̂r, of the representative points,
x̂r, is given by:
Thereafter, with the block-to-point estimation, the representative points are prepared to describe the object surface.
2.3. Fine Registration
In the transformation model in step 1, the registration errors between the scanner stations and the reference coordinate system are adjusted, but the errors among the multiple scans are not included. In order to detect and decrease the distance discrepancies among these multiple scans, a rigid body movement and a similarity transformation are used. The ICP method is commonly used to estimate the rigid body movement parameters. For a comparison with the rigid body movement, the similarity transformation is applied with the help of the SVD algorithm.
The red points in
Figure 2 are the estimated representative point correspondence of two blocks from two scan datasets,
X1 and
X2. If there are no systematic or random errors, the distance discrepancy between these two representative points would be zero. In order to adjust this kind of error, fine registration is proposed either with the help of the six-parameter transformation or with the help of the seven-parameter transformation (see step 3 in
Figure 1). The ICP method estimates six parameters (three translations and three rotations), while the SVD algorithm estimates seven parameters (three translations, three rotations and one scale). The performance of each of the algorithms depends on the quality of observations and the measurement environment. If the scale factor is important in a given case, seven-parameter transformation would perform better in the processing of fine registration. Otherwise, if the distances among the scans and the object surface are quite close to each other, the scale factor may be insignificant.
4. Conclusions and Future Work
This paper was focused on block-to-point fine registration in Terrestrial Laser Scanning (TLS), which is used for combined evaluation of point clouds from TLS and of geodetic networks observed using total stations. Of specific interest were those methods that utilize corresponding representative points, as they employ the original point clouds for registration (e.g., [
14]). The block-to-point fine registration approach applies the quadratic form estimation and segmentation methods to guarantee the explicitly corresponding representative points between the scans. Thereafter, the six-parameter transformation was performed for fine registration to correct the distance discrepancies of representative points between the scans by applying the Iterative Closest Point (ICP) method. For comparisons, fine registration with seven-parameter transformation was also applied with the help of the Singular Value Decomposition (SVD) algorithm.
Comparing fine registration with the six-parameter transformation and with the seven-parameter transformation, the seven-parameter transformation worked better, mainly due to the scale factor. Not only the registration errors from the scans to the reference coordinate system were considered, but also the distance discrepancies between the scans were treated. This is meaningful in real applications. In this paper with real TLS data from the Oker dam surface, the mean standard deviation of the distance discrepancies between the scans was decreased by approximately 3.2 mm based on the comparison between seven-parameter fine registration and [
14]. This means that the systematic errors were decreased by approximately 60%.
The proposed block-to-point fine registration can be used in applications for objects with regular shapes. For complicated object surfaces that cannot be described by several parameters, other modeling methods are suggested to describe the object surfaces, such as the Non-Uniform Rational B-Spline (NURBS) model. Future work will try to apply the proposed method to other applications, such as tunnel monitoring tasks.