Next Article in Journal
Novel Modification of Integrated Optimization Method for Sensor’s Communication in Wi-Fi Public Networks
Previous Article in Journal
Design and Implementation of a Linear Active Disturbance Rejection Control-Based Position Servo Control System of an Electromotive Valve for Exhaust Gas Recirculation
Previous Article in Special Issue
Vehicle–Bridge Interaction Modelling Using Precise 3D Road Surface Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Template-Guided Hierarchical Multi-View Registration Framework of Unordered Bridge Terrestrial Laser Scanning Data

1
Key Laboratory of New Technology for Construction of Cities in Mountain Area (Ministry of Education), Chongqing University, Chongqing 400045, China
2
School of Civil Engineering, Chongqing University, Chongqing 400045, China
3
Chongqing Academy of Surveying and Mapping, Chongqing 401121, China
4
Technology Innovation Center for Spatio-Temporal Information and Equipment of Intelligent City, Ministry of Natural Resources, Chongqing 401121, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2024, 24(5), 1394; https://doi.org/10.3390/s24051394
Submission received: 20 January 2024 / Revised: 18 February 2024 / Accepted: 20 February 2024 / Published: 21 February 2024
(This article belongs to the Special Issue Sensors in Civil Structural Health Monitoring)

Abstract

:
The registration of bridge point cloud data (PCD) is an important preprocessing step for tasks such as bridge modeling, deformation detection, and bridge health monitoring. However, most existing research on bridge PCD registration only focused on pairwise registration, and payed insufficient attention to multi-view registration. In addition, to recover the overlaps of unordered multiple scans and obtain the merging order, extensive pairwise matching and the creation of a fully connected graph of all scans are often required, resulting in low efficiency. To address these issues, this paper proposes a marker-free template-guided method to align multiple unordered bridge PCD to a global coordinate system. Firstly, by aligning each scan to a given registration template, the overlaps between all the scans are recovered. Secondly, a fully connected graph is created based on the overlaps and scanning locations, and then a graph-partition algorithm is utilized to construct the scan-blocks. Then, the coarse-to-fine registration is performed within each scan-block, and the transformation matrix of coarse registration is obtained using an intelligent optimization algorithm. Finally, global block-to-block registration is performed to align all scans to a unified coordinate reference system. We tested our framework on different bridge point cloud datasets, including a suspension bridge and a continuous rigid frame bridge, to evaluate its accuracy. Experimental results demonstrate that our method has high accuracy.

1. Introduction

According to the Ministry of Transport of the People’s Republic of China, the total number of highway bridges surpassed one million at the end of 2022 in China [1]. These bridges include large-scale structures that span seas, rivers, and lakes, as well as bridges of various sizes. As the operation time of these bridges increases, issues related to their structural health and carrying capacity have become increasingly prominent [2,3]. Regular operation and maintenance management of bridges are essential to ensure traffic safety, extend the lifespan of bridges, and reduce maintenance costs. There are several methods available to obtain structural dynamic parameters that can reflect the health condition of bridges. One commonly used approach is the eigen perturbation method [4], which enables the extraction of modal parameters. In addition to analytical methods, manual inspections remain essential for certain tasks such as crack inspections, settlement detection, bridge alignment checks. However, traditional manual maintenance methods are expensive, time-consuming, and subjective and relies on inspectors to obtain all information accurately [5].
Recently, the three-dimensional (3D) laser scanning technique has opened up new possibilities for bridge health monitoring and maintenance [6]. The terrestrial laser scanner (TLS) can acquire 3D point cloud data (PCD) of bridges in a short period of time. PCD provides a comprehensive and accurate representation of bridges, allowing for a detailed assessment of its current condition. By analyzing the acquired PCD, valuable insights can be gained to guide maintenance efforts and ensure the longevity and safety of the bridge structure. Furthermore, PCD can be leveraged to create building information modelling (BIM) models. A BIM model provides a digital representation of a bridge and can be valuable for bridge maintenance and digital twins [5,7]. To obtain PCD of the entire bridge, it is necessary to perform multiple scans from different locations because of the limited view of TLS. However, each scan has its own coordinate system, making it essential to align all the scans into a unified coordinate system in order to generate an accurate 3D model of the bridge. Because the accuracy of the 3D model is heavily influenced by the registration accuracy, PCD registration therefore becomes a crucial preprocessing step for TLS-based bridge health monitoring and maintenance.
Although some researchers have studied registration methods of bridge PCD, they mainly focused on pairwise registration [8,9]. However, to align multiple scans of a bridge to a unified coordinate system, the utilization of multi-view registration techniques is essential. In engineering practice, artificial markers such as target spheres [10] and target papers [11] have been commonly used to assist registration and improve registration accuracy when dealing with multiple scans. Although the registration results using these artificial markers are highly reliable, placing artificial markers is time-consuming and costly [12], especially for large infrastructure like bridges. Extensive research has been conducted on marker-free multi-view registration techniques. However, when it comes to applying these techniques to the registration of bridge PCD, there are two notable challenges that need to be addressed. The first challenge involves recovering the overlaps between unordered multiple scans. This typically requires extensive pairwise matching and the creation of a fully connected graph that encompasses all the scans [13]. However, when deal with TLS scans of large bridges, the process of extensive pairwise matching can become quite time-consuming. The second challenge is related to determining the merging order of the scans. Typically, the creation of a fully connected graph that encompasses all the scans is necessary to establish the merging order [14]. However, when dealing with a large number of scans, creating such a fully connected graph can be inefficient. To tackle this challenge, Wu et al. [12] proposed a method to simplify this heavy registration task by subdividing all the scans into several scan-blocks. In this way, the number of pairwise registration procedures is significantly reduced. However, they constructed scan-blocks according to the scanning order, which is sometimes impractical in bridge scenes. Because bridges are essential components of urban transportation infrastructure, minimizing the disruption to urban traffic during the scanning process is crucial. As a result, the PCD of bridge scans cannot always be registered in a sequential order. Manual specification of the registration sequence becomes necessary in such cases. However, for large bridges, where dozens of scans may be acquired, the manual specification process can become quite tedious and inefficient.
To address these issues, this paper proposes a template-guided hierarchical multi-view registration framework that can register unordered bridge terrestrial laser scanning data without any artificial targets. Firstly, the overlaps of multiple scans are recovered using a template-guided initial pose estimation method and extensive pairwise matching is avoided. Next, all the scans are partitioned into different scan-blocks based on their locations and overlaps. Subsequently, the pairwise coarse registration is conducted in each scanning block, and the transformations are obtained using an intelligent optimization algorithm. After that, the fine registration is then performed to further refine the poses within each block. Finally, the above steps are repeated between blocks until all scans are merged into a unified common coordinate system. The main contributions of the proposed method are as follows:
(1)
A marker-free multi-view registration framework is proposed to hierarchically align unordered bridge terrestrial laser scanning data.
(2)
A template-based initial pose estimation method is proposed to recover the overlaps of unordered PCD, which avoids extensive pairwise matching and improves the efficiency.
(3)
To group scans with high overlaps into the same block, a graph partition algorithm based on the overlaps and scanning locations is utilized to construct scan-blocks.

2. Research Background

2.1. Registration of Bridge PCD

Several researchers have conducted studies on the registration techniques for bridge PCD. For example, Data et al. [8] proposed a method for pairwise registration of bridge PCD. They relied on straight-line edges as linear features and used the random sample consensus (RANSAC) algorithm with a hash table of line pairs to match line pairs. The experiments in three bridge datasets demonstrated its accuracy and efficiency. Zhao et al. [15] proposed an efficient local descriptor for pairwise registration and validated its accuracy using a dataset of a large-scale high-pier concrete bridge. Based on a long-span suspension bridge dataset, Zhang et al. [9] compared and evaluated three iterative closest point (ICP) registration methods from various aspects including convergence rate, execution time, and accuracy. They also provided suggestions for efficiency optimization and accuracy improvement to enhance long-span bridge deformation analysis. Deng et al. [16] introduced a novel dual-purpose target for total station and laser scanner applications. Their experimental results on a long-span arch bridge demonstrated that the using this target can improve the accuracy of PCD registration. While these studies have made significant contributions to the field of bridge PCD registration, it is important to note that they primarily focused on pairwise registration and did not explore multi-view registration methods. However, multi-view registration is more crucial in determining the speed and accuracy of bridge PCD registration.

2.2. Multi-View Registration

For unordered multiple scans, multi-view registration first obtains the scanning order of all the scans, and then merges them into a single PCD. The methods for obtaining the scanning order can be roughly classified to marker-based methods and marker-free methods. Marker-based methods involve placing unique markers in the scanning scene and identifying these markers from the scanned PCD to assist with registration. For example, Singh et al. [17] developed 3D unique identifiers along with a 3D registration workflow for mapping and monitoring applications in underground mines. Ge et al. [18] presented an automated marker-based approach for the registration of unordered scans in complex forest scenes, which can automatically detect the artificial markers and build a geometric network to judge their connectivity. However, since installing markers is time-consuming and expensive, these methods are only suitable for complex or single repetitive scenarios. For marker-free methods, Dong et al. [13] utilized multi-level descriptors, including a local descriptor called binary shape context (BSC) and a global descriptor named the vector of locally aggregated descriptor (VLAD), to describe each scan. They evaluated the similarity between every pair of scans based on the distance between their features. Then, they iteratively merged two scans until all the scans were combined into a single coordinate system. However, the performance of VLAD is directly influenced by the distinctiveness of the descriptor. To enhance the robustness of descriptors in urban scenarios dominated by planar structures, Ge et al. [19] proposed a new descriptor by combining panoramic images and the planar features of PCD. However, due to the diverse types of components in bridge scenes, it is challenging to segment them using planar structures, making this method difficult to apply in such scenarios. In addition, the extensive pairwise matching is time-consuming when dealing with a large number of scans.
The methods for merging scans can be roughly divided into sequential registration methods and joint registration methods. Sequential registration methods utilize pairwise registration to iteratively merge two scans until all the scans are merged into a common coordinate system [20]. To determine the merging order, Weber et al. [14] created a fully connected graph and extracted the minimum spanning tree (MST) to define the merging order from root to leaf. The weights of edges in the fully connected graph were defined based on the number of true correspondences between two scans. However, when applied to large-scale scans, searching for corresponding point pairs can be time-consuming. To improve the efficiency of the MST-based methods, some researchers have adopted shape growing algorithms to merge scans, which have higher computational efficiency [19,21]. Other researchers have employed hierarchical merging-based methods to enhance the robustness and efficiency of scans with limited overlaps. For example, Wu et al. [12] proposed a framework for aligning ordered multi-view PCDs. Their approach involves partitioning all the scans into blocks based on the scanning order and then merging the scans within each block to achieve alignment. However, this method is not suitable for unordered scans. Joint registration-based methods consider the multi-view registration problem as a graph optimization problem by using all pairwise registration results to minimize global residuals. Theiler et al. [22] constructed an energy function by considering all loop consistency constraints in the graph and then utilized the lazy flipper algorithm to minimize the energy function, thereby reducing the overall registration error. Tang et al. [23] proposed a method to minimize the sum of squared distances between point pairs generated by the ICP algorithm. Their approach includes a constraint that ensures the product of all transformations within each loop results in an identity matrix. This constraint helps to optimize the alignment process and improve the accuracy of the registration results. However, these methods only redistribute the registration errors on the constructed graph without updating correspondences, and they cannot reduce the total registration error [20].

3. Methodology

The proposed method aims to addresses the challenge of automatically aligning unordered scans of bridges without any makers. The overview of the registration framework is plotted in Figure 1. For unordered scans, a template-guided approach is employed in Section 3.1 to recover the overlaps. Subsequently, all the scans are partitioned into different scan-blocks based on their scanning locations and overlaps in Section 3.2. Within each block, pairwise coarse registration is conducted using an intelligent optimization algorithm, as described in Section 3.3. Finally, the fine registration is performed to further refine the poses within each block in Section 3.4. These steps are iteratively executed between blocks until all scans are merged into a unified common coordinate system.

3.1. Template-Guided Initial Pose Estimation

Typically, bridges are long and symmetrical structures, and their geometric features on the side view are distinct and can be leveraged to approximate the location of each scan. By comparing the geometric features of the side view of an individual scan of the entire bridge, the approximate position of that scan along the bridge’s length can be determined. Therefore, using the as-designed side view as the template, an initial pose estimation method is proposed.

3.1.1. Acquisition of Side View Geometric Features

The side view geometric features of each scan can be obtained by projecting the PCD into a binary image along the direction perpendicular to the traffic direction of a bridge. First, the traffic direction of a bridge can be determined using the principal component analysis (PCA) algorithm. Considering the non-uniformity of scans, only the points that form the minimum convex hull of two-dimensional (2D) PCD on the xOy plane are utilized in the PCA algorithm. Second, the PCD is rotated so that its traffic direction aligns with the x axis, as shown in Figure 2a. Then, the 2D PCD can be obtained by projecting the PCD into the xOz plane. Next, the 2D PCD is voxelized using the grid sizes δx and δz (their calculation methods will be introduced later), and then it is converted to a binary image with each grid corresponding to a pixel. If the number of points within a grid is more than one, the grey value of the corresponding pixel is set to 1; Otherwise, it is set to 0. The binary image of the PCD in Figure 2a is shown in Figure 2b.

3.1.2. Unified Scale of the Template and Binary Images

We employ an as-designed side view as the registration template. To ensure that the scale of the template is same as the binary images converted by PCD, the grid sizes δx and δz are determined using the following two formulas.
δ x = L N x
δ z = H N y
where L and H represent the length and height of the bridge, respectively; Nx and Ny are the width and height of the minimum bounding box of the non-blank region in the template.

3.1.3. Image Matching

Image matching algorithms can be utilized to obtain the approximate positions of scans relative to the bridge. Based on the utilized primitives, image-matching algorithms can be classified as area-based matching (ABM) or feature-based matching (FBM) algorithms [24]. ABM is based on the idea that grey values of pixels of conjugate points have similar radiometric characteristics [25], while FBM is based on feature extraction, feature description, and correspondence feature matching. We opted for the ABM for image matching, which has high efficiency and reliability [26]. Given two images, ABM considers one as the reference image, the other one as the matching image. The matching image slides within a search window on the reference image, and a similarity measure is calculated at each position. The location of the matching image is assumed to be the position of the best agreement [27]. In our study, the registration template is set to the reference image and the search window is the entire image. Binary images are the matching images. An example of image matching is shown in Figure 3a. Considering that the image matching step can only determine the positions of PCDs but cannot adjust their orientations, the matching score between the horizontally mirrored image of the binary image and the registration template is also calculated. The case with the higher matching score is selected as the final result. Based on the matching relationship between each scan and the registration template, the relative positional relationships on the xOz plane between different scans can be obtained, as shown in Figure 3b. The alignment in the y-axis direction is determined based on the centres of the scans.

3.1.4. Correction of False Matching

In some cases, same geometric patterns may appear multiple times in a bridge, which may result in false matching. As shown in Figure 4, the geometric elements within the red boxes appear four times. When the coverage area of a scan is small, these repeated patterns may cause false matching.
To address this issue, a method for identifying and correcting false matching is proposed. A good scanning plan should provide a uniform coverage of a bridge while avoiding over-coverage in any particular area. When a scan is matched to a false location, the scan is considered “redundant” for the false location because it results in an over-coverage for the false location. Therefore, we can quantify the redundancy of a scan by comparing it with its adjacent scans. For the ith scan, its projected image is denoted as Ii, and its areas intersecting with all the other scans are calculated based on the locations obtained in the template matching step. The n scans with the largest intersection areas are selected, their projected images denoted as I1, I2, …, In are utilized to calculate the redundancy score Si for the ith scan using Formula (3).
S i = I 1 I 2 I n I i a = 1 A b = 1 B I i a , b
where A and B are the height and width of Ii, respectively; ( ) denotes the Hadamard product. The scores for all the scans are then statistically analyzed to calculate the coefficient of variation. A higher coefficient of variation indicates a greater discreteness of scans, meaning false matching may be present. The top 20% scans with the highest scores are selected as potentially mismatched PCDs, awaiting further validation.
Considering that there may be differences between the as-designed side view and the actual state of the bridge, the similarity between a scan with its adjacent scans in the right location will be higher than the similarity between the scan with its adjacent scans in the false location. Therefore, further validation of whether a scan is mismatched can be determined based on its positional relationship with adjacent scans. For the ith scan which is labelled as potentially mismatched, three locations with the top three highest matching scores in the template matching step are considered as candidate positions. All the candidate positions are the pixels with the highest matching scores in a certain neighbourhood area. For each candidate position, the intersection areas with other scans that have not been labelled as mismatched are calculated. Three scans with the largest intersection areas are selected, and their projected images are added to the registration template with a proportion of 5% to modify the template, as shown in Figure 5. The matching score between the modified registration template and the ith scan is computed at each candidate position. The candidate position with the highest matching score is selected as the final position for the ith scan.

3.2. Overlap-Based Scan-Block Construction

The main focus of our work lies in the efficiency of processing large bridges, such as suspension bridges, which may require dozens of TLS scans. The fully connected graph that arises from all these scans can be quite large, making the registration task computationally intensive. To address this challenge, a hierarchical registration strategy can be employed. This strategy can recursively subdivide and fuse the heavy registration task into smaller, more manageable subsets. By breaking down the registration process into hierarchical levels, the number of pairwise registration procedures can be significantly reduced [12]. In more detail, the hierarchical registration strategy involves partitioning all the scans into different scan-blocks and then locally aligns the scans in each scan-block, and finally performs the global block-to-block registration. The local registration accuracy of each scan-block plays a vital role in determining the overall registration accuracy of all the scans. To improve the registration accuracy, it is crucial to ensure significant overlaps within each block. Therefore, the block partitioning step plays a critical role in improving both registration precision and efficiency.
Partitioning scan-blocks must take into account both the scanning locations and overlaps between scans, and the overlaps of two scans can be qualified by the overlaps of their bounding boxes. The normalized cut (Ncut) algorithm [28] is utilized in our work for scan-block construction. First, a fully connected graph is created with each scan treated as a node. The edge weight wij between the nodes vi and vj is defined as
w i j = I o U Δ x + Δ z
where Δx represents the ratio between the difference in x-coordinate values of the scanner positions and the width of the bridge, and Δz the ratio between the interpolated z-coordinate values and the height of the bridge. IoU is defined as the intersection ratio of the bounding boxes of two PCDs, i.e.,
I o U = B i B j B i B j
Bj and Bi in Equation (5) represent the bounding boxes of two scans.
The objective of the Ncut algorithm is to minimize the cut between different blocks while maximizing the sum of the edge weights within each block. The objective function can be formulated as
N c u t A 1 , A 2 , , A k = i = 1 k c u t A i , A ¯ i v o l A i
where Ai represents the ith scan-block, A ¯ i the complement of Ai, k the total number of blocks, and cut(Ai, Aj) the cut between Ai and jth scan-block Aj, i.e.,
c u t A i , A j = i A i , j A j w i j
vol(Ai) represents the sum of degrees of each node within Ai, i.e.,
v o l A i = i A d i
Based on the normalized Laplacian matrix, the objective function can be simplified, and the Rayleigh–Ritz theorem can be employed for the solution. Finally, the partitioning results can be obtained using the k-means clustering algorithm [29]. Using the TLS scans of a suspension bridge as an example, the scan-block construction result is shown in Figure 6, where the numbers within the purple circles indicate the scanner positioned above the bridge deck, while the numbers within the black circles represent the scanner positioned below the bridge deck. It can be observed that our block construction method is capable of grouping scans with closer proximity and higher overlap into the same block.

3.3. Pairwise Coarse Registration by Optimization Algorithms

The pairwise coarse registration is commonly formulated as the maximum consensus set (MCS) problem [12]. Given two scans, P and Q, the point pi in P is paired with its nearest neighbour qj in Q, forming a point pair (pi, qj). The set of point pairs is represented as H = {(pi, qj)}1k, where k is the number of point pairs. If the distance between pi and qj is less than a threshold δ, the point pair is considered as a true correspondence. The MCS problem [30] aims to find the transformation matrix corresponding to the maximum number of true correspondences, i.e.,
max R , t , I H I ,   s u b j e c t   t o   p i R q j + t < δ ,   p i , q j I
The feature-based coarse registration method is the most commonly used coarse registration method [20]. It first extracts key points from two PCDs, and computes features of key points. Then a subset of key points is randomly selected from one PCD, and their corresponding points are searched in the other PCD using the feature similarity. After that, a transformation matrix can be obtained based on the matched point pairs, along with the number of true correspondences, which is called a consensus set. This process is usually iterated multiple times, and the transformation matrix associated with the highest number of true correspondences, which is called the maximum consensus set, will be output as the final transformations. However, existing feature extraction methods are susceptible to variations in point density and noise, and it makes the above methods less robust [20]. Considering that all the scans already have rough relative position relationships after the template matching step, performing pairwise coarse registration on the basis of this can significantly reduce the search space. We employ the particle swarm optimization (PSO) algorithm to search for the maximum consensus set.
The PSO algorithm is an iterative intelligent optimization algorithm that relies on collaboration and information sharing among particles to search for the optimal solution. During the searching process, each particle records its current position as well as its historical best solution, and the population also records its historical best solution. Based on the historical best solutions, the positions and velocities of the particles are updated, enabling the particle swarm to iteratively evolve and converge towards the global optimum. In this study, the particles denote the combinations of R and T, and each particle is described by speed vij and position xij, which are updated by Equations (10) and (11).
v i j ( t + 1 ) = w v i j ( t ) + c 1 × r 1 × ( p i j x i j ( t ) ) + c 2 × r 2 × ( p g j x i j ( t ) )
x i j ( t + 1 ) = x i j ( t ) + v i j ( t + 1 )
where i is the particle number and j is a variable dimension; vij(t) and vij(t + 1) are the particle speeds at times t and t + 1, respectively; xij(t) and xij(t + 1) are the particle positions at times t and t + 1, respectively; pij is the historical optimal solution of the current particle; and pgj is the historical optimal solution of the swarm; c1 and c2 are acceleration constants and both set to 2; w is the inertia weight and is set to 0.8; r1 and r2 are random numbers in the closed interval [0,1].
For the initial pose shown in Figure 7a, the PSO algorithm can obtain an accurate transformation, as shown in Figure 7b. Although the PSO algorithm can avoid becoming trapped in local optima, there may still be cases where two scans within the same scan-block cannot obtain the correct transformation matrix due to low overlap or slightly larger scanning distances. To identify false matches, Wu et al. [12] proposed a method based on the loop closure constraint and proved to be effective. This method is utilized in this study to reject false scan-to-scan matches obtained by the PSO algorithm.

3.4. Fine Registration and Pose Optimization

To obtain the optimal merging order of scans within each scan-block, the minimum spanning tree (MST) is utilized to extract a cycle-free and well-pairwise-registered graph [30]. The MST relies on edge weights to define the shortest path. In this study, the weight of each edge is defined as the number of true correspondences after the pairwise coarse registration, and the weights of the edges with false matches are set to 0. Along the edges of the MST, all the other nodes in a scan-block can be merged into the root node using the coarse registration matrix calculated by the PSO algorithm and the fine registration matrix computed by the ICP algorithm [31]. However, the errors will accumulate along the edges from the root to the leaf node since the MST is a cycle-free graph [13]. To address the issue of error accumulation, the Lu–Milios algorithm [32] is utilized to further optimize the pose.
Treating each scan-block as a new scan, the registration between scan-blocks can be accomplished using the same methods described in Section 1, Section 2, Section 3 and Section 4. This process is repeated until all scans are merged into a single scan. To avoid excessive point growth, PCD down-sampling is performed after the scan merging step.

4. Experiments and Analysis

4.1. Datasets Description and Evaluation Criteria

We evaluated the performance of the proposed method using two bridge point cloud datasets, including a suspension bridge and a continuous rigid frame bridge. The suspension bridge named Cuntan Yangtze River Bridge is located in Chongqing with a total length of approximately 1.6 km and a width of about 38 m (Figure 8a). It was scanned using Leica P40 with a ranging error of 1.2 mm + 10 ppm and an angular accuracy of 8″ [33]. The continuous rigid frame bridge named Huanghuayuan Jialing River Bridge is located in Chongqing with a total length of approximately 1.2 km and a width of about 31 m. It was scanned using Faro S350 with a ranging error of 1 mm between 10 m to 25 m and an angular accuracy of 19″ [34]. More details about the two datasets are listed in Table 1.
For the Huanghuayuan Jialing River Bridge, only the area below the bridge deck was scanned using TLS, while the bridge deck area was scanned using a mobile scanning system. Prior to inputting the data into the algorithm, background points were roughly removed. In addition, the ground truth result was obtained through manual registration.
The performance of the method is evaluated by the axis-angle rotation error er and translation error et of all transformations among the scans [35], which are as follows:
e r = arccos t r R g R e T 1 2
e t = t e t g
where Re and te represent the estimated rotation matrix and translation vector, respectively, and Rg and tg are the those of the ground truth. In addition, the successful registration rate (SRR) is also utilized and defined by
SRR = N s N 1
where N represents the total number of scans, and Ns the number of successfully aligned scans. A scan is considered successfully aligned when its rotation error and translation error are both less than the specified thresholds σr and σt, respectively.

4.2. Results of Template-Guided Initial Pose Estimation

Based on the bridge datasets, we evaluated the accuracy and efficiency of the template-guided initial pose estimation method. The proposed method is implemented in Python through an Intel Core i7-7700K CPU (Intel, Santa Clara, CA, USA). The initial pose estimation results of the Cuntan Yangtze River Bridge and Huanghuayuan Jialing River Bridge are shown in Figure 9a and Figure 9b, respectively. It can be seen that this step can roughly align the scans and obtain the relative scanning positions.
The average rotation error, average translation error and running time for two bridges are listed in Table 2. For the Cuntan Yangtze River Bridge, the average rotation error is 19.1 mdeg, and the average translation error in the three coordinate axis directions are 0.66 m, 3.62 m, and 1.88 m, respectively. For the Huanghuayuan Jialing River Bridge, the average rotation error is 19.7 mdeg, and the average translation error in the three coordinate axis directions are 0.64 m, 5.16 m, and 2.35 m, respectively. Due to the alignment in the y-axis direction being solely based on the centres of the PCD, translation errors are more pronounced in the y-coordinate.
The running times of the two bridges are 6.08 min and 4.43 min, respectively. Compared with extensive pairwise matching, template-guided initial pose estimation has a higher efficiency, primarily due to two key reasons. Firstly, in our method, only N template matching is required to create a graph, where N represents the number of scans. In contrast, extensive pairwise matching necessitates C N 2 = N N 1 / 2 matches. This reduction in the number of matches significantly improves the efficiency of our approach. Secondly, in template-guided initial pose estimation, the problem scale in each individual match is determined by the number of points in the scan, denoted as n. The fundamental operation is mapping each point to a 2D grid. As a result, the time complexity for processing a single scan in the initial pose estimation step is O(n), and the time complexity for processing all scans is O(Nn). The quadratic time complexity ensures that our method has high efficiency and is adaptable to larger datasets. In summary, the data showed that our method can achieve a relatively high alignment accuracy, provide good initial poses for subsequent steps, and it helps to improve efficiency by avoiding extensive pairwise matching.

4.3. End-to-End Performance Evaluation

During the experiments, we set the average number of scans in a scan-block to 5, the down-sampling voxel size to 0.1 m, and σr and σt to 100 mdeg and 100 mm, respectively. The experimental setup utilized a hybrid programming approach, with the template-guided initial pose estimation step implemented in Python and the remaining parts in C++. Our method was tested on a laptop with 32 GB RAM and an Intel Core i7-7700K CPU.
The registration results for the two bridges are plotted in Figure 10. The points are coloured according to the difference in height. It can be seen that the mismatching in Figure 9 has been removed. This demonstrates that our method can deal with PCD of large bridges, and obtain a good performance in accuracy.
The accuracy and efficiency of the proposed method was evaluated, and the rotation and translation errors with the root mean square error (RMSE), SRR, and running time are listed in Table 3. The average rotation errors of our method are 0.96, 0.74 (mdeg) and average translation errors are 28.04, 43.25 (mm), with the SSR of 100% and 100%. The results show that our method achieves relatively small rotation errors, but slightly larger translation errors. Compared with the Cuntan Yangtze River Bridge, the Huanghuayuan Jialing River Bridge has a relatively simple geometric shape and lacks complex features and structures, resulting a larger translation error. In conclusion, the registration results listed above prove that our method performs well in registering the TLS scans of varying bridges, and the accuracy can satisfy the requirements of component extraction and 3D reconstruction.

4.4. Discussion

Although our method has demonstrated good performance, it still has two limitations. Firstly, prior to registration, the manual removal of background points is required. This step can become cumbersome when dealing with a large number of scans. Therefore, future research will focus on developing automatic background-point removal techniques. Secondly, our method is currently only applicable to bridges, which limits its application. Future research will aim to expand this method to more diverse scenarios and environments.

5. Conclusions

In this paper, we present a template-guided hierarchical multi-view registration framework for aligning unordered bridge terrestrial laser scanning data without any markers. The proposed framework incorporates two distinct features: template-guided initial pose estimation and overlap-based scan-block generation. The former enables the rapid recovery of the scanning sequence for multiple unordered bridge scans, while the latter leverages the overlap between scans and a graph segmentation algorithm to partition the scans into different scan-blocks. Experimental results on two different bridge datasets demonstrate the accuracy and efficiency of the proposed method.

Author Contributions

Conceptualization, N.C. and J.L.; Methodology, G.X., N.C., Y.Z. and H.X.; Software, N.C. and H.C.; Validation, G.X.; Formal analysis, N.C., Y.Z. and H.C.; Investigation, G.X. and N.C.; Resources, J.L.; Data curation, G.X., N.C., Y.Z. and C.H.; Writing—original draft, N.C. and H.X.; Writing—review & editing, H.X.; Visualization, C.H.; Supervision, J.L.; Project administration, J.L. and H.X.; Funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 52130801).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is unavailable due to privacy.

Conflicts of Interest

We declare that we do not have any commercial or associative interests that represents a conflict of interest in connection with the work submitted.

References

  1. Ministry of Transport of the People’s Republic of China. Statistical Bulletin on the Development of the Transport Industry in 2022; Ministry of Transport of the People’s Republic of China: Beijing, China, 2023. (In Chinese)
  2. Yang, Y.B.; Chen, J.; Xu, H. Normalized formula for removing damping effect in recovering bridge mode shapes using a moving and a stationary vehicle. J. Sound Vib. 2023, 573, 118219. [Google Scholar] [CrossRef]
  3. Xu, H.; Liu, Y.H.; Chen, J.; Yang, D.S.; Yang, Y.B. Novel formula for determining bridge damping ratio from two wheels of a scanning vehicle by wavelet transform. Mech. Syst. Signal Pract. 2024, 208, 111026. [Google Scholar] [CrossRef]
  4. Bhowmik, B.; Tripura, T.; Hazra, B.; Pakrashi, V. Real time structural modal identification using recursive canonical correlation analysis and application towards online structural damage detection. J. Sound Vib. 2020, 468, 115101. [Google Scholar] [CrossRef]
  5. Kaewunruen, S.; Sresakoolchai, J.; Ma, W.; Phil-Ebosie, O. Digital twin aided vulnerability assessment and risk-based maintenance planning of bridge infrastructures exposed to extreme conditions. Sustainability 2021, 13, 2051. [Google Scholar] [CrossRef]
  6. Rashidi, M.; Mohammadi, M.; Sadeghlou Kivi, S.; Abdolvand, M.M.; Truong-Hong, L.; Samali, B. A decade of modern bridge monitoring using terrestrial laser scanning: Review and future directions. Remote Sens. 2020, 12, 3796. [Google Scholar] [CrossRef]
  7. Kaewunruen, S.; AbdelHadi, M.; Kongpuang, M.; Pansuk, W.; Remennikov, A.M. Digital Twins for Managing Railway Bridge Maintenance, Resilience, and Climate Change Adaptation. Sensors 2023, 23, 252. [Google Scholar] [CrossRef]
  8. Date, H.; Yokoyama, T.; Kanai, S.; Hada, Y.; Nakao, M.; Sugawara, T. Efficient registration of laser-scanned point clouds of bridges using linear features. Int. J. Autom. Technol. 2018, 12, 328–338. [Google Scholar] [CrossRef]
  9. Zhang, H.; Zhu, Y.; Xiong, W.; Cai, C.S. Point Cloud Registration Methods for Long-Span Bridge Spatial Deformation Monitoring Using Terrestrial Laser Scanning. Struct. Control Health Monit. 2023, 2023, 2629418. [Google Scholar] [CrossRef]
  10. Cheng, G.; Liu, J.; Cui, N.; Hu, H.; Xu, C.; Tang, J. Virtual trial assembly of large steel members with bolted connections based on point cloud data. Automat. Constr. 2023, 151, 104866. [Google Scholar] [CrossRef]
  11. Liu, J.; Cui, N.; Zhou, X.; Li, D.; Cheng, G.; Zeng, Y.; Cao, Y. Intelligent inspection method for dimensional quality of houses based on 3D laser scanning. J. Archit. Civ. Eng. 2022, 71–80. (In Chinese) [Google Scholar]
  12. Wu, H.; Yan, L.; Xie, H.; Wei, P.; Dai, J. A hierarchical multiview registration framework of tls point clouds based on loop constraint. ISPRS J. Photogramm. 2023, 195, 65–76. [Google Scholar] [CrossRef]
  13. Dong, Z.; Yang, B.; Liang, F.; Huang, R.; Scherer, S. Hierarchical registration of unordered TLS point clouds based on binary shape context descriptor. ISPRS J. Photogramm. 2018, 144, 61–79. [Google Scholar] [CrossRef]
  14. Weber, T.; Hänsch, R.; Hellwich, O. Automatic registration of unordered point clouds acquired by Kinect sensors using an overlap heuristic. ISPRS J. Photogramm. 2015, 102, 96–109. [Google Scholar] [CrossRef]
  15. Zhao, L.; Xiang, Z.; Chen, M.; Ma, X.; Zhou, Y.; Zhang, S.; Hu, C.; Hu, K. Establishment and Extension of a Fast Descriptor for Point Cloud Registration. Remote Sens. 2022, 14, 4346. [Google Scholar] [CrossRef]
  16. Deng, X.; Wang, S.; Zhuang, Y.; Fan, Y.; Zhou, Y. Evaluation of new target structure and recognition for point cloud registration and coordinates transformation of China’s large double-span bridge. J. Eng. Appl. Sci. 2023, 70, 135. [Google Scholar] [CrossRef]
  17. Singh, S.K.; Banerjee, B.P.; Raval, S. Three-dimensional unique-identifier-based automated georeferencing and coregistration of point clouds in underground mines. Remote Sens. 2021, 13, 3145. [Google Scholar] [CrossRef]
  18. Ge, X.; Zhu, Q. Target-based automated matching of multiple terrestrial laser scans for complex forest scenes. ISPRS J. Photogramm. 2021, 179, 1–13. [Google Scholar] [CrossRef]
  19. Ge, X.; Hu, H.; Wu, B. Image-guided registration of unordered terrestrial laser scanning point clouds for urban scenes. IEEE Trans. Geosci. Remote Sens. 2019, 57, 9264–9276. [Google Scholar] [CrossRef]
  20. Dong, Z.; Liang, F.; Yang, B.; Xu, Y.; Zang, Y.; Li, J.; Wang, Y.; Dai, W.; Fan, H.; Hyyppä, J. Registration of large-scale terrestrial laser scanner point clouds: A review and benchmark. ISPRS J. Photogramm. 2020, 163, 327–342. [Google Scholar] [CrossRef]
  21. Guo, Y.; Sohel, F.; Bennamoun, M.; Wan, J.; Lu, M. An accurate and robust range image registration algorithm for 3D object modeling. IEEE Trans. Multimed. 2014, 16, 1377–1390. [Google Scholar] [CrossRef]
  22. Theiler, P.W.; Wegner, J.D.; Schindler, K. Keypoint-based 4-points congruent sets–automated marker-less registration of laser scans. ISPRS J. Photogramm. 2014, 96, 149–163. [Google Scholar] [CrossRef]
  23. Tang, Y.; Feng, J. Hierarchical multiview rigid registration. Comput. Graph. Forum 2015, 34, 77–87. [Google Scholar] [CrossRef]
  24. Remondino, F.; Spera, M.G.; Nocerino, E.; Menna, F.; Nex, F. State of the art in high density image matching. Photogramm. Rec. 2014, 29, 144–166. [Google Scholar] [CrossRef]
  25. Schenk, T. Digital Photogrammetry; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  26. Fras, M.K.; Vezocnik, R.; Gvozdanovic, T.; Kogoj, D. Complete automation of the relative orientation of a stereopair. Geod. Vestn. 2008, 52, 241–253. [Google Scholar]
  27. Tjahjadi, M.E. Photogrammetric area-based least square image matching for surface reconstruction. Elektrika 2017, 1, 22–26. [Google Scholar]
  28. Shi, J.; Malik, J. Normalized cuts and image segmentation. IEEE Trans. Pattern Anal. Mach. Intell. 2000, 22, 888–905. [Google Scholar]
  29. Anon. sklearn.cluster. KMeans. Available online: https://scikit-learn.org/stable/modules/generated/sklearn.cluster.KMeans.html (accessed on 14 January 2024).
  30. Chin, T.-J.; Suter, D. The Maximum Consensus Problem: Recent Algorithmic Advances; Springer Nature: Berlin/Heidelberg, Germany, 2022. [Google Scholar]
  31. Besl, P.J.; McKay, N.D. Method for registration of 3-D shapes. Proc. SPIE 1992, 1611, 586–606. [Google Scholar] [CrossRef]
  32. Lu, F.; Milios, E. Globally consistent range scan alignment for environment mapping. Auton. Robot. 1997, 4, 333–349. [Google Scholar] [CrossRef]
  33. Anon. Leica ScanStation P40/P30—High-Definition 3D Laser Scanning Solution. Available online: https://leica-geosystems.com/products/laser-scanners/scanners/leica-scanstation-p40--p30 (accessed on 14 January 2024).
  34. Anon. User Manuals and Quick Start Guides for the Focus Laser Scanner. Available online: https://knowledge.faro.com/Hardware/Focus/Focus/User_Manuals_and_Quick_Start_Guides_for_the_Focus_Laser_Scanner (accessed on 14 January 2024).
  35. Wei, P.; Yan, L.; Xie, H.; Huang, M. Automatic coarse registration of point clouds using plane contour shape descriptor and topological graph voting. Autom. Constr. 2022, 134, 104055. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed method.
Figure 1. Flowchart of the proposed method.
Sensors 24 01394 g001
Figure 2. An example of the binary image: (a) Traffic direction estimation. (b) Binary image.
Figure 2. An example of the binary image: (a) Traffic direction estimation. (b) Binary image.
Sensors 24 01394 g002
Figure 3. An example of image matching: (a) The relative positions of a scan with respect to the registration template. (b) The relative positions between different scans.
Figure 3. An example of image matching: (a) The relative positions of a scan with respect to the registration template. (b) The relative positions between different scans.
Sensors 24 01394 g003
Figure 4. Repeated geometric patterns in a bridge.
Figure 4. Repeated geometric patterns in a bridge.
Sensors 24 01394 g004
Figure 5. Candidate positions of a scan.
Figure 5. Candidate positions of a scan.
Sensors 24 01394 g005
Figure 6. An example of the scan-block construction result.
Figure 6. An example of the scan-block construction result.
Sensors 24 01394 g006
Figure 7. An example of coarse registration by the PSO algorithm: (a) Before PSO. (b) After PSO.
Figure 7. An example of coarse registration by the PSO algorithm: (a) Before PSO. (b) After PSO.
Sensors 24 01394 g007
Figure 8. Two bridges in the bridge datasets: (a) Cuntan Yangtze River Bridge. (b) Huanghuayuan Jialing River Bridge.
Figure 8. Two bridges in the bridge datasets: (a) Cuntan Yangtze River Bridge. (b) Huanghuayuan Jialing River Bridge.
Sensors 24 01394 g008aSensors 24 01394 g008b
Figure 9. Initial pose estimation results of two bridges: (a) Cuntan Yangtze River Bridge. (b) Huanghuayuan Jialing River Bridge.
Figure 9. Initial pose estimation results of two bridges: (a) Cuntan Yangtze River Bridge. (b) Huanghuayuan Jialing River Bridge.
Sensors 24 01394 g009
Figure 10. Registration results of the bridge dataset: (a) Cuntan Yangtze River Bridge. (b) Huanghuayuan Jialing River Bridge.
Figure 10. Registration results of the bridge dataset: (a) Cuntan Yangtze River Bridge. (b) Huanghuayuan Jialing River Bridge.
Sensors 24 01394 g010
Table 1. The detailed description of the two datasets.
Table 1. The detailed description of the two datasets.
DatasetScannersScansPts (Billion)
Cuntan Yangtze River BridgeLeica P40
(Leica, Wetzlar, Germany)
291.52
Huanghuayuan Jialing River BridgeFaro S350
(Faro, Lake Mary, FL, USA)
201.68
Table 2. Evaluation on template-guided initial pose estimation.
Table 2. Evaluation on template-guided initial pose estimation.
DatasetRotation Error (mdeg)Translation Error (m)Time (min)
ΔxΔyΔz
Cuntan Yangtze River Bridge19.10.663.621.886.08
Huanghuayuan Jialing River Bridge19.70.645.162.354.43
Table 3. Evaluation of the proposed method.
Table 3. Evaluation of the proposed method.
DatasetRotation Error (mdeg)Translation Error (mm)SSR (%)Time (min)
AverageRMSEAverageRMSE
Cuntan Yangtze River Bridge0.960.6728.0414.1310050.6
Huanghuayuan Jialing River Bridge0.740.6643.2528.4810040.5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xiong, G.; Cui, N.; Liu, J.; Zeng, Y.; Chen, H.; Huang, C.; Xu, H. Template-Guided Hierarchical Multi-View Registration Framework of Unordered Bridge Terrestrial Laser Scanning Data. Sensors 2024, 24, 1394. https://doi.org/10.3390/s24051394

AMA Style

Xiong G, Cui N, Liu J, Zeng Y, Chen H, Huang C, Xu H. Template-Guided Hierarchical Multi-View Registration Framework of Unordered Bridge Terrestrial Laser Scanning Data. Sensors. 2024; 24(5):1394. https://doi.org/10.3390/s24051394

Chicago/Turabian Style

Xiong, Guikai, Na Cui, Jiepeng Liu, Yan Zeng, Hanxin Chen, Chengliang Huang, and Hao Xu. 2024. "Template-Guided Hierarchical Multi-View Registration Framework of Unordered Bridge Terrestrial Laser Scanning Data" Sensors 24, no. 5: 1394. https://doi.org/10.3390/s24051394

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop