Next Article in Journal
3D Wavelet Finite-Element Modeling of Frequency-Domain Airborne EM Data Based on B-Spline Wavelet on the Interval Using Potentials
Next Article in Special Issue
Estimation of Inter-System Biases between BDS-3/GPS/Galileo and Its Application in RTK Positioning
Previous Article in Journal
Absence of Surface Water Temperature Trends in Lake Kinneret despite Present Atmospheric Warming: Comparisons with Dead Sea Trends
Previous Article in Special Issue
LEO-Constellation-Augmented BDS Precise Orbit Determination Considering Spaceborne Observational Errors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Parallel Approach for Multi-GNSS Ultra-Rapid Orbit Determination

1
School of Geomatics, Liaoning Technical University, Fuxin 123000, China
2
State Key Laboratory of Geo-Information Engineering, Xi’an 710054, China
3
Institute of Space Science, Shandong University, Weihai 214209, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(17), 3464; https://doi.org/10.3390/rs13173464
Submission received: 9 July 2021 / Revised: 29 August 2021 / Accepted: 30 August 2021 / Published: 1 September 2021
(This article belongs to the Special Issue Beidou/GNSS Precise Positioning and Atmospheric Modeling)

Abstract

:
Global Navigation Satellite System (GNSS) ultra-rapid orbit is critical for geoscience and real-time engineering applications. To improve the computational efficiency and the accuracy of predicted orbit, a parallel approach for multi-GNSS ultra-rapid orbit determination is proposed based on Message Passing Interface (MPI)/Open Multi Processing (OpenMP). This approach, compared with earlier efficient methods, can improve the efficiency of multi-GNSS ultra-rapid orbit solution without changing the original observation data and retaining the continuity and consistency of the original parameters to be estimated. To obtain high efficiency, three steps are involved in the approach. First and foremost, the normal equation construction is optimized in parallel based on MPI. Second, equivalent reduction of the estimated parameters is optimized using OpenMP parallel method. Third, multithreading is used for parallel orbit extrapolation. Thus, GNSS ultra-rapid orbit determination is comprehensively optimized in parallel, and the computation efficiency is greatly improved. Based on the data from MGEX and IGS stations, experiments are carried out to analyze the performance of the proposed approach in computational efficiency, accuracy and stability. The results show that the approach greatly improves the efficiency of satellite orbit determination. It can realize 1-h update frequency for the multi-GNSS ultra-rapid orbit determination using 88 stations with four-system observations. The accuracy of the GPS, GLONASS, Galileo and BDS ultra-rapid orbit with 1-h update frequency using the parallel approach is approximately 33.4%,31.4%,40.1% and 32.8% higher than that of the original orbit, respectively. The root mean squares (RMS) of GPS, GLONASS, Galileo and BDS predicted orbit are about 3.2 cm, 5.1 cm, 5.6 cm and 11.8 cm. Moreover, the orbit provided by the proposed method has a better stability. The precision loss of all parallel optimization can be negligible and the original correlation between the parameters is fully preserved.

1. Introduction

Global Navigation Satellite System (GNSS) ultra-rapid orbit is the basis of geodetic research and real-time applications, such as real-time atmosphere monitoring [1,2], real-time PPP [3,4,5], real-time satellite clock estimation [6], and real-time orbit determination of LEO satellite [7]. With the development of Galileo and implementation of BDS global service plan, the number of satellites in the constellation has increased dramatically to more than 130 [8]. Currently, the International GNSS Service (IGS), the Multi-GNSS Experiment and Pilot Project (MGEX) and the International GNSS Monitoring and Assessment System (iGMAS) have provided more than 500 global distributed stations. A combined multi-GNSS solution can effectively reduce systematic errors and improve the accuracy and reliability of satellite orbit. However, on the other hand, the increase in the number of satellite constellations and stations makes the processing model more complicated and results in larger computational burden [9], which decrease the computational efficiency of multi-GNSS ultra-rapid orbit determination [10]. Real-time or quasi-real-time users have higher requirements for the accuracy, stability, and update frequency of GNSS ultra-rapid orbits. Therefore, further research on the rapid determination of multi-GNSS predicted orbits is critical for its high-precision applications.
GNSS orbit products that have been widely used are gradually developing toward multi-system, higher accuracy and higher update frequency. Since its establishment in 1994, IGS has vigorously carried out research on GNSS precise orbit determination, and the accuracy and reliability of the products have been continuously improved [11,12]. At present, the accuracy of the ultra-rapid orbit is better than 5 cm, which is provided for real-time users. The details of the accuracy and latency of orbit products are shown in Table 1 (http://www.igs.org/products, accessed January 2019). IGS also generates experimental final combined multi-GNSS orbits, which include GPS, GLONASS, Galileo, BDS-2, BDS-3, and QZSS [13]. The IGS called for participation in MGEX in 2012 (http://mgex.igs.org, accessed January 2019). Subsequently, five analysis centers (ACs) provided multi-system orbit products, including CNES, TUM, GFZ, CODE and WUH. The later three ACs can provide the orbit product of GPS, GLONASS, Galileo and BDS, as well as QZSS. In recent years, ACs have increased the update frequency of ultra-rapid orbit for real-time applications (ftp://cddis.gsfc.nasa.gov/pub/gps/products/mgex/, accessed January 2019). IGS increased the update frequency from 12-h to 6-h. The frequency of ultra-rapid orbit has been further improved to 3-h by the German Research Centre for Geosciences (GFZ) (ftp.gfz-potsdam.de, accessed January 2019). The analysis center at Wuhan university is the first to provide a 1-h update frequency ultra-rapid orbit product.
The algorithm of precise orbit determination (POD) mainly uses dynamics and reduced dynamics methods [14,15,16,17,18]. The methods above are first applied to GPS satellite orbit determination and are widely used in GLONASS, Galileo, BDS and other satellite systems, as well as the combination of them with the development of GNSS [19,20,21,22,23,24]. Based on this, the predicted orbit of ultra-rapid products is generally obtained by integration of initial orbit and force model parameters [25]. The accuracy of orbit for real-time applications is directly related to the accuracy of satellite orbit initial values and the length of extrapolation time. In order to obtain accurate solutions and achieve higher computational efficiency, many scholars have optimized the theories and methods of orbit determination by dividing subnets, grouping combined systems, dividing arc length of observation, and reducing observations. Beutler (1999), for the first time, proposed the method of merging the past three consecutive single-day arcs to orbit determination in 1999 [26]. Yao (2008) further used a dynamic smoothing algorithm to smooth the orbit which can improve the accuracy and reliability of the comprehensive algorithm solution. In order to improve the efficiency of ultra-rapid orbit determination [27], Lou (2008) deeply studied the problem of ambiguity processing in real-time orbit determination and proposed a sliding window method [28]. Through the superposition of short-arc normal equations, the orbit can be calculated with higher efficiency, and the accuracy is improved when the ambiguities are fixed. In addition, TUM used a two-step approach for generating Galileo and QZSS products [29]. Furthermore, a three-step approach was applied by WHU to orbit product calculation: GPS and GLONASS orbit and ERPs are solved first, then station coordinates, zenith tropospheric delay (ZTD) and receiver clock are estimated, and the orbit of other satellites are solved in the last step. Then, a block recursive least squares method was proposed to improve the accuracy of GNSS ultra-rapid orbits. The update frequency of ultra-rapid orbits was increased from 6-h to 1-h, and the accuracy was significantly improved [30]. Li (2019) realized 1-h update multi-GNSS ultra-rapid orbit determination by multi-threading and changed observation intervals [31]. The accuracy of predicted orbit was improved and the kinematic PPP solutions could achieve a shorter time to first fix and higher positioning accuracy.
To promote the accuracy of the predicted multi-GNSS orbit, we proposed a new approach for multi-GNSS hourly ultra-rapid orbit determination. This approach improves the efficiency of multi-GNSS ultra-rapid orbit solution without changing the original observation data or processing interval, and retains the continuity and consistency of the original parameters to be estimated using the multi-process and multi-thread method based on Message Passing Interface (MPI)/Open Multi Processing (OpenMP). MPI is a message passing library specification which provides a powerful and portable way for expressing parallel programs. OpenMP is a set of compiler directives as well as an API for programs that provides support for parallel programming in shared-memory environments. Therefore, MPI/OpenMP can achieve higher accuracy for multi-GNSS ultra-rapid orbit as the initial orbit is precise, but the length of extrapolation time is shortened significantly. To realize that, first and foremost, the construction of the normal equation for initial value of satellite orbit is divided into multi-arc parallel processing by the multi-process method based on MPI, and each process achieves consistency through processing the arc data and interacting matrix. In addition, the multi-thread method is used in parameter elimination and orbit extrapolation to further improve the efficiency. In this approach, the GPS+GLONASS+Galileo+BDS observations are processed together in one common procedure with the double-differenced network model. The promotion of computational efficiency, accuracy and stability is validated based on the data from MGEX and IGS stations.
In the following sections, the principle and main steps of the modified parallel multi-GNSS obit determination approach will be introduced. Then, some processing strategies are presented. Afterward, using MGEX/IGS data and related products, the improvement of computational efficiency of the new approach is analyzed. At last, the accuracy of multi-GNSS obit using the proposed approach is evaluated, and the results indicate that the proposed parallel approach is viable.

2. Methodology

The double-differenced network solution is used in the proposed parallel multi-GNSS ultra-rapid orbit determination. After the data preprocessing, the basic procedure of multi-GNSS ultra-rapid orbit determination mainly includes three steps. First, the double-differenced normal equations are constructed based on un-differenced observation equations. Second, the normal equations are solved after eliminating the redundant parameters. Finally, the predicted orbits are obtained by extrapolation on the basis of initial orbit and model parameters using orbit integration. We optimize the above core processing to achieve higher computational efficiency using a parallel computing method based on MPI/OpenMP. Since the original observational model and stochastic model have not been changed, the approach retains the original correlation of the estimated parameters and the consistency of the ambiguity parameters.

2.1. Multi-Arc Parallel Normal Equations Construction Based on MPI

In the multi-GNSS ultra-rapid orbit determination, the sequential batch processing method is adopted. The double-differenced carrier observation equation can be expressed as:
L C k l i j ( t ) = ρ k l i j ( t ) + Δ T k l i j t r o p ( t ) + N k l i j λ c + ε Φ
where t represents epoch, L C k l i j represents the ionospheric-free phase observation of station k, l and satellite i, j, ρ k l i j denotes pseudorange, Δ T k l   t r o p i j denotes tropospheric delay, N k l i j is ambiguity, λ is the wavelength of ionospheric-free observations, and ε Φ is the observation noise. The observation equation (1) is linearized and expressed as:
L C k l i j ( t ) = ρ k l 0 i j ( t , X r 0 , X s ( t ) ) + B δ X r + H δ X 0 s + Δ T k l i j t r o p ( t ) + N k l i j λ c + ε Φ
where B = ρ k l i j ( t ) / X r is the partial derivative of the measured value with respect to the state parameters, Ψ ( t , t 0 ) is the transition matrix which can be obtained by solving the variational equation, X r 0 , δ X r are the approximate coordinates of the station and their corrections, respectively, X s ( t ) is the state of satellite at time t, which is obtained by solving the differential equation of satellite motion, ρ k l 0 i j ( t , X r 0 , X s ( t ) ) is the approximate value of the double-differenced satellite-ground geometric distance at time t, and δ X 0 s is the correction of the satellite state at the reference time. Other symbols have the same meaning as in (1). Based on the observation (2), the GNSS satellite reduce-dynamic precise orbit determination (RDPOD) can be performed. Table 2 shows the main parameter settings of the RDPOD strategy. In this paper, we have no a priori box-wing model used for the orbit determination, which may generate orbit errors, especially for Galileo satellites due to the elongated shape of the satellite buses. The a priori box-wing model may be a solution for twice-per-revolution empiri-cal parameters from the CODE’s new solar radiation pressure model are used for GNSS orbit determination [32]. In the paper, we adopt the one-step orbit determination of the GPS, GLONASS, Galileo and BDS in order to study that the parallel algorithm improves the efficiency of satellite orbit prediction, which does not affect the computational efficiency of the parallel algorithm. For BDS, no official metadata for a priori box-wing models are available. However, estimated parameters can be still used [33].
Equation (2) is the basic observation equation for ultra-rapid orbit determination. In order to determine the ultra-rapid orbit, the traditional sequential batch processing method traverses the observation data in all arcs by epoch to construct the normal equations. In this study, a parallel method based on MPI is used in this process. The main idea is to establish the multi-arc normal equation in the orbit determination by adopting the multi-processing method. Assuming that the single-day data is divided into three arcs, the length of each arc is 8 h. Each sub-process processes the measured data of the current arc and then calculates the observation equation coefficients. Based on the current time period, sub-normal equations including initial satellite state, force model parameters, earth rotation parameters, atmospheric parameters, station coordinates and ambiguity parameters are generated. Each process sends the sub-normal equation matrices and auxiliary parameters to the main process to obtain the total normal equations.
The parameter sequence of each period is consistent with that of the entire original period. The mathematical model and the parameters have not changed, so there is no additional conversion and connection. Supposing that the entire arc is divided into k sub-arcs, the least squares problem can be solved by accumulating the normal equations. The sub-normal equations and the entire normal equation can be expressed as:
{ N 1 δ x = W 1 N 2 δ x = W 2 N k δ x = W k
( N 1 + N 2 + + N k ) δ x = ( W 1 + W 2 + + W k )
where Nk and Wk are the related matrix of sub-normal equations and δ x is unknown parameter vector. Based on this entire normal equation, the elimination of parameters and least squares (LS) estimation can be then carried out. It is worth noting that the number of time periods are designed according to the computer hardware and time requirements.

2.2. Parallel Parameter Elimination Based on OpenMP

The parameters to be estimated in the GNSS observation equation of POD include station coordinate parameters, satellite orbit parameters, tropospheric parameters, ambiguity parameters, and earth rotation parameters. Usually, only some of these parameters are of concern to us, and the other parameters can be eliminated to improve the computational efficiency [41]. Equivalent elimination is a classic algorithm that can reduce the dimensionality of normal equations without a precision loss (Yao 2004). The error equation can be expressed as:
v = L ( A B ) [ x ^ 1 x ^ 2 ] , P ,
where L is the observation vector of dimension m, A and B are the design matrices of dimension m × ( n r ) and m × r , x ^ 1 and x ^ 2 are unknown parameters of dimension n r and r , v is the residual error matrix of dimension m , n is the number of unknowns, m is the number of observations, and P is the weight matrix of dimension m × m .
The parameters in the equation are decomposed into [ x ^ 1 x ^ 2 ] as needed, taking into account the correlation between the observations. The normal equations are:
[ N 11 N 12 N 21 N 22 ] [ x ^ 1 x ^ 2 ] = [ b 1 b 2 ]
where
[ N 11 N 12 N 21 N 22 ] = [ A T P A A T P B B T P A B T P B ]
b 1 = A T P L , b 2 = B T P L
If only x 1 is concerned, x 2 can be eliminated in order to reduce the dimension of the normal equation; [ E Z 0 E ] is defined, where E is the identity matrix, 0 is the zero matrix, Z = N 12 N 22 1 , and N 22 1 is the inverse matrix of N 22 . The left and right sides of (8) are multiplied by the matrix as follows:
[ E Z 0 E ] [ N 11 N 12 N 21 N 22 ] [ x ^ 1 x ^ 2 ] = [ E Z 0 E ] [ b 1 b 2 ]
[ N 1 0 N 21 N 22 ] [ x ^ 1 x ^ 2 ] = [ R 1 b 2 ]
where N 1 = N 11 N 12 N 22 1 N 21 , R 1 = b 1 N 12 N 22 1 b 2 . Then the normal equation after elimination is:
N 1 x ^ 1 = R 1
The solution of the normal equation is:
x ^ 1 = N 1 1 R 1
Q x ^ 1 x ^ 1 = ( N 11 N 12 N 22 1 N 21 ) 1 = N 11 1 + N 11 1 N 21 ( N 22 N 12 N 11 1 N 12 ) 1 N 12 N 11
Equation (12) is substituted into (5):
x ^ 2 = N 22 1 ( b 2 N 21 x ^ 1 )
Q x ^ 2 x ^ 2 = N 22 1 + N 22 1 N 21 ( N 11 N 12 N 22 1 N 21 ) 1 N 12 N 22 1
When the number of observation equations and parameter dimensions is very large, the parameter elimination is time-consuming. In order to further improve computational efficiency, parallel optimization using OpenMP technology for parameter elimination is adopted.
When the initial orbit parameters and the force model parameters are determined by using LS estimation, the orbit of the entire arc will be updated in the order of epoch and satellite by using orbit integration. As with the the above-mentioned, we use the parallel orbit integration method to update the final orbit sequence.

2.3. Parallel Orbit Determination Process

The parallel approach proposed for the ultra-rapid orbit determination of the double-differenced network solution mainly adopts the MPI/OpenMP hybrid parallel programming model. The optimized orbit determination process is shown in Figure 1. After reading observation files, navigation ephemeris files, and some auxiliary model correction files, data preprocessing includes gross error detection, basic model correction, forming double-differenced baseline, and floating solution. Then, the entire time period to be calculated is divided into n segments, and multiple processes are developed based on MPI to read and construct sub-normal equations in parallel. Each process read from the first epoch to the last epoch of the current period so that the consistency of the current arc parameter sequence can be maintained. The above multiple sub-normal equations are accumulated to obtain the total normal equation. Next, in order to further improve the efficiency of orbit determination, the OpenMP multithreading model is used for the parameter elimination process. Finally, after obtaining the precise dynamic parameters by the LS estimation, the same multi-thread model is used to realize the parallel orbit integration to obtain the final satellite POD results. The entire parallel optimization does not destroy the correlation between the parameters without the ambiguity segmentation problem. The input and output of orbit determination do not need to be adjusted greatly, which is easy to implement. It can be applied to single satellite POD as well as multi-GNSS POD.

3. Computation and Comparison

In order to verify the feasibility and effectiveness of the proposed GNSS parallel ultra-rapid orbit determination approach, multiple sets of experiments are carried out. The precision loss of the parallel POD, the improvement of computational efficiency, and the orbit accuracy after using the parallel approach are analyzed in detail.
Two sets of tracking stations from IGS/MGEX networks around the world are selected to analyze the accuracy and efficiency of multi-GNSS parallel ultra-rapid POD. The networks are shown in Figure 2. Relevant data from day of year (DOY) 182 to 189, 2018 are selected for the processing. The approach is implemented by Fortran and intel MPI/OpenMP parallel library. The experiments are carried out on a Linux computer cluster.

3.1. Impact of Parallel Approach on the Orbit Solution

Before further analyzing the efficiency and accuracy of the parallel approach, it is necessary to analyze the accuracy of the non-parallel GNSS POD method. The observation data from DOY 182 to 188, 2018 of the above 131 IGS/MGEX stations are used to carry out the multi-GNSS POD. Analysis of the errors in the radial, along and cross directions of the multi-GNSS orbits using the non-parallel method is performed. As a result, the along track direction accuracy is the lowest, followed by the cross track, and the radial is the highest. Except for slightly larger fluctuations in the BDS in the sense of the orbit RMS, the other systems have little change and the accuracy is relatively stable. The accuracy of each system is further analyzed, which is shown in Table 3. It can be seen from the table that the 1D mean RMS ( 1 D R M S = ( σ R 2 + σ A 2 + σ C 2 ) / 3 ) of GPS, GLONASS, Galileo and BDS are 1.3 cm, 2.7 cm, 2.7 cm, and 5.4 cm, respectively, which are at the same accuracy level as reported by international institutions [42]. The results show the correctness of the non-parallel POD processing, which lays the foundation for the proposed parallel POD.
It is necessary to further analyze the impact of the parallel approach on the covariance matrix and solution of the normal equations, as well as the precision loss of orbit. Based on MPI, the entire 480 epochs are divided into 2 periods. Statistics of the difference between the matrix elements of the parallel method and those of the non-parallel method is shown in Table 4. It can be seen from these figures that the error impact of the parallel approach on the covariance matrix and unknown parameters are approximately 10−13 and 10−6, respectively.
Using the same set of data, the normal equation coefficients of the above-mentioned parallel approach and the normal equation coefficients of the non-parallel approach are saved. After the procedures of parameter estimation, orbit update, and result output, the two sets of orbits are obtained. The comparison of the two sets of orbits show that their difference is below the millimeter level, which indicates that the precision loss caused by parallel method can be negligible.

3.2. Improvement of Computational Efficiency of Parallel Approach

In order to verify the characteristics of the parallel approach in terms of efficiency, the parallel approach is compared with the non-parallel method. The efficiency improvement of the approach in many aspects is further analyzed, such as multi-system, number of threads, and number of processes. In addition, the computational time changing with the number of stations is analyzed, which reflects the advantage of the parallel approach in the improvement of efficiency and provides a reference for the selection of the number of stations in the 1-h update frequency ultra-rapid POD.

3.2.1. Improvement of Efficiency with Different Number of Satellite Systems

The computational efficiency of GPS, GPS+GLONASS and GPS+GLONASS+Galileo+BDS parallel POD are compared to analyze the influence of the number of satellites. Two processes and ten threads are used in the parallel scheme. The statistical results of computational time are shown in Table 5, and the comparison between non-parallel and parallel approach is shown in Figure 3.
It can be seen from Table 5 that the more satellite systems that are involved, the longer it takes for POD. After using MPI parallelism, the calculation time of this process is significantly reduced, indicating that it is reasonable to optimize this process in parallel.
It can be seen from Figure 3 that the efficiency of GPS, GPS+GLONASS and GPS+GLONASS+Galileo+BDS POD have been obviously improved by the proposed parallel approach. As the number of satellites increases, the efficiency of the parallel POD has been improved more significantly, and the efficiency of the four-system combined POD has been improved the most. The parameters to be estimated and the dimension of the equation both increase nonlinearly, so the advantages of the parallel POD are more obvious.

3.2.2. Improvement of Efficiency with Different Number of Threads

The computational efficiency of the parallel parameter elimination process with a different number of threads based on OpenMP is further analyzed. According to the hardware configuration of the platform, the number of threads is sequentially increased from 2 to 16 at an interval of 2, and the computational time for the parameter elimination process of each parallel scheme is compared, which is shown in Figure 4.
It can be seen from Figure 4 that as the number of threads increases, the computational time of parameter elimination in GNSS POD gradually decreases and the scheme with 16 threads has the shortest computational time. However, the speed of efficiency improvement gradually decreases, most likely because when the number of threads is close to or greater than the number of cores, queuing round calculations occur, and the opening and ending of threads take more time. Therefore, further increasing the number of threads beyond the hardware core is of little significance to the improvement of parallel efficiency.

3.2.3. Improvement of Efficiency with Different Number of Processes

In order to further reflect the improvement of efficiency by the parallel POD approach, the influence of the number of processes on efficiency is analyzed. Multi-process is mainly used to construct normal equations in parallel. All schemes use the same number of threads, the same processing strategy, and select the number of processes as the only control variable to analyze the effect of increasing the number of processes on the computational time. The time consumption is shown in Table 6, and the speedup is shown in Figure 5.
As can be seen from Figure 5 and Table 6, compared to a single-process, the speedup ratios of two-process, three-process, four-process, and five-process are 1.788, 2.525, 3.102, and 3.406, respectively. As the number of processes increases, the speedup ratio increases. When five processes are used, the speedup in POD reaches the maximum. The speedup shows a nonlinear trend probably because the cost of opening and closing processes and inter-process communication increases as the number of processes increases, and depends on the proportion of parallel programs to the total program. In general, if the hardware can support more processes, the computational efficiency of POD can be further improved, which is similar to multithreading.

3.2.4. Improvement of Efficiency with Different Number of Stations

Observation data from 40 to 110 IGS and MGEX stations are selected to analyze the improvement of efficiency of parallel multi-GNSS POD. The results are shown in Figure 6.
It can be seen from Figure 6 that the parallel POD is obviously more efficient than the non-parallel POD. As the number of IGS/MGEX stations increases, the serial calculation time of non-parallel and parallel POD both gradually increases. The speedup ratio of parallel POD increases with the increase in the number of stations. The efficiency of POD using 80 to 90 stations has been significantly improved. Moreover, the parallel POD of four satellite system can be controlled within one hour, which provides the possibility to increase the update frequency to 1-h for multi-GNSS parallel POD.

3.3. Analysis of the Accuracy of Parallel Ultra-Rapid POD

Aiming at the problem that the predicted orbit accuracy of GNSS ultra-rapid products becomes worse as the extrapolation time is longer, we apply the parallel approach to the multi-GNSS ultra-rapid POD. Based on the results of the improvement of efficiency with different number of stations in the previous section, 88 IGS/MGEX stations are selected for multi-GNSS ultra-rapid POD. The distribution of stations is shown in Figure 2.

3.3.1. Improvement of Accuracy with Increased Update Frequency by the Parallel POD

The update frequency of the original ultra-rapid orbit is 6 h, and the update frequency can be increased to 1 h after applying the parallel POD. The accuracy of orbit with 6-h update frequency and 1-h update frequency are compared in this section. Observation data from DOY 182 to 189, 2018 are selected for multi-GNSS ultra-rapid POD. The accuracy of GPS, GLONASS, Galileo, and BDS predicted orbit with parallel and non-parallel POD are analyzed, which are shown in Figure 7, Figure 8, Figure 9 and Figure 10.
From Figure 7, Figure 8, Figure 9 and Figure 10, it can be seen that the accuracy of the predicted orbit of the 6-h update frequency fluctuates at different extrapolation times. The accuracies of orbit at 0 o’clock, 6 o’clock, 12 o’clock and 18 o’clock are the highest, and then gradually decrease. After the update frequency is increased to 1-h by parallel approach, the accuracy of predicted orbit is higher and more stable at different extrapolation time. The statistical results of the accuracy of GPS/GLONASS/Galileo/BDS orbit are shown in Table 7. Selecting the precise orbit of IGS and CODE analysis center as reference, the along, normal, radial and 1D mean RMS values of predicted orbit of 1-h update frequency and 6-h update frequency are compared.
It can be seen from the table that the accuracy of GPS, GLONASS, Galileo, and BDS 1-h updated orbits are improved in normal, radial and along direction, and the radial direction is improved the most compared to the 6-h updated orbit. The RMS of the predicted orbit of GPS, GLONASS, Galileo and BDS are 33.4%, 31.4%, 40.1% and 32.8% higher than those of the original 6-h results, respectively.

3.3.2. Evaluation of Predicted Orbit Accuracy of Ultra-Rapid POD by Parallel Approach

In order to further analyze the performance of the parallel approach, the accuracy of the ultra-rapid predicted orbit is evaluated. The GPS/GLONASS precise orbit of IGS and the Galileo/BDS precise orbit of CODE are selected as reference values. The experiment includes three groups. The first one uses the proposed parallel POD approach with a 1-h updated result (denoted as SDU-1h). The second one uses the non-parallel POD approach with a 6-h updated result (denoted as SDU-6h), and the third one is the 1-h updated result of Wuhan University analysis center (denoted as WHU-1h). The statistical results of accuracy are shown in Figure 11.
From Figure 11, it can be seen that as the extrapolation time increases, the orbit accuracy of SDU-6h decreases in the orbit sequence. The phenomenon that the orbit accuracy fluctuates with extrapolation time does not exist in WUH-1h and SDU-1h. The accuracy of WUH-1h is better than that of SDU-6h result, and it is relatively stable except for occasionally large outliers. The predicted orbits of GPS, GLONASS, Galileo, and BDS of SDU-1h have the highest accuracy and the most stable performance. Combined with the statistical results in Table 8, the 1D mean RMS of the ultra-rapid predicted orbit with the proposed parallel approach is 3.2 cm, 5.1 cm, 5.6 cm and 11.8 cm, respectively.
The accuracy of the predicted orbit of each satellite from GPS/GLONASS/Galileo/BDS are plotted in Figure 12, and BDS-2 MEO and BDS-2 IGSO of the BDS system are used in this processing. It can be seen from these figures that the predicted orbits of GPS, GLONASS, Galileo and BDS determined by the parallel approach with 1-h update frequency have the highest accuracy. The accuracies of GPS, GLONASS, Galileo and BDS predicted orbit are mostly better than 3 cm, 5 cm, 6 cm, and 12 cm, respectively. One of the GPS satellites, G18, has poor accuracy, which may be due to its large orbit error during the shadow period. In general, the accuracy of the ultra-rapid predicted orbit of each GNSS satellite with the parallel approach and 1-h updated frequency have greatly improved.

4. Conclusions

A parallel ultra-rapid POD approach has been proposed in this study. The parallel approach properly makes use of parallel computing to improve the efficiency of multi-GNSS orbit determination. The construction of normal equations, parameter elimination and orbit extrapolation in the GNSS POD algorithm are optimized based on MPI/OpenMP. The following conclusions can be drawn:
(1) The precision loss and feasibility of parallel POD are verified. The results show that 1D mean RMS of GPS, GLONASS, Galileo and BDS orbit in non-parallel mode are 1.3 cm, 2.7 cm, 2.7 cm, and 5.4 cm, respectively. Moreover, the precision loss of using the parallel approach is at the sub-millimeter level, which can be ignored, thus verifying the feasibility of the parallel method.
(2) The improvement of efficiency by the parallel approach is analyzed in detail, including multiple factors such as the number of processes, the number of threads, the number of stations and the number of satellite systems. The results show that with the increase of the number of threads and processes, the efficiency of the parallel approach increases more obviously. The approach is more applicable to large reference networks and multiple satellite systems.
(3) The superiority of the parallel approach used in multi-GNSS ultra-rapid POD is verified. The improvement of accuracy of predicted orbit with an increased update frequency of 1-h is analyzed in detail. The 1-h updated orbit determined by the parallel approach is compared with the 6-h updated orbit determined by the non-parallel approach, and the results show that the accuracy of GPS, GLONASS, Galileo and BDS are 33.4%, 31.4%, 40.1% and 32.8% higher than the 6-h results of the non-parallel approach, respectively. Additionally, the accuracy of the 1-h updated orbit is evaluated by comparing with IGS and CODE precise orbit products. The 1-h orbit product of Wuhan University is selected for comparison. The results show that the orbit with 1-h update frequency by the parallel approach has higher accuracy and stability. Further, statistical results show that the 1D mean RMS of the ultra-rapid predicted orbit of GPS, GLONASS, Galileo, and BDS with the proposed parallel algorithm is 3.2 cm, 5.1 cm, 5.6 cm and 11.8 cm, respectively.
In summary, the proposed parallel approach can greatly improve the efficiency of multi-GNSS ultra-rapid orbit determination and enhance the update frequency from 6-h to 1-h, which can help for the improvement of the ultra-rapid orbit accuracy and stability. The proposed parallel approach makes it possible to use a denser network of reference stations for multi-GNSS ultra-rapid POD with high update frequency for real-time users.
It should be noted that we adopted the ECOM1 SRP model for the one-step orbit determination of the GPS, GLONASS, Galileo and BDS in this study. For the specific elongated shape of Galileo satellites, orbit accuracy will be improved by using a prior model. In addition, clock determination is not related in this study, making it difficult to get the Signal In Space Range Error (SISRE). In the next studies, we will use an a priori model to reduce the orbit perturbations for Galileo, and we will combine the orbit and clocks through the processing in order to calculate clocks and orbital SISRE for the analysis of the algorithm.

Author Contributions

All authors contributed to the design of the study. C.J. and T.X. came up with the idea of the study. C.J. carried out the experiments and performed the data analysis. C.J., W.N., Z.F. and S.W. discussed the results. C.J. wrote the manuscript with support of T.X., W.N., Z.F., S.W. and A.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was founded by the National Key Research Program and Development of China (2020YFB0505800, 2016YFB0501900, and 2016YFC0803102), the National Natural Science Foundation of China (41874032,42074012 and 42030109), the Liaoning Key Research and Development Program (2020JH2/10100044).

Data Availability Statement

The GNSS observation data can be downloaded from https://cddis.nasa.gov/archive/gnss/data/daily/, accessed January 2019. The orbit products can be downloaded from https://cddis.nasa.gov/archive/gps/products/, accessed January 2019, ftp://igs.gnsswhu.cn/pub/gnss/products/mgex/, accessed January 2019 and http://ftp.aiub.unibe.ch/CODE_MGEX/, accessed January 2019. Questions and requests for materials can be addressed to [email protected] or [email protected].

Acknowledgments

The authors would like to extend their sincere gratitude to IGS/MGEX for providing the relevant data. The numerical calculation in this study have been done on the supercomputing system in the Supercomputing Center, Shandong University, Weihai.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, C.; Li, X.; Ge, M.; Heinkelmann, R.; Nilsson, T.; Soja, B.; Schuh, H. Estimation and evaluation of real-time precipitable water vapor from GLONASS and GPS. GPS Solut. 2016, 20, 703–713. [Google Scholar] [CrossRef] [Green Version]
  2. Ding, W.; Teferle, F.N.; Kazmierski, K.; Laurichesse, D.; Yuan, Y. An evaluation of realtime troposphere estimation based on GNSS Precise Point Positioning. J. Geophys. Res. Atmos. 2017, 122, 2779–2790. [Google Scholar] [CrossRef]
  3. Ge, M.; Chen, J.; Gendt, G. EPOS-RT: Software for real-time GNSS data processing. Geophys. Res. Abstr. 2009, 11, EGU2009-8933. [Google Scholar]
  4. Zhang, B.; Chen, Y.; Yuan, Y. PPP-RTK based on undifferenced and uncombined observations: Theoretical and practical aspects. J. Geod. 2018, 93, 1011–1024. [Google Scholar] [CrossRef]
  5. Zhang, B.; Teunissen, P.J.G.; Odijk, D. A Novel Un-differenced PPP-RTK Concept. J. Navig. 2011, 64, S180–S191. [Google Scholar] [CrossRef] [Green Version]
  6. 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]
  7. Hauschild, A.; Tegedor, J.; Montenbruck, O.; Visser, H.; Markgraf, M. Precise onboard orbit determination for LEO satellites with real-time orbit and clock corrections. In Proceedings of the ION GNSS 2016, Portland, OR, USA, 12–16 September 2016; pp. 3715–3723. [Google Scholar]
  8. Hein, G.W. Status, perspectives and trends of satellite navigation. Satell. Navig. 2020, 1, 22. [Google Scholar] [CrossRef]
  9. Cui, Y.; Chen, Z.; Li, L.; Zhang, Q.; Luo, S.; Lu, Z. An efficient parallel computing strategy for the processing of large GNSS network datasets. GPS Solut. 2021, 25, 36. [Google Scholar] [CrossRef]
  10. Montenbruck, O.; Steigenberger, P.; Prange, L.; Deng, Z.; Zhao, Q.; Perosanz, F.; Romero, I.; Noll, C.; Stürze, A.; Weber, G.; et al. The Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS)–achievements, prospects and challenges. Adv. Space Res. 2017, 59, 1671–1697. [Google Scholar] [CrossRef]
  11. Kouba, J.; Héroux, P. Precise Point Positioning Using IGS Orbit and Clock Products. GPS Solut. 2001, 5, 12–28. [Google Scholar] [CrossRef]
  12. Guo, F.; Li, X.; Zhang, X.; Wang, J. Assessment of precise orbit and clock products for Galileo, BeiDou, and QZSS from IGS Multi-GNSS Experiment (MGEX). GPS Solut. 2016, 21, 279–290. [Google Scholar] [CrossRef]
  13. Sośnica, K.; Zajdel, R.; Bury, G.; Bosy, J. Quality assessment of experimental IGS multi-GNSS combined orbits. GPS Solut. 2020, 24, 54. [Google Scholar] [CrossRef] [Green Version]
  14. Lin, C.S.; Cole, T.R. Dynamic model for global positioning system block IIR space vehicle. J. Spacecr. Rocket. 1997, 34, 354–359. [Google Scholar] [CrossRef]
  15. Han, S.C.; Kwon, J.H.; Jekeli, C. Accurate absolute GPS positioning through satellite clock error estimation. J. Geod. 2001, 75, 33–43. [Google Scholar] [CrossRef]
  16. IGS Central Bureau. International GNSS Service Progress Report 2008–2012; Jet Propulsion Laboratory California Institute of Technology: Pasadena, CA, USA, 2013. [Google Scholar]
  17. Dach, R.; Lutz, S.; Walser, P.; Fridez, P. Bernese GNSS Software Version 5.2; Astronomical Institute, University of Bern: Bern, Switzerland, 2015. [Google Scholar]
  18. Zeng, T.; Sui, L.; Ruan, R.; Jia, X.; Feng, L. Uncombined precise orbit and clock determination of GPS and BDS-3. Satell. Navigation. 2020, 1, 19. [Google Scholar] [CrossRef]
  19. Ge, M. Study on Theory and Software of Precise Orbit Determination for GPS Satellite. Ph.D. Thesis, Wuhan University, Wuhan, China, 1995. (In Chinese). [Google Scholar]
  20. Li, X.; Xu, Y.; Wang, L. Undifferenced Precise Satellite Clock Error Estimation and Precision Analysis. Geomat. Inf. Sci. Wuhan Univ. 2010, 35, 661–664. [Google Scholar]
  21. Steigenberger, P.; Hugentobler, U.; Loyer, S.; Perosanz, F.; Montenbruck, O. Galileo orbit and clock quality of the IGS Multi-GNSS Experiment. Adv. Space Res. 2014, 55, 269–281. [Google Scholar] [CrossRef]
  22. Guo, J.; Xu, X.; Zhao, Q.; Liu, J. Precise orbit determination for quad-constellation satellites at Wuhan University: Strategy, result validation, and comparison. J. Geod. 2015, 90, 1–7. [Google Scholar] [CrossRef]
  23. Liu, Y.; Ge, M.; Shi, C.; Lou, Y.; Wickert, J.; Schuh, H. Improving integer ambiguity resolution for GLONASS precise orbit determination. J. Geod. 2016, 90, 715–726. [Google Scholar] [CrossRef]
  24. Li, J.; Zeng, G.; Zhu, J.; Yang, J. precise orbit determination and clock error estimation for multi-GNSS[C]. In China Satellite Navigation Conference; Springer: Singapore, 2017; pp. 50–56. [Google Scholar]
  25. Springer, T.A.; Hugentobler, U. IGS ultra rapid products for (near-) real-time applications. Phys. Chem. Earth Part A Solid Earth Geod. 2001, 26, 623–628. [Google Scholar] [CrossRef]
  26. Beutler, G.; Rothacher, M.; Schaer, S.; Springer, T.A.; Kouba, J.; Neilan, R.E. The International GPS Service (IGS): An interdisciplinary service in support of Earth sciences. Adv. Space Res. 1999, 23, 631–653. [Google Scholar] [CrossRef]
  27. Yao, Y. Theory and realization of GPS orbit integration. Geospat. Inf. Sci. 2008, 11, 1–5. [Google Scholar] [CrossRef] [Green Version]
  28. Lou, Y.; Shi, C.; Ge, M.; Zhao, Q. GPS Real Time Orbit Determination and Initial Results Analysis. Geomat. Inf. Sci. Wuhan Univ. 2008, 33, 815–817. [Google Scholar]
  29. Steigenberger, P.; Hugentobler, U.; Montenbruck, O.; Hauschild, A. Precise orbit determination of GIOVE-B based on the CONGO network. J. Geod. 2011, 85, 357–365. [Google Scholar] [CrossRef]
  30. Zhao, Q.; Xu, X.; Ma, H.; Liu, J. Real-Time Precise Orbit Determination of BDS/GNSS:Method and Service. Geomat. Inf. Sci. Wuhan Univ. 2018, 43, 2157–2166. [Google Scholar]
  31. Li, X.; Chen, X.; Ge, M.; Schuh, H. Improving multi-GNSS ultra-rapid orbit determination for real-time precise point positioning. J. Geod. 2019, 93, 45–64. [Google Scholar] [CrossRef]
  32. Arnold, D.; Meindl, M.; Beutler, G.; Dach, R.; Schaer, S.; Lutz, S.; Prange, L.; Sośnica, K.; Mervart, L.; Jäggi, A. CODE’s new solar radiation pressure model for GNSS orbit determination. J. Geod. 2015, 89, 775–791. [Google Scholar] [CrossRef] [Green Version]
  33. Duan, B.; Hugentobler, U.; Selmke, I.; Marz, S.; Killian, M.; Rott, M. BeiDou satellite radiation force models for precise orbit determination and geodetic applications. TechRxiv 2021. [Google Scholar] [CrossRef]
  34. Böhm, J.; Niell, A.; Tregoning, P.; Schuh, H. Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data. Geophys. Res. Lett. 2006, 33, L07304. [Google Scholar] [CrossRef] [Green Version]
  35. Chen, G.; Herring, T.A. Effects of atmospheric azimuthal asymmetry on the analysis of space geodetic data. J. Geophys. Res. 1997, 102, 20489–20502. [Google Scholar] [CrossRef]
  36. Springer, T.A.; Beutler, G.; Rothacher, M. A new solar radiation pressure model for the GPS satellites. Adv. Space Res. 1999, 23, 673–676. [Google Scholar] [CrossRef]
  37. He, L.; Ge, M.; Wang, J.; Wickert, J.; Schuh, H. Experimental Study on the Precise Orbit Determination of the BeiDou Navigation Satellite System. Sensors 2013, 13, 2911–2928. [Google Scholar] [CrossRef] [Green Version]
  38. Pavlis, N.; Holmes, S.; Kenyon, S.; Factor, J. The development and evaluation of the earth gravitational model 2008 (EGM2008). J. Geophys. Res. 2012, 117, B04406. [Google Scholar] [CrossRef] [Green Version]
  39. Petit, G.; Luzum, B. IERS Conventions. IERS Conventions Centre. France. 2010. Available online: https://iers-conventions.obspm.fr/content/tn36.pdf (accessed on 1 September 2010).
  40. Lyard, F.; Lefevre, F.; Letellier, T.; Francis, O. Modelling the global ocean tides: Modern insights from FES2004. Ocean Dyn. 2006, 56, 394–415. [Google Scholar] [CrossRef]
  41. Xu, G.; Xu, Y. GPS: Theory, Algorithms and Applications; Springer Publishing Company: Potsdam, Germany, 2016. [Google Scholar]
  42. Kazmierski, K.; Sośnica, K.; Hadas, T. Quality assessment of multi-GNSS orbits and clocks for real-time precise point positioning. GPS Solut. 2018, 22, 1–12. [Google Scholar] [CrossRef] [Green Version]
Figure 1. GNSS parallel precision orbit determination flowchart.
Figure 1. GNSS parallel precision orbit determination flowchart.
Remotesensing 13 03464 g001
Figure 2. Distribution of GNSS-tracking stations. Red triangles denote stations used for GNSS precise orbit determination and blue pentagrams denote stations used for GNSS ultra-rapid orbit determination.
Figure 2. Distribution of GNSS-tracking stations. Red triangles denote stations used for GNSS precise orbit determination and blue pentagrams denote stations used for GNSS ultra-rapid orbit determination.
Remotesensing 13 03464 g002
Figure 3. Computational time of parallel and non-parallel orbit determination.
Figure 3. Computational time of parallel and non-parallel orbit determination.
Remotesensing 13 03464 g003
Figure 4. Computational time for different number of threads.
Figure 4. Computational time for different number of threads.
Remotesensing 13 03464 g004
Figure 5. Speedup of satellite orbit determination with multi-process parallel computing.
Figure 5. Speedup of satellite orbit determination with multi-process parallel computing.
Remotesensing 13 03464 g005
Figure 6. Serial and parallel computational time for satellite orbit determination with different number of stations.
Figure 6. Serial and parallel computational time for satellite orbit determination with different number of stations.
Remotesensing 13 03464 g006
Figure 7. Accuracy of GPS ultra-rapid predicted orbit updated per 6 h and GPS ultra-rapid predicted orbit updated per 1 h.
Figure 7. Accuracy of GPS ultra-rapid predicted orbit updated per 6 h and GPS ultra-rapid predicted orbit updated per 1 h.
Remotesensing 13 03464 g007
Figure 8. Accuracy of GLONASS ultra-rapid predicted orbit updated per 6 h and GLONASS ultra-rapid predicted orbit updated per 1 h.
Figure 8. Accuracy of GLONASS ultra-rapid predicted orbit updated per 6 h and GLONASS ultra-rapid predicted orbit updated per 1 h.
Remotesensing 13 03464 g008
Figure 9. Accuracy of Galileo ultra-rapid predicted orbit updated per 6 h and Galileo ultra-rapid predicted orbit updated per 1 h.
Figure 9. Accuracy of Galileo ultra-rapid predicted orbit updated per 6 h and Galileo ultra-rapid predicted orbit updated per 1 h.
Remotesensing 13 03464 g009
Figure 10. Accuracy of BDS ultra-rapid predicted orbit updated per 6 h and BDS ultra-rapid predicted orbit updated per 1 h.
Figure 10. Accuracy of BDS ultra-rapid predicted orbit updated per 6 h and BDS ultra-rapid predicted orbit updated per 1 h.
Remotesensing 13 03464 g010
Figure 11. RMS of GPS, GLONASS, Galileo and BDS (from top to bottom) orbit of SDU-1h, SDU-6h and WHU-1h compared with CODE precise orbit.
Figure 11. RMS of GPS, GLONASS, Galileo and BDS (from top to bottom) orbit of SDU-1h, SDU-6h and WHU-1h compared with CODE precise orbit.
Remotesensing 13 03464 g011
Figure 12. RMS of predicted orbit of each satellite from GPS, GLONASS, Galileo and BDS (from top to bottom) using SDU-1h, SDU-6h and WUH-1h strategy.
Figure 12. RMS of predicted orbit of each satellite from GPS, GLONASS, Galileo and BDS (from top to bottom) using SDU-1h, SDU-6h and WUH-1h strategy.
Remotesensing 13 03464 g012
Table 1. Accuracy and latency of each IGS orbit products.
Table 1. Accuracy and latency of each IGS orbit products.
Type AccuracyLatencyUpdatesSample Interval
Broadcastorbits~100 cmreal-time-daily
Ultra-rapidorbits~5 cmreal-time4 times 15 min
Ultra-rapidorbits~3 cm3–9 h4 times15 min
Rapidorbits~2.5 cm17–41 hdaily15 min
Finalorbits~2.5 cm12–18 daysweekly15 min
Table 2. Multi-GNSS ultra-rapid orbit determination strategies.
Table 2. Multi-GNSS ultra-rapid orbit determination strategies.
ItemDescriptionRemark
ObservationsDouble-differenced phase and pseudo-range measurements External a
Satellite systemGPS, GLONASS, BDS, GalileoExternal
Sampling rate 30 sExternal
Elevation cut-off angleExternal
Ionosphere modelIonosphere-free combination with higher-order-ionosphere correctionExternal
Tropospheric modelGMF [34] with 2-h resolution of station-specific troposphere parameters, gradients [35] with daily resolutionEstimated b
PCO/PCVIGS14.atxExternal
A priori orbitMGEX broadcast ephemerisExternal
Stations coordinatesIGS SINEX productsExternal
Orbital model3 position and 3 velocity coordinates at the start of the orbital arc, without pseudo-stochastic pulsesEstimated
Solar Radiation Parameter (SRP) modelECOM1 model with five parameters [36]External
ERPxp, yp, xp-rate, yp-rate and LOD are estimated with tight constraint of UT1-UTC [37]Estimated
Gravity field modelEGM2008, 12 × 12 degree [38]External
N-body perturbationJPL DE405 ephemerisExternal
Solid earth tides, pole tides/Semidiurnal tide IERS Conventions 2010 [39]External
Ocean Tidal modelFES2004 [40]External
Precession and nutation modelsIAU2006/IAU2000AExternal
Relativistic effectsIERS Conventions 2010External
a The remark of External means the data source from models or products provided by other agencies; b The remark of Estimated means the data source generated by our software.
Table 3. Accuracy statistics of multi-GNSS orbit using the non-parallel method (cm).
Table 3. Accuracy statistics of multi-GNSS orbit using the non-parallel method (cm).
Satellite SystemRACRMS
GPS1.01.51.21.3
GLONASS1.63.22.62.7
Galileo2.72.82.42.7
BDS5.95.34.25.4
Table 4. Statistics of the difference between the matrix elements of the parallel method and those of the non-parallel method.
Table 4. Statistics of the difference between the matrix elements of the parallel method and those of the non-parallel method.
MaxMinMeanSTD
Difference of ATPA1.019 × 10−9−1.019 × 10−9−9.633 × 10−181.365 × 10−13
Difference of ATPL1.023 × 10−12−1 × 10−11−6.233 × 10−161.22 × 10−13
Difference of unknown parameters2.913 × 10−5−2.525 × 10−57.519 × 10−83.274 × 10−6
Table 5. Time-consuming statistical results of GNSS parallel and non-parallel POD (unit: min).
Table 5. Time-consuming statistical results of GNSS parallel and non-parallel POD (unit: min).
Satellite SystemNon-Parallel PODInitial Orbit EstimationParallel Initial Orbit EstimationParallel POD
GPS2812623
GPS + GLONASS60462246
GNSS(G/R/E/C)17816046122
Table 6. Time consumption of multi-process parallel computing.
Table 6. Time consumption of multi-process parallel computing.
Number of Process12345
Computational Time824.0s460.8s326.3s265.6s241.9s
Table 7. Accuracy of ultra-rapid predicted Multi-GNSS orbit updated per 6 h.
Table 7. Accuracy of ultra-rapid predicted Multi-GNSS orbit updated per 6 h.
R (cm)A (cm)C (cm)1D RMS (cm)Improvement of Accuracy
GPS_62.36.82.74.8
GPS_12.33.92.23.233.4%
GLO_63.210.44.47.4
GLO_13.16.43.65.131.4%
Gal_66.312.44.49.3
Gal_15.85.23.65.640.1%
BDS_614.321.06.916.9
BDS_112.212.05.611.332.8%
Table 8. Statistics of the accuracy of ultra-rapid predicted orbit (cm).
Table 8. Statistics of the accuracy of ultra-rapid predicted orbit (cm).
Satellite SystemScheme 1Scheme 2Scheme 3
GPS3.24.84.3
GLONASS5.17.36.5
Galileo5.69.211.2
BDS11.816.916.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jiang, C.; Xu, T.; Nie, W.; Fang, Z.; Wang, S.; Xu, A. A Parallel Approach for Multi-GNSS Ultra-Rapid Orbit Determination. Remote Sens. 2021, 13, 3464. https://doi.org/10.3390/rs13173464

AMA Style

Jiang C, Xu T, Nie W, Fang Z, Wang S, Xu A. A Parallel Approach for Multi-GNSS Ultra-Rapid Orbit Determination. Remote Sensing. 2021; 13(17):3464. https://doi.org/10.3390/rs13173464

Chicago/Turabian Style

Jiang, Chunhua, Tianhe Xu, Wenfeng Nie, Zhenlong Fang, Shuaimin Wang, and Aigong Xu. 2021. "A Parallel Approach for Multi-GNSS Ultra-Rapid Orbit Determination" Remote Sensing 13, no. 17: 3464. https://doi.org/10.3390/rs13173464

APA Style

Jiang, C., Xu, T., Nie, W., Fang, Z., Wang, S., & Xu, A. (2021). A Parallel Approach for Multi-GNSS Ultra-Rapid Orbit Determination. Remote Sensing, 13(17), 3464. https://doi.org/10.3390/rs13173464

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