Next Article in Journal
Guava Detection and Pose Estimation Using a Low-Cost RGB-D Sensor in the Field
Next Article in Special Issue
Design and DOF Analysis of a Novel Compliant Parallel Mechanism for Large Load
Previous Article in Journal
Chemical Source Searching by Controlling a Wheeled Mobile Robot to Follow an Online Planned Route in Outdoor Field Environments
Previous Article in Special Issue
Modelling and Design of MEMS Piezoresistive Out-of-Plane Shear and Normal Stress Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

HEALPix-IA: A Global Registration Algorithm for Initial Alignment

State Key Laboratory of Robotics and System, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(2), 427; https://doi.org/10.3390/s19020427
Submission received: 12 December 2018 / Revised: 16 January 2019 / Accepted: 18 January 2019 / Published: 21 January 2019
(This article belongs to the Special Issue Sensors for MEMS and Microsystems)

Abstract

:
Methods of point cloud registration based on ICP algorithm are always limited by convergence rate, which is related to initial guess. A good initial alignment transformation can sharply reduce convergence time and raise efficiency. In this paper, we propose a global registration method to estimate the initial alignment transformation based on HEALPix (Hierarchical Equal Area isoLatitude Pixelation of a sphere), an algorithm for spherical projections. We adopt EGI (Extended Gaussian Image) method to map the normals of the point cloud and estimate the transformation with optimized point correspondence. Cross-correlation method is used to search the best alignment results in consideration of the accuracy and robustness of the algorithm. The efficiency and accuracy of the proposed algorithm were verified with created model and real data from various sensors in comparison with similar methods.

1. Introduction

A point cloud is a set of data points in space. A 3D point cloud is generally produced by sensors such as 3D scanner, 3D camera, Light Detection and Ranging equipment (LiDAR) and so on. Although accuracy and density from the sensors above are different [1], similar analysis methods, such as conversion to 3D surfaces, alignment and registration, are processing in different fields. 3D point cloud registration or alignment is a fundamental problem in computer and robotic vision. The main task is to find the transformation describing rigid motion in SE(3) (Lie group for 3D Rigid transformations) between two given sets of the same or similar point clouds through the analysis of local or global feature. This technique is widely applied to object pose estimation [2], simultaneous localization and mapping (SLAM) [3], medical image processing [4], industrial manufacturing and inspection [5,6], culture heritage reconstruction [7,8,9], topography change detection [10], etc. With the development of 3D scanning sensors and reconstruction technology, the accuracy of point cloud data has improved while the cost of scanning sensor has dropped. The point cloud registration technique is gradually entering low cost application fields, especially industrial areas. For the production mode of small lots and multiple categories, it is important to apply industrial robots and other flexible manufacturing equipments into grinding and polishing in consideration of poor working condition, increasing human cost and frequent product changes. However, the robot system needs a whole process with information such as CAD model, machining allowance (MA), tool path planning and process, and process inspection. Among these processes, 3D point cloud registration technique can be applied to MA analysis and process inspection. When a workpiece is scanned, point cloud data are acquired to align with the CAD model. The redundant parts are calculated as machining allowance, which is taken as input for tool path planning and process. On the other hand, registration is applied to calculate the error between workpiece and CAD model in process inspection.
The application mentioned above in industrial field has used coordinate measuring machine (CMM) to acquire sparse patterns, which are used for localization and MA analysis in early studies [5,11,12], and even is still the mainstream method today. Nevertheless, the CMM method is limited by shape of the workpiece. Complex curve surfaces need plenty of time for measuring and recording or even cannot be dealt with. With the development of laser scanners, stereo cameras and other same type of sensors, the density of points is sharply increasing, as we call the dataset “point cloud”. The improved accuracy of points and reconstruction makes it more precise for registration, but more complex. Researchers have proposed different methods for MA registration with dense point cloud. Li, X. et al. [13] introduced a novel two-step rough–precise automatic method. A PCA-based (principal component analysis) rough alignment algorithm is utilized with the precise registration method, which only considers plane feature matching. Wang, H. et al. [14] introduced a heuristic algorithm into PCA/ICP registration method with multi-objective optimization considering envelopment, wall thickness and plane alignment. Sun, Y. et al. [15] combined the multipliers method and the BFGS algorithm to handle many constraints. The objective function is driven by Euclidean oriented distance with the assumption of a small initial rigid transformation and only envelopment constraint is considered. Whatever the optimization constraints, the methods above have been processed by two steps, rough alignment and optimal registration. Since dense point cloud data are acquired, Sac-IA [16], PCA [17,18,19], 3D-NDT [20] and EGI-based alignment algorithms [21,22,23,24,25] are common methods for rough alignment in different applications. Optimal registration methods are always based on ICP algorithm with different process requirements. Generally, the optimal registration algorithms mentioned above are all local methods, which could be trapped into local minima easily [26]. Thus, a good initial transformation could accelerate the convergence and improve accuracy of the optimization method significantly.
We propose a fast global registration method based on EGI mapping to resolve the contradiction between efficiency and accuracy of dense point cloud registration. The main algorithm includes several steps: (1) equal-area pixelization; (2) EGI mapping; (3) point correspondence searching; (4) transformation estimating; and (5) correlation optimization. The point correspondence searching and correlation optimization could achieve a high accuracy within several iterations in our experiments. The proposed method was compared with other initial alignment registration methods, such as Sac-IA, PCA, 3D-NDT and EGI-based methods. Several contributions are made in this paper: First, a new division method is suggested to improve the division efficiency. Second, a point corresponding optimization is proposed to improve the accuracy. Finally, the application of the algorithm to MA analysis is presented.
The rest of the paper is organized as follows: First, the related works on other initial alignment registration methods are introduced in Section 2. The problem is expressed in Section 3. The main strategy, HEALPix-IA Algorithm, is discussed in Section 4. In Section 5, experiments are presented to demonstrate the advantages of this algorithm over other rough registration methods. Section 6 concludes the paper.

2. Related Works

Various rough registration methods are suggested for initial alignment to solve the local minima problem within a couple of seconds or even more quickly.
PCA-based algorithm can perform object alignment in real-time and without constraints on the three registration parameters (i.e., translation, rotation, and scaling). Principal Components Analysis (PCA) computes the mutually orthogonal directions of maximum variance in a collection of d-dimensional data and measures the global features by eigenvalues. It is widely used for image registration [27,28] and 3D alignment [17,18,19]. However, PCA-based methods are sensitive to noise [29], which limits accuracy.
The Sample Consensus-Initial Alignment (SAC-IA) algorithm was introduced by Rusu et al. [16,30], using 16-dimensional point feature histograms (FPFH) that describe the local surface structure. FPFH gives a good discriminative power for point correspondence search. Nevertheless, the robustness to outliers, invariance to pose, sampling density and measurement noise are not mentioned and the extensive computational steps result in time-consuming process.
Another 3D point set registration algorithm is the normal distributions transform (NDT) [20] that represents the underlying scene geometry as a Gaussian probability distribution. 3D-NDT method generates disjoint voxels in space and represents points within the voxels as a probability density function (PDF). Its benefit is to give piece-wise smooth spatial representations; however, the division results in discontinuities in the cost function that could be trapped in local minima [31].
EGI-based algorithms are another branch of 3D point registration, which describes the global features by Extended Gaussian Images (EGI) of laser scan. The intensity image mapping from normals is invariant to translation, thus such algorithms have the advantage that pose parameters can be recovered in two separate steps. Latitude and longitude tessellation [21], regular polyhedron tessellation [22], and other complex EGI division methods [23,24,25] have been suggested to raise the accuracy. A voting/correlation procedure, such as spherical correlation, cross correlation and kernel correlation [32] is generated to search the best alignment. However, the accuracy and efficiency of these algorithms are limited by the division method [33] and it is not sensitive to constant EGI (such as a sphere) [34].
In this paper, an EGI-based registration algorithm is introduced with new division method, which is equal area pixelization. We suggest a special procedure for point correspondence optimization to reduce the error introduced by low level division method. In addition, several conditions are configured to avoid the mismatch caused by constant EGI.

3. Problem Formulation

In the standard point-to-point ICP algorithm, registration problem is defined with two 3D point-sets X = { x i } , i = 1 , , N and Y = { y j } , j = 1 , , M , where x i , y i R are point coordinates. We call X model point cloud and Y data. The registration problem is to estimate a rigid motion transformation with translation t R 3 and rotation R S O ( 3 ) (Lie group for 3D rotations), which minimizes the error E:
E ( R , t ) = i = 1 N e i ( R , t ) 2 = i = 1 N R x i + t y j * 2 ,
where e i ( R , t ) 2 is the per-point residual error for x i , and · denotes the Euclidean norm. Given r and t , the point y j * Y is denoted as the optimal correspondence of x i , which is the closest point to the transformed x i in Y , i.e.,
j * = arg min j { 1 , , M } R x i + t y j * .
According to ICP algorithm, two prominent cost functions are used to measure “resemblance”: Hausdorff distance and root mean square distances (RMS) [32]. In this paper, we discuss the registration application that data and model are not totally resembled because of machining allowance. Hausdorff distance shows the measure of maximal error and RMS, which is more suitable for the “under-resembled” registration problem, shows the average measure. We used the notation
R M S ( R , t ) = 1 N i = 1 N R x i + t y j * 2 = 1 N E ( R , t ) ,
for the evaluation of various methods in experiments.
In this paper we also define the registration problem with X , Y , R and t . The proposed method is applied to the rough alignment stage of the whole registration process, so that global measure is selected in consideration of convergence region, noise sensitivity and robustness. The histogram is brought in as a global feature notation to simplify the local features of the point set, as the notation H ( X ) denotes the feature histogram of data point cloud. As correlation-based method is demonstrated to outperform ICP algorithm [32], we defined the objective function to minimize the following cost function,
COST ( X , Y , R , t ) = C o r r ( H ( X ) , H ( T ( Y , R , t ) ) ,
where C o r r ( · ) is the correlation function. Since two point clouds are resembled, t is the difference of the centroids, which are easily calculated with the coordinates of the points. Ignoring the slight error of t , the problem is transformed into solving the value of R .

4. HEALPix-IA Algorithm

We call the algorithm HEALPix-IA because it is an initial alignment method based on HEALPix. To minimize the cost function in Equation (4), several steps need to be followed in order. In this section, HEALPix is introduced to reduce the dimensions by mapping the local features to 2-sphere first. After getting the histogram, corresponding points are estimated and optimized. The rotation can be estimated with the corresponding points and optimized with the correlation function.

4.1. Pixelization, Projection and Indexing

HEALPix (Hierarchical Equal Area isoLatitude Pixelization of a 2-sphere) is an algorithm for pixelization of the 2-sphere, and the associated class of map projections [35].
It was first used for satellite missions to measure the cosmic microwave background anisotropy. The HEALPix map was created for the construction of full-sky maps of the microwave sky at an angular resolution of a few arcminutes. The sphere is hierarchically tessellated into curvilinear quadrilaterals. The lowest resolution partition is comprised of 12 base pixels. Resolution of the tessellation increases by division of each pixel into four new ones. Figure 1 illustrates (clockwise from upper-left to bottom-left) the resolution increase by three steps from the base level (i.e., the sphere is partitioned, respectively, into 12, 48, 192, and 768 pixels). We use the notation
m = 12 × 4 n 1 , n = 1 , 2 , 3 , ,
where m denotes the number of the pixel under nth level of the division.
The scheme has three advantages for our method: (1) Equal-area mapping can reduce the normal feature error generated in different posture (SO(3)). (2) Areas of all pixels at a given resolution are identical. (3) It has a less complexity computation of integrals over individual spherical harmonics. This gives a strong base for high speed performance of the proposed algorithm.
For the point sets data and model, let the scale be equalled by voxel filter to simplify the projection. In consideration of the efficiency and accuracy, let n = 5 and the tessellation be divided into 3072 grids. As we obtain the EGI or normal sphere (NS) of the point cloud, it can be projected to HEALPix sphere (HS). The HealpixLib is utilized to generate sampling points and convert them from Spherical Coordination to Cartesian Coordination. The projection is equivalent to finding the nearest sample point of the normal vector using k-nearest neighbors algorithm. Thus, the projection from a single normal to the corresponding Healpix sample point is defined as
P x = P ( N s ) .
In this paper, we give hierarchic projection method to accelerate accessing data, which is different from other EGI methods. Let n = 1 in Equation (5); we call the 12 pixels basic girds. In every basic grid, there are 256 subdivided grids. Here, we create the index using the notation
H S = { H 1 S , , H i S , , H 12 S } , H T = { H 1 T , , H i T , , H 12 S } ,
where H i S and H i T are the histograms of the 12 basic grids. Each histogram has 256 values, which stand for the intensity of subdivided grids. The mapping relation from the normals to histogram is defined as
H x = H ( N s ) ,
where N s denotes a group normals and H x denotes the intensive histogram of the corresponding Healpix grids.

4.2. Point Correspondence Searching

Corresponding points are used to estimate the rotation. Generally, we use the peaks in the histogram to estimate at least two corresponding points. However, two situations would make this method fail. One is the histogram is relatively even in some range. The other is the two corresponding points are opposite (nearly 180 degree) in EGI. To tackle this problem, we use a relative feature instead of the absolute coordinate of the corresponding points.
First, finding the maximal grid in H i S and H i T , we get
P S = { P 1 S , , P i S , , P 12 S } , a n d P T = { P 1 T , , P i T , , P 12 S } ,
where P i S and P i T are the indices of the maximal intensity and sorted in intensity descending order. Then, we let the center point of the sphere be the origin and get the vectors V i S , V j S , V i T and V j T point to the sample point of P i S , P j S , P i T and P j T , where P j S and P j T are the indices of the maximal intensity in H j S and H j T . For P i S , P j S must meet the condition
arccos ( V i S , V j S ) < 0.94 , i < j 12 ,
or P j S is excluded. The parameter 0.94 is an empirical value less than 1, as the two vectors are not absolutely coincident. The same condition is met with V i T and V j T . Distance between the rest grids is calculated and denote as
D S = { d i j S P i S , P j S P S } , a n d D T = { d i j T P i T , P j T P T } ,
where d i j S and d m n T are the distance between each two rest grids. If i , j , let
ε = d i j S d i j T < 0.05 ,
then subdivided grids P i S and P j S can be aligned with P i T and P j T initially. The error 0.05 is an empirical value according to the difference and accuracy of the two data.

4.3. Point Correspondence Optimization

The Healpix method is limited by its accuracy for the pixelization. Although the sphere is divided into over 3000 grids, there is still intolerable error up to four degrees. To improve the corresponding point alignment accuracy, we optimize the method to calculate the peak of the intensive region in geometry instead of the sample points of Healpix grids.
In the last subsection, the two pairs of corresponding points are P 1 S , P 2 S , P 1 T and P 2 T . The normal set is defined as
N p = { N q P ( N q ) { P 0 , , P n } }
where P n are the n neighbours of the subdivided Healpix grid P 0 . P 0 is one of the corresponding points mentioned above. Then, a plane A is created through the origin of the Healpix sphere, perpendicular to the normal vector of P 0 . The projection of N p is defined on the plane A, as N p R 2 . In this optimal method, the intensity peak in plane A is searched by KD-tree using a circle with fixed radius, which is denoted as
r = max x i x j τ ,
where τ is the parameter to control radius, and numerator is the range of the normals in x axis. In the range of radius r, we mesh the region with grid in fixed step, and search for the maximum intensity grid with KD-tree method. Assuming ( x p , y p ) is the peak point, it is projected back to Healpix sphere, denoted as P 0 , which is the modified corresponding point. Thus, we modify the corresponding points as P 1 S , P 2 S , P 1 T and P 2 T . Algorithm 1 describes how to estimate and optimize the corresponding points.

4.4. Rotation Estimation

To align two EGIs, we assume that the origins coincide, so that only two rotation angles are needed for entire alignment. Defining the normal vectors of the corresponding points V i S , and V i T , the rotation vector is
V i = V i T × V i S V i T × V i S ,
and rotation angle is
θ i = arccos ( V i T × V i S ) V i T · V i S .
According to Rodrigues’ rotation formula, the rotation matrix is
R i = cos ( θ i ) I + ( 1 cos ( θ i ) ) V i V i H + sin ( θ i ) 0 V i z V i z V i z 0 V i x V i y V i x 0 ,
where I is identity matrix, V i H is the transposed matrix, and V i x , V i y , and V i z are the components of V i . Then, we calculate the final rotation matrix with Equation (17),
R = R 1 · R 2 .
Algorithm 1: Estimate and optimize the corresponding points.
input: ModelS and DataT with normals of size N × 6 (N is the scale of Model/Data)
output:  2 pairs of corresponding points P 1 S , P 2 S , P 1 T and P 2 T
1 extract normals N S and N T from S and T
2 H S H ( N S ) ; H T H ( N T ) ;// see Equations (7) & (8)
3for i 1 to 12 do
4 H S ( i ) , I S ( i ) max ( H i S ) ;// maximum intensity and the index
5 H T ( i ) , I T ( i ) max ( H i T )
6 N S ( i ) I S ( i ) ; N T ( i ) I T ( i ) ;// normals of the max
7 P i S P ( N S ( i ) ) ; P i T P ( N T ( i ) ) ;// see Equation (6)
8 construct Struct Q ( i ) = { H ( i ) , I ( i ) , P i }
9end for
10 Q S sortrows ( H S ; Q T ) sortrows ( H T );// descending sort
11 i = 2
12 P t e m p S = P 2 S
13while arccos ( V i S , V t e m p S ) > 0.94 do// see Eqution (10)
14 P t e m p S = P i + 1 S
15 i i + 1
16end while
17 P 2 S = P t e m p S
18for i 2 to 12 do
19 P t e m p T = P i T
20if d 1 t e m p S d 1 t e m p T < 0.05 then// see Equation (12)
21   P 2 T = P t e m p T
22  break;
23end if
24end for
25let n = 4 then calculate N 1 S , N 2 S , N 1 T , N 2 T ;// see Equation (13)
26 a , b , c = V 1 S
27create e 1 = ( 1 , 1 , a + b c ) in plane A then e 2 = e 1 × V 1 S
28 N 1 S = ( e i ; e 2 · ( N 1 S ) T ) T ;// project normals to A
29calculate r 1 S ;// see Equation (14)
30 mesh r 1 S × r 1 S region into 100 grids;
31search ( x p , y p ) with the peak intensity in meshed grids;
32calculate inverse projection point P 1 S
33calculate P 2 S , P 1 T , P 2 T

4.5. Correlation Searching

According to Equation (4), cost function is employed to maximize the correlation of the model and transformed data. In this method, we use cross correlation method and simply define the correlation function as
C o r r = H S , H T H S · H T .
When C o r r 1 , the model and transformed data are most correlative. Iterative search is carried out with halt conditions in this method. Moreover, parameter correlation and iteration are recommended as 0.98 and 50.
In the process of the optimizing search, random rotation angle θ is given in rotation matrix form
T r a n d = R ω ( θ ) R ϕ ( θ ) R κ ( θ ) ,
where ω , ϕ , κ ( π , π ) . The basic rotations follow the right-hand rule. Algorithm 2 describes the search strategy of transformation estimation.
Algorithm 2: Optimal transformation search.
input: ModelS and DataT with normals of size N × 6 (N is the scale of Model/Data)
output: transformation T, correlation c o r r , iteration i t e r
1initialize i t e r = 0 ; s t a l l = 0 ; c o r r = 0 ; T = e y e ( 3 )
2extract normals N S and N T from S and T
3 H S H ( N S ) ; // see Equation (8)
4while c o r r < 0.98 and i t e r < 50 do
5 estimate initial transformation T t e m p with Equations (15)∼(18)
6 H T H ( ( T t e m p · ( N T ) T ) T )
7 calculate c o r r t e m p ;// see Equation (19)
8 s t a l l s t a l l + 1
9if c o r r < c o r r t e m p then
10  if c o r r t e m p > c o r r × 1.005 then
11    s t a l l s t a l l 1
12  end if
13   c o r r c o r r t e m p
14   T T t e m p · T
15   N T ( T t e m p · ( N T ) T ) T
16end if
17if s t a l l 3 and c o r r t e m p < 0.98 then
18   T T r a n d · T ; // see Equation (20)
19   N T ( T r a n d · ( N T ) T ) T
20   s t a l l 0
21end if
22 i t e r i t e r + 1
23end while

5. Results and Evaluation

The proposed algorithm was verified with both model simulation and real data. It was compared with other EGI-based methods and other sophisticated rough registration algorithms.

5.1. Simulation

Three tests were carried out in simulation. The test method was self-registration with random rotations. Different division methods and point correspondence optimization were verified in the first test. The performance of the proposed algorithm was compared with other rough registration methods in the second test. How different features affect on the method was tested in the end. The four CAD models shown in Figure 2 were prepared to test HEALPix-IA algorithm. They were all casting parts that need grinding after molding.

5.1.1. Compared with EGI-Based Methods

In the proposed algorithm, new division method and point correspondence method are suggested. To measure the effectivity of the algorithm, root mean square distances (RMS) and running time are two key evaluation criteria. We used several different types of models to test the registration results with two EGI-based algorithms. Latitude and longitude tessellation (Lati-Longi, in Ameesh’s algorithm [21]) and regular polyhedron tessellation (Reg-Poly, in Christoph’s algorithm [22]) were both tested in this experiment. In addition, point correspondence optimization (PC-Opt) was also separated. Six trials included HEALPix tessellation, Lati-Longi tessellation, Reg-Poly tessellation, HEALPix-IA, and Lati-Longi, respectively, with point correspondence optimization and Reg-Poly with point correspondence optimization, carried out with two models and 50 samples. The results are shown in Figure 3, where RMS and running time are both in logarithm. The tests were configured with the same halt conditions: 0.02 correlation value and 50 iterations.
The two selected models were water turbine blade and casting parts 2 with typical features friendly to EGI-based methods. The first three groups without PC-Opt in the boxplot show HEALPix division method won in both RMS and running time, while the three groups with PC-Opt performed more distinctly. The comparison between each division method with or without PC-Opt showed that PC-Opt gave better performance for HEALPix and Lati-Longi. HEALPix and Lati-Longi overwhelmed Reg-Poly in running time because the letter one uses an iteration division method that costs time for every generation. Furthermore, the higher level the sphere is divided, the more time is cost. This is the main limit of the regular polyhedron tessellation. The methods with PC-Opt ran more accurately than the ones with only divide tessellations. With the same RMS for two methods with PC-Opt, the running time of the proposed algorithm had an advantage over other methods.

5.1.2. Compared with Other Rough Registration Methods

With PCL library, PCA, 3D-NDT, and Sac-IA were applied, respectively, on the registration tasks of four CAD models together with HEALPix-IA (C + + ). One hundred trials were carried out for each model. The results are shown in Figure 4 and Figure 5. The RMS values are shown in logarithm.
From the boxplot, we can see HEALPix-IA performed the best in both accuracy and efficiency and PCA performed the worst. For accuracy, HEALPix-IA and Sac-IA obtained about the same RMS; however, Sac-IA had a more stable performance because it had over 1000 searching iterations, which resulted in long running time. In addition, HEALPix-IA won by almost an order of magnitude in running time over Sac-IA. Nevertheless, 3D-NDT’s results depended on its complex parameter configurations, which would hinder its application in automatic process. Although it ran quickly, the accuracy was quite low. Furthermore, 3D-NDT is not safe because it is a local method that is easily trapped into local minima. The three other global methods are better choices to carry out an automatic rough registration task.

5.1.3. Tests on Real Data with Different Features via Various Sensors

We chose nine different data, shown in Figure 6, to test the proposed algorithm. Three kinds of sensors were used in this test: 3D laser scanner, 3D cameras and 3D lidar (Figure 7). The 3D laser scanner can reach the volumetric accuracy up to 30 μm and the 3D camera with image stitching technique can reach about 0.8 mm. All data obtained from the sensors had their noise eliminated. The different features are explained in Table 1 and the test results are shown in Figure 8.
In this test, we raised the search iteration to 50 generations to adapt to different types of features. The boxplot shows that HEALPix-IA performed poorly with “human arm dummy” and “laser cut rocket”. The reason is that they both have the constant EGI features similar to sphere surface, which is a disadvantage of all EGI-based methods. The “casting parts 1” with thin-walled plane led to a fast convergence in searching iteration, so that it performed a quick registration with unstable accuracy. Besides, the two other sensors achieved better results in both accuracy and efficiency. In summary, HEALPix-IA is applicable with most features and different sensors, yet its stability needs to be improved further.

5.2. MA Analysis Test

In this paper, we suggest HEALPix-IA is applicable to MA analysis for manufacturing. In this experiment, two real workpieces were scanned by 3D scanner (Creaform HandyScan 700) and aligned to the CAD models. As the rough registration is the first step for MA analysis, another evaluation criterion, optimization time for MA analysis process, was used. The results are shown in Table 2 with Sac-IA and HEALPix-IA. Here, we abandoned PCA and 3D-NDT for stability and local method. HEALPix-IA was verified to perform relatively better in accuracy (RMS) and efficiency (time) for MA analysis application with almost the same stability though the variance values. The MA analysis results are shown in Figure 9 which can be used for follow-up manufacturing processes such as debugging and grinding.

6. Conclusions

In this paper, a new rough registration method—HEALPix-IA—is proposed especially for MA analysis application. HEALPix-IA is based on EGI, and has the same features with other EGI-based algorithms, such as global method, robust to noise, and sensitive to constant EGI. Furthermore, we suggest a new division method, which was proven to be more effective. Point correspondence optimization is also implemented to improve the accuracy. HEALPix-IA was compared with other sophisticated rough registration algorithms, PCA, 3D-NDT and Sac-IA, on RMS and running time, utilizing nine real scanned data and two pairs of datasets for MA analysis, which were obtained from three different kinds of sensors. The experiments verified HEALPix-IA shows a better performance on accuracy and efficiency.

Author Contributions

Conceptualization, Y.G.; methodology, Y.G.; software, W.X.; validation, M.L.; investigation, Y.G.; writing—original draft preparation, Y.G.; writing—review and editing, W.D.; supervision, Z.D.; and project administration, Z.D.

Funding

This research was funded by State Administration for Science, Technology and Industry for National Defense grant number JSCG2015603B003.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study.

Abbreviations

The following abbreviations are used in this manuscript:
ICPIterative Closest Point
EGIExtended Gaussian Image
SLAMSimultaneous Localization And Mapping
MAMachining Allowance
CMMcoordinate measuring machine
PCAprincipal component analysis
FPFHFast Point Feature Histograms
Sac-IASample Consensus-Initial Alignment
NDTNormal Distributions Transform
RMSRoot Mean Square distances

References

  1. Balsa-Barreiro, J.; Avariento, J.P.; Lerma, J.L. Airborne light detection and ranging (LiDAR) point density analysis. Sci. Res. Essays 2012, 7, 3010–3019. [Google Scholar] [CrossRef]
  2. Choi, C.; Taguchi, Y.; Tuzel, O.; Liu, M.Y.; Ramalingam, S. Voting-based pose estimation for robotic assembly using a 3D sensor. In Proceedings of the IEEE International Conference on Robotics and Automation, St. Paul, MN, USA, 14–18 May 2012; pp. 1724–1731. [Google Scholar]
  3. Holz, D.; Behnke, S. Sancta Simplicitas—On the efficiency and achievable results of SLAM using ICP-based incremental registration. In Proceedings of the 2010 International Conference on Robotics and Automation, Anchorage, AK, USA, 4–8 May 2010; pp. 1380–1387. [Google Scholar]
  4. Miao, S.; Zhang, W.J.; Liao, R. A CNN regression approach for real-time 2D/3D registration. IEEE Trans. Med. Imaging 2016, 35, 1352–1363. [Google Scholar] [CrossRef] [PubMed]
  5. Li, Z.; Gou, J.; Chu, Y. Geometric algorithms for workpiece localization. IEEE Trans. Robot. Autom. 1998, 14, 864–878. [Google Scholar]
  6. Tan, G.; Zhang, L.; Liu, S.; Zhang, W. A fast and differentiated localization method for complex surfaces inspection. Int. J. Precis. Eng. Manuf. 2015, 16, 2631–2639. [Google Scholar] [CrossRef]
  7. Shanoer, M.M.; Abed, F.M. Evaluate 3D laser point clouds registration for cultural heritage documentation. Egypt. J. Remote Sens. Space Sci. 2018, 21, 295–304. [Google Scholar] [CrossRef]
  8. Balsa-Barreiro, J.; Fritsch, D. Generation of visually aesthetic and detailed 3D models of historical cities by using laser scanning and digital photogrammetry. Digit. Appl. Archaeol. Cult. Herit. 2018, 8, 57–64. [Google Scholar] [CrossRef]
  9. Owda, A.; Balsa-Barreiro, J.; Fritsch, D. Methodology for digital preservation of the cultural and patrimonial heritage: Generation of a 3D model of the Church St. Peter and Paul (Calw, Germany) by using laser scanning and digital photogrammetry. Sens. Rev. 2018, 38, 282–288. [Google Scholar] [CrossRef]
  10. Lague, D.; Brodu, N.; Leroux, J. Accurate 3D comparison of complex topography with terrestrial laser scanner: Application to the Rangitikei canyon (N-Z). ISPRS J. Photogramm. Remote Sens. 2013, 82, 10–26. [Google Scholar] [CrossRef] [Green Version]
  11. Chu, Y.X.; Gou, J.B.; Li, Z.X. On the Hybrid Workpiece Localization/Envelopment Problem. IEEE Int. Conf. Robot. Autom. 1998, 4, 1–25. [Google Scholar]
  12. Gou, J.B.; Chu, Y.X.; Xiong, Z.H.; Li, Z.X. A Geometric Method for Computation of Datum Reference Frames. IEEE Trans. Robot. Autom. 2000, 16, 797–806. [Google Scholar] [CrossRef]
  13. Li, X.; Li, W.; Jiang, H.; Zhao, H. Automatic evaluation of machining allowance of precision castings based on plane features from 3D point cloud. Comput. Ind. 2013, 64, 1129–1137. [Google Scholar] [CrossRef]
  14. Wang, H.; Zhou, M.X.; Zheng, W.Z.; Shi, Z.B.; Li, H.W. 3D machining allowance analysis method for the large thin-walled aerospace component. Int. J. Precis. Eng. Manuf. 2017, 18, 399–406. [Google Scholar] [CrossRef]
  15. Sun, Y.W.; Xu, J.T.; Guo, D.M.; Jia, Z.Y. A unified localization approach for machining allowance optimization of complex curved surfaces. Precis. Eng. 2009, 33, 516–523. [Google Scholar] [CrossRef]
  16. Rusu, R.B.; Blodow, N.; Beetz, M. Fast Point Feature Histograms (FPFH) for 3D registration. In Proceedings of the 2009 IEEE International Conference on Robotics and Automation, Kobe, Japan, 12–17 May 2009; pp. 3212–3217. [Google Scholar]
  17. Weingarten, J.; Gruener, G.; Siegwart, R. A Fast and Robust 3D Feature Extraction Algorithm for Structured Environment Reconstruction; EPFL Scientific Publications: Geneva, Switzerland, 2003; pp. 1–6. [Google Scholar]
  18. Prasad, G. Segmentation of 3D MR Images of the Brain Using a PCA Atlas and Nonrigid Registration; University of California: Los Angeles, CA, USA, 2010; pp. 1–63. [Google Scholar]
  19. Schneider, D.C.; Eisert, P. Automatic and robust semantic registration of 3D head scans. In Proceedings of the 5th European Conference on Visual Media Production, London, UK, 26–27 November 2008. [Google Scholar] [CrossRef]
  20. Magnusson, M. The Three-Dimensional Normal-Distributions Transform—An Efficient Representation for Registration, Surface Analysis, and Loop Detection. Ph.D. Thesis, Örebro University, Örebro, Sweden, 2009. [Google Scholar]
  21. Makadia, A.A.; Patterson, V.; Daniilidis, K. Fully automatic registration of 3D point clouds. In Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, New York, NY, USA, 17–22 June 2006; pp. 1297–1304. [Google Scholar]
  22. Dold, C. Extended Gaussian images for the registration of terrestrial scan data. In Proceedings of the ISPRS Workshop “Laser Scanning 2005”, Enschede, The Netherlands, 12–14 September 2005; pp. 180–185. [Google Scholar]
  23. Kang, S.B.; Ikeuchi, K. The Complex EGI: A New Representation for 3-D Pose Determination. IEEE Trans. Pattern Anal. Mach. Intell. 1993, 15, 707–721. [Google Scholar] [CrossRef]
  24. Hilaga, M.; Shinagawa, Y.; Kohmura, T.; Kunii, T.L. Topology matching for fully automatic similarity estimation of 3D shapes. In Proceedings of the 28th Annual Conference on Computer Graphics and Interactive Techniques—SIGGRAPH’01, Los Angeles, CA, USA, 12–17 August 2001; pp. 203–212. [Google Scholar]
  25. Hu, Z.; Chung, R.; Fung, K.S. EC-EGI: Enriched complex EGI for 3D shape registration. Mach. Vis. Appl. 2010, 21, 177–188. [Google Scholar] [CrossRef]
  26. Ezra, E.; Sharir, M.; Efrat, A. On the performance of the ICP algorithm. Comput. Geom. 2008, 41, 77–93. [Google Scholar] [CrossRef] [Green Version]
  27. Huizinga, W.; Poot, D.H.; Guyader, J.M.; Klaassen, R.; Coolen, B.F.; Van Kranenburg, M.; Van Geuns, R.J.; Uitterdijk, A.; Polfliet, M.; Vandemeulebroucke, J.; et al. PCA-based groupwise image registration for quantitative MRI. Med. Image Anal. 2016, 29, 65–78. [Google Scholar] [CrossRef] [PubMed]
  28. Sasatani, S.; Han, X.H.; Chen, Y.W. Image Registration Using PCA And Gradient Method For Super-resolution Imaging. In Proceedings of the 2nd International Conference on Software Engineering and Data Mining, Chengdu, China, 23–25 June 2010; pp. 631–634. [Google Scholar]
  29. Zhan, Y.; Yin, J. Robust local tangent space alignment via iterative weighted PCA. Neurocomputing 2011, 74, 1985–1993. [Google Scholar] [CrossRef]
  30. Rusu, R.B.; Blodow, N.; Marton, Z.C.; Beetz, M. Aligning point cloud views using persistent feature histograms. In Proceedings of the 2008 IEEE/RSJ International Conference on Intelligent Robots and Systems, Nice, France, 22–26 September 2008; pp. 3384–3391. [Google Scholar]
  31. Magnusson, M.; Duckett, T. A Comparison of 3D Registration Algorithms for Autonomous Underground Mining Vehicles. In Proceedings of the European Conference on Mobile Robotics (ECMR 2005), Ancona, Italy, 7–10 September 2005; pp. 86–91. [Google Scholar]
  32. Tsin, Y.; Kanade, T. A Correlation-Based Approach to Robust Point Set Registration. In European Conference on Computer Vision; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  33. Straub, J.; Campbell, T.; How, J.P.; Fisher, J.W. Efficient global point cloud alignment using Bayesian nonparametric mixtures. In Proceedings of the 30th IEEE Conference on Computer Vision and Pattern Recognition—CVPR 2017, Honolulu, HI, USA, 21–26 July 2017; pp. 2403–2412. [Google Scholar]
  34. Ahuja, S.; Waslander, S.L. 3D scan registration using curvelet features. In Proceedings of the Conference on Computer and Robot Vision—CRV 2014, Montréal, QC, Canada, 7–9 May 2014; pp. 77–83. [Google Scholar]
  35. Calabretta, M.R.; Roukema, B.F. Mapping on the HEALPix grid. Mon. Not. R. Astron. Soc. 2007, 381, 865–872. [Google Scholar] [CrossRef] [Green Version]
Figure 1. HEALPix: The resolution increase by three steps from the base level (the figure is cited from https://healpix.jpl.nasa.gov/).
Figure 1. HEALPix: The resolution increase by three steps from the base level (the figure is cited from https://healpix.jpl.nasa.gov/).
Sensors 19 00427 g001
Figure 2. CAD models for HEALPix-IA tests.
Figure 2. CAD models for HEALPix-IA tests.
Sensors 19 00427 g002
Figure 3. The results of different EGI-based methods for HEALPix-IA tests.
Figure 3. The results of different EGI-based methods for HEALPix-IA tests.
Sensors 19 00427 g003
Figure 4. RMS results in logarithm of four rough registration methods contrast tests.
Figure 4. RMS results in logarithm of four rough registration methods contrast tests.
Sensors 19 00427 g004
Figure 5. Running time results of four rough registration methods contrast tests.
Figure 5. Running time results of four rough registration methods contrast tests.
Sensors 19 00427 g005
Figure 6. Other types of models for HEALPix-IA tests.
Figure 6. Other types of models for HEALPix-IA tests.
Sensors 19 00427 g006
Figure 7. Various sensors used for tests.
Figure 7. Various sensors used for tests.
Sensors 19 00427 g007
Figure 8. Results of different types of models for HEALPix-IA tests.
Figure 8. Results of different types of models for HEALPix-IA tests.
Sensors 19 00427 g008
Figure 9. Results of real workpiece MA analysis using HEALPix-IA.
Figure 9. Results of real workpiece MA analysis using HEALPix-IA.
Sensors 19 00427 g009
Table 1. Features of nine data for tests.
Table 1. Features of nine data for tests.
Model (Abbreviation)FeaturesSensors
Water Turbine Blade (WaTB)complex Curve Surface3D laser scanner
Wind Turbine Blade (WiTB)complex curve Surface3D laser scanner
Casting Parts 1 (CP1)thin-walled, large plane3D laser scanner
Casting Parts 2 (CP2)mirror symmetry3D laser scanner
Carbon Fiber Seat (CFS)thin-walled, curve surface3D laser scanner
Human Arm Dummy (HAD)unstructureed features3D laser scanner
Laser Cut Rocket (LCR)rotational symmetry3D laser scanner
Star Destroyer (SD)mirror symmetry3D Camera
3D Map In Door (Map)low dense data3D LiDAR
Table 2. MA analysis with different rough registration methods.
Table 2. MA analysis with different rough registration methods.
Models and Criteria Sac-IAHEALPix-IA
WaTB-RMSAve2.9 × 10−23.4 × 10−3
Var8.1 × 10−44.9 × 10−5
WaTB-time (ms)Ave4932.13830.8
Var6.3 × 1032.7 × 104
WaTB-opt time (ms)Ave7222.710,110.5
Var1.6 × 1072.2 × 107
CP2-RMSAve1.2 × 10−27.2 × 10−3
Var3.5 × 10−51.1 × 10−6
CP2-time (ms)Ave5953.24818.3
Var8.7 × 1034.5 × 103
CP2-opt time (ms)Ave9049.98066.5
Var2.65 × 1071.4 × 107

Share and Cite

MDPI and ACS Style

Gao, Y.; Du, Z.; Xu, W.; Li, M.; Dong, W. HEALPix-IA: A Global Registration Algorithm for Initial Alignment. Sensors 2019, 19, 427. https://doi.org/10.3390/s19020427

AMA Style

Gao Y, Du Z, Xu W, Li M, Dong W. HEALPix-IA: A Global Registration Algorithm for Initial Alignment. Sensors. 2019; 19(2):427. https://doi.org/10.3390/s19020427

Chicago/Turabian Style

Gao, Yongzhuo, Zhijiang Du, Wei Xu, Mingyang Li, and Wei Dong. 2019. "HEALPix-IA: A Global Registration Algorithm for Initial Alignment" Sensors 19, no. 2: 427. https://doi.org/10.3390/s19020427

APA Style

Gao, Y., Du, Z., Xu, W., Li, M., & Dong, W. (2019). HEALPix-IA: A Global Registration Algorithm for Initial Alignment. Sensors, 19(2), 427. https://doi.org/10.3390/s19020427

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