Next Article in Journal
Research Progress on Power Visual Detection of Overhead Line Bolt Defects Based on UAV Images
Next Article in Special Issue
An AI-Based Deep Learning with K-Mean Approach for Enhancing Altitude Estimation Accuracy in Unmanned Aerial Vehicles
Previous Article in Journal
Advancing mmWave Altimetry for Unmanned Aerial Systems: A Signal Processing Framework for Optimized Waveform Design
Previous Article in Special Issue
Distributed Control for Multi-Robot Interactive Swarming Using Voronoi Partioning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integration of Deep Sequence Learning-Based Virtual GPS Model and EKF for AUV Navigation

1
Ocean College, Jiangsu University of Science and Technology, Zhenjiang 212003, China
2
Shenzhen Institute for Advanced Study, University of Electronic Science and Technology of China, Shenzhen 518110, China
3
Jiangsu Marine Technology Innovation Center, Nantong 226199, China
*
Author to whom correspondence should be addressed.
Drones 2024, 8(9), 441; https://doi.org/10.3390/drones8090441
Submission received: 9 August 2024 / Revised: 26 August 2024 / Accepted: 27 August 2024 / Published: 29 August 2024

Abstract

:
To address the issue of increasing navigation errors in low-cost autonomous underwater vehicles (AUVs) operating without assisted positioning underwater, this paper proposes a Virtual GPS Model (VGPSM) based on deep sequence learning. This model is integrated with an Extended Kalman Filter (EKF) to provide a high-precision navigation solution for AUVs. The VGPSM leverages the time-series characteristics of data from sensors such as the Attitude and Heading Reference System (AHRS) and the Doppler Velocity Log (DVL) while the AUV is on the surface. It learns the relationship between these sensor data and GPS data by utilizing a hybrid model of Long Short-Term Memory (LSTM) and Bidirectional Long Short-Term Memory (Bi-LSTM), which are well-suited for processing and predicting time-series data. This approach constructs a virtual GPS model that generates virtual GPS displacements updated at the same frequency as the real GPS data. When the AUV navigates underwater, the virtual GPS displacements generated using the VGPSM in real-time are used as measurements to assist the EKF in state estimation, thereby enhancing the accuracy and robustness of underwater navigation. The effectiveness of the proposed method is validated through a series of experiments under various conditions. The experimental results demonstrate that the proposed method significantly reduces cumulative errors, with navigation accuracy improvements ranging from 29.2% to 69.56% compared to the standard EKF, indicating strong adaptability and robustness.

1. Introduction

With the rapid development of marine science and technology, the demand for autonomous underwater vehicles (AUVs) in civil, military, and scientific research fields is increasing [1]. As the core component of AUV autonomy and intelligence, the positioning accuracy of the navigation system directly affects the success and efficiency of AUV underwater missions [2,3]. However, the special characteristics of the underwater environment bring great challenges to the navigation and positioning of AUVs [4,5]. Therefore, it is of great research significance to improve the navigation accuracy of underwater vehicles.
Currently, AUV underwater navigation sensors mainly include the inertial navigation system (INS) [6], Doppler Velocimetry (DVL) [7], a Depth System (PS) [8], and a Global Positioning System (GPS) [9]. In general, INS can provide AUVs with acceleration, angular velocity, angle, velocity, and position information. However, since the calculation of velocity and position relies on the integral accumulation of acceleration and angular velocity, it requires the high accuracy of gyroscopes and accelerometers. Conventional AUVs generally use an INS that cannot provide velocity and position information, such as inertial measurement units (IMUs) and attitude and heading reference systems (AHRSs). DVL provides velocity relative to the bottom through acoustic Doppler frequency shift, which assists inertial navigation in generating position information. PS provides the AUV with depth information between the current position and the water surface. GPS can provide accurate position information on the water surface, but the signal cannot be transmitted in underwater environments [10]. An ultra-short baseline (USBL) [11] can provide underwater positioning; however, it is expensive and relies on external equipment. Therefore, the integrated navigation technology of INS/DVL/PS/GPS has become a research hotspot for low-cost AUV underwater navigation [2,3], which achieves a higher accuracy of navigation through fusing data from multiple sensors. However, it is still a challenge to optimize the sensor data fusion algorithm in the underwater complex environment.
In practical applications, commonly used navigation algorithms include dead reckoning (DR) [12] and the extended Kalman filter (EKF) [13]. DR estimates the position of the AUV by continuously calculating its motion state, but is easily affected by cumulative errors. EKF effectively reduces the cumulative error and provides more accurate navigation information through multiple filtering and the fusion of sensor data. However, the performance of EKF depends on accurate sensor models and environmental models, and the accuracy is still limited without the assistance of external position information [14]. In recent years, machine learning [15,16], deep learning [17,18], and other neural network-based methods in AUV navigation [19] have gained widespread attention.
More and more scholars are exploring the impact of machine learning and deep learning methods on improving the accuracy of AUV navigation systems [20,21,22,23]. He et al. summarized how deep neural networks (DNN) can suppress the drift error of navigation systems by calibrating inertial sensor noise and fusing navigation sensor information [24]. Xie et al. proposed a neural network-based approach to explore the time-varying relationship between the acceleration and attitude angle, using acoustic localization and a neural network to estimate the pitch angle and reduce the impact of waves on inertial sensors [25]. Huang et al. proposed a multi-error model-assisted extended Kalman filter method for attitude estimation, which improves the performance of MEMS sensors and enhances the accuracy of underwater glider navigation systems [26]. Li et al. constructed a nonlinear regressive neural network model for the supervised learning of the velocity increments of DVL and INS, assisting the navigation when the DVL fails [27]. Sabet et al. proposed a dynamic surge identification model based on low-cost MEMS to assist the dead reckoning, which improves the navigation accuracy by reducing the effects of external acceleration and electromagnetic interference on attitude [28]. Topini et al. used four deep learning models, including MLP [29], LSTM [30], CNN [31], and CNN-LSTM [32], to estimate DVL velocity and combined dead reckoning for navigation [33]. The above methods are mainly based on supervised learning models outputting information such as attitude and velocity to assist in improving navigation and positioning accuracy.
Since GPS can provide the navigation system with position data without cumulative errors when the AUV is sailing on the surface, the divergence of the navigation system error is suppressed. However, the GPS information is not available due to the attenuation of radio signals underwater, and the USBL data update rate is low, which makes it impossible to control the error dispersion of the navigation system. Scholars have begun to study the direct relationship between the inertial navigation system, the Doppler velocity log, the depth sensor, and GPS positioning data. The most common method is to achieve end-to-end position output through a neural network [34]. Lv et al. proposed an adaptive GPS filtering method to address the GPS drift caused by the GPS antenna being easily covered when the AUV is sailing near the surface. By comparing the length of time without GPS data and the distance difference between GPS position and navigation position with the threshold value, the GPS drift data are filtered out and the confidence of GPS position is improved [35]. Then, an AUV navigation dataset with GPS correction is constructed, a hybrid recurrent neural network is selected to learn the AUV navigation data, and a position correction model is constructed for the position output of the navigation system [36]. Guo et al. estimated the AUV position based on an online process model by reconstructing state information between valid data from an ultrashort baseline positioning system [37]. Mu et al. trained GPS and other navigation data directly through LSTM and bidirectional long- and short-term memory, and obtained the AUV position data [38] through an end-to-end approach [39].
Although the above navigation technologies have suppressed the divergence rate of navigation system errors to a certain extent, they still have certain limitations [40]. For example, deep learning methods rely on a large amount of high-quality data and computing resources. To overcome the above problems, this paper proposes a navigation method that combines the traditional EKF and the virtual GPS model (VGPSM) based on deep sequence learning to provide high-precision and high-update rate navigation information for an autonomous underwater vehicle. The main contributions of this paper are summarized as follows:
  • According to the high-quality navigation sensor data, when GPS data are available on the surface, based on the time-series characteristics of the navigation system data, a deep sequence learning model integrating LSTM and Bi-LSTM is constructed to learn the relationship between the GPS data and the DVL and AHRS data in the surface environment, and forms a virtual GPS model, which solves the problem of no GPS signal underwater.
  • Observation vectors are constructed for two cases, depending on whether the AUVs have positional aids when performing the mission. The virtual GPS model is updated at the same frequency as the real GPS data. When there is an update in the surface GPS or the underwater VGPSM, the corresponding position is added to the observation vector. Otherwise, the observation vector does not contain position information, which effectively solves the limitations of traditional navigation methods when there is no external position assistance.
  • Through experiments on real navigation data, the performance and practical application effects of the proposed method in different situations are verified, which provides new ideas and references for the development of underwater vehicle navigation technology.
The structure of this paper is organized as follows: Section 2 introduces the basic knowledge of navigation system. Section 3 describes the specific implementation process of a virtual GPS model based on deep sequence learning and the proposed methodology. Section 4 conducts experimental design and results, followed by Section 5, which offers an in-depth analysis and discussion of the experimental findings. Finally, Section 6 summarizes the key conclusions and contributions of this study. Through these contents, this study hopes to provide new theoretical and technical support for the development of AUV navigation systems.

2. Navigation System

The AUV navigation system contains functions such as sensor data acquisition and navigation positioning. The data collected using each sensor will be sent to the data center (MOOSDB). The navigation algorithm module, after acquiring the sensor data at the current moment from the MOOSDB, calculates a stable and reliable position, velocity, attitude, and other information in real time, and then sends it to the MOOSDB for the other systems of the AUV to acquire in order to perform the task more accurately. This section focuses on the two coordinate systems and the main onboard sensors of the AUV navigation system in this paper.

2.1. Coordinate System for AUV Navigation

The navigation system in this paper uses both the navigation coordinate system and the body coordinate system, as shown in Figure 1. In general, the position, velocity, attitude, acceleration, and angular velocity of the AUV come from the navigation coordinate system by default, but the velocity and acceleration of the AUV during the control system’s manipulation of the AUV’s movement come from the AUV coordinate system.
The letter n is generally used to denote the navigation coordinate system, whose origin O n is a point fixed to the Earth, generally starting from a stable position obtained after the navigation system is switched on. The O n X n axis points geographical north along the tangent direction of the reference ellipsoid meridian, the O n Y n axis points geographical east along the normal direction of the reference ellipsoid, and the O n Z n axis forms a right-handed coordinate system with the O n X n axis and the O n Y n axis in the local horizontal plane, which is also known as the North-East-Down (NED) coordinate system due to its orientation.
The letter b is generally used to denote the carrier coordinate system, whose origin O b is fixed at the center of mass and moves with the movement of the AUV. The O b X b axis is in the same direction as the roll axis of the AUV’s angular motion and points to the front of the AUV. The O b Y b axis is in the same direction as the pitch axis of the AUV’s angular motion and points to the right of the AUV. The O b Z b axis is in the same direction as the heading axis of the AUV’s angular motion, and forms a right-handed coordinate system with the O b X b axis and the O b Y b axis, also known as the Forward-Right-Down (FRD) coordinate system.

2.2. Onboard Sensors for AUV Navigation

This section focuses on the onboard sensors of the navigation system as shown in Figure 1, as well as their performance parameters and output data types, which are shown in Table 1.
The attitude and heading reference system (AHRS), model Ellipse-A AHRS, provides heading ( ξ z b ), pitch ( ξ y b ), roll ( ξ x b ), three-axis acceleration ( a x b , a y b , a z b ,), and three-axis angular velocity ( w x b , w y b , w z b ) for the navigation system. The data of AHRS are mainly obtained based on three-axis gyroscopes, three-axis accelerometers, and three-axis magnetometers. Due to the use of the magnetometer, the AHRS needs to be calibrated for the local magnetic field every time the AUV is tested in a new environment.
The Doppler velocity log (DVL) sensor, model NavQuest 600 micro DVL, is used to provide three-dimensional velocity ( v x b , v y b , v z b ) and bottom height ( a l z b ) for the navigation system. The transducer array of DVL adopts a four-beam JANUS configuration, which is able to reduce the adverse effects of AUV motion on DVL measurement, and has a high accuracy in the velocity measurement. The bottom height can assist in determining the validity of DVL data.
The global positioning system (GPS), model LEA-M8T GPS, provides latitude and longitude with no accumulated errors for the navigation system. This sensor can provide high-precision data near the water surface where its antenna is not covered, but it is unable to provide a position underwater due to the rapid decay of radio signals, and this paper focuses on improving this.
The intelligent pressure system (IPS), model miniIPS, provides depth ( p z n ) for navigation system. The sensor obtains the current depth with high accuracy based on the effect of water pressure on the resistance value, and can provide it to the AUV continuously.

3. Proposed Method for AUV Navigation

As a commonly used method for AUV navigation, EKF can efficiently handle nonlinear systems with good real-time performance and robustness, but it still inevitably suffers from accumulated errors. Sensor noise and model uncertainty lead to gradual deviations in state estimation over time, resulting in cumulative errors. Therefore, these errors need to be corrected by periodic correction to ensure the long-term accuracy and reliability of the navigation system. The navigation system can use GPS position data for error correction while sailing on the surface, but it cannot be performed underwater. We propose a virtual GPS model based on deep sequence learning to address this issue, which outputs the virtual GPS position during AUV underwater navigation for error correction. This section will describe the proposed method in detail, focusing on the proposed virtual/real GPS-assisted EKF and the virtual GPS model.

3.1. Proposed Navigation Method Based on Virtual/Real GPS

The flowchart of the proposed navigation method based on virtual/real GPS is shown in Figure 2. First, a virtual GPS model is constructed and trained based on sensor data, which is consistent with the update frequency of the GPS data. Then, when the AUV navigates near the water surface, its position is assigned to the measurement vector Z 2 k R 9 if there is a GPS data update, corresponding to “Case2” in Figure 2; otherwise, the measurement vector Z 1 k R 7 has no two-dimensional (2D) position, corresponding to “Case1”. When the AUV has no GPS data underwater, the VGPSM position is assigned to the measurement vector Z 2 k R 9 if there is an update of the virtual GPS model; otherwise, there is no 2D position in the measurement vector Z 1 k R 7 . Finally, the AUV performs navigation positioning based on the prediction and update steps of the EKF.
The motion model and measurement model of AUV are shown in Equation (1) and Equation (2), respectively.
X k = f ( X k 1 , U k , W k )
Z k = h ( X k , V k )
W k N ( 0 , Q k )
V k N ( 0 , R k )
where f ( · ) and h ( · ) represent the motion function and measurement function, respectively, X k R 9 is the state vector at time k as shown in Equation (5), W k R 9 is Gaussian noise at time k with a mean zero and covariance matrix Q k as shown in Equation (3), U k R 6 is the control vector at time k as shown in Equation (6), V k is Gaussian noise at time k with a mean zero and covariance matrix R k as shown in Equation (4), and Z k is the measurement vector at time k. When there is no update from GPS or VGPSM, Z k is expressed as Z 1 k as shown in Equation (7); otherwise, Z k is expressed as Z 2 k , as shown in Equation (8).
X k = p x , k n , p y , k n , p z , k n , v x , k b , v y , k b , v z , k b , ξ x , k , ξ y , k , ξ z , k T
U k = a x , k b , a y , k b , a z , k b , w x , k , w y , k , w z , k T
Z 1 k = p z , k I P S , v x , k D V L , v y , k D V L , v z , k D V L , ξ x , k A H R S , ξ y , k A H R S , ξ z , k A H R S T
Z 2 k = p x , k G P S / V G P S M , p y , k G P S / V G P S M , p z , k I P S , v x , k D V L , v y , k D V L , v z , k D V L , ξ x , k A H R S , ξ y , k A H R S , ξ z , k A H R S T
where
  • p x , k n , p y , k n , p z , k n : the position of the x, y, and z axes in the navigation coordinate system.
  • v x , k b , v y , k b , v z , k b : the velocity of the x, y, and z axes in the body coordinate system.
  • ξ x , k , ξ y , k , ξ z , k : the three-dimensional (3D) attitude of the AUV for roll, pitch, and heading.
  • a x , k b , a y , k b , a z , k b : the 3D acceleration of the AUV in the body coordinate system.
  • w x , k , w y , k , w z , k : the 3D angular velocities of the AUV.
  • p x , k G P S / V G P S M , p y , k G P S / V G P S M : the position of the x and y axes in the navigation coordinate system from GPS or VGPSM.
  • p z , k I P S : the position of the z axes in the navigation coordinate system from IPS.
  • v x , k D V L , v y , k D V L , v z , k D V L : the 3D velocity in the body coordinate system from DVL.
  • ξ x , k A H R S , y , k A H R S , ξ z , k A H R S : the 3D attitude of the AUV for roll, pitch, and heading from AHRS.
The pseudo-code of the proposed method is shown in Algorithm 1. When the navigation system is initialized, the state vector X 0 and covariance matrix P 0 are set up at the initial moment. At this stage, the current GPS position ( l a t , l o n ) is recorded and converted into Universal Transverse Mercator (UTM) coordinates ( p x , 0 U T M , p y , 0 U T M ) . During the AUV navigation process, the system primarily relies on two stages: the time update (prediction) and the measurement update (correction) stages. The correction stage is divided into two scenarios: one where the measurement vector Z 1 k R 7 does not include the two-dimensional position, corresponding to Case 1, and another where the measurement vector Z 2 k R 9 includes the two-dimensional position, corresponding to Case 2. The detailed procedures for these two scenarios are described in Section 3.1.1 and Section 3.1.2, respectively.
Algorithm 1: The pseudo-code of proposed method
Drones 08 00441 i001

3.1.1. The Prediction Stage of Proposed Navigation Method

The prediction and update of the state vector X k | k 1 and the covariance matrix P k | k 1 are performed at each time step k according to the update frequency of the navigation system. The prediction steps are consistent with Figure 2, as shown in Equations (9) and (10).
X k | k 1 = f ( X k 1 , U k , 0 ) = p x , k 1 n + C 11 δ p x , k b + C 12 δ p y , k b + C 13 δ p z , k b p y , k 1 n + C 21 δ p x , k b + C 22 δ p y , k b + C 23 δ p z , k b p z , k 1 n + C 31 δ p x , k b + C 32 δ p y , k b + C 33 δ p z , k b v x , k 1 b + a x , k b δ t v y , k 1 b + a y , k b δ t v z , k 1 b + a z , k b δ t ξ x , k 1 + w x , k δ t ξ y , k 1 + w y , k δ t ξ z , k 1 + w z , k δ t
P k | k 1 = F k 1 P k 1 F k 1 T + Q k
where δ p x , k b , δ p y , k b , δ p z , k b denote the displacement of the x, y, and z axes in the body coordinate system between time steps k 1 and k as shown in Equation (11), δ t denotes the time between time step k 1 and k, F k 1 is the Jacobi matrix of the motion function f ( · ) on the state vector X k 1 as shown in Equation (12), and C i j ( 1 i 3 , 1 j 3 ) denotes the data in the i-th row and j-th column of the rotation matrix C b n that transform from the body coordinate system to the navigation coordinate system, as shown in Equation (13).
Δ P k b = δ p x , k b δ p y , k b δ p z , k b = v x , k 1 b + 1 2 a x , k b δ t 2 v y , k 1 b + 1 2 a y , k b δ t 2 v z , k 1 b + 1 2 a z , k b δ t 2
F k 1 = f ( X k 1 , U k , 0 ) X k 1
C b n = cos ξ y cos ξ z cos ξ y sin ξ z sin ξ y cos ξ x sin ξ y cos ξ z + cos ξ x sin ξ z sin ξ x sin ξ y sin ξ z + cos ξ x cos ξ z sin ξ y cos ξ y cos ξ x sin ξ y cos ξ z + sin ξ x sin ξ z cos ξ x sin ξ y sin ξ z + sin ξ x cos ξ z cos ξ x cos ξ y

3.1.2. The Correction Stage of Proposed Navigation Method

After the prediction stage, the measurement update stage follows. When the AUV is sailing on the surface, if the GPS data are updated at this time, the current GPS position ( l a t k G P S , l o n k G P S ) is converted to the UTM coordinate position ( p x , k U T M , p y , k U T M ) ; then, the measured position ( p x , k G P S , p y , k G P S ) , as shown in Equation (14), is obtained and assigned to Z 2 k , and the measured values of IPS, DVL, and AHRS are also assigned to Z 2 k , as shown in Equation (15).
p x , k G P S p y , k G P S = p x , k U T M p x , 0 U T M p y , k U T M p y , 0 U T M
Z 2 k = p x , k G P S , p y , k G P S , p z , k I P S , v x , k D V L , v y , k D V L , v z , k D V L , ξ x , k A H R S , ξ y , k A H R S , ξ z , k A H R S T
The measurement value Z 2 k | k 1 , as shown in Equation (16), is predicted based on the measurement function h 2 ( · ) , where H 2 k is a nine dimensional identity matrix as shown in Equation (17).
Z 2 k | k 1 = h 2 ( X k | k 1 , 0 ) = H 2 k X k | k 1
H 2 k = h 2 ( X k | k 1 , 0 ) X k | k 1
The Kalman gain K k at time step k can be calculated using Equation (18). Then the state vector X k and its covariance Matrix P k are corrected as shown in Equations (19) and (20).
K k = P k | k 1 H 2 k T [ H 2 k P k | k 1 H 2 k T + R k ] 1
X k = X k | k 1 + K k [ Z 2 k Z 2 k | k 1 ]
P k = [ I K k H 2 k ] P k | k 1
If the virtual GPS model is updated when the AUV is sailing underwater, the measurement position ( p x , k V G P S M , p y , k V G P S M ) , as shown in Equation (21), is obtained based on the displacement output by VGPSM ( δ p x , k V G P S M , δ p y , k V G P S M ) and the AUV position ( p x , k m n , p y , k m n ) at time step k m , which is assigned to Z 2 k together with the measurement data of IPS, DVL, and AHRS, as shown in Equation (22).
p x , k V G P S M p y , k V G P S M = p x , k m n + δ p x , k V G P S M p y , k m n + δ p y , k V G P S M
Z 2 k = p x , k V G P S M , p y , k V G P S M , p z , k I P S , v x , k D V L , v y , k D V L , v z , k D V L , ξ x , k A H R S , ξ y , k A H R S , ξ z , k A H R S T
where m is the number of time steps within one GPS update cycle. Then the updated state vector X k and its covariance matrix P k can be obtained based on Equations (18)–(20).
The measurement vector Z 1 k , as shown in Equation (7), does not contain a 2D position if there is no GPS data on the surface and no VGPSM data underwater, which is consistent with the situation of “Case 1” in Figure 2. The measurement value Z 1 k | k 1 , as shown in Equation (23), is predicted based on the measurement function h 1 ( · ) , where the 7 × 9 matrix H 1 k , as shown in Equation (24), is the Jacobian matrix of the measurement function h 1 ( · ) with respect to the predicted state vector X k | k 1 .
Z 1 k | k 1 = h 1 ( X k | k 1 , 0 ) = H 1 k X k | k 1
H 1 k = h 1 ( X k | k 1 , 0 ) X k | k 1
The Kalman gain K k at time step k can be calculated using Equation (25). Then the state vector X k and its covariance Matrix P k are corrected, as shown in Equations (26) and (27).
K k = P k | k 1 H 1 k T [ H 2 k P k | k 1 H 1 k T + R k ] 1
X k = X k | k 1 + K k [ Z 1 k Z 1 k | k 1 ]
P k = [ I K k H 1 k ] P k | k 1

3.2. Virtual GPS Model

Aiming at the problem of an increasing cumulative error in the navigation system without underwater positioning assistance, we propose to use a virtual GPS model, as shown in Figure 3, for the underwater position assistance of the AUV navigation system, in order to slow down the increase in cumulative error.

3.2.1. The Construction of Virtual GPS Model

VGPSM is trained and generated based on deep sequence learning, which comprehensively utilizes long short-term memory and bidirectional long short-term memory to effectively capture long-term dependencies in time-series navigation data, and, when combined with a two-layer neural network to learn the nonlinear relationship between navigation sensor data and GPS data, generates virtual GPS displacement.
The construction process of the virtual GPS model is as follows:
  • First, the construction of high-quality datasets. Deep sequence learning models can efficiently handle the effects of nonlinearity and noise, but their training process relies on a large amount of high-quality data. Therefore, based on the data from AHRS, DVL, and GPS, we combine improved incremental smoothing and mapping (iSAM) and an adaptive GPS filter to generate position data that are consistent with the update frequency of the navigation system, and then construct them into a high-quality navigation dataset.
  • Then, the structural design of a deep sequence model. Based on the time-series characteristics of navigation data, LSTM, which can capture long-term dependencies, and Bi-LSTM, which considers both forward and backward sequence data, are selected to process navigation sequence data. After that, data integration prediction is carried out by combining a 2-layer fully-connected layer, a Dropout layer, which prevents overfitting, and a Regression layer. Finally, the predicted virtual GPS position can be obtained using the output layer.
  • Finally, the training and prediction process of the virtual GPS model. In total, 80% of the navigation dataset constructed above is used for training and 20% is used for testing. The input and output of the dataset are normalized before the training begins. The training process uses the mean square error (MSE) as the loss function, the optimizer is selected as Adam, and the initial learning rate is set to 0.001. After the DVL and AHRS data are processed with the same normalization parameters, they are input into the trained virtual GPSM model to generate the virtual GPS displacement.
Before constructing the navigation dataset, it is necessary to collect and save the data of the same AUV under various trajectories when there is an update of GPS on the surface, including the 3D attitude of DVL, the 2D position of GPS, the 3D acceleration, the 3D angular velocity, and the 3D attitude of AHRS. The update frequency of the navigation system is 10 Hz, and the update frequency of the GPS is 1 Hz. Since the GPS antenna may be covered by waves on the water surface causing drift errors, the above data are put through improved iSAM and the adaptive GPS filter for the navigation computation, and the navigation data with the update frequency of 10 Hz are generated uniformly.
Based on the described procedures, we can construct a dataset ‘NavD’ specifically for AUV navigation. Considering the navigation update frequency, the dataset is sampled at 0.1 s intervals. Taking the j-th row as an example, the input vector I j n shown in Equation (28) consists of sensor data including DVL, AHRS, IPS, and GPS with a feature set of 15 in every time step j.
I j n = v x , j D V L , v y , j D V L , v z , j D V L , ξ x , j A H R S , ξ y , j A H R S , ξ z , j A H R S , a x , j A H R S , a y , j A H R S , a z , j A H R S , w x , j A H R S , w y , j A H R S , w z , j A H R S , p z , j I P S , l a t j G P S , l o n j G P S
The output vector O j n of ‘NavD’ shown in Equation (29) includes the northward and eastward positions ( p x , j n , p y , j n ) , as well as the northward and eastward displacements ( p x , j n p x , j 1 n , p y , j n p y , j 1 n ) at time step j.
O j n = p x , j n , p y , j n , p x , j n p x , j 1 n , p y , j n p y , j 1 n
Based on the existing dataset ‘NavD’, the relevant portions are selected to construct a specialized dataset “VGD” for the training and testing of the virtual GPS model. As shown in Figure 3, the input sequence x ( i ) at time step i is shown in Equation (30).
x ( i ) = I m × ( i 1 ) + 1 n ( 1 : 12 ) , I m × ( i 1 ) + 2 n ( 1 : 12 ) , I m × ( i 1 ) + 3 n ( 1 : 12 ) , , I m × i n ( 1 : 12 )
where I i × m n ( 1 : 12 ) is the DVL and AHRS data at time step i × m , and m is the number of time steps within one GPS update cycle.
The output vector y ( i ) of ‘VGD’ shown in Equation (31) is the northward and eastward displacements of m time steps in the navigation coordinate system, and its output frequency is consistent with GPS.
y ( i ) = O m × i n ( 1 : 2 ) O m × ( i 1 ) n ( 1 : 2 )
The training of the virtual GPS model is based on the input x ( i ) and output y ( i ) from the dataset ‘VGD’. The model starts with an input layer, where each sequence contains 10 time steps, with 12 features per time step. The hidden layer is a hybrid gated recurrent neural network (HGRNN) composed of an LSTM and a bidirectional LSTM (Bi-LSTM). The Bi-LSTM layer has 128 units, followed by a unidirectional LSTM layer with 64 units. The output is then passed through a fully connected layer (FC1) with 64 units, followed by a Dropout layer with a dropout rate of 0.5 to prevent overfitting. Afterward, the data go through another fully connected layer (FC2) with 32 units, and finally, the regression layer outputs two features. The model is trained with a batch size of 64, a learning rate of 0.001, and for 100 epochs. The loss function used is Mean Squared Error (MSE), optimized using the Adam optimizer.
The virtual GPS model is imported into the AUV navigation system for the virtual GPS displacement output y o ( i ) , as shown in Equation (32), during navigation.
y o ( i ) = δ p x , i V G P S M , δ p y , i V G P S M
where δ p x , i V G P S M and δ p y , i V G P S M denote the northward and eastward displacement over m time steps.

3.2.2. Bi-LSTM and LSTM

A hybrid recurrent gated neural network consisting of Bi-LSTM and LSTM is used in the virtual GPS model. The Bi-LSTM, as shown in Figure 4, consisting of forward and backward LSTM, can learn the bidirectional long-term dependencies of sequence navigation data. The structure of the LSTM cell can integrate long-term memory and short-term memory based on gating with different functions.
LSTM is a recurrent neural network (RNN) that can effectively solve the gradient vanishing problem in traditional RNN, and is suitable for time-series-based data prediction. As shown in Figure 4, the LSTM cell consists of an input gate, a forget gate, and an output gate.
As described in the yellow dashed box in Figure 4, the forget gate processes the output s ( t 1 ) of the previous moment and the input x ( t ) of the current moment through a sigmoid layer to obtain the activation vector f ( t ) , shown in Equation (33).
f ( t ) = σ ( W f s ( t 1 ) + U f x ( t ) + b f )
As shown in the blue dashed box, the input gate combines s ( t 1 ) and x ( t ) through the sigmoid layer to obtain the input gate’s activation vector i ( t ) , as shown in Equation (34), and through the tanh layer to obtain the cell input activation vector i t , as shown in Equation (35).
i ( t ) = σ ( W i s ( t 1 ) + U i x ( t ) + b i )
c 1 ( t ) = t a n h ( W c s ( t 1 ) + U c x ( t ) + b c )
The light green line represents long-term memory. The forget gate determines the degree to which the long-term memory c ( t 1 ) at t 1 moment is forgotten, and the input gate selects the degree to which the short-term memory is increased, which together determine the long-term memory c ( t ) of the current moment, as shown in Equation (36).
c ( t ) = g i c 1 ( t ) + f ( t ) c ( t 1 )
As indicated by the pink dashed box, the combination of s ( t 1 ) and x ( t ) is passed through the sigmoid layer to obtain the output gate’s activation vector o ( t ) , as shown in Equation (37). Combined with the tanh layer’s processing of c ( t ) , the output gate’s activation vector s ( t ) , as shown in Equation (38), is obtained.
o ( t ) = σ ( W o s ( t 1 ) + U o x ( t ) + b o )
s ( t ) = o ( t ) t a n h ( c ( t ) )
where W and U denote the weight matrices, and b denotes the bias vector parameter.

4. Experimental Results

This section will mainly present the performance of the virtual GPS model and its correction of the navigation system in combination with the real GPS position.
To validate the proposed algorithm, we conducted a series of carefully controlled experiments. The experimental data were collected from the Sailfish AUV, shown in Figure 2, under optimal sea conditions. The datasets used for this validation are ‘NavD’ and ‘VGD,’ as described in Section 3.2.1. The ‘NavD’ dataset, with a time interval of 0.1 s, serves as the foundational dataset from which the ‘VGD’ dataset was derived. ‘VGD’ was constructed by subsampling ‘NavD’ at 1-second intervals to match the GPS update frequency.
The sensor data incorporated in these experiments include 3D attitude, 3D acceleration, and 3D angular velocity from AHRS, 3D velocity from DVL, latitude and longitude from GPS, and depth measurements from IPS. The sensors used in these experiments are detailed in Table 1, which lists their models and key specifications. Furthermore, the primary parameters employed in the experiments are summarized in Table 2.

4.1. The Results of Virtual GPS Model

The virtual GPS model provides position information for the AUV underwater without GPS assistance to improve the navigation system accuracy. Figure 5 shows the result of the virtual GPS model. The northward and eastward displacement of GPS is taken as the true value, and that of the virtual GPS model is compared in terms of displacement waveform, root mean square error (RMSE), and single-step error.
In Figure 5, the red and blue lines represent the true values of northward and eastward displacement, respectively, and the cyan and magenta lines represent the northward and eastward displacement of the virtual GPS model, respectively. (a) and (b) represent the northward and eastward displacement comparison between the virtual GPS model and the GPS truth. (c) and (d) represent the northward and eastward root mean square error between the virtual GPS model and the ground truth, respectively. (e) and (f) compare the error between the virtual GPS model and the truth in two ways, where (e) shows the error at each step, and (f) shows the statistics of all errors.
In order to analyze the performance of the virtual GPS model specifically, we formed a table, as shown in Table 3, which mainly analyzes the performance of the virtual GPS model in terms of the mean, standard deviation, and median of the errors, as well as the root mean square error.
As shown in Figure 5e, since there are positive and negative errors, we conducted two kinds of analyses on the mean and standard deviation of the errors. One is to take the absolute value of all errors and then solve for the mean (|Error-M|), as shown in Equation (39), and the standard deviation (|Error|-S), as shown in Equation (40). The other is to directly solve for the mean of the errors (Error-M), as shown in Equation (41), and the standard deviation (Error-S), as shown in Equation (42).
| E r r o r | M = 1 T m = 1 T | x m x ^ m |
| E r r o r | S = 1 T m = 1 T ( | x m x ^ m | | E r r o r | M ) 2
E r r o r M = 1 T m = 1 T ( x m x ^ m )
E r r o r S = 1 T m = 1 T ( x m x ^ m E r r o r M ) 2
The calculation of RMSE is shown in Equation (43).
R M S E = 1 T m = 1 T ( x m x ^ m ) 2
where x m denotes the true value of northward or eastward displacement at time step m, x ^ m denotes the predicted value of northward or eastward displacement output using the virtual GPS model at time step m, and T denotes the total time steps.

4.2. The Results of Proposed Algorithm

To evaluate the effectiveness of the proposed algorithm, we have meticulously designed three experimental scenarios.
Scenario 1, as depicted in Figure 6, provides continuous GPS assistance throughout the entire operation. In this scenario, both EKF-GPS and the proposed algorithm receive position updates directly from the GPS, while EKF-VGPSM utilizes a virtual velocity model for position correction.
Scenario 2, illustrated in Figure 7, involves GPS assistance only during the first third of the operation, followed by a period without GPS support during the middle third. Specifically, in this scenario, EKF-GPS and the proposed algorithm are aided by GPS during the initial phase. During the subsequent two-thirds of the operation, the proposed algorithm relies on a virtual GPS model for position updates, while EKF-VGPSM depends entirely on the virtual GPS model for the entire duration.
Scenario 3, shown in Figure 8, operates without any GPS assistance for the entire duration. In this case, EKF-GPS behaves identically to a standard EKF without GPS, while the proposed algorithm functions equivalently to EKF-VGPSM, both relying solely on the virtual GPS model for positional updates.
These three scenarios are designed to comprehensively assess the performance of the proposed algorithm under varying conditions, particularly focusing on its robustness and accuracy in the absence of real-time GPS data, which is critical for underwater navigation systems.
In Figure 6, Figure 7 and Figure 8, the black line represents the ground truth, the red line denotes the trajectory estimated using the standard extended Kalman filter (EKF), the blue line illustrates the trajectory estimated using the EKF augmented with the virtual GPS model (VGPSM), the green line depicts the trajectory estimated using the EKF with GPS positional assistance (EKF-GPS), and the magenta line represents the trajectory generated using the proposed algorithm. Subfigure (a) shows a trajectory comparison among the aforementioned algorithms against the ground truth. Subfigure (b) compares the root mean square error (RMSE) of the different methods. Subfigure (c) contrasts the single-step error across the algorithms, while Subfigure (d) breaks down the single-step error into the northward and eastward components.
To enhance the clarity of the results presented in Figure 6, Figure 7 and Figure 8, we summarize the key findings in Table 4. This table primarily includes the trajectory length, the root mean square error (RMSE), the endpoint error, the northward error (Error-N), and the eastward error (Error-E) at the endpoint, and the overall navigation accuracy for each scenario. The RMSE represents the root mean square of the displacement errors, calculated as shown in Equation (44). The displacement error at time step i is determined according to Equation (45), and the navigation accuracy is computed using Equation (46). Additionally, the optimal values achieved by each algorithm across the three scenarios are highlighted in bold to facilitate a clearer comparison of the strengths and weaknesses of each method.
R M S E = 1 T i = 1 T E r r o r i 2
E r r o r i = ( x i x ^ i ) 2 + ( y i y ^ i ) 2
A c c u r a c y = R M S E D i s t a n c e
where x i and y i denote the ground truth of the displacement at time step i in the north and east directions, respectively, x ^ i and y ^ i denote the outputs of the virtual GPS model at time step i in the north and east directions, respectively, T denotes the total number of time steps, and D i s t a n c e denotes the total distance of this set of data.

5. Experimental Analysis

5.1. The Performance Analysis of Virtual GPS Model

In order to verify the effectiveness of the virtual GPS model, we first carried out the comparison between the virtual GPS model and the true GPS data in the northward and eastward displacements, the results of which are shown in Figure 5 and Table 3.
From the displacement comparison in (a) and (b), it can be seen that the northward and eastward displacements output by the virtual GPS model are very close to the true displacement values of GPS. From the comparison of the root mean square error in the subfigures (c) and (d), it can be seen that the RMSE of northward and eastward displacements is basically stable at around 0.05 m. Combined with Table 3, the root mean square errors of the virtual GPS model in the north and east directions throughout the entire process are 0.0538 m and 0.0502 m, respectively. From the error comparisons in (e) and (f), it can be seen that the single-step error fluctuates around 0.1 m to 0.1 m, and the median error is close to 0 m. Combined with Table 3, it can be seen that the mean and covariance of the errors between the northward displacement and the true value of the virtual GPS model are 0.0019 m and 0.0538 m, respectively, and the mean and covariance of the errors between the virtual GPS model and the true value in northward displacement are 0.0019 m and 0.0538 m, respectively, and those in eastward displacement are −0.0059 m and 0.0498 m, respectively. To show the distribution of the data, the mean of the absolute value of the error is higher than that of the error, and the mean squared error of the absolute value of the error is lower than that of the error.

5.2. The Performance Analysis of Proposed Method

To evaluate the effectiveness of the proposed algorithm, we conducted three experimental scenarios under different conditions. In the first scenario, depicted in Figure 6, both the EKF-GPS and the proposed algorithm operated with continuous GPS assistance throughout the entire process. In this setup, EKF-VGPSM relied on the position data provided by the output of the virtual GPS model. The results clearly show that the navigation accuracy of EKF-VGPSM was significantly better than that of the standard EKF. Moreover, the trajectories produced by both EKF-GPS and the proposed algorithm coincided and were noticeably more accurate than those of EKF-VGPS alone. As summarized in Table 4, the proposed algorithm and EKF-GPS exhibited the smallest root mean square error (RMSE), endpoint error, and eastward error at the endpoint. The proposed algorithm improved navigation accuracy by 90.77 % compared to EKF, primarily due to the continuous GPS support. Additionally, EKF-VGPSM, leveraging the virtual GPS model, achieved a 69.56 % improvement in accuracy over EKF, further validating the effectiveness of the proposed method.
In the second scenario, illustrated in Figure 7, both EKF-GPS and the proposed algorithm were initially assisted by GPS for the first one-third of the experiment. However, during the remaining two-thirds of the time, the EKF-GPS lost its positional assistance, while the proposed algorithm continued to benefit from the virtual GPS model’s output. In the first phase, the errors of both EKF-GPS and the proposed algorithm were minimal due to GPS assistance. However, in the second phase, the error of EKF-GPS began to increase as it lost GPS support. In contrast, the proposed algorithm, aided by the virtual GPS model, maintained significantly lower errors and outperformed all other algorithms. According to Table 4, the proposed algorithm demonstrated the smallest RMSE, endpoint error, and errors in both the northward and eastward directions. The navigation accuracy of the proposed algorithm was improved by 37.77 % compared to EKF, showcasing its robustness even with partial GPS assistance.
In the third scenario, shown in Figure 8, no GPS assistance was available. Under these conditions, EKF-GPS lost its available GPS position data, and its navigation results aligned with those of the standard EKF, as indicated by the overlapping red and green trajectories in the figure. Both EKF-VGPSM and the proposed algorithm relied entirely on the virtual GPS model for positional data, resulting in consistent navigation outcomes, as reflected by the overlapping magenta and blue lines. This scenario effectively compares the navigation performance of EKF and EKF-VGPSM. The results clearly indicate that EKF-VGPSM and the proposed algorithm, both supported by the virtual GPS model, achieved significantly better navigation accuracy than the standard EKF and EKF-GPS. Table 4 confirms that the proposed algorithm reduced the root mean square error to 4.51 m compared to 6.37 m for EKF, resulting in a 29.2% improvement in navigation accuracy.
The effectiveness of the proposed method has been validated through a series of experiments under varying conditions. The results demonstrate that the proposed approach effectively reduces cumulative errors, with navigation accuracy improvements ranging from 29.2% to 69.56% compared to the conventional EKF. Additionally, the virtual GPS model’s output accurately predicts real GPS displacement changes, with root mean square errors of 0.0538 m northward and 0.0502 m eastward, respectively. These experimental results robustly demonstrate the proposed method’s effectiveness in improving navigation accuracy.

5.3. Error Analysis

The primary sources of error stem from sensor inaccuracies and model limitations. AHRS, while essential for orientation data, can introduce drift over time due to inherent biases and noise in gyroscope and accelerometer readings. DVL, used to measure velocity, can suffer from inaccuracies in low-quality bottom lock or turbulent conditions, leading to velocity estimation errors. GPS, although accurate, can introduce errors due to multipath effects or signal loss, especially in challenging underwater environments.
In the virtual GPS model, errors may arise from the integration of these noisy sensor data, where slight inaccuracies in velocity or orientation measurements can accumulate over time, leading to larger discrepancies. The proposed algorithm mitigates some of these errors through its innovative use of deep neural networks, which likely improves the sensor fusion process. However, the residual errors observed, such as the fluctuations in single-step errors and the slight biases in displacement predictions, indicate that the model could still be susceptible to overfitting or insufficient generalization in diverse conditions.
Overall, while the virtual GPS model and the proposed method significantly enhance navigation accuracy, particularly when GPS data are sparse or unavailable, the errors observed underscore the need for continued refinement in sensor fusion algorithms and model robustness to achieve even higher precision in underwater navigation.

5.4. Computational Complexity Analysis

The computational complexity of the extended Kalman filter primarily hinges on the matrix operations performed during the prediction and correction steps. Let n denote the dimension of the state vector, and o the dimension of the observation vector. The prediction step involves matrix multiplication and addition, leading to a complexity of O ( n 2 ) . The correction step requires matrix inversion and multiplication, with a complexity of O ( n 3 ) . Therefore, the overall complexity for one EKF iteration is O ( n 3 ) .
The virtual GPS model, which integrates Bi-LSTM and LSTM networks, has its computational complexity determined by the number of LSTM units, the sequence length, and the structure of the fully connected layers. Let d represent the input feature dimension, T the sequence length, u the number of LSTM units, and p and q the number of units in the fully connected layers. Given the bidirectional nature of Bi-LSTM, the complexity for this layer is O ( 2 m d u 2 ) . For the LSTM layer, processing the sequence unidirectionally yields a complexity of O ( m d u 2 ) . The fully connected layer computes the output from the LSTM units to the final regression layer, with a complexity of O ( u p + p q + 2 q ) . Thus, the total complexity for processing a sequence through VGPSM is O ( m d u 2 ) , assuming the fully connected layers and dropout do not significantly affect the complexity.
In practical AUV navigation applications, real-time performance is of paramount importance. The computational complexity of the EKF under standard conditions is O ( n 3 ) . When invoking the VGPSM every m time step, the complexity increases by O ( m d u 2 ) . While the introduction of VGPSM enhances navigational accuracy, it simultaneously imposes additional demands on computational resources. However, with the support of modern hardware, this level of computational complexity remains within an acceptable range. Specifically, by employing extended time steps m and leveraging efficient parallel computing frameworks, the method can meet the demands of most real-time AUV navigation tasks. Consequently, this approach successfully maintains high accuracy while still offering a significant real-time advantage.

6. Conclusions

To address the issue of increasing navigation errors in low-cost AUV without assisted positioning underwater, this paper introduces a virtual GPS model (VGPSM) based on deep sequence learning. By integrating VGPSM with the extended Kalman filter (EKF), the proposed approach provides a high-precision navigation solution for AUVs. The model leverages time-series data from sensors such as AHRS and DVL during surface navigation. It learns the relationship between these sensor inputs and GPS data using a hybrid architecture of LSTM and Bi-LSTM networks, which are particularly suited for processing and predicting time-series data. The VGPSM generates virtual GPS displacements updated at the same frequency as real GPS data. During underwater navigation, these virtual GPS displacements are utilized as measurements to assist the EKF in state estimation, significantly improving the accuracy and robustness of the navigation system.
However, this study presents certain limitations that warrant consideration. Firstly, the navigation datasets used for training and validation were obtained from a single AUV operating under favorable sea conditions. As a result, the generalizability of the model across different AUV platforms and in more challenging or variable marine environments remains uncertain. Additionally, the model’s computational efficiency requires further optimization to ensure its viability for deployment in real-world scenarios, where computational resources may be limited, and real-time processing demands are more stringent.
Future work will focus on addressing these limitations to enhance the robustness and applicability of the proposed method. Specifically, we aim to diversify the training datasets by incorporating data from multiple AUVs operating under varied environmental conditions, thereby improving the model’s adaptability and performance across diverse scenarios. Furthermore, we will explore advanced optimization techniques to reduce computational overhead, enabling the system to operate more efficiently across a broader range of hardware platforms.

Author Contributions

Conceptualization, P.-F.L.; methodology, P.-F.L.; software, P.-F.L.; validation, P.-F.L. and L.-X.X.; formal analysis, P.-F.L.; investigation, P.-F.L. and Z.-C.H.; resources, P.-F.L.; data curation, P.-F.L.; writing—original draft preparation, P.-F.L.; writing—review and editing, P.-F.L. and J.-Y.L.; visualization, P.-F.L. and Z.-C.H.; supervision, P.-F.L. and L.-X.X.; project administration, P.-F.L. and L.-X.X.; funding acquisition, P.-F.L. and L.-X.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Research of Jiangsu Higher Education Institutions of China (22KJB570003), Jiangsu Marine Technology Innovation Center (1012992301), and National Key Research and Development Program of China (2022YFC2806600, 2022YFC2806604).

Data Availability Statement

Data set available on request from the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhou, J.; Si, Y.; Chen, Y. A review of subsea AUV technology. J. Mar. Sci. Eng. 2023, 11, 1119. [Google Scholar] [CrossRef]
  2. Leonard, J.J.; Bahr, A. Autonomous underwater vehicle navigation. In Handbook of Ocean Engineering; Springer: Berlin/Heidelberg, Germany, 2016; pp. 341–358. [Google Scholar]
  3. Stutters, L.; Liu, H.; Tiltman, C.; Brown, D.J. Navigation technologies for autonomous underwater vehicles. IEEE Trans. Syst. Man, Cybern. Part C Appl. Rev. 2008, 38, 581–589. [Google Scholar] [CrossRef]
  4. Zhang, B.; Ji, D.; Liu, S.; Zhu, X.; Xu, W. Autonomous Underwater Vehicle navigation: A review. Ocean Eng. 2023, 273, 113861. [Google Scholar] [CrossRef]
  5. Maurelli, F.; Krupiński, S.; Xiang, X.; Petillot, Y. AUV localisation: A review of passive and active techniques. Int. J. Intell. Robot. Appl. 2022, 6, 246–269. [Google Scholar] [CrossRef]
  6. El-Sheimy, N.; Youssef, A. Inertial sensors technologies for navigation applications: State of the art and future trends. Satell. Navig. 2020, 1, 2. [Google Scholar] [CrossRef]
  7. Bian, H.; Li, A.; Ma, H.; Wang, R. Velocity Measurement. In Essentials of Navigation: A Guide for Marine Navigation; Springer: Berlin/Heidelberg, Germany, 2024; pp. 149–168. [Google Scholar]
  8. Zhang, L.; Wu, K.; Zhou, G.; Lan, Z. Design, Preparation, and Experiment of the Titanium Alloy Thin Film Pressure Sensor for Ocean Depth Measurements. IEEE Sens. J. 2024, 24, 14042–14049. [Google Scholar] [CrossRef]
  9. Kussat, N.H.; Chadwell, C.D.; Zimmerman, R. Absolute positioning of an autonomous underwater vehicle using GPS and acoustic measurements. IEEE J. Ocean. Eng. 2005, 30, 153–164. [Google Scholar] [CrossRef]
  10. Taraldsen, G.; Reinen, T.A.; Berg, T. The underwater GPS problem. In Proceedings of the Oceans 2011 IEEE, Santander, Spain, 6–9 June 2011; IEEE: Piscataway, NJ, USA, 2011; pp. 1–8. [Google Scholar]
  11. Zhang, L.; Hsu, L.T.; Zhang, T. A novel INS/USBL integrated navigation scheme via factor graph optimization. IEEE Trans. Veh. Technol. 2022, 71, 9239–9249. [Google Scholar] [CrossRef]
  12. Brossard, M.; Barrau, A.; Bonnabel, S. AI-IMU dead-reckoning. IEEE Trans. Intell. Veh. 2020, 5, 585–595. [Google Scholar] [CrossRef]
  13. Konrad, T.; Gehrt, J.J.; Lin, J.; Zweigel, R.; Abel, D. Advanced state estimation for navigation of automated vehicles. Annu. Rev. Control 2018, 46, 181–195. [Google Scholar] [CrossRef]
  14. Potokar, E.R.; Norman, K.; Mangelson, J.G. Invariant extended kalman filtering for underwater navigation. IEEE Robot. Autom. Lett. 2021, 6, 5792–5799. [Google Scholar] [CrossRef]
  15. Mahesh, B. Machine learning algorithms-a review. Int. J. Sci. Res. 2020, 9, 381–386. [Google Scholar] [CrossRef]
  16. Ma, D.; Chen, X.; Ma, W.; Zheng, H.; Qu, F. Neural network model-based reinforcement learning control for auv 3-d path following. IEEE Trans. Intell. Veh. 2023, 9, 893–904. [Google Scholar] [CrossRef]
  17. LeCun, Y.; Bengio, Y.; Hinton, G. Deep learning. Nature 2015, 521, 436–444. [Google Scholar] [CrossRef] [PubMed]
  18. Zhao, Z.; Chen, W.; Wu, X.; Chen, P.C.; Liu, J. LSTM network: A deep learning approach for short-term traffic forecast. IET Intell. Transp. Syst. 2017, 11, 68–75. [Google Scholar] [CrossRef]
  19. Hou, X.; Wang, J.; Bai, T.; Deng, Y.; Ren, Y.; Hanzo, L. Environment-aware AUV trajectory design and resource management for multi-tier underwater computing. IEEE J. Sel. Areas Commun. 2022, 41, 474–490. [Google Scholar] [CrossRef]
  20. Lv, P.F.; He, B.; Guo, J.; Shen, Y.; Yan, T.H.; Sha, Q.X. Underwater navigation methodology based on intelligent velocity model for standard AUV. Ocean Eng. 2020, 202, 107073. [Google Scholar] [CrossRef]
  21. Vandavasi, B.N.J.; Gidugu, A.R.; Venkataraman, H. Deep Learning Aided Magnetostatic Fields Based Real-Time Pose Estimation of AUV for Homing Applications. IEEE Sens. Lett. 2023, 7, 1–4. [Google Scholar] [CrossRef]
  22. Cohen, N.; Klein, I. BeamsNet: A data-driven approach enhancing Doppler velocity log measurements for autonomous underwater vehicle navigation. Eng. Appl. Artif. Intell. 2022, 114, 105216. [Google Scholar] [CrossRef]
  23. Qiao, Y.; Yin, J.; Wang, W.; Duarte, F.; Yang, J.; Ratti, C. Survey of deep learning for autonomous surface vehicles in marine environments. IEEE Trans. Intell. Transp. Syst. 2023, 24, 3678–3701. [Google Scholar] [CrossRef]
  24. He, Q.; Yu, H.; Fang, Y. Deep Learning-Based Inertial Navigation Technology for Autonomous Underwater Vehicle Long-Distance Navigation—A Review. Gyroscopy Navig. 2023, 14, 267–275. [Google Scholar]
  25. Xie, Y.X.; Liu, J.; Hu, C.Q.; Cui, J.H.; Xu, H. AUV dead-reckoning navigation based on neural network using a single accelerometer. In Proceedings of the 11th International Conference on Underwater Networks & Systems, Shanghai, China, 24–26 October 2016; pp. 1–5. [Google Scholar]
  26. Huang, H.; Chen, X.; Zhang, B.; Wang, J. High accuracy navigation information estimation for inertial system using the multi-model EKF fusing adams explicit formula applied to underwater gliders. ISA Trans. 2017, 66, 414–424. [Google Scholar] [CrossRef] [PubMed]
  27. Li, W.; Chen, M.; Zhang, C.; Zhang, L.; Chen, R. A Novel Neural Network-Based SINS/DVL Integrated Navigation Approach to Deal with DVL Malfunction for Underwater Vehicles. Math. Probl. Eng. 2020, 2020, 2891572. [Google Scholar] [CrossRef]
  28. Sabet, M.T.; Daniali, H.M.; Fathi, A.; Alizadeh, E. A low-cost dead reckoning navigation system for an AUV using a robust AHRS: Design and experimental analysis. IEEE J. Ocean. Eng. 2017, 43, 927–939. [Google Scholar] [CrossRef]
  29. Taud, H.; Mas, J. Multilayer Perceptron (MLP). In Geomatic Approaches for Modeling Land Change Scenarios; Camacho Olmedo, M.T., Paegelow, M., Mas, J.F., Escobar, F., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 451–455. [Google Scholar]
  30. Gonzalez, J.; Yu, W. Non-linear system modeling using LSTM neural networks. IFAC-PapersOnLine 2018, 51, 485–489. [Google Scholar] [CrossRef]
  31. Luo, Z.; Shi, D.; Ji, J.; Shen, X.; Gan, W.S. Real-time implementation and explainable AI analysis of delayless CNN-based selective fixed-filter active noise control. Mech. Syst. Signal Process. 2024, 214, 111364. [Google Scholar] [CrossRef]
  32. Xia, H.; Zhu, R.; Yuan, H.; Song, C. Rapid quantitative analysis of cotton-polyester blended fabrics using near-infrared spectroscopy combined with CNN-LSTM. Microchem. J. 2024, 200, 110391. [Google Scholar] [CrossRef]
  33. Topini, E.; Fanelli, F.; Topini, A.; Pebody, M.; Ridolfi, A.; Phillips, A.B.; Allotta, B. An experimental comparison of Deep Learning strategies for AUV navigation in DVL-denied environments. Ocean Eng. 2023, 274, 114034. [Google Scholar] [CrossRef]
  34. Prieto, A.; Prieto, B.; Ortigosa, E.M.; Ros, E.; Pelayo, F.; Ortega, J.; Rojas, I. Neural networks: An overview of early research, current frameworks and new challenges. Neurocomputing 2016, 214, 242–268. [Google Scholar] [CrossRef]
  35. Lv, P.F.; He, B.; Guo, J. Position correction model based on gated hybrid RNN for AUV navigation. IEEE Trans. Veh. Technol. 2021, 70, 5648–5657. [Google Scholar] [CrossRef]
  36. Lv, P.; Guo, J.; Song, Y.; Sha, Q.; Jiang, J.; Mu, X.; Yan, T.; He, B. Autonomous navigation based on iSAM and GPS filter for AUV. In Proceedings of the 2017 IEEE Underwater Technology (UT), Busan, Republic of Korea, 21–24 February 2017; pp. 1–4. [Google Scholar]
  37. Guo, J.; Li, D.; He, B. Intelligent collaborative navigation and control for AUV tracking. IEEE Trans. Ind. Inform. 2020, 17, 1732–1741. [Google Scholar] [CrossRef]
  38. Mu, X.; He, B.; Zhang, X.; Song, Y.; Shen, Y.; Feng, C. End-to-end navigation for autonomous underwater vehicle with hybrid recurrent neural networks. Ocean Eng. 2019, 194, 106602. [Google Scholar] [CrossRef]
  39. Staudemeyer, R.C.; Morris, E.R. Understanding LSTM–a tutorial into long short-term memory recurrent neural networks. arXiv 2019, arXiv:1909.09586. [Google Scholar]
  40. Greff, K.; Srivastava, R.K.; Koutník, J.; Steunebrink, B.R.; Schmidhuber, J. LSTM: A search space odyssey. IEEE Trans. Neural Networks Learn. Syst. 2016, 28, 2222–2232. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The coordinate system and onboard sensors for the AUV navigation system.
Figure 1. The coordinate system and onboard sensors for the AUV navigation system.
Drones 08 00441 g001
Figure 2. The framework of the proposed navigation method for AUVs.
Figure 2. The framework of the proposed navigation method for AUVs.
Drones 08 00441 g002
Figure 3. The flowchart of virtual GPS model based on deep sequence learning.
Figure 3. The flowchart of virtual GPS model based on deep sequence learning.
Drones 08 00441 g003
Figure 4. The structure of Bi-LSTM and LSTM. (a) Bi-LSTM. (b) LSTM.
Figure 4. The structure of Bi-LSTM and LSTM. (a) Bi-LSTM. (b) LSTM.
Drones 08 00441 g004
Figure 5. The result of virtual GPS model (a) Northward Displacement. (b) Eastward Displacement. (c) RMSE of Northward Displacement. (d) RMSE of Northward Displacement. (e,f) Error of Northward and Eastward Displacement.
Figure 5. The result of virtual GPS model (a) Northward Displacement. (b) Eastward Displacement. (c) RMSE of Northward Displacement. (d) RMSE of Northward Displacement. (e,f) Error of Northward and Eastward Displacement.
Drones 08 00441 g005
Figure 6. The comparison between proposed method and other algorithms in Scenario 1. (a) AUV trajectory. (b) Position RMSE. (c) Position error. (d) Northward and Eastward error.
Figure 6. The comparison between proposed method and other algorithms in Scenario 1. (a) AUV trajectory. (b) Position RMSE. (c) Position error. (d) Northward and Eastward error.
Drones 08 00441 g006
Figure 7. The comparison between proposed method and other algorithms in Scenario 2. (a) AUV trajectory. (b) Position RMSE. (c) Position error. (d) Northward and Eastward error.
Figure 7. The comparison between proposed method and other algorithms in Scenario 2. (a) AUV trajectory. (b) Position RMSE. (c) Position error. (d) Northward and Eastward error.
Drones 08 00441 g007
Figure 8. The comparison between proposed method and other algorithms in Scenario 3. (a) AUV trajectory. (b) Position RMSE. (c) Position error. (d) Northward and Eastward error.
Figure 8. The comparison between proposed method and other algorithms in Scenario 3. (a) AUV trajectory. (b) Position RMSE. (c) Position error. (d) Northward and Eastward error.
Drones 08 00441 g008
Table 1. Performance parameters and output data types of onboard sensors.
Table 1. Performance parameters and output data types of onboard sensors.
SensorsData TypeParameterValue
AHRS
(Ellipse-A)
w x b , w y b , w z b ,
a x b , a y b , a z b ,
ξ x b , ξ y b , ξ z b
Heading1° RMS
Pitch0.2° RMS
Roll0.2° RMS
Update Rate10 Hz (up to 200 Hz)
DVL
(NavQuest 600)
v x b , v y b , v z b ,
a l z b
Velocity Accuracy1 m/s
Altitude Range0.3–110 m
Maximum Velocity20 knot
Update Rate1 Hz
GPS
(LEA-M8T)
l a t , l o n Position Accuracy2.5 m CEP
AcquisitionGPS & BeiDou
Update Rate1 Hz
Cold Starts25 s
IPS
(Valeport)
p z n Accuracy0.01 % Range
Update Rate1 Hz
Table 2. Primary parameters employed in the experiments.
Table 2. Primary parameters employed in the experiments.
Parameter NameSymbolValueDescription
Time Step δ t 0.1Time interval for each time step.
Time Intervalm10Number of time steps in one GPS update cycle
Initial State Vector X 0 0 9 × 1 The initial 3D position, velocity, and attitude in navigation coordinate system.
Initial State Covariance P 0 0.1 I 3 , 0 3 × 6 0 3 × 3 , 0.01 I 3 , 0 3 × 3 0 3 × 6 , 0.001 I 3 Covariance matrix of initial state error.
Jacobian matrix 1 H 1 k 0 7 × 2 , I 7 Jacobian matrix of measurement function h 1 ( · )
Jacobian matrix 2 H 2 k I 9 Jacobian matrix of measurement function h 2 ( · )
Table 3. Numerical performance of virtual GPS model.
Table 3. Numerical performance of virtual GPS model.
Type|Error|-M (m)|Error|-S (m)Error-M (m)Error-S (m)Error-Med (m)RMSE (m)
Northward0.04330.0320.00190.05380.00210.0538
Displacement
Eastward0.03930.0312−0.00590.04980.00460.0502
Displacement
Table 4. Numerical performance of different algorithms for AUV navigation.
Table 4. Numerical performance of different algorithms for AUV navigation.
DataDistance (m)AlgorithmRMSE (m)Error (m)Error-N (m)Error-E (m)Accuracy
Scenario 1287.33EKF5.426.140.276.130.01886
EKF-VGPSM1.654.160.074.150.00574
EKF-GPS0.50.950.940.130.00174
Proposed0.50.950.940.130.00174
Scenario 2641.48EKF6.8310.756.048.90.01065
EKF-VGPSM6.069.54.458.390.00945
EKF-GPS4.457.765.225.740.00694
Proposed3.86.674.035.310.00592
Scenario 3457EKF6.377.355.065.330.01394
EKF-VGPSM4.514.943.313.670.00987
EKF-GPS6.377.355.065.330.01394
Proposed4.514.943.313.670.00987
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.

Share and Cite

MDPI and ACS Style

Lv, P.-F.; Lv, J.-Y.; Hong, Z.-C.; Xu, L.-X. Integration of Deep Sequence Learning-Based Virtual GPS Model and EKF for AUV Navigation. Drones 2024, 8, 441. https://doi.org/10.3390/drones8090441

AMA Style

Lv P-F, Lv J-Y, Hong Z-C, Xu L-X. Integration of Deep Sequence Learning-Based Virtual GPS Model and EKF for AUV Navigation. Drones. 2024; 8(9):441. https://doi.org/10.3390/drones8090441

Chicago/Turabian Style

Lv, Peng-Fei, Jun-Yi Lv, Zhi-Chao Hong, and Li-Xin Xu. 2024. "Integration of Deep Sequence Learning-Based Virtual GPS Model and EKF for AUV Navigation" Drones 8, no. 9: 441. https://doi.org/10.3390/drones8090441

APA Style

Lv, P. -F., Lv, J. -Y., Hong, Z. -C., & Xu, L. -X. (2024). Integration of Deep Sequence Learning-Based Virtual GPS Model and EKF for AUV Navigation. Drones, 8(9), 441. https://doi.org/10.3390/drones8090441

Article Metrics

Back to TopTop