Abstract

A three-dimensional (3D) laser scanner with characteristics such as acquiring huge point cloud data and noncontact measurement has revolutionized the surveying and mapping industry. Nonetheless, how to guarantee the 3D laser scanner precision remains the critical factor that determines the excellence of 3D laser scanners. Hence, this study proposes a 3D laser scanner error analysis and calibration-method-based D-H model, applies the D-H model method in the robot area to the 3D laser scanner coordinate for calculating the point cloud data and creatively derive the error model, comprehensively analyzes six external parameters and seven inner structure parameters that affect point cloud coordinator error, and designs two calibration platforms for inner structure parameters. To validate the proposed method, we used SOKKIA total station and BLSS-PE 3D laser scanner to attain the center coordinate of the testing target sphere and then evaluate the external parameters and modify the point coordinate. Based on modifying the point coordinate, comparing the point coordinate that considered the inner structure parameters with the point coordinate that did not consider the inner structure parameters, the experiment revealed that the BLSS-PE 3D laser scanner’s precision enhanced after considering the inner structure parameters, demonstrating that the error analysis and calibration method was correct and feasible.

1. Introduction

Three-dimensional (3D) laser scanning technology [1] is a new spatial data acquisition technology that has revolutionized the surveying and mapping industry. Compared with the 3D photogrammetry and triangle photogrammetry technology, 3D laser scanning technology offers advantages like high measurement accuracy, simple data processing, and broad range. Meanwhile, 3D laser scanning technology is a cost-effective and practical solution for 3D measurements [2, 3] because of a low-energy laser pulse, strong ability to resist outside light interference, and no effect on the human eye.

While using a 3D laser scanner, the precision and index of the instrument have strict requirements, and how to ensure the precision of the 3D laser scanner that can fulfill the actual use of requirements is imperative. Currently, many scientific research institutions are researching the theory and technology of 3D laser scanning and have made some achievements. A study [4] explored the factors affecting the performance of a 3D laser scanner from the aspects of measuring distance, the color of the object to be measured, and laser incident angle. Another study [5] assessed the quality of the data collected by the ground 3D laser scanner through experiments and obtained the plane coordinates and elevation accuracy of the scanning point of the system. In addition, a study [6] performed theoretical analysis on multiple factors affecting the measurement error of a 3D laser scanner. Furthermore, the measurement accuracy of the Riegl LMS-Q140I-80 3D laser scanner was evaluated in another study [7], and the experimental results revealed that the actual measurement accuracy corroborated the nominal measurement accuracy. All the studies mentioned above primarily investigated the measurement performance and precision of a 3D laser scanner.

In recent years, some scholars started using the self-checking calibration method to validate the system error of 3D laser scanners. A study [8] assumed that the system error type of the ground 3D laser scanner was similar to that of the total station and, thus, proposed the concept of the 3D laser scanner calibration model, including nine parameters such as rotation angle, translation between three shafts of a 3D laser scanner, and instrument addition and multiplication constant. In addition, a study [9] proposed an additional seven-parameter model and the self-checking process of the 3D free network by using the theodolite error model, as well as executed the checking of the Faro 880 3D laser scanner. Another study [10] enhanced the error model previously constructed and reconstructed the mathematical model with 19 additional parameters. A study [11] explored the impact of several crucial factors and their control methods on the accuracy of Leica’s HDS 3D laser scanner. A study [12] analyzed the influence of six error factors on the coordinates of scanning points through a theoretical simulation test and then validated the detection model with simulated numerical values. Moreover, a study [13] built a scanner self-checking model with 10 parameters, examined three 3D laser scanners, and attained the ranging error, horizontal angle error, and vertical angle error of 3D laser scanners. A study [14] analyzed the influence of the positioning error of an airborne LIDAR system on the positioning accuracy and determined the errors of the positioning parameters among the GPS, laser scanner, and inertial navigation unit by using the self-checking calibration method. Furthermore, a study [15] summarized the determination method of the 3D laser scanning system error and proposed a self-check calibration method to calibrate the instrument system error. To date, no uniform standard exists for the parameters in the self-checking method.

Starting from the working principle and instrument structure of a 360° laser scanner developed by SureStar Technology Co., Ltd., a study [16] comprehensively analyzed the source of angle measurement error of such equipment, emphatically analyzed the impact of installation eccentricity of the photoelectric encoder on the system’s measurement accuracy, and deduced and established the angle error model. Based on the correlation between position and angle of each sensor in the vehicle-mounted 3D laser scanner system developed by Wuhan University, a study [17] proposed the concept of combined check and calibration of the vehicle-mounted mobile measurement system. The combined check and calibration of the laser ranging sensor and IMU are primarily implemented by using the check and calibration method based on the control point and feature surface. To date, some studies have considered the impact of the internal structure of the 3D laser scanner or sensor installation on measurement error; however, there is a lack of an overall and essential error analysis model.

The error analysis of the 3D laser scanner primarily analyzed the impact of the external influence factors on the measurement error, using the method of mathematical modeling, self-calibration, or combination calibration method for system calibration—all these methods [18, 19] can enhance the system accuracy to a certain extent; however, these methods just partially considered the influence on the equipment’s accuracy due to the system’s internal structure. This study innovatively introduces the D-H modeling method of robotics into the measurement error analysis [20, 21] and derives the error influencing factors of the 3D laser scanner, including the influence of external parameters and internal structure parameters on the measurement accuracy. To ensure the measurement accuracy of the system to the maximum extent, two sets of internal parameter calibration platforms are designed. To validate the system’s accuracy improvement effect, the high-precision measurement method is used for verification, and the test results demonstrate the method’s efficacy.

2. D–H Error Model Establishment and Analysis

2.1. 3D Laser Scanning Principle

The fundamental principle of the 3D laser scanner is that, first, a drive scanner is used to rotate in the axial direction by the axial motor and the axial rotation angle α is obtained by using an angle encoder and, second, a drive scanner is used to rotate in the radial direction by the radial motor and the radial rotation angle β is obtained by using an angle encoder; at the same time, the laser sensor measures the fly time and evaluates the distance S to the target. The point P space coordinates (shown in Figure 1) calculating formula is as follows:

Distance calculating formula is , where is the speed at which a laser travels through the atmosphere and is the round-trip time of laser flying.

2.2. Error Formula Derivation
2.2.1. D-H Model Establishment Method

D–H modeling method is a standard method to represent robots and model robot motions. Assumedly, the robot comprises a series of joints and linkages, which could be sliding or rotating and can be placed in any order and any plane. First, a reference coordinate system is assigned to each joint. Then, the steps for the transformation from one joint to the next are determined. Second, all transformations from the base to the first joint and then from the first joint to the second joint until to the last joint are combined to obtain the total transformation matrix of the robot. Figure 2 shows the D-H modeling method.

2.2.2. Calculating Model Establishment of Point Cloud

The 3D laser scanning principle presents the formula of how to calculate the coordinates of spatial point clouds. The formula is derived on the basis of the simplified model, which only considers the measurement distance S, axial scanning angle α, and radial scanning angle β, but many internal structural factors are not considered. The D-H modeling method is introduced to entirely derive the coordinate solution formula of the point cloud space.

The following definitions are given to fulfill the D-H model establishment principle:Point P is assumed as the target point of the measured objectPoint Z is assumed to be zero of the laser ranging sensorPoint O0 is assumed as the intersection point of the end face of the cross target and target laser axisPoint O1 is assumed as the intersection point of the end face of the cross target and axial motor axisPoint O2 is assumed as the intersection point of the axial motor axis and radial motor axisPoint O3 is assumed as the intersection point of the radial motor axis and measurement face of the laser sensorO0, O1, O2, and O3 are assumed as the origin, and coordinate system 0, coordinate system 1, coordinate system 2, and coordinate system 3 are established successivelyThe measuring distance of the laser ranging module is assumed as d, axial rotation angle as α, radial rotation angle as β, the distance from the laser ranging distance sensor zero point to the YZ plane of the coordinate system as L1, the distance from the laser ranging distance sensor zero point to the XZ plane of the coordinate system as L2, the distance from the coordinate system 1 origin O1 to the coordinate system 3 origin O3 as L3, the distance from the coordinate system 0 origin O0 to the coordinate system 2 origin O2 as L4, and the distance from the coordinate system 0 origin O0 to the coordinate system 1 origin O1 as L5ai−1 is assumed as the measurement distance from Zi−1 to Zi in the direction Xi−1, αi−1 as the rotation angle from Zi−1 to Zi in the direction, Xi−1, di as the measurement distance from Xi−1 to Xi in the direction Zi, and ϑi as the rotation angle from Xi−1 to Xi in the direction Zi

Figure 3 shows the internal structure of the 3D laser scanner and the relevant definitions of the D-H modeling auxiliary coordinate system.

Table 1 lists the D-H modeling parameters.

The homogeneous coordinate transformation matrix between coordinate system {i} and coordinate system {i − 1} is

We assume any point in the coordinate system 3, and then, the 3D coordinate calculating formula of any point P0 of coordinate system 0 is as follows:

, , , , and are the fixed distance, known quantity; d is the measurement distance of the laser sensor, α is the rotation angle of the radial motor, and β is the rotation angle of the axial motor. Table 2 shows the D-H modeling parameter.

The coordinate of point P is in coordinate system 3, which is the laser sensor coordinate system.

The coordinate of point P in coordinate system 0, which is the 3D laser scanner coordinate system, is shown as follows:

2.2.3. Point Cloud Error Formula Derivation

The 3D laser scanner point cloud calculation formula contains two types of parameters. One is the transformation parameter of the scanner coordinate system and the external coordinate system, which is called the external orientation parameter and comprises three translation parameters and three rotation parameters . The other is the internal error parameters of the instrument; the corresponding distance measurement error of parameter is , the corresponding distance measurement error of parameter is , the corresponding distance measurement error of parameter is , the corresponding distance measurement error of parameter is , the corresponding distance measurement error of parameter is , the corresponding rotation angle error of parameter α is , and the corresponding rotation angle error of parameter β is . In addition, set and are the coordinates of any point in the external coordinate system and the coordinates in the scanner coordinate system, respectively. Then, the following formula is satisfied between the external coordinate points and the coordinate points in the scanner coordinate system:

According to formula (4) and system error of the 3D laser scanner, the coordinate calculating formula of point P in the scanner coordinate system is shown follows:

3. Initial Value Determination for an Unknown Parameter

3.1. Initial Value Determination for an External Parameter

The external orientation parameter contains three translation parameters and three rotation parameters. To evaluate the initial value of the external orientation parameter, we pick two pairs of homonymy points , and , , which should satisfy formula (5):

The result using the abovementioned equation minus the equation below is

Based on the characteristics of the orthogonal matrix and formula (8), calculating the three angle elements , the translation can be evaluated by substituting them into any formula in equation (7), and the parameter can be taken as the initial value of the external orientation parameter.

3.2. Initial Value Determination for an Internal Parameter
3.2.1. Calculating the Internal Distance Parameter

Based on formula (6), there are primarily seven internal parameters of the 3D laser scanner, , , , , , , and , and the seven internal parameters include distance error and angle error. The distance parameter value was obtained according to the design drawing in the theory case; however, the actual situation is the scanner processing structure and assembly process inevitably introduced errors, so using the value of the design drawing was not desirable. In this study, an internal distance parameter calibration platform was designed, including a high-precision clamping tool and a digital level (Figure 4). The Vernier caliper was used to precisely measure the real value , , , , and of parameters , , , , and . Then, using the equation to calculate the distance error, the result was , , , , and .

3.2.2. Calculating the Internal Angle Parameter

In this study, we designed the internal angle parameter calibration platform, including high-precision linear guide, high-precision clamping tool, and calibration coordinate paper (as shown in Figure 5). By moving the 3D laser scanner on the linear guide rail, the 2n sets coordinate values and of laser indicators corresponding to the axial direction and radial direction on the calibration coordinate paper were collected, and two lines and were fitted by the least-square formula:

The internal angle parameters and were precisely obtained by calculating the intersection angle between fitting the line , (as shown in Table 3), and the vertical line of the calibrated coordinate paper.

4. The Result of Error Verification

4.1. Calculating the External Orientation Parameter

We used a mine 3D laser scanner named BLSS-PE to calculate external orientation parameters; six standard test target balls were set up in the laboratory, the entire laboratory scene was scanned with the 3D laser scanner to extract the cloud data of test target ballpoints, and the central coordinates of the target ball were fitted with the software of the 3D laser scanner (as shown in Figure 6). Meanwhile, the total station was used to measure the coordinates of six standard test target balls, and the center coordinates of the target balls fitted by the 3D laser scanner and the coordinates measured by the total station were used as homonyms to calculate the external orientation parameters (as shown in Figure 7).

Table 4 shows the measurement coordinates of the total station and the 3D laser scanner.

Using the coordinate values of the first four targets and then using formulas (7) and (8), the external orientation parameters are calculated as follows:

Using the calculated external orientation parameters, the scan coordinates of the other two target balls were converted to the coordinate system of the total station, and the difference between them is calculated (Table 5).

The error in each direction is as follows:

The mean square error of measurement points is as follows:

4.2. Correcting Error

By evaluating the external directional parameters and combining the internal parameters determined by the calibration platform of the 3D laser scanner, the coordinates of the standard test target ball center obtained by the 3D laser scanner can be converted into the coordinate system of the total station, and the difference between the measured coordinates and the corrected coordinates can be calculated (Table 6).

The error in each direction is as follows:

The mean square error of measurement points is as follows:

Table 7 reflects the point position accuracy comparison before and after the error correction of the 3D laser scanner.

After the correction of the instrument system error, the errors in all directions of the test point were improved compared with those before the correction. Meanwhile, the error in the point position was increased from 67 mm before correction to 20 mm, and the accuracy was increased by approximately 70.1%.

5. Conclusions

To ensure that the accuracy of the 3D laser scanner fulfills the requirements, this study innovatively uses the D-H modeling method in the robotics field to derive the error model of the 3D laser scanner, along with a detailed analysis of the six external orientation parameters and seven internal structure parameters. All these parameters are involved in the coordinate conversion error of the device, machining error, electrical installation error, and tooling error. In addition, this study integrally analyzes the error sources and whether the concrete measurement results are correct; experiments revealed that the D-H error model and calibration method proposed could enhance the instrument’s accuracy significantly. Of note, the method proposed in this study analyzes most error sources, but does not consider the environmental temperature, air pressure, reflectivity, and other factors. How to integrate these factors to further enhance the accuracy of the system would be the next research direction.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially supported by the National Key Research and Development Program of China (grant nos. 2018YFE0121000 and 2020YFE0202800).