[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Influence of Electrostatic Forces on the Vibrational Characteristics of Resonators for Coriolis Vibratory Gyroscopes
Previous Article in Journal
Automatic Indoor as-Built Building Information Models Generation by Using Low-Cost RGB-D Sensors
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 20
Issue 1
10.3390/s20010294
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Positioning Error Compensation Method for a Mobile Measurement System Based on Plane Control

by
Bo Shi
1,
Fan Zhang
1,
Fanlin Yang
1,*,
Yanquan Lyu
1,
Shun Zhang
1 and
Guoyu Li
2
1
College of Geomatics, Shandong University of Science and Technology, Qingdao 266000, China
2
Qingdao Xiushan Mobile Surveying Co., Ltd., Qingdao 266000, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(1), 294; https://doi.org/10.3390/s20010294
Submission received: 19 November 2019 / Revised: 30 December 2019 / Accepted: 2 January 2020 / Published: 4 January 2020
(This article belongs to the Section Remote Sensors)
Browse Figures
Figure 1
<p>Smooth schematic diagram during the satellite signal loss-of-lock.</p> "> Figure 2
<p>Comparison of POS error after 1-min smooth solution of satellite signal loss of lock. (<b>a</b>) Position error comparison; (<b>b</b>) attitude error comparison.</p> "> Figure 3
<p>Comparison of POS errors after 3-min smooth solution of satellite signal loss-of-lock. (<b>a</b>) Position error comparison; (<b>b</b>) attitude error comparison.</p> "> Figure 4
<p>Positioning error compensation based on plane control.</p> "> Figure 5
<p>Comparison of position error curves and Gaussian curves after GNSS signal loss-of-lock.</p> "> Figure 6
<p>VSurs-E mobile measurement system.</p> "> Figure 7
<p>The J1 teaching building experimental area control plane.</p> "> Figure 8
<p>The restaurant experiment area control plane.</p> "> Figure 9
<p>Position error curve of satellite loss of lock in the experimental area of J1 teaching building. (<b>a</b>) loss of lock for 90 s; (<b>b</b>) satellite loss of lock for 180 s.</p> "> Figure 10
<p>Compensated position error curve in the experimental area of J1 teaching building. (<b>a</b>) Satellite loss of lock for 90 s; (<b>b</b>) satellite loss of lock for 180 s.</p> "> Figure 11
<p>The satellite loss of lock 180 s positioning error compensation before and after the graph in the experimental area of restaurant. (<b>a</b>) Satellite loss of lock for 180 s; (<b>b</b>) satellite loss of lock 180 s compensation.</p> "> Figure 12
<p>3D point cloud image of the experimental area.</p> "> Figure 13
<p>(<b>a</b>) The same proportion of the national flag point cloud; (<b>b</b>) the same proportion of the stone ball top view point cloud; (<b>c</b>) the same proportion of the JIL and SZSF plane point cloud top view point cloud; (<b>d</b>) the same proportion of the window point cloud.</p> ">
Versions Notes

Abstract

:
Global navigation satellite system (GNSS)/inertial navigation system (INS) navigation technology is one of the core technologies in a mobile measurement system and can provide real-time geo-referenced information. However, in the environment measurements, buildings cover up the GNSS signal, causing satellite signals to experience loss-of-lock. At this time errors of INS independent navigation accumulate rapidly, so it cannot meet the needs of the mobile measurement system. In this paper, a positioning error compensation method based on plane control is proposed by analyzing the error characteristics of position and orientation caused by satellite signal loss-of-lock in the urban environment. This method control uses planar features and the laser point cloud positioning equation to establish an adjustment model that ignores the influence of the attitude error and finds the positioning error at the middle point of the GNSS signal loss-of-lock time period, and then compensates for the error at other points by using the characteristics of the positioning error. The experimental results show that the accuracy of the compensated laser point cloud has been significantly improved, and the feasibility of the method is verified. Meanwhile, the method can rely on the existing building plane so the method is adaptable and easy to implement.

1. Introduction

A mobile mapping system (MMS) integrates measurement sensors such as a position and orientation system (POS), a laser scanning system and a charge-coupled device (CCD) camera on a mobile carrier [1,2], enabling the fast and efficient acquisition of spatial data on surrounding targets, which is widely used in three-dimensional modeling, digital city, street view map, digital mine and other fields [3,4,5,6]. In the measurement process, MMS needs the high-precision position and orientation information obtained by POS to realize the conversion from a laser scanner coordinates system to the specified mapping coordinate system. At present, the POS of the mobile measurement system mainly adopts the Global Navigation Satellite System (GNSS)/inertial navigation system (INS) combined navigation [7]. When MMS is measured in an urban environment, GNSS signals can easily be blocked by the tall buildings in the city which cause the satellite signals to experience loss of lock. An integrated navigation system can only rely on INS independent navigation. Since INS is a continuous integration system, if navigation is performed for a long time, the inertial device errors will continue to accumulate. In this case, the POS system cannot provide accurate position and orientation information, and the accuracy of the laser scanning system will be reduced from the sub-meter to the meter level [8,9].
Domestic and foreign researchers have conducted a great deal of research on how to improve the accuracy of MMS during satellite signal loss-of-lock. The research results can be summarized into two strategies: correcting POS information by means of constraints and smoothing, or compensating for laser point cloud accuracy with feature points. The former fixes the navigation data error by adding constraints and constructing a filter equation of additional motion constraints by using the motion law of the carrier [10,11,12]. However, this method usually constructs constrained observation vectors with fixed observations. The actual observation value is not a constant, and is greatly affected by topographical and carrier motion changes, so the navigation accuracy after improvement cannot be guaranteed. The performance of the GNSS/INS integrated navigation system can also be improved by adding auxiliary sensors such as odometers, visual sensors, and LiDAR (light detection and ranging) to compensate for the INS errors [13,14]. However, the initial value of the odometer scale might be affected by temperature and tire pressure in this method is not clear, the error of the odometer scale accumulates linearly with driving distance and is greatly affected by the speed of the vehicle. The visual sensor is susceptible to light and so we cannot guarantee POS accuracy. LiDAR is less affected by the environment and is inexpensive, so it is a sensor suitable for positioning assistance. However, this method mostly uses a Kalman filter (KF) or other filter models to fuse the navigation information of LiDAR and INS to achieve positioning and orientation. Among them, for LiDAR, the estimation of carrier position and attitude is mainly obtained by scanning matching on continuously-scanned images. When there is no overlapping area between consecutively-scanned images or the image surface geometry is not good, it will seriously affect the accuracy of scanning matching and lead to a significant reduction in POS accuracy. In addition, in the post-processing of GNSS/INS combined navigation data, smoothing is used to improve the accuracy of combined navigation [15,16,17,18,19,20]. The smoothing process can directly correct the POS error in the loss-of-lock period by directly using the GNSS data from before and after without adding other observations. The effect of improving the attitude error after smoothing is significant, but positioning error increases with satellite signal loss-of-lock time, so it is difficult to ensure the accuracy of the MMS measurement. It is possible to obtain an external positioning update at a given time by calculating the coordinate difference between the pre-deployed feature points with precise coordinates and the coordinate points of the same-named feature points in the point cloud [21,22]. Feature points can also be extracted within the trajectory, and the improved Iterative Closest Point (ICP) algorithm is used to analyze the incorrect correspondence of points. Then the least squares method is used to calculate the corrected position and attitude of the trajectory feature points. Finally, the positional rotation parameters of the corresponding trajectory points are calculated by linear interpolation [23]. However, this method requires a large number of feature points to be set and requires precise feature points calculated by the traditional method, which has a significant workload and low efficiency. Moreover, extracting feature points of the same name mainly relies on manual operation. The accuracy of the extraction is highly correlated with the point cloud density. In addition, this method can only correct the error of the laser point cloud; it does not improve the accuracy of the POS, and it is difficult to obtain accurate trajectory information, so the method has certain limitations in terms of its practical application.
In this paper, a positioning error compensation method based on plane control is proposed. The method makes full use of the geographical environment in the satellite loss-of-lock area and uses the existing building plane as the control plane to calculate the positioning error at an intermediate time. Based on the analysis of characteristics for the positioning error after smoothing, the positioning error at other times can be compensated for. This not only avoids the difficulty of feature point extraction and the inaccuracy of extraction, but also obtains accurate track information. In addition, the method overcomes the limitations of the hardware and environment brought about by constraint conditions and auxiliary sensors. At the same time, it improves the work efficiency in practical applications.

2. Analysis of Position and Orientation Error Characteristics

In this paper, the proposed method compensates for the POS error by using the relationship between the smoothed POS error and the GNSS loss-of-lock time. Therefore, the characteristics of the smoothed POS error need to be analyzed. Taking a traditional forward-backward smoothing filtering algorithm as an example, the forward and backward filtering achieve smooth estimation of the period loss-of-lock of satellite signals by using the data before and after the loss of the GNSS signal [24], as shown in Figure 1.
Figure 1 shows that the smoothed results present a trend of first rising and then falling, similar to the shape of the Gaussian function distribution. The related literature research has yielded similar results [15,16,17,18,19]. This article uses a set of measured data from a mobile measurement system as an example. Navigation data were calculated using IE8.6, Novatel’s commercial software (The software is manufactured by Novatel Corporation of Alberta, Canada.), and GNSS satellite signals were set to cover a loss-of-lock of 1 min or 3 min. The results were different from those obtained without the loss-of-lock. The results are shown in Figure 2 and Figure 3, and the position and attitude error data are shown in Table 1.
The satellite loss of lock lasts 1 min, as shown in Figure 2.
The satellite loss of lock lasts for 3 min, as shown in Figure 3.
From the trend of the curve in Figure 2 and Figure 3, it can be seen that, after the smoothing process, the attitude error has a certain fluctuation. The error value is always within a thousandth of a degree, which meets the requirements of the mobile measurement accuracy. The positioning error shows a trend of increasing from zero to the maximum error and then decreasing from the maximum error to zero in the satellite loss-of-lock period. This trend is consistent with the shape proposed in this paper, which is similar to the Gaussian function distribution. When the GNSS signal experiences loss-of-lock for 1 min, the positioning error is still on the decimeter level. Although it definitely experiences an improvement, it still cannot meet the needs of mobile mapping.

3. Positioning Error Compensation Based on Plane Control

The positioning accuracy cannot meet the needs of mobile mapping after the smoothing solution. In order to improve the positioning accuracy, this paper proposes a compensation method based on plane control. In this paper, the positioning error at the intermediate moment of GNSS loss-of-lock was first calculated in the area where the GNSS signal loss-of-lock using the control plane. Then, according to a relationship similar to a Gaussian function between the smoothed positioning error and the GNSS loss-of-lock time, the positioning error at other times is compensated for. As shown in Figure 4, first of all, a total station is used to measure the coordinates of 15–30 measuring points which are evenly distributed on the space plane, and then the space plane equation of the plane is calculated. At the same time, the laser point cloud data of the same plane measured by the mobile measurement system and the smoothed position and attitude information are used to construct the laser point cloud positioning equation. Then, Based on the spatial plane equation and the laser point cloud positioning equation, a position error solving model based on plane control is established. Finally, by analyzing the characteristics of the smoothed position error and using the calculated position error to construct a Gaussian function error compensation model, the position error is compensated.

3.1. Laser Point Cloud Localization Equation

The basic steps of a mobile measurement system positioning are based on the polar diameter and angle information of the target point measured by the laser scanner. It is the first to calculate the coordinates of the target point in the laser scanner coordinate system, and then goes through a series of coordinate transformations to achieve the conversion of the target point from a laser scanner coordinate system to the WGS-84 coordinate system. The positioning equation of the laser point cloud in the WGS-84 coordinate system is:
[ X Y Z ] W G S 84 = [ X Y Z ] o e + R l e R b l [ R s b [ X Y Z ] s + [ l x l y l z ] ] ,
where [ X , Y , Z ] S T is the point cloud coordinate in the laser scanner frame; [ l x , l y , l z ] T is the amount of translation of the laser scanner frame to the body frame; R s b is the rotation matrix of the laser scanner frame to the body frame; R b l is the rotation matrix of the body frame to the local level frame; R l e is the rotation matrix of the local level frame to the WGS-84 coordinate frame; [ X , Y , Z ] oe T are the coordinates of the local horizontal coordinate origin in the WGS-84 coordinate frame; and [ X , Y , Z ] e T is the coordinates of the laser point cloud in the WGS-84 coordinate frame.

3.2. Positioning Error Solving Principle

The errors that affect the measurement accuracy of the mobile measurement system can be divided into three major components: POS error, laser scanner measurement error, and placement error after the integration of each sensor. With the development of laser scanner hardware, we must ensure that the measuring accuracy of the point position meets the specific measurement needs. The placement error between each sensor can reduce the measurement error by calibration. Therefore, the method in this paper ignores the influence of the laser scanner measurement error and the placement error after the integration of each sensor. Equation (1) can be expressed as follows:
[ X Y Z ] W G S 84 = [ X Y Z ] o e + R l e R b l [ X Y Z ] b ,
The POS error includes three positioning errors Δ X oe , Δ Y oe , Δ Y oe and three attitude errors Δ r , Δ p , Δ y , which are added to Equation (2) to obtain the POS error:
[ X Y Z ] W G S 84 = [ X Y Z ] o e + Δ [ X Y Z ] o e + R l e * × ( I + Ω l * l ) R b l * × [ X Y Z ] b ,     Ω l * l = [ 0 Δ y Δ r − Δ y 0 − Δ p − Δ r Δ p 0 ] ,
in which R l e ∗ is the rotation matrix of the local level frame to the WGS-84 coordinate frame when there is a positioning error; Ω l * l is the attitude error matrix; and R b l * is the rotation matrix of the body frame to the local level frame when there is an attitude error.
Through the propagation of error, the error model of Equation (3) can be obtained:
{ σ x = [ ( ∂ X ∂ r ) 2 σ r 2 + ( ∂ X ∂ p ) 2 σ p 2 + ( ∂ X ∂ y ) 2 σ y 2 + ( ∂ X ∂ X o e ) 2 σ X o e 2 + ( ∂ X ∂ Y o e ) 2 σ Y o e 2 + ( ∂ X ∂ Z o e ) 2 σ Z o e 2 ] 1 2 σ Y = [ ( ∂ Y ∂ r ) 2 σ r 2 + ( ∂ Y ∂ p ) 2 σ p 2 + ( ∂ Y ∂ y ) 2 σ y 2 + ( ∂ Y ∂ X o e ) 2 σ X o e 2 + ( ∂ Y ∂ Y o e ) 2 σ Y o e 2 + ( ∂ Y ∂ Z o e ) 2 σ Z o e 2 ] 1 2 σ Z = [ ( ∂ Z ∂ r ) 2 σ r 2 + ( ∂ Z ∂ p ) 2 σ p 2 + ( ∂ Z ∂ y ) 2 σ y 2 + ( ∂ Z ∂ X o e ) 2 σ X o e 2 + ( ∂ Z ∂ Y o e ) 2 σ Y o e 2 + ( ∂ Z ∂ Z o e ) 2 σ Z o e 2 ] 1 2 }
In this equation, σ X , σ Y , σ Z are the mean square errors of the laser point cloud in the X, Y, and Z directions, respectively, in the WGS-84 coordinate frame; σ r 2 , σ p 2 , σ y 2 are the mean square errors of the roll angle, pitch angle, and heading angle, respectively, and σ X o e 2 , σ Y o e 2 , σ Z o e 2 are the mean square errors of the positioning in the X, Y, and Z directions, respectively, in the WGS-84 coordinate frame.
According to the official data provided by Novatel, after IE8.6 software processing, the mean square error of the pitch angle and roll angle is 0.005°, the mean square error of heading angle is 0.008°, the mean square error of the horizontal direction is 0.01 m, and the mean square error of the vertical direction is 0.015 m when the satellite signal is not loss-of-lock. The mean square error of pitch angle and roll angle is 0.006°, the mean square error of the heading angle is 0.006°, the mean square error of the heading angle is 0.010°, the mean square error of the horizontal direction is 0.011 m, and the mean square error of the vertical direction is 0.03 m when the satellite signal loss-of-lock lasts for 60 s.
It can be seen that as the satellite signal loses lock time increases, the positioning error increases significantly, while the attitude error varies little. In addition, when the error model is used to calculate the accuracy of the laser point cloud, the effect of the positioning error on the accuracy of the laser point cloud is much greater than the attitude error. The attitude error can be smoothed to meet the needs of movement measurement. Therefore, the influence of attitude error is not considered in the calculation, so Equation (3) can be simplified as follows:
[ X Y Z ] W G S 84 = [ X Y Z ] o e + Δ [ X Y Z ] o e + R l e × R b l * × [ X Y Z ] b ,
The general form of the space plane is:
a x + b y + c z − d = 0 ,
where a , b , c is the unit normal vector of the plane and d is the distance from the origin of the coordinate to the plane. Then the reference plane parameter a , b , c , d can be obtained using the eigenvalue method [25].
The laser point cloud data on the control surface satisfy the plane equation of the control surface, so:
[ a p b p c p ] [ X Y Z ] e − d p = 0 , [ a p b p c p ] [ [ X Y Z ] o e + Δ [ X Y Z ] o e + R l e × R b l * × [ X Y Z ] b ] − d p = 0
In Equation (7), a p , b p , c p , d p represent the plane parameter of the p-th control plane. [ X Y Z ] b T is the point cloud data coordinates in the body frame. Δ [ X Y Z ] o e is the obtained positioning error.

3.3. Similar to the Gaussian Distribution Function Error Compensation Model

By analyzing the characteristics of the positioning error after smoothing, it can be observed that the relationship between the positioning error after smoothing and time is similar to the Gaussian function feature. Taking the GNSS signal loss-of-lock for 90 s as an example, the relationship between the smoothed position error curve and the Gaussian function curve is shown in Figure 5. The red curve is the smoothed positioning error curve, and the green curve is the Gaussian function curve. It can be seen that the smoothed position error curve is similar to the Gaussian function curve.
It is assumed that the satellite signal is occluded at time t 0 ~ t n and the navigation information calculated after smoothing reaches the maximum value at the intermediate time t m i d . The error and time in the satellite loss-of-lock period have the following Gaussian relationship:
y = b × e − ( t − t m i d ) 2 2 σ 2 σ = ( t n − t 0 ) 6
where b is the positioning error value at time t m i d , which can be obtained by the above method. The positioning error at other time points is calculated by a model similar to the Gaussian function, thereby compensating for the positioning error of the entire satellite signal during the loss-of-lock period.

4. Experimental Analysis

In order to verify the accuracy of the positioning error compensation method based on plane control, this paper chooses a VSurs-E mobile measurement system as the experimental platform. The feasibility of this method is demonstrated by comparing the positioning error before and after compensation and the accuracy of the laser point cloud.

4.1. Experimental Platform and Control Plane Selection

The VSurs-E road information collection and inspection system was jointly developed by Shandong University of Science and Technology and Qingdao Xiushan Mobile Measurement Co., Ltd., as shown in Figure 6. The system is composed of Z+F PROFILER 9012 laser scanner (The Z + F Profiler 9012 manufactured by Zoller & Fröhlich GmbH Corporation of Wangen im Allgäu, Germany), a SPAN-LCI (The equipment is manufactured by Novatel Corporation of Alberta, Canada.) integrated navigation system, industrial cameras and other sensors. The system uses a car as a carrier to quickly collect three-dimensional spatial data such as road facilities and road topography. The parameters of VSurs-E’s laser scanner and integrated navigation system are given in Table 2.
When the vehicle mobile measurement system is measured in an urban canyon environment, it can collect rich plane features, and there is a rich plane in the teaching building area. Therefore, the teaching building area is selected as the experimental site. In the J1 teaching building area, a total of seven groups of planes, which are the building’s smooth wall surfaces different orientations are selected as the control plane as shown in Figure 7. In addition, in the school restaurant, three planes with flat surfaces and different orientations are selected as the control plane as shown in Figure 8. The measurement of the experimental sites with the VSurs-E mobile measurement system and the reference planes were accurately extracted from the laser point cloud data. At the same time, based on accurate measurement of the control points, a high-precision total station is used to measure 15–30 points for each reference plane. These acquired data points provide data support for the parameters that compute the plane equations.

4.2. Positioning Error Simulation Analysis

In order to verify the accuracy of the method in this paper and analyze the influence of attitude error on positioning error, a simulation experiment was performed by using a set of navigation data that did not involve loss-of-lock. Firstly, the original positioning error in the navigation data was calculated using the method of this paper. Then we set the constant positioning error and the orientation error in the period when the measuring vehicle passes through the area where the control plane is laid, and the positioning error was calculated again using the above method. The positioning error at this time include the original positioning error and the set constant value positioning error. Finally, after removing the positioning error, compare it with the set constant error.
Experimental statistics are shown in Table 3. Setting different values and different positive and negative positioning errors means that the positioning error obtained by our method is different from the simulated constant positioning error in the millimeter range. It can be seen from Equation (3) that due to the influence of the attitude error, the position corresponding to each point changes when the coordinate system is rotated, so the positioning error obtained in this paper is different from the simulation constant positioning error. With the increase in attitude error, the difference between the positioning error obtained by the method and the simulated constant positioning error is gradually increasing. However, when the attitude error increases to 0.010 degrees, the difference is still in the millimeter range, which has little effect on the final solution.

4.3. Positioning Error Compensation Experiment

First, the original data of GNSS was processed. The GNSS information for the experimental area of the J1 teaching building was deleted for 90 s and 180 s, thus causing an artificial loss-of-lock. The difference between the smoothed loss-of-lock data and the normalized unlocked data is shown in Figure 9. When the satellite signal experienced loss-of-lock for 90 s, the positioning error reached 0.13 m in the ECEF-Z direction and 0.09 m in the ECEF-Y direction. When the satellite loss-of-lock time increases to 180 s, the maximum positioning error in the ECEF-Z direction was 0.15 m and in the ECEF-Y direction it was 0.15 m. Since these data come from the calibration data of mobile measurement, the measuring car needs to be stationary for a period of time before and after passing through the control plane area, thus the inertial navigation error accumulates slowly. Therefore, the positioning error is still on the decimeter level when the satellite experiences loss-of-lock for 180 s.
The positioning error obtained after compensation is shown in Figure 10. When the satellite experiences loss of lock for 90 s, the compensated positioning error ECEF-X direction is reduced from 0.03 m to 0.01 m, and in the ECEF-Y direction it is reduced from 0.09 m to 0.04 m. The ECEF-Z direction is reduced from 0.13 m to 0.02 m. When the satellite experiences loss of lock for 180 s, the compensated positioning error ECEF-X direction is reduced from 0.03 m to 0.02 m and in the ECEF-Y direction it is reduced from 0.15 m to 0.02 m, while the ECEF-Z direction is reduced from 0.15 m to 0.01 m. It can be concluded that with the increase in the time of satellite loss of lock, the compensation effect of this method becomes more obvious.
In addition, in order to verify the reliability of the method, the 180 s GNSS observation information of the restaurant experiment area was deleted, thus causing an artificial loss-of-lock. The difference between the data positioning error before and after the GNSS loss-of-lock compensation and the normal unlocked data is shown in Figure 11.
The compensated positioning error in the ECEF-X direction is reduced from 0.15 m to 0.02 m, and in the ECEF-Y direction it is reduced from 0.21 m to 0.03 m. The ECEF-Z direction is reduced from 0.13 m to 0.01 m. It can be seen that the method has good repeatability and can provide useful reference for practical engineering applications.
In order to visually show the compensation effect of laser point cloud positioning accuracy. The navigation data before and after the compensation are merged with the point cloud data measured by the laser scanner, and the coordinates of the laser foot point in the ECEF coordinate system are obtained. Then the value is plugged into the corresponding equation of the control plane to calculate the residual of each laser foot point, as follows:
ω i = ( a x i + b y i + z c i ) − d ,
The residual means square error statistics of residues before and after compensation in each plane are given in Table 4 and Table 5. The residual means square error of residues of seven planes are shown in Table 6 and Table 7. It is easy to see that the accuracy of the laser point cloud after compensation has been significantly improved.
Figure 12 is a three-dimensional point cloud image of the experimental area. a , b , c , and d marked in this figure corresponds to those in Figure 13.
Figure 13 shows some details of the three-dimensional point cloud map of the experimental area. Among them, the red point cloud indicates the laser point cloud when GNSS is not experiencing loss-of-lock. The green laser point cloud on the left indicates the laser point cloud before compensation, and the green laser point cloud on the right indicates the compensated laser point cloud. The bottom line in the figure is the size of the scale, which is used to show that the scale of the laser point cloud before and after compensation is the same. From the above experiment, it can be seen that the point cloud before compensation does not coincide with the point cloud calculated normally, and the point cloud compensated for by the method in this paper has an extremely obvious effect of overlapping with the point cloud calculated normally.
It can be seen from the conclusion of the whole experiment that the positioning error compensation method of the mobile measurement system, based on the plane control can effectively reduce the positioning error of the integrated navigation after the loss-of-lock. In addition, the accuracy of the compensated laser point cloud has been significantly improved.

5. Conclusions

In this paper, by analyzing the POS error characteristics of satellites in the unlocking state of GNSS, a method of positioning error compensation for mobile measurement systems based on plane control is proposed. The method can calculate the positioning error at an intermediate point and compensate for the entire GNSS loss-of-lock period. This not only overcomes the problem of feature point extraction and precision guarantee, but also improves the work efficiency, at the same time obtaining accurate track information. In our study, not only by comparing the accuracy of the moving measurement system before and after compensation, but also by comparing the positioning error before and after compensation and the accuracy of the mobile measurement system, we prove that this method is a reliable and effective.
In the process of solving the positioning error, when all the point cloud data are used for solving, the solving speed is slow due to the data volume being so large. In this study, the positioning error obtained may not be the optimal solution because only partial points are used to solve the solution. So, it is necessary to further study how to make full use of data to improve the accuracy of the solution. The idea behind this method can also be extended to other features. In future work, we will study the positioning error compensation of mobile measurement systems based on wider features, such as surface features suitable for tunnels.

Author Contributions

B.S. proposed the original idea, performed the preliminary experiments, drafted the manuscript, and was responsible for the arrangement of the entire work. F.Z. contributed to the main results and code implementation, and collected and analyzed the experimental data. F.Y. provided ideas to improve the quality of the paper and also revised the manuscript. Y.L. and S.Z. helped in data cross-checking, manuscript formatting and editing. G.L. designed and performed the experiments. All authors have read and agreed to the published version of the manuscript.

Funding

This work was jointly funded by the National Natural Science Foundation of China (grant no. 41930535), the National Key R and D Program of China (grant nos. 2018YFF0212203, 2018YFC1405900, 2017YFC1405006, 2016YFB0501705, and 2016YFC1401210), the Shandong Provincial Key R and D Program (grant no. 2019JZZY010429), the SDUST Research Fund (grant no. 2019TDJH103) and a Project of the Shandong Province Higher Educational Youth Innovation Science and Technology Program (grant no. 2019KJH007).

Acknowledgments

The authors thank Xiaokai Li who polished the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wan, R.; Huang, Y.; Xie, R.; Ma, P. Combined Lane Mapping Using a Mobile Mapping System. Remote Sens. 2019, 11, 305. [Google Scholar] [CrossRef] [Green Version]
  2. Cui, T.; Ji, S.; Shan, J.; Gong, J.; Liu, K. Line-Based Registration of Panoramic Images and LiDAR Point Clouds for Mobile Mapping. Sensors 2017, 17, 70. [Google Scholar] [CrossRef] [PubMed]
  3. Nivedita, S.; Sudhagar, N.; Scott, O. Development of mobile mapping system for 3D road asset inventory. Sensors 2016, 16, 367. [Google Scholar]
  4. Borja, R.C.; Silverio, G.C.; Celestino, O.; Maria, C.A. An approach to detect and delineate street curbs from MLS 3D point cloud data. Autom. Constr. 2015, 51, 103–112. [Google Scholar]
  5. Fryskowska, A.; Wroblewski, P. Mobile laser scanning accuracy assessment for the purpose of base-map updating. Geod. Cartogr. 2018, 67, 35–55. [Google Scholar]
  6. Handel, P.; Ohlsson, J.; Ohlsson, M. Smartphone-Based Measurement Systems for Road Vehicle Traffic Monitoring and Usage-Based Insurance. IEEE Syst. J. 2014, 8, 1238–1248. [Google Scholar] [CrossRef] [Green Version]
  7. Yu, M.J. INS/GPS integration system using adaptive filter for estimating measurement noise variance. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 1786–1792. [Google Scholar] [CrossRef]
  8. Puente, I.; González-Jorge, H.; Martínez-Sánchez, J.; Arias, P. Review of mobile mapping and surveying technologies. Measurement 2013, 46, 2127–2145. [Google Scholar] [CrossRef]
  9. Barber, D.; Mills, J.; Smith-Voysey, S. Geometric validation of a ground-based mobile laser scanning system. ISPRS J. Photogramm. Remote Sens. 2008, 63, 128–141. [Google Scholar] [CrossRef]
  10. Peng, K.Y.; Lin, C.A.; Chiang, K.W. The performance analysis of an AKF based tightly-coupled INS/GPS integrated positioning and orientation scheme with odometer and non-holonomic constraints. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, XXXIX-B7, 481–486. [Google Scholar] [CrossRef] [Green Version]
  11. Angrisano, A.; Gaglione, S.; Gioia, C. Performance assessment of GPS/GLONASS single point positioning in an urban environment. Acta Geod. Geophys. Hung. 2013, 48, 1–13. [Google Scholar] [CrossRef]
  12. Klein, I.; Filin, S.; Toledo, T. Vehicle constraints enhancement for supporting INS navigation in urban environments. Navigation 2011, 58, 7–15. [Google Scholar] [CrossRef]
  13. Vagle, N.; Broumandan, A.; Lachapelle, G. Multi-Antenna GNSS and Inertial Sensors/Odometer Coupling for Robust Vehicular Navigation. IEEE Internet Things J. 2018, 5, 4816–4828. [Google Scholar] [CrossRef]
  14. Broumandan, A.; Lachapelle, G. Spoofing detection using GNSS/INS/Odometer coupling for vehicular navigation. Sensors 2018, 18, 1305. [Google Scholar] [CrossRef] [Green Version]
  15. Gong, X.; Zhang, R.; Fang, J. Application of unscented R–T–S smoothing on INS/GPS integration system post processing for airborne earth observation. Measurement 2013, 46, 1074–1083. [Google Scholar] [CrossRef]
  16. Shi, B.; Lu, X.; Chen, Y.F. Improving position and attitude precision of GPS/INS by applying EKF smoothing algorithm. J. Geomat. Sci. Technol. 2012, 29, 334–338. [Google Scholar]
  17. Qiu, R.; Cheng, Y. The interpolation application of interval quartering algorithm of singular spectrum analysis iterative in GPS coordinate time series. J. Geod. Geodyn. 2015, 35, 1031–1035. [Google Scholar]
  18. Liu, S.; Sun, F.; Li, H. Forward-backward-smoothing algorithm with application to tightly coupled PPP/INS data post-processing. J. Chin. Inert. Technol. 2015, 23, 85–91. [Google Scholar]
  19. Liu, H.; Nassar, S.; El-Sheimy, N. I Two-Filter Smoothing for Accurate INS/GPS Land-Vehicle Navigation in Urban Centers. IEEE Trans. Veh. Technol. 2010, 59, 4256–4267. [Google Scholar] [CrossRef]
  20. Nassar, S.; Niu, X.; El-Sheimy, N. Land-Vehicle INS/GPS accurate positioning during GPS signal blockage periods. J. Surv. Eng. 2007, 133, 134–143. [Google Scholar] [CrossRef]
  21. Takai, S.; Date, H.; Kanai, S.; Niina, Y.; Oda, K.; Ikeda, T. Accurate registration of MMS point clouds of urban areas using trajectory. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2013, II-5/W2, 277–282. [Google Scholar] [CrossRef] [Green Version]
  22. Schaer, P.; Vallet, J. Trajectory adjustment of mobile laser scan data in GPS denied environments. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2016, XL-3/W4, 61–64. [Google Scholar] [CrossRef]
  23. Han, J.Y.; Chen, C.S.; Lo, C.T. Time-Variant registration of point clouds acquired by a mobile mapping system. IEEE Geosci. Remote Sens. Lett. 2013, 11, 196–199. [Google Scholar] [CrossRef]
  24. Kai-Wei, C.; Trung, D.T.; Liao, J.K. On-Line smoothing for an integrated navigation system with Low-Cost MEMS inertial sensors. Sensors 2012, 12, 17372–17389. [Google Scholar]
  25. Sharkarji, C.M. Least-Squares fitting algorithms of the NIST algorithm testing system. J. Res. Natl. Inst. Stand. Technol. 1998, 103, 633–641. [Google Scholar] [CrossRef] [PubMed]
Sensors 20 00294 g001 550
Figure 1. Smooth schematic diagram during the satellite signal loss-of-lock.
Figure 1. Smooth schematic diagram during the satellite signal loss-of-lock.
Sensors 20 00294 g001
Sensors 20 00294 g002 550
Figure 2. Comparison of POS error after 1-min smooth solution of satellite signal loss of lock. (a) Position error comparison; (b) attitude error comparison.
Figure 2. Comparison of POS error after 1-min smooth solution of satellite signal loss of lock. (a) Position error comparison; (b) attitude error comparison.
Sensors 20 00294 g002
Sensors 20 00294 g003 550
Figure 3. Comparison of POS errors after 3-min smooth solution of satellite signal loss-of-lock. (a) Position error comparison; (b) attitude error comparison.
Figure 3. Comparison of POS errors after 3-min smooth solution of satellite signal loss-of-lock. (a) Position error comparison; (b) attitude error comparison.
Sensors 20 00294 g003
Sensors 20 00294 g004 550
Figure 4. Positioning error compensation based on plane control.
Figure 4. Positioning error compensation based on plane control.
Sensors 20 00294 g004
Sensors 20 00294 g005 550
Figure 5. Comparison of position error curves and Gaussian curves after GNSS signal loss-of-lock.
Figure 5. Comparison of position error curves and Gaussian curves after GNSS signal loss-of-lock.
Sensors 20 00294 g005
Sensors 20 00294 g006 550
Figure 6. VSurs-E mobile measurement system.
Figure 6. VSurs-E mobile measurement system.
Sensors 20 00294 g006
Sensors 20 00294 g007 550
Figure 7. The J1 teaching building experimental area control plane.
Figure 7. The J1 teaching building experimental area control plane.
Sensors 20 00294 g007
Sensors 20 00294 g008 550
Figure 8. The restaurant experiment area control plane.
Figure 8. The restaurant experiment area control plane.
Sensors 20 00294 g008
Sensors 20 00294 g009 550
Figure 9. Position error curve of satellite loss of lock in the experimental area of J1 teaching building. (a) loss of lock for 90 s; (b) satellite loss of lock for 180 s.
Figure 9. Position error curve of satellite loss of lock in the experimental area of J1 teaching building. (a) loss of lock for 90 s; (b) satellite loss of lock for 180 s.
Sensors 20 00294 g009
Sensors 20 00294 g010 550
Figure 10. Compensated position error curve in the experimental area of J1 teaching building. (a) Satellite loss of lock for 90 s; (b) satellite loss of lock for 180 s.
Figure 10. Compensated position error curve in the experimental area of J1 teaching building. (a) Satellite loss of lock for 90 s; (b) satellite loss of lock for 180 s.
Sensors 20 00294 g010
Sensors 20 00294 g011 550
Figure 11. The satellite loss of lock 180 s positioning error compensation before and after the graph in the experimental area of restaurant. (a) Satellite loss of lock for 180 s; (b) satellite loss of lock 180 s compensation.
Figure 11. The satellite loss of lock 180 s positioning error compensation before and after the graph in the experimental area of restaurant. (a) Satellite loss of lock for 180 s; (b) satellite loss of lock 180 s compensation.
Sensors 20 00294 g011
Sensors 20 00294 g012 550
Figure 12. 3D point cloud image of the experimental area.
Figure 12. 3D point cloud image of the experimental area.
Sensors 20 00294 g012
Sensors 20 00294 g013 550
Figure 13. (a) The same proportion of the national flag point cloud; (b) the same proportion of the stone ball top view point cloud; (c) the same proportion of the JIL and SZSF plane point cloud top view point cloud; (d) the same proportion of the window point cloud.
Figure 13. (a) The same proportion of the national flag point cloud; (b) the same proportion of the stone ball top view point cloud; (c) the same proportion of the JIL and SZSF plane point cloud top view point cloud; (d) the same proportion of the window point cloud.
Sensors 20 00294 g013
Table
Table 1. Position and attitude error data after GNSS/INS smoothing solution.
Table 1. Position and attitude error data after GNSS/INS smoothing solution.
TimeECEF-X (m)ECEF-Y (m)ECEF-Z (m)Roll (deg)Pitch (deg)Yaw (deg)
1-min0.1580.0590.0330.00450.00160.0011
3-min0.4001.0050.2100.00420.00620.0065
Table
Table 2. Technical parameters.
Table 2. Technical parameters.
SensorTechnical Parameters
Laser scanner
(Z+F PROFILER 9012)		Laser class1
Ambiguity distance (m)≤0.001
Angular resolution (deg)0.0088°
POS(SPAN-LCI)Position Accuracy/Vertical Position Accuracy (m)0.01/0.015
Pitch/Roll Angle Accuracy (deg)0.005
Heading Angle Accuracy (deg)0.008
Table
Table 3. Simulation experiment results.
Table 3. Simulation experiment results.
Simulation Constant ErrorThis Paper Obtained the Positioning Error
ECEF-X (m)ECEF-Y (m)ECEF-Z (m)Roll (deg)Pitch (deg)Yaw (deg)ECEF-X (m)ECEF-Y (m)ECEF-Z (m)
0.050.050.050.0050.0050.0050.04920.04880.0498
0.100.100.100.0050.0050.0050.09920.09880.0998
0.150.150.150.0050.0050.0050.14920.14880.1498
0.050.050.050.0080.0080.0080.04870.04810.0504
0.050.050.050.0100.0100.0100.04840.04760.0495
0.050.050.05−0.005−0.005−0.0050.05080.05120.0502
−0.05−0.05−0.05−0.005−0.005−0.005−0.0492−0.0488−0.0498
Table
Table 4. Compensation for satellite loss of lock of 90 s.
Table 4. Compensation for satellite loss of lock of 90 s.
PlaneBefore Compensation (m)After Compensation (m)
EJ1E0.0830.010
J1D0.0300.009
J1DN0.0260.009
J1L0.1230.010
J1M0.0780.008
J1R0.0440.016
SZSF0.0980.006
Table
Table 5. Compensation for satellite loss of lock of 180 s.
Table 5. Compensation for satellite loss of lock of 180 s.
PlaneBefore Compensation (m)After Compensation (m)
EJ1E0.0920.004
J1D0.0740.015
J1DN0.0690.015
J1L0.1820.005
J1M0.1060.005
J1R0.0620.014
SZSF0.0950.014
Table
Table 6. Compensation for satellite loss of lock of 90 s.
Table 6. Compensation for satellite loss of lock of 90 s.
Check the Number of FacesBefore Compensation (m)After Compensation (m)
70.0730.011
Table
Table 7. Compensation for satellite loss of lock of 180 s.
Table 7. Compensation for satellite loss of lock of 180 s.
Check the Number of FacesBefore Compensation (m)After Compensation (m)
70.1100.011

Share and Cite

MDPI and ACS Style

Shi, B.; Zhang, F.; Yang, F.; Lyu, Y.; Zhang, S.; Li, G. A Positioning Error Compensation Method for a Mobile Measurement System Based on Plane Control. Sensors 2020, 20, 294. https://doi.org/10.3390/s20010294

AMA Style

Shi B, Zhang F, Yang F, Lyu Y, Zhang S, Li G. A Positioning Error Compensation Method for a Mobile Measurement System Based on Plane Control. Sensors. 2020; 20(1):294. https://doi.org/10.3390/s20010294

Chicago/Turabian Style

Shi, Bo, Fan Zhang, Fanlin Yang, Yanquan Lyu, Shun Zhang, and Guoyu Li. 2020. "A Positioning Error Compensation Method for a Mobile Measurement System Based on Plane Control" Sensors 20, no. 1: 294. https://doi.org/10.3390/s20010294

APA Style

Shi, B., Zhang, F., Yang, F., Lyu, Y., Zhang, S., & Li, G. (2020). A Positioning Error Compensation Method for a Mobile Measurement System Based on Plane Control. Sensors, 20(1), 294. https://doi.org/10.3390/s20010294

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop