Next Article in Journal
Fine-Scale Improved Carbon Bookkeeping Model Using Landsat Time Series for Subtropical Forest, Southern China
Previous Article in Journal
Remote Sensing Monitoring of Grasslands Based on Adaptive Feature Fusion with Multi-Source Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A GRU and AKF-Based Hybrid Algorithm for Improving INS/GNSS Navigation Accuracy during GNSS Outage

1
GNSS Research Center, Wuhan University, Wuhan 430079, China
2
School of Electronics and Information Engineering, Hubei University of Science and Technology, Xianning 437100, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(3), 752; https://doi.org/10.3390/rs14030752
Submission received: 10 January 2022 / Revised: 28 January 2022 / Accepted: 29 January 2022 / Published: 6 February 2022
(This article belongs to the Topic GNSS Measurement Technique in Aerial Navigation)

Abstract

:
The integrated navigation system consisting of an inertial navigation system (INS) and Global Navigation Satellite System (GNSS) provides continuous high-accuracy positioning whereas the navigation accuracy during a GNSS outage inevitably degrades owing to INS error divergence. To reduce such degradation, a gated recurrent unit (GRU) and adaptive Kalman filter (AKF)-based hybrid algorithm is proposed. The GRU network, which has advantages of high accuracy and efficiency, is constructed to predict the position variations during GNSS outage. Furthermore, this paper takes the GRU-predicted error accumulation into consideration, and introduces AKF as a supplementary methodology to improve the navigation performance. The proposed hybrid algorithm is trained and tested by practical road datasets and compared with four algorithms, including the standard KF, Multi-Layer Perceptron (MLP)-aided KF, Long Short Time Memory (LSTM) aided KF, and GRU-aided KF. Periods of 180 and 120 s GNSS outage are employed to test the performance of the proposed algorithm in different time scales. The comparison result between the standard KF and neural network-aided KF indicates that the neural network is an effective methodology for bridging GNSS outages. The performance comparison between three kinds of neural networks demonstrate that both recurrent neural networks surpass the MLP in prediction position variation, and the GRU transcends the LSTM in prediction accuracy and training efficiency. Furthermore, it is concluded that the adaptive estimation theory is an effective complement to neural network-aided navigation, as the GRU-aided AKF reduced the horizontal error of GRU-aided KF by 31.71% and 16.12% after 180 and 120 s of GNSS outage, respectively.

Graphical Abstract

1. Introduction

INS and GNSS are two of the most widely used navigation techniques in both civilian and military fields. GNSS provides high accuracy position and velocity information with a relatively stable noise level in open-sky outdoor environments [1,2]. Nevertheless, it suffers from the shortcoming of signal vulnerability, which leads to accuracy degradation in complex urban environment, including overpasses, boulevards, and urban canyons, etc. On the other hand, INS is a self-contained navigation system that estimates the position, velocity, and attitude with a high update frequency. Although, INS suffers from a drawback in that its error accumulates over time [3]. Since INS and GNSS have complementary characteristics, they are combined as an INS/GNSS integrated navigation system which surpasses both stand-alone systems [4,5,6,7].
Although INS error could be estimated and compensated for by the KF-based INS/GNSS integration algorithm [8], low-cost INS error accumulates rapidly [9] during GNSS outage. Generally, there are, in general, two categories of techniques used to enhance low-cost INS /GNSS integrated navigation performance during GNSS outage. The first one is adopting auxiliary sensors, such as a camera [10,11,12,13], Lidar [14,15], magnetic sensor [16], or odometer [17,18], to collect various information to aid INS. The drawbacks of this kind of solution include increases in hardware cost, equipment volume, and power consumption. Another category of techniques is suppressing INS error accumulation by algorithm improvement without extra sensors. For instance, non-holonomic constraint (NHC) [17,19], assuming zero velocity at the land vehicle laterally and vertically, adopts pseudo velocity measurement to suppress INS velocity error divergence, and zero velocity update (ZUPT) [20] corrects INS error when the navigation system is stationary. Recently, the artificial intelligent (AI) algorithm has been employed in the navigation field due to its advantages of solving the non-linear problem, and AI-aided navigation systems have shown impressive performance. However, the AI algorithm also has the disadvantages of a long training time and the predicted value inevitably containing error. The former disadvantage could be partly solved by simplifying the neural network structure while the latter one could be suppressed using the adaptive algorithm. Thus, this paper constructed a simple and effective GRU network to aid INS during GNSS outage, and adopts the AKF algorithm to estimate pseudo position measurement noise, thus optimizing the integrated navigation accuracy.
The novelties of this paper are summarized as follows: (1)A GRU-based pseudo GNSS position variation prediction neural network is proposed in this work. The proposed network is verified to have the advantages of a high prediction accuracy and high training efficiency. (2) Meanwhile, different from existing related AI methodology, this work takes the AI prediction error into consideration. An adaptive filtering strategy is adopted as the supplement of the GRU network, and the Sage-Husa AKF is introduced to estimate the state vector and pseudo GNSS position noise simultaneously. The rationale behind this methodology is that assigning proper weight to the GRU-predicted information during the INS and pseudo GNSS data fusion process improves the navigation performance.
The structure of the remaining part is organized as follow. The rationale of the GRU and AKF hybrid algorithm is presented in Section 2. Then, the detailed GRU network implement method is described in Section 3, followed by the test and result analysis in Section 4. Finally, the discussion and conclusions are put forward in Section 5 and Section 6, respectively.

2. Rationale of the Proposed Algorithm

2.1. AI Module Input and Output Parameters

It is crucial to select proper input and output parameters for the AI module, since these parameters directly determine the navigation accuracy and training efficiency. All existing AI models applied in bridging GNSS outage can be divided into three categories. The first kind is the OINS-δPGNSS, INS model [21,22,23,24,25], which predicts δPGNSS, INS standing for the position error between GNSS measurement and INS estimation with the input parameters OINS (including the position, velocity, attitude estimated by INS, and angular velocity, specific force measured by the inertial measurement unit). Another category is OINS-X model [26], which estimates the state vector X of the integrated navigation system directly. The third type of model OINS-ΔPGNSS [27,28,29], which calculates the GNSS position variation ΔPGNSS, is adopted in this paper. The GNSS position variation error is expressed in Equation (1), where Δ   P ^ GNSS ( t k ) is the estimated position variance between epochs t k and t k 1 , and P ^ GNSS ( t k ) and   P ^ GNSS ( t k 1 ) represent the estimated position at t k and t k 1 , respectively. The estimated position   P ^ GNSS ( t k ) consists of the truth value P GNSS t k and the error term δ   P GNSS t k .     P ^ GNSS ( t k 1 ) similarly consists of the truth value P GNSS t k 1 and the error term δ   P GNSS t k 1 . It can be deduced from the formula that differential arithmetic between adjacent epochs could remove the GNSS measurement noise, which has time correlativity. In other words, the OINS-ΔPGNSS model has lower prediction error than the former two models:
Δ P ^ GNSS ( t k ) = P ^ GNSS ( t k ) P ^ GNSS ( t k 1 ) = ( P GNSS t k + δ   P GNSS t k ) P GNSS t k 1 + δ   P GNSS t k 1 = Δ P GNSS t k δ   P GNSS t k δ   P GNSS t k 1
The position variation Δ P GNSS ( t k )   can also be expressed as the difference in the INS derived position between adjacent epochs as shown in Equation (2):
Δ P GNSS = C b n f ib b t 2 ω ie n t + ω en n t × V n t + G n dtdt
where b, i, e, and n stand for the body frame (b-frame), inertial frame (i-frame), earth frame (e-frame), and navigation frame (n-frame), respectively. C b n is the directional cosine matrix from the n-frame to b-frame, f ib b is the specific force acquired by the inertial measurement unit (IMU), ω ie n is the angular rate of the e-frame relative to the i-frame, ω en n is the angular rate of the n-frame relative to the e-frame, V n is the INS-derived velocity, and G n is the gravity vector. Since ω ie n ,   ω en n , and G n are related to the longitude and latitude, which vary slowly in land vehicles, these three parameters can be neglected to improve the training efficiency. The input parameters of the AI model are the IMU data, including f ib b and ω ib b and the INS-estimated velocity ( V n ) and attitude ( A INS ) .

2.2. Neural Network Model: GRU

Various neural network models are adopted in AI-aided integrated navigation. Assad adopted the radial basis function neural network (RBF) to improve the attitude estimation accuracy in GPS-denied environments [30]. Yao estimated the pseudo GPS position with the MLP network when GPS signal was unavailable [28]. However, these methods employ a static neural network, which only deals with current data but neglects the information of historical data. Neglecting the historical data may result in a low prediction accuracy. Thus, the recurrent neural network that is widely used in time series signals processing is a better choice for integrated navigation. Fang suggested a method for bridging GNSS outage that used a kind of recurrent neural network called the LSTM, and obtained a satisfying effect [29]. The GRU, which is also a recurrent neural network, with the same estimation accuracy and higher training efficiency compared with LSTM [31], is adopted in this paper. As shown in Figure 1, the GRU unit consists of an update gate and a reset gate. The output of the update gate (zt) controls the extent of the impact on the current state forced by previous states; that is, the smaller zt is, the less previous information is retained. The function of the reset gate(rt) is to determine the forgetfulness degree of the hidden state information h t ˜ ; that is, the smaller the value is, the more past state information is discarded.
The GRU forward propagation formula is shown in Equation (3):
r t = σ ( W xr x t + W hr h t 1 + b r ) z t = σ ( W xz x t + W hz h t 1 + b z ) h ˜ t h t = z t h t 1 + ( 1 z t )   h ˜ t = tanh ( W xh x t + W hh ( r t h t 1 ) + b h ) h t = z t h t 1 + ( 1 z t )   h ˜ t
where, W represents the weight matrix between different structures; subscripts x, r, h, and z denote the input, reset gate, hidden state, and update gate; and b r , b z , and b h   are the offset vectors of the corresponding structure.

2.3. Proposed GRU-aided AKF Algorithm

The proposed INS/GNSS integrated navigation algorithm employs an AKF-based loosely coupled integration strategy with a GRU module, which works in training mode when GNSS signal is available and prediction mode during a GNSS outage.
According to the above analysis, the proposed scheme consisting of the training mode and predicting mode is illustrated in Figure 2. As shown in Figure 2a, the GRU module operates in training mode when GNSS signal is available. With the input parameters of f ib b , ω ib b , V n , and A INS , the GRU module outputs the predicted position variance ΔPGNSS, which is compared with the high precision ΔPtrue representing the position variation derived from the real GNSS signal. The GRU module adjusts the internal weight continuously until the loss function drops below the threshold. Meanwhile, the KF module works as a traditional loosely coupled integrated navigation system, providing high accuracy position, velocity, attitude, and IMU error information in real time. When GNSS signal is unavailable, the GRU module operates in the predicting mode shown in Figure 2b. The GRU module predicts the position variation ΔPGNSS, which is integrated as a pseudo GNSS position measurement. Because ΔPGNSS inevitably includes error, the pseudo position measurement also includes error, which increases during the ΔPGNSS integral process. However, all existing literature [23,24,25] that adopted the OINS-ΔPGNSS model did not take this phenomenon into consideration.
In this paper, the 15-dimentional state vector X is defined as X = [δr, δv, φ, bg, ba] T, where δr and δv are the 3-dimensional INS position error and velocity error in the n-frame, φ stands for the misalignment angles, and bg and ba are the gyroscope and accelerometer bias in the x axis, y axis, and z axis. The three-dimensional measurement vector Z is defined as the difference between the INS-estimated position and the GRU-estimated position. The discrete-time Kalman state function and measurement function are defined as Equation (4):
X k = Φ k | k 1 X k 1 + Γ k 1 ω k 1 Z k = H k X k + V k
where Xk is the state vector at time k, Φk|k−1 is the one-step state transition matrix from time k − 1 to k, Γk−1 is the system noise-input mapping matrix, ωk−1 is the system noise vector presenting a zero mean Gauss distribution with variance Qk−1, Z is the measurement vector, Hk represents the measurement matrix, and Vk is the measurement noise vector subjected to a zero mean Gauss distribution with variance Rk−1. Theoretically, KF could acquire an optimal estimation when noise statistical parameters, such as Qk−1 and Rk−1, are precisely known [32]. However, the ΔPGNSS predicted by the GRU module inevitably contains error, which is accumulated when ΔPGNSS is integral to PGNSS. All the existing literatures utilizing the OINS-ΔPGNSS model [27,28,29] did not take this phenomenon into consideration.
Standard KF includes the time update and measurement update process. The time update process is formulized as Equation (5):
X ^ k | k 1 = Φ k | k 1   X ^ k 1 P k | k 1 = Φ k | k 1 P k 1 Φ k | k 1 T + Γ k 1 Q k 1 Γ k 1 T
where   X ^ k | k 1 is the prediction of the state vector at time k,   X ^ k 1 is the optimal estimation at time k − 1, P k | k 1 is the prediction covariance matrix at time k, and P k 1 is the optimal covariance matrix estimation at k − 1.
The KF measurement update process is illustrated in Equation (6):
K k = P k | k 1 H k T H k P k | k 1 H k T + R k 1 X ^ k = X ^ k | k 1 + K k ( Z k H k X ^ k | k 1 ) P k = ( I K k H k ) P k | k 1 I K k H k T + K k R k R k T
where K k is the KF gain, and   X ^ k and P k are the optimal estimated state and its covariance matrix, respectively.
In order to furtherly improve the accuracy, the pseudo noise increase (namely the time-varying R) of the OINS-ΔPGNSS model should be taken into consideration. The Sage-Husa filter, which was proposed with a standard KF framework with the capability of estimating the system state and noise parameters simultaneously, is a solution for the increasing pseudo measurement noise. Its essence is matching the theoretical covariance of the innovation or residual sequence with its practical equivalent [33]. The former methodology is also known as the IAE (innovation-based adaptive estimation) filter and the latter one as the RAE (residual-based adaptive estimation) filter. The RAE can only estimate noise covariance with residuals prior to the current epoch, and once the current epoch measurement precision is inconsistent with the historical precision, it is hard for this methodology to guarantee estimation reliability [34]. Therefore, the IAE filter is adopted in this paper.
The innovation is defined as Equation (7) shows:
  Z ˜ k | k 1 = Z k   Z ^ k | k 1 = H k X k + V k H k   X ^ k | k 1 = H k   X ˜ k | k 1 + V k
The state prediction error   X ˜ k | k 1   is unbiased as the initial value is unbiased, and the mean of the measured noise Vk is assumed to be zero. It can be concluded that the innovation is zero mean. Since   X ˜ k | k 1 and Vk are unrelated, the variance of   Z ˜ k | k 1 can be expressed as Equation (8):
E Z ˜ k / k 1   Z ˜ k / k 1 T = H k P k / k 1 H k T + R k
The measurement noise matrix Rk is presented as:
R k = E Z ˜ k / k 1 Z ˜ k / k 1 T H k P k / k 1 H k T
In theory, E   Z ˜ k / k 1   Z ˜ k / k 1 T can be replaced by the innovation average value over time, and the R k equally weighted recursive estimation equation can be constructed as Equation (10):
R ^ k = 1 k i = 1 k   Z ˜ i | i 1   Z ˜ i | i 1 T H i P i | i 1 H i T = 1 k k 1   R ^ k 1 +   Z ˜ k | k 1   Z ˜ k | k 1 T H k P k | k 1 H k T = 1 1 k   R ^ k 1 + 1 k   Z ˜ k | k 1   Z ˜ k | k 1 T H k P k | k 1 H k T
The Sage-Husa AKF is fully constructed while the   R ^ k   is substituted into Equation (6), and the filter can estimate the noise matrix in time and adjust the KF gain.

3. Implementation of the Proposed Algorithm

3.1. Training and Testing Dataset Description

In order to train and verify the performance of the proposed algorithm, a dataset was acquired with equipment mounted on a land vehicle shown in Figure 3. The data collection equipment consisted of micro-electromechanical systems (MEMSs) IMU whose crucial specifications are summarized in Table 1 and a single-frequency GNSS receiver chip. Meanwhile, a set of reference data was gathered with NovAtel SPAN CPT.experimental platform.
The dataset was collected in Wuhan City, Hubei Province, China, and the dataset collection trajectory is shown in Figure 4. The total length of the dataset was 5200 s (GPS time: 198,900–204,100). Among them, data from GPS time 198,900 to 201,900 s were utilized as the GRU network training dataset, and data from GPS time 201,900 to 204,100 s were used as the testing dataset.

3.2. Training Process

The GRU module was trained with 12 input parameters, including the 3-axis gyroscope data, 3-axis accelerometer data, and 3-axis velocity and attitude data. Meanwhile, the GNSS position variation ΔPGNSS was taken as the training label. To balance the training efficiency and predicting accuracy, the GRU model comprised a GRU layer and a fully connected layer activated by the sigmoid function. Furthermore, two significant hyperparameters were utilized, namely the number of neurons and the step time. In order to accelerate the training convergence rate and prevent overfitting, experiments were carried out for the adjustment of both hyperparameters, and the prediction performance is shown in Table 2. Four kinds of time step (1, 2, 4, 8) and hidden neurons of GRU (32, 64, 128, 256) were evaluated. As is shown in Table 2, the training time of the network increases with the increase of both hyperparameters, whereas the prediction accuracy is not positively correlated with them. The position accuracy predicted by 128 neurons surpasses that predicted by 32, 64, and 256 neurons. An excessive neuron number could result in the overfitting phenomenon and degradation of the generalization performance while an insufficient neuron number cannot express the relationship between OINS and ΔPGNSS adequately. After setting the hidden neurons as 128, it can be seen in Table 2 that the prediction accuracy with a time step of 4 surpasses the other two groups. Therefore, the neural number and time step were set as 128 and 4, respectively.

4. Test and Result Analysis

4.1. GRU Prediction Accuracy

The GRU network prediction accuracy was verified with the 2200 s (GPS time: 201,900–204,100 s) test dataset. Meanwhile, an LSTM network and MLP network were constructed as comparison groups. The LSTM network consisted of an LSTM layer with 128 neurons and a fully connected layer activated by the sigmoid function, and the time step was set as 4. The MLP network consisted of 128 neurons but without the time step parameter. The truth value was acquired by SPAN-CPT, and the predicted position variation root mean square error (RMSE) values of the different networks are summarized in Table 3. The prediction RMSE values of the GRU, LSTM, and MLP network in the North direction are 3.05 × 10−8, 3.24 × 10−8, and 6.49 × 10−7, and the RMSE values in the East direction are 3.04 × 10−8, 4.96 × 10−8, and 6.81 × 10−7. It can be concluded that both recurrent networks (LSTM and GRU) have a similar prediction accuracy and visibly surpass the performance of the MLP network. On the other hand, the training time consumption per epoch of GRU and LSTM is 1.90 and 2.19 s and this result verifies that the GRU network surpasses LSTM in efficiency. The comparison results between LSTM and GRU are also consistent with the conclusion of [27].
The North and East position variations in the GRU prediction, LSTM prediction, and MLP prediction and truth value are shown in Figure 5a,b, respectively. As Figure 5 indicates, the predicted value of GRU and LSTM identifies the truth value while the MLP-predicted value obviously differs from the truth value.
Figure 6 shows the GRU-, LSTM-, and MLP-predicted horizontal trajectory, which is an integral value of the predicted position variation. More specifically, the below trajectory reflects the change in the position of each epoch relative to the starting point (GPS time: 201,900 s, latitude: 30.44681283°, longitude: 114.461966°). As shown in Figure 6, the starting points of the predicted curves are consistent with that of the reference curve, and then the prediction position error accumulates gradually. The reason for this phenomenon is that the predicted position variation, namely ΔPGNSS, inevitably includes error, and the prediction error, which is positively related to the measurement noise, can accumulate while the ΔPGNSS is integrated into pseudo PGNSS, thus it is necessary to adopt AKF to estimate the pseudo GNSS measurement noise matrix.

4.2. Integrated Navigation Accuracy

In order to verify the performance of the GRU-aided AKF algorithm, a period of 180 s (GPS time:202,100–202,280 s) and a period of 120 s (GPS time: 202,300–202,420) were selected from the test dataset to simulate a GNSS outage. The position information derived from the proposed GRU network replaces the real GNSS information in the measurement update process of AKF and KF. Meanwhile, the LSTM-aided KF, MLP-aided KF, and not-aided KF were also adopted as control groups.
As is shown in Figure 7, the MLP, LSTM, and GRU network could suppress the velocity error accumulation in both the North and East directions. The maximum velocity errors in the North and East directions estimated by the above-mentioned algorithm are summarized in Table 4. Compared with the not-aided KF, MLP-aided and LSTM-aided KF, and GRU-aided KF methods, the proposed GRU-aided AKF reduced the maximum North velocity error by 7.071%, 6.105%, 3.922%, and 5.431%, and the maximum east velocity error by 36.686%, 34.896%, 26.87%, and 25.595%.
Figure 8 shows a comparison of the North, East, and horizontal errors of different algorithms during the 180 s GNSS outage. It is illustrated in Figure 8 that the LSTM-aided KF and GRU-aided KF have almost the same position error during the 180 s GNSS outage, and the GRU network is a better choice due to its higher efficiency. Meanwhile, the GRU-aided AKF’s prediction accuracy exceeds that of the GRU-aided KF after the 120 s GNSS outage, and it can be concluded that the proposed GRU-aided AKF is better suited to a long period of GNSS outage.
The maximum position errors and position error RMS during the 180 s GNSS outage are summarized in Table 5. The quantitative analysis indicates that the proposed GRU-aided AKF algorithm reduced the maximal horizontal error by 76.03%, 70.10%, 32.30%, and 31.76%, and the horizontal error RMS by 83.03%, 76.55%, 33.40%, and 31.99% compared with the Not-aided, MLP-aided, LSTM-aided, and GRU-aided KF, respectively. Meanwhile, by comparing the data of the MLP-aided, LSTM-aided, and GRU-aided KF, it can be concluded that both recurrent neural networks, namely the LSTM and the GRU, have a similar positioning accuracy and surpass the performance of MLP.
Figure 9 shows a comparison of the North, East, and horizontal errors of the different algorithms during the 120 s GNSS outage. It is obvious that all of the AI-aided algorithms surpass the no-aided one and all of the recurrent neural network-aided algorithms surpassed the MLP-aided one regarding the positioning accuracy. Besides, it can be observed that the proposed GRU-aided AKF position error curve is lower than that of the control groups.
The maximum position errors and position error RMS during the 120 s GNSS outage are summarized in Table 6. The proposed GRU-aided AKF algorithm reduced the maximal horizontal error by 82.65%, 50.36%, 12.68%, and 16.12% and the horizontal error RMS by 75.39%, 42.17%, 8.56%, and 12.22% of the not-aided, MLP-aided, LSTM-aided, and GRU-aided KF, respectively. By comparing the not-aided KF and the proposed GRU-aided AKF horizontal error RMS in Table 5 and Table 6, it can observed that the proposed algorithm reduced the horizon error RMS by 75.39% and 83.03% under the 120 and 180 s GNSS outage respectively. This indicates that the proposed algorithm performs better during longer periods of GNSS outage.

5. Discussion

The main advantages of the GRU and AKF hybrid algorithm can be summarized as follows: (1) A high efficiency recurrent neural cell with a simple structure, namely the GRU, was chosen to estimate the pseudo position variation. Benefiting from proper neural cell selection and hyperparameter setting, the proposed algorithm possesses superior navigation accuracy and training efficiency. (2) This paper combined the modern artificial intelligence methodology and conventional adaptation theory to improve the navigation accuracy during GNSS outage. The experiment results under different time scales indicate that the hybrid algorithm surpasses the stand-along adaptive and AI methods.
Meanwhile, further improvements of this work are summarized as follows: (1) The algorithm is currently run on PC and cannot satisfy the real-time dynamic navigation requirement. A potential solution is to integrate the AI chip and the embedded platform. The integrated computational platform could support the proposed method, and has lower power and less volume advantages compared with PC. (2) The training and testing dataset could not cover all the application scenarios, and the GRU network should be trained with more different types of data. As the amount of training data grows, the proposed algorithm will acquire better applicability and adaptability.

6. Conclusions

Aiming to improve the positioning accuracy of INS/GNSS integrated navigation during GNSS outage, this study proposed a GRU and AKF-based integrated navigation algorithm. Owing to proper neural cell selection and hyperparameter setting, the constructed GRU network has a higher prediction accuracy and higher training efficiency. Furthermore, the proposed algorithm introduces adaptive KF to reduce the position error by estimating the GRU prediction error dynamically. Compared with the standard KF-based integrated navigation algorithm, the test results show that the hybrid algorithm reduces the root mean square error by 83.03% and 75.39% and the maximum position error by 76.03% and 82.65% during the 180 and 120 s GNSS outage, respectively. The algorithm can improve the navigation accuracy at different time scales, and has real-time navigation application potential due to its high efficiency.

Author Contributions

Conceptualization, J.J. and Y.T. (Yanan Tang); methodology, J.J., Y.T. (Yanan Tang) and Y.T. (Yifeng Tao); software, Y.T. (Yanan Tang); validation, Y.T. (Yanan Tang) and Y.T. (Yifeng Tao); formal analysis, Y.T. (Yanan Tang); investigation, J.J. and Y.T. (Yanan Tang); resources, Y.T. (Yanan Tang) and Y.T. (Yifeng Tao); writing—original draft preparation, J.J. and Y.T. (Yanan Tang); writing—review and editing, Y.T. (Yanan Tang), J.L. (Jianghua Liu), P.Y. and J.L. (Jingnan Liu); visualization, J.L. (Jianghua Liu); supervision, J.J. and Y.T. (Yanan Tang); project administration, J.L. (Jianghua Liu); funding acquisition, J.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2018YFB0505200 and 2018YFB0505201), the Fundamental Research Funds for the Central Universities (2042018kf0253).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Thanks for the financial support of the National Key Research and Development Program of China (2018YFB0505200 and 2018YFB0505201). Many thanks for professor Jingnan Liu, Professor Jinguang Jiang and his team.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Godha, S.; Cannon, M.E. GPS/MEMS INS Integrated System for Navigation in Urban Areas. GPS Solut. 2007, 11, 193–203. [Google Scholar] [CrossRef]
  2. Xiao, Y.; Luo, H.; Zhao, F.; Wu, F.; Gao, X.; Wang, Q.; Cui, L. Residual Attention Network-Based Confidence Estimation Algorithm for Non-Holonomic Constraint in GNSS/INS Integrated Navigation System. IEEE Trans. Veh. Technol. 2021, 70, 11404–11418. [Google Scholar] [CrossRef]
  3. Grewal, M.S.; Weill, L.R.; Andrews, A.P. Global Positioning Systems, Inertial Navigation, and Integration; John Wiley: Hoboken, NJ, USA, 2001. [Google Scholar]
  4. Jiang, H.; Shi, C.; Li, T.; Dong, Y.; Li, Y.; Jing, G. Low-cost GPS/INS Integration with Accurate Measurement Modeling Using an Extended State Observer. GPS Solut. J. Glob. Navig. Satell. Syst. 2021, 25, 17. [Google Scholar] [CrossRef]
  5. Wang, D.; Dong, Y.; Li, Z.; Li, Q.; Wu, J. Constrained MEMS-Based GNSS/INS Tightly Coupled System with Robust Kalman Filter for Accurate Land Vehicular Navigation. IEEE Trans. Instrum. Meas. 2020, 69, 5138–5148. [Google Scholar] [CrossRef]
  6. Wang, S.Z.; Zhan, X.Q.; Zhai, Y.W.; Chi, C.; Liu, X.Y. Ensuring High Navigation Integrity for Urban Air Mobility Using Tightly Coupled Gnss/Ins System. J. Aeronaut. Astronaut. Aviat. 2020, 52, 429–441. [Google Scholar]
  7. Kim, Y.; An, J.; Lee, J. Robust Navigational System for a Transporter Using GPS/INS Fusion. IEEE Trans. Ind. Electron. 2018, 65, 3346–3354. [Google Scholar] [CrossRef]
  8. Basar, T. A New Approach to Linear Filtering and Prediction Problems. In Control Theory: Twenty-Five Seminal Papers (1932–1981); IEEE: Piscataway, NJ, USA, 2001; pp. 167–179. [Google Scholar]
  9. Priyanka, A. Mems-Based Integrated Navigation; Gnss Technology and Applications Series; Artech House: Boston, MA, USA, 2010. [Google Scholar]
  10. Yue, Z.; Lian, B.; Gao, Y. Robust Adaptive Filter Using Fuzzy Logic for Tightly-Coupled Visual Inertial Odometry Navigation System. IET Radar Sonar Navig. 2020, 14, 364–371. [Google Scholar] [CrossRef]
  11. Li, T.; Zhang, H.; Gao, Z.; Niu, X.; El-Sheimy, N. Tight Fusion of a Monocular Camera, MEMS-IMU, and Single-Frequency Multi-GNSS RTK for Precise Navigation in GNSS-Challenged Environments. Remote Sens. 2019, 11, 610. [Google Scholar] [CrossRef] [Green Version]
  12. Sun, R.; Yang, Y.; Chiang, K.-W.; Duong, T.-T.; Lin, K.-Y.; Tsai, G.-J. Robust IMU/GPS/VO Integration for Vehicle Navigation in GNSS Degraded Urban Areas. IEEE Sens. J. 2020, 20, 10110–10122. [Google Scholar] [CrossRef]
  13. Wang, M.; Tayebi, A. Hybrid Nonlinear Observers for Inertial Navigation Using Landmark Measurements. IEEE Trans. Autom. Control 2020, 65, 5173–5188. [Google Scholar] [CrossRef] [Green Version]
  14. Li, N.; Guan, L.; Gao, Y.; Du, S.; Wu, M.; Guang, X.; Cong, X. Indoor and Outdoor Low-Cost Seamless Integrated Navigation System Based on the Integration of INS/GNSS/LIDAR System. Remote Sens. 2020, 12, 3271. [Google Scholar] [CrossRef]
  15. Schütz, A.; Sánchez-Morales, D.E.; Pany, T. Precise Positioning through a Loosely-Coupled Sensor Fusion of Gnss-Rtk, Ins and Lidar for Autonomous Driving. In Proceedings of the 2020 IEEE/ION Position, Location and Navigation Symposium (PLANS), Portland, OR, USA, 20–23 April 2020. [Google Scholar]
  16. Batista, P.; Petit, N.; Silvestre, C.; Oliveira, P. Further Results on the Observability in Magneto-Inertial Navigation. In Proceedings of the 2013 American Control Conference, Washington, DC, USA, 17–19 June 2013; pp. 2503–2508. [Google Scholar]
  17. Zhu, K.; Guo, X.; Jiang, C.; Xue, Y.; Li, Y.; Han, L.; Chen, Y. MIMU/Odometer Fusion with State Constraints for Vehicle Positioning during BeiDou Signal Outage: Testing and Results. Sensors 2020, 20, 2302. [Google Scholar] [CrossRef] [Green Version]
  18. Ran, C.; Cheng, X. A Direct and Non-Singular Ukf Approach Using Euler Angle Kinematics for Integrated Navigation Systems. Sensors 2016, 16, 1415. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Dissanayake, G.; Sukkarieh, S.; Nebot, E.; Durrant-Whyte, H. The Aiding of a Low-Cost Strapdown Inertial Measurement Unit Using Vehicle Model Constraints for land Vehicle Applications. IEEE Trans. Robot. Autom. 2001, 17, 731–747. [Google Scholar] [CrossRef] [Green Version]
  20. Sun, W.; Yang, Y. BDS PPP/INS Tight Coupling Method Based on Non-Holonomic Constraint and Zero Velocity Update. IEEE Access 2020, 8, 128866–128876. [Google Scholar] [CrossRef]
  21. Abdolkarimi, E.S.; Mosavi, M.-R. A Low-Cost Integrated MEMS-Based INS/GPS Vehicle Navigation System with Challenging Conditions Based on an Optimized IT2FNN in Occluded Environments. GPS Solut. 2020, 24, 108. [Google Scholar] [CrossRef]
  22. Zhang, Y.; Wang, L. A Hybrid Intelligent Algorithm DGP-MLP for GNSS/INS Integration during GNSS Outages. J. Navig. 2019, 72, 375–388. [Google Scholar] [CrossRef]
  23. Doostdar, P.; Keighobadi, J.; Hamed, M.A. INS/GNSS Integration Using Recurrent Fuzzy Wavelet Neural Networks. GPS Solut. 2019, 24, 29. [Google Scholar] [CrossRef]
  24. Al Bitar, N.; Gavrilov, A. A New Method for Compensating the Errors of Integrated Navigation Systems Using Artificial Neural Networks. Measurement 2021, 168, 108391. [Google Scholar] [CrossRef]
  25. Dai, H.-F.; Bian, H.-W.; Wang, R.-Y.; Ma, H. An INS/GNSS Integrated Navigation in GNSS Denied Environment Using Recurrent Neural Network. Def. Technol. 2020, 16, 334–340. [Google Scholar] [CrossRef]
  26. Noureldin, A.; El-Shafie, A.; Bayoumi, M. GPS/INS Integration Utilizing Dynamic Neural Networks for Vehicular Navigation. Inf. Fusion 2011, 12, 48–57. [Google Scholar] [CrossRef]
  27. Lu, S.; Gong, Y.; Luo, H.; Zhao, F.; Li, Z.; Jiang, J. Heterogeneous Multi-Task Learning for Multiple Pseudo-Measurement Estimation to Bridge GPS Outages. IEEE Trans. Instrum. Meas. 2021, 70, 1–16. [Google Scholar] [CrossRef]
  28. Yao, Y.; Xu, X.; Zhu, C.; Chan, C.-Y. A Hybrid Fusion Algorithm for GPS/INS Integration during GPS Outages. Measurement 2017, 103, 42–51. [Google Scholar] [CrossRef]
  29. Fang, W.; Jiang, J.; Lu, S.; Gong, Y.; Tao, Y.; Tang, Y.; Yan, P.; Luo, H.; Liu, J. A LSTM Algorithm Estimating Pseudo Measurements for Aiding INS during GNSS Signal Outages. Remote Sens. 2020, 12, 256. [Google Scholar] [CrossRef] [Green Version]
  30. Assad, A.; Khalaf, W.; Chouaib, I. Radial Basis Function Kalman Filter for Attitude Estimation in GPS-Denied Environment. IET Radar Sonar Navig. 2020, 14, 736–746. [Google Scholar] [CrossRef]
  31. Jozefowicz, R.; Zaremba, W.; Sutskever, I. An Empirical Exploration of Recurrent Network Architectures. In Proceedings of the 32nd International Conference on International Conference on Machine Learning, Lille, France, 6–11 July 2015; Volume 37, pp. 2342–2350. [Google Scholar]
  32. Yan, G.; Weng, J. Principle of Strapdown Inertial Navigation Algorithm and Integrated Navigation; Northwestern Polytechnical University Press: Xi’an, China, 2019. [Google Scholar]
  33. Kottath, R.; Poddar, S.; Das, A.; Kumar, V. Window Based Multiple Model Adaptive Estimation for Navigational Framework. Aerosp. Sci. Technol. 2016, 50, 88–95. [Google Scholar] [CrossRef]
  34. Yang, Y.; He, H.; Xu, T. Adaptive Robust Filtering for Kinematic Gps Positioning. Acta Geod. Cartogr. Sin. 2001, 30, 293–298. [Google Scholar]
Figure 1. The structure of GRU.
Figure 1. The structure of GRU.
Remotesensing 14 00752 g001
Figure 2. Structure of GRU-aiding INS/GNSS integrated navigation. (a) Training mode; (b) predicting mode
Figure 2. Structure of GRU-aiding INS/GNSS integrated navigation. (a) Training mode; (b) predicting mode
Remotesensing 14 00752 g002
Figure 3. Data acquisition experimental platform.
Figure 3. Data acquisition experimental platform.
Remotesensing 14 00752 g003
Figure 4. Dataset collection trajectory.
Figure 4. Dataset collection trajectory.
Remotesensing 14 00752 g004
Figure 5. Position variation between GPS time 201,900 and 204,100 s. (a) North position variation; (b) East position variation.
Figure 5. Position variation between GPS time 201,900 and 204,100 s. (a) North position variation; (b) East position variation.
Remotesensing 14 00752 g005
Figure 6. Predicted and real horizontal trajectories.
Figure 6. Predicted and real horizontal trajectories.
Remotesensing 14 00752 g006
Figure 7. Velocity errors of 180 s outages estimated by different algorithms. (a) North velocity error; (b) East velocity error.
Figure 7. Velocity errors of 180 s outages estimated by different algorithms. (a) North velocity error; (b) East velocity error.
Remotesensing 14 00752 g007
Figure 8. Position errors of 180 s outages estimated by different algorithms. (a) Latitude position error; (b) longitude position error; (c) horizontal position error.
Figure 8. Position errors of 180 s outages estimated by different algorithms. (a) Latitude position error; (b) longitude position error; (c) horizontal position error.
Remotesensing 14 00752 g008
Figure 9. Position errors of 120 s outages estimated by different algorithms. (a) Latitude position error; (b) longitude position error; (c) horizontal position error.
Figure 9. Position errors of 120 s outages estimated by different algorithms. (a) Latitude position error; (b) longitude position error; (c) horizontal position error.
Remotesensing 14 00752 g009aRemotesensing 14 00752 g009b
Table 1. Hardware specifications.
Table 1. Hardware specifications.
Performance Parameter
Gyroscope bias stability8°/h
Gyroscope angle random walk0.12°/sqrt(h)
Accelerometer bias stability0.2 mg
Accelerometer velocity random walk0.09m/s/sqrt(h)
Table 2. GRU performance with different hyperparameters.
Table 2. GRU performance with different hyperparameters.
Time StepHidden
Neurons
North RMSE/°East RMSE/°Training Time/s
(per Epoch)
1324.05 × 10−86.16 × 10−80.269
2323.96 × 10−85.20 × 10−80.324
4323.99 × 10−87.18 × 10−80.486
8324.36 × 10−85.26 × 10−80.861
1644.12 × 10−84.52 × 10−80.610
2644.86 × 10−84.18 × 10−80.741
4647.83 × 10−81.06 × 10−70.991
8644.35 × 10−85.08 × 10−81.679
11283.66 × 10−83.56 × 10−81.017
21283.16 × 10−83.61 × 10−81.554
41283.05 × 10−83.04 × 10−81.900
81284.96 × 10−85.44 × 10−82.781
12563.04 × 10−83.71 × 10−82.63
22563.54 × 10−84.48 × 10−83.171
42564.41 × 10−84.75 × 10−83.903
82563.63 × 10−86.28 × 10−86.645
Table 3. The predicted position variation RMSE of the different networks.
Table 3. The predicted position variation RMSE of the different networks.
North RMSE/radEast RMSE/rad
GRU3.05 × 10−83.04 × 10−8
LSTM3.24 × 10−84.96 × 10−8
MLP6.49 × 10−76.81 × 10−7
Table 4. Maximum velocity errors after 180 s GNSS outage.
Table 4. Maximum velocity errors after 180 s GNSS outage.
North Velocity ErrorEast Velocity Error
Not-aided KF3.691 m/s2.764 m/s
MLP-aided KF3.653 m/s2.688 m/s
LSTM-aided KF3.570 m/s2.393 m/s
GRU-aided KF3.627 m/s2.352 m/s
GRU-aided AKF3.430 m/s1.750 m/s
Table 5. Position errors after the 180 s GNSS outage.
Table 5. Position errors after the 180 s GNSS outage.
Maximal
North Error
Maximal
East Error
Maximal
Horizontal Error
Horizontal Error
RMS
Not-aided KF303.866 m227.599 m379.652 m157.182 m
MLP-aided KF281.452 m116.134 m304.471 m113.758 m
LSTM-aided KF120.761 m59.109 m134.452 m40.062 m
GRU-aided KF117.314 m63.462 m133.379 m39.232 m
GRU-aided AKF78.949 m45.294 m91.019 m26.680 m
Table 6. Position errors after the 120 s GNSS outage.
Table 6. Position errors after the 120 s GNSS outage.
Maximal
North error
Maximal
East Error
Maximal
Horizontal Error
Horizontal Eerror
RMS
Not-aided KF109.639 m105.480 m152.141 m64.258 m
MLP-aided KF28.513 m30.177 m53.285 m27.342 m
LSTM-aided KF19.990 m26.512 m30.291 m17.294 m
GRU-aided KF21.227 m26.467 m31.532 m18.015 m
GRU-aided AKF21.331 m20.544 m26.450 m15.813 m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tang, Y.; Jiang, J.; Liu, J.; Yan, P.; Tao, Y.; Liu, J. A GRU and AKF-Based Hybrid Algorithm for Improving INS/GNSS Navigation Accuracy during GNSS Outage. Remote Sens. 2022, 14, 752. https://doi.org/10.3390/rs14030752

AMA Style

Tang Y, Jiang J, Liu J, Yan P, Tao Y, Liu J. A GRU and AKF-Based Hybrid Algorithm for Improving INS/GNSS Navigation Accuracy during GNSS Outage. Remote Sensing. 2022; 14(3):752. https://doi.org/10.3390/rs14030752

Chicago/Turabian Style

Tang, Yanan, Jinguang Jiang, Jianghua Liu, Peihui Yan, Yifeng Tao, and Jingnan Liu. 2022. "A GRU and AKF-Based Hybrid Algorithm for Improving INS/GNSS Navigation Accuracy during GNSS Outage" Remote Sensing 14, no. 3: 752. https://doi.org/10.3390/rs14030752

APA Style

Tang, Y., Jiang, J., Liu, J., Yan, P., Tao, Y., & Liu, J. (2022). A GRU and AKF-Based Hybrid Algorithm for Improving INS/GNSS Navigation Accuracy during GNSS Outage. Remote Sensing, 14(3), 752. https://doi.org/10.3390/rs14030752

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