# Research on Mapping Error Control of Underground Space Mobile LiDAR Constrained by Cooperative Targets

Authors: Huiyun Liu; Yongqiang Li; Baoshan Wang; Zengzhang Guo
Journal: Journal of Sensors
 (2022)

---

## Abstract

In the underground environment, it is difficult to obtain spatial three-dimensional data because of occlusion and its complexity. Mobile light detection and ranging (LiDAR) measurement technology has the ability to obtain three-dimensional spatial information quickly and accurately, but in the underground environment, because of the lack of global navigation satellite system (GNSS) signal in the integrated navigation system, the measurement accuracy decreases with the increase of time. In this paper, an extended Kalman filter-based loose mode is constructed using the real coordinates of the center of gravity of the target and the measurement information of the laser foothold. It provides an additional space-time reference for the integrated navigation system of the mobile LiDAR system and restricts the errors. The simulation results show that the proposed method can effectively improve the measurement accuracy of mobile LiDAR system and enhance the 3D spatial perception ability of underground space.


---

## Body

## 1. Introduction

Traditional underground space measurement mainly adopts engineering measurement methods such as total station, measuring robot, artificial static and discrete observation [1–3], and cloud point/image [4–6]. Among them, total station, measuring robot has high accuracy and good reliability, which can reach the millimeter level. The accuracy of photogrammetry can reach decimeter level. However, when traditional measurement methods obtain data, they face many problems such as low efficiency, poor safety of surveyors, and high cost of manpower and material resources, the accuracy of information acquisition is not high, and the details are incomplete [3]. The acquisition method of underground space information is inconsistent with the current demand for reliable, accurate, and fast acquisition of three-dimensional underground space information. Mobile light detection and ranging (LiDAR) technology has the advantages of fast data acquisition speed and high precision and has been widely used in 3D spatial information acquisition. Mobile LiDAR technology adopts the active measurement method, which does not depend on illumination and can efficiently obtain 3D data of underground space. However, the mobile LiDAR system relies on the accurate pose information of position and orientation system (POS). When outdoor global navigation satellite system (GNSS) signals are good, POS system can provide centimeter-level positioning accuracy [7]. In the underground space, the GNSS signals are weak or blocked, resulting in mobile LiDAR measurement data quality over time cannot meet the measurement accuracy requirements [8–10].Although mobile LiDAR technology has been applied in underground spatial information measurement, there are still challenges for large-scale 3D data acquisition in long tunnels and indoor spaces without GNSS signals due to the lack of significant features or complex environments. The current research on the rapid acquisition and mapping of underground space information can be divided into four categories: mobile measurement system for measuring and its improvement [11–21], laser based on simultaneous localization and mapping (SLAM) algorithm [22–25], monocular/binocular vision sensor [26], and RGBD (Red, Green, Blue, Depth map) vision sensor mobile robot system [27–30] and its corresponding unmanned aerial vehicle (UAV) [31–35] platform.The improved underground space information acquisition equipment based on the traditional mobile measurement system has the advantages of high relative accuracy, no need to form a closed-loop condition, and small error accumulation of long-time data acquisition. However, it also has the disadvantages of low automation, no real-time processing, and extra work. The acquisition of spatial data in underground and indoor areas without GNSS signal based on SLAM system has the advantages of high automation, good real-time data, and fast acquisition speed. However, this method also has the disadvantages of large relative errors of data and large accumulation of long-time acquisition errors, especially in special environments with no closed loop for a long time, insignificant features, or dynamic changes in space.In GNSS-denied environment, the traditional vehicle robots used for underground space positioning have high precision, but it is difficult to achieve high dynamic continuous positioning [36]. Although the radio frequency identification (RFID) [37], ZigBee technology [38], infrared technology [39], Bluetooth [40], Ultra Wideband Technology (UWB) [41], Wi-Fi, and other wireless positioning method developed in recent years can solve the problem of dynamic positioning, it has low precision, small coverage, and high cost. Another underground space positioning method, pseudostar positioning can achieve centimeter-level accuracy when the signal is good, but the signal is easily interfered, and it is not suitable for positioning in long linear space such as subway tunnel [42].In this paper, the authors analyzed the measurement error characteristics of the mobile LiDAR and then constructed the GNSS/IMU (inertial measurement unit)/target observation based on the current research. The cooperative targets are arranged to restrain the error drift of the mobile measurement system. The simulation experiment is analyzed and verified, and the appropriate target spacing is given for reference.


## 2. Theory of Method

Kalman filter (KF) is an optimal estimation algorithm based on hidden Markov model, which is simple and easy, and has been widely used in navigation, control, and data fusion. The recursive method is adopted to update the state estimate through the weighted average of the observation data and the prior values, without storing the measured values during the whole observation process. Therefore, KF is suitable for dynamic modeling. In the discrete Kalman filter, its state equation and measurement equation are shown in(1)Xk=Φk/k−1Xk−1+Γk/k−1Wk−1,Zk=HkXk+Vk,whereXk is the state vector, the dimension is n×1; Zk is the measurement vector, the dimension is m×1; ∅k/k−1, Γk/k−1, and Hk are, respectively, the state transition matrix n×n, system noise allocation matrix n×l, and measurement matrix m×n at the corresponding moment; Wk−1 and Vk are, respectively, the system noise vector l×l and the measurement noise vector m×l; they are Gaussian white noises which are not correlated with each other, satisfying the expected normal distribution of zero.In the mobile LiDAR system, the accuracy of the point cloud is mainly determined by the spatiotemporal reference provided by POS. To solve the problem of accurate acquisition of point cloud data in the environment of poor GNSS signal or missing, three-dimensional cooperative target is used to provide the spatiotemporal reference for the trajectory, and the extended Kalman filter (EKF) algorithm is used for data fusion to improve the trajectory calculation accuracy and reduce the error drift. However, for the navigation system composed of GNSS, IMU, and target, the state vector transfer and the conversion between observation quantity and state quantity are strongly nonlinear. The Kalman filter algorithm is only applicable to the linear data and cannot be directly used to process the nonlinear data.For the extended Kalman filter in nonlinear space, the premise is to assume a model for the discrete state space and replace the system state matrix and observation matrix of the standard model with a nonlinear function.(2)Xk=fXk−1+Γk−1Wk−1,Zk=hXk+Vk,where Wk−1 and Vk are system noise and observation noise; they both obey Gaussian white noise with a mean value of zero. Γk/k−1 system noise allocation matrix n×l. fXk and hXk are m- and n-dimensional nonlinear vector functions, respectively.The system transfer functionsfx and hx in the system state equation are expanded by first-order Taylor series, respectively, and the first-order terms are reserved for approximation. A new linear system and an indirect filtering system are formed by state deviation and measurement deviation, and then, the best state is estimated by Kalman filter. The filtering model is as follows:
(3)∆Xk−=Φk/k−1∆Xk−1+,Pk−=Φk/k−1Pk−1+Φk/k−1T+Γk/k−1Qk−1Γk/k−1T,Kk=Pk−HkTHkPk−HkT+Rk−1,∆Xk+=∆Xk−+KkZk−Hk∆Xk−,Pk+=I−KkHkPk−,where∆Xk− is the difference of the state vector estimated value from ∆Xk−1+ to ∆Xk−; Φk/k−1 is the Jacobian matrix corresponding to the nonlinear vector equation fx. Pk− is the covariance matrix from Pk−1+ to Pk−; Kk is the Kalman filter gain matrix; Hk is the Jacobian matrix corresponding to the nonlinear vector equation hx. ∆Xk+ is the state estimation value difference updated from ∆Xk− to ∆Xk+; Pk+ is the error covariance matrix from Pk− to Pk+ update. Qk−1 and Rk are the system noise covariance matrix and the observation noise covariance matrix at the corresponding filter moment, respectively.According to the optimal estimation of state deviation and initial state value at timeK, the optimal state estimation value at time K can be calculated as
(4)Xk=∆Xk+Xk−,or
(5)Xk=∆Xk+fXk−10.


## 3. Error Analysis of Mobile LiDAR for Underground Space

The mobile LiDAR system is a measurement system integrated with multiple sensors such as GNSS, IMU, and laser scanners. When the mobile LiDAR system is used to obtain scene information such as underground space and indoors, there are measurement errors due to many factors. The error sources are mainly divided into laser measurement error, spatiotemporal synchronization error, target system error, and integrated navigation error. Laser scanner has relatively high measurement accuracy, and the measurement error is generally within 2 cm, within the limited measurement distance such as indoor and underground space, which has limited impact on the overall measurement. Laser scanners can be corrected by laser ranging so that their errors become random noise of a smaller magnitude [42]. The spatiotemporal synchronization error includes the spatial synchronization error caused by the eccentricity and the installation angle error of the multisensor integration and the time synchronization error corresponding to the time reference of each sensor. In many literatures and researches, the error of spatial synchronization is treated as a kind of systematic error [43] and calibrated through self-calibration [44–48] or control point-based calibration methods [49–52]. Spatial synchronization can obtain the precise geometric transformation relationship among the sensors through the calibration method, which can reduce or even eliminate the error. The time synchronization error is caused by the performance difference and sampling frequency difference of various sensors, which can be effectively contained by the precise time synchronization and frequency difference. The spatiotemporal synchronization error can generally be controlled within 1-2 cm, which has a limited influence on the overall measurement. The positioning and orientation system (POS) error is the most difficult one to deal with, which is closely related to the environment, like denied the GNSS signals [42]. Cooperative target error is an extra error introduced in underground and indoor spaces to restrict spatial positioning and provide spatial reference. This error does not include a system error, and its error is considered as an accidental error. Moreover, the target data has high accuracy and the navigation system related data has low accuracy; the final result can meet the accuracy requirements of the underground navigation system by introducing target data into the navigation system to control the overall navigation error. The cooperative target adopts sensitive materials with good reflectivity to laser pulses and selects easy-to-recognize graphics, which is conducive to automatically identify the center position of the target from the laser point cloud. A large number of studies have been carried out in this area, among which the robust target center positioning algorithm based on M estimation [13] can achieve the target center positioning of 1-2 mm. In addition, the absolute position of the target is accurately measured by the precision total station, and the error is generally less than 2 cm, so it can be seen that the cooperative target error has limited influence on the overall measurement. In the mobile LiDAR system, in the loosely coupled mode of GNSS/IMU integrated navigation system, smoothing the POS trajectory can achieve relatively accurate spatial positioning and attitude determination, and the error presents a certain periodicity (Figure 1(a)), but the overall error is relatively small, generally within 1-5 cm. However, in the underground and indoor blind environment, GNSS signals are missing, and the pose of the mobile LiDAR system is independently provided by IMU, which will rapidly bring very serious accumulated errors with the increase of time (Figure 1(b)). It can be seen that in the underground blind environment, the effective containment and correction of IMU estimation errors have become the primary task to solve the underground space positioning.Figure 1
POS integrated navigation solution error characteristics.
(a)
Combined calculation error
(b)
IMU alone calculating error


## 4. Trajectory Correction Model of Mobile LiDAR with Cooperative Target Constraints

In the underground space without GNSS signals, the mobile LiDAR system has no significant difference in point cloud acquisition error and spatiotemporal synchronization error with the outdoor GNSS signal in good condition. However, in the integrated navigation error, because there is no GNSS, it only relies on IMU for trajectory estimation, and the error cannot be effectively controlled and increases sharply. When the mobile LiDAR system obtains underground space information, it usually adopts a method similar to the ground environment, which requires the participation of GNSS. However, after entering the blind environment, GNSS signal disappears. When the data is solved, GNSS will no longer participate in the relevant calculations due to the absence of reliable GNSS signal. This study is based on the extended Kalman filter algorithm using the loose coupling mode and derives the trajectory error solution model based on the graphical target constraint (Figure2).Figure 2
POS location loose coupling model with graphic target constraints.


### 4.1. IMU Trajectory State Model

Among the measurement errors of gyroscope and accelerometer in the IMU, the installation errors and the scale factor errors can be measured and compensated. Considering the simplicity of the system, only the random errors that cannot be compensated are described.Differential equation of attitude error:(6)ϕ̇=ϕ×ωinn+δωinn−δωibb,where∅ is the attitude error function. δωinn is the error between the calculated navigation coordinate system of the strapdown inertial navigation and the actual navigation coordinate system, and δωibb is the error generated by the gyroscope measurement.In the navigation coordinate system, the differential equation of velocity error is(7)δv̇n=fsfn×ϕ+vn×2δωien+δωenn−ωien+ωenn×δvn+δfsfb+δgn,where δfsfb is the error generated by the acceleration measurement difference, δωien is the calculation error of the earth rotation speed, δωenn is the calculation error of the navigation system rotation relative to the earth-fixed coordinate system, and δgn is the gravity model error. The above error calculation formulas are as follows:
(8)δfsfb=f~sfb−fsfb,δωien=ω~ien−ωien,δωenn=ω~enn−ωenn,δg=g~n−gn.The position error differential equation has different forms in different directions, which are(9)δL̇=1RM+hδvn−vNRM+h2δh,δλ=secLRN+ḣδvE+vEsecLtanLRN+hδL−vEsecLRN+h2δh,δh=δvH,where L and λ are longitude and latitude and vN, vE, and vH, respectively, are the projections of velocity in the three directions of the local coordinate system. The IMU sensor error model in this paper includes gyroscope drift error and accelerometer error, and its error differential equation is
(10)ε̇b=0,ε̇r=−1Tεr+ωr,∆̇=−1T∆+ωa,where εb is the random drift of the gyroscope, εr is the first-order Markov process of the gyroscope error, T is the correlation time, and ∆ is the accelerometer measurement error.According to the IMU attitude error equation (Equation (6)), velocity error equation (Equation (7)), position error equation (Equation (9)), and IMU sensor error (Equation (10)), the motion state equation of Kalman filter is constructed. In the equation of state, attitude error, velocity error, position error, and IMU sensor error, a total of 18-dimensional parameters are the components of the state parameters, as shown in
(11)X=φvpεbεr∆gT,(12)φ=φEφNφH,v=vEvNvH,p=pEpNpH,εb=εxεyεz,εr=εxεyεz,∆g=∆x∆y∆z.The general expression of IMU motion state equation is(13)Ẋ=FtXt+GtWt.In the expression,Ft is the state transition matrix corresponding to the state equation, for any moment t,
(14)F=FNFs0FM,where FN is the system dynamic matrix with 9 parameters of attitude, speed, and position. Fs is a matrix of 9∗9, and its form in the strapdown inertial navigation system is
(15)Fs=CbtCbt03×3Cbt06×603×39×9,where Cbt is the attitude transformation matrix between the local geographic coordinate system and the platform coordinate system. Fs is a matrix of 9∗9, and its form in the strapdown inertial navigation system is
(16)Fs=CbtCbt03×3Cbt06×603×39×9,whereFM is the state transition parameter of the inertial navigation model, and its specific form is a diagonal matrix
(17)FM=diag01×3−1Tgx−1Tgy−1Tgz−1Tax−1Tay−1Taz,where Tgx,Tgy,Tgz is the time of the gyroscope in the x,y,z direction and Tax,Tay,Taz is the time of the acceleration in the x,y,z direction.G
t is the control vector function of system noise in the state matrix, and at any time t is
(18)G=Cbt03×303×303×303×303×303×3I3×303×303×303×3I3×318×9.W
t is the noise matrix
(19)W=ωgxωgyωgzωbxωbyωbzωaxωayωazT,whereωgx,ωgy,ωgz,ωbx,ωby,ωbz,ωax,ωay,ωaz is the white noise of the gyroscope, sensor platform, and accelerometer in the x,y,z direction.


### 4.2. Mobile LiDAR Trajectory Observation Model

GNSS and IMU form an integrated navigation system. In the indirect Kalman filter mode, the IMU track estimation results are taken as the initial observation, and GNSS observation results are taken as the information to calculate the track estimation errors, and then, the navigation results are obtained. GNSS observations and IMU measurements are carried out independently of each other, and there is no correlation. Meanwhile, due to the limitation of physical conditions, the antenna center of GNSS receiver does not coincide with the coordinate center of IMU, so it is necessary to calculate the result of GNSS difference decomposition into the IMU coordinate system, to participate in the combined solution of GNSS/IMU/target. The target center coordinates are obtained in other ways and correspond to the UTC time of POS trajectory, but there is still a difference between the target center point cloud and the coordinate system of POS measurement center. Therefore, it needs to be calculated inversely according to the LiDAR positioning equation to obtain the corresponding IMU position at the moment, which participate in the GNSS/IMU/target combination calculation.


#### 4.2.1. GNSS Observation

The position results of GNSS differential observations are normalized to the IMU coordinate system through an arm effect, as shown in Figure3.Figure 3
GNSS arm lever effect reduction correction.The reduced equation is(20)δl=rGNSSe−RIMUe,where δl is the arm length vector from IMU to GNSS, RIMUe is the position vector of IMU in the geocentric earth-fixed coordinate system, and rGNSSe is the position vector of GNSS in the geocentric earth-fixed coordinate system.There is a difference in the length of the arm between GNSS and IMU, and there is also a difference in the speed of the two when the carrier moves. To integrate GNSS and IMU data, GNSS information needs to be transformed to unify the reference benchmark.In the geocentric earth-fixed coordinate system, the position vector of GNSS can be expressed asrGNSSe, and the position vector of IMU is rIMUe. Then, the position difference between GNSS and IMU, namely, the length of the space arm, is
(21)δl=rGNSSe−RIMUe.In the local navigation coordinate system, there are(22)δln=δLEδLNδLHT=Cnbδlb.Therefore, the relationship between GNSS observation result reduction and correction to the inertial navigation coordinate system is(23)LIMU=LGNSS−δLNRMh,λIMU=λGNSS−δLEsecLRNh,hIMU=hGNSS−δLH.In general, the installation positions of GNSS receiving antenna and IMU do not change, which can be regarded as a constant value. Then, the derivative of the relative geocentric earth-fixed coordinate system on both sides of Equation (21) is obtained:
(24)dδldt=rGNSSedt−drIMUedt.NotedrIMUe/dt as venIMU and rGNSSe/dt as venGNSS. The above equation projects to the inertial navigation coordinate system, where
(25)δvln=venGNSS−venIMU=Cnbωebb×δlb.


#### 4.2.2. Target Observation

In the initial calculation of the mobile LiDAR, although the long-term lack of GNSS signals caused a severe drift in navigation trajectory information, the trajectory maintained a relatively high accuracy within a short time interval within the target range. According to the graphic characteristics of the mobile LiDAR graphical target, the center coordinates of the target and the corresponding UTC time can be obtained semiautomatically with manual assistance. Since the obtained characteristic coordinate points of the graphical target have both the UTC time reference and the spatial reference of the global coordinate system after precision measurement, it can provide constraints on the spatial reference of the mobile LiDAR system and improve the accuracy of obtaining spatial information of the system.After calculating the point cloud for the first time, according to the target center positioning method in the literature [13], firstly, the point cloud scene is manually viewed and the range of the graphic target is selected. Secondly, the target center point and its corresponding UTC time are determined according to the iterative center positioning estimated by M. Finally, according to other measurement methods, the actual precise coordinates of the graphical target in the scene are obtained.According to the LiDAR positioning geometric model, the laser scanner rotates at high speed to obtain the coordinateXL of the point cloud in the scanning coordinate system. The coordinates of the laser foot point in the WGS-84 coordinate system can be measured according to the placement parameters, the estimated IMU center position, and the estimated attitude by POS system. The positioning geometric equation is shown in
(26)XWGS‐84=RNRMXL+XIL+XIMU,where XWGS‐84 is the coordinates of the WGS-84 coordinate system of the center foot point of the target, XIMU is the measurement center of inertial measurement unit, RN is the attitude estimation of POS system, and RM and XIL are the placement rotation matrix and eccentricity of the placement parameters, respectively.The target center coordinateXWGS‐84′ was obtained according to the traditional measurement method, and then, the inertial measurement center coordinate could be obtained by reverse calculation, as shown in
(27)XIMU=XWGS‐84′−RNRMXL+XIL.By integrating the GNSS and target measurement information, the measurement equation of Kalman filter in the combined system can be obtained as follows:(28)Z=vIMU−vGNSSXIMU−XGNSSXIMU−XT.


## 4.1. IMU Trajectory State Model

Among the measurement errors of gyroscope and accelerometer in the IMU, the installation errors and the scale factor errors can be measured and compensated. Considering the simplicity of the system, only the random errors that cannot be compensated are described.Differential equation of attitude error:(6)ϕ̇=ϕ×ωinn+δωinn−δωibb,where∅ is the attitude error function. δωinn is the error between the calculated navigation coordinate system of the strapdown inertial navigation and the actual navigation coordinate system, and δωibb is the error generated by the gyroscope measurement.In the navigation coordinate system, the differential equation of velocity error is(7)δv̇n=fsfn×ϕ+vn×2δωien+δωenn−ωien+ωenn×δvn+δfsfb+δgn,where δfsfb is the error generated by the acceleration measurement difference, δωien is the calculation error of the earth rotation speed, δωenn is the calculation error of the navigation system rotation relative to the earth-fixed coordinate system, and δgn is the gravity model error. The above error calculation formulas are as follows:
(8)δfsfb=f~sfb−fsfb,δωien=ω~ien−ωien,δωenn=ω~enn−ωenn,δg=g~n−gn.The position error differential equation has different forms in different directions, which are(9)δL̇=1RM+hδvn−vNRM+h2δh,δλ=secLRN+ḣδvE+vEsecLtanLRN+hδL−vEsecLRN+h2δh,δh=δvH,where L and λ are longitude and latitude and vN, vE, and vH, respectively, are the projections of velocity in the three directions of the local coordinate system. The IMU sensor error model in this paper includes gyroscope drift error and accelerometer error, and its error differential equation is
(10)ε̇b=0,ε̇r=−1Tεr+ωr,∆̇=−1T∆+ωa,where εb is the random drift of the gyroscope, εr is the first-order Markov process of the gyroscope error, T is the correlation time, and ∆ is the accelerometer measurement error.According to the IMU attitude error equation (Equation (6)), velocity error equation (Equation (7)), position error equation (Equation (9)), and IMU sensor error (Equation (10)), the motion state equation of Kalman filter is constructed. In the equation of state, attitude error, velocity error, position error, and IMU sensor error, a total of 18-dimensional parameters are the components of the state parameters, as shown in
(11)X=φvpεbεr∆gT,(12)φ=φEφNφH,v=vEvNvH,p=pEpNpH,εb=εxεyεz,εr=εxεyεz,∆g=∆x∆y∆z.The general expression of IMU motion state equation is(13)Ẋ=FtXt+GtWt.In the expression,Ft is the state transition matrix corresponding to the state equation, for any moment t,
(14)F=FNFs0FM,where FN is the system dynamic matrix with 9 parameters of attitude, speed, and position. Fs is a matrix of 9∗9, and its form in the strapdown inertial navigation system is
(15)Fs=CbtCbt03×3Cbt06×603×39×9,where Cbt is the attitude transformation matrix between the local geographic coordinate system and the platform coordinate system. Fs is a matrix of 9∗9, and its form in the strapdown inertial navigation system is
(16)Fs=CbtCbt03×3Cbt06×603×39×9,whereFM is the state transition parameter of the inertial navigation model, and its specific form is a diagonal matrix
(17)FM=diag01×3−1Tgx−1Tgy−1Tgz−1Tax−1Tay−1Taz,where Tgx,Tgy,Tgz is the time of the gyroscope in the x,y,z direction and Tax,Tay,Taz is the time of the acceleration in the x,y,z direction.G
t is the control vector function of system noise in the state matrix, and at any time t is
(18)G=Cbt03×303×303×303×303×303×3I3×303×303×303×3I3×318×9.W
t is the noise matrix
(19)W=ωgxωgyωgzωbxωbyωbzωaxωayωazT,whereωgx,ωgy,ωgz,ωbx,ωby,ωbz,ωax,ωay,ωaz is the white noise of the gyroscope, sensor platform, and accelerometer in the x,y,z direction.


## 4.2. Mobile LiDAR Trajectory Observation Model

GNSS and IMU form an integrated navigation system. In the indirect Kalman filter mode, the IMU track estimation results are taken as the initial observation, and GNSS observation results are taken as the information to calculate the track estimation errors, and then, the navigation results are obtained. GNSS observations and IMU measurements are carried out independently of each other, and there is no correlation. Meanwhile, due to the limitation of physical conditions, the antenna center of GNSS receiver does not coincide with the coordinate center of IMU, so it is necessary to calculate the result of GNSS difference decomposition into the IMU coordinate system, to participate in the combined solution of GNSS/IMU/target. The target center coordinates are obtained in other ways and correspond to the UTC time of POS trajectory, but there is still a difference between the target center point cloud and the coordinate system of POS measurement center. Therefore, it needs to be calculated inversely according to the LiDAR positioning equation to obtain the corresponding IMU position at the moment, which participate in the GNSS/IMU/target combination calculation.


### 4.2.1. GNSS Observation

The position results of GNSS differential observations are normalized to the IMU coordinate system through an arm effect, as shown in Figure3.Figure 3
GNSS arm lever effect reduction correction.The reduced equation is(20)δl=rGNSSe−RIMUe,where δl is the arm length vector from IMU to GNSS, RIMUe is the position vector of IMU in the geocentric earth-fixed coordinate system, and rGNSSe is the position vector of GNSS in the geocentric earth-fixed coordinate system.There is a difference in the length of the arm between GNSS and IMU, and there is also a difference in the speed of the two when the carrier moves. To integrate GNSS and IMU data, GNSS information needs to be transformed to unify the reference benchmark.In the geocentric earth-fixed coordinate system, the position vector of GNSS can be expressed asrGNSSe, and the position vector of IMU is rIMUe. Then, the position difference between GNSS and IMU, namely, the length of the space arm, is
(21)δl=rGNSSe−RIMUe.In the local navigation coordinate system, there are(22)δln=δLEδLNδLHT=Cnbδlb.Therefore, the relationship between GNSS observation result reduction and correction to the inertial navigation coordinate system is(23)LIMU=LGNSS−δLNRMh,λIMU=λGNSS−δLEsecLRNh,hIMU=hGNSS−δLH.In general, the installation positions of GNSS receiving antenna and IMU do not change, which can be regarded as a constant value. Then, the derivative of the relative geocentric earth-fixed coordinate system on both sides of Equation (21) is obtained:
(24)dδldt=rGNSSedt−drIMUedt.NotedrIMUe/dt as venIMU and rGNSSe/dt as venGNSS. The above equation projects to the inertial navigation coordinate system, where
(25)δvln=venGNSS−venIMU=Cnbωebb×δlb.


### 4.2.2. Target Observation

In the initial calculation of the mobile LiDAR, although the long-term lack of GNSS signals caused a severe drift in navigation trajectory information, the trajectory maintained a relatively high accuracy within a short time interval within the target range. According to the graphic characteristics of the mobile LiDAR graphical target, the center coordinates of the target and the corresponding UTC time can be obtained semiautomatically with manual assistance. Since the obtained characteristic coordinate points of the graphical target have both the UTC time reference and the spatial reference of the global coordinate system after precision measurement, it can provide constraints on the spatial reference of the mobile LiDAR system and improve the accuracy of obtaining spatial information of the system.After calculating the point cloud for the first time, according to the target center positioning method in the literature [13], firstly, the point cloud scene is manually viewed and the range of the graphic target is selected. Secondly, the target center point and its corresponding UTC time are determined according to the iterative center positioning estimated by M. Finally, according to other measurement methods, the actual precise coordinates of the graphical target in the scene are obtained.According to the LiDAR positioning geometric model, the laser scanner rotates at high speed to obtain the coordinateXL of the point cloud in the scanning coordinate system. The coordinates of the laser foot point in the WGS-84 coordinate system can be measured according to the placement parameters, the estimated IMU center position, and the estimated attitude by POS system. The positioning geometric equation is shown in
(26)XWGS‐84=RNRMXL+XIL+XIMU,where XWGS‐84 is the coordinates of the WGS-84 coordinate system of the center foot point of the target, XIMU is the measurement center of inertial measurement unit, RN is the attitude estimation of POS system, and RM and XIL are the placement rotation matrix and eccentricity of the placement parameters, respectively.The target center coordinateXWGS‐84′ was obtained according to the traditional measurement method, and then, the inertial measurement center coordinate could be obtained by reverse calculation, as shown in
(27)XIMU=XWGS‐84′−RNRMXL+XIL.By integrating the GNSS and target measurement information, the measurement equation of Kalman filter in the combined system can be obtained as follows:(28)Z=vIMU−vGNSSXIMU−XGNSSXIMU−XT.


## 4.2.1. GNSS Observation

The position results of GNSS differential observations are normalized to the IMU coordinate system through an arm effect, as shown in Figure3.Figure 3
GNSS arm lever effect reduction correction.The reduced equation is(20)δl=rGNSSe−RIMUe,where δl is the arm length vector from IMU to GNSS, RIMUe is the position vector of IMU in the geocentric earth-fixed coordinate system, and rGNSSe is the position vector of GNSS in the geocentric earth-fixed coordinate system.There is a difference in the length of the arm between GNSS and IMU, and there is also a difference in the speed of the two when the carrier moves. To integrate GNSS and IMU data, GNSS information needs to be transformed to unify the reference benchmark.In the geocentric earth-fixed coordinate system, the position vector of GNSS can be expressed asrGNSSe, and the position vector of IMU is rIMUe. Then, the position difference between GNSS and IMU, namely, the length of the space arm, is
(21)δl=rGNSSe−RIMUe.In the local navigation coordinate system, there are(22)δln=δLEδLNδLHT=Cnbδlb.Therefore, the relationship between GNSS observation result reduction and correction to the inertial navigation coordinate system is(23)LIMU=LGNSS−δLNRMh,λIMU=λGNSS−δLEsecLRNh,hIMU=hGNSS−δLH.In general, the installation positions of GNSS receiving antenna and IMU do not change, which can be regarded as a constant value. Then, the derivative of the relative geocentric earth-fixed coordinate system on both sides of Equation (21) is obtained:
(24)dδldt=rGNSSedt−drIMUedt.NotedrIMUe/dt as venIMU and rGNSSe/dt as venGNSS. The above equation projects to the inertial navigation coordinate system, where
(25)δvln=venGNSS−venIMU=Cnbωebb×δlb.


## 4.2.2. Target Observation

In the initial calculation of the mobile LiDAR, although the long-term lack of GNSS signals caused a severe drift in navigation trajectory information, the trajectory maintained a relatively high accuracy within a short time interval within the target range. According to the graphic characteristics of the mobile LiDAR graphical target, the center coordinates of the target and the corresponding UTC time can be obtained semiautomatically with manual assistance. Since the obtained characteristic coordinate points of the graphical target have both the UTC time reference and the spatial reference of the global coordinate system after precision measurement, it can provide constraints on the spatial reference of the mobile LiDAR system and improve the accuracy of obtaining spatial information of the system.After calculating the point cloud for the first time, according to the target center positioning method in the literature [13], firstly, the point cloud scene is manually viewed and the range of the graphic target is selected. Secondly, the target center point and its corresponding UTC time are determined according to the iterative center positioning estimated by M. Finally, according to other measurement methods, the actual precise coordinates of the graphical target in the scene are obtained.According to the LiDAR positioning geometric model, the laser scanner rotates at high speed to obtain the coordinateXL of the point cloud in the scanning coordinate system. The coordinates of the laser foot point in the WGS-84 coordinate system can be measured according to the placement parameters, the estimated IMU center position, and the estimated attitude by POS system. The positioning geometric equation is shown in
(26)XWGS‐84=RNRMXL+XIL+XIMU,where XWGS‐84 is the coordinates of the WGS-84 coordinate system of the center foot point of the target, XIMU is the measurement center of inertial measurement unit, RN is the attitude estimation of POS system, and RM and XIL are the placement rotation matrix and eccentricity of the placement parameters, respectively.The target center coordinateXWGS‐84′ was obtained according to the traditional measurement method, and then, the inertial measurement center coordinate could be obtained by reverse calculation, as shown in
(27)XIMU=XWGS‐84′−RNRMXL+XIL.By integrating the GNSS and target measurement information, the measurement equation of Kalman filter in the combined system can be obtained as follows:(28)Z=vIMU−vGNSSXIMU−XGNSSXIMU−XT.


## 5. Experiment Analysis

To verify the effectiveness of the proposed method in this paper, two sets of experiments were adopted for verification. In the first experiment, the authors obtain the real trajectory in the open-air environment, and the GNSS signal was removed artificially to verify the degree of conformity between the calculated trajectory and the real trajectory under target constraints. In the second experiment, the authors collected the data in the subway tunnels without GNSS signal and then compared and analyzed the measured values based on target constraints and the total station observations, and the degree of coincidence of the two test results was verified.


### 5.1. Experimental Scheme

Experiment 1: the test area is located in the newly built road environment in Jiaozuo High-tech Zone, Henan Province, with a length of 400 m. The portable flatbed vehicle is used to carry the SSW-3 mobile LiDAR system. Data are collected along nonmotorized lanes on both sides of the road. The test area is open over the sky, and there are no surrounding environments such as tall buildings and water areas that have adverse factors on GNSS signals. Plane graphical targets were set at 25 m intervals along the edge of the road, and the coordinates of the target center were accurately measured by the total station with the measurement error less than 2 cm. The time of data acquisition was October 2018, GNSS frequency was 10 Hz, IMU frequency was 200 Hz, data acquisition time was 1878 s, and the trajectory length was 836 m. Figure4 shows the point cloud image and target layout in the test area.Figure 4
Point cloud and target arrangement in test area 1.Experiment 2: the test area is located in a section of Beijing Metro Line 15. A rail car was used to carry the same type of mobile LiDAR system. The data was collected by manually dragging along the track. Cross-shaped plane targets are arranged at 50 m intervals on the tunnel wall. The center point of each plane target is precisely measured by the total station, and the measurement error is less than 2 cm. The time of data acquisition was August 2014, GNSS frequency was 10 Hz, IMU frequency was 200 Hz, the data acquisition time was 270 s, and the track length was 1700 m. Figure5 shows the data collection status and target cloud in the test area.Figure 5
Data acquisition and target arrangement in test area 2.The data solutions of the two test areas are similar, so this study focuses on the data of the test one for analysis and research.


### 5.2. Data Calculation

According to the trajectory processing method, it mainly includes POS trajectory precalculation, point cloud preliminary calculation and target point cloud center location, POS trajectory calculation with cooperative target constraint, and point cloud generation.


#### 5.2.1. POS Trajectory Precalculation

In the first experiment, the POS system acquired complete GNSS information and IMU information, and this paper used Kalman filter to obtain the complete trajectory as a reference. The maximum errors ofE (east), N (north), and H (height) in the trajectory calculation were 0.003 m, 0.004 m, and 0.010 m, respectively. Then, GNSS information is ignored to obtain the trajectory calculated only by IMU, which is used as the simulated trajectory of underground space without target constraints. In the simulation environment, because the drift error in POS system cannot be effectively suppressed, the E, N, and H errors of trajectory calculation are 8.523 m, 6.522 m, and 1.122 m, respectively. The results are shown in Figure 6. It can be seen that there is a large deviation between the above two scenarios over time.Figure 6
Simulated trajectory and real trajectory.


#### 5.2.2. Point Cloud Precalculation and Target Center Location

In the environment of simulated GNSS signal absence, the calculated point cloud trajectory is seriously deviated due to the lack of GNSS correction for IMU drift and other errors, and the calculated point cloud based on the pose information provided by IMU has a great deviation. However, the data acquisition time of each cooperative target data is very short. In a limited period of time, the pose calculated by IMU data still has high relative accuracy, and each target cloud has high overall consistency. According to the robust target center extraction algorithm, relatively reliable target center coordinates can be obtained at this time.


#### 5.2.3. Trajectory Solution of Cooperative Target Constraint

According to the target cloud extracted from the precalculated point cloud, the center coordinates of the target and the measurement results of the laser scanner at the corresponding UTC time were determined. Then, according to the observation model, the known reference points of the trajectory are calculated and the trajectory is drift constrained. The blue line in Figure6 is the corrected track. It can be seen that the corrected track has a high degree of coincidence with the real track, which greatly improves the reliability of the track.


#### 5.2.4. Point Cloud Solution

According to the pose information provided by the corrected trajectory, the mobile LiDAR point cloud was recalculated to obtain reliable point cloud information without GNSS signal. In this study, the length of the test area was 400 m, and the targets were arranged on both sides of the road. A total of 37 targets were set, and each target interval was 25 m. Finally, the integrated navigation trajectory is solved according to the intervals of 25 m, 50 m, 100 m, and 200 m, and the target point cloud data is finally obtained to evaluate the accuracy of the constrained point cloud results.


### 5.3. Constraint Results and Analysis

In the simulation of GNSS signal-free state, the error of pose estimation by IMU appears to drift over time. The deviations between the calculated coordinates and the measured coordinates of the 33 cooperative target centers are shown in Table1. The mean point-to-point errors in x, y, and z directions are 0.564 m, 10.698 m, and 0.010 m, respectively. The maximum point deviation is 19.925 m. The position deviation of each point is shown in Figure 7. The test track adopts a closed route. The farther away from the starting point, the greater the point error caused by IMU drift.Table 1
Accuracy statistics of underground simulated state.
ID
1
2
3
4
5
6
7
8
9
10
11
Δ
x
0.041
-0.026
0.189
0.163
0.378
0.305
0.379
0.211
-0.131
-0.268
-0.399
Δ
y
0.014
-0.001
0.167
0.478
1.051
1.983
3.150
5.095
6.999
9.158
9.259
Δ
z
0.000
-0.013
0.024
0.051
0.084
0.109
0.145
0.156
0.178
0.176
0.181
Δ
p
0.043
0.029
0.253
0.508
1.120
2.009
3.176
5.102
7.002
9.164
9.269
ID
12
13
14
15
16
17
18
19
20
21
22
Δ
x
-0.668
-0.557
-0.955
-0.835
-1.064
-1.394
-1.295
0.293
0.426
0.303
0.413
Δ
y
11.277
11.431
13.751
13.799
14.738
17.892
18.478
19.922
19.900
18.898
18.837
Δ
z
0.172
0.194
0.209
0.210
0.217
0.229
0.215
0.127
0.121
0.071
0.057
Δ
p
11.298
11.446
13.786
13.826
14.778
17.948
18.525
19.925
19.905
18.901
18.842
ID
23
24
25
26
27
28
29
30
31
32
33
Δ
x
0.393
0.252
0.355
0.236
0.288
0.423
0.480
0.371
0.601
0.517
0.683
Δ
y
17.668
17.606
16.048
14.392
12.499
12.375
11.054
10.595
9.004
8.707
6.823
Δ
z
0.030
0.035
0.003
-0.032
-0.062
-0.072
-0.076
-0.076
-0.045
-0.083
-0.078
Δ
p
17.672
17.608
16.052
14.394
12.502
12.382
11.065
10.602
9.024
8.723
6.858
ID: the target number;Δx is the east error; Δy is the north error; Δz is the height error; Δp is the point error.Figure 7
Accuracy results of underground simulated state.Target intervals of 25 m, 50 m, 100 m, and 200 m were, respectively, used to constrain the trajectory, and the point position errors obtained are shown in Figures8–11.Figure 8
Point error under the target constraint of 25 m interval.Figure 9
Point error under the target constraint of 50 m interval.Figure 10
Point error under the target constraint of 100 m interval.Figure 11
Point error under the target constraint of 200 m interval.From the analysis of the two test results, in the absence of GNSS signal, the trajectory of the mobile LiDAR system accumulates with time, and the measurement error increases rapidly. Based on the algorithm in this paper, the measurement error decreases rapidly under the constraint of cooperative target, and the two test results are consistent. The smaller the cooperative target interval is, the smaller the measurement error value is. When the cooperative target interval is about 50 m, the measurement error can meet the requirements of general indoor space measurement. With the increase of the interval, the accuracy and instability of the measurement error gradually increase.In the second experiment, the point cloud was solved. Without target constraint, the coordinates of 34 target points were compared with the measured values, and the mean point median errors in thex, y, and z directions were 4.076 m, 12.099 m, and 9.343 m, respectively. The proposed method was used to eliminate some unqualified targets and constrain the trajectory based on the measured target values. After recalculation, the median errors of the point cloud coordinates and the measured values were 0.051 m, 0.162 m, and 0.096 m in the x, y, and z directions, respectively.


## 5.1. Experimental Scheme

Experiment 1: the test area is located in the newly built road environment in Jiaozuo High-tech Zone, Henan Province, with a length of 400 m. The portable flatbed vehicle is used to carry the SSW-3 mobile LiDAR system. Data are collected along nonmotorized lanes on both sides of the road. The test area is open over the sky, and there are no surrounding environments such as tall buildings and water areas that have adverse factors on GNSS signals. Plane graphical targets were set at 25 m intervals along the edge of the road, and the coordinates of the target center were accurately measured by the total station with the measurement error less than 2 cm. The time of data acquisition was October 2018, GNSS frequency was 10 Hz, IMU frequency was 200 Hz, data acquisition time was 1878 s, and the trajectory length was 836 m. Figure4 shows the point cloud image and target layout in the test area.Figure 4
Point cloud and target arrangement in test area 1.Experiment 2: the test area is located in a section of Beijing Metro Line 15. A rail car was used to carry the same type of mobile LiDAR system. The data was collected by manually dragging along the track. Cross-shaped plane targets are arranged at 50 m intervals on the tunnel wall. The center point of each plane target is precisely measured by the total station, and the measurement error is less than 2 cm. The time of data acquisition was August 2014, GNSS frequency was 10 Hz, IMU frequency was 200 Hz, the data acquisition time was 270 s, and the track length was 1700 m. Figure5 shows the data collection status and target cloud in the test area.Figure 5
Data acquisition and target arrangement in test area 2.The data solutions of the two test areas are similar, so this study focuses on the data of the test one for analysis and research.


## 5.2. Data Calculation

According to the trajectory processing method, it mainly includes POS trajectory precalculation, point cloud preliminary calculation and target point cloud center location, POS trajectory calculation with cooperative target constraint, and point cloud generation.


### 5.2.1. POS Trajectory Precalculation

In the first experiment, the POS system acquired complete GNSS information and IMU information, and this paper used Kalman filter to obtain the complete trajectory as a reference. The maximum errors ofE (east), N (north), and H (height) in the trajectory calculation were 0.003 m, 0.004 m, and 0.010 m, respectively. Then, GNSS information is ignored to obtain the trajectory calculated only by IMU, which is used as the simulated trajectory of underground space without target constraints. In the simulation environment, because the drift error in POS system cannot be effectively suppressed, the E, N, and H errors of trajectory calculation are 8.523 m, 6.522 m, and 1.122 m, respectively. The results are shown in Figure 6. It can be seen that there is a large deviation between the above two scenarios over time.Figure 6
Simulated trajectory and real trajectory.


### 5.2.2. Point Cloud Precalculation and Target Center Location

In the environment of simulated GNSS signal absence, the calculated point cloud trajectory is seriously deviated due to the lack of GNSS correction for IMU drift and other errors, and the calculated point cloud based on the pose information provided by IMU has a great deviation. However, the data acquisition time of each cooperative target data is very short. In a limited period of time, the pose calculated by IMU data still has high relative accuracy, and each target cloud has high overall consistency. According to the robust target center extraction algorithm, relatively reliable target center coordinates can be obtained at this time.


### 5.2.3. Trajectory Solution of Cooperative Target Constraint

According to the target cloud extracted from the precalculated point cloud, the center coordinates of the target and the measurement results of the laser scanner at the corresponding UTC time were determined. Then, according to the observation model, the known reference points of the trajectory are calculated and the trajectory is drift constrained. The blue line in Figure6 is the corrected track. It can be seen that the corrected track has a high degree of coincidence with the real track, which greatly improves the reliability of the track.


### 5.2.4. Point Cloud Solution

According to the pose information provided by the corrected trajectory, the mobile LiDAR point cloud was recalculated to obtain reliable point cloud information without GNSS signal. In this study, the length of the test area was 400 m, and the targets were arranged on both sides of the road. A total of 37 targets were set, and each target interval was 25 m. Finally, the integrated navigation trajectory is solved according to the intervals of 25 m, 50 m, 100 m, and 200 m, and the target point cloud data is finally obtained to evaluate the accuracy of the constrained point cloud results.


## 5.2.1. POS Trajectory Precalculation

In the first experiment, the POS system acquired complete GNSS information and IMU information, and this paper used Kalman filter to obtain the complete trajectory as a reference. The maximum errors ofE (east), N (north), and H (height) in the trajectory calculation were 0.003 m, 0.004 m, and 0.010 m, respectively. Then, GNSS information is ignored to obtain the trajectory calculated only by IMU, which is used as the simulated trajectory of underground space without target constraints. In the simulation environment, because the drift error in POS system cannot be effectively suppressed, the E, N, and H errors of trajectory calculation are 8.523 m, 6.522 m, and 1.122 m, respectively. The results are shown in Figure 6. It can be seen that there is a large deviation between the above two scenarios over time.Figure 6
Simulated trajectory and real trajectory.


## 5.2.2. Point Cloud Precalculation and Target Center Location

In the environment of simulated GNSS signal absence, the calculated point cloud trajectory is seriously deviated due to the lack of GNSS correction for IMU drift and other errors, and the calculated point cloud based on the pose information provided by IMU has a great deviation. However, the data acquisition time of each cooperative target data is very short. In a limited period of time, the pose calculated by IMU data still has high relative accuracy, and each target cloud has high overall consistency. According to the robust target center extraction algorithm, relatively reliable target center coordinates can be obtained at this time.


## 5.2.3. Trajectory Solution of Cooperative Target Constraint

According to the target cloud extracted from the precalculated point cloud, the center coordinates of the target and the measurement results of the laser scanner at the corresponding UTC time were determined. Then, according to the observation model, the known reference points of the trajectory are calculated and the trajectory is drift constrained. The blue line in Figure6 is the corrected track. It can be seen that the corrected track has a high degree of coincidence with the real track, which greatly improves the reliability of the track.


## 5.2.4. Point Cloud Solution

According to the pose information provided by the corrected trajectory, the mobile LiDAR point cloud was recalculated to obtain reliable point cloud information without GNSS signal. In this study, the length of the test area was 400 m, and the targets were arranged on both sides of the road. A total of 37 targets were set, and each target interval was 25 m. Finally, the integrated navigation trajectory is solved according to the intervals of 25 m, 50 m, 100 m, and 200 m, and the target point cloud data is finally obtained to evaluate the accuracy of the constrained point cloud results.


## 5.3. Constraint Results and Analysis

In the simulation of GNSS signal-free state, the error of pose estimation by IMU appears to drift over time. The deviations between the calculated coordinates and the measured coordinates of the 33 cooperative target centers are shown in Table1. The mean point-to-point errors in x, y, and z directions are 0.564 m, 10.698 m, and 0.010 m, respectively. The maximum point deviation is 19.925 m. The position deviation of each point is shown in Figure 7. The test track adopts a closed route. The farther away from the starting point, the greater the point error caused by IMU drift.Table 1
Accuracy statistics of underground simulated state.
ID
1
2
3
4
5
6
7
8
9
10
11
Δ
x
0.041
-0.026
0.189
0.163
0.378
0.305
0.379
0.211
-0.131
-0.268
-0.399
Δ
y
0.014
-0.001
0.167
0.478
1.051
1.983
3.150
5.095
6.999
9.158
9.259
Δ
z
0.000
-0.013
0.024
0.051
0.084
0.109
0.145
0.156
0.178
0.176
0.181
Δ
p
0.043
0.029
0.253
0.508
1.120
2.009
3.176
5.102
7.002
9.164
9.269
ID
12
13
14
15
16
17
18
19
20
21
22
Δ
x
-0.668
-0.557
-0.955
-0.835
-1.064
-1.394
-1.295
0.293
0.426
0.303
0.413
Δ
y
11.277
11.431
13.751
13.799
14.738
17.892
18.478
19.922
19.900
18.898
18.837
Δ
z
0.172
0.194
0.209
0.210
0.217
0.229
0.215
0.127
0.121
0.071
0.057
Δ
p
11.298
11.446
13.786
13.826
14.778
17.948
18.525
19.925
19.905
18.901
18.842
ID
23
24
25
26
27
28
29
30
31
32
33
Δ
x
0.393
0.252
0.355
0.236
0.288
0.423
0.480
0.371
0.601
0.517
0.683
Δ
y
17.668
17.606
16.048
14.392
12.499
12.375
11.054
10.595
9.004
8.707
6.823
Δ
z
0.030
0.035
0.003
-0.032
-0.062
-0.072
-0.076
-0.076
-0.045
-0.083
-0.078
Δ
p
17.672
17.608
16.052
14.394
12.502
12.382
11.065
10.602
9.024
8.723
6.858
ID: the target number;Δx is the east error; Δy is the north error; Δz is the height error; Δp is the point error.Figure 7
Accuracy results of underground simulated state.Target intervals of 25 m, 50 m, 100 m, and 200 m were, respectively, used to constrain the trajectory, and the point position errors obtained are shown in Figures8–11.Figure 8
Point error under the target constraint of 25 m interval.Figure 9
Point error under the target constraint of 50 m interval.Figure 10
Point error under the target constraint of 100 m interval.Figure 11
Point error under the target constraint of 200 m interval.From the analysis of the two test results, in the absence of GNSS signal, the trajectory of the mobile LiDAR system accumulates with time, and the measurement error increases rapidly. Based on the algorithm in this paper, the measurement error decreases rapidly under the constraint of cooperative target, and the two test results are consistent. The smaller the cooperative target interval is, the smaller the measurement error value is. When the cooperative target interval is about 50 m, the measurement error can meet the requirements of general indoor space measurement. With the increase of the interval, the accuracy and instability of the measurement error gradually increase.In the second experiment, the point cloud was solved. Without target constraint, the coordinates of 34 target points were compared with the measured values, and the mean point median errors in thex, y, and z directions were 4.076 m, 12.099 m, and 9.343 m, respectively. The proposed method was used to eliminate some unqualified targets and constrain the trajectory based on the measured target values. After recalculation, the median errors of the point cloud coordinates and the measured values were 0.051 m, 0.162 m, and 0.096 m in the x, y, and z directions, respectively.


## 6. Conclusion

Aiming at the measurement problem of mobile LiDAR in underground space without GNSS signal, this paper analyzes the error sources of mobile LiDAR system and concludes that the dead position and attitude error calculated by IMU is the main error source of underground mobile LiDAR. In the loose coupling mode, the trajectory error solving model based on the graphical target constraint is derived. The algorithm in this paper is verified by simulating the environment without GNSS signal and the real subway tunnel scene, respectively. The research shows that the estimated trajectory error of mobile LiDAR increases rapidly with time in the absence of GNSS signal. With the constraint of target, the measurement error can be effectively controlled. When the distance between the graphic targets is about 50 m, a more stable and reliable measurement result can be obtained. Based on the error correction of the target constraint, although it can achieve effective correction of the pose to a certain extent, each graphical target only can be used as an effective constraint point, which has limited restraint ability on the estimated trajectory and less restraint on the posture. In the later study, the cooperative target and the corresponding correction model will be expanded to apply to more complex environment. And a more robust correction model that can simultaneously correct the pose will be adopted.



---

*Source: 8690532-2022-10-10.xml*