Next Article in Journal
Land Use Change and Climate Variation in the Three Gorges Reservoir Catchment from 2000 to 2015 Based on the Google Earth Engine
Previous Article in Journal
Application of the TDR Soil Moisture Sensor for Terramechanical Research
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Sub-Regional Calibration Method That Can Accomplish Error Compensation for Photoelectric Scanning Measurement Network

State Key Laboratory of Precision Measuring Technology and Instruments, Tianjin University, Tianjin 300072, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(9), 2117; https://doi.org/10.3390/s19092117
Submission received: 15 April 2019 / Revised: 30 April 2019 / Accepted: 3 May 2019 / Published: 7 May 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
In the measurement process of photoelectric scanning measurement network, the laser surface edge area has lower measurement accuracy than the middle area due to the geometrical distortions of the laser surface of the transmitter. This paper presents a sub-regional calibration method that can accomplish error compensation for the measurement system. Unlike the camera sub-regional calibration, the regional division and identification of the laser surface are more difficult. In this paper, the pitch angle in the transmitter coordinate frame of the spatial point was used as the basis for the division and identification of the laser surface. In the calibration process, the laser surface of the transmitter was divided into different regions and each region was calibrated independently, so that an intrinsic parameters database containing the intrinsic parameters of different regions could be established. Based on the database, the region identification and error compensation algorithm were designed, and comparison experiments were carried out. With the novel calibration method, the measurement accuracy of the system had an obvious upgrade, especially at the edges of the laser surface within a certain measurement area, which could enlarge the effective measurement area of the transmitter and would broaden and deepen the application fields of photoelectric scanning measurement network.

1. Introduction

Large-scale metrology, as an indispensable role for modern industry, is routinely used in large-scale equipment manufacturing industry and precision engineering. The measurement systems for large-scale metrology can be classified by means of their system topology into centralized and distributed systems [1]. The centralized measurement system is a stand-alone unit that can independently measure spatial coordinates. The common centralized measurement systems include laser trackers [2], laser radars [3], tachymeters, laser scanners [4], and more. Laser scanners and laser radars are laser based range finders which enable a measurement of non-cooperative targets and featureless surfaces. For some laser scanners, the 3D coordinate measurement is based on triangulation. The distributed measurement system consists of several, separate and independent units whose separately gathered measurement information needs to be jointly processed in order for the system to determine the coordinates of a point. It is also possible to use centralized systems in combination as a distributed system, such as laser trackers. Compared with the centralized measurement system, the distributed system has advantages of measurement range scalability, which is very practical in large-scale metrology. There have been several distributed measurement systems provided for large-scale measurement, including theodolite networks [5,6], digital photogrammetry [7] and laser trackers. However, the theodolite networks cannot measure automatically in real time, the digital photogrammetry has a limited field of view and the laser tracker is unable to measure more than one point concurrently. Compared with these previous methods, photoelectric scanning measurement network, having concurrent measurement capability, has developed rapidly and shown great potential in the large-scale manufacturing industry [8]. As a representative, the workshop measurement positioning system (wMPS) has been widely advocated in the manufacturing industry, and intensive research efforts have been conducted on this system [9,10,11,12,13,14].
The wMPS consists of several transmitters, receivers, a signal processor and a terminal computer. As a node of the system, the calibration for intrinsic parameters of the transmitter is a crucial link in the whole measurement process, which can directly affect the final measurement accuracy [15]. However, there are relatively few methods for calibrating the intrinsic parameters of the transmitter. The earliest calibration method was to directly measure the point coordinates in the laser surfaces of the transmitter by using the theodolite measuring system and the white screen with the grid line, and then fitting the laser surface by the least squares method, which was very inefficient and had a large manual error. Zhao et al. [16] designed an automatic calibration method by using the laser tracker and the precise control network, which improved calibration efficiency and reduced manual error. In both existing methods, the laser surfaces were represented by plane equations. In the measurement process, we discovered that the laser surface edge area had lower measurement accuracy than the middle area, which reduced the effective measurement area of the transmitter. There are many factors that cause this phenomenon, the most important one of which is that the laser surface is not an ideal plane. The distortion of a cylindrical lens in a laser module and assembly errors of a component can introduce the geometrical distortions of the laser surface.
In order to reduce the error caused by laser surface distortions, Wang et al. [17] proposed a new laser surface mathematical model. Based on the quadratic interpolation principle, the quadratic surface equation was used to represent the laser surface, and the corresponding calibration and measurement algorithms were designed [17]. This method was computationally intensive, and its improvement in measurement accuracy was limited due to the existence of the Runge’s phenomenon. Therefore, in this paper we used sub-regional interpolation to build the laser surface mathematical model, which could better represent the laser surface and had a smaller amount of calculation. The sub-regional calibration has been proposed for camera calibration in photogrammetry [18], but for the wMPS the sub-regional calibration is innovative, which is different from the camera sub-regional calibration. The camera sub-regional calibration is to divide the imaging plane of the camera into different areas, and respectively calibrate the corresponding camera parameters and distortion parameters. Because the imaging plane is rectangular, its area division and recognition are more direct. The laser surface of the transmitter is a fan-shaped surface, so the division and identification of the surface are more difficult. In this paper, the pitch angle in transmitter coordinate frame of the spatial point was used as the basis for the division and identification of the laser surface.
In the proposed method, a high-accuracy 3-axis turntable and a precise 3-D coordinate control network were introduced. The transmitter for calibration needs be fixed on the turntable so that different regions of the laser surface can scan the spatial fixed control points along with the rotation of the turntable. The key to this process is to establish the corresponding coordinate system of each instrument, which is utilized to obtain the coordinates of the spatial points in the coordinate system of the transmitter under different rotation angles of the turntable. The laser surface can be divided into different regions, and each region will have a group of intrinsic parameters so that we can create a database of intrinsic parameters. In the following measurement process, the region identification arithmetic and error compensation algorithm were designed based on the database of intrinsic parameters. A series experiments were been carried out and confirmed that the novel calibration method can improve the measurement accuracy, especially at the edges of the laser surface. This paper is organized as follows: Section 2 gives a description of the sub-regional calibration mechanism. Section 3 presents the mathematical solution, including the sub-regional calibration algorithm and the error compensation algorithm. The experimental platform is set up, and experiments are prepared in Section 4 to verify the validity of the novel calibration method. Finally, in Section 5, we present concluding remarks and potential future improvements.

2. The Principle of Calibration

The novel calibration process includes the establishment of the corresponding coordinate frames and the sub-regional calibration. The detailed calibration process is as follow.

2.1. The Establishment of the Corresponding Coordinate Frames

Establishing the corresponding coordinate frame of each instrument is the key to the novel calibration method. A high-accuracy 3-axis turntable and a precise 3-D coordinate control network are introduced in this method. As shown in Figure 1, there are three coordinate frames introduced, including the turntable coordinate frame, the reference coordinate frame and the transmitter coordinate frame. The three coordinate frames are respectively defined as follows. The precise 3-D coordinate control network is established by the laser tracker, so the coordinate frame of laser tracker can be treated as the reference coordinate frame, which is defined as Ol-xlylzl (CL). The turntable coordinate frame Or-xryrzr (CR) is defined as follows. The middle shaft of the turntable is defined as X axis while the inner and outer shafts are treated as YZ axis based on the right-hand rule. The origin O is the rotation center. For the transmitter coordinate frame Ot-xtytzt (CT), the rotating shaft of the transmitter is defined as Z axis. The origin O is the intersection of laser plane 1 and Z axis. X axis is in laser plane 1 at the initial time and perpendicular to Z axis. Y axis is determined according to the right-hand rule.
Before the sub-regional calibration, a precise calibration with the reference frame Ol-xlylzl (CL) is carried out to obtain the spatial relationship between the transmitter and the turntable. As shown in Figure 1, we first need to establish a precise control network with the laser tracker. We then utilize the laser tracker to obtain the spatial relationship between the turntable and the reference coordinate frame based on fitting algorithms. What calls for special attention, is that the position of laser tracker needs to stay fixed during this process. As illustrated in Figure 2, with the help of laser tracker, we can acquire different circular locus’ centers of the three rotation shafts based on circle fitting, and then fit the three direction vectors together based on the collinear constraints and vertical constraints utilizing the least square method. The result of the fitting process can be expressed as a rotation matrix defined as R L R The coordinates of the turntable rotation center can then be obtained by the sphere fitting method, which is defined as T L R . Based on transformation principle of coordinate system, we have
{ R R L = ( R L R ) 1 T R L = ( R L R ) 1 T L R
where R R L and T R L can represent the spatial rotation and translation vector from the turntable coordinate frame Or-xryrzr (CR) to the reference coordinate frame Ol-xlylzl (CL).
Still keeping the position of the laser tracker unchanged, we then fixed the transmitter on the turntable. Several spatial points were placed in the space and the coordinates of each point were measured by the laser tracker. For each spatial point, there were two corresponding scanning angles in the transmitter. As illustrated in Figure 2, with the coordinates of spatial fixed points and corresponding scan angles of each point, plane constraint equations could be established. Based on space resection, the spatial relationship between the transmitter coordinate frame and the reference coordinate frame could then be calculated. The rotation matrix and translation vector from the transmitter coordinate frame Ot-xtytzt (CT) to the reference frame Ol-xlylzl (CL) are defined as R T L and T T L . Finally, with the rotation matrixes and translation vectors of the turntable and transmitter with respect to the reference coordinate frame, the rotation matrix and translation vector from the transmitter to the turntable could be obtained, which are defined as R T R and T T R , and we have
{ R T R = R T L ( R R L ) 1 T T R = R T L ( R R L ) 1 T R L + T T L

2.2. The Sub-regional Calibration Mechanism

With the corresponding coordinate frames established and the spatial relationship between the transmitter and the turntable acquired, the sub-regional calibration process can be carried out in the follow process. First we need keep the transmitter and the turntable unchanged and utilize the precise control network to calibrate the transmitter preliminarily, which can be called the pre-calibration. For each control point, there will be two corresponding scanning angles. With the coordinate values and corresponding scan angles of the control points, multiple plane constraint equations can be acquired. Based on the space resection, we can acquire a group of intrinsic parameters. We define this group of data as the initial intrinsic parameters, which will be utilized for the region identification of target points. Then with the help of the turntable, we begin to execute the formal calibration. With the rotation of the turntable, we can ensure that the whole laser surface can scan the control network. Based on the theory of relative motion, the clockwise rotation of the turntable along the X axis can be imagined as a counterclockwise rotation of the control network along the X axis of the turntable. As shown in Figure 3a, one control point can be transformed into different virtual points due to the change in rotation angle. With the spatial relationship between the transmitter and the reference frame and the coordinates of control points under the reference frame, we can acquire the pitch angle of the virtual points in the transmitter coordinate frame, which will be the standard of region division. For each virtual point, there will be two corresponding scanning angles. In each region, the coordinates of virtual points in the reference coordinate frame and the corresponding scanning angles can be converted to the plane constraint equations. Based on space resection, the intrinsic parameters in corresponding region can be calculated with the L-M optimization algorithm [17]. In this way, we can build a database of intrinsic parameters, which will be utilized in the error compensation process.
Then in the following measurement process, the initial intrinsic parameters are first utilized to acquire the pitch angle of the target point, which is the evidence of region identification. Based on forward intersection, the laser surfaces that have scanned the target point can establish the multiple constraints and we can then calculate the initial coordinates of target point, which can help us acquire the pitch angle of the target point. We can then carry out the error compensation based on the database of intrinsic parameters we have established. The pitch angle of the target point can decide which one of those groups of intrinsic parameters will be selected for error compensation. Finally utilizing the corresponding intrinsic parameters, we acquire a more accurate set of coordinate values for the target point based on forward intersection. The whole process of the sub-regional calibration and the error compensation is explicated in the Figure 3b.

3. The Mathematical Solution

3.1. The Sub-regional Calibration Algorithm

First we need calculate the initial intrinsic parameters, which is called the pre-calibration. Suppose there are M points in the control network. Based on space resection, utilizing the precise control network we can set up multi-plane constraints as follows:
d ij = ( a i 0 b i 0 c i 0 d i 0 ) ( R T ( θ i ) 0 0 1 ) ( R 0 T 0 0 1 ) ( x j y j z j 1 ) 2 ( i = 1 , 2 1 j M )
The i is the index of the laser surface of the transmitter, the R 0 and T 0 is the initial orientation parameters of the transmitter in the reference coordinate frame. The R T ( θ i ) is the rotation matrix determined by the scanning angle θ i . The R 0 and T 0 can be expressed as r01, r02, r03, r04, r05, r06, r07, r08, r09, tx0, ty0, tz0 and based on properties of the rotation identity matrix we have the equations below:
{ f 1 = r 01 2 + r 02 2 + r 03 2 1 = 0 f 2 = r 04 2 + r 05 2 + r 06 2 1 = 0 f 3 = r 07 2 + r 08 2 + r 09 2 1 = 0 f 4 = r 01 r 04 + r 02 r 05 + r 03 r 06 = 0 f 5 = r 01 r 07 + r 02 r 08 + r 03 r 09 = 0 f 6 = r 04 r 07 + r 05 r 08 + r 06 r 09 = 0
With Equations (3) and (4), we introduce 14 unknown coefficients, including the initial intrinsic parameters and orientation parameters of the transmitter, while 2M + 6 constraints are also introduced. By means of the Lagrange multiplier approach, the objective function can be written as:
F = i = 1 2 j = 1 M [ ( d i j ) 2 ] + λ i = 1 6 f i 2
Utilizing the L-M algorithm we can acquire the initial intrinsic parameters, which are expressed as a10, b10, c10, d10, a20, b20, c20, d20 [16].
We then start the formal calibration process. According to sub-regional calibration mechanism, one control point can be transformed to different virtual points due to the change in rotation angle. With the spatial relationship between the transmitter and the reference frame, we can acquire the coordinates of virtual points in the transmitter coordinate frame by the equation bellows:
X T = R T R R X ( R R L X L + T R L ) + T T R
XT ( x t y t z t ) T is the coordinates of virtual point in the transmitter coordinate frame and XL ( x l y l z l ) T is the coordinates of the corresponding control point in the reference coordinate frame. The ( R T R T T R ) , ( R R L T R L ) can be obtained in the process of establishing the corresponding coordinate frames and the R X is the clockwise rotation matrix along the X axis of the turntable.
With the coordinates of the virtual point XT ( x t y t z t ) T , we can calculate its pitch angle in the corresponding transmitter coordinate frame, which is defined as β and the equation is as follows:
β = arctan ( z t x t 2 + y t 2 )
The regional division of the laser surface is as follows. The fan angle ϕ of the laser surface, the angle φ between the laser surface and the axis of the transmitter are structural parameters during the manufacture of the transmitter, which are known during the calibration process. Based on the minimum angle principle, we have
c o s ϕ = c o s φ c o s γ
where the γ is the pitch angle range of laser surface in the transmitter coordinate frame, supposing that the laser surface is divided into n regions. The γ n will be the angle range of each region. Based on the pitch angle β and the angle range of each region, the virtual points can be divided into different regions.
After the region division, we can make an independent calibration in each region with the corresponding virtual points. Take one of the regions as an example. Suppose there are M virtual points in this region based on the standard of region division. The coordinates of virtual point XV ( x v j y v j z v j ) T in the turntable coordinate frame can be acquired as follows:
X V = R X ( R R L X L + T R L )
For each virtual point, there will be two corresponding scanning angles (defined as θ 1 and θ 2 ). The coordinates of the virtual point in the reference coordinate frame and the corresponding scanning angles can be converted to the plane constraint equations, which is described as follows:
d i j = ( a i cos θ i + b i sin θ i a i sin θ i b i cos θ i c i d i ) T ( R T 0 1 ) ( x v j y v j z v j 1 ) 2 ( i = 1 , 2 1 j M )
where the i is the index of the scanning optical surface of the transmitter, and the ai, bi, ci, di are the intrinsic parameters of the laser surface in this region. The θ i is the rotation angle of the laser surface. The R and T are the orientation parameters of the transmitter in the turntable coordinate frame, which can be expressed as r1, r2, r3, r4, r5, r6, r7, r8, r9, tx, ty, tz and they also follow the Equation (4).
In each region there are eight unknown coefficients for the intrinsic parameters of two laser surfaces. Considering the orientation parameters of the transmitter, there are a total of 14 unknown coefficients in this optimization. There also are 2M + 6 constraints introduced. By means of the Lagrange multiplier approach, the objective function can be written as:
F = i = 1 2 j = 1 M [ ( d i j ) 2 ] + λ i = 1 6 f i 2
This group of intrinsic parameters can then be calculated with the L-M optimization algorithm. The calibrations for the other regions are the same. After the whole regions have been calibrated, a database of intrinsic parameters has been established.

3.2. The Error Compensation Algorithm

In the measurement process, we first need to acquire the orientation parameters of each transmitter written as R t and T t , which is called the establishment of the measurement network [19]. Then, with the initial intrinsic parameters and the orientation parameters of each transmitter, we can acquire the pitch angle of the target point. Based on forward intersection, the laser surfaces that have scanned the target point can establish the constraint equation as follows:
( a i 0 b i 0 c i 0 d i 0 ) ( R T ( θ i ) 0 0 1 ) ( R t j T t j 0 1 ) ( x y z ) = 0 ( i = 1 , 2 1 j N )
where the ai0, bi0, ci0, di0 are the initial intrinsic parameters, and the i is the index of the laser surface of the transmitter. The R T ( θ i ) is the rotation matrix determined by the scanning angle θ i . The R t j and T t j are the orientation parameters of the transmitter. The j is the index of the transmitter, supposing there are N transmitters in the measurement network. The x, y, z are the initial coordinates of the target point. So, there are a total of 2N equations and only 3 unknown coefficients. With the Equation (12), we can acquire the initial coordinates of the target point and then we can calculate the pitch angle of the target point in each transmitter coordinate frame as follows:
{ ( x t 0 y t 0 z t 0 ) T = R t j ( x y z ) T + T t j β = arctan ( z t 0 x t 0 2 + y t 0 2 )
As we know, the β is the evidence of region identification. After region identification, we can carry out the error compensation based on the database of intrinsic parameters. In the error compensation process, we use the corresponding group of intrinsic parameters instead of the initial coordinates of the target point. Based on forward intersection we can then acquire new coordinates of the target point with the constraint equation as follows
( a c i b c i c c i d c i ) ( R T ( θ i ) 0 0 1 ) ( R t j T t j 0 1 ) X = 0
where the aci, bci, cci, dci are the corresponding group of intrinsic parameters from the database. The R T ( θ i ) is the rotation matrix determined by the scanning angle θ i . The R t j and T t j are the orientation parameters of the transmitter. There are a total of 2N equations and only 3 unknown coefficients as before. Finally, we can calculate the coordinates of the target point written as Xnew ( x y z ) .

4. Experiments

4.1. Set-Up of the Verification Platform

A high-accuracy 3-axis turntable and a precise 3-D coordinate control network were used to carry out the novel calibration method. The angle position accuracy of 3-axis turntable is within 1″ along each rotation axis. The control network was established by the Leica AT901-LR laser tracker, whose measurement accuracy is 15 μm + 6 μm/m. The transmitter was designed and manufactured by Tianjin University, China. The angle measurement accuracy of the transmitter is 2.6″ (based on three standard deviations) [20]. The platform for calibration using the Leica AT901-LR laser tracker and the 3-axis turntable is shown in Figure 4.

4.2. Verification of the Shaft Rotation Stability

The wMPS is a type of distributed measurement system based on the principle of angle intersection, so the angle-measurement performance of the single transmitter is the core factor of this system, which is evaluated by the angle-measurement repeatability error. Before the formal calibration, we need verify the stability of angle-measurement performance of the transmitter with the change of transmitter posture. With the help of the turntable, we tested the angle-measurement performance of the transmitter at different inclination angles respectively, and the results of the comparison tests are shown in Figure 5. The ranges of the measurement repeatability error of the scanning angles 1 and 2 were respectively 0.38″ and 0.39″, which indicated that within a certain range, the angle-measurement performance of the transmitter had good stability with the change of inclination angle.

4.3. The Sub-Regional Calibration

After the verification platform had been set up, a precise calibration was first carried out to obtain the spatial relationship between the transmitter and the turntable. Then, the sub-regional calibration was carried out and as a contrast we also utilized the conventional calibration method [21] to calibrate the intrinsic parameters of the transmitter. Based on the measurement principle of the wMPS, we need calibrate at least two transmitters to form an effective measurement network. As shown in the Table 1, the intrinsic parameters of the two transmitters were acquired by utilizing the conditional method, which could also be treated as the initial intrinsic parameters of the sub-regional calibration.
Then, utilizing the sub-regional calibration, we established the databases of intrinsic parameters for the two transmitters, which are shown in the Table 2 and Table 3 respectively.

4.4. Verification of the Error Compensation

To verify the viability and effectiveness of the proposed error compensation algorithm, a series of comparison experiments were designed to compare the measurement precision between the conventional calibration method and the novel calibration method. The Leica AT901-LR laser tracker, whose accuracy is 15 μm + 6 μm/m, was required to supply standard coordinates. Two transmitters with a distance of 5 m from each other were set up to form a measurement network. As shown in the Figure 6a, we set up 8 sets of comparison experiments in the eight areas A, B, C, D, E, F, G, H. Based on the theory of optimal measurement area proposed by Xiong et al. [22], the intersection angle of the system with the best measurement effect was 110 degrees, that is, the A area. The measurement distances of the B, C, D and E areas were close to those of the A area, but the intersection angles of them gradually changed. The measurement distances of the F, G, and H areas were gradually increased relative to the A area, and were utilized to explore the compensation effect of the novel calibration method under the different measurement distances. As shown in Figure 6b, in each area, we set up a test point every 100 mm in height direction. There were a total of 30 test points for each set of experiments, which were numbered inclusively from 1 through 30. Each point was measured by the wMPS and the laser tracker, respectively. The measured results of the wMPS were resolved twice: the proposed error compensation algorithm and the conventional algorithm. The coordinates measured by the wMPS were translated into the laser tracker coordinate system and then directly compared with the measured results of the laser tracker.
Figure 7 shows the results of the 8 sets of comparison experiments. To reject the possibility of random error appearance, we chose the average error of each region as the evaluation index of measurement performance. Some experimental results can not be seen directly from the graph, but can be obtained through data analysis. Figure 7a–h are the corresponding experimental results for areas A to H, respectively. As can be seen from the Figure 7a, in the A area for the both calibration methods, the measurement errors of the points at the edges of the laser surface were larger than those points at the central areas. The error compensation algorithm of the novel method could reduce the measurement errors of the wMPS, especially for those points at the edges of the laser surface. From point 1 to point 10, the novel calibration method reduced the average error from 0.85 mm to 0.40 mm and from point 21 to point 30 the average error was reduced from 1.01 mm to 0.60 mm. Meanwhile, points from 11 to 20, which were at the central areas of the laser surface, had a little decrease of average error from 0.41 mm to 0.3 0mm. At the points 10 and 17, the measurement errors of the novel method were larger than those of the conventional method. The reason may be that for the sub-regional interpolation algorithm, the interpolation errors at the intersection of different regions are larger, which may be far from the true value.
As shown in Figure 7b,c, compared with the A area, the global measurement errors of the two methods in the B and C areas were larger, which was caused by the change of the intersection angle. Same as the A area, the measurement errors of points at the edges of the laser surface were still larger than those points at the central areas in this two areas. The error compensation algorithm of the novel method could reduce the measurement errors of the wMPS, but the effect was not as obvious as the A area. In the B area, from point 1 to point 10, the novel calibration method reduced the average error from 1.18 mm to 1.04 mm and from point 21 to point 30 the average error was reduced from 1.43 mm to 1.14 mm. The points from 11 to 20 had a little decrease of average error from 0.88 mm to 0.82 mm. In the C area, the average error of the points from 1 to 10 was reduced from 1.32 mm to 1.00 mm and from point 21 to point 30 the average error was reduced from 1.40 mm to 1.27 mm. The points from 11 to 20 had a little decrease of average error from 0.85 mm to 0.80 mm. As shown in Figure 7d,e, in the D and E areas, the measurement results of the wMPS in the two areas were very poor, no matter which calibration method was used. The measurement errors of the points at the edges of the laser surface were also larger than those points at the central areas, but the error compensation algorithm of the novel method had no significant improvement over the conventional method.
We can see from Figure 7a,f–h that as the measurement distance increased, the global measurement error increased gradually, and the measurement errors of the points at the edges of the laser surface were larger than those points at the central areas in these three areas. In the three areas A, F, G, the error compensation algorithm of the novel method could reduce the measurement errors of the wMPS, but the error compensation effect deteriorated as the measurement distance increased. In the F area, from point 1 to point 10, the novel calibration method reduced the average error from 1.15 mm to 0.97 mm, and from point 21 to point 30 the average error was reduced from 1.32 mm to 1.25 mm, and from point 11 to point 20 the average error reduced from 0.85 mm to 0.83 mm, which is very inconspicuous. In the G area, from point 1 to point 10, the novel calibration method reduced the average error from 1.38 mm to 1.27 mm and the points from 21 to 30 had a decrease of average error from 1.37 mm to 1.35 mm. From point 11 to point 20 the average error was reduced from 1.00 mm to 0.98 mm, which was also not obvious. In the H area, the measurement errors of the two calibration methods were close, and the novel method had not improved significantly.
From the above 8 sets of experiments, in the A, B, C, F, and G areas of Figure 7a, the novel calibration method could achieve error compensation, but in the D, E, and H areas, there was no improvement in measurement accuracy compared with the conventional method. In the A, B, C, F, and G areas, the A area had the best error compensation effect, whose intersection angle was the best in the measurement network. The error compensation effect was worse in the B and C areas than in the A area. It is possible that the optical signal recognition error of the receiver photosensitive unit was increased due to the change of the intersection angle, thereby affecting the error compensation effect. The F and G areas were the areas with the worst error compensation, but there was still a certain compensation effect. Moreover, the error compensation effect of the F area was better than that of the G area. There were two factors that affected the final error compensation effect in the two areas: the increase of the optical signal recognition error of the receiver, and the decrease of the length measurement accuracy of the transmitter due to the increase of the measurement distance. In the D, E, and H areas, the novel calibration method had no advantage over the conventional method. This might be due to the fact that the measurement error introduced by the laser surface distortion was not the main error in the areas with poor intersection angle or the far measurement areas. Other error sources such as the transmitter’s own angle measurement error, the receiver optical signal recognition error, the disturbance of air, the uneven refractive index of air, etc., might have a greater impact on the experimental results, which required further experimental exploration. In summary, the novel calibration method has an error compensation effect in the area where the intersection angle is close to the optimal junction angle and the measurement distance is not far, and the error compensation effect is better in the area closer to the optimal measurement area.
The results from the experiment verified the viability and effectiveness of the error compensation algorithm and also indicated that the novel calibration method could improve the measurement accuracy of the wMPS, especially at an edge of the laser surface within a certain measurement area, which would broaden and deepen the fields of the wMPS application.

5. Conclusions

A novel intrinsic parameters calibration method that can accomplish error compensation for photoelectric scanning measurement network has been presented. This method uses the sub-regional calibration to establish the intrinsic parameters database and then resolves the target values with error compensation based on the database. This calibration method introduces a high-accuracy 3-axis turntable and a precise 3-D coordinate control network. The calibration principle is expounded in detail, including the establishment of the corresponding coordinate frames and the sub-regional calibration. Moreover, the sub-regional calibration algorithm and error compensation algorithm are also discussed in this paper. The viability and effectiveness of the novel calibration method are verified through the comparison experiments of coordinate measurement. Experimental results reveal that the novel calibration method has an error compensation effect in the area where the intersection angle is close to the optimal junction angle and the measurement distance is not far, and that the error compensation effect is better in the area closer to the optimal measurement area. The results clearly demonstrate that the novel calibration method could improve the measurement accuracy of the wMPS, especially at the edge of the laser surface within a certain measurement area. This will broaden and deepen the fields of the wMPS application.
With the expanded high-accuracy measuring volume of the transmitter, the wMPS will be able to measure larger objects while it still keeps a high accuracy. The enlarged measurement volume also gives a hand to the improvement of the dynamic measurement accuracy of this system. However, this method still has certain limitations, that is, it has an error compensation effect within a certain measurement area. In future research, we will try to realize the error compensation in the areas of the poor intersection angle and the long measurement distance, or in this case, find out the main error source and design a new error compensation method, which can be used together with the method proposed in this paper.

Author Contributions

Conceptualization, methodology, experiment design and original draft, Z.Z. and Y.R.; experiment and analysis, Z.Z. and S.S.; software, L.Y. and J.L.; supervision and funding acquisition, J.Z.; review and editing, Z.Z.

Funding

This research was funded by National Key Research and Development Program of China (Grant No. 2017YFF0204802), National Natural Science Foundation of China (Grant Nos. 51775380, 51835007, 51721003) and Natural Science Foundation of Tianjin (Grant No. 18JCYBJC19400).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schmitt, R.H.; Peterek, M.; Morse, E.; Knapp, W.; Galetto, M.; Hartig, F.; Goch, G.; Hughes, B.; Forbes, A.; Estler, W.T. Advances in large-scale metrology—review and future trends. CIRP Ann. 2016, 65, 643–665. [Google Scholar] [CrossRef]
  2. Muralikrishnan, B.; Phillips, S.; Sawyer, D. Laser trackers for large-scale dimensional metrology: A review. Precis. Eng. 2016, 44, 13–28. [Google Scholar] [CrossRef]
  3. Molebny, V.; McManamon, P.; Steinvall, O.; Kobayashi, T.; Chen, W.B. Laser radar: historical prospective-from the East to the West. Opt. Eng. 2017, 56, 14–16. [Google Scholar] [CrossRef]
  4. Bouzakis, K.; Pantermalis, D.; Mirisidis, I.; Grigoriadou, M.; Varitis, E.; Sakellaridou, A.; Bouzaki, A.; Efstathiou, K.; Theodoridou, S.; Papacharisi, A.; et al. 3D-Laser Scanning of the Parthenon West Frieze Blocks and Their Digital Assembly Based on Extracted Characteristic Geometrical Details. J. Archaeol. Sci. 2016, 6, 94–108. [Google Scholar]
  5. Franceschini, F.; Galetto, M.; Maisano, D.; Mastrogiacomo, L. Large-Scale Dimensional Metrology (LSDM): from Tapes and Theodolites to Multi-Sensor Systems. Int. J. Precis. Eng. Manuf. 2014, 15, 1739–1758. [Google Scholar] [CrossRef]
  6. Zhang, X.H.; Zhu, Z.K.; Yuan, Y.; Li, L.C.; Sun, X.Y.; Yu, Q.F.; Ou, J.L. A universal and flexible theodolite-camera system for making accurate measurements over large volumes. Opt. Laser Eng. 2012, 50, 1611–1620. [Google Scholar] [CrossRef]
  7. Cuypers, W.; Van, G.N.; Voet, A.; Kruth, J. Optical measurement techniques for mobile and large-scale dimensional metrology. Opt. Laser Eng. 2009, 47, 292–300. [Google Scholar] [CrossRef] [Green Version]
  8. Franceschini, F.; Maisano, D. The evolution of large-scale dimensional metrology from the perspective of scientific articles and patents. Int. J. Adv. Manuf. Technol. 2014, 70, 887–909. [Google Scholar] [CrossRef]
  9. Norman, A.R.; Schönberg, A.; Gorlach, I.A.; Schmitt, R. Validation of IGPS as an external measurement system for cooperative robot positioning. Int. J. Adv. Manuf. Technol. 2013, 64, 427–446. [Google Scholar] [CrossRef]
  10. Xue, B.; Zhu, J.G.; Yang, L.H.; Zhao, Z.Y.; Ye, S.H. The application of the wMPS in airplane level measurement. Opto-Electron. Eng. 2014, 41, 22–26. [Google Scholar]
  11. Xue, B.; Zhu, J.G.; Zhao, Z.Y.; Wu, J.; Liu, Z.X.; Wang, Q. Validation and mathematical model of workspace Measuring and Positioning System as an integrated metrology system for improving industrial robot positioning. Proc. Ins. Mech. Eng. 2014, 228, 422–440. [Google Scholar] [CrossRef]
  12. Ren, Y.; Lin, J.R.; Zhu, J.G.; Sun, B.; Ye, S.H. Coordinate Transformation Uncertainty Analysis in Large-Scale Metrology. IEEE Trans. Instrum. Meas. 2015, 64, 2380–2388. [Google Scholar] [CrossRef]
  13. Guo, S.Y.; Lin, J.R.; Ren, Y.J.; Shi, S.D.; Zhu, J.G. Application of a self-compensation mechanism to a rotary-laser scanning measurement system. Meas. Sci. Technol. 2017, 8, 115007. [Google Scholar] [CrossRef]
  14. Shi, S.D.; Yang, L.H.; Lin, J.R.; Long, C.Y.; Deng, R.; Zhang, Z.Y.; Zhu, J.G. Dynamic measurement error modeling and analysis in a photoelectric scanning measurement network. Appl. Sci. 2019, 9, 62. [Google Scholar] [CrossRef]
  15. Xiong, Z.; Zhu, J.G.; Zhao, Z.Y.; Yang, X.Y.; Ye, S.H. Workspace measuring and positioning system based on rotating laser planes. Mechanics 2012, 18, 94–98. [Google Scholar] [CrossRef]
  16. Zhao, Z.Y.; Zhu, J.G.; Lin, J.R.; Yang, L.H.; Xue, B.; Xiong, Z. Transmitter parameter calibration of the workspace measurement and positioning system by using precise three-dimensional coordinate control network. Opt. Eng. 2014, 53, 81–83. [Google Scholar] [CrossRef]
  17. Wang, Q.; Zhu, J.G.; Xue, B.; Zhao, Z.Y. The fitting and evaluation technique of scanning light surface for photoelectric scanning measurement network. Chin. J. Opto-Electron. Eng. 2013, 6, 84–90. [Google Scholar]
  18. Shen, X.L.; Wang, Z.; Liu, C.J.; Fu, L.H. A new technique for high precision sub-regional camera calibration based on checkerboard pattern. Chin. J. Meas. Sci. Instrum. 2016, 7, 342–349. [Google Scholar]
  19. More, J.J. The Levenberg–Marquardt algorithm: implementation and theory. Proc. Bienn. Conf. Numer. Anal. 1978, 105–116. Available online: https://www.researchgate.net/publication/313149172_The_Levenberg-_Marquardt_algorithm_implementation_and_theory (accessed on 1 January 1977).
  20. Xiong, Z.; Zhu, J.G.; Geng, L.; Ren, Y.; Yang, X.Y.; Ye, S.H. Verifcation of angle measuring uncertainty for workspace measuring and positioning system. Chin. J. Sens. Actuat. 2012, 25, 229–235. [Google Scholar]
  21. Zhao, Z.Y.; Zhu, J.G.; Xue, B.; Yang, L.H. Optimization for calibration of large-scale optical measurement positioning system by using spherical constraint. JOSA A 2014, 31, 1427–1435. [Google Scholar] [CrossRef] [PubMed]
  22. Xiong, Z.; Zhu, J.G.; Ren, Y.J.; Lao, D.B.; Gen, L.; Yang, X.Y.; Ye, S.H. Analysis and design of the best layout based on the network measurement of wMPS. In Proceedings of the 10th International Symposium on Measurement and Quality Control (ISMQC), Osaka, Japan, 5–9 September 2010. [Google Scholar]
Figure 1. The schematic of the sub-regional calibration.
Figure 1. The schematic of the sub-regional calibration.
Sensors 19 02117 g001
Figure 2. The spatial calibration between the transmitter and the turntable.
Figure 2. The spatial calibration between the transmitter and the turntable.
Sensors 19 02117 g002
Figure 3. The sub-regional calibration mechanism and the error compensation algorithm. (a) The experimental schematic; (b) The algorithm structure.
Figure 3. The sub-regional calibration mechanism and the error compensation algorithm. (a) The experimental schematic; (b) The algorithm structure.
Sensors 19 02117 g003
Figure 4. The experiment platform for the novel calibration method.
Figure 4. The experiment platform for the novel calibration method.
Sensors 19 02117 g004
Figure 5. The results of the angle-measurement performance tests.
Figure 5. The results of the angle-measurement performance tests.
Sensors 19 02117 g005
Figure 6. The schematic of the comparison experiments of coordinate measurement. (a) The space layout; (b) The distribution of test points in each group.
Figure 6. The schematic of the comparison experiments of coordinate measurement. (a) The space layout; (b) The distribution of test points in each group.
Sensors 19 02117 g006
Figure 7. The measurement errors of the two calibration methods (the blue line represents the conventional method and the red line represents the novel method). (a) The experimental results of A area; (b) The experimental results of B area; (c) The experimental results of C area; (d) The experimental results of D area; (e) The experimental results of E area; (f) The experimental results of F area; (g) The experimental results of G area; (h) The experimental results of H area.
Figure 7. The measurement errors of the two calibration methods (the blue line represents the conventional method and the red line represents the novel method). (a) The experimental results of A area; (b) The experimental results of B area; (c) The experimental results of C area; (d) The experimental results of D area; (e) The experimental results of E area; (f) The experimental results of F area; (g) The experimental results of G area; (h) The experimental results of H area.
Sensors 19 02117 g007
Table 1. The initial intrinsic parameters of the two transmitters.
Table 1. The initial intrinsic parameters of the two transmitters.
Transmittera1b1c1d1a2b2c2d2
10−0.6957450.7182890−0.779335−0.0024870.626602−0.119252
20−0.6910380.7228180−0.695678−0.0105810.7182751.32968
Table 2. The database of intrinsic parameters for the transmitter 1.
Table 2. The database of intrinsic parameters for the transmitter 1.
Regiona11b11c11d11a12b12c12d12
β < 12 0−0.6956730.7183580−0.779260−0.0024880.626696−0.119247
12 β < 4 0−0.6957090.7183240−0.779297−0.0024850.626649−0.119249
4 β 4 0−0.6957570.7182770−0.779341−0.0024900.626599−0.119251
4 < β 12 0−0.6957810.7182540−0.779373−0.0024910.626555−0.119254
β > 12 0−0.6958170.7182190−0.779411−0.0024870.626508−0.119257
Table 3. The database of intrinsic parameters for the transmitter 2.
Table 3. The database of intrinsic parameters for the transmitter 2.
Regiona21b21c21d21a22b22c22d22
β < 12 0−0.6911110.7227490−0.695606−0.0105800.7183451.329637
12 β < 4 0−0.6910740.7227830−0.695642−0.0105860.7183101.329663
4 β 4 0−0.6910240.7228320−0.695687−0.0105840.7182661.329661
4 < β 12 0−0.6910020.7228520−0.695714−0.0105780.7182411.329710
β > 12 0−0.6909640.7228870−0.695749−0.0105820.7182051.329731

Share and Cite

MDPI and ACS Style

Zhang, Z.; Ren, Y.; Yang, L.; Lin, J.; Shi, S.; Zhu, J. A Sub-Regional Calibration Method That Can Accomplish Error Compensation for Photoelectric Scanning Measurement Network. Sensors 2019, 19, 2117. https://doi.org/10.3390/s19092117

AMA Style

Zhang Z, Ren Y, Yang L, Lin J, Shi S, Zhu J. A Sub-Regional Calibration Method That Can Accomplish Error Compensation for Photoelectric Scanning Measurement Network. Sensors. 2019; 19(9):2117. https://doi.org/10.3390/s19092117

Chicago/Turabian Style

Zhang, Zhenyu, Yongjie Ren, Linghui Yang, Jiarui Lin, Shendong Shi, and Jigui Zhu. 2019. "A Sub-Regional Calibration Method That Can Accomplish Error Compensation for Photoelectric Scanning Measurement Network" Sensors 19, no. 9: 2117. https://doi.org/10.3390/s19092117

APA Style

Zhang, Z., Ren, Y., Yang, L., Lin, J., Shi, S., & Zhu, J. (2019). A Sub-Regional Calibration Method That Can Accomplish Error Compensation for Photoelectric Scanning Measurement Network. Sensors, 19(9), 2117. https://doi.org/10.3390/s19092117

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