Next Article in Journal
Development of an Unmanned Aerial Vehicle-Borne Crop-Growth Monitoring System
Next Article in Special Issue
LiDAR-IMU Time Delay Calibration Based on Iterative Closest Point and Iterated Sigma Point Kalman Filter
Previous Article in Journal
Domain 4 (D4) of Perfringolysin O to Visualize Cholesterol in Cellular Membranes—The Update
Previous Article in Special Issue
Utilization of a Terrestrial Laser Scanner for the Calibration of Mobile Mapping Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Unified Model for BDS Wide Area and Local Area Augmentation Positioning Based on Raw Observations

1
National Time Service Center, Chinese Academy of Sciences, Shu Yuan Road, Xi’an 710600, China
2
Key Laboratory of Precision Navigation and Timing Technology, Chinese Academy of Sciences, Shu Yuan Road, Xi’an 710600, China
3
College of earth sciences, University of Chinese Academy of Sciences, Yu Quan Road, Beijing 100049, China
4
Germany Research Center for Geosciences, Telegrafenberg, Potsdam 14473, Germany
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(3), 507; https://doi.org/10.3390/s17030507
Submission received: 7 January 2017 / Revised: 26 February 2017 / Accepted: 1 March 2017 / Published: 3 March 2017
(This article belongs to the Special Issue Multi-Sensor Integration and Fusion)

Abstract

:
In this study, a unified model for BeiDou Navigation Satellite System (BDS) wide area and local area augmentation positioning based on raw observations has been proposed. Applying this model, both the Real-Time Kinematic (RTK) and Precise Point Positioning (PPP) service can be realized by performing different corrections at the user end. This algorithm was assessed and validated with the BDS data collected at four regional stations from Day of Year (DOY) 080 to 083 of 2016. When the users are located within the local reference network, the fast and high precision RTK service can be achieved using the regional observation corrections, revealing a convergence time of about several seconds and a precision of about 2–3 cm. For the users out of the regional reference network, the global broadcast State-Space Represented (SSR) corrections can be utilized to realize the global PPP service which shows a convergence time of about 25 min for achieving an accuracy of 10 cm. With this unified model, it can not only integrate the Network RTK (NRTK) and PPP into a seamless positioning service, but also recover the ionosphere Vertical Total Electronic Content (VTEC) and Differential Code Bias (DCB) values that are useful for the ionosphere monitoring and modeling.

1. Introduction

Relative positioning and Precise Point Positioning (PPP) are two basic approaches in precise Global Navigation Satellite System (GNSS) data processing. For the relative positioning, data from a network of stations are analyzed simultaneously to estimate the receiver positions and other parameters such as zenith troposphere delays and the integer-cycle phase ambiguities [1,2]. It shows an advantage in effectively eliminating the common errors derived from the multiple satellites and stations, including the GNSS receiver clock and satellite clock errors. Thus, the integer-cycle phase ambiguities can be easily resolved to their correct integer values to achieve mm-cm level precision in a local reference frame [3]. However, it requires the simultaneous observations for the whole network and the selection of reference stations.
In terms of the PPP technology [4], it relies on fixed satellite orbits and satellite clocks to estimate the positions of receivers, phase ambiguities, zenith troposphere delays, and receiver clock parameters. As no specific reference station is needed, PPP is more flexible and has been widely used in many high-precision areas [5,6,7]. After an initialization of half an hour or longer, the dynamic PPP can obtain a centimeter precision in the global reference. Furthermore, the convergence time and positioning accuracy can be improved by the newly developed PPP ambiguity resolution methods [8,9,10,11]. In addition, the latest development on speeding up PPP initialization can be resolving un-differenced ambiguities for multi-GNSS simultaneously [12,13,14,15].
The Real-Time Kinematic (RTK) positioning is also of great interest for the GNSS community. Generally, the ambiguity fixed solutions can be achieved using a few seconds of measurements. The positioning accuracy is at the centimeter level when the baselines are shorter than a few tens of kilometers [16,17,18]. In the case of longer baselines, a network of Continuously Operating Reference Stations (CORS) can be established to spatially interpolate the atmospheric corrections, and thus the positioning accuracy is significantly improved compared to the RTK [19,20,21]. This service mode is named as Network RTK (NRTK). The deficiencies of NRTK are also obvious. On one hand, the interpolated atmospheric corrections are only effective within the area of the CORS network, as the accuracy of atmospheric corrections is rapidly degraded for users located out of the CORS network [22,23]. On the other hand, due to the high cost of deploying the CORS networks, NRTK can only cover a regional area and hence can hardly evolve into a precise positioning service on the global scale [24]. Accordingly, some new positioning algorithms were developed. Ge et al. proposed the Un-differenced network RTK (URTK), which utilized Double Differenced (DD) observations for the server solutions and Un-differenced (UD) observations for the user solutions [25]. It provided the rovers a rapid and precise absolute positioning service. Teunissen et al. proposed a CORS based PPP-RTK, which made use of UD observations for both the server and rover solutions [26]. URTK is different from the CORS based PPP-RTK, as URTK does not need clock corrections and ionospheric delays for the reference network solution, and URTK’s rover can use precise or broadcast ephemeris products. It is also different from regional reference network-augmented precise point positioning as the DD-ambiguity resolution is more rapid than fixing the integer ambiguity in PPP for the reference network [27], thus shortening the initialization time of modeling the corrections. The main differential corrections between NRTK and URTK are also different. As the NRTK usually uses the Virtual Reference Station (VRS) technique, it will form a virtual observation and broadcast it to the corresponding user, then the user using the standard RTK for the fast positioning [21]. And for the URTK, it usually broadcasts the UD Observed Minus Computed (OMC) corrections from the reference stations. The users then selected and interpolated these corrections and using PPP for the fast positioning [25]. Meanwhile, the URTK is a one-way communication, with the need of receiving the server’s corrections. In addition, it is possible to integrate URTK and PPP into a seamless positioning service [28,29].
Nowadays, some researchers have proposed a new PPP model using the raw observations [11,30,31,32,33], which shows several advantages compared to the ionospheric-free observations [34]. Firstly, it will reduce the observing noise and increase the number of the observations. Secondly, the ionosphere and Differential Code Bias (DCB) information can be recovered. Thirdly, it can be used for the single-frequency users. For example, Gu et al. discussed the BeiDou Navigation System (BDS) PPP with raw observations and tried to resolve the un-differenced ambiguities [35]. However, as most of the previous studies mainly focused on (Global Positioning System) GPS, GLObalnaya NAvigatsionnaya Sputnikovaya Sistema (GLONASS) and BDS PPP, further investigations concerning the BDS NRTK using raw observations are still needed.
So far, 14 BDS satellites (five Geostationary Earth Orbit (GEO), five Inclined Geosynchronous Orbit (IGSO), four Mid Earth Orbit (MEO)) are in orbit and providing positioning service for the Asian-Pacific region [33,36,37,38]. Besides, the BDS ground-based augmentation system is also under construction. It consists of 150 base stations for providing the NRTK and/or the PPP service in China. Therefore, the study of the unified model for BDS wide area and local area augmentation positioning based on raw observations is of crucial significance, not only for the integration of PPP and NRTK but also for geosciences application.
In this contribution, we proposed a unified model for BDS wide area and local area augmentation positioning based on raw observations. First of all, the separated PPP and/or URTK model with raw observations are introduced. Secondly, the integrated model of PPP and/or URTK with the raw observations is presented. Then, the experimental BDS data are applied for the purpose of evaluation and validation. Finally, the conclusions and discussion are provided.

2. The PPP Solution Model Based on Raw Observations

Generally, the ionospheric-free combinations with dual-frequency observations are applied to eliminate the first-order ionospheric delays, which, however, will lead to the increase of the measurement noise. In this study, as the raw observations are used, the ionospheric delays are estimated as unknown parameters by adding constrains of a priori information of their spatial and temporal characteristics at each epoch.

2.1. A priori Constraint

Usually, the ionosphere models such as the International Reference Ionosphere (IRI) model [39], Bent model [40] and Klobuchar model [41] can be used for the priori constraint. In this study, the Global Ionospheric Model (GIM) provided by the Center of Orbit Determination in Europe (CODE) is used. The ionospheric delays provided by the priori ionospheric model can be expressed as:
i o n i , s p r i o r p = M I i , s p × [ 40.28 / f s 2 × V T E C i , s p r i o r p ]
where, i o n i , s p r i o r p is the priori constrains for ionospheric delays, M I i , s p is the mapping function of the ionospheric delays. V T E C i , s p r i o r p represents the total electron content in the vertical component by the prior ionospheric model, p denotes the satellite, i denotes the receiver, s is the signal channel, and f is the signal frequency.

2.2. Spatial Constraint

The ionosphere VTEC’s spatial behavior can be expressed as a polynomial function of the latitude differences and longitude differences in the sun-fixed reference frame [11,42]:
i o n i , s s p a c e p = M I i , s p × [ 40.28 f s 2 × ( n 0 = 0 n m 0 = 0 m E n 0 m 0 ( φ φ 0 ) n 0 ( λ λ 0 ) m 0 ) ]
where, i o n i , s s p a c e p represents the value of ionospheric space constraint, m , n are the order numbers of the polynomial, the value of the order number can be set as two. φ , λ are latitude and longitude of IPP in the sun-fixed reference frame; and φ 0 , λ 0 describe latitude and longitude of the station in the sun-fixed reference frame; E n 0 m 0 is the estimated polynomial coefficient. The coefficients can be estimated epoch wise by Equation (3) when the number of visible satellites is more than four for the BDS system. If the satellite in view is less than four, the coefficients of the last epoch are used [42].
n 0 = 0 n m 0 = 0 m E n 0 m 0 ( φ φ 0 ) n 0 ( λ λ 0 ) m 0 ) = V T E C i , p r i o r p

2.3. Temporal Constraint

As the temporal variation of ionospheric VTEC is very slow, the epoch-wise variation of ionosphere delays can be modeled by a stochastic process. In this study, the temporal correlations of the DD ionospheric delay residuals are considered as a random walk process [43]. The temporal constraint of each satellite-receiver pair is described as follows:
i o n i , s t e m p p = I i , s p ( l a s t ) + Δ δ I i , s p
where, i o n i , s t e m p p is the value of ionospheric temporal constraint, Δ δ I i , s p represents the epoch-wise variation of the ionospheric delay.

2.4. The PPP Model

In our PPP model, the ionospheric delays are estimated as unknown parameters by adding constraint of a priori information, i.e., the spatial and temporal characteristics. Concerning the tropospheric delays, the Saastamoinen (SAAS) troposphere model and Global Mapping Function (GMF) are used to calculate the tropospheric error corrections, and the tropospheric residuals are estimated as random walk processes [44]. The satellite DCB is corrected by DCB files provided by the International GNSS Service (IGS), and the receiver DCB is estimated as an unknown parameter at each epoch. The satellite ephemeris errors are corrected by real-time corrections provided by the International GNSS Monitoring and Assessment Service (IGMAS) analysis center, and the receiver clock errors are estimated as white noises. The antenna parameters, Earth rotation parameters and DCB parameters are obtained from the IGS. The antenna phase center offsets (PCO) and variations (PCV) for satellites and receivers are corrected by the absolute antenna model [45]. In addition, the effects including solid tide, ocean tide, relativity effect, and phase wind up are also corrected following the descriptions in Dach et al. [44]. The PPP model can be written as follows:
P i , s p = ρ p + [ t r o p + i o n i , s p D C B i , s c δ t i ] + [ c δ t p + D C B i , s p + M P p ] + ε P , i , s p
ε P , i , s p ~ N ( 0 , σ P , i , s 2 )
Φ i , s p = ρ p + [ t r o p i o n i , s p c δ t i + λ s A m b s ] + [ c δ t p + M Φ p ] + ε Φ , i , s p
ε Φ , i , s p ~ N ( 0 , σ Φ , i , s 2 )
i o n p r i o r p = i o n i , s p + ε p r i o r
ε p r i o r ~ N ( 0 , σ p r i o r 2 )
i o n t e m p p = i o n i , s p + ε t e m p
ε t e m p ~ N ( 0 , σ t e m p 2 )
i o n s p a c e p = i o n i , s p + ε s p a c e
ε s p a c e ~ N ( 0 , σ s p a c e 2 )
where, P , Φ are the pseudo-range and carrier phase measurements, respectively. ρ is the geometric distance from satellite to the receiver. p represents the satellite, i represents the receiver, s denotes the frequency, λ and A m b are the wave length and phase ambiguity, respectively. M expresses the combined correction of other un-modeled errors. i o n p r i o r , i o n t e m p , i o n s p a c e are the virtual observations of the ionospheric priori information constraint, i.e., temporal constraint and spatial constraint. ε is the observation noise, δ is the standard deviation. Based on Equations (5)–(14), the least square is employed for data processing. When N satellites are tracked for one station at each epoch, there are 3 × N + 6 parameters needed to be estimated. These parameters are summarized in Table 1.
In parameters estimation, the receiver DCB, troposphere residuals and ionosphere parameters are estimated as random walk processes at each epoch. The phase ambiguity is treated as a constant in a continuous arc, which should be re-initialization when a cycle slip occurs. The measurement noise level and the satellite elevations are used to determine the observation weights. And the weights of the ionospheric virtual observations are determined by local time, noise level and satellite elevations [32,42].

3. The URTK Solution Model Based on Raw Observations

For the URTK with raw observations, the key issues include baseline solution of the reference stations, retrieval and interpolation of the un-differenced observation corrections, the rover’s positioning and ambiguity resolution.

3.1. The Baseline Solution Model for the Reference Stations with Raw Observations

For the baseline solution of the reference station, we used the DD observations by estimating the baseline components ( d X j , d Y j , d Z j ), DD ambiguity ( A m b L s , i j p q ) and DD atmospheric delay ( δ I i j p q , δ T i j p q ) parameters. The linearized DD pseudo-range and carrier phase observation equations are expressed as follows:
P P s , i j p q = ρ i j p q + ( l j q l j p ) d X j + ( m j q m j p ) d Y j + ( n j q n j p ) d Z j + f 1 2 f s 2 δ I i j p q + M T δ T i j p q + ε P s i j p q
ε P s i j p q ~ N ( 0 , σ P s i j p q 2 )
Φ L s , i j p q = ρ i j p q + ( l j q l j p ) d X j + ( m j q m j p ) d Y j + ( n j q n j p ) d Z j f 1 2 f s 2 δ I i j p q + M T δ T i j p q + λ s A m b L s , i j p q + ε L s i j p q
ε L s i j p q ~ N ( 0 , σ L s i j p q 2 )
where the subscripts i and j represent the base and user receivers, respectively, and the superscripts p and q represent the pair of satellites while p is the reference one. s represents the signal channel, “1” is the L1 channel. The combined subscript i j means differencing between two receivers, and the combined superscript p q denotes the differencing between two satellites. P , Φ are the pseudo-range and carrier phase observables, respectively, and λ denotes the wavelength of the carrier. The baseline components d X , d Y , d Z are to be solved, and the symbols l , m and n are the unit vectors on the line-of sight from receiver to satellite. A m b is the integer ambiguity of the carrier phase, and δ I and δ T are the ionospheric delay and tropospheric delay, respectively. f is the carrier phase frequency, ε denotes the measurement noise, σ expresses the standard deviation of the noise, and ρ is the geometric distance from satellite to the receiver. Moreover, the mapping function M T for the wet components of tropospheric delay is defined as follows:
M T = [ M T i p + M T i q + M T j q + M T j q ] / 4.0
In Equations (15)–(19), the least square is used for data processing. When N common satellites are tracked for one pair station at each epoch, there will be 3 × ( N − 1) + 4 parameters needed to be estimated. These parameters are listed in Table 2.
For the parameters estimation, the station coordinates are tightly constrained. The troposphere and ionosphere parameters are estimated as random walk processes. The DD phase ambiguity is treated as a constant in a continuous period, which should be re-initialization when a cycle slip occurs. And the observation weights are determined by the measurement noise level and satellite elevation [32,42].

3.2. The Retrieval of the Un-Differenced Observation Corrections

In order to retrieve the UD observation corrections, the DD ambiguities need to be fixed firstly. The LAMBDA method is utilized here for the integer ambiguity searching [46]. The retrieval of the UD observation corrections can be divided into two steps. The first step is to recover the UD ambiguity, and the second part is to retrieve the UD observation corrections.
Based on the constraint that the DD ambiguity is integer and has been fixed, the other un-differenced ambiguities which formed the DD ambiguity can be recovered using Equation (20) when the integer part of some un-differenced ambiguities are defined as the datum in one pair of baseline:
[ A m b i p A m b i q A m b j p A m b j q ] = [ 1 0 0 1 0 1 0 1 0 0 1 1 0 0 0 1 ] 1 [ b i p b i q b j p A m b i j p q ]
Here, b i p , b i q , b j p are the selected UD integer ambiguity datum. A m b i p , A m b i q , A m b j p , A m b j q are the recovered UD integer ambiguities. A m b i j p q is the fixed DD integer ambiguity.
By transferring the datum ambiguities of the reference satellites through different satellite pairs, the other un-differenced ambiguities of this baseline can be recovered. Then, by the transfer of datum ambiguities of the reference stations through different baselines, all the un-differenced ambiguities of the whole network can be recovered. It is noteworthy that the selected UD ambiguities datum must be consistent for the same network.
After all UD ambiguities are recovered, the UD observation corrections can be calculated using the following equations:
o m c P i , s p = P i , s p ρ i , s p c δ t i , s + c δ t p + M o d e l P i , s p
o m c Φ i , s p = λ Φ i , s p ρ i , s p + λ A m b i , s p c δ t i , s + c δ t p + M o d e l Φ i , s p
Here, M o d e l describes the modeled errors, such as the tropospheric errors, relativity errors, phase center offset error and phase wind up. o m c is the comprehensive correction, it contains the orbit residual, clock residual, Uncalibrated Phase Delay (UPD), DCB and some other un-modeled errors.

3.3. The Interpolation of the UD Observation Corrections

At the user end, three nearby reference stations are selected as the augmentation stations. Based on these UD observation corrections from the augmentation, we can calculate the rover’s corrections using Equation (23):
o m c s , u p = a ( o m c s , 2 p o m c s , 1 p ) + b ( o m c s , 3 p o m c s , 1 p ) + o m c s , 1 p
a = ( x 3 x 1 ) ( y u y 1 ) ( x u x 1 ) ( y 3 y 1 ) ( x 3 x 1 ) ( y 2 y 1 ) ( x 2 x 1 ) ( y 3 y 1 )
b = ( x u x 1 ) ( y 2 y 1 ) ( x 2 x 1 ) ( y u y 1 ) ( x 3 x 1 ) ( y 2 y 1 ) ( x 2 x 1 ) ( y 3 y 1 )
Here, 1 , 2 , 3 represents the augmentation stations, u refers to the rover station, x , y are the plane coordinates of each station.

3.4. The URTK Model

At the user end, the observation corrections can be interpolated using Equation (23) at first. Then the UD observations are corrected and the positioning is performed. The positioning model at the user end can be written as:
P i , s p = ρ i p + l i p X i + m i p Y i + n i p Z i + [ D C B F c δ t i ] + [ o m c P i , s p ] + M o d e l P i , s p + ε P i , s p
ε P i , s p ~ N ( 0 , σ P s 2 )
Φ i , s p = ρ i p + l i p X i + m i p Y i + n i p Z i + [ λ s A m b s p c δ t i ] + [ o m c Φ i , s p ] + M o d e l Φ i , s i p + ε Φ i , s p
ε Φ i , s p ~ N ( 0 , σ Φ s 2 )
Using Equations (26) and (29), the least square is employed for the estimation. The estimated unknown parameters are the positioning components ( X i , Y i , Z i ), ambiguity ( A m b s , i p ) and the combined receiver clock error δ t i . The clock error is estimated at each epoch as a white noise. The observation corrections include the DCB and UPD of both satellites and receivers at the reference stations. For the user positioning, the satellite UPD can be cancelled when using the observation corrections, but the receiver DCB still need to be estimated. Meanwhile, the receiver UPD cannot be eliminated, thus the ambiguity is not integer. Fortunately, the Single-Differenced (SD) projection can be applied for the ambiguity resolution to realize the fast convergence and high-precision positioning.

4. The Integration of PPP and URTK Based on Raw Observations

4.1. An Integrated Model

In order to integrate PPP and URTK based on raw observations into a seamless positioning service, the key point is to develop a unified model of PPP and URTK. The unified model of PPP and URTK can be written as follows:
P i , s p = ρ p + [ E s t i m a t i o n ] + [ C o r r e c t i o n ] + ε P , i , s p
ε P , i , s p ~ N ( 0 , σ P , i , s 2 )
Φ i , s p = ρ p + [ E s t i m a t i o n ] + [ C o r r e c t i o n ] + ε Φ , i , s p
ε Φ , i , s p ~ N ( 0 , σ Φ , i , s 2 )
V i r t u a l = [ E s t i m a t i o n ] + ε v , i , s p
ε v , i , s p ~ N ( 0 , σ v , i , s 2 )
Here, E s t i m a t i o n consists of the unknown parameters such as the coordinates of the receiver, phase ambiguity A m b , troposphere delay t r o p , ionosphere delay i o n s p , hardware delay D C B s and the receiver clock error δ t i . C o r r e c t i o n represents the corresponding correction parts in Equations (5)–(8), (21) and (22). Moreover, in terms of the integrated model for the URTK solution, the ionosphere estimates are the residuals. Thus, the values of ionosphere virtual observations can be set to zero for a regional area, i.e., V i r t u a l = 0 .

4.2. Data Processing Flow

The data processing is shown in Figure 1, which consists of two parts: the server and the user end. At the server end, the global reference stations are processed in real-time to generate the precise orbit, clock and ionosphere products, while the regional reference stations are processed in real-time to generate the OMC corrections. Both the global and regional products are broadcasted to the users by the caster.
At the user end, when they are global users, the satellite orbit and clock errors are corrected by the received corrections, and the other modeling errors are corrected by using the empirical models. Moreover, the un-differenced and un-combined solution models for the absolute positioning are used to get the PPP solutions. For the regional users, all the observation errors are considered as the comprehensive errors, which are corrected in the observation domain. Besides, the un-differenced and un-combined solution models for the absolute positioning are applied to achieve the URTK solutions.

5. Validations and Results

In order to assess the proposed algorithm for BDS wide area and local area augmentation positioning, a series of experiments were carried out in Xi’an. The data processing was carried out using the GNSS software package developed by our research group, which is capable of providing both network solution and PPP/RTK services [32]. Although its old version is only applicable to the GPS measurements, this software has been recently developed and is now able to handle with the BDS data. In the following, the results based on these experiments will be presented.

5.1. Data Description

The BDS observations are collected on day of year (DOY) 080 to 083 in 2016.The antenna type is HG-GOYH7151 and the receiver type is UR380, which is capable of tracking the BDS B1, B2 and B3 signals. The sampling rate is 1 s, observing from UTC 0:0:00 to 23:59:59. The data are processed in a simulated real-time mode, which indicates that all the raw observations were input to the software package epoch by epoch similar to the real-time stream, and the station positions are estimated and output epoch by epoch. The observation signals of B1 and B2 are used for the dual-frequency solution, while those of the B1 are used for the single-frequency solution. The broadcast ephemeris products are utilized here. The stations are located in the city of Xi’an in Shan Xi Province. Figure 2 shows the distribution of the collected stations, where three stations marked in red are selected as the reference stations (i.e., NT02, NT03, NT04) and the station marked in black is selected as the rover station (i.e., NT01). The distances for the three baseline pairs (NT04-NT03, NT04-NT02 and NT03-NT02) are 30, 33 and 27 km, respectively. The precise satellite orbit, clock, ionosphere and DCB products are provided by the IGMAS and used in a simulated real-time mode. The antenna PCO and PCV parameters, as well as the Earth Rotation Parameters (ERP) files are provided by the IGS.

5.2. Accuracy Validation

The collocated data were processed using the unified model with four different strategies in simulated real-time modes: (a) Using the dual-frequency (B1 and B2) observations and real-time orbit, clock products, processed in a static mode (It is equal to static PPP solution); (b) Using the dual-frequency (B1 and B2) observations and real-time orbit, clock products, processed in a dynamic mode (It is equal to dynamic PPP solution); (c) Using the dual-frequency (B1 and B2) observations and real-time UD OMC corrections, processed in a static mode (It is equal to static URTK solution); (d) Using the dual-frequency (B1 and B2) observations and real-time UD OMC corrections, processed in a dynamic mode (It is equal to dynamic URTK solution). The reference coordinates of these stations are calculated by using the static GNSS-PPP of three-day observations, the coordinate accuracy is better than 5 mm in horizontal and 1 cm in vertical [42].
Figure 3 shows the comparison of positioning errors for different solution strategies. We can notice that the results of dynamic model show larger noise than the static model, especially in the vertical component. In Table 3, the statistics of the positioning accuracy for four-day solutions is presented. For the PPP, the Root Mean Square (RMS) is (8, 6, 15) mm and (25, 16, 56) mm in East, North and Up components for the static and dynamic solutions, respectively. For the URTK, the RMS is (6, 5, 14) mm and (21, 17, 30) mm in East, North and Up components for the static and dynamic solutions, respectively.

5.3. Convergence Time

In order to make a test on the convergence time, the estimation of ambiguity parameter is reinitialized every half an hour for the PPP solution and every five seconds for the URTK solution. The criterion for the convergence is defined as that the bias of each horizontal component is smaller than 10 cm for PPP and smaller than 3 cm for URTK.
Figure 4 shows the solution convergence of the test. It can be seen that the convergence speed of URTK model is much faster than that of the PPP model. It may be due to the fact that the common errors are precisely corrected by the observation corrections and the satellite to satellite single-differenced ambiguity can be easily fixed by URTK. In addition, Table 4 summarizes the average convergence time for the four days. The mean convergence time is about 1500 and 1.7 s for the PPP and URTK model, respectively.

5.4. Observation Residuals

Figure 5 shows the comparison of observation residuals. For the URTK model, the residuals are at the level of several mm. But for the PPP model, the residuals are a little dispersive. The Standard Deviation (STD) is about 1–2 cm, which means that the errors correction in the URTK model is much better than in the PPP model. Thus the position errors for the PPP model are also much larger than for the URTK mode.

5.5. Ionospheric Delay

In this study, the ionospheric delay is treated as an unknown parameter and estimated together with other positioning parameters. It is also of benefit for retrieving the ionospheric VTEC. For each station, the estimated ionospheric delay can be used to establish a polynomial ionospheric model like Equation (2), thus the VTEC at the station can be obtained. As shown in Figure 6a, the VTEC gained by the PPP solution agree well with the Grid Ionosphere Model (GIM), showing the same variation trend. Figure 6b shows the differences of the VTEC between the grid ionosphere model and the PPP solution. The largest difference of the two solutions is smaller than 2 TECU. Compared to the traditional approach which usually uses the pseudo-range or carrier phase smoothed pseudo-range, our proposed method directly makes use of the carrier phase observations based on the easily operated PPP technology. It provides an alternative way for the ionosphere monitoring and modeling.

5.6. Analysis of Receiver DCB

For precise positioning based on raw observations, both satellite and receiver DCB must be considered. In this study, the satellite DCB are corrected using the published DCB files, and the receiver DCB are estimated as unknown parameters at each epoch.
Figure 7 shows time series of the estimated receiver DCB. It can be seen that the daily variation of the receiver DCB is very small, with a mean STD of 0.15 ns and 0.18 ns for P1C1 and P1P2, respectively. The DCB corrections for the P1C1 are very small, but are much larger for the P1P2. As its magnitude can reach several nanosecond or even larger, the DCB correction plays a significant role for the single-frequency PPP initialization that utilizes the raw observations, and it also critical for the ionosphere estimation.

6. Discussion and Conclusions

In this contribution, we have proposed a unified model for BDS wide area and local area augmentation positioning based on raw observations. In this model, the users can realize both PPP and RTK service. When the users are within the local reference network, the regional observation corrections can be utilized in order to achieve the fast and high precision RTK service, with a convergence time of few seconds and a precision of about 2–3 cm. When the users are out of the local reference network, the global broadcast SSR corrections can be applied to realize the PPP service, achieving a convergence time of 25 min and a positioning accuracy of 10 cm.
The proposed method was realized and validated with the BDS data collocated at four regional stations from DOY 080 to 083 in 2016. In this unified model, the PPP and NRTK can be integrated into a seamless positioning service, which is able to provide an accuracy of about 10 cm anywhere and upgrade to a few centimeters within several seconds with the support of a regional network. Meanwhile, the ionosphere VTEC and DCB values that are important for the ionosphere monitoring and modeling can be recovered.
Up to now, the multi-system and multi-frequency satellite navigation systems are under rapid development. In future, we may focus on studying the combined PPP/RTK algorithm utilizing BDS, GPS, GLONASS and Galileo systems, so as to improve the accuracy, continuity and reliability of the precise positioning. Meanwhile, applying the estimated ionosphere information to establish a high-precision and high-resolution ionosphere model is also one of the focuses for our future work. The influences of the solar cycle, seasonal and local-time dependent ionospheric variability are also needed to be studied later.
Furthermore, this study is focusing on real-time positioning, not only for PPP but also for URTK. Some previous studies have proposed the “Virtual RINEX”, which is basically NRTK for post-processing applications [47,48]. Virtual RINEX data is of potential benefits to both everyday and upmarket boutique users. For the everyday users, it can provide a vital backup that enables workable results in those annoying communications black spot areas. Besides, it also provides the potential of forward and reverses post-processing to minimize the effects of GNSS data gaps and maximize the performance. Thus, the precision, availability and reliability achieved with the Virtual RINEX will be much better than the rea-time users. For the further study, it is also meaningful to use our approach to generate the Virtual RINEX for the post-processing applications.

Acknowledgments

The work was partly supported by the program of National Natural Science Foundation of China (Grant Nos.: 41504006, 41674034), National Key Research and Development Plan of China (Grant No.: 2016YFB0501804), and Chinese Academy of Sciences (CAS) programs of “Pioneer Hundred Talents”, and “The Frontier Science Research Project” (Grant No.: QYZDB-SSW-DQC028).

Author Contributions

Rui Tu and Jinhai Liu conceived and designed the experiments; Rui Zhang and Jinhai Liu performed the experiments; Rui Tu and Pengfei Zhang analyzed the data; Rui Tu and Cuixian Lu contributed reagents/materials/analysis tools; Rui Tu, Cuixian Lu and Xiaochun Lu wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dong, D.N.; Bock, Y. Global Positioning System network analysis with phase ambiguity resolution applied to crustal deformation studies in California. J. Geophys. Res. 1989, 94, 3949–3966. [Google Scholar] [CrossRef]
  2. Blewitt, G. Carrier phase ambiguity resolution for the Global Positioning System applied to geodetic baselines up to 2000 km. J. Geophys. Res. 1989, 94, 10187–10203. [Google Scholar] [CrossRef]
  3. Bock, Y.; Melgar, D.; Crowell, B.W. Real-time strong-motion broadband displacements from collocated GPS and accelerometers. Bull. Seism. Soc. Am. 2011, 101, 2904–2925. [Google Scholar] [CrossRef]
  4. Zumberge, J.F.; Heflin, M.B.; Jefferson, D.C.; Watkins, M. Precise Point Positioning for the efficient and robust analysis of GPS data from large networks. J. Geophys. Res. 1997, 102, 5005–5017. [Google Scholar] [CrossRef]
  5. Bar-Sever, Y.; Blewitt, G.; Gross, R.S.; Hammond, W.C.; Hsu, V.; Hudnut, K.W.; Khachikyan, R.; Kreemer, C.W.; Meyer, R.; Plag, H.P.; et al. A GPS Real Time Earthquake and Tsunami (GREAT) alert system. In Proceedings of the AGU Fall Meeting Abstracts, San Francisco, CA, USA, 14–18 December 2009.
  6. Wright, T.J.; Houli´e, N.; Hildyard, M.; Iwabuchi, T. Real-time, reliable magnitudes for large earthquakes from 1 Hz GPS Precise Point Positioning: The 2011 Tohoku-Oki (Japan) earthquake. Geophys. Res. Lett. 2012, 39, L12302. [Google Scholar] [CrossRef]
  7. Li, X.; Dick, G.; Ge, M.; Heise, S.; Wickert, J.; Bender, M. Real-time GPS sensing of atmospheric water vapor: Precise point positioning with orbit, clock, and phase delay corrections. Geophys. Res. Lett. 2014, 41, 3615–3621. [Google Scholar] [CrossRef]
  8. 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]
  9. Feng, Y. GNSS three carrier ambiguity resolution using ionosphere-reduced virtual signals. J. Geod. 2008, 82, 847–862. [Google Scholar] [CrossRef]
  10. Geng, J.; Meng, X.; Dodson, A.; Ge, M.; Teferle, F. Rapid reconvergences to ambiguity-fixed solutions in precise point positioning. J. Geod. 2010, 84, 705–714. [Google Scholar] [CrossRef]
  11. Li, X.; Ge, M.; Zhang, X.H.; Wickert, J. A method for improving uncalibrated phase delay estimation and ambiguity-fixing in real-time precise point positioning. J. Geod. 2013, 87, 405–416. [Google Scholar] [CrossRef]
  12. Geng, J.; Teferle, F.N.; Meng, X.; Dodson, A.H. Towards PPP-RTK: Ambiguity resolution in real-time precise point positioning. Adv. Space Res. 2011, 47, 1664–1673. [Google Scholar] [CrossRef]
  13. Geng, J.; Shi, C.; Ge, M.; Dodson, A.H.; Lou, Y.; Zhao, Q.; Liu, J. Improving the estimation of fractional-cycle biases for ambiguity resolution in precise point positioning. J. Geod. 2012, 86, 579–589. [Google Scholar] [CrossRef]
  14. Geng, J.; Shi, C. Rapid initialization of real-time PPP by resolving undifferenced GPS and GLONASS ambiguities simultaneously. J. Geod. 2016. [Google Scholar] [CrossRef]
  15. Liu, Y.; Ye, S.; Song, W.; Lou, Y.; Chen, D. Integrating GPS and BDS to shorten the initialization time for ambiguity-fixed PPP. GPS Solut. 2016. [Google Scholar] [CrossRef]
  16. Takasu, T.; Yasuda, A. Kalman filter based integer ambiguity resolution strategy for long-baseline RTK with ionosphere and troposphere estimation. In Proceedings of the 23rd International Technical Meeting of the Satellite Division of the Institute of Navigation, Portland, OR, USA, 21–24 September 2010; pp. 161–171.
  17. Yanase, T.; Tanaka, H.; Ohashi, M.; Kubo, Y.; Sugimoto, S. Long baseline relative positioning with estimating ionosphere and troposphere gradients. In Proceedings of the 23rd International Technical Meeting of the Satellite Division of the Institute of Navigation, Portland, OR, USA, 21–24 September 2010; pp. 196–206.
  18. Dai, L.; Eslinger, D.; Sharpe, T. Innovative algorithms to improve long range RTK reliability and availability. In Proceedings of the ION National Technical Meeting, San Diego, CA, USA, 28–30 January 2007; pp. 860–872.
  19. Fotopoulos, G.; Cannon, M.E. An overview of multi-reference station methods for cm-level positioning. GPS Solut. 2001, 4, 1–10. [Google Scholar] [CrossRef]
  20. Musa, T.A.; Lim, S.; Rizos, C. GPS network-based approach to mitigate residual tropospheric delay in low latitude area. In Proceedings of the ION GNSS 18th International Technical Meeting of the Satellite Division, Long Beach, CA, USA, 13–16 September 2005; pp. 2606–2613.
  21. Snay, R.A.; Soler, T.M. Continuously Operating Reference Station (CORS): History, Applications, and Future Enhancements. J. Surv. Eng. 2008, 134, 95–104. [Google Scholar] [CrossRef]
  22. Grejner-Brzezinska, D.A.; Kashani, I.; Wielgosz, P. On accuracy and reliability of instantaneous network RTK as a function of network geometry, station separation, and data processing strategy. GPS Solut. 2005, 9, 212–225. [Google Scholar] [CrossRef]
  23. Grejner-Brzezinska, D.A.; Arslan, N.; Wielgosz, P.; Hong, C.K. Network calibration for unfavorable reference-rover geometry in network-based RTK: Ohio CORS case study. J. Surv. Eng. ASCE 2009, 135, 90–100. [Google Scholar] [CrossRef]
  24. Rizos, C. Alternatives to current GPS-RTK services and some implications for CORS infrastructure and operations. GPS Solut. 2007, 11, 151–158. [Google Scholar] [CrossRef]
  25. Ge, M.; Zou, X.; Dick, G.; Jiang, W.P.; Wickert, J.; Liu, J.N. An alternative network RTK approach based on undifferenced observation corrections. In Proceedings of the ION GNSS 2010, Institute of Navigation, Portland, OR, USA, 21–24 September 2010; pp. 21–24.
  26. Teunissen, P.J.G.; Odijk, D.; Zhang, B. PPP-RTK: Results of CORS network-based PPP with integer ambiguity resolution. J. Aeronaut. Astronaut. Aviat. Ser. A 2010, 42, 223–230. [Google Scholar]
  27. Li, X.; Zhang, X.H.; Ge, M.R. Regional reference network augmented precise point positioning for instantaneous ambiguity resolution. J. Geod. 2011, 85, 151–158. [Google Scholar] [CrossRef]
  28. Zou, X.; Ge, M.; Tang, W.; Shi, C.; Liu, J. URTK: Undifferenced network RTK positioning. GPS Solut. 2013, 17, 283–293. [Google Scholar] [CrossRef]
  29. Zou, X.; Tang, W.; Shi, C.; Liu, J.N. Instantaneous ambiguity resolution for URTK and its seamless transition with PPP-AR. GPS Solut. 2015, 19, 559–567. [Google Scholar] [CrossRef]
  30. Schonemann, E.; Becker, M.; Springer, T. A new approach for GNSS analysis in a multi-GNSS and multi-signal environment. J. Geod. 2011, 1, 204–214. [Google Scholar] [CrossRef]
  31. Zhang, B.; Ou, J.; Yuan, Y.; Li, Z. Calibration of slant total electron content and satellite-receiver’s differential code biases with uncombined precise point positioning (PPP) technique. Acta Geodaeticaet Cartographica Sinica 2011, 40, 447–453. [Google Scholar]
  32. Tu, R.; Ge, M.; Zhang, H.; Huang, G. The realization and convergence analysis of combined PPP based on raw observation. Adv. Space Res. 2013, 52, 211–221. [Google Scholar] [CrossRef]
  33. Lou, Y.; Zheng, F.; Gu, S.; Wang, C.; Guo, H.; Feng, Y. Multi-GNSS precise point positioning with raw single-frequency and dual-frequency measurement models. GPS Solut. 2016, 20, 849–862. [Google Scholar] [CrossRef]
  34. Li, X.; Ge, M.; Douša, J.; Wickert, J. Real-time precise point positioning regional augmentation for large GPS reference networks. GPS Solut. 2014, 18, 61–71. [Google Scholar] [CrossRef]
  35. Gu, S.; Lou, Y.; Shi, C.; Liu, J. BeiDou phase bias estimation and its application in precise point positioning with triple-frequency observable. J. Geod. 2015, 89, 979–992. [Google Scholar] [CrossRef]
  36. Yang, Y.; Li, J.; Wang, A. Preliminary assessment of the navigation and positioning performance of BeiDou regional navigation satellite system. Sci. China Earth Sci. 2014, 57, 144–152. [Google Scholar] [CrossRef]
  37. Li, X.; Ge, M.; Dai, X.; Ren, X.; Fritsche, M.; Wickert, J.; Schuh, H. Accuracy and reliability of multi-GNSS real-time precise positioning: GPS, GLONASS, BeiDou, and Galileo. J. Geod. 2015, 89, 607–635. [Google Scholar] [CrossRef]
  38. Liu, T.; Yuan, Y.; Zhang, B.; Wang, N.; Tan, B.; Chen, Y. Multi-GNSS precise point positioning (MGPPP) using raw observations. J. Geod. 2017, 91, 253–268. [Google Scholar] [CrossRef]
  39. Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 2001, 36, 261–275. [Google Scholar] [CrossRef]
  40. Bent, R.; Llewellyn, S.; Schmid, P.A. Highly Successful Empirical Model for the Worldwide Ionospheric Electron Density Profile; DBA Systems: Melbourne, FL, USA, 1972. [Google Scholar]
  41. Klobuchar, J.A. Ionospheric time-delay algorithm for single frequency GPS users. IEEE Trans. Aerosp. Electron. Syst. AES 1987, 23, 325–331. [Google Scholar] [CrossRef]
  42. Tu, R.; Zhang, H.; Ge, M.; Huang, G. A real-time ionospheric model based on GNSS Precise Point Positioning. Adv. Space Res. 2013, 52, 1125–1134. [Google Scholar] [CrossRef]
  43. Bock, H.; Hugentobler, U.; Beutler, G. Kinematic and dynamic determination of trajectories for low Earth satellites using GPS. In First CHAMP Mission Results for Gravity, Magnetic and Atmospheric Studies; Reigber, C., Luhr, H., Schwintzer, P., Eds.; Springer: Heidelberg, Germany, 2003; pp. 65–69. [Google Scholar]
  44. Dach, R.; Hungentobler, U.; Fridez, P.; Michael, M. Bernese GPS Software Version 5.0; Astronomical Institute, University of Bern: Bern, Switzerland, 2007. [Google Scholar]
  45. Schmid, R.; Steigenberger, P.; Gendt, G.; Ge, M.; Rothacher, M. Generation of a Consistent Absolute Phase Center Correction Model for GPS Receiver and Satellite Antennas; Springer: Berlin, Germany, 2007. [Google Scholar]
  46. Teunissen, P.J.G. The least squares ambiguity decorrelation adjustment: A method for fast GPS integer estimation. J. Geod. 1995, 70, 65–82. [Google Scholar] [CrossRef]
  47. Janssen, V.; McElroy, S. Virtual RINEX: Science or fiction? Position 2013, 67, 38–41. [Google Scholar]
  48. Dabove, P.; Cina, A.; Manzino, A.M. How reliable is a Virtual RINEX? In Proceedings of the 2016 IEEE/ION Position, Location and Navigation Symposium (PLANS), Savannah, GA, USA, 11–14 April 2016; pp. 255–262.
Figure 1. The processing flow of integrated PPP/RTK algorithm (“Rov.” represents the rover stations, “Obs” represents the observation domain, “SSR” denotes the state space represents, “OMC” is the observed minus the computed corrections, “Caster” is used to receive and broadcast the global and the regional products).
Figure 1. The processing flow of integrated PPP/RTK algorithm (“Rov.” represents the rover stations, “Obs” represents the observation domain, “SSR” denotes the state space represents, “OMC” is the observed minus the computed corrections, “Caster” is used to receive and broadcast the global and the regional products).
Sensors 17 00507 g001
Figure 2. The distribution of the collected reference and rover stations (The reference stations are shown in red, and the rover station is shown in black).
Figure 2. The distribution of the collected reference and rover stations (The reference stations are shown in red, and the rover station is shown in black).
Sensors 17 00507 g002
Figure 3. The positioning errors of the BDS PPP/URTK ((a,b) represents the static and dynamic results by PPP mode; and (c,d) represents the static and dynamic results by URTK mode).
Figure 3. The positioning errors of the BDS PPP/URTK ((a,b) represents the static and dynamic results by PPP mode; and (c,d) represents the static and dynamic results by URTK mode).
Sensors 17 00507 g003
Figure 4. The convergence of the PPP/RTK ((a,b) represents the results of PPP and URTK, respectively).
Figure 4. The convergence of the PPP/RTK ((a,b) represents the results of PPP and URTK, respectively).
Sensors 17 00507 g004
Figure 5. The comparison of observation residuals. (a) The residuals of the L1 and L2-frequency carrier phase observations by PPP mode; (b) The residuals of the L1 and L2-frequency carrier phase observations by URTK mode.
Figure 5. The comparison of observation residuals. (a) The residuals of the L1 and L2-frequency carrier phase observations by PPP mode; (b) The residuals of the L1 and L2-frequency carrier phase observations by URTK mode.
Sensors 17 00507 g005
Figure 6. (a) The time series of VTEC at station “NT01” on day of 080 for the grid ionosphere model (black) and the PPP solution (red); (b) The differences of VTEC between the grid ionosphere model and the PPP solution for the four selected stations.
Figure 6. (a) The time series of VTEC at station “NT01” on day of 080 for the grid ionosphere model (black) and the PPP solution (red); (b) The differences of VTEC between the grid ionosphere model and the PPP solution for the four selected stations.
Sensors 17 00507 g006
Figure 7. The DCB time series of (a) P1C1 and (b) P1P2 for the four selected stations.
Figure 7. The DCB time series of (a) P1C1 and (b) P1P2 for the four selected stations.
Sensors 17 00507 g007
Table 1. The classify of the estimated parameters for PPP based on raw observations.
Table 1. The classify of the estimated parameters for PPP based on raw observations.
SystemGroup 1Group 2Group 3No.
BDS X , Y , Z , T , A m b L 1 ( N ) I ( N ) 3 × N + 6
D C B , c l k A m b L 2 ( N )
Table 2. The classify of the estimated parameters for DD baseline solution.
Table 2. The classify of the estimated parameters for DD baseline solution.
SystemGroup 1Group 2Group 3No.
BDS d X , d Y , d Z A m b L 1 ( N 1 ) δ I ( N 1 ) 3 × ( N 1 ) + 4
δ T A m b L 2 ( N 1 )
Table 3. Statistics of the positioning accuracy.
Table 3. Statistics of the positioning accuracy.
StrategyComponentsStatic (m)Dynamic (m) RMS
PPPNorth0.0080.025
East0.0060.016
Up0.0150.056
URTKNorth0.0060.021
East0.0050.017
Up0.0140.030
Table 4. The statistics of the average convergence time (Unit: s).
Table 4. The statistics of the average convergence time (Unit: s).
DOY080081082083 Average
PPP1460152015401485
RTK1.61.71.81.7

Share and Cite

MDPI and ACS Style

Tu, R.; Zhang, R.; Lu, C.; Zhang, P.; Liu, J.; Lu, X. A Unified Model for BDS Wide Area and Local Area Augmentation Positioning Based on Raw Observations. Sensors 2017, 17, 507. https://doi.org/10.3390/s17030507

AMA Style

Tu R, Zhang R, Lu C, Zhang P, Liu J, Lu X. A Unified Model for BDS Wide Area and Local Area Augmentation Positioning Based on Raw Observations. Sensors. 2017; 17(3):507. https://doi.org/10.3390/s17030507

Chicago/Turabian Style

Tu, Rui, Rui Zhang, Cuixian Lu, Pengfei Zhang, Jinhai Liu, and Xiaochun Lu. 2017. "A Unified Model for BDS Wide Area and Local Area Augmentation Positioning Based on Raw Observations" Sensors 17, no. 3: 507. https://doi.org/10.3390/s17030507

APA Style

Tu, R., Zhang, R., Lu, C., Zhang, P., Liu, J., & Lu, X. (2017). A Unified Model for BDS Wide Area and Local Area Augmentation Positioning Based on Raw Observations. Sensors, 17(3), 507. https://doi.org/10.3390/s17030507

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