Next Article in Journal
Short-Term and Long-Term Forecasting for the 3D Point Position Changing by Using Artificial Neural Networks
Previous Article in Journal
Generative Street Addresses from Satellite Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Multiple GNSS Real-Time SSR Products from Different Analysis Centers

1
Academy of Opto-Electronics, Chinese Academy of Sciences, 9 Dengzhuang South Road, Haidian District, Beijing 100094, China
2
University of Chinese Academy of Sciences, 19A Yuquan Road, Shijingshan District, Beijing 100049, China
3
School of Environment and Spatial Informatics, China University of Mining and Technology, 1 Daxue Road, Xuzhou 221116, China
*
Authors to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2018, 7(3), 85; https://doi.org/10.3390/ijgi7030085
Submission received: 15 January 2018 / Revised: 24 February 2018 / Accepted: 7 March 2018 / Published: 8 March 2018

Abstract

:
The real-time State Space Representation (SSR) product of the GNSS (Global Navigation Satellite System) orbit and clock is one of the most essential corrections for real-time precise point positioning (PPP). In this work, the performance of current SSR products from eight analysis centers were assessed by comparing it with the final product and the accuracy of real-time PPP. Numerical results showed that (1) the accuracies of the GPS SSR product were better than 8 cm for the satellite orbit and 0.3 ns for the satellite clock; (2) the accuracies of the GLONASS (GLObalnaya NAvigatsionnaya Sputnikovaya Sistema) SSR product were better than 10 cm for orbit RMS (Root Mean Square) and 0.6 ns for clock STD (Standard Deviation); and (3) the accuracies of the BDS (BeiDou Navigation Satellite System) and Galileo SSR products from CLK93 were about 14.54 and 4.42 cm for the orbit RMS and 0.32 and 0.18 ns for the clock STD, respectively. The simulated kinematic PPP results obtained using the SSR products from CLK93 and CLK51 performed better than those using other SSR products; and the accuracy of PPP based on all products was better than 6 and 10 cm in the horizontal and vertical directions, respectively. The real-time kinematic PPP experiment carried out in Beijing, Tianjin, and Shijiazhuang, China indicated that the SSR product CLK93 from Centre National d’Etudes Spatiales (CNES) had a better performance than CAS01. Moreover, the PPP with GPS + BDS dual systems had a higher accuracy than those with only a GPS single system.

1. Introduction

Precise Point Positioning (PPP) is one of the most widely-used approaches for high-precision real-time positioning with the development of multi-frequency global navigation satellite systems (GNSS). However, the PPP approach relies heavily on the availability of the high-precision satellite orbit and clock corrections [1,2,3,4,5,6,7,8,9]. Currently, the International GNSS Service (IGS) agency and various analysis centers (ACs) provide users with precise satellite orbit and clock products through FTP (File Transfer Protocol) in three forms: ultra-rapid, rapid, and final [10,11,12]. The rapid and final orbit/clock products are available after around 17 h after the end of the previous UTC (Coordinated Universal Time) day and 13 days after the end of the solution week, respectively, which mean that they cannot be used for real-time applications [13]. Although ultra-rapid products are available for real-time applications, its accuracy is not good enough for high-precision PPP. Each ultra-rapid orbit file usually covers 48 h, but only the first 24 h of the orbit are generated using actual observations and the second 24 h are extrapolated using the first 24-h orbit.
To meet the growing needs for real-time high-precision positioning and application, IGS founded a Real-Time Working Group in 2002 committed to the construction of infrastructure and to set up standards as well as technical specifications related to high-precision real-time GNSS [14]. In 2007, IGS started the Real-Time Pilot Project (RTPP) and has extended its capability to support applications requiring real-time access to IGS products since 2013 by providing GPS and GLONASS dual-system orbit and clock corrections based on RTCM (Radio Technical Commission for Maritime Services) and NTRIP (Networked Transport of RTCM via Internet Protocol) [15,16,17,18,19,20,21]. Multi-GNSS real-time orbit and clock products are also making headways with the development of BDS and Galileo. Currently, there are a wide collection of real-time orbit and clock products for either GPS or GPS + GLONASS, developed by ACs such as BKG (Bundesamt für Kartographie und Geodäsie), CNES (Centre National d’Etudes Spatiales), DLR (Deutsches Zentrum für Luft- und Raumfahrt), ESA (European Space Agency), GFZ (Deutsches GeoForschungsZentrum), and GMV (GMV Aerospace and Defense). CNES was the first to provide RTS for all four systems (GPS/GLONASS/BDS /GALILEO) since 2015 (IGS MAIL 7183).
Moreover, China started the construction of the international GNSS Monitoring and Assessment System (iGMAS) in 2012. The main task of iGMAS was to (1) establish a worldwide near-real-time tracking network for BDS, GPS, GLONASS, and Galileo; (2) build an information service platform for data collection, storage, analysis, management, and publication; (3) monitor and assess the operation status and key performance indicators of all GNSS [22]. At present, there are 30 global tracking stations, three data centers and eight analysis centers that can provide precise products to support satellite navigation technology testing, monitoring assessment, scientific research, and various applications [23].The Institute of Geodesy and Geophysics (IGG) of the Chinese Academy of Sciences (CAS), one of the iGMAS analysis centers, is also beginning to develop multi-system real-time orbit and clock correction products with the advance of iGMAS.

2. The Acquisition of Real-Time Observation Data and State Space Representation Product

It is critical for the real-time PPP to access real-time data and SSR products in an efficient way [8,11,16,18,21,24]. Accessing GNSS data via the Internet based on NTRIP has been widely used in many applications. For instance, it has been adopted in data transmission between CORS (Continuously Operating Reference Stations) servers and receivers. The NTRIP agreement, which officially became an RTCM standard in November 2004, is used for sending data streams in the format of RTCM 2.0 and 3.0. Real-time orbit and clock correction data generated by ACs in IGS-RTPP are released in SSR (State Space Representation) format in compliance with the RTCM standard and broadcasted via NTRIP (RTCM 2011). Figure 1 shows the broadcasting, receiving, and precise positioning process of real-time GNSS data/products. The GNSS data transmission system based on NTRIP generally consists of four parts: the data source, server (NtripServer), broadcaster (NtripCaster), and client terminal (NtripClient). The table of data sources generated by NTRIP broadcasters contains general information about data sources including their ID, RTCM version, data type, etc. One can access this table via the Internet on a client terminal, and select proper mount points to obtain raw data or corrections from NTRIP data sources with a short latency [25].
BKG Ntrip Client (BNC) is one of the most widely-used software package for obtaining real-time data and products [20], but it only supports data decoding in the RTCM format, making it difficult to broadcast real-time products in the latest format. In light of this, the GNSS research group at the IGG of Chinese Academy of Sciences has developed an alternative software named IGG-Ntrip, which presents the following features:
  • Supports both RTCM and iGMAS formats;
  • Supports real-time data and products for four systems (GPS/GLONASS/BDS/GALILEO) with multiple addresses and mount points;
  • Provides a data sharing mechanism based on sharing memory and Socket;
  • A user-friendly graphic interface that allows users to select stations on a map;
All the real-time data and products used in this work were obtained via IGG-Ntrip.

3. Real-Time Precise Orbit and Clock Recovery

The real-time data streams from IGS/ACs provide corrections of orbits and clocks to broadcast ephemeris. As above-mentioned, these corrections are essential for obtaining high-precision orbits and clocks for precise real-time PPP [10,26,27,28].
Real-time orbit corrections are provided in radial, along-track, and cross-track directions in a satellite-fixed coordinate system. Thus, it was necessary to first convert orbit corrections into an Earth-Fixed reference frame (ECEF) system, which was adopted for the positioning. RTCM-SSR corrections included the following parameters [29,30,31]:
Δ ssr ( t 0 , I O D E ) = ( δ O r , δ O a , δ O c , δ O r , δ O a , δ O c , C 0 , C 1 , C 2 ) ,
where IODE is the issue of data; δ O r , δ O a , δ O c , δ O r , δ O a , δ O c are the corrections and rate of change at time t 0 in the radial, along-track, and cross-track directions, respectively; C 0 , C 1 , C 2 are the polynomial coefficients for calculating clock corrections, respectively.
The orbit correction δ = [ δ r δ a δ c ] T at epoch t can be calculated using Equation (2) based on the SSR product
δ = [ δ r δ a δ c ] = [ δ O r δ O a δ O c ] + [ δ O r δ O a δ O c ] ( t t 0 )
After converting the corrections from the satellite-fixed system to the ECEF system in X, Y, and Z directions, the precise orbit R can be calculated using Equation (3),
R = [ X Y Z ] = [ x b y b z b ] [ e r e a e c ] [ δ r δ a δ c ]
where e r , e a , e c are the unit vectors in the radial, along-track, and cross-track directions, respectively; x = [ x b y b z b ] T is the satellite orbit calculated from the broadcast ephemeris.
It should be noted that there are two types of reference points for satellite position corrections in the SSR messages provided by NTRIP: APC (antenna phase center) and COM (center of mass). The data source table provided by the NTRIP casters indicates which reference point is used. If APC is adopted, the antenna phase bias correction shall be taken into account to obtain the coordinates of the satellite’s center of mass under ITRF.
Regarding the recovery of precise clocks, Equation (4) is used to determine the clock correction Δ k at epoch t with the polynomial coefficients C 0 , C 1 , C 2 given at the reference epoch t 0 .
Δ k = C 0 + C 1 ( t t 0 ) + C 2 ( t t 0 ) 2
Then, the precise satellite clock Δ t s at epoch t can be calculated with the following equation,
Δ t s = Δ t b Δ k V c
where V c is the speed of light in the vacuum; and Δ t b is the clock correction generated from the broadcast ephemeris.

4. Assessment of the Real-Time Orbit and Clock Corrections for Multi-GNSS

The precision of orbit and clock corrections is a critical issue for high-precision positioning using PPP. In this study, the satellite orbits and clocks were calculated epoch-by-epoch using the above-mentioned method with the multi-GNSS SSR products provided by eight selected ACs. The results were then compared against the final precise products released by IGS/ACs to assess the accuracy of these multi-system real-time orbit and clock corrections.
The real-time products used in this study were developed by IGS, BKG, DLR, ESA, GFZ, GMV, CNES, and CAS, as listed in Table 1 in detail. It can be seen from Table 1 that all products supported GPS, while the products CLK93 and CAS01 supported BDS and Galileo. The data used in this experiment were collected from 12 to 18 September 2017 with a sampling interval of 30 s. The final products released by IGS and ESA were selected as references for assessing the performance of real-time SSR products for GPS and GLONASS systems, respectively, while for BDS and Galileo, the final products of GBM (Geodetic Benchmark) released by GFZ were used as references.
The RMS of differences between orbits calculated using real-time corrections and final products was calculated for the radial (R), along-track (A), and cross-track (C) directions in a satellite-fixed coordinate system with the equation
RMS = i = 1 n Δ i 2 n
where Δ i represents the orbit mutual differences in nodes; and n is the number of mutual differences.
For clock comparison, the RMS and STD of the differences between the real-time clock and the reference clock are generally calculated for the assessment after unifying statistical criteria, where RMS reflects the compliance of clock correction and pseudo-range, while STD represents the real resolution precision of clock correction which has a great influence on processing phase data [21,32]. Meanwhile, we considered the mean of the clocks for all the satellites as zero as the datum instead of selecting a reference satellite as the datum [3,28], which may cause the loss of precision for the reference satellite. The adopted mathematical approach is as follows:
{ s u m r t i = j = 1 j = s c l k r t , i j s u m p t i = j = 1 j = s c l k p t , i j Δ t i = s u m r t i s u m p t i s { Δ c l k i j = c l k r t , i j c l k p t , i j Δ t i Δ j ¯ = ( i = 1 i = n Δ c l k i j ) / n R M S j = i = 1 i = n ( Δ c l k i j ) 2 n S T D j = i = 1 i = n ( Δ c l k i j Δ j ¯ ) 2 n
where i and j are the epoch and the satellite number; n is the sum of epochs; s is the number of satellites; c l k r t , i j and c l k p t , i j are the real-time and post-processing reference clock corrections; Δ t i is the datum difference between the two clock correction products with the selected datum; R M S j and S T D j are the statistical indicators of RMS and STD, respectively.

4.1. Validation Results of GPS Real-Time Orbit and Clock Products

Using the IGS final product as the reference, the RMS of the differences of the GPS real-time orbit products from the eight selected ACs are shown in Figure 2, where the RMS of the orbit differences were calculated in the R, A, and C directions, respectively, as denoted by the red, green, and blue bars; the RMS and STD of the differences of the GPS real-time clock products from the eight selected ACs are shown by the blue and red bars in Figure 3. The average accuracies of the orbit and clock corrections for each product over all the GPS satellites are given in Table 2.
It can be seen from Table 2 that the accuracies of the GPS orbit product in the radial direction was much better than that in the along-track and cross-track directions. The CLK51 performed the best in the 1D-RMS comparison, i.e., 2.09, 3.29, and 2.74 cm in the radial, along-track, and along-cross directions, respectively, while the IGS03 was the worst, i.e., 5.79, 8.04, and 7.02 cm, again in those three directions. Regarding the GPS real-time clock product, all the RMS and STD were better than 0.5 and 0.2 ns, except for CLK10 and CAS01. The real-time GPS clock product from CLK80 was the best with an RMS of 0.31 ns and STD of 0.17 ns, while the CLK10 was the worst one with an RMS of 0.79 ns and STD of 0.29 ns. Overall, the real-time GPS orbit and clock product from CLK51 and CLK93 basically had the same accuracy, i.e., the 1D-RMS of the orbit product was better than 3.0 cm and the STD of the clock product was better than 0.15 ns. The differences in accuracies between the different ACs may have resulted because of the different distributions of contributed stations and the different strategies adopted by CLK51 and CLK93; however, these differences generally do not affect the real-time PPP achieving a positioning with about a 5–10 cm accuracy. Aside from the accuracy of the products from different ACs, it should be noted that the latency of the real-time SSR corrections is also an important aspect that should be taken into account in the quality assessment since the positioning accuracy will decrease with the increasing latency of the real-time SSR corrections [33]. It can be seen from Table 3, that the CLKxx products generally had a latency of about 4–10 s, but the latency for the IGS03 products was around 30 s. Meanwhile, we also analyzed different generations of GPS satellites (BLOCKIIR-A, BLOKII-B, BLOKIIR-M, BLOKIIF) using the same product, but we did not find anything particularly regular. Thus, we were able to conclude that the accuracy of orbits and clock products had no obvious correlation with different satellite generations.

4.2. Validation Results of GLONASS Real-Time Orbit and Clock Products

Different from the real-time GPS orbit and clock products, the GLONASS products are only broadcasted in IGS03, CLK20, CLK70, CLK80, CLK93, and CAS01. The performance of GLONASS products is validated by comparing them with the ESA final product. The accuracies of GLONASS orbits are also given in the radial, along-track, and along-cross directions shown in Figure 4, and the accuracies of the GLONASS clock products from the selected ACs are given in Figure 5. The average accuracies of the GLONASS satellite orbit and clock products during the test period are illustrated in Table 4. It should be noted that the satellite orbits of R12 and R14 were absent in the products from IGS03, CLK70, CLK 80, CLK93, and CAS01 in the test period, and the result of R8 and R12 were also absent in the products from CLK20 and CLK80, respectively.
It can be seen from Table 4 that the average accuracy (1D-RMS) of the real-time GLONASS satellite orbit from all the ACs was almost at the same level (about 5–8 cm) except for the product from CAS01. Except for IGS03 and CAS01, the STDs of the real-time GLONASS satellite clock products from the other ACs were about 0.2–0.3 ns. Regarding the RMS of the GLONASS satellite clock, although the different references in the clock products from different ACs was considered before comparison, all the products performed much worse than that of the GPS satellite except for the products from CLK93, leaving large room for improvement. In particular, the RMS of IGS03 was 9.38 ns, and even higher for R10 (26 ns). The reason may be that GLONASS adopts FDMA and different ACs adopt different strategies in dealing with pseudo-range IFB in resolving clock corrections [34,35]. In addition, the GLONASS product from CAS01 is calculated by the Institute of Geodesy and Geophysics and broadcasted by the Academy of Opto-Electronics. However, it is currently in the initial and test phase and the performance of the GLONASS product will be further improved. At the same time, the different generations of GLONASS satellite (GLONASS-M and GLONASS-KI) were analyzed by the same product, and we also found nothing regular between them.

4.3. Validation Results of BDS and Galileo Real-Time Orbit and Clock Products

The real-time BDS and Galileo orbit and clock products are currently only broadcast by CLK93. The real-time BDS and Galileo products were validated by comparing them with the final precise product “GBM” released by GFZ. The validation results of the BDS/Galileo orbit and clock product are shown in Figure 6 and Figure 7, respectively; and their average accuracy during the test period is illustrated in Table 5. In the test period, the BDS constellation consisted of five GEO, six IGSO- and three MEO-satellites, and the accuracies of the real-time or final orbit and clock product for those GEO satellites were generally much worse than those of other satellites due to the poor geometry of the GEO satellites tracked by the ground stations in orbit and clock determination. Thus, the results for the BDS IGSO (C06–C10, C13) and MEO (C11, C12, C14) satellites were only considered in our experiment. In view of the Galileo constellation, there are 17 satellites in orbit, of which 15 satellites are marked as available in the ephemeris, and CLK93 estimates the orbit and clock only for those available satellites.
It can be seen from Figure 6 that the real-time Galileo orbit was more accurate than that of BDS, but a bit worse than that of GPS, and the accuracy of the BDS MEO satellites was better than that of the BDS IGSO satellites. Regarding the real-time clock accuracy shown in Figure 7, the Galileo satellites also outperformed the BDS satellites. It can be found in Table 5 that the differences of BDS satellite orbits between the CLK93 real-time product and the GBM final product in the radial, along-track, and along-cross directions were individually 7.83, 14.75, and 18.85 cm. The average STD of real-time BDS satellite clock product was 0.32 ns, and its average RMS reached 3.00 ns where the RMS of MEO satellites was larger than that of the IGSO satellites. Moreover, the differences of the Galileo satellite orbits between the real-time CLK93 product and the GBM final product in the radial, along-track, and along-cross directions were 3.21, 5.39, and 4.41 cm, respectively. The average STD of clock correction was 0.18 ns, while the average RMS was 0.39 ns. All in all, the performance of the Galileo satellite orbits and clocks was much better than that of BDS in the CLK93 product. The primary reason is that the number of contributed ground stations that can track BDS satellites is much smaller than that of Galileo.

5. Positioning Results of Real-Time PPP with IGS Stations

The GNSS satellite orbit and clock products are generally used for the real-time precise point positioning. In this section, we further validated the performance of real-time products by the PPP technique. The PPP software used in this experiment, named RTPosNavi_AOE (version 1.0.17a, Beijing, China), was developed by the Academy of Opto-Electronics (AOE) based on the open-source software RTKLIB (version 2.4.3, Japan) [19]. The PPP rover used for this validation was conducted in static and kinematic modes, respectively, but the data processing strategy was all designed as kinematic mode, i.e., the estimated parameters of rover coordinates in two neighboring epochs were considered as independent. The PPP observation model was the ionosphere-free (IF) combination of dual-frequency carrier phase and pseudo-range observations, and the Kalman filter was adopted for parameter estimation. The specific processing strategy is listed in Table 6. To evaluate the accuracy of real-time PPP, the real-time data stream in RTCM format from 10 stations from the IGS/MGEX network was selected in this study. The distribution of the selected stations for the experiment is shown in Figure 8.

5.1. Result of Real-Time PPP in Static Mode

In this section, the selected reference stations were tested for static mode. RTPosNavi_AOE received the real-time raw data and the GNSS satellite orbit and clock products via IGG-Ntrip software carried out the PPP result. Taking the case of the CEDU station as an example, our experiment began at 00:00:00 (UTC) on 16 February 2018 and lasted for 20 h. The time series of the differences between the PPP-estimated and IGS-released coordinates in the E (East-West), N (North-South), and U (Up-Down) directions are shown in Figure 9, where the convergence time was usually less than 30 min, and the positioning accuracy after finishing the convergence process was almost better than 10 cm in the E, N, and U directions.
Table 7 shows the average positioning accuracies of the static mode in the E, N, and U directions during all test periods (without the convergence time) for all selected stations. Regarding 3D positioning accuracy, CLK93 had the best performance (about 4.18 cm), followed by CLK51, IGS03, CAS01, CLK70, CLK80, CLK10, and CLK20. Thus, we concluded that the real-time PPP could achieve a positioning with an accuracy of about 6 cm, provided one was able to obtain the real-time orbit and clock products from the current eight selected ACs.

5.2. Result of Real-Time PPP in Simulated Kinematic Mode

In this section, the selected reference stations were tested for simulated kinematic mode. Since the processing strategy was designed as a kinematic mode, the positioning result with the data from a static station can be generally considered as the best performance the real kinematic mode can achieve. This is because the effect of cycle slip, multi-path, and the interruption of signals can be reduced as much as possible.
RTPosNavi_AOE receives the real-time raw data and the GNSS satellite orbit and clock products via the IGG-Ntrip software and carries out the PPP result. Taking CEDU station as an example, our experiment began at 00:00:00 (UTC) on 12 September 2017 and lasted for 20 h. IGS has released the known coordinates of every station and we directly used them as the reference for validating the positioning accuracy. The time series of the differences between the PPP-estimated and IGS-released coordinates in the E, N, and U directions are shown in Figure 10, where the results based on the different real-time GNSS satellite orbit and clock products are individually illustrated by each sub-figure. The ranges of the vertical axis of each sub-figure were set as the same for convenience during comparison, and the horizontal axis was the hour from the experiment start time. With the exception of the real-time product from IGS03 and CLK80, the convergence time was usually less than 30 min, and the positioning accuracy after finishing the convergence process was almost better than 15 cm in the E, N, and U directions.
Table 8 shows the average positioning accuracies of kinematic mode in the E, N, and U directions during all test periods (without the convergence time). The positioning result based on CNES CLK93 and ESA CLK51 had the best accuracy in the horizontal (about 4 cm) and vertical (about 6 cm) directions, respectively. Regarding 3D positioning accuracy, the CLK51 and CLK93 also had the best performance (about 6.90 cm), followed by CLK10, CLK20, IGS03, CLK80, CLK70, and CAS01. It should be noted that the number of systems provided by different ACs was different, only the product from CNES CLK93 included all GPS, BDS, GLONASS, and Galileo systems. Nevertheless, we concluded that real-time PPP could achieve a positioning with an accuracy of about 10 cm provided one was able to obtain the real-time orbit and clock products from the current eight selected ACs.

6. Positioning Result of Real-Time Kinematic PPP in Urban Experiment

In this section, we conducted a real kinematic test in the cities of Beijing, Tianjin, and Shijiazhuang in China for a better understanding of the impact of real-time orbit and clock products from different ACs on the accuracy of PPP. Based on the analysis result, we found that the performance of the real-time GPS and BDS orbit (CLK93) from CNES was the best one among the eight selected ACs and is the only center that can simultaneously broadcast GPS, BDS, GLONASS, and Galileo products. Moreover, the BDS orbit and clock product began broadcasting by CAS01 from October 2017 in the test phase. Thus, the real-time satellite orbit and clock product from CLK93 and CAS01 were only selected in this experiment, as the test result could indicate the different performances of the products from those two ACs.
The terminal used for this real-time kinematic PPP experiment was assembled by integrating a NovAtel GNSS receiver board 618 and an Advanced RISC Machines (ARM) for the PPP calculation; a 4G mobile communication module was used for receiving the real-time orbit and clock product from the servers. The NovAtel GNSS receiver board can track the GPS, BDS, GLONASS, and Galileo signals and output the raw observations with an interval of 1 s. The PPP software RTPosNavi_AOE was installed on the corresponding ARM and the positioning result was sent to our servers for showing on a precise map. The terminal and precise map are given in Figure 11. Moreover, we used two terminals in our test, one was only for the GPS mode and the other for the GPS+BDS mode. Two copies of RTPosNavi_AOE were run on the ARM to use the real-time products from CLK93 and CAS01, respectively. To obtain the real coordinates for the positioning accuracy validation, the IGS final orbit and clock product were introduced for PPP resolution in a post-processing mode.
The test was carried out in Beijing, Tianjin, and Shijiazhuang individually during the periods of 10:00:00 to 12:30:00 (UTC) on 22 October 2017, 15:40:00 to 17:00:00 (UTC) on 24 October 2017 and 13:00:00 to 14:40:00 (UTC) on 26 October 2017. The maximum speed of the experimental vehicle was about 80 km/h. The experiment routes in these three cities included the urban, expressway, and flyover environment and are individually displayed in Figure 12. The distances of these three routes were about 23.5, 34.2, and 26.1 km, respectively.
Figure 13 shows the time series of the differences between the real-time PPP-estimated and post PPP at the different cities in the E (East-West), N (North-South), and U (Up-Down) directions. The results were based on the real-time products from CLK93 where the convergence time was usually also less than 30 min, and the positioning accuracy after finishing the convergence process was almost better than 0.5 m in all directions. The accuracy of the real-time kinematic PPP result based on the products from CLK93 and CAS01 in GPS and GPS + BDS modes are illustrated in Table 9, where the result during the convergence period is not removed. It was found that the positioning accuracy of the simulated kinematic PPP decreased significantly with respect to the static result in Section 5.2 regardless of which product was used. The average accuracies of the real-time kinematic PPP based on the CLK93 product were about 0.45–0.65 m and 0.75–1.5 m in the horizontal and vertical directions, respectively, whereas they were about 0.50–0.65 m and 0.70–1.20 m for CAS01, except for the horizontal accuracy in Tianjin. The reason why the positioning result based on CAS01 in Tianjin was much worse than that in other cities needs to be further investigated.
Moreover, comparing the GPS result with that of GPS + BDS, the positioning accuracies were all slightly improved in the horizontal and vertical directions regardless of which product was used. The average accuracies in the Beijing, Tianjin, and Shijiazhuang areas are also given in Table 9. We found that the real-time kinematic PPP could achieve the positioning with an accuracy of better than 0.6 and 1.0 m in the horizontal and vertical directions based on the real-time orbit and clock product from CLK93, while they were 0.95 and 1.0 m for CAS01. The product from CAS needs to be improved further.

7. Conclusions

The SSR product including the real-time orbit and clock corrections is an essential and indispensable product for the GPS/BDS/GLONASS/Galileo PPP user achieving precise point positioning in real-time mode. In this study, a brief introduction of the methods to use real-time products for correcting satellite orbit and clock errors was summarized and the accuracy of real-time orbit and clock products from eight selected ACs were validated by using IGS, ESA, and GBM precise products as references. Further insights into PPP performance were provided by using static data to simulate kinematic positioning and real-time kinematic tests. The main conclusions are presented as follows:
  • A comparison of real-time precise products developed by eight selected ACs, i.e., BKG, DLR, ESA, GFZ, GMV, CNES, CAS, and IGS, with respect to the final precise products, showed that the RMS of differences for GPS between real-time orbit and IGS final orbit was better than 8 cm, whereas the STD of GPS clock differences was better than 0.3 ns. The RMS of GLONASS orbit differences between the real-time product and ESA final product was less than 10 cm, and clock STD was better than 0.6 ns, except for CAS01. The 1DRMS of orbit differences for BDS and Galileo in the CLK93 product were 14.54 and 4.42 cm, respectively, in contrast to 0.32 and 0.18 ns for clock STD.
  • Static mode and simulated kinematic mode PPP results from 10 stations in the IGS/MGEX network indicated that the PPP results using real-time products developed by BKG, DLR, ESA, GFA, GMV, CNES, and CAS needed nearly the same convergence time (about 20–30 min) to reach about 10 cm and 15 cm level accuracy as using IGS final precise ephemeris and clock products.
  • CNES and CAS products were selected for the real-time kinematic PPP tests carried out in Beijing, Tianjin, and Shijiazhuang, China. Results indicated that CLK93 products showed a higher accuracy than CAS01 in real-time kinematic PPP, horizontally and vertically, when being tested under the same conditions. Moreover, when using the same product for real-time PPP solution, the GPS + BDS dual system showed higher accuracy than the GPS single system.
  • The real-time SSR product from the eight selected ACs was advantageous to aid the real-time PPP user achieving positioning with sub-meter level accuracy (including the result during the convergence or re-convergence period) in the horizontal direction and could be used for lane identification in the near future.
It must be pointed out that the tests conducted in this paper were only based on one-week data from real-time products. Data from a longer timeframe will be used for the follow-up research, where real-time data from reference stations in different locations and different periods will be collected for systematic assessment and analysis of real-time product performance.

Acknowledgments

Partially funded by the National Natural Science Foundation of China (41674043, 41704038, 41730109), the National Key Research Program of China “Collaborative Precision Positioning Project” (No. 2016YFB050190), Beijing Nova program (xx2017042), Beijing Natural Science Foundation (8184092) and Beijing youth talent support program (2017000021223ZK13). The third author acknowledges the financial support from the China Scholarship Council. Thanks also to the IGG-Ntrip software provided by the Institute of Geodesy and Geophysics, Chinese Academy of Sciences, the original observations provided by IGS, and the multi-system SSR products provided by the analysis centers such as BKG, GMV, and CNES.

Author Contributions

Zishen Li and Liang Wang conceived and designed the experiments; Zhiyu Wang performed the experiments, analyzed the data, and wrote the paper; Liang Wang also helped to analyze the data; Zishen Li, Liang Wang, Xiaoming Wang, and Hong Yuan helped to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abdi, N.; Ardalan, A.A.; Karimi, R.; Rezvani, M.H. Performance assessment of multi-gnss real-time ppp over iran. Adv. Space Res. 2017, 59, 2870–2879. [Google Scholar] [CrossRef]
  2. Barsever, Y.; Muellerschoen, R. NASA’s global differential gps system. In Proceedings of the EGS-AGU-EUG Joint Assembly, Nice, France, 6–11 April 2003. [Google Scholar]
  3. Ge, M.; Chen, J.; Douša, J.; Gendt, G.; Wickert, J. A computationally efficient approach for estimating high-rate satellite clock corrections in realtime. GPS Solut. 2012, 16, 9–17. [Google Scholar] [CrossRef]
  4. Hauschild, A.; Montenbruck, O. Kalman-filter-based gps clock estimation for near real-time positioning. GPS Solut. 2009, 13, 173–182. [Google Scholar] [CrossRef]
  5. Huang, G.; Zhang, Q.; Li, H.; Fu, W. Quality variation of gps satellite clocks on-orbit using igs clock products. Adv. Space Res. 2013, 51, 978–987. [Google Scholar] [CrossRef]
  6. Yu, X.; Gao, J. Kinematic precise point positioning using multi-constellation global navigation satellite system (gnss) observations. Int. J. Geo-Inf. 2017, 6, 6. [Google Scholar] [CrossRef]
  7. Asari, K.; Suzuki, T.; Kubo, N.; Shima, M.; Saito, M. Ssr ppp for asia-pacific services. In Proceedings of the 27th International Technical Meeting of the Satellite Division of the Institute of Navigation, ION GNSS 2014, Tampa, FL, USA, 8–12 September 2014. [Google Scholar]
  8. Bar-Sever, Y.; Bell, B.; Dorsey, A.; Srinivasan, J. Space Applications of the Nasa Global Differential GPS System; Jpl Technical Reports; Institute of Navigation: Portland, OR, USA, 2003. [Google Scholar]
  9. Mervart, L.; Weber, G. Real-time combination of gnss orbit and clock correction streams using a kalman filter approach. In Proceedings of the International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2011), Portland, OR, USA, 20–23 September 2011; pp. 707–711. [Google Scholar]
  10. Dixon, K. Starfire™: A global sbas for sub-decimeter precise point positioning. In Proceedings of the International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2006), Fort Worth, TX, USA, 26–29 September 2006; pp. 2286–2296. [Google Scholar]
  11. Kouba, J. A Guide to Using International Gnss Service (Igs) Products; IGS Central Bureau: Pasadena, CA, USA, 2009. [Google Scholar]
  12. Melgard, T.; Vigen, E.; Jong, K.D.; Lapucha, D.; Visser, H.; Oerpen, O. G2-the first real-time gps and glonass precise orbit and clock service. In Proceedings of the International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2009), Savannah, GA, USA, 22–25 September 2009; pp. 1885–1891. [Google Scholar]
  13. IGS International GNSS Service. Available online: http://www.igs.org/products (accessed on 7 February 2018).
  14. Caissy, M.; Agrotis, L.; Weber, G.; Hernandez-Pajares, M.; Hugentobler, U. Coming Soon: The International Gnss Real-Time Service. GPS World 2012, 23, 52–58. [Google Scholar]
  15. Chen, L.; Zhao, Q.; Hu, Z.; Jiang, X.; Geng, C.; Ge, M.; Shi, C. Gnss global real-time augmentation positioning: Real-time precise satellite clock estimation, prototype system construction and performance analysis. Adv. Space Res. 2018, 61, 367–384. [Google Scholar] [CrossRef]
  16. Colombo, O. Real-time, wide-area, precise kinematic positioning using data from internet ntrip streams. In Proceedings of the 21st International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS 2008), Network Rtk Session, Savannah, GA, USA, 16–19 September 2008. [Google Scholar]
  17. Li, Z.; Zhang, J.; Li, T.; He, X.; Wu, M. Analysis of static and dynamic real-time precise point positioning and precision based on ssr correction. In Proceedings of the IEEE International Conference on Information and Automation, Macau, China, 18–20 July 2017; pp. 2022–2027. [Google Scholar]
  18. Liu, Z.Q.; Wang, J.X. Realization and analysis of real-time precise point positioning based on ssr broadcast ephemeris corrections. Sci. Surv. Mapp. 2014, 39, 15–19. [Google Scholar]
  19. Wang, L.; Li, Z.; Ge, M.; Neitzel, F.; Wang, Z.; Yuan, H. Validation and assessment of multi-gnss real-time precise point positioning in simulated kinematic mode using igs real-time service. Remote Sens. 2018, 10, 337. [Google Scholar] [CrossRef]
  20. Weber, G.; Mervart, L.; Lukes, Z.; Rocken, C.; Dousa, J. Real-time clock and orbit corrections for improved point positioning via Ntrip. In Proceedings of the ION GNSS 20th International Technical Meeting of the Satellite Division, Fort Worth, TX, USA, 25–28 September 2007; pp. 1992–1998. [Google Scholar]
  21. Zhang, B. Three methods to retrieve slant total electron content measurements from ground-based gps receivers and performance assessment. Radio Sci. 2016, 51, 972–988. [Google Scholar] [CrossRef]
  22. Jiao, W. International GNSS Monitoring and Assessment System (iGMAS) and latest progress. In Proceedings of the China Satellite Navigation Conference (CSNC), Nanjing, China, 21–23 May 2014. [Google Scholar]
  23. International GNSS Monitoring & Assessment System. Available online: http://www.igmas.org (accessed on 7 February 2018).
  24. Zumberge, J.F.; Heflin, M.B.; Jefferson, D.C.; Watkins, M.M.; Webb, F.H. 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]
  25. Gebhard, H.; Weber, G. NTRIP: Networked transport of rtcm via internet protocol—Internet radio technology for real-time gnss purposes. In Proceedings of the AGU Fall Meeting, San Francisco, CA, USA, 8–12 December 2003. [Google Scholar]
  26. Collins, P.; Gao, Y.; Lahaye, F.; Héroux, P.; Macleod, K.; Chen, K. Accessing and processing real-time gps corrections for precise point positioning … some user considerations. In Proceedings of the International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2005), Long Beach, CA, USA, 13–16 September 2005; pp. 1483–1491. [Google Scholar]
  27. Pérez, J.; Agrotis, L.; Fernández, J.; Garcia, C.; Dow, J. Esa/esoc real time data processing. In Proceedings of the IGS workshop 2006, Darmstadt, Germany, 8–11 May 2006. [Google Scholar]
  28. Zhang, X.; Li, X.; Guo, F. Satellite clock estimation at 1 hz for realtime kinematic ppp applications. GPS Solut. 2011, 15, 315–324. [Google Scholar] [CrossRef]
  29. El-Mowafy, A.; Deo, M.; Kubo, N. Maintaining real-time precise point positioning during outages of orbit and clock corrections. GPS Solut. 2017, 21, 937–947. [Google Scholar] [CrossRef]
  30. 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]
  31. Hadas, T.; Bosy, J. Igs Rts Precise Orbits and Clocks Verification and Quality Degradation over Time; Springer: New York, NY, USA, 2015; pp. 93–105. [Google Scholar]
  32. Ge, M.; Zhang, H.; Jia, X.; Song, S.; Wickert, J. What is achievable with the current compass constellation? GPS World 2012, 1, 29–34. [Google Scholar]
  33. Martín, A.; Hadas, T.; Dimas, A.; Anquela, A.B.; Berné, J.L. Influence of real-time products latency on kinematic ppp results. In Proceedings of the ESA 5th International Colloquium on Scientific and Fundamental Aspects of the Galileo Program, Braunschweig, Germany, 27–29 October 2015. [Google Scholar]
  34. Chen, Y.; Yuan, Y.; Ding, W.; Zhang, B.; Liu, T. Glonass pseudorange inter-channel biases considerations when jointly estimating gps and glonass clock offset. GPS Solut. 2017, 21, 1525–1533. [Google Scholar] [CrossRef]
  35. Song, W.; Yi, W.; Lou, Y.; Shi, C.; Yao, Y.; Liu, Y.; Mao, Y.; Xiang, Y. Impact of glonass pseudorange inter-channel biases on satellite clock corrections. GPS Solut. 2014, 18, 323–333. [Google Scholar] [CrossRef]
Figure 1. Broadcasting, receiving, and precise positioning process of real-time Global Navigation Satellite System (GNSS) data/products.
Figure 1. Broadcasting, receiving, and precise positioning process of real-time Global Navigation Satellite System (GNSS) data/products.
Ijgi 07 00085 g001
Figure 2. The RMS of the differences of real-time GPS orbit between the eight selected ACs’ products and the IGS final product.
Figure 2. The RMS of the differences of real-time GPS orbit between the eight selected ACs’ products and the IGS final product.
Ijgi 07 00085 g002
Figure 3. The Root Mean Square (RMS) (blue bar) and Standard Deviation (STD) (red bar) of the differences of real-time Global Positioning System (GPS) clock between the eight selected analysis centers’ (ACs’) products and the International GNSS Service (IGS) final product.
Figure 3. The Root Mean Square (RMS) (blue bar) and Standard Deviation (STD) (red bar) of the differences of real-time Global Positioning System (GPS) clock between the eight selected analysis centers’ (ACs’) products and the International GNSS Service (IGS) final product.
Ijgi 07 00085 g003
Figure 4. The RMS of the differences of real-time GLONASS orbit between the six selected ACs’ products and the European Space Agency (ESA) final product.
Figure 4. The RMS of the differences of real-time GLONASS orbit between the six selected ACs’ products and the European Space Agency (ESA) final product.
Ijgi 07 00085 g004
Figure 5. The RMS (blue bar) and STD (red bar) of the differences of real-time GLONASS clock between the six selected ACs’ products and the IGS final product.
Figure 5. The RMS (blue bar) and STD (red bar) of the differences of real-time GLONASS clock between the six selected ACs’ products and the IGS final product.
Ijgi 07 00085 g005
Figure 6. The RMS of the differences of the BDS and Galileo orbit between the CLK93 real-time product and the GBM final product.
Figure 6. The RMS of the differences of the BDS and Galileo orbit between the CLK93 real-time product and the GBM final product.
Ijgi 07 00085 g006
Figure 7. The RMS (blue bar) and STD (red bar) of the differences of BDS (top) and Galileo (bottom) between the CLK 93 real-time product and the GBM final product.
Figure 7. The RMS (blue bar) and STD (red bar) of the differences of BDS (top) and Galileo (bottom) between the CLK 93 real-time product and the GBM final product.
Ijgi 07 00085 g007
Figure 8. The distribution of the selected stations for the real-time PPP test experiment.
Figure 8. The distribution of the selected stations for the real-time PPP test experiment.
Ijgi 07 00085 g008
Figure 9. Time series of the differences between PPP-estimated (static mode) and IGS-released coordinates in E (East-West), N (North-South), U (Up-Down) directions at CEDU station. The results based on the real-time products from IGS03, CLK10, CLK20, CLK51, CLK70, CLK80, CLK93, and CAS01 are individually illustrated by each sub-figure.
Figure 9. Time series of the differences between PPP-estimated (static mode) and IGS-released coordinates in E (East-West), N (North-South), U (Up-Down) directions at CEDU station. The results based on the real-time products from IGS03, CLK10, CLK20, CLK51, CLK70, CLK80, CLK93, and CAS01 are individually illustrated by each sub-figure.
Ijgi 07 00085 g009
Figure 10. Time series of the differences between PPP-estimated (simulated kinematic mode) and IGS-released coordinates in E (East-West), N (North-South), U (Up-Down) directions at CEDU station. The results based on the real-time products from IGS03, CLK10, CLK20, CLK51, CLK70, CLK80, CLK93, and CAS01 are individually illustrated by each sub-figure.
Figure 10. Time series of the differences between PPP-estimated (simulated kinematic mode) and IGS-released coordinates in E (East-West), N (North-South), U (Up-Down) directions at CEDU station. The results based on the real-time products from IGS03, CLK10, CLK20, CLK51, CLK70, CLK80, CLK93, and CAS01 are individually illustrated by each sub-figure.
Ijgi 07 00085 g010
Figure 11. The positioning terminal (left) and the precise map (right) used for the real-time kinematic PPP test.
Figure 11. The positioning terminal (left) and the precise map (right) used for the real-time kinematic PPP test.
Ijgi 07 00085 g011
Figure 12. The road selected for our experiment in Beijing (left), Tianjin (middle), and Shijiazhuang (right).
Figure 12. The road selected for our experiment in Beijing (left), Tianjin (middle), and Shijiazhuang (right).
Ijgi 07 00085 g012
Figure 13. Time series of the differences between real-time PPP-estimated and post PPP at Beijing (left), Tianjin (middle), and Shijiazhuang (right) in the E, N, and U directions. The results of GPS (top) and GPS + BDS (bottom) based on the real-time products from CLK93 are only shown here.
Figure 13. Time series of the differences between real-time PPP-estimated and post PPP at Beijing (left), Tianjin (middle), and Shijiazhuang (right) in the E, N, and U directions. The results of GPS (top) and GPS + BDS (bottom) based on the real-time products from CLK93 are only shown here.
Ijgi 07 00085 g013
Table 1. The details of the real-time products used in this work.
Table 1. The details of the real-time products used in this work.
ACsIGSBKGDLRESAGFZGMVCNESCAS
ProductsIGS03CLK10CLK20CLK51CLK70CLK80CLK93CAS01
SystemG + RGG + RGG + RG + RG + R + C + EG + R + C + E
Table 2. The average accuracies of GPS real-time orbit and clock products from the eight selected ACs with respect to the IGS final products during the test period.
Table 2. The average accuracies of GPS real-time orbit and clock products from the eight selected ACs with respect to the IGS final products during the test period.
ProductsOrbits (cm)Clocks (ns)
RAC1DRMSRMSSTD
IGS035.798.047.027.010.350.19
CLK102.764.173.563.540.790.29
CLK202.483.972.943.190.430.12
CLK512.093.292.742.750.460.13
CLK705.147.466.116.310.320.18
CLK805.057.195.956.130.310.17
CLK932.193.502.792.880.450.14
CAS015.016.775.675.860.690.19
Table 3. The latency of the real-time State Space Representation (SSR) corrections from different ACs.
Table 3. The latency of the real-time State Space Representation (SSR) corrections from different ACs.
ProductsOrbit/Clock Update IntervalLatency of Orbit/Clock
CLK1060 s/5 s5 s
CLK205 s/5 s4 s
CLK515 s/5 s6 s
CLK7010 s/5 s6 s
CLK805 s/5 s6 s
CLK935 s/5 s8 s
CAS011 s/1 s8 s
IGS0360 s/10 s28 s
Table 4. The average accuracies of GLObalnaya NAvigatsionnaya Sputnikovaya Sistema (GLONASS) orbit and clock products from the six selected ACs compared with the IGS final products during the test period.
Table 4. The average accuracies of GLObalnaya NAvigatsionnaya Sputnikovaya Sistema (GLONASS) orbit and clock products from the six selected ACs compared with the IGS final products during the test period.
ProductsOrbits (cm)Clocks (ns)
RAC1DRMSRMSSTD
IGS035.348.746.957.159.380.50
CLK204.457.095.695.848.290.31
CLK706.218.557.737.561.740.17
CLK806.449.027.467.710.820.33
CLK934.157.165.325.682.970.21
CAS0110.3334.8420.5724.115.491.86
Table 5. The average accuracies of the BeiDou Navigation Satellite System (BDS)/Galileo orbit and clock products from CLK93 when compared with the GBM final product during the test period.
Table 5. The average accuracies of the BeiDou Navigation Satellite System (BDS)/Galileo orbit and clock products from CLK93 when compared with the GBM final product during the test period.
SystemOrbits (cm)Clocks (ns)
RAC1DRMSRMSSTD
BDS7.8314.7518.8514.543.000.32
Galileo3.215.394.414.420.390.18
Table 6. The processing strategy for precise point positioning (PPP) used for validating the performance of the real-time GNSS orbit and clock product.
Table 6. The processing strategy for precise point positioning (PPP) used for validating the performance of the real-time GNSS orbit and clock product.
IndexItemsProcessing Strategies
1. Rover coordinatesReal-time estimation; initial value determined by SPP
2. Ionospheric delayIonosphere-free combination
3. Satellite orbit and clockSSR corrections + Broadcast Ephemeris
4. Receiver clockReal-time estimation
5. Tropospheric delaySaastamoinen model + real-time estimation
6. Antenna PCO/PCVigs14.atx
7. Cut-off elevation angle
8. AmbiguityFloat
9. Cycle SlipDetected by MW and GF
Table 7. Static mode positioning RMS of the PPP-estimated result using the orbit and clock products from different ACs.
Table 7. Static mode positioning RMS of the PPP-estimated result using the orbit and clock products from different ACs.
StationsDirection (cm)IGS
IGS03
BKG
CLK10
DLR
CLK20
ESA
CLK51
GFZ
CLK70
GMV
CLK80
CNES
CLK93
CAS
CAS01
CEDUE1.943.753.092.633.162.912.043.82
N3.623.614.753.134.222.982.810.12
U2.912.812.572.150.333.141.201.69
3D5.045.916.224.625.285.213.674.17
AZU1E2.013.823.153.613.023.192.112.66
N3.643.654.780.993.024.212.853.16
U2.882.792.622.653.110.381.182.17
3D5.065.976.304.595.285.303.744.67
RVDIE1.953.83.163.583.173.112.152.56
N3.593.714.721.143.363.112.953.29
U2.682.732.482.352.452.971.132.67
3D4.895.976.204.435.235.313.824.95
ORIDE1.813.883.123.533.053.22.012.59
N3.493.554.31.043.034.063.263.25
U2.432.182.322.713.090.692.032.54
3D4.625.695.804.575.295.224.334.87
STJOE1.773.63.223.623.023.262.482.7
N3.253.414.131.113.683.693.513.37
U1.982.452.152.612.483.141.932.44
3D4.205.535.664.605.375.844.714.96
SALUE2.153.413.373.563.123.092.722.46
N3.913.714.360.953.213.823.443.52
U1.742.372.022.423.153.321.362.27
3D4.795.575.874.415.475.934.594.86
CAS1E1.663.253.123.473.063.152.252.32
N3.663.834.171.324.24.363.333.61
U1.852.282.222.530.910.771.452.53
3D4.425.525.664.495.285.434.274.98
CHOFE1.733.023.213.333.163.012.072.23
N3.713.994.021.253.693.143.023.85
U1.922.022.052.811.833.181.552.55
3D4.525.405.544.535.195.393.985.13
RIO2E1.873.013.123.673.193.221.942.44
N3.553.74.211.043.783.643.513.59
U1.681.932.132.361.913.021.642.47
3D4.355.155.664.495.305.724.334.99
HARBE1.913.323.173.513.163.21.832.32
N3.463.513.661.233.453.343.573.45
U1.821.732.002.872.012.951.772.29
3D4.325.175.264.705.095.504.364.79
AverageE1.883.493.173.453.113.132.162.61
N3.593.674.311.323.563.643.233.12
U2.192.332.262.552.132.361.522.36
3D4.625.595.824.545.285.484.184.84
Table 8. Simulated kinematic mode Positioning RMS of the PPP-estimated result using the orbit and clock products from different ACs.
Table 8. Simulated kinematic mode Positioning RMS of the PPP-estimated result using the orbit and clock products from different ACs.
StationsDirection (cm)IGS
IGS03
BKG
CLK10
DLR
CLK20
ESA
CLK51
GFZ
CLK70
GMV
CLK80
CNES
CLK93
CAS
CAS01
CEDUE3.362.742.243.786.024.891.755.03
N3.862.511.891.963.762.601.672.57
U5.956.687.005.087.128.316.369.53
3D7.857.657.596.6210.059.996.8110.86
AZU1E5.563.432.483.845.813.722.013.95
N5.842.021.962.641.862.601.622.28
U6.446.086.874.987.327.536.138.53
3D10.317.277.566.829.538.796.659.67
RVDIE4.833.392.882.564.994.194.454.81
N4.712.341.861.783.242.582.033.22
U6.026.716.545.996.997.765.338.56
3D9.047.877.386.759.189.197.2310.33
ORIDE5.172.682.492.815.124.133.034.11
N4.061.932.041.903.992.542.522.48
U5.985.997.135.887.418.494.899.67
3D8.886.847.826.799.859.786.2810.79
STJOE4.413.162.373.783.845.661.824.81
N4.511.701.982.072.643.791.552.85
U5.797.336.966.018.737.426.439.34
3D8.568.167.347.409.8910.076.8610.89
SALUE5.472.223.732.754.914.323.485.16
N4.061.852.861.883.672.402.332.63
U6.076.786.445.637.418.615.338.98
3D9.127.377.976.549.629.926.7810.68
CAS1E4.773.192.454.016.134.942.124.26
N3.622.211.781.882.652.591.832.52
U5.116.936.895.688.037.816.169.83
3D7.877.947.537.2010.449.606.7711.01
CHOFE5.023.684.573.135.243.992.863.91
N5.941.742.051.992.842.061.712.00
U6.746.895.725.988.018.576.7310.06
3D10.298.007.607.049.989.677.5111.71
RIO2E4.332.892.923.235.534.732.114.74
N5.152.101.762.342.351.831.971.74
U6.546.767.095.467.247.946.888.67
3D9.387.647.876.769.419.427.4610.03
HARBE4.683.142.484.015.885.012.125.65
N5.512.241.961.913.693.082.012.25
U6.026.406.796.087.377.766.098.84
3D9.417.477.497.5210.129.996.7510.72
AverageE4.763.052.863.395.354.562.584.64
N4.732.062.022.043.072.611.922.45
U6.076.656.745.677.568.026.038.84
3D9.077.627.626.949.819.646.9110.67
Table 9. The accuracies of real-time kinematic PPP using the CLK93 and CAS01 in GPS and GPS + BDS modes.
Table 9. The accuracies of real-time kinematic PPP using the CLK93 and CAS01 in GPS and GPS + BDS modes.
CityProductSystemAccuracy (m)Available Epochs
HorizontalVertical
BeijingCLK93GPS0.470.788005
GPS + BDS0.460.738023
CAS01GPS0.600.778005
GPS + BDS0.560.748023
TianjinCLK93GPS0.631.064615
GPS + BDS0.510.504758
CAS01GPS1.380.924615
GPS + BDS1.290.894758
ShijiazhuangCLK93GPS0.621.145500
GPS + BDS0.530.835591
CAS01GPS0.601.165500
GPS + BDS0.531.085591
AverageCLK93GPS0.581.0118,120
GPS + BDS0.500.7018,372
CAS01GPS0.940.9618,120
GPS + BDS0.870.9118,372

Share and Cite

MDPI and ACS Style

Wang, Z.; Li, Z.; Wang, L.; Wang, X.; Yuan, H. Assessment of Multiple GNSS Real-Time SSR Products from Different Analysis Centers. ISPRS Int. J. Geo-Inf. 2018, 7, 85. https://doi.org/10.3390/ijgi7030085

AMA Style

Wang Z, Li Z, Wang L, Wang X, Yuan H. Assessment of Multiple GNSS Real-Time SSR Products from Different Analysis Centers. ISPRS International Journal of Geo-Information. 2018; 7(3):85. https://doi.org/10.3390/ijgi7030085

Chicago/Turabian Style

Wang, Zhiyu, Zishen Li, Liang Wang, Xiaoming Wang, and Hong Yuan. 2018. "Assessment of Multiple GNSS Real-Time SSR Products from Different Analysis Centers" ISPRS International Journal of Geo-Information 7, no. 3: 85. https://doi.org/10.3390/ijgi7030085

APA Style

Wang, Z., Li, Z., Wang, L., Wang, X., & Yuan, H. (2018). Assessment of Multiple GNSS Real-Time SSR Products from Different Analysis Centers. ISPRS International Journal of Geo-Information, 7(3), 85. https://doi.org/10.3390/ijgi7030085

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