1. Introduction
Currently, many countries are vigorously developing intelligent connected vehicle-road cooperative technologies. High-precision dynamic positioning of highway vehicles serves as the foundation and prerequisite for achieving intelligent connected transportation.
In general road environments, vehicles can maintain good continuous positioning performance based on the combination of GNSS (Global Navigation Satellite System) and INS (Inertial Navigation System) [
1,
2,
3]. However, in strongly shielded environments such as long tunnels, effectively suppressing the divergence speed of INS is a important issue that urgently needs to be addressed [
4,
5]. When GNSS signals are unavailable because of the obstruction, multipath effects, and others, the GNSS/INS/OD(odometer) combination can effectively mitigate the accumulation of INS errors, to a certain extent, achieving high-precision stable navigation and positioning of the vehicle [
6,
7].
To simulate the positioning performance of tunnel vehicles using the GNSS/INS/OD combination, the authors conducted experiments on a certain section of the highway where all GNSS signals were artificially removed. Only inertial navigation and odometer were used for dead reckoning, and the GNSS/INS data fusion results were used as the reference trajectory (truth value) for comparison. The experimental results indicate that as the dead reckoning distance increases, the positioning results deviate more and more from the true trajectory. At a distance of 545 meters, the positioning error is approximately 0.88 meters. This cannot meet the practical requirements for tunnel vehicle positioning.
The reason behind this lies in the fact that under vehicle motion, the output of in-ertial devices contains significant noise. How to distinguish between valid data and noise, and reduce the impact of vibration noise on system stability and positioning accuracy, remains to be studied [
8,
9,
10]. In different scenarios, when the effectiveness of sensors in the integrated navigation system changes, corresponding adjustments need to be made to the state space model [
11,
12]. Multi-sensor systems are often heterogeneous due to their different principles, making it a challenge to construct reliable multi-sensor fusion function models that can accommodate real-time addition and removal of sensor types [
13,
14].
Based on this, this paper uses high-precision digital lane maps, as well as the ad-vantages of LiDAR in distance measurement and odometry in mileage counting, to propose a combined GNSS/INS/OD positioning method for tunnel vehicles based on lateral distance measurement and lane line constraints. It mainly consists of two steps: (1) lateral distance measurement of vehicles inside the tunnel based on laser radar point cloud data; (2) map matching combined positioning based on lateral distance measurement, odometer, and lane markings. The prerequisite for this positioning method is that a high-precision electronic map of the tunnel has been previously measured or obtained.
2. Research on Lateral Distance Measurement for Tunnel Vehicle Localization
The principle of a laser radar, also known as LiDAR (Light Detection and Ranging), is based on emitting an infrared laser beam from the radar itself and receiving the reflected infrared beam from the surrounding environment. By measuring the time difference between the emission and reception of the laser pulses, the distance can be calculated. Taking the example of a VLP-16 LiDAR (
Figure 1), it consists of 16 laser beams distributed within a 30-degree angle. The beams are evenly spaced at a 2-degree interval, with 15 degrees covering the upward and downward directions.
Using LiDAR for lateral distance measurement of tunnel vehicles, the specific method is as follows:
(1) Classify LiDAR points into surface feature points and line feature points based on curvature, and remove the line feature points.
By applying the LiDAR-measured point cloud data of the tunnel, calculate the curvature(
) of each point using the following formula:
In the equation, represents the set of five continuous points obtained by the same laser beam, and represent the coordinates of the j-th and i-th points in the k-th frame in the LiDAR coordinate system. Typically, the points with the highest curvature are selected as line feature points, while those with the lowest curvature are considered as surface feature points.
(2) Calculate the horizontal distance from the LiDAR to obstacles when LiDAR is horizontal, and select the m smallest values. These values correspond to the point cloud data that forms a set .
Assuming the maximum elevation angle of the LiDAR is
, select the point cloud data with an elevation angle of
obtained from step (1). Calculate the horizontal distance between the LiDAR and obstacles at different horizontal rotation angles when the elevation angle is
, using the following formula:
In the equation,
represents the horizontal rotation angle of the LiDAR,
represents the horizontal distance between the LiDAR and obstacles at an elevation angle of
and a horizontal rotation angle of
, and
represents the slope distance between the LiDAR and obstacles at an elevation angle of
and a horizontal rotation angle of
.
represents the set of point cloud data from the m LiDARs with the smallest
, and m represents the length of
. Refer to
Figure 2 for illustration.
(3) Calculate the horizontal distance from the LiDAR to obstacles when it is tilted.
Uneven road surfaces or tilting of the LiDAR may cause it to be non-horizontal, as shown in
Figure 3.
Select the point cloud data obtained from the LiDAR with an elevation angle of
from step (1). For each LiDAR point in
described in step (2), select the point cloud data with the same horizontal rotation angle and an elevation angle of
, then calculate the horizontal distance
between the LiDAR and obstacles at that horizontal rotation angle. The calculation formula is as follows:
In the equation, represents the elevation angle of the LiDAR, where . represents the distance between a LiDAR point with an elevation angle of and a horizontal rotation angle of and another LiDAR point with an elevation angle of and a horizontal rotation angle of . represents the slope distance between the LiDAR and obstacles measured at an elevation angle of and a horizontal rotation angle of . represents the horizontal distance between the LiDAR and obstacles measured at a horizontal rotation angle .
Figure 3.
Distance measurement when LiDAR is tilted.
Figure 3.
Distance measurement when LiDAR is tilted.
(4) Calculate the precise horizontal distance between the LiDAR and obstacles.
Calculate based on the parameter
described in step (3), using the following formula:
In the equation, represents the precise distance value between the LiDAR and obstacles, where is the minimum value among .
(5) Calculate the distance from the LiDAR to the other side of obstacles.
Select the distances
measured between the LiDAR and obstacles at an elevation angle of
and a horizontal rotation angle of
from step (1), as well as the distances
measured at an elevation angle of
and a horizontal rotation angle of
. Calculate the distance between the LiDAR and obstacles on the other side using the following formula:
In the equation, represents the distance between the LiDAR point with an elevation angle of and a horizontal rotation angle of and the LiDAR point with an elevation angle of and a horizontal rotation angle of , and represents the horizontal distance between the LiDAR and obstacles measured at a horizontal rotation angle .
(6) Correction of lateral distance.
In tunnels, emergency stopping zones are provided at regular intervals. Additionally, there are situations where large and high vehicles run parallel to the test vehicle. As shown in
Figure 4. Under these two special circumstances, the calculated lateral distance values are wrong and require correction. Let
and
represent the horizontal distances from the onboard LiDAR to the left and right obstacles, respectively, and
represents the typical width of a tunnel segment. In normal conditions, horizontal distances of the previous time
from the LiDAR to the left and right obstacles are represented by
and
, respectively. Considering that the width of the typical tunnel segment remains unchanged, the correction method for the lateral distance is as follows:
In the equation, within the judgment condition, an error threshold can be added, for example, by setting .
Figure 4.
Distance scheme for special situations.
Figure 4.
Distance scheme for special situations.
3. Combined Positioning of Tunnel Vehicles Constrained by Lane Lines and Lateral Distance
(1) Determine the starting point of mileage constraint. Before entering the tunnel, according to the GNSS RTK positioning results of the satellite, the corresponding lane line points are geometrically matched in the vertical direction, and are recorded as the starting point of the mileage constraint.
(2) Calculate accumulated vehicle mileage. According to the wheel speedometer, the cumulative calculation of vehicle mileage is accumulated by the encoder pulse count of the cumulative mileage meter. Assuming that the cumulative number of pulses of the odometer at a certain time is n, the total cumulative mileage at this time is:
Where,denotes the circumference of the wheel (prior value, unit: m), and represents the perimeter error coefficient in the process of vehicle travel, which is estimated by real-time joint solution calculation with GNSS RTK / INS.
Calculate the distance between the vehicle and the lane line. In the tunnel, assuming that the distance between the lane line and the tunnel wall remains unchanged (there will be a small change in reality) is
, as shown in
Figure 5, the distance between the vehicle and the lane line can be obtained as
Determine matching points. As shown in
Figure 5, the method of matching the corresponding points of lane lines according to the total mileage
is
Where, indicates the previous recent spot number of the match, and indicates the nearest point of the latter match; indicates the distance between lane line points (the value in this paper is fixed at 1m), and indicates the start dot number.
Calculate the coordinates of vertical points. According to the principle of distance division, the coordinates of vertical point
can be obtained as
Where,, indicates the hole down.
According to the lateral distance of the vehicle from the lane line and the azimuth of the point tangent line, the coordinates of the vehicle point
(
) are obtained, as shown in
Figure 6 and formula (14).
According to the above mileage matching and lane line constraint positioning method, higher lateral positioning accuracy can be obtained. However, in the longitudinal direction, due to the odometer calibration error, it is easy to produce mileage accumulation error when driving for a long time, so in the long-distance application environment, it is necessary to fully consider the impact of the odometer calibration error. Therefore, this paper proposes a positioning idea for further location combination of mileage matching positioning results and INS/OD. In other words, the above positioning results are further integrated with INS/OD combination, and the mileage matching positioning results are used to correct and feedback the INS zero bias error and odometer scale error. The odometer calibration error feedback calculated by INS/OD combination corrects the matching position of odometer and lane line, and realizes the mutual correction and complementary advantages between sensors and different positioning means.
4. Case Study of Lateral Distance Measurement
4.1. Static Testing
During the static test, the situation in the tunnel is simulated, the experimental environment is set in the corridor, the LiDAR is placed on the horizontal ground, the algorithm result value is measured by LiDAR, and the measured result value is measured by steel tape. The experimental results of multiple measurements are shown in
Table 1. As can be seen from the table, in this measurement, the result of the algorithm is in the centimeter level accuracy, and the fluctuation range is within 1 centimeter, and the lateral positioning results with centimeter-level accuracy in the tunnel can be obtained.
4.2. Dynamic Test
In a tunnel section, the vehicle lateral location measurement experiment is carried out, and the dynamic and practicability of the algorithm is verified by the actual data. The vehicle drove in the left lane, and the test obtained that the distance of the vehicle from the left wall was about 5.74 meters, the distance from the right wall was about 6.21 meters, and the total width of the tunnel was about 11.96 meters. According to experience, the design width of the tunnel is about 12 meters, so the accuracy of the distance measured from the LiDAR to the walls on both sides of the tunnel meets the requirements. The distance between the vehicle and the walls on both sides cannot be unchanged in the process of driving, so the distance between the vehicle and the walls on both sides fluctuates within a certain range.At the same time, the measured tunnel width is relatively stable, and the fluctuation is also in the centimeter-level range.
5. Case Study of Combined Positioning in Tunnels
The positioning performance in a real tunnel environment was verified by field tests in the Yudushan Tunnel in Beijing. The tunnel has a total length of 9.2 km, with open sections at the entrances and exits that can receive GNSS signals; however, the interior of the tunnel is a closed environment where GNSS data cannot be received. Vehicles travel at a speed of approximately 60 km/h inside the tunnel. The sensors used for the tests include GNSS RTK, INS, odometer, a laser radar for lateral distance measurement inside the tunnel, and lane map data. Translate GNSS RTK differential data to virtual reference station data provided by Qianxun Location FindCM service.
The tests involved driving from Beijing towards Zhangjiakou (designated as the forward direction), then returning in the opposite direction, making a total of two round trips with four passages through the tunnel. Three methods were used in each test:①GNSS/INS combination;②GNSS/INS/OD combination;③GNSS/INS/OD combination with lateral distance measurement and lane constraints.
(1)Comparison of Three Positioning Methods
Figure 7 shows the positioning trajectories of three positioning methods in the tunnel. The green points represent the GNSS/INS combined positioning results, the yellow points represent the GNSS/INS/OD combined positioning results, the red points represent the GNSS/INS/OD combined positioning results constrained by lateral distance and lane lines, and the blue points represent the actual lane lines. From the figure, it can be seen that the GNSS/INS combination deviates the most from the lane lines with the largest positioning error. The GNSS/INS/OD combination follows, although it has significantly improved compared to the GNSS/INS combination, the positioning deviation is still considerable over long distances. The GNSS/INS/OD combination with lateral distance lane line constraints provides the best positioning results.
Due to the moving status of vehicles in the tunnel and the lack of high-precision dynamic positioning assessment means on the roadside, it is difficult to quantitatively evaluate the absolute positioning accuracy in the tunnel. Therefore, in this study, in the open area outside the tunnel, the positioning data is compared with the RTK positioning results obtained at the same time as the first fixed solution to calculate the positioning accuracy when the vehicle exits the tunnel.
Figure 8,
Figure 9 and
Figure 10 show the positioning results outside the tunnel in the opposite direction for the first time. The yellow points represent the first fixed RTK results outside the tunnel, which serve as a reference due to the centimeter-level accuracy of RTK fixed solutions; the red points represent the positioning points extended from the tunnel positioning method (excluding the GNSS signals received when exiting the tunnel). The distances between corresponding points (at the same time) in the figures represent the error between the tunnel positioning points and the RTK reference points.
Table 2 provides the corresponding error statistics for two tunnel driving scenarios. It can be seen from the table that in the two sets of experiments, the positioning errors of GNSS/INS and GNSS/INS/OD combinations at the RTK fixed points outside the tunnel are over 100 meters, indicating that relying solely on INS, odometers, and other dead reckoning sensors in severe GNSS-denied environments in long tunnels it is difficult to maintain high-precision positioning and navigation capabilities. The GNSS/INS/OD combination with lateral distance and lane markings constraint provides positioning results closest to the RTK results, with a horizontal positioning accuracy between 1.2m and 1.5m.
(2)Analysis of Lateral Positioning Error
In order to quantitatively analyze the lateral accuracy of GNSS/INS/OD combined positioning constrained by lane lines inside the tunnel, the experiment results in the opposite direction are analyzed and statistically processed below. In this experiment, the vehicle's driving speed inside the tunnel is approximately around 60 km/h, as shown in
Figure 11.
Figure 12 provides the lateral error of a 9.2 km tunnel, with the elliptical marked portion representing a 4.6 km half-length tunnel. The calculation method of lateral error is as follows: at each comparison moment, calculate the vertical distance between the positioning point and the lane line, then subtract the actual distance of the vehicle from the lane line at the corresponding moment (this distance is obtained by subtracting the lane line distance from the tunnel wall distance measured by LiDAR).
As can be seen from
Figure 12, with the increase of driving distance in the tunnel, the positioning error shows a gradually diverging trend. Especially in the second half of the tunnel, the positioning error diverges severely, reaching a maximum of 1m. The positioning error results for the 4.6km midpoint of the tunnel are summarized in
Table 3. The table includes the arithmetic mean, absolute error mean, root mean square error, maximum error, minimum error, and the 95th percentile of the absolute error. In
Table 3, a total of 61,001 positioning points were analyzed. It can be observed that the arithmetic mean of the lateral error is 0.013m, and the absolute error mean is 0.294m, which characterizes the absolute lateral accuracy of vehicle positioning. The maximum error and minimum error are 0.671m and -0.865m, respectively, with the 95th percentile of the absolute error being 0.612m, indicating that further improvement is needed in the stability of vehicle positioning in the tunnel.
(3)Analysis of Longitudinal Positioning Errors
For the longitudinal positioning errors inside the tunnel, due to the lack of reference points for dynamic positioning evaluation within the tunnel, the longitudinal accuracy at the 9.2 km tunnel entrance based on RTK positioning results is used to infer the longitudinal accuracy over the 4.6 km tunnel length. Similarly, the analysis is conducted using experimental results in the reverse direction. When the vehicle first obtains a fixed RTK solution upon exiting the tunnel, the comparison between the positioning points extended by the tunnel's positioning method and the RTK positioning points is shown in
Figure 10. At this point, the total horizontal positioning error is 1.43m, and after subtracting the lateral positioning error of that point (0.86m), the longitudinal positioning error is calculated to be 1.14m. Based on the divergence characteristics of INS based on navigation estimation sensors and the cumulative nature of odometer errors with distance, it is known that the positioning error diverges faster with increasing driving distance, and the divergence speed in the latter half of the journey is significantly faster than in the first half, as indicated in
Figure 12, which also shows a similar phenomenon in lateral positioning results. Therefore, it can be conservatively estimated that the longitudinal positioning error at the 4.6 km mark of the tunnel is not greater than 0.6m. It should be noted that this longitudinal error is only an inferred value based on the results at the tunnel exit, representing a general accuracy level. The positioning errors at different positions within the tunnel may still fluctuate due to various measurement errors, as illustrated in
Figure 12.
6. Conclusions
Building upon the GNSS/INS/OD combination positioning, this paper proposes a lane-constrained positioning method integrating odometer mileage calculation and laser radar lateral distance measurement. The constrained positioning results are used to update the measurements, achieving fusion between lane-constrained positioning within the tunnel and INS/OD combination positioning. First, the paper introduces a lateral distance measurement method inside tunnels based on laser radar. By constraining and extracting point clouds, the closest distances from the laser radar to the walls on both sides of the tunnel are calculated. Based on the horizontal rotation angle of these points, the left and right attributes of the walls are further determined. Additionally, by assuming a constant total width of the tunnel, it is solved for interference of dynamic vehicles, emergency parking spaces, emergency exits and other structural change scenarios on lateral distance measurement, enabling continuous lateral distance measurement. Furthermore, the positioning points calculated from odometer mileage and lateral distance are fused with the INS/OD combination positioning. On one hand, this approach replaces GNSS measurement updates with positioning measurement points, effectively suppressing the divergence of the INS/OD errors. Simultaneously, the real-time calibration errors estimated by the odometer are fed back into the odometer mileage calculation, enhancing the accuracy of lane-constrained positioning. Through the fusion of these two positioning models, errors between different positioning models are mutually compensated and suppressed, thereby enhancing the accuracy and robustness of tunnel positioning.
The vehicle test results conducted at speeds around 60 km/h inside tunnels demonstrate that the lane-constrained GNSS/INS/OD combination positioning model significantly improves the positioning accuracy of vehicles in long tunnels compared to GNSS/INS combination and GNSS/INS/OD combination alone. Through statistical analysis of 61,001 positioning points collected along a 4.6 km tunnel, the average absolute error of lateral positioning is 0.294 meters. In terms of longitudinal accuracy, comparing with the precision at the tunnel exit (approximately 9.2 km) and RTK fixed solution, the calculated longitudinal error at the midpoint of the 4.6 km tunnel is less than 0.6 meters.
Author Contributions
Formal analysis, Hongbin Zhang; Methodology, Hongbin Zhang; Writing – original draft, Hongbin Zhang. Writing-review & editing, Xu Zhang.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
No new data were created.
Acknowledgments
Acknowledge the tunnel administration for providing the test environment.
Conflicts of Interest
The authors declare no conflicts of interest.
References
- Zhang, H.; Xia, X.; Nitsch, M.; Abel, D. Continuous-Time Factor Graph Optimization for Trajectory Smoothness of GNSS/INS Navigation in Temporarily GNSS-Denied Environments. IEEE Robotics and Automation Letters. 2022, 7(4):9115-9122. [CrossRef]
- He, Y.Li, J.Liu, J. Research on GNSS INS & GNSS/INS Integrated Navigation Method for Autonomous Vehicles: A Survey, EEE Access Access, IEEE 2023,11:79033-79055. [CrossRef]
- Chen, Z.; Li, Z.; Chen, W.; Sun, Y.; Ding, X.. Wasserstein Distance-Assisted Variational Bayesian Kalman Filter With Statistical Similarity Measure for GNSS/INS Integrated Navigation. IEEE Sensors Journal. 2024(6):8807-8820. [CrossRef]
- Gao, Wei.Zhan, Xingqun.Yang, Rong INS-aiding information error modeling in GNSS/INS ultra-tight integration, GPS Solutions: The Journal of Global Navigation Satellite Systems. 2024, 28(1). [CrossRef]
- SHANGGUAN Wei; CHEN Jingjing; XIE Chaoxi; JIANG Wei. Magnetometer-enhanced Combination Positioning Method Based on IMU Error Compensation. Journal of the China Railway Society. 2022, 44(7):80-90.
- Ke, Y.; Lv, Z.; Zhang, C.; Deng, X.; Zhou, W.; Song, D. Tightly Coupled GNSS/INS Integration Spoofing Detection Algorithm Based on Innovation Rate Optimization and Robust Estimation. IEEE Access. 2022,10:72444-72457. [CrossRef]
- ZHU Zhoulin; LU Yishi; XIE Zhilong. GNSS/INS/Wheel Speed Integrated Positioning Post-Processing Algorithm Based on Iterative Smoothing. Automobile Parts. 2022(z1):41-47.
- QIN Honglei; DU Yansong. An integrated positioning algorithm of Iridium opportunity signals and MEMS-INS. Journal of Navigation and Positioning. 2023, 11(3):45-52.
- Wang, Junwei; Chen, Xiyuan; Shi, Chunfeng; Liu, Jianguo. An Improved Robust Estimation Method for GNSS/SINS under GNSS-Challenged Environment. 2022 11th International Conference on Control, Automation and Information Sciences (ICCAIS):79-83 Nov, 2022.
- Lu, W.; Teng, Y.; Yan, X.; Jia, X.; Du, Y. Improved adaptive SRCKF algorithm for GNSS/SINS integrated navigation based on measurement characteristics. Journal of Chinese Inertial Technology, April 2023, 31(4):327-334.
- Chen, Wei; Yang, Gongliu; Tu, Yongqiang. A Digital Track Map-Assisted SINS/OD Fusion Algorithm for Onboard Train Localization. Applied Sciences. December, 2023, Vol. 14 Issue 1; MDPI AG. [CrossRef]
- Chen, He; Zhang, Zhili; Zhou, Zhaofa; Liu, Pengpeng; Guo, Qi. SINS/ODIntegrated Navigation Algorithm Based on Body Frame Position Increment for Land Vehicles. Mathematical Problems in Engineering. 8/13/2018, p1-11. [CrossRef]
- YF,wang. Research on multi-frequency GNSS/INS/OD fusion high-precision positioning technology in urban environment[D]. Southeast University, 2021.
- LIU Wanke; NONG Qi; TAO Xianlu; ZHU Feng; HU Jie. OD/SINS adaptive integrated navigation method with non-holonomic constraints. Acta Geodaetica et Cartographica Sinica, Vol 51, Iss 1, p9-17, 2022. [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).