1. Introduction
A self-rotating, linear-structured-light (LSL) scanning, 3D reconstruction system mounted on a mobile robot can be applied to the 3D reconstruction of the internal dimensions of structures such as pipelines [
1,
2,
3], tunnels [
4], tanks, indoor environments [
5,
6,
7], nuclear reactor internals [
8,
9], underwater environments [
10], etc. That can help identify internal defects, damages, leakages, etc. As shown in
Figure 1, the self-rotating, LSL scanning, 3D reconstruction system is composed of an LSL module and a rotating module. The coordinate system of the self-rotating, LSL scanning, 3D reconstruction system is usually established on the rotation center of the rotating module. However, the position parameters between the coordinate system of the LSL module and the rotation center of the rotating module are difficult to determine by the mechanical installation position. Inaccurate position parameters will lead to a decrease in system measurement accuracy [
11]. Therefore, it is necessary to establish a set of calibration methods to calibrate the position parameters between the coordinate system of the LSL module and the coordinate system of the self-rotating, LSL scanning, 3D reconstruction system to improve the measurement accuracy of the system.
Scholars have proposed some calibration methods for solving the transformation matrix between the measurement coordinate system and the world coordinate system. However, most of them are the calibration methods for rotating 2D light detection and ranging (LIDAR). Zeng et al. [
12] established the measurement coordinate system when the 2D LIDAR rotation angle is zero. When installing, make each axis parallel to each axis of the motion module coordinate system as much as possible. The origins of the two coordinates coincide in the vertical direction. Since it is difficult to ensure that the coordinate axis is completely parallel, a calibration method for the deflection angle of the center line and the coordinate axis was established. They scanned a calibration plane that is perpendicular to the Z axis of the measurement coordinate system and divided the scanned data into two, then used the projection of the measured distance of each point on the plane to the Z axis. The distances are equal; the first data are used to find the offset angle, and the second data are used to verify the results, and finally, the optimal solution is obtained. However, this method only considers the skew problem and does not calibrate the position parameters of the two coordinate origins. Yu Qiqi et al. [
13] obtained the distance between the origins of the two coordinates directly through the installation and only calibrated the angle between the lines of the two origins with respect to the horizontal plane through the parallel planes at different distances from the device. Therefore, the direct position parameters of the origins of the two coordinate systems are not fully calibrated.
Alismail et al. [
14] proposed an algorithm to automatically calibrate the rotation center for the problem that the 2D LIDAR rotation center and the driving device center may not overlap, and this method does not require a special target and can be used in general measurement scenarios. This method is based on the assumption that the surface in the local neighborhood of the 2D LIDAR point cloud can be well approximated to a plane. Since the full scan is symmetrical, the full scan is divided into two half scans of the same object, and the distance from the plane to the center in the range of the corresponding angle laser point is equal, and the transformation matrix that minimizes the dissimilarity of the two data is calculated so that the calibration of the center of rotation is completed. This method does not require special calibration objects, but the algorithm is more complicated. Huang Fengshan et al. [
15] used a rectangular parallelepiped with a triangular prism in the middle channel as the calibration piece to complete the calibration. The calibration piece is large and requires special customization. Cai et al. [
16] proposed using a special calibration plate with hollow holes for calibration, This method needs to fit the coordinates of each center according to the measured data and then complete the calibration according to the position relationship of each center. The calibration methods for rotating 2D LIDAR can be used as references for calibrating self-rotating, LSL scanning, 3D reconstruction systems, but they cannot be directly applied.
Li and Xi [
17] and Xiao et al. [
18] presented a rotational laser scanner by mounting the LSL sensor to a turntable, but they did not mention the calibration of the center of rotation. Manakov et al. [
19], Wissel et al. [
20], and chi et al. [
10] designed a galvanometric laser scanner with a laser projector and proposed a model-driven calibration method, but they all assumed that the galvanometer rotation axis coincides with the line intersected by the mirror of the galvanometer and the laser plane. Yang M et al. [
21] proposed a flexible plane-constraint-based calibration method for the galvanometric laser scanner, which is effective in the vision inspection system, but this method needs to be combined with a lookup table.
Aiming at the above problems, this paper proposes a calibration method for a self-rotating, LSL scanning, 3D reconstruction system based on plane constraints. By collecting the point cloud data of the plane target, based on the basic principle that the points on the plane target should be constrained by the plane equation, position parameters between the coordinate system of the LSL module and the coordinate system of the self-rotating, LSL scanning, 3D reconstruction system were calculated.
2. Transformation between Coordinate Systems
In the establishment of the coordinate system of a self-rotating, LSL scanning, 3D reconstruction system,
, the axis
is coincident with the rotation axis of a high-precision rotating turntable, and the direction is upward. The plane
is established based on the right-hand rule. Additionally, the plane
and the plane of the coordinate system of the LSL module are set coplanar when the rotation angle is zero, while the axis
is parallel to the axis Y and the axis
is parallel to the X axis. As shown in
Figure 2,
is the coordinate system of the self-rotating, LSL scanning, 3D reconstruction system, while
is the coordinate system of the high-precision rotating turntable.
is the coordinate system of the LSL module when the rotation angle is zero.
is the coordinate system of the LSL module when the axis rotates clockwise by an angle
.
Point P is a measurement point on the laser plane when the rotation angle of the self-rotating, LSL scanning, 3D reconstruction system. Since the laser plane coincides with the plane , the coordinate of point P in the coordinate system is . The coordinate of the point P in the coordinate system is , since the Z axes of the two coordinate systems are parallel and in the same direction, so .
When the rotation angle of the self-rotating, LSL scanning, 3D reconstruction system is zero, the coordinate value of the origin O of the LSL module in the coordinate system be . Pass point O to make a vertical line perpendicular to the axis , and point A is the vertical foot. When the rotation angle of the self-rotating, LSL scanning, 3D reconstruction system is , the origin of the coordinate system of the LSL module shifts from point O to point , and the corresponding position of point A shifts to point Ai.
From the geometric relationship,
where
,
,
. According to the geometric relationship in
Figure 2, we can obtain
Therefore, the coordinate value
of the point P in the coordinate system
can be expressed as follows [
18]:
It can be seen from Equation (7) that the coordinate value of the point P in the coordinate system is related to the parameters M, N, and the coordinate value of the point P in the coordinate system . The angle can be measured by the encoder and the coordinate value of the point P in the coordinate system can be obtained by the LSL module, while the position parameters M and N can only be a rough estimate of their value range, and further calibration is required to obtain their precise values.
3. Calibration Based on Plane Constraints
This paper proposes a calibration method for self-rotating, LSL scanning, 3D reconstruction systems based on plane constraints to achieve the precise acquisition of the position parameters M and N described in
Section 2. The specific method is as follows:
(1) First, according to the design parameters and installation positions of the LSL module and high-precision rotating turntable, the position parameters are preliminarily estimated a ≤ M ≤ b, c ≤ N ≤ d.
(2) As shown in
Figure 3, the calibration plane is placed within the measurement range of the self-rotating, LSL scanning, 3D reconstruction system. The high-precision rotary turntable is rotated so that the LSL module sweeps across the plane, and at the same time, the data of the plane profile are collected in the LSL module at various angles; then, the transformation method of Equation (7) is used to unify the point cloud data to the coordinate system of the self-rotating, LSL scanning, 3D reconstruction system.
Assume that the three-dimensional coordinates of the measurement point on the calibration plane in the self-rotating, LSL scanning, 3D reconstruction system are
,
……
. These measurement points are all located on the same plane, so these points are constrained to a plane. The plane equation is
, and there is
which is
Solve the equation
and obtain the fitted plane by the least square method, the least-square solution of the equation is
The distance from the points to the fitting plane is
, where
, so the root-mean-squared error of the distance from the measuring point to the fitting plane is
Change the calibration plane position, as shown in
Figure 4, and collect the point cloud of the planes at different positions. According to the above calculation method, we can obtain the root-mean-squared error of the fitting plane at different positions
,
……
. The root-mean-squared error of each fitting plane and the average value is
Take Equation (15) as the objective function and search for the best value of M and N in the interval of a ≤ M ≤ b, c ≤ N ≤ d with a certain step length. Minimize the value of S, and use this optimized value of M and N as the result of positional parameter calibration. The Algorithm 1 is as the follows:
Algorithm 1: Calibration of self-rotating, LSL scanning, 3D reconstruction system based on plane constraints. |
Input: 1. The collected point cloud data of m groups of planes that are not parallel to each other; 2. The value range of parameter M is [a, b]. The value range of parameter N is [c, d]. |
Algorithm steps: 3. M = a, N = c, S = 100, step size is δ; 4. While M≯b do 5. While N≯c do 6. Based on the least square method to get the fitting planes, and solve the root mean square error of each group of point clouds to each fitting plane ,……; 7. If then 8. ; ; ; 9. end if 10. N = N + δ; 11. end while 12. M = M + δ; 13. end while Output: 14. The minimum S value and position parameters and . |
4. Experiment and Results
As shown in
Figure 5, the plane of the aluminum alloy plate was placed within the measurement range of the self-rotating, LSL scanning, 3D reconstruction system, and the self-rotating, LSL scanning, 3D reconstruction system was controlled to scan and collect the point cloud data of the aluminum alloy plate. The collected data were the data in the coordinate system of the LSL module. The plane position of the aluminum alloy plate was adjusted, and the point cloud data were collected three times.
According to the self-rotating, LSL scanning, 3D reconstruction system designed in
Section 2, it was preliminarily estimated that the M value would be in the interval (110, 130), and the N value would be in the interval (395, 425). Using the calibration method described in
Section 3 to process the collected data, the relationship between the parameters M, N, and the root-mean-squared error S was obtained, as shown in
Figure 6. The position parameters Mopt and Nopt corresponding to the minimum root-mean-squared error were the calibration results of the parameters M and N, as shown in
Table 1.
In order to verify the accuracy of the parameters Mopt and Nopt, obtained in the above calibration process, the test piece, as shown in
Figure 7, was processed, its characteristic dimensions were measured by a coordinate measurement machine, and the size measured by the coordinate machine was used as the true value.
In
Figure 8, the test piece was placed within the measurement range of the self-rotating, LSL scanning, 3D reconstruction system, and the self-rotating, LSL scanning, 3D reconstruction system was controlled to scan and collect the point cloud of the test piece’s surface.
To convert the collected data into the data in the coordinate system of the self-rotating, LSL scanning, 3D reconstruction system, M = 123 and N = 409, as shown in
Figure 9; the relevant feature size was then measured.
As shown in
Table 2, when M = 123 and N = 409, where M and N are the optimal value, the measurement accuracy of each feature size on the test piece is better than others. When the value of M and N is smaller than the optimal value, the measured value of the feature size along the scanning direction on the test piece is smaller than the measured value when M and N take the optimal value, while the measured value of the feature size perpendicular to the scanning direction is equivalent when M and N take the optimal value. When the value of M and N is greater than the optimal value, the measured value of the feature size along the scanning direction of the test piece is greater than the measurement value when M and N take the optimal value, and the measured value of the feature size perpendicular to the scanning direction is equivalent when M and N take the optimal value.
5. Conclusions
The coordinate system of the LSL module does not coincide with the coordinate system of the self-rotating, LSL scanning, 3D reconstruction system, and the positional parameters between them are difficult to measure through mechanical installation. This paper first analyzed the transformation relationship between the coordinate system of the LSL module and the coordinate system of the self-rotating, LSL scanning, 3D reconstruction system, then proposed a calibration method for the self-rotating, LSL scanning, 3D reconstruction system based on plane constraints, and finally, designed the validation experiment. It was confirmed by the experiment that the proposed method can accurately calculate the positional parameters between the coordinate system of the LSL module and the coordinate system of the self-rotating, LSL scanning, 3D reconstruction system, and can significantly improve the measurement accuracy of the measurement system.
Author Contributions
Conceptualization, J.Z. and G.C.; methodology, J.Z. and Y.C.; software, S.H.; validation, L.L.; formal analysis, G.W.; investigation, J.Z. and L.Y.; resources, C.F. and L.Y.; data curation, S.H.; writing—original draft preparation, J.Z.; writing—review and editing, J.Z. and L.Y.; visualization, Y.C. and L.Y. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Informed consent was obtained from all subjects involved in the study.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Zhang, J.; Pan, L.I.; Fan, S.H.; Wang, S.; Lu, L. Research and Application of 3D Reconstruction Techniques of Pipeline Based on Line Scanning. Radio Eng. 2017, 47, 66–70. (In Chinese) [Google Scholar]
- Wang, Y.; Zhang, R. In-pipe surface circular structured light 3D vision inspection system. Infrared Laser Eng. 2014, 43, 891–896. [Google Scholar]
- Tang, X.; Pang, Y.J.; Zhang, H.; Zeng, W.J.; Li, Y. Underwater Pipeline Detection by AUV Based on Monocular Vision. Robot 2010, 32, 592–600. [Google Scholar]
- Roh, T.H.; Kim, J.S.; Lee, Y.D. Development of Cross Section Management System in Tunnel using Terrestrial Laser Scanning Data. J. Korean Assoc. Geogr. Inf. Stud. 2008, 11, 90–104. [Google Scholar]
- Zhang, L.; Sun, J.; Yin, G.; Zhao, J.; Han, Q. A Cross Structured Light Sensor and Stripe Segmentation Method for Visual Tracking of a Wall Climbing Robot. Sensors 2015, 15, 13725–13751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Liang, X.; Chen, H.; Li, Y.; Liu, Y. Visual laser-SLAM in large-scale indoor environments. In Proceedings of the 2016 IEEE International Conference on Robotics and Biomimetics (ROBIO), Qingdao, China, 3–7 December 2016. [Google Scholar]
- Chen, H.; Huang, H.; Qin, Y.; Li, Y.; Liu, Y. Vision and laser fused SLAM in indoor environments with multi-robot system. Assem. Autom. 2019, 39, 297–307. [Google Scholar] [CrossRef]
- Zhao, J.P.; Feng, C.; Cai, G.; Zhang, R.; Chen, Z.; Cheng, Y.; Xu, B. Three-dimensional reconstruction and measurement of fuel assemblies for sodium-cooled fast reactor using linear structured light. Ann. Nucl. Energy 2021, 160, 108397. [Google Scholar] [CrossRef]
- Zavyalov, P.S.; Senchenko, E.S.; Finogenov, L.V.; Khakimov, D.R. A structured-light method for the measurement of deformations in fuel assemblies in the cooling ponds of nuclear power plants. Russ. J. Nondestruct. Test. 2012, 48, 705–711. [Google Scholar] [CrossRef]
- Chi, S.; Xie, Z.; Chen, W. A laser line auto-scanning system for underwater 3D reconstruction. Sensors 2016, 16, 1534. [Google Scholar] [CrossRef] [Green Version]
- Wang, C.M.; Huang, F.S.; Xue, Z. Lidar measurement system of loading robot and its calibration method. Optoelectron. Eng. 2019, 46, 190002-1–190002-8. [Google Scholar]
- Zeng, Y.; Yu, H.; Dai, H.; Song, S.; Lin, M.; Sun, B.; Meng, M.Q.H. An Improved Calibration Method for a Rotating 2D LIDAR System. Sensors 2018, 18, 497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Yu, Q.Q.; Cui, Z.S. A calibration method of scanning system based on 2D lidar. Laser Infrared 2017, 47, 1234–1237. [Google Scholar]
- Alismail, H.; Browning, B. Automatic Calibration of Spinning Actuated Lidar Internal Parameters. J. Field Robot. 2015, 32, 723–747. [Google Scholar] [CrossRef]
- Huang, F.S.; Ma, S.; Xue, Z. Rotating two-dimensional lidar measurement system and its calibration method. J. Opt. Laser 2018, 29, 67–75. [Google Scholar]
- Cai, J.; Zhao, Y.; Li, Y.H.; Xie, Y.M. Design and parameter calibration of a three-dimensional laser scanning system. J. Beijing Univ. Aeronaut. Astronaut. 2018, 44, 2208–2216. [Google Scholar]
- Li, L.; Xi, J. Free and global pose calibration of a rotating laser monocular vision sensor for robotic 3D measurementsystem. In Proceedings of the International Conference on Optics in Precision Engineering and Nanotechnology (icOPEN2013), Singapore, 22 June 2013. [Google Scholar]
- Xiao, J.; Hu, X.; Lu, W.; Ma, J.; Guo, X. A new three-dimensional laser scanner design and its performance analysis. Optik 2015, 126, 701–707. [Google Scholar] [CrossRef]
- Manakov, A.; Seidel, H.P.; Ihrke, I. A mathematical model and calibration procedure for galvanometric laser scanning systems. In Vision, Modeling, and Visualization Workshop 2011; The Eurographics Association: Berlin, Germany, 2011; pp. 207–214. [Google Scholar]
- Wissel, T.; Wagner, B.; Stuber, P.; Schweikard, A.; Ernst, F. Data-driven learning for calibrating galvanometric laser scanners. IEEE Sens. J. 2015, 15, 5709–5717. [Google Scholar] [CrossRef]
- Mao, Y.; Zeng, L.C.; Jiang, J.; Yu, C. Plane-constraint-based calibration method for a galvanometric laser scanner. Adv. Mech. Eng. 2018, 10, 168781401877367. [Google Scholar] [CrossRef] [Green Version]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).