Next Article in Journal
Assessment of Machine Learning Techniques for Oil Rig Classification in C-Band SAR Images
Next Article in Special Issue
Performance of Multi-GNSS in the Asia-Pacific Region: Signal Quality, Broadcast Ephemeris and Precise Point Positioning (PPP)
Previous Article in Journal
Fast Approach for SAR Imaging of Ground-Based Moving Targets Based on Range Azimuth Joint Processing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Assessment of BDS-3 PPP-B2b/INS Loosely Coupled Integration

1
College of Oceanography and Space Informatics, China University of Petroleum, Qingdao 266580, China
2
Shengli Branch, Sinopec Geophysical Corporation, Dongying 257100, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(13), 2957; https://0-doi-org.brum.beds.ac.uk/10.3390/rs14132957
Submission received: 27 May 2022 / Revised: 15 June 2022 / Accepted: 19 June 2022 / Published: 21 June 2022
(This article belongs to the Special Issue Precise Point Positioning with GPS, GLONASS, BeiDou, and Galileo)

Abstract

:
The BeiDou global navigation satellite system (BDS-3) has been officially providing a real-time precise point positioning (PPP) augmentation service, known as the PPP-B2b service, since 2020. Decimeter-level positioning accuracy is expected to be achieved based on the PPP-B2b service. It shows great potential for global navigation satellite system (GNSS) real-time applications, including, for example, vehicle positioning on land. However, the application of the PPP-B2b service is still full of challenges in the urban environment because of GNSS signal blockage. The inertial navigation system (INS) is a popular technology which can provide continuous positions under GNSS challenging scenarios. In this study, we constructed a BDS-3 PPP-B2b/INS loosely coupled integration system for vehicle positioning and evaluated its performance through two automotive experiments. In the first experiment, four periods of 30 s GNSS outages were simulated to evaluate the performance of PPP-B2b/INS loosely coupled integration during GNSS outages. During the simulated GNSS outages, PPP-B2b positioning did not work. Nevertheless, PPP-B2b/INS loosely coupled integration provided continues solution through INS mechanization. The averaged positioning errors at the last epoch of outages were 300.6/498.0/41.0 cm for PPP-B2b/MEMS-IMU and 18.6/21.8/6.1 cm for PPP-B2b/Tactical-IMU loosely coupled integration, in the east, north and up directions, respectively. In the second experiment, we drove the land vehicle in a complex urban environment for 15 min. During this period, two GNSS signal interruptions occurred due to the occlusion of bridges, lasting 15 s and 5 s, respectively. The results show that the improvement of positioning accuracy in the east, north, and up components were 64.1%, 77.8%, and 73.8% respectively for PPP-B2b/MEMS-IMU loosely coupled integration, and 63.9%, 79.5%, and 74.4% respectively for PPP-B2b/Tactical-IMU loosely coupled integration, as compared to the positioning accuracy of PPP-B2b only.

1. Introduction

Real-time kinematic (RTK) and precise point positioning (PPP) are usually conducted to satisfy the high-precision positioning requirements of autonomous cars and unmanned aerial vehicles. RTK technology can archive centimeter-level positioning accuracy with short initialization time [1,2], but it needs regional corrections from a reference station or reference network. The concept of PPP technology was proposed in 1997 [3,4], and it could obtain decimeter- even to centimeter-level positioning accuracy using standalone global navigation satellite system (GNSS) equipment [5]. Over the decades, PPP has become a widely used high-precision GNSS positioning technology. Precise GNSS satellite products are indispensable for the implementation of PPP. Motivated by the requirements of GNSS real-time applications, the international GNSS service (IGS) has been providing real-time service (RTS) through Internet communication since 2013 [6,7]. The RTS corrections, including precise orbit correction, clock offset correction, and code biases, are sent to users based on the Networked Transport of RTCM via Internet Protocol (NTRIP) [8]. By applying these corrections, the users can compute precise satellite orbits and clocks, and then carry out real-time PPP. However, RTS may be interrupted when the Internet is unavailable [9].
The commissioning of the BeiDou global navigation satellite system (BDS-3) was announced on 31 July 2020 [10,11]. Several featured services are provided by BDS-3, including short message communication, international search and rescue, and the PPP-B2b service, in addition to global positioning, navigation and time (PNT) services [12]. The BDS-3 PPP-B2b service can support real-time PPP based on PPP-B2b corrections broadcast by 3 BDS-3 geostationary earth orbit (GEO) satellites [13,14]. The PPP-B2b corrections include satellite orbit correction, clock correction, and code bias correction. At present, the BDS-3 PPP-B2b service only provides corrections for BDS-3/GPS satellites. The BDS-3 PPP-B2b service can be obtained free of charge. In addition, unlike IGS-RTS, Internet communication is not required since the BDS-3 PPP-B2b service is a satellite-based service. Therefore, the PPP-B2b service has great potential for real-time GNSS applications. Recently, publications have presented the precision of PPP-B2b corrections and positioning performance based on the PPP-B2b service [15,16,17,18,19]. These results indicate that the signal-in-space accuracy of precise ephemeris calculated with PPP-B2b corrections is at the decimeter- even to centimeter-level. In term of positioning performance, centimeter-level accuracy in the static mode and decimeter-level in the kinematic mode could be achieved. However, positioning performance based on the BDS-3 PPP-B2b service would be significantly degraded in severe environments, due to GNSS signal blockages. To overcome this weakness, GNSS should be integrated with other sensors [20]. The inertial navigation system (INS) is an autonomous and spontaneous navigation system, which has potential to overcome degradation of GNSS positioning performance when GNSS signal outage occurs. Many valuable studies in PPP/INS integration have been published in the past decades, including loosely coupled (LC) model and tightly coupled (TC) integration [21,22,23,24,25]. In these publications, the precise ephemeris derived from IGS final precise products or from the IGS RTS service is adopted for PPP/INS integration processing. However, there is a considerable delay obtaining IGS final precise products. As for the IGS RTS service, the PPP/INS integration is restricted when there is a lack of Internet communication. It is worth mentioned that the above disadvantages can be overcome when the BDS-3 PPP-B2b service has implemented PPP/INS integration.
To our knowledge, there has until now been no published research on BDS-3 PPP-B2b/INS integration. In this paper, we focus on performance assessment of BDS-3 PPP-B2b/INS loosely coupled integration. Two vehicle kinematic experiments were carried out to assess the performance of PPP-B2b/INS loosely coupled integration. The rest of the paper is organized as follows. In Section 2, the methodology is introduced, including PPP positioning based on PPP-B2b corrections, the INS model, and their loosely coupled integration. In Section 3, the experimental set-ups and data processing strategies are presented in detail. Section 4 presents the performance of PPP-B2b/INS loosely coupled integration in urban environments. Finally, conclusions and perspectives are illustrated in Section 5.

2. Methodology

BDS-3 PPP-B2b/INS loosely coupled integration can be divided into three parts, including PPP positioning based on the PPP-B2b service, the INS model, and the integration of BDS-3 PPP-B2b and INS. In this section, these three parts are described in detail.

2.1. PPP Positioning Based on PPP-B2b Service

As is known, precise satellite orbit and clock products play a key role in PPP processing. For PPP-B2b service, the satellite orbit corrections and clock offset corrections are broadcast by BDS-3 GEO satellites. Up to now, the PPP-B2b service has provided corrections only for BDS-3/GPS satellites. According to the PPP-B2b Interface Control Document (ICD) from China Satellite Navigation Office (CSNO) [26], the PPP-B2b orbit corrections are given in the radial ( δ O r a d i a l ) , along-track ( δ O a l o n g ) , and cross-track ( δ O c r o s s ) directions. Hence, the orbit correction vector δ O B 2 b = [ δ O r a d i a l δ O a l o n g δ O c r o s s ] T should first be transformed to the Earth-Center Earth-Fixed (ECEF) frame. This is because the satellite positions derived from CNAV1/LNAV broadcast ephemeris are based on the ECEF frame. This transformation can be described as follows:
δ X s a t = [ e r a d i a l e a l o n g e c r o s s ] · δ O B 2 b
with
{ e r a d i a l = r | r | e c r o s s = r × r ˙ | r × r ˙ | e a l o n g = e c r o s s × e r a d i a l
where δ X s a t = [ δ O x δ O y δ O z ] T represents the PPP-B2b orbit correction vector in the ECEF frame. r and r ˙ are the satellite position and velocity vectors derived from broadcast ephemeris. Then, the precise satellite position vector X p r e c can be computed by
X p r e c = X b r d c δ X s a t
where X b r d c is the satellite position vector derived from broadcast ephemeris. As described in ICD document, the PPP-B2b clock correction parameter is defined as offset to the broadcast ephemeris clock in meters. The precise satellite clock offset is given by
d t p r e c s a t = d t b r d c s a t C 0 c
where C 0 represents the PPP-B2b clock correction parameter; d t b r d c s a t is the satellite clock offset derived from broadcast ephemeris; d t p r e c s a t denotes the precise PPP-B2b clock offset; c is the velocity of light in a vacuum.
The GNSS raw code and carrier-phase measurements can read as [27,28]
{ P i = ρ + c · ( d t r d t s ) + T + f 1 2 f i 2 · I 1 + B r , i B i s + ε P i L i = ρ + c · ( d t r d t s ) + T f 1 2 f i 2 · I 1 + λ i ( N i + b r , i b i s ) + ε L i
where i represents the frequency number; P i and L i are raw code and phase measurements; ρ denotes the geometric distance from satellite to receiver; d t r and d t s are the clock offsets at the receiver and satellite, respectively; f i is the frequency value, T is the tropospheric delay, and I 1 denotes the ionospheric delay for L 1 ; B r , i and B i s represent the code bias at the receiver/satellite end; b r , i and b i s denote the phase bias at the receiver/satellite end; ε P i and ε L i are the unmodelled errors of code/phase measurements. The relativistic effect, Sagnac effect, Shapiro time delay [29], site displacements [30], and phase windup [31] should be corrected according to the corresponding models.
The ionosphere-free code and phase combinations are usually adopted by PPP to eliminate the first-order ionospheric delay. The ionospheric-free code and phase combinations read as [32]:
{ P I F = α · P i + ( 1 α ) · P j = ρ + c · d t ¯ r c · d t s B I F s + T + ε P I F L I F = α · L i + ( 1 α ) · L j = ρ + c · d t ¯ r c · d t s + T + λ I F · N ¯ I F + ε L I F
where α = f i 2 / ( f i 2 f j 2 ) and λ I F · N ¯ I F = λ I F N I F b I F s + b r , I F B r , I F , b I F s and b r , I F are phase bias of ionospheric-free phase combination at the satellite/receiver end, B I F s and B r , I F represent code bias of ionospheric-free code combination at the satellite/receiver end, λ I F and N I F are the ionospheric-free wavelength and the ionospheric-free ambiguity, respectively; d t ¯ r is the recombined receiver clock offset, which absorbs the receiver code bias of ionospheric-free code combination; ε P I F and ε L I F denote the unmodelled errors of ionospheric-free code and phase combinations. It is noted that the ionospheric-free code bias at the satellite end should be corrected by applying PPP-B2b differential code bias corrections [15].
When the recovered PPP-B2b precise satellite orbits and clock offsets have been applied, the satellite orbit and clock errors are considered eliminated. The tropospheric delay can be divided into the dry and wet parts. The Saastamoinen model is usually used to correct the dry part, and the wet part must be estimated as unknown. Then the ionospheric-free code and phase combinations can be linearized as
{ p I F = e · x + c · d t ¯ r + M W · z w d + ε P I F l I F = e · x + c · d t ¯ r + M W · z w d + λ I F · N ¯ I F + ε L I F
where e , x represent the unit vector from receiver to satellite, and the vector of position increment, respectively; z w d denotes the zenith wet delay and M W is the corresponding mapping function. In this equation, the remaining unknown parameters include only the receiver position increment vector x , the receiver clock offset d t ¯ r , the zenith wet delay z w d and the ionosphere-free ambiguity λ I F N ¯ I F . By adopting a Kalman filter, the unknown parameters can be exactly estimated. It should be noted that in this paper the ionospheric-free Doppler combination is used to derive the velocity [23].

2.2. INS Model

The mechanization of INS in the navigation frame ( n -frame) can be expressed as an integral process for the following equation [33,34]:
[ p ˙ I N S n v ˙ I N S n C ˙ b n ] = [ D 1 v I N S n C b n f b ( 2 ω i e n + ω e n n ) × v I N S n + g n C b n ( ω i b b × ) ( ω i e n + ω e n n ) × C b n ]
where D 1 is the diagonal matrix to transform the rectangular coordinates into geodetic coordinates; p I N S n and v I N S n denote the INS position vector and INS velocity vector in n -frame; C b n represents the transformation matrix from body-frame ( b -frame) to n -frame; f b is the inertial measurement unit (IMU) specific force measurement projected in b -frame; ω i b b denotes the IMU angular rate measurement expressed in b -frame; ω i e n is the vector of earth rotation rate expressed in n -frame; g n is the gravity vector presented in n -frame; ω e n n represents the rotation rate of ECEF frame relative to that of n -frame expressed in n -frame; The symbol “ × ” represents the cross-product operator. Obviously, the INS update of position, velocity, and attitude can be implement based on Equation (8).
The perturbation of position, velocity and attitude must be considered in the INS data processing. In this paper, the error model of attitude, velocity, and position is defined as the psi-angle error model, which reads as
[ δ p ˙ I N S n δ v ˙ I N S n ϕ ˙ ] = [ ω e n n × δ p I N S n + δ v I N S n ( 2 ω i e n + ω e n n ) × δ v I N S n + C b n f b × ϕ + δ g n + C b n δ f b ( ω i e n + ω e n n ) × ϕ C b n δ ω i b b ]
where δ p I N S n and δ v I N S n are the error vector of INS position and velocity, respectively; ϕ indicates the vector of misalignment angles; δ g n is the gravity error; δ f b is error vector of specific force measurement; δ ω i b b denotes error vector of the IMU angular rate. Generally, the error vectors of specific force and angular rate measurement include scale factor errors, bias errors, and white noise [35]. Here, only bias errors are considered and the bias errors of IMU sensors are described as the first Gauss-Markov procedure [34]:
[ δ B ˙ a δ B ˙ g ] = 1 τ [ δ B a δ B g ] + [ w a w g ]
where δ B g and δ B a represent the error vector of gyro biases and accelerometer biases, respectively; τ is the correction time; and w a , w g represents the driving noise for accelerometer and gyro, respectively.

2.3. PPP-B2b/INS Loosely Coupled Integration Model

A typical state vector of 15 states is used in PPP-B2b/INS loosely coupled integration, which can be expressed as
X k = [ δ p I N S n δ v I N S n ϕ δ B a δ B g ] T
The system model of PPP-B2b/INS loosely coupled integration in the discrete form can be simplified as
X k = ( I + F · Δ t ) X k 1 + Γ k 1 w k 1 , w k 1 ~ N ( 0 , Q k 1 )
where X k 1 and X k are the vector of state parameter at epoch k 1 and k , respectively; I is the identify matrix; Δ t denotes the time interval between two adjacent epochs; w k 1 represents the vector of the process noise at epoch k 1 , and Γ k 1 denotes the matrix of the noise distribution; Q k 1 represents the covariance matrix of process noise; F is the dynamic matrix, which can be derived from the INS error model as follows:
F = [ F r r F r v 0 0 0 F v r F v v ( f b × ) C b n 0 F ϕ r F ϕ v ( ω i e n + ω e n n × ) 0 C b n 0 0 0 1 τ I 0 0 0 0 0 1 τ I ]
where 0 denotes the zero matrix; F r r , F r v , F v r , F v v , F ϕ r , and F ϕ v are sub-matrices of dynamic matrix, which can be derived from Equation (9). Then the Kalman time update can be expressed as [36]
{ X k = Φ k , k 1 X k 1 + P k = Φ k , k 1 P k 1 + Φ k , k 1 T + Q k 1
where the superscript “ + ” and “ ” denote the estimated and predicted information; Φ k , k 1 = ( I + F · Δ t ) represents the state transition matrix; P k represents the covariance matrix of the states.
When the PPP-B2b positioning solution and INS predicted solution are available at the same epoch, the Kalman filter measurement update can be operated. However, the positions and velocities from INS are referred to the IMU center, while those from PPP-B2b positioning are typically based on the GNSS antenna phase center. Therefore, the corresponding lever-arm offsets must be compensated. Here, the accurate lever-arm offsets are measured and directly applied to the position and velocity from INS. Then, the observation model for PPP-B2b/INS loosely coupled integration can be expressed as
Z k = H k X k + [ ε p ε v ] , [ ε p ε v ] ~ N ( 0 , R k )
with
Z k = [ p G N S S n v G N S S n ] [ p I N S n v I N S n ]
H k = [ I 3 × 3 0 3 × 3 ( C b n l b × ) 0 3 × 3 0 3 × 3 0 3 × 3 I 3 × 3 H ϕ 0 3 × 3 C b n ( l b × ) ]
where H ϕ = [ ( ω i e n + ω e n n ) × C b n ( l b × ) + C b n ( l b × ω i b b ) × ] ; p G N S S n , v G N S S n denote the vector of position and velocity obtained from PPP-B2b positioning; [ ε p ε v ] T denotes the vector of measurement noise of position and velocity and R k represents its related covariance matrix; l b represents the vector of lever-arm offsets; H k is the design matrix; Z k is defined as the innovation vector, which consists of the difference of positions and velocities between PPP-B2b positioning solution and INS predicted solution. Finally, the Kalman measurement update can be given as
{ K k = P k H k T ( H k P k H k T + R k ) 1 X k + = X k + K k ( Z k H k X k ) P k + = ( I K k H k ) P k
where K k is the Kalman gain.

3. Description of Experiments

In order to assess the performance of BDS-3 PPP-B2b/INS loosely coupled integration, two land vehicle experiments in open-sky scenery and urban canyon scenery (hereinafter referred to as Experiment A and Experiment B, respectively) were carried out in Qingdao, China. In these two experiments, a tactical-grade IMU (ISA100C) [37] and a MEMS-grade IMU (ADIS-16505) [38] were equipped to evaluate the impact of IMU grade on BDS-3 PPP-B2b/INS loosely coupled integration. Two GNSS antennas were mounted on the roof of the land vehicle, and connected to a NovAtel PwrPak7 receiver [39] and a FRII-PLUS PPP-B2b receiver (http://www.femtomes.com, accessed on 20 October 2021). The FRII-PLUS receiver was only used to collected PPP-B2b messages, while all GNSS observations used in the evaluation were collected by the PwrPak7 receiver. The sampling rate of GNSS observations was 1 Hz. Figure 1 shows the installation of experimental equipment. Relevant specific information for the IMU sensors is listed in Table 1.
During these two experiments, a GNSS base station was established on the roof of Engineering-C building in the China University of Petroleum (East China). It should be noted that the distance between the base and rover stations was less than 5 km. Using GNSS observations collected by PwrPak7, and IMU data from ISA100C, together with GNSS measurements at the base station, the smoothed RTK/INS tightly coupled solution was obtained by commercial Inertial Explorer software (IE 8.90) from the NovAtel company, and this solution was used as reference.
In PPP-B2b positioning, the B1C and B2a signals of BDS-3, L1 and L2 signals of GPS were selected to form the ionospheric-free combination, and the cut-off elevation angle was set to 7 degrees. The remaining processing strategies are shown in detail in Table 2.

4. Result and Discussion

In this section, the performance of BDS-3 PPP-B2b/INS loosely coupled integration was investigated in both open-sky and complex urban environments.

4.1. Experiment A

Experiment A was conducted from 16:30:40 to 16:59:59 on 3 December 2021 in GPS time. The trajectory of the land vehicle is shown in Figure 2.
During Experiment A, the velocities of the land vehicle were within ± 20 m/s in the east and north directions and within ± 1 m/s in the vertical direction. The number of visible BDS-3/GPS satellites and the position dilution of precision (PDOP) values are shown in Figure 3. The number of visible satellites ranged from 9 to 15 with an average value of 12.5, and the PDOP varied from 1.28 to 2.70 with an average value of 1.59.
Figure 4 shows the positioning errors of PPP-B2b only, PPP-B2b/MEMS-IMU loosely coupled integration, and PPP-B2b/tactical-IMU loosely coupled integration, respectively. As with traditional PPP processing, a convergence period is indispensable in PPP-B2b positioning and PPP-B2b/INS loosely coupled integration. In this study, the convergence time was defined for horizontal/vertical positioning accuracy better than 30 cm/60 cm, with such a positioning accuracy for at least 60 continuous epochs. The statistics of the positioning results are summarized in Table 3. It should be noted that positioning results before full convergence were excluded when calculating positioning accuracy. The biases in the east, north and up directions were 6.7/23.7/26.5 cm for PPP-B2b only, and the corresponding values for PPP-B2b/tactical-IMU and PPP-B2b/MEMS-IMU loosely coupled integration schemes were 6.6/23.3/26.0 cm and 6.2/23.5/25.8 cm, respectively. In terms of the root mean square (RMS) value of positioning errors, the improvement of PPP-B2b/INS loosely coupled integration was not significant compared to PPP-B2b only. This is because Experiment A was carried out within an open-sky environment and thus PPP-B2b on its own was already able to obtain high-precision positioning accuracy.
It is recognized that PPP-B2b only is unable to obtain a high-precision solution during GNSS signal outages, while PPP-B2b/INS loosely coupled integration continues to work by INS mechanization. Four periods of GNSS signal outage were simulated to present the superiority of PPP-B2b/INS loosely coupled integration. The periods of GNSS outage were simulated from the 600th/900th/1200th/1500th epoch, and each period lasted 30 s. The corresponding positioning errors of PPP-B2b/MEMS-IMU and PPP-B2b/tactical-IMU loosely coupled integration are plotted in Figure 5 and Figure 6, respectively. Clearly, the positioning accuracy become progressively worse with increased time of GNSS outages in both PPP-B2b/MEMS-IMU and PPP-B2b/tactical-IMU loosely coupled integration schemes. When the GNSS signal outage time increased to 30 s, the positioning errors in the east, north, and up directions decreased to 300.0 cm, 498.0 cm, and 41.0 cm, respectively, for PPP-B2b/MEMS-IMU loosely coupled integration. The results for the PPP-B2b/tactical-IMU loosely coupled integration scheme dropped to 18.6 cm, 21.8 cm, and 6.1 cm in the east, north, and up components, respectively. Table 4 lists the statistics results, showing that the positioning performance of PPP-B2b/tactical-IMU loosely coupled integration is significantly better than that of PPP-B2b/MEMS-IMU loosely coupled integration.

4.2. Experiment B

To fully demonstrate the performance of BDS-3 PPP-B2b/INS loosely coupled integration in an urban scenario, we conducted Experiment B in a complex urban environment. Experiment B was carried out from 09:00:00 to 09:49:59 on 5 December 2021 in GPS time. The trajectory of the land vehicle is shown in Figure 7.
Figure 8 presents the number of visible BDS-3/GPS satellites and the values of PDOP. The number of available BDS-3/GPS satellites dropped frequently, and at the same time the PDOP value significantly increased as several typical city features such as tall buildings, trees, and over-bridges appeared along the vehicle route.
During the first 35 min, there were enough satellites for the implementation of PPP-B2b positioning. The positioning errors of PPP-B2b only, PPP-B2b/MEMS-IMU loosely coupled integration and PPP-B2b/tactical-IMU loosely coupled integration are plotted in Figure 9. It can be found that all three schemes provided similar positioning accuracy in this 35-min period. The RMS values of positioning errors were 20.7 cm, 17.3 cm, and 47.6 cm in the east, north, and up directions for PPP-B2b only. For PPP-B2b/INS loosely coupled integration, the corresponding RMS values were 20.8 cm, 17.4 cm, and 47.5 cm for the PPP-B2b/MEMS-IMU scheme, and 20.3 cm, 17.4 cm, and 47.5 cm for the PPP-B2b/tactical-IMU scheme. It should be noted that a convergence period of 15 min was removed when we computed the positioning accuracy. Clearly, the equipped IMU sensors could not improve the performance of the PPP-B2b service when there was no GNSS signal blockage.
From 09:35:00 to the end of Experiment B, the land vehicle crossed two over-bridges, as shown in Figure 10. Two GNSS signal outages appeared when the land vehicle crossed the bridges, lasting 15 s and 5 s, respectively. During these two GNSS signal outages, PPP-B2b positioning could be implemented, because the number of visible satellites was less than four. Therefore, the PPP-B2b positioning was re-converged after crossing the bridges. Figure 11 shows the positioning errors of PPP-B2b only, PPP-B2b/MEMS-IMU loosely coupled integration, and PPP-B2b/tactical-IMU loosely coupled integration. It can be seen that the positioning accuracy of PPP-B2b/INS loosely coupled integration exhibited better performance than PPP-B2b positioning during the re-convergence phase. The overall positioning accuracies for the three schemes were calculated from 09:35:00 to 09:49:59 and are presented in Table 5. The RMS values of positioning errors for the east, north, and up components were 114.2 cm, 200.0 cm, and 244.3 cm for PPP-B2b only. These figures improved to 41.0 cm, 44.5 cm, and 64.0 cm for PPP-B2b/MEMS-IMU loosely coupled integration, and to 41.2 cm, 41.0 cm, and 62.5 cm for PPP-B2b/tactical-IMU loosely coupled integration. Obviously, the performance of PPP-B2b/INS loosely coupled integration was significantly improved compared to PPP-B2b only in the last 15-min period of Experiment B. The positioning accuracy of PPP-B2b/tactical-IMU loosely coupled integration was slightly better than that of PPP-B2b/MEMS-IMU loosely coupled integration.

5. Conclusions

Since July 2020, PPP-B2b has served as a featured service of BDS-3 to support real-time PPP. Compared to the IGS RTS service, PPP-B2b is a satellite-based service, which is not limited by Internet communication. Its high-precision positioning performance shows great potential in real-time GNSS applications. However, the BDS-3 PPP-B2b service continues to encounter challenges in urban environments due to GNSS signal outages. In order to overcome the drawbacks of the BDS-3 PPP-B2b service in urban environments, we set up a BDS-3 PPP-B2b/INS loosely coupled integration system and evaluated its performance in different urban scenarios.
The experimental results indicate that PPP-B2b/INS loosely coupled integration cannot improve the positioning performance in open-sky environments where there are enough GPS/BDS-3 satellites for PPP-B2b positioning. However, PPP-B2b/INS loosely coupled integration can show its superiority when GNSS signal outage occurs. During GNSS signal outages, the INS mechanization can provide continuous positioning. Therefore, the performance of PPP-B2b positioning can be improved by adopting PPP-B2b/INS loosely coupled integration in GNSS blockage environments. In the first experiment, we simulated four 30-s periods of GNSS signal outages. At the last epoch of the simulated outages, the averaged positioning errors in the east, north, and up directions were 300.0 cm, 498.0 cm, and 41.0 cm, respectively, for PPP-B2b/MEMS-IMU loosely coupled integration, and 18.6 cm, 21.8 cm, and 6.1 cm, respectively, for PPP-B2b/tactical-IMU loosely coupled integration. In addition, we also evaluated the performance of PPP-B2b/INS loosely coupled integration in a real complex urban environment. When the land vehicle crossed bridges, two GNSS signal interruptions appeared that lasted 15 s and 5 s, respectively. Compared to PPP-B2b only, the respective improvement of positioning accuracy in the east, north, and up components was 64.1%, 77.8%, and 73.8% for PPP-B2b/MEMS-IMU loosely coupled integration, and 63.9%, 79.5%, and 74.4% for PPP-B2b/Tactical-IMU loosely coupled integration.
It can be concluded that when there are sufficient valid satellites for PPP-B2b positioning, the positioning accuracy of PPP-B2b/INS loosely coupled integration depends mainly on the accuracy of PPP-B2b. In urban environments, positioning performance can be significantly improved compared to PPP–B2b alone, by adopting PPP-B2b/INS loosely coupled integration especially in GNSS signal blockage environments.

Author Contributions

X.X.; Conceptualization, methodology, software, visualization, and writing—original draft preparation. Z.N.; Conceptualization, software, writing—review and editing, and funding acquisition. Z.W.; writing—review and editing, supervision, project administration, and funding acquisition. B.W.; validation and formal analysis. Q.D.; validation and formal analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (No. 42104011), the Shandong Provincial Natural Science Foundation (No. ZR2021QD069), the National Key Research and Development Program of China (No. 2019YFC1509205), and Sinopec Science and Technology Entry Program (JP210004).

Data Availability Statement

The datasets used in this study are managed by the College of Oceanography and Space Informatics, China University of Petroleum, Qingdao, China and can be available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bock, Y.; Nikolaidis, R.M.; de Jonge, P.J.; Bevis, M. Instantaneous geodetic positioning at medium distances with the Global Positioning System. J. Geophys. Res. Solid Earth 2020, 105, 28223–28253. [Google Scholar] [CrossRef]
  2. Rizos, C. Network RTK research and implementation-a geodetic perspective. J. Glob. Posit. Syst. 2009, 1, 144–151. [Google Scholar] [CrossRef] [Green Version]
  3. Kouba, J.; Héroux, P. Precise point positioning using IGS orbit and clock products. GPS Solut. 2001, 5, 12–28. [Google Scholar] [CrossRef]
  4. Zumberge, J.; Heflin, M.; Jefferson, D.; Watkins, M.; Webb, F. Precise point positioning for the efficient and robust analysis of GPS data from large networks. J. Geophys. Res. Solid Earth 1997, 102, 5005–5017. [Google Scholar] [CrossRef] [Green Version]
  5. Teunissen, P.; Montenbruck, O. Springer Handbook of Global Navigation Satellite Systems; Springer International Publishing: New York, NY, USA, 2017. [Google Scholar]
  6. Hadas, T.; Bosy, J. IGS RTS precise orbits and clocks verification and quality degradation over time. GPS Solut. 2015, 19, 93–105. [Google Scholar] [CrossRef] [Green Version]
  7. Nie, Z.; Gao, Y.; Wang, Z.; Ji, S.; Yang, H. An approach to GPS clock prediction for real-time PPP during outages of RTS stream. GPS Solut. 2018, 22, 14. [Google Scholar] [CrossRef]
  8. Elsobeiey, M.; Al-Harbi, S. Performance of real-time Precise Point Positioning using IGS real-time service. GPS Solut. 2016, 20, 565–571. [Google Scholar] [CrossRef]
  9. Nie, Z.; Wang, B.; Wang, Z.; He, K. An offshore real-time precise point positioning technique based on a single set of BeiDou short-message communication devices. J. Geod. 2020, 94, 1–11. [Google Scholar] [CrossRef]
  10. Cai, H.; Meng, Y.; Geng, C.; Gao, W.; Zhang, T.; Li, G.; Shao, B.; Xin, J.; Lu, H.; Mao, Y. BDS-3 performance assessment: PNT, SBAS, PPP, SMC and SAR. Acta Geod. Cartogr. Sinica. 2021, 50, 427–435. [Google Scholar]
  11. Yang, Y.; Liu, L.; Li, J.; Yang, Y.; Zhang, T.; Mao, Y.; Sun, B.; Ren, X. Featured services and performance of BDS-3. Sci. Bull. 2021, 66, 2135–2143. [Google Scholar] [CrossRef]
  12. China Satellite Navigation Office (CSNO). BeiDou Navigation Satellite System Open Service Performance Standard (Version 2.0). Available online: http://www.beidou.gov.cn/xt/gfxz/201812/P020181227529449178798.pdf (accessed on 30 September 2020).
  13. Liu, C.; Gao, W.; Liu, T.; Wang, D.; Yao, Z.; Gao, Y.; Nie, X.; Wang, W.; Li, D.; Zhang, W.; et al. Design and implementation of a BDS precise point positioning service. Navig. J. Instit. Navig. 2020, 67, 875–891. [Google Scholar] [CrossRef]
  14. Yang, Y.; Gao, W.; Guo, S.; Mao, Y.; Yang, Y. Introduction to BeiDou-3 navigation satellite system. Navigation 2019, 66, 7–18. [Google Scholar] [CrossRef] [Green Version]
  15. Nie, Z.; Xu, X.; Wang, Z.; Du, J. Initial Assessment of BDS PPP-B2b Service: Precision of Orbit and Clock Corrections, and PPP Performance. Remote Sens. 2021, 13, 2050. [Google Scholar] [CrossRef]
  16. Zhang, W.; Lou, Y.; Song, W.; Sun, W.; Zou, X.; Gong, X. Initial assessment of BDS-3 precise point positioning service on GEO B2b signal. Adv. Space Res. 2021, 69, 690–700. [Google Scholar] [CrossRef]
  17. Ren, Z.; Gong, H.; Peng, J.; Tang, C.; Huang, X.; Sun, G. Performance assessment of real-time precise point positioning using BDS PPP-B2b service signal. Adv. Space Res. 2021, 68, 3242–3254. [Google Scholar] [CrossRef]
  18. Tao, J.; Liu, J.; Hu, Z.; Zhao, Q.; Chen, G.; Ju, B. Initial Assessment of the BDS-3 PPP-B2b RTS compared with the CNES RTS. GPS Solut. 2021, 25, 1–16. [Google Scholar] [CrossRef]
  19. Xu, Y.; Yang, Y.; Li, J. Performance evaluation of BDS-3 PPP-B2b precise point positioning service. GPS Solut. 2021, 25, 1–14. [Google Scholar] [CrossRef]
  20. Cox, D.B. Integration of GPS with Inertial Navigation Systems. Navig. J. Inst. Navig. 1978, 25, 236–245. [Google Scholar] [CrossRef]
  21. Abd Rabbou, M.; El-Rabbany, A. Tightly coupled integration of GPS precise point positioning and MEMS-based inertial systems. GPS Solut. 2015, 19, 601–609. [Google Scholar] [CrossRef]
  22. Li, X.; Li, X.; Huang, J.; Shen, Z.; Wang, B.; Yuan, Y.; Zhang, K. Improving PPP-RTK in urban environment by tightly coupled integration of GNSS and INS. J. Geod. 2021, 95, 1–18. [Google Scholar] [CrossRef]
  23. Liu, S.; Sun, F.; Zhang, L.; Li, W.; Zhu, X. Tight integration of ambiguity-fixed PPP and INS: Model description and initial results. GPS Solut. 2016, 20, 39–49. [Google Scholar] [CrossRef]
  24. Roesler, G.; Martell, H. Tightly coupled processing of precise point positioning (PPP) and INS data. In Proceedings of the 22nd International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS 2009), Savannah, GA, USA, 22–25 September 2009; pp. 1898–1905. [Google Scholar]
  25. Zhang, X.; Zhu, F.; Zhang, Y.; Mohamed, F.; Zhou, W. The improvement in integer ambiguity resolution with INS aiding for kinematic precise point positioning. J. Geod. 2019, 93, 993–1010. [Google Scholar] [CrossRef]
  26. China Satellite Navigation Office (CSNO). BeiDou Navigation Satellite System Signal in Space Interface Control Document Precise Point Positioning Service Signal PPP-B2b (Version 1.0). Available online: http://www.beidou.gov.cn/xt/gfxz/202008/P020200803362062482940.pdf (accessed on 30 September 2020).
  27. Leick, A.; Rapoport, L.; Tatarnikov, D. GPS Satellite Surveying, 4th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2015. [Google Scholar]
  28. Zhang, B.; Hou, P.; Liu, T. PPP-RTK functional models formulated with undifferenced and un combined GNSS observations. Satell. Navig. 2022, 3, 1–15. [Google Scholar] [CrossRef]
  29. Ashby, N. Relativity in the Global Positioning System Imprint/Terms of Use. Living Rev. Relativ. 2003, 6, 1–42. [Google Scholar] [CrossRef] [Green Version]
  30. Petit, G.; Luzum, B. IERS Conventions (2010); Bureau International Des Poids et Mesures Sevres: Frankfurt, Germany, 2010. [Google Scholar]
  31. Wu, J.-T.; Wu, S.C.; Hajj, G.A.; Bertiger, W.I.; Lichten, S.M. Effects of antenna orientation on GPS carrier phase. In Proceedings of the AAS/AIAA Astrodynamics Conference, Durango, CO, USA, 19–22 August 1991; pp. 1647–1660. [Google Scholar]
  32. Geng, J.; Meng, X.; Dodson, A.; Teferle, F. Integer ambiguity resolution in precise point positioning: Method comparison. J. Geod. 2010, 84, 569–581. [Google Scholar] [CrossRef] [Green Version]
  33. Niu, X.; Zhang, Q.; Gong, L.; Liu, C.; Zhang, H.; Shi, C.; Wang, J.; Coleman, M. Development and evaluation of GNSS/INS data processing software for position and orientation systems. Surv. Rev. 2015, 47, 87–98. [Google Scholar] [CrossRef]
  34. Shin, E.H. Estimation Techniques for Low-Cost Inertial Navigation. Ph.D. Thesis, University of Calgary, Calgary, AB, Canada, 2005. [Google Scholar]
  35. Groves, P.D. Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems, 2nd ed.; Artech House: Norwood, MA, USA, 2013. [Google Scholar]
  36. Yang, Y. Adaptive Navigation and Kinematic Positioning; Surveying and Mapping Press: Beijing, China, 2006. [Google Scholar]
  37. NovAtel Inc. IMU-ISA-100C Product Sheet. Available online: https://hexagondownloads.blob.core.windows.net/public/Novatel/assets/Documents/Papers/IMUISA100CD19539v2/IMU-ISA-100C Product Sheet.pdf (accessed on 15 October 2021).
  38. Analog Devices Inc. ADIS16505 Data Sheet. Available online: https://www.analog.com/media/en/technical-documentation/data-sheets/adis16505.pdf (accessed on 20 October 2021).
  39. NovAtel Inc. PwrPak7 Product Sheet. Available online: https://hexagondownloads.blob.core.windows.net/public/Novatel/assets/Documents/Papers/PwrPak7-Product-Sheet/PwrPak7-Product-Sheet.pdf (accessed on 15 October 2021).
  40. Ge, M.; Gendt, G.; Rothacher, M.; Shi, C.; Liu, J. Resolution of GPS carrier-phase ambiguities in precise point positioning (PPP) with daily observations. J. Geod. 2008, 82, 389–399. [Google Scholar] [CrossRef]
  41. Saastamoinen, J. Atmospheric correction for the troposphere and stratosphere in radio ranging satellites. J. Geophys. Res. Atmosph. 1972, 15, 247–251. [Google Scholar] [CrossRef]
  42. Böhm, J.; Niell, A.; Tregoning, P.; Schuh, H. Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data. Geophys. Res. Lett. 2006, 33, 1–4. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Experimental equipment mounted on the vehicle.
Figure 1. Experimental equipment mounted on the vehicle.
Remotesensing 14 02957 g001
Figure 2. The trajectory of land vehicle in Experiment A.
Figure 2. The trajectory of land vehicle in Experiment A.
Remotesensing 14 02957 g002
Figure 3. Number of visible BDS-3/GPS satellites and PDOPs in Experiment A.
Figure 3. Number of visible BDS-3/GPS satellites and PDOPs in Experiment A.
Remotesensing 14 02957 g003
Figure 4. Positioning errors of PPP-B2b only (top panel), PPP-B2b/MEMS-IMU loosely coupled integration (middle panel), and PPP-B2b/tactical−IMU loosely coupled integration (bottom panel), for Experiment A.
Figure 4. Positioning errors of PPP-B2b only (top panel), PPP-B2b/MEMS-IMU loosely coupled integration (middle panel), and PPP-B2b/tactical−IMU loosely coupled integration (bottom panel), for Experiment A.
Remotesensing 14 02957 g004
Figure 5. Positioning errors of the PPP-B2b/MEMS−IMU loosely coupled integration during the simulated GNSS signal outages.
Figure 5. Positioning errors of the PPP-B2b/MEMS−IMU loosely coupled integration during the simulated GNSS signal outages.
Remotesensing 14 02957 g005
Figure 6. Positioning errors of the PPP-B2b/tactical−IMU loosely coupled integration during the simulated GNSS signal outages.
Figure 6. Positioning errors of the PPP-B2b/tactical−IMU loosely coupled integration during the simulated GNSS signal outages.
Remotesensing 14 02957 g006
Figure 7. The trajectory of the land vehicle in Experiment B.
Figure 7. The trajectory of the land vehicle in Experiment B.
Remotesensing 14 02957 g007
Figure 8. Number of visible BDS-3/GPS satellites and PDOPs in Experiment B.
Figure 8. Number of visible BDS-3/GPS satellites and PDOPs in Experiment B.
Remotesensing 14 02957 g008
Figure 9. Positioning errors of PPP-B2b only (top panel), PPP-B2b/MEMS-IMU loosely coupled integration (middle panel), and PPP-B2b/tactical−IMU loosely coupled integration (bottom panel) during the period from 09:00:00 to 09:34:59.
Figure 9. Positioning errors of PPP-B2b only (top panel), PPP-B2b/MEMS-IMU loosely coupled integration (middle panel), and PPP-B2b/tactical−IMU loosely coupled integration (bottom panel) during the period from 09:00:00 to 09:34:59.
Remotesensing 14 02957 g009
Figure 10. Two bridge scenarios during the period from 09:35:00 to 09:49:59.
Figure 10. Two bridge scenarios during the period from 09:35:00 to 09:49:59.
Remotesensing 14 02957 g010
Figure 11. Positioning errors of PPP-B2b only (top panel), PPP-B2b/MEMS-IMU loosely coupled integration (middle panel), and PPP-B2b/tactical−IMU loosely coupled integration (bottom panel) during the period from 09:35:00 to 09:49:59.
Figure 11. Positioning errors of PPP-B2b only (top panel), PPP-B2b/MEMS-IMU loosely coupled integration (middle panel), and PPP-B2b/tactical−IMU loosely coupled integration (bottom panel) during the period from 09:35:00 to 09:49:59.
Remotesensing 14 02957 g011
Table 1. Parameters of the IMU sensors used in the experiments.
Table 1. Parameters of the IMU sensors used in the experiments.
IMU SensorGradeSampling RateBiasRandom Walk
Gyro.
(°/h)
Acc.
(mGal)
Angular
(°/ h )
Velocity
( m / s / h )
IAS100CTactical200 Hz0.7510000.030.1
ADIS-16505MEMS100 Hz1015000.340.18
Table 2. Strategies of PPP-B2b positioning.
Table 2. Strategies of PPP-B2b positioning.
ItemModel/Strategy
Satellite orbit and clockApplying the PPP-B2b precise products
Weight methodElevation-dependent weight [40]
Phase windupCorrected [31]
Sagnac effectCorrected [29]
Special relativistic effectCorrected [29]
Shapiro time delayCorrected [29]
Tidal effectsCorrected according to IERS Conventions 2010 [30]
TroposphereSaastamoinen model [41] is adopted to correct hydrostatic delay, the zenith wet delay is estimated for each epoch as a random walk noise, and Global Mapping Function [42] is used.
Phase ambiguity Estimated as a constant for each ambiguity arc
EstimatorKalman filter
Table 3. Statistics of the positioning errors (unit: cm).
Table 3. Statistics of the positioning errors (unit: cm).
SchemesBiasRMS
EastNorthUpEastNorthUp
PPP-B2b6.723.726.58.124.227.5
PPP-B2b/MEMS-IMU6.623.326.07.924.027.0
PPP-B2b/tactical-IMU6.223.525.87.123.826.5
Table 4. Statistics of positioning errors for PPP-B2b/INS loosely coupled integration during GNSS outages (unit: cm).
Table 4. Statistics of positioning errors for PPP-B2b/INS loosely coupled integration during GNSS outages (unit: cm).
SchemesDirection5 s10 s15 s20 s25 s30 s
PPP-B2b/tactical-IMUEast0.51.83.97.712.718.6
North0.82.55.39.314.621.8
Up0.41.32.53.75.06.1
PPP-B2b/MEMS-IMUEast3.817.445.797.4180.5300.6
North8.035.287.3174.6308.9498.0
Up1.87.215.724.733.441.0
Table 5. The positioning accuracy during the period from 09:35:00 to 09:49:59 (unit: cm).
Table 5. The positioning accuracy during the period from 09:35:00 to 09:49:59 (unit: cm).
SchemesBiasRMS
EastNorthUpEastNorthUp
PPP-B2b14.948.112.9114.2200.0244.3
PPP-B2b/MEMS-IMU31.725.125.241.044.564.0
PPP-B2b/tactical-IMU32.124.225.141.241.062.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, X.; Nie, Z.; Wang, Z.; Wang, B.; Du, Q. Performance Assessment of BDS-3 PPP-B2b/INS Loosely Coupled Integration. Remote Sens. 2022, 14, 2957. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14132957

AMA Style

Xu X, Nie Z, Wang Z, Wang B, Du Q. Performance Assessment of BDS-3 PPP-B2b/INS Loosely Coupled Integration. Remote Sensing. 2022; 14(13):2957. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14132957

Chicago/Turabian Style

Xu, Xiaofei, Zhixi Nie, Zhenjie Wang, Boyang Wang, and Qinghuai Du. 2022. "Performance Assessment of BDS-3 PPP-B2b/INS Loosely Coupled Integration" Remote Sensing 14, no. 13: 2957. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14132957

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