Next Article in Journal
Experimental Demonstration of Single-Channel EEG Signal Using 32 × 32 Pixel OLED Screen and Camera
Next Article in Special Issue
Multivariate Temporal Convolutional Network: A Deep Neural Networks Approach for Multivariate Time Series Forecasting
Previous Article in Journal
A Minimum Mean-Square-Error (MMSE) Decoder for Quasi-Orthogonal Space–Time Block Code
Previous Article in Special Issue
Modeling and Analysis of Adaptive Temperature Compensation for Humidity Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization and Correction of the Geometric Errors in Using Confocal Microscope for Extended Topography Measurement. Part I: Models, Algorithms Development and Validation

1
Department of Mechanical Engineering, Chemical and Industrial Design, ETS of Engineering and Industrial Design, Technical University of Madrid, 28012 Madrid, Spain
2
Department of Precision Mechanical Engineering, Shanghai University. No.333, Nanchen Rd., Shanghai 200444, China
*
Author to whom correspondence should be addressed.
Submission received: 2 June 2019 / Revised: 25 June 2019 / Accepted: 26 June 2019 / Published: 27 June 2019

Abstract

:
This work presents a method for characterizing and correcting the geometric errors of the movement of the lateral stage of Imaging Confocal Microscope (CM) in extended topography measurement. For an extended topography measurement, a defined number of 2D images are taken and stitched by correlation methods. Inaccuracies due to linear displacement, vertical and horizontal straightness errors, angular errors, and squareness errors based on the assumption of the rigid body kinematics are described. A mathematical model for the scale calibration of the X- and Y- coordinates is derived according to the system kinematics, the axis chain vector of CM, and the geometric error functions and their approximations by Legendre polynomials. The correction coefficients of the kinematic modelling are determined by the measured and certified data of a dot grid target standard artefact. To process the measurement data, algorithms for data partitions, fittings of cylinder centers, and determinations of coefficients are developed and validated. During which methods such as form removal, K-means clustering, linear and non-linear Least Squares are implemented. Results of the correction coefficients are presented in Part II based on the experimental studies. The mean residual reduces 29.6% after the correction of the lateral stage errors.

1. Introduction

The calibration of measuring machines is both important for test acceptance and error compensation [1]. Use of software techniques started from the very beginning for the correction of the systematic errors of measurement instruments, as the mechanical accuracy is expensive while repeatability cost little [2]. Over the last several decades, measurement accuracy and error compensation have been an area of intensive investigation [3,4,5].
Lateral calibration of the X- and Y- coordinates serves as a calibration of the magnification of the X- and Y-axis scale [6]. Lateral calibration/correction of confocal microscope (CM) can be classified mainly into two groups: imaging system calibration and machine system calibration [7]. Calibrations of the imaging system studies aberrations of refractive systems of objective lenses [8,9,10], axial distortion by point spread function or refractive-index mismatches [11,12], etc. The machine system calibrations investigate geometric errors generated by the movement of probes [2,13,14]. The first study of using kinematic geometric errors for error compensation of coordinate measurement traces back to the work of G Zhang et al. [15]. From then on, characterization and correction of kinematic geometric errors has mainly focused on coordinate measuring machines (CMM) [16,17], and has seldom been applied for optical measurement instruments.
However, areal measurement by CM is realized by both imaging systems and machine systems. CM scan in the lateral and vertical directions using different physical mechanisms [18]. To obtain a surface topography, CM captures a series of two-dimensional images by stepping either the specimen holding base or the objectives along the Z-stage [19]. For a single topography measurement, the 2D images are acquired in lateral directions either all at once by a CCD array (Imaging CM) or by a raster scan (Laser Scanning CM) [20,21]. For an extended topography measurement, a defined number of 2D images are taken and stitched by correlation methods. Accordingly, calibrating the imaging system only is not sufficient for a CM, as its measurement is realized by both imaging and machine systems. Lateral calibration or machine system calibration is necessary and important. It is worth mentioning that some works have been carried out on characterization of metrological characteristics defined in ISO 25178-600 [22], such as amplification coefficient, linearity deviation, x-y perpendicularity deviation, etc. [23,24,25,26]. These metrological characteristics imply the whole instrumentation system’s characteristics to some extent.
In this work, we present a method to characterize and correct the geometric errors of the lateral stage of the CM, with implementation of a standard artefact of dot grid target. Extended topography measurements are influenced by kinematic errors of the measurand holding base, which moves horizontally in the X- and Y-directions. This work first introduces the theory and mathematics of the 21 rigid body geometric errors applied generally in the study of error compensations for CMM and machine tools in Section 2 Afterwards, we develop our own kinematic modelling based on the theory of kinematic and geometric errors for the lateral stage movement of CM in the same Section. Section 3 introduces the methodology of the study. Mathematical models and algorithms for measurement data processing are presented in Section 4. Section 5 validates our developed models and algorithms using synthetic data. Section 6 draws the conclusions. Experimental studies of the calibration and uncertainty evaluation are presented in Part II, published as another paper in the same journal.

2. Mathematical Model for the X- and Y-Scale Calibration

For each of the three axes there exist three translational deviations, i.e. the linearity deviation, the straightness deviation and the ortoganality deviation. The straightness deviation is defined as the displacement orthogonal to the axis of motion, with two straightness deviation functions for each of the orthogonal directions. Rotational movements of stages are roll, yaw and pitch, with roll denoting a screw like movement, and pitch and yaw denoting angular deviation functions that describe movements within a plane [1,15,27]. Some parameters of mechanical systems, such as bearing spacing and guideway geometric errors, determine axis motion errors of machine axes [28]. T. O. Ekinci et al. categorized machine errors into three levels, i.e., geometric errors, joint kinematic errors, and volumetric errors [29]. This work clarifies the terminologies and highlights the usefulness of a machine error modeling approach.
For the measurement of a specimen using the 3D Imaging CM, the workpiece is placed on the X- and Y-stage above the supporting base, which is usually made of granite. As shown in Figure 1a, to focus the laser beam onto the measurand surface, two movements are carried ou, i.e., movement along the X- and Y-axis, controlled by the lateral stage, and movement along the Z-axis, controlled by the vertical stage. For example, to measure the point P on the measurand surface, its position can be denoted as P ( P x , P y , P z ) with respect to its coordinate origin O ( 0 , 0 , 0 ) . For the focus of the probe and the point of workpiece, the first movement include translation of X carriage followed by the translation of Y carriage; the second movement is the translation of Z carriage. It can be suggested that vectors X , Y and Z represent the translations of X, Y and Z carriages, the rotation matrices R ( X ) , R ( Y ) and R ( Z ) denote the angular motion errors caused by the translations of X, Y and Z carriages, the vector T represents the X, Y and Z abbe offsets of the probe with respect to the carriage to which the probe is attached. For our Imaging CM, the probe is attached to the Z-stage. Therefore, it is obvious that the translation of the X and Y carriage, with an offset of the position P , stops at the same point reached by the translation of the Z carriage with a probe offset T , which can be understood as the focus of the probe with the measurand. According to this conclusion, a diagram of the axis chain vector of 3D imaging microscopy is plotted and shown in Figure 1b.
As shown in Figure 1b, there are two kinematical paths from the coordinate origin O ( 0 , 0 , 0 ) to the laser beam focus, i.e., Z T and X Y P . In each kinematical path, the actual movement of that axis is affected by the rotational error of its predecessors. The movement of T is affected by the rotational error of the movement of Z . Y is affected by the rotational error of X . P is affected by the rotational errors of X and Y . Therefore, the chain vector shown in Figure 1b for the two kinematical paths, from the same origin reaching the same laser beam focus, can be expressed by Equation (1).
Z + R - 1 ( Z ) T = X + R - 1 ( X ) Y + R - 1 ( X ) R - 1 ( Y ) P
After rearrangement of the above equation, vector P can be expressed by all the other vectors and matrices, as the same result presented by G. Zhang et al. [1].
The rotation can be expressed by the infinitesimal rotation matrix. The vectors X , Y , Z and T can be substituted by their position matrix. This has been widely accepted and implemented by many former investigators [1,15,27,30].
In the 3D Imaging CM, there are no real X- and Y-carriages, the workpiece is carried by a lateral stage mounted on a fixed base [13,27], the Z-stage is mounted on another fixed base. In addition to that, the beam focus is considered to be the probe. As the focus is always on the workpiece surface and the information is acquired exactly at that point, there is no Abbe offset on x , y and z , nor is there an angular term. This means T = 0 and R ( Z ) = 0 .
In this work, only geometric errors along the X- and Y-directions are studied. After dropping out the Z-stage variables, the geometric error functions for X- and Y-stage are:
P x = x + δ x ( x ) + δ x ( y ) y σ z ( y )
P y = y + δ y ( x ) + x σ z ( x ) + δ y ( y ) + x σ z ( y )
Those errors in Equations (2) and (3) are functions, which can be simulated by Legendre polynomials. Implementation of Legendre polynomials is computationally simpler and provides a reduction of around 2% higher than using Chebyshev polynomials [31]. After substituting the error functions in Equations (21) and (22) by using the approximation of Legendre polynomials, the final mathematical models for the beam focus point P ( P x , P y , 0 ) with error corrections can be acquired:
P x = x + a 1 x + 3 2 a 2 x 2 + 3 2 a 3 ( 5 x 3 3 x ) + + b 1 y + 3 2 b 2 y 2 + 3 2 b 3 ( 5 y 3 3 y ) y ( c 1 y + 3 2 c 2 y 2 + 3 2 c 3 ( 5 y 3 3 y ) )
P y = y + d 1 x + 3 2 d 2 x 2 + 3 2 d 3 ( 5 x 3 3 x ) + + x ( e 1 x + 3 2 e 2 x 2 + 3 2 e 3 ( 5 x 3 3 x ) ) + + f 1 y + 3 2 f 2 y 2 + 3 2 f 3 ( 5 y 3 3 y ) + + x ( c 1 y + 3 2 c 2 y 2 + 3 2 c 3 ( 5 y 3 3 y ) )

3. Methodology and a Brief Introduction of Experimental Design

To calibrate the kinematic geometric errors of the lateral stage of a dual core 3D CM, a standard artefact with dot grid target on glass is measured and a series of algorithms are developed for the process of the raw measurement data, which include separations of cylinders and flat, fitting of cylinder centers, determination of coefficients, etc. Figure 2 illustrates the standard, which was provided by Max Levy Autograph, Inc. Those dots are deposited on the substrate, which is glass with a thickness of 1.5 mm. Table 1 gives an indication of the size, the dot spacing, the diameter, the accuracy, etc., of the dot pattern.
This work aims at characterizing and correcting the kinematic geometric errors of the movement of the lateral stage, by calculating and applying the coefficients of Equations (4) and (5). To calculate the coefficients, the introduced dot grid target standard is measured, the measurement data are processed for obtaining the coordinate values ( x i , y i ) of the dots’ centers, the coordinate values ( x i , y i ) and their corresponding certified values are put into Equations (4) and (5), the coefficients are obtained by solving the equations using non-linear Least Squares method. The obtained coefficients and Equations (4) and (5) are implemented for correcting the new measurement data. Another different area of the introduced dot grid target standard is measured, and the measurement data is corrected. By comparing the residuals before correction and after correction, the significance and practicality of our model for kinematic geometric error correction can be observed.
The algorithms used for measurement data processing and for calculating the coefficients are introduced in Section 4. Validations of those algorithms are introduced in Section 5. The details of the experimental procedures, parameters, and analysis are presented in Part II, which is published in another paper in the same journal.

4. Algorithms and Procedures for Measurement Data Processing

4.1. Algorithm and Procedure for the Separation of Flats and Cylinders

The raw measurement data usually contains defects, outliers, unmeasured points, etc. [32]. The measured surface often inclines, as it is impossible to locate the measurand surface perfectly perpendicular to the Z-axis. Figure 3 gives an example of the measurement of the dot grid target, implementing an objective of magnification 50×, numerical aperture of 0.90. The acquisition parameter of measurement area was defined as the topography stitching measurement, with 8 × 8 extended topographies, covering an area of 1.78 × 1.33 mm2. The parameter of the overlapping area is 25%, and the correlation takes the XYZ option. The level of resolution is 2, and the measured extended topographies contain 2673 × 2003 pixels.
This measurement has many unmeasured points in the right-top corner due to the Z-range parameter setting. The first step is to choose an area of this raw surface with appropriate size and location, targeting at minimizing the influence of defects and removing outliers [32] by selecting valid grid positions with loss of the uniformity of the grid. Figure 4 shows a surface reconstruction by trimming the raw data with 10% edges along the X- and Y-directions. As there are too many measurement defects when y > 800   μ m , this part is also abandoned.
After obtaining an appropriate area for analysis, it is necessary to eliminate its form. The principal axes of the distribution of data points of the measured point cloud are determined by solving the eigen value problem. The data are rotated by principal axis transformation to align them in the direction of the eigen vectors. Figure 5 gives an example of the rotated surface, which is rotated from the surface shown in Figure 4. Figure 5a was plotted by the Matlab® function ‘plot3′, while Figure 5b was plotted by the function ‘surf’. It can be found that after rotation, the range of the Z-axis is much smaller, and many details of the surface can be observed. It is obvious that this surface is constructed using point clouds, which form a flat plane and many cylinders perpendicular to this flat.
This rotated surface is now ready for the separation of flats and cylinders. Figure 6 gives an indication of the distribution of the heights of the surface. Figure 7 and Figure 8 give examples of the separations of the rotated surface shown in Figure 5.
After separating the cylinders from the flat, the next step is to obtain the three coordinates of each cylinder point cloud. Figure 9 gives an example of the separation of cylinders into individual clusters, characterizing the cylinders with different colors. Here the squared Euclidean distance metric is applied for distance calculation, as shown in Equation (6). Moreover, each centroid of the point cloud is the mean of the coordinates’ values of the points.
d ( x , c ) = ( x c ) ( x c ) T
where x is the initial observations, c is the centroid.
According to the previous investigations, there are several efficient algorithms for data partitioning, such as the K-means clustering algorithm [33,34] and Lloyd’s algorithm [35]. These two algorithms are iterative data-partitioning algorithms, which assign the initial observations into k clusters according to the squared Euclidean distance metric and the centroids. The number of clusters k is defined before data partition. In our work, the Lloyd’s algorithm is used for finding the centroids of k clusters while the K-means clustering algorithm is used for partitioning the initial observations into k clusters according to the squared Euclidean distance metric with respect to the centroid of each cluster. The K-means clustering algorithm requires the input as a matrix of M points in N dimensions as well as a matrix of K initial cluster centers in N dimensions. Denoting the number of points in cluster L is N C ( L ) , and the Euclidean distance between point I and cluster L is D ( I , L ) , the general procedure of the K-means clustering algorithm is to find a K-partition of clusters with locally optimal within-cluster sum of squares by changing the initial observations from one cluster to another [33], as indicated by Equation (7).
min D 2 i = 1 K j = 1 N C ( L ) D 2 ( I , L )
The procedures of the algorithm for partitioning the filtered cylinders into individuals are shown in Figure 10, which consists of 6 steps. It is worth mentioning that the data are subjected to dimensionality reduction before data partitioning. These cylinders are a point cloud of three dimensions. The aim of dimensionality reduction is to produce a compact low-dimensional encoding of a given high-dimensional data set [36]. The dimensionality can be reduced in two ways [37]. The first is by only keeping the most relevant variables from the original dataset, which is called feature selection. The second way is to exploit the redundancy of the input data and to find a smaller set of new variables, each being a combination of the input variables containing basically the same information as the input variables. In this work, the dimensionality of the data is reduced to 2D from 3D by projection along Z-axis.
Since the equations characterizing lateral linearity deviations, Equations (4) and (5) are independent of the vertical axis z, the dimension of the 3D point cloud is reduced to the two dimensions of the x-y-plane.

4.2. Algorithm for Determination of Coefficients

This section aims at the determination of the coefficients, a 1 , a 2 , a 3 , b 1 …, f 3 , defined in Equations (4) and (5). The objective is to fit the 18 coefficients to obtain the least squared residuals between the corrected X- and Y- coordinates and the nominal position. That is:
i = 1 n P n o m i ( x , y ) P i c o r r ( x , y ) 2 2 = min
where, P n o m i ( x , y ) represents the i t h nominal position defined by the X- and Y- coordinates, P c o r r i ( x , y ) denotes the i t h corrected position obtained from the measured data and Equations (4) and (5).
Therefore, the objective function is
min F ( x , y ) | a 1 , a 2 , , g 3 = m i n i = 1 n P n o m i ( x , y ) P i c o r r ( x , y ) 2 2 = min ( ( x n o m i x c o r r i ) 2 + ( y n o m i y c o r r i ) 2 ) = min ( ( x n o m i x m e a s i + a 1 x m e a s i + 3 2 a 2 ( x m e a s i ) 2 + 3 2 a 3 ( 5 ( x m e a s i ) 3 3 x m e a s i ) + + b 1 y m e a s i + 3 2 b 2 ( y m e a s i ) 2 + 3 2 b 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) y m e a s i ( c 1 y m e a s i + 3 2 c 2 ( y m e a s i ) 2 + 3 2 c 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) ) ) 2 + + ( y n o m i y m e a s i + d 1 x m e a s i + 3 2 d 2 ( x m e a s i ) 2 + 3 2 d 3 ( 5 ( x m e a s i ) 3 3 x m e a s i ) + + x m e a s i ( e 1 x m e a s i + 3 2 e 2 ( x m e a s i ) 2 + 3 2 e 3 ( 5 ( x m e a s i ) 3 3 x m e a s i ) ) + + f 1 y m e a s i + 3 2 f 2 ( y m e a s i ) 2 + 3 2 f 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) + + x m e a s i ( c 1 y m e a s i + 3 2 c 2 ( y m e a s i ) 2 + 3 2 c 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) ) ) 2 )
The minimum value of F arrives when the gradient is zero [38,39], i.e.,:
F β i = 0
Where F is the objective function, β i represents the parameters of coefficients to be determined. There are 18 coefficients in this part of the work, and hence there are 18 partial derivatives:
F a 1 = i = 1 n 2 ( x n o m i x m e a s i + a 1 x m e a s i + 3 2 a 2 ( x m e a s i ) 2 + 3 2 a 3 ( 5 ( x m e a s i ) 3 3 x m e a s i ) + + b 1 y m e a s i + 3 2 b 2 ( y m e a s i ) 2 + 3 2 b 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) y m e a s i ( c 1 y m e a s i + 3 2 c 2 ( y m e a s i ) 2 + 3 2 c 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) ) ) x m e a s i = 0
F a 2 = i = 1 n 2 ( x n o m i x m e a s i + a 1 x m e a s i + 3 2 a 2 ( x m e a s i ) 2 + 3 2 a 3 ( 5 ( x m e a s i ) 3 3 x m e a s i ) + + b 1 y m e a s i + 3 2 b 2 ( y m e a s i ) 2 + 3 2 b 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) y m e a s i ( c 1 y m e a s i + 3 2 c 2 ( y m e a s i ) 2 + 3 2 c 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) ) ) 3 2 ( x m e a s i ) 2 = 0
F f 3 = i = 1 n 2 ( y n o m i y m e a s i + d 1 x m e a s i + 3 2 d 2 ( x m e a s i ) 2 + 3 2 d 3 ( 5 ( x m e a s i ) 3 3 x m e a s i ) + + x m e a s i ( e 1 x m e a s i + 3 2 e 2 ( x m e a s i ) 2 + 3 2 e 3 ( 5 ( x m e a s i ) 3 3 x m e a s i ) ) + + f 1 y m e a s i + 3 2 f 2 ( y m e a s i ) 2 + 3 2 f 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) + + x m e a s i ( c 1 y m e a s i + 3 2 c 2 ( y m e a s i ) 2 + 3 2 c 3 ( 5 ( y m e a s i ) 3 3 y m e a s i ) ) ) 3 2 ( 5 ( y m e a s i ) 3 3 y m e a s i ) = 0
These equations can first be solved by solving the equation for the initial solutions. The initial solutions might be far from the correct solutions, as there are too many more constraints than unknown coefficients. Those initial solutions can be used as the starting point for iteratively finding the nonlinear least squares solutions. Extensive work has been done on the nonlinear least squares algorithms [40,41]. Let the model for data fitting be [41]:
E ( y ) = f ( x 1 , x 2 , , x m ; b 1 , b 2 , , b k ) = f ( x , b )
where x 1 , x 2 , , x m are independent variables, b 1 , b 2 , , b k are k parameters of coefficients to be determined, E ( y ) is the expected value of the dependent variable y . Denote the data points as: ( Y i , X 1 i , X 2 i , X m i ) , where, i = 1 , 2 , , n . The objective is to calculate the k parameters of coefficients which will minimize the squares of the residuals:
min i = 1 n ( Y i Y ^ i ) 2 = min Y - Y ^ 2
This problem can be written as an objective function which aims at optimizing the coefficients of each function Y i Y ^ i , denoting the vector of the parameters to be optimized as t , where t = [ a 1 , a 2 , , f 3 ] in this work. The objective function is:
min F ( t ) = i = 1 k f i 2 ( t )
where, f i ( t ) are functions of the to-be-optimized vector t . When f i ( t ) are nonlinear functions with respect to t , this problem is about nonlinear optimization. The solution is to use Taylor expansion to convert f i ( t ) into linear functions. As only the first order of the Taylor series is linear, here we approximate it by the first order:
φ i ( t ) = f i ( t ( k ) ) + f i ( t ( k ) ) T ( t t ( k ) ) = f i ( t ( k ) ) T t f i ( t ( k ) ) T t ( k ) + f i ( t ( k ) )
where f i ( t ( k ) ) is the value of the first derivative of f i ( t ( k ) ) on vector t evaluated at the point t ( k ) .
Substituting f i ( t ) in Equation (16) with its Taylor Expansion approximation indicated by Equation (17), the approximation of F ( t ) is:
min ϕ ( t ) = i = 1 k φ i 2 ( t ) = i = 1 k ( f i ( t ( k ) ) T t f i ( t ( k ) ) T t ( k ) + f i ( t ( k ) ) ) 2
To simplify the above equation, use A k to represent f i ( t ( k ) ) T , and b to represent f i ( t ( k ) ) T t ( k ) f i ( t ( k ) ) . Equation (18) is simplified into:
min ϕ ( t ) = i = 1 k ( A k b ) 2
where,
A k = [ f 1 ( t ( k ) ) T f m ( t ( k ) ) T ] = [ f 1 ( t ( k ) ) t 1 f 1 ( t ( k ) ) t 2 f 1 ( t ( k ) ) t n f m ( t ( k ) ) t 1 f m ( t ( k ) ) t 2 f m ( t ( k ) ) t n ]
b = [ f 1 ( t ( k ) ) T t ( k ) f 1 ( t ( k ) ) f m ( t ( k ) ) T t ( k ) f m ( t ( k ) ) ] = A k t ( k ) f ( k )
f ( k ) = [ f 1 ( t ( k ) ) f 2 ( t ( k ) ) f m ( t ( k ) ) ]
Therefore, it can be obtained:
ϕ ( t ) = ( A k t b ) T ( A k t b )
The solution of this equation is:
t ( k + 1 ) = t ( k ) ( A k T A k ) 1 A k T f ( k )
This solution can be simplified as:
t ( k + 1 ) = t ( k ) H k 1 F ( t ( k ) )
where, H k is the Hessian matrix:
H k = 2 A k T A k
F ( x ( k ) ) = 2 A k T f ( k )
Equation (25) can be solved iteratively by starting with an initial solution.

5. Validation of the Algorithms with Synthetic Data

5.1. Validation of the Algorithm for Determination of Dots’ Centers and Distance

To validate the algorithm for determination of dots’ centers, synthetic data with known cylinder centers and point cloud distributions is generated. This generated data contains 35 cylinders of point cloud, distributed as 5 × 7 along the X-axis and Y-axis, as shown in Figure 11.
For each point cloud, the radius of the cylinder is 31.25 μm, with a uniform distribution of r ~ U (−0.1, 0.1) μm. The values of the Z-coordinate of each point follows a standard normal distribution z ~ N (0, 1) μm. The centers of each cylinder, as indicated by red asterisks in Figure 11b, have an interval of 125 μm. Their X-coordinates are [0, 125, 250, 375, 500] μm and the Y-coordinates are [0, 125, 250, 375, 500, 625, 750] μm. The generated data contains X-, Y-, and Z-coordinates’ values. Those data are saved in a Matlab® file with the suffix name ‘.mat’ for processing.
The generated data are processed using our developed algorithm for the determination of dots’ centers. The first step is to separate the data into individual point clouds. The results of the point cloud separation are shown in Figure 12, with each individual being indicated by numbers and different colors. In this step, the algorithm not only separates the point cloud into individuals, but also gives an indication of the initial centroids of each cloud.
After the separation, each individual point cloud is processed by our developed algorithm in order to fit the circle centers. The calculated centers are exactly the same with the generated ones. Figure 13 makes a plot of the calculated centers of the synthetic data generated for algorithm validation.
The results indicate that the algorithm developed in this work for the determination of the centers of cylinder point cloud is valid.

5.2. Validation of the Algorithm for Determination of Coefficients

Here, we use two methods for the validation of the algorithm for the determination of the coefficients.
The first method is to create a dataset containing points with intervals of 125 μm as certified positions, while making a little displacement with each point. Calculate the correction coefficients and correct the points with the displacement using the algorithm defined in Section 4.2. The points with the displacement with respect to the certified points represent the measured points in the experiments. Therefore, two residuals can be obtained, i.e., the residual between the certified points and the points with displacement as well as the residual between the certified points and the corrected points. Compare those two residuals to check whether the correction is meaningful.
The algorithm first creates the two data described above. Then it rotates the points with displacement to the certified points, aiming at making the axis parallel. This arises from the methodology of this work, which requires locating the standard artefact that is parallel with the axis. As it is impossible to realize it manually, the coordinate axes are adjusted mathematically. After that, the error coefficients are calculated by the algorithm. With those coefficients, the values of the X- and Y-coordinates of the corrected points can be obtained. Figure 14 shows a comparison of the positions of the certified, measured and corrected points. The contour of the errors between certified and measured points is displayed in Figure 15. The contour of the errors between the certified and corrected points is shown in Figure 16. It is obvious that almost all of the errors are smaller after correction.
The error is calculated by the Euclidean distance between the certified point and its corresponding point with displacement. The same goes with the error between the certified point and the corrected point. Table 2 makes a comparison of the error before correction and after correction. The difference is calculated as:
D i f f ( i ) = E r r c ( i ) E r r m ( i )
where, D i f f ( i ) is the difference between the error after correction and before correction. When D i f f ( i ) is negative, it means the error is smaller after correction. E r r c ( i ) is the Euclidean distance between the certified point and corrected point. E r r m ( i ) is the Euclidean distance between the certified point and the measured point. The mean residual after correction is 5.60 μm, while the mean error between the measured and the certified point is 9.65 μm. The square residual of all the 35 corrected points is 1587.99 μm2, while it was 3684.85 μm2 before correction. It can be found that some points are more distorted after correction. This is because the fitting method used is the nonlinear least squares, the objective function of which aims at finding the sum of least squares of the residuals of all points, but does not ensure that every residual is smaller after fitting.
The second method creates another 10 data by varying the displacement from the certified positions based on this data. In these 10 simulations, the certified positions are all the same with those certified positions introduced above. The measured points’ positions are assigned with displacements, the values of which are shown in Table 3.
The 10 data are processed by this algorithm. The results of each simulation are compared and analyzed. Figure 17 indicates the positions of the certified, measured, and corrected points of those 10 simulations. It is obvious that almost all of the corrected points are closer than the measured points to the certified points.
Table 4 lists two mean residual and two squared residuals. One of the mean residuals arises from the displacements of the measured points from the certified points, while the other mean residual comes from the corrected points from the certified points. The mean residual is calculated by Equations (29) and (30):
Re s m e a s M = 1 n i = 1 n P c e r t ( i ) P m e a s ( i ) 2
Re s c o r r M = 1 n i = 1 n P c e r t ( i ) P c o r r ( i ) 2
where, Re s m e a s M and Re s c o r r M represent the mean residual of the measured points and the mean residual of the corrected points, respectively, P c e r t ( i ) represents the position of the i t h point of the certified positions, P c o r r ( i ) represents the position of the i t h point of the corrected positions.
The two squared residuals include the squared residual of the measured points to the certified points and the squared residual of the corrected points to the certified points. The squared residual is calculated by:
Re s m e a s S = i = 1 n ( P c e r t ( i ) P m e a s ( i ) 2 ) 2
Re s c o r r S = i = 1 n ( P c e r t ( i ) P c o r r ( i ) 2 ) 2
Figure 18 and Figure 19 plot the mean residuals and squared residuals of the 10 simulations. It can be observed that both the mean residual curve and the squared residual curve of the measured points are concave in the middle and convex at the two sides. The two residuals of the corrected points are much smoother. The larger the original residuals, the more the correction can be realized.
The results of the Simulation No. 10 are shown. Figure 20 compares the positions of the certified, measured, and corrected points of Simulation No. 10. Figure 21 and Figure 22 presents contours of the mean residual and the squared residual between certified and measured points of Simulation No. 10. Figure 23 shows the error vectors between the measured and corrected points of Simulation No. 10.
As the results above show, all the measured values are improved after geometric error correction. It can be concluded that our developed algorithm for the determination of correction coefficients is valid.

6. Conclusions

The work presented herein is Part I of the whole work, describing our methods for characterization and correction of lateral stage geometric errors. With implementation of the 21 parametric errors based on the assumption of the rigid body kinematics for calculating the machine volumetric errors of a three-axis machine, a mathematical model for correcting the lateral stage deviations of the CM was developed.
According to the measurement principle of CM, a diagram of the axis chain vector was drawn. Equation (1) was developed according to the kinematic path of the chain vector. After the rearrangement of Equation (1) and its simplification by dropping Z and T , this equation was further simplified due to the study specifications and characteristics of the CM. The function errors were approximated by Legendre polynomials. Finally, Equations (4) and (5) for corrections of kinematic geometric errors were obtained. The methodology for calculation and application of the corrections coefficients defined in Equations (4) and (5) are introduced afterwards. Experiments on a dot grid target standard artefact will be measured and studied in Part II, published in the same journal.
The algorithms that will be implemented for the measurement data processing are presented and validated. These algorithms apply a lot of methods, including form removal, K-means clustering, linear and non-linear Least Squares. Their mathematical models were also introduced. The results of the validation based on synthetic data imply that our mathematical models and algorithms are reliable.

Author Contributions

Conceptualization and Methodology, C.W. and E.G.; Algorithms, Experiments, Data, Manuscript writing, C.W.; Resources, Funding, Supervision, E.G. and Y.Y.; Review and Editing, C.W., E.G. and Y.Y.

Funding

This work is funded by the Spanish State Programme of Promotion of Scientific Research and Technique of Excellence, State Sub-programme of Generation of Knowledge. Project DPI2016-78476-P “Desarrollo Colaborativo de Patrones de Software y Estudios de Trazabilidad e Intercomparación en la Caracterización Metrológica de Superficies”, belonging to the 2016 call for R & D Projects. The authors acknowledge the Chinese Scholarship Council (CSC) for funding the first author’s doctoral study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CMconfocal microscope
CCDcharge-coupled device
CMMcoordinate measuring machine

References

  1. Zhang, G.; Ouyang, R.; Lu, B.; Hocken, R.; Veale, R.; Donmez, A. A Displacement Method for Machine Geometry Calibration. CIRP Ann. 1988, 37, 515–518. [Google Scholar] [CrossRef]
  2. Hermann, G. Geometric Error Correction in Coordinate Measurement. Acta Polytech. Hung. 2007, 4, 47–62. [Google Scholar]
  3. Krolczyk, G.M.; Maruda, R.W.; Nieslony, P.; Wieczorowski, M. Surface morphology analysis of Duplex Stainless Steel (DSS) in Clean Production using the Power Spectral Density. Measurement 2016, 94, 464–470. [Google Scholar] [CrossRef]
  4. Krolczyk, G.M.; Krolczyk, J.B.; Maruda, R.W.; Legutko, S.; Tomaszewski, M. Metrological changes in surface morphology of high-strength steels in manufacturing processes. Measurement 2016, 88, 176–185. [Google Scholar] [CrossRef]
  5. Wang, C.; D’Amato, R.; Gómez, E. Confidence Distance Matrix for outlier identification: A new method to improve the characterizations of surfaces measured by confocal microscopy. Measurement 2019, 137, 484–500. [Google Scholar] [CrossRef]
  6. Harding, K. (Ed.) Handbook of Optical Dimensional Metrology; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2013. [Google Scholar]
  7. Corle, T.R.; Kino, G.S. Confocal Scanning Optical Microscopy and Related Imaging Systems; Academic Press: Cambridge, MA, USA, 1996. [Google Scholar]
  8. Thomas, D.; Sugimoto, A. Parametric surface representation with bump image for dense 3d modeling using an rbg-d camera. Int. J. Comput. Vis. 2017, 123, 206–225. [Google Scholar] [CrossRef]
  9. Bailey, D.G. A new approach to lens distortion correction. Proc. Image Vis. Comput. N. Z. 2002, 2002, 59–64. [Google Scholar]
  10. Wang, J.; Shi, F.; Zhang, J.; Liu, Y. A new calibration model of camera lens distortion. Pattern Recognit. 2007, 41, 607–615. [Google Scholar] [CrossRef]
  11. Besseling, T.H.; Jose, J.; Van Blaaderen, A. Methods to calibrate and scale axial distances in confocal microscopy as as function of refractive index. J. Microsc. 2015, 257, 142–150. [Google Scholar] [CrossRef]
  12. Cole, R.W.; Jinadasa, T.; Brown, C.M. Measuring and interpreting point spread functions to determine confocal microscope resolution and ensure quality control. Nat. Protoc. 2011, 6, 1929–1941. [Google Scholar] [CrossRef]
  13. Wang, C. Current issues on 3D volumetric positioning accuracy: Measurement, Compensation and Definition. Proc. SPIE 2008, 7128. [Google Scholar] [CrossRef]
  14. Barakat, N.A.; Elbestawi, M.A.; Spence, A.D. Kinematic and geometric error compensation of a coordinate measuring machine. Int. J. Mach. Tools Manuf. 2000, 40, 833–850. [Google Scholar] [CrossRef]
  15. Zhang, G.; Veale, R.; Charlton, T.; Borchardt, B.; Hocken, R. Error Compensation of Coordinate Measuring Machines. CLRP Ann. 1985, 34, 445–448. [Google Scholar] [CrossRef]
  16. Schwenke, H.; Knapp, W.; Haitjema, H.; Weckenmann, A.; Schmitte, R.; Delbressine, F. Geometric error measurement and compensation of machines—An update. CIRP Ann. 2008, 57, 660–675. [Google Scholar] [CrossRef]
  17. Umetsu, K.; Furutnani, R.; Osawa, S.; Takatsuji, T.; Kurosawa, T. Geometric calibration of a coordinate measuring machine using a laser tracking system. Meas. Sci. Technol. 2005, 16, 2466–2472. [Google Scholar] [CrossRef]
  18. Semwogerere, D.; Weeks, E.R. Confocal Microscopy; Taylor & Francis: Abingdon, UK, 2005. [Google Scholar]
  19. Jensen, K.E.; Weitz, D.A.; Spaepen, F. Note: A three-dimensional calibration device for the confocal microscope. Rev. Sci. Instrum. 2013, 84, 016108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Senin, N.; Leach, R. Information-rich surface metrology. Procedia CIRP 2018, 75, 19–26. [Google Scholar] [CrossRef]
  21. Claxton, N.S.; Fellers, T.J.; Davidson, M.W. Laser Scanning Confocal Microscopy; Technical Report; Department of Optical Microscopy and Digital Imaging, Florida State University: Tallahassee, FL, USA, 2006. [Google Scholar]
  22. ISO 25178-600. Geometrical Product Specifications (GPS)—Surface Texture: Areal-Part 600: Metrological Characteristics for Areal-Topography Measuring Methods; International Organization for Standardization: Geneva, Switzerland, 2019. [Google Scholar]
  23. Leach, R.K.; Giusca, C.L.; Haitjema, H.; Evans, C.; Jiang, X. Calibration and verification of areal surface texture measuring instruments. CIRP Ann. 2015, 64, 797–813. [Google Scholar] [CrossRef]
  24. Alburayt, A.; Syam, W.P.; Leach, R. Lateral scale calibration for focus variation microscopy. Meas. Sci. Technol. 2018, 29, 065012. [Google Scholar] [CrossRef]
  25. Leach, R. Optical Measurement of Surface Topography; Springer: Berlin, Germany, 2011. [Google Scholar]
  26. Nouira, H.; El-Hayek, N.; Yuan, X.; Anwer, N.; Salgado, J. Metrological characterization of optical confocal sensors measurements (20 and 350 travel ranges). J. Phys. Conf. Ser. 2015, 483, 012015. [Google Scholar] [CrossRef]
  27. Tong, K.; Lehtihet, E.A.; Joshi, S. Parametric error modeling and software error compensation for rapid prototyping. Rapid Prototyp. J. 2003, 9, 301–313. [Google Scholar] [CrossRef]
  28. Ekinci, T.O.; Mayer, J.R.R.; Cloutier, G.M. Investigation of accuracy of aerostatic guideways. Int. J. Mach. Tools Manuf. 2009, 49, 478–487. [Google Scholar] [CrossRef]
  29. Ekinci, T.O.; Mayer, J.R.R. Relationships between straightness and angular kinematic errors in machines. Int. J. Mach. Tools Manuf. 2007, 47, 1997–2004. [Google Scholar] [CrossRef]
  30. Yang, J.; Ren, Y.; Wang, C.; Liotto, G. Theoretical derivations of 4 body diagonal displacement errors in 4 machine configurations. In Proceedings of the LAMDAMAP Conference, Cransfield, UK, 27–30 June 2005. [Google Scholar]
  31. Aguado, S.; Samper, D.; Santolaria, J.; Aguilar, J.J. Towards an effective identification strategy in volumetric error compensation of machine tools. Meas. Sci. Technol. 2012, 23, 065003. [Google Scholar] [CrossRef]
  32. Wang, C.; Caja, J.; Gomez, E. Comparison of methods for outlier identification in surface characterization. Measurement 2018, 117, 312–325. [Google Scholar] [CrossRef]
  33. Hartigan, J.A.; Wong, M.A. A K-means Clustering Algorithm. J. R. Stat. Soc. 1979, 28, 100–108. [Google Scholar]
  34. Jain, K. Data clustering: 50 years beyond K-means. Pattern Recognit. Lett. 2010, 31, 651–666. [Google Scholar] [CrossRef]
  35. Lloyd, S.P. Least Squares Quantization in PCM. IEEE Trans. Inf. Theory 1982, 28, 129–137. [Google Scholar] [CrossRef]
  36. Ghodsi, A. Dimensionality Reduction—A Short Tutorial; Technical Report; University of Waterloo: Waterloo, ON, Canada, 2006. [Google Scholar]
  37. Sorzano, S.; Vargas, J.; Montano, A.P. A survey of dimensionality reduction techniques, Technical Report. arXiv 2014, arXiv:1403.2877. [Google Scholar]
  38. Miller, S.J. The Method of Least Squares; Brown University: Providence, RI, USA, 2006. [Google Scholar]
  39. Weisstein, E.W. Least Squares Fitting, From MathWorld–A Wolfram Web Resource. Available online: http://mathworld.wolfram.com/LeastSquaresFitting.html (accessed on 1 October 2018).
  40. Dennis, J.E.; Gay, D.M.; Welsch, R.E. An adaptive nonlinear least-squares algorithm. ACM Trans. Math. Softw. 1981, 7, 348–368. [Google Scholar] [CrossRef]
  41. Marquardt, D.W. An algorithm for Least-Squares estimation of nonlinear parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
Figure 1. Manifestation of the vectors in 3D Imaging CM: (a) 3D view; (b) diagram of the axis chain vector.
Figure 1. Manifestation of the vectors in 3D Imaging CM: (a) 3D view; (b) diagram of the axis chain vector.
Electronics 08 00733 g001
Figure 2. Illustration of the standard artefact of the dot grid target.
Figure 2. Illustration of the standard artefact of the dot grid target.
Electronics 08 00733 g002
Figure 3. Reconstruction of the measured surface with its raw measurement data: (a) contour plot; (b) surf plot.
Figure 3. Reconstruction of the measured surface with its raw measurement data: (a) contour plot; (b) surf plot.
Electronics 08 00733 g003
Figure 4. Reconstruction of the surface with trimmed raw data: (a) contour plot; (b) areal plot.
Figure 4. Reconstruction of the surface with trimmed raw data: (a) contour plot; (b) areal plot.
Electronics 08 00733 g004
Figure 5. Reconstruction of the rotated surface: (a) point plot; (b) surfplot using Z for the color data.
Figure 5. Reconstruction of the rotated surface: (a) point plot; (b) surfplot using Z for the color data.
Electronics 08 00733 g005
Figure 6. Histogram of the Z values of the rotated surface.
Figure 6. Histogram of the Z values of the rotated surface.
Electronics 08 00733 g006
Figure 7. The separated flat plane of the rotated surface: (a) 3D view; (b) view from top down.
Figure 7. The separated flat plane of the rotated surface: (a) 3D view; (b) view from top down.
Electronics 08 00733 g007
Figure 8. The separated cylinders of the rotated surface: (a) 3D view; (b) view from top down.
Figure 8. The separated cylinders of the rotated surface: (a) 3D view; (b) view from top down.
Electronics 08 00733 g008
Figure 9. Separation of cylinders into individuals.
Figure 9. Separation of cylinders into individuals.
Electronics 08 00733 g009
Figure 10. Flowchart of the algorithm for data partition.
Figure 10. Flowchart of the algorithm for data partition.
Electronics 08 00733 g010
Figure 11. Plot of the created data for algorithm validation: (a) 3D view of the point cloud; (b) view from the Z-axis.
Figure 11. Plot of the created data for algorithm validation: (a) 3D view of the point cloud; (b) view from the Z-axis.
Electronics 08 00733 g011
Figure 12. The separation of cylinders into individual point clouds.
Figure 12. The separation of cylinders into individual point clouds.
Electronics 08 00733 g012
Figure 13. Plots of the calculated centers of the generated data.
Figure 13. Plots of the calculated centers of the generated data.
Electronics 08 00733 g013
Figure 14. Comparison of the positions of certified, measured and corrected points.
Figure 14. Comparison of the positions of certified, measured and corrected points.
Electronics 08 00733 g014
Figure 15. Contour of the errors between the certified and measured points.
Figure 15. Contour of the errors between the certified and measured points.
Electronics 08 00733 g015
Figure 16. Contour of the errors between the certified and corrected points.
Figure 16. Contour of the errors between the certified and corrected points.
Electronics 08 00733 g016
Figure 17. Comparison of the certified, measured, and corrected positions of the 10 simulations.
Figure 17. Comparison of the certified, measured, and corrected positions of the 10 simulations.
Electronics 08 00733 g017
Figure 18. Comparison of the mean residuals before correction and after correction.
Figure 18. Comparison of the mean residuals before correction and after correction.
Electronics 08 00733 g018
Figure 19. Comparison of the squared residuals before correction and after correction.
Figure 19. Comparison of the squared residuals before correction and after correction.
Electronics 08 00733 g019
Figure 20. Comparison of the positions of the certified, measured, and corrected points of Simulation No. 10.
Figure 20. Comparison of the positions of the certified, measured, and corrected points of Simulation No. 10.
Electronics 08 00733 g020
Figure 21. Contour of the mean residual between certified and measured points of Simulation No. 10.
Figure 21. Contour of the mean residual between certified and measured points of Simulation No. 10.
Electronics 08 00733 g021
Figure 22. Contour of the mean residual between certified and corrected points of Simulation No. 10.
Figure 22. Contour of the mean residual between certified and corrected points of Simulation No. 10.
Electronics 08 00733 g022
Figure 23. Comparison of the error vectors between the measured and corrected points of Simulation No. 10.
Figure 23. Comparison of the error vectors between the measured and corrected points of Simulation No. 10.
Electronics 08 00733 g023
Table 1. Detailed information about the dot grid target standard artefact.
Table 1. Detailed information about the dot grid target standard artefact.
Width [mm]Dot Diameter [mm]Dot Spacing [mm]X and Y Axis Accuracy [mm]Dot Array
XYTotal
250.06250.125±0.00120120140401
Table 2. The values of the error of each individual point.
Table 2. The values of the error of each individual point.
Point No. iErrc(i) [μm]Errm(i) [μm]Diff(i) [μm]Point No. iErrc(i) [μm]Errm(i) [μm]Diff(i) [μm]
1000194.667.65−2.99
22.924.35−1.43204.7212.64−7.93
34.879.20−4.31214.1613.09−8.94
46.3612.48−6.12228.2613.13−4.87
55.5513.25−7.702312.879.053.82
65.2713.13−7.862413.5213.92−0.41
72.639.64−7.01253.569.58−6.01
83.085.16−2.082614.1711.123.04
92.865.85−2.99273.397.24−3.86
105.3610.66−5.31281.618.03−6.43
114.6311.28−6.65297.4713.86−6.39
1211.9011.320.59303.7313.14−9.41
134.177.38−3.20312.888.69−5.81
143.8211.99−8.173215.9313.782.15
157.1910.38−3.20336.5914.88−8.28
160.623.85−3.23344.158.68−4.53
173.998.23−4.25354.493.820.67
184.727.36−2.64-
mean error [μm]5.609.65
Sum of Squared error [μm2]1587.993684.85
Table 3. Further displacements of the measured points of the 10 simulations.
Table 3. Further displacements of the measured points of the 10 simulations.
SimulationNo. 1No. 2No. 3No. 4No. 5No. 6No. 7No. 8No. 9No. 10
X [μm]−20−16−12−8−40481216
Y [μm]−15−12−9−6−3036912
Table 4. Mean residuals and squared residuals of the measured and corrected points to the certified points.
Table 4. Mean residuals and squared residuals of the measured and corrected points to the certified points.
No. SimulationMean Residuals [μm]Squared Residuals [μm2]
MeasuredCorrectedMeasuredCorrected
117.948.1112155.653287.68
213.138.456961.493269.41
38.568.133517.332986.80
45.815.231823.171206.33
56.425.961879.011288.29
69.655.603684.851587.99
713.936.407240.691903.75
818.538.2012546.533101.98
923.319.6619602.374179.53
1028.1811.0628408.215389.67

Share and Cite

MDPI and ACS Style

Wang, C.; Gómez, E.; Yu, Y. Characterization and Correction of the Geometric Errors in Using Confocal Microscope for Extended Topography Measurement. Part I: Models, Algorithms Development and Validation. Electronics 2019, 8, 733. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics8070733

AMA Style

Wang C, Gómez E, Yu Y. Characterization and Correction of the Geometric Errors in Using Confocal Microscope for Extended Topography Measurement. Part I: Models, Algorithms Development and Validation. Electronics. 2019; 8(7):733. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics8070733

Chicago/Turabian Style

Wang, Chen, Emilio Gómez, and Yingjie Yu. 2019. "Characterization and Correction of the Geometric Errors in Using Confocal Microscope for Extended Topography Measurement. Part I: Models, Algorithms Development and Validation" Electronics 8, no. 7: 733. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics8070733

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