Research on Accelerating Single-Frequency Precise Point Positioning Convergence with Atmospheric Constraint

: An increasing number of researchers have conducted in-depth research on the advantages of low-cost single-frequency (SF) receivers, which can e ﬀ ectively use ionospheric information when compared to dual-frequency ionospheric-free combination. However, SF observations are bound to increase the unknown parameters and prolong the convergence time. It is desirable if the convergence time can be reduced by external information constraints, for example atmospheric constraints, which include ionosphere- or troposphere constraints. In this study, ionospheric delay constraints, tropospheric delay constraints, and their dual constraints were considered. Additionally, a total of 18,720 test experiments were performed. First, the nearest-neighbor extrapolation (NENE), bilinear-(BILI), bicubic- (BICU), and Junkins weighted-interpolation (JUNK) method of Global Ionospheric Map (GIM) grid products were analyzed. The statistically veriﬁed BILI in the percentage of convergence time, average convergence time, and computation time consumption of them shows a good advantage. Next, the inﬂuences of global troposphere- and ionosphere-constrained on the convergence time of SF Precise Point Positioning (PPP) were analyzed. It is veriﬁed that the ionosphere-constrained (TIC2) has signiﬁcant inﬂuence on the convergence time in the horizontal and vertical components, while the troposphere-constrained (TIC1) has better e ﬀ ect on the convergence time in the vertical components within some thresholds. Of course, the dual constraint (TIC3) has the shortest average convergence time, which is at least 46.5% shorter in static mode and 5.4% in kinematic mode than standard SF PPP (TIC0).


Introduction
The concept and technology of Precise Point Positioning (PPP) were first proposed and implemented for the Global Positioning System (GPS) by the American Jet Propulsion Laboratory (JPL) in the late 1990s [1]. PPP has attracted significant interest over the intervening years due to its high accuracy without needing a specific reference station, providing correctional information, simple operations, and cost effectiveness due to reductions in labor and equipment costs. Therefore, it has been extensively used in the areas, SF PPP cannot be combined or differenced to be eliminated or attenuated by a part of the error, like dual-frequency PPP. If there are some information related to the parameters, and this information is accurate enough, the precision of traditional SF PPP can be largely improved and the convergence time can be shortened. Zhang et al. [31] studied real-time GIM and its application in SF positioning, Aggrey and Bisnath [32] studied the effect of atmospheric-constrained on the convergence time of dual-and triple-frequency PPP, and Gao et al. [33] applied the Inertial Navigation System (INS) to the ionosphere-constrained PPP to overcome the drawbacks that accompany unexpected and unavoidable substandard observation environments.
This study uses the GIM products and the tropospheric zenith path delays from the IGS as the constrained information for SF PPP. The organization of this study is as follows: the next step details the methods for mathematical models of standard SF GPS PPP, troposphere-constrained SF GPS PPP, and ionosphere-constrained SF GPS PPP, and the details of the four interpolation methods of GIM products. Afterwards, in Section 3 we compare and analysis the influence of interpolation methods and constraint methods on SF GPS PPP convergence time, and finally in Section 4 we draw conclusions. Later, we will study ionospheric delay and tropospheric delay prediction models to provide virtual atmospheric delay observations for real-time PPP and also provide a priori information for the constraint processing.

Methods
The GPS PPP observation models are first derived. Afterwards, different interpolation methods for GIM products are introduced.

GPS PPP Observation Models
This study uses an undifferenced and uncombined GPS PPP model as compared with the traditionally used IF-PPP model. These models include a standard SF GPS PPP model, troposphere-constrained SF GPS PPP models, ionosphere-constrained SF GPS PPP models, and SF GPS PPP models with troposphere-and ionosphere-constrained.

Standard Single-frequency GPS PPP Observation Model
The distance from the satellite to the receiver can be measured while using pseudorange and carrier phase observations, with the following expression [34]: where, P s r and L s r are the original pseudorange and the carrier phase Φ s r is multiplied by wavelength λ 1 for the specific receiver r and satellite s; ρ s r is the geometric distance from the satellite to the receiver; c is the speed of light in vacuum; dt r and dt s are the clock error of the receiver and satellite; I s r is the slant ionospheric delay on the first frequency f 1 ; M w is the wet mapping function; Z w is the zenith wet delay; M h is the hydrostatic mapping function; Z h is the zenith hydrostatic delay; d r is the frequency-dependent receiver uncalibrated code delay (UCD) with respect to satellite s; d s is the frequency-dependent satellite UCD; N s r is the integer phase ambiguity; b r and b s are the uncalibrated phase delays (UPDs) for the receiver and satellite, which is frequency-dependent; and, ε s r and ξ s r are the sum of measurement noise caused by the pseudorange and carrier phase observations and the error caused by the multipath effect. Other errors have been modeled in advance.
The original pseudorange and carrier phase observations equations are linearized, and the receiver clock offset only absorbs common part of frequency dependent receiver UCDs. Supposing that m satellites are tracked simultaneously at a certain epoch by the receiver r, the standard SF PPP model can be written, as follows [35]: with where µ r is the unit vector of the coordinate component between the receiver and the satellite; x is the vector of the receiver position increments relative to the priori position; and, I t is a unit vector of 2 · m rows and one column, corresponding to the receiver clock parameter dt r . In matrix K, the elements for the corresponding p s r and l s r are 1 and −1, respectively, corresponding to the ionospheric parameter I s r . R is the matrix corresponding to the ambiguity parameters N s r , and the elements for the corresponding p s r and l s r are 0 and 1, respectively. Q L is the stochastic model of the observed minus computed (OMC) observables. IGS precise satellite clock was calculated by the IF combination observables and the IF combination of satellite UCDs was absorbed by the satellite clock offsets, as follows [2]:

Troposphere-Constrained Single-Frequency GPS PPP Observation Model
According to the tropospheric delay product published by IGS, a virtual observation is used in the observation equation, and a troposphere-constrained SF GPS PPP is then constructed. The observation equation is expressed, as follows: where I T is a unit vector of m rows and one column, corresponding to the zenith wet delay parameter Z w ; ∧ T r is derived from external tropospheric products, the products for each IGS station at 5 min intervals provided by the IGS Analysis Center; ε r,trop is the corresponding noise; Q T denotes the stochastic model of virtual troposphere observables; and, 0 is a zero matrix. The ionosphere-constrained SF GPS PPP can add the GIM product as a virtual observation to constrain the ionospheric parameters as compared with the standard SF GPS PPP. The constraint equation for the observation is as follows [35]: where ∧ I s r is derived from external ionospheric products with the corresponding noise ε r,ion ; I I is a unit vector of m rows and one column; and, Q I denotes the stochastic model of virtual ionospheric observables.

Troposphere-and Ionosphere-Constrained Single-Frequency GPS PPP Observation Model
If the external troposphere and ionospheric products are both used as virtual observations to constrain the observation equation, the following observation equation is obtained:

Different Interpolation Methods for GIM Products
It is necessary to interpolate the ionospheric grid point from GIM products to obtain the observation station VTEC before performing the ionosphere-constrained modeling. In this study, a comparative analysis of four interpolation methods is performed. These methods include nearest-neighbor extrapolation [36], bilinear interpolation [18], Junkins weighted interpolation, and bicubic interpolation. The principles of the Junkins weighted and bicubic interpolation methods are described below. Figure 1 shows the schematic of Junkins interpolation. In Figure 1, the latitude and longitude of the interpolation point are b and l, that is, p(b, l), and the values of the four GIM grid points adjacent to it are p 11 (b 1 , l 1 ), p 12 (b 1 , l 2 ), p 21 (b 2 , l 1 ), and p 22 (b 2 , l 2 ). It is assumed that the TEC values at p 11 , p 12 , p 21 , and p 22 are TEC 11 , TEC 12 , TEC 21 , and TEC 22 , respectively. The VTEC of point p can be interpolated according to the VTEC of the surrounding four grid points. The interpolation formula used is as follows:

Junkins Weighted Interpolation
The weighted function w n (x, y), n = 1, 2, 3, 4 is as follows: where We can use nearest-neighbor extrapolation when the p point only has two proximate GIM grid points.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 20 We can use nearest-neighbor extrapolation when the p point only has two proximate GIM grid points.

Bicubic Interpolation
The bicubic interpolation [37] method uses 16 adjacent points for interpolation ( Figure 2). Bicubic interpolation is an extension of cubic interpolation in a two-dimensional space. The cubic interpolation kernel is an approximation of the convolution interpolation of the ideal sampling sinc function at [−2,2] that is based on the cubic polynomial [38]. Typically, cubic interpolation produces interpolation coefficients that are based on a third-order polynomial: where q represents the distance between the interpolation point and the reference point, w is a tunable parameter, and the best result is obtained by verifying 0.5 w = - [38] with a large amount of data. Under these conditions, Equation (18) is simplified to:

Bicubic Interpolation
The bicubic interpolation [37] method uses 16 adjacent points for interpolation ( Figure 2). Bicubic interpolation is an extension of cubic interpolation in a two-dimensional space. The cubic interpolation kernel is an approximation of the convolution interpolation of the ideal sampling sinc function at [−2,2] that is based on the cubic polynomial [38]. Typically, cubic interpolation produces interpolation coefficients that are based on a third-order polynomial: where q represents the distance between the interpolation point and the reference point, w is a tunable parameter, and the best result is obtained by verifying w = −0.5 [38] with a large amount of data. Under these conditions, Equation (18) is simplified to: Bicubic interpolation is similar to bilinear interpolation. It is also decomposed into two one-dimensional interpolations: horizontal interpolation and vertical interpolation. Figure 2 shows the interpolation process. First, vertical interpolation is performed to obtain four virtual values, and then the values of the p point are obtained via interpolating based on the four virtual values. The calculation process is as follows: namely where  We can use bilinear interpolation when the number of GIM grid points around point p is less than 16.
Tropospheric delay interpolation is similar to inverse distance-weighted interpolation [19], and it is performed according to the relationship between the tropospheric delay value and time in the tropospheric products.

Experimental Data Set and Analysis
The PPP results based on raw observations are evaluated with 78 global IGS stations for the whole month of September 2017 to verify the improved troposphere-and ionosphere-constrained for GPS SF PPP performance. For every station, the 24-h observations were divided into eight 3-h sessions to evaluate the performance of PPP. Therefore, a total of 18720 experiments were tested and compared. Figure 3 shows the distribution of these stations in the world. Figure 4 show a flowchart of our study procedure, while Table 1 lists the models and strategies that were used in this study.
The PPP performance in terms of convergence time in the horizontal and vertical components is evaluated at a different confidence level (the 68% and 95%) in two modes (static and kinematic). The criterion for judging convergence is when the positioning error in the horizontal and vertical components is less than 0.5 m (95% level) and 0.3 m (68% level), respectively [11].  We can use bilinear interpolation when the number of GIM grid points around point p is less than 16.
Tropospheric delay interpolation is similar to inverse distance-weighted interpolation [19], and it is performed according to the relationship between the tropospheric delay value and time in the tropospheric products.

Experimental Data Set and Analysis
The PPP results based on raw observations are evaluated with 78 global IGS stations for the whole month of September 2017 to verify the improved troposphere-and ionosphere-constrained for GPS SF PPP performance. For every station, the 24-h observations were divided into eight 3-h sessions to evaluate the performance of PPP. Therefore, a total of 18720 experiments were tested and compared. Figure 3 shows the distribution of these stations in the world. Figure 4 show a flowchart of our study procedure, while Table 1 lists the models and strategies that were used in this study.
The PPP performance in terms of convergence time in the horizontal and vertical components is evaluated at a different confidence level (the 68% and 95%) in two modes (static and kinematic). The criterion for judging convergence is when the positioning error in the horizontal and vertical components is less than 0.5 m (95% level) and 0.3 m (68% level), respectively [11]. We can use bilinear interpolation when the number of GIM grid points around point p is less than 16.
Tropospheric delay interpolation is similar to inverse distance-weighted interpolation [19], and it is performed according to the relationship between the tropospheric delay value and time in the tropospheric products.

Experimental Data Set and Analysis
The PPP results based on raw observations are evaluated with 78 global IGS stations for the whole month of September 2017 to verify the improved troposphere-and ionosphere-constrained for GPS SF PPP performance. For every station, the 24-h observations were divided into eight 3-h sessions to evaluate the performance of PPP. Therefore, a total of 18720 experiments were tested and compared. Figure 3 shows the distribution of these stations in the world. Figure 4 show a flowchart of our study procedure, while Table 1 lists the models and strategies that were used in this study.
The PPP performance in terms of convergence time in the horizontal and vertical components is evaluated at a different confidence level (the 68% and 95%) in two modes (static and kinematic). The criterion for judging convergence is when the positioning error in the horizontal and vertical components is less than 0.5 m (95% level) and 0.3 m (68% level), respectively [11].

The Effect of Interpolation Method on Ionosphere-constrained GPS Single-frequency PPP
For the convenience of description, Table 2 shows the abbreviations, according to the four interpolation methods that were introduced above. For GIM grid products, this study discuss the  For the convenience of description, Table 2 shows the abbreviations, according to the four interpolation methods that were introduced above. For GIM grid products, this study discuss the convergence time of the four interpolation methods under ionosphere-constrained and the average time of calculation of 78 stations.  Figures 5 and 6 show the percentage of convergence time. From these two pictures, whether in the 95% or 68% level, and whether it is static mode or kinematic mode, the convergence time is less than 10 min, the percentage of the BILI method calculation result is slightly lower than the other three methods, and in the subsequent convergence time less than 20 min, 30 min, etc., the percentage of BILI method calculation results has been higher than the other three methods, as shown in the green histogram. That is to say, the convergence time of the BILI interpolation method is the shortest when compared to the other three methods. Figure 7 shows the statistical results of the root mean square (RMS) error in the horizontal and vertical components in the static and kinematic mode, respectively. The RMS errors are based on the statistics of the last 15 min of the position solution error [44]. In the static mode, the percentage of BILI method calculations exceeds 40% in the horizontal and vertical components when the RMS is less than 5 cm, and the other three methods are lower than 40% and, in the subsequent statistical stages, the percentage of the BILI method is higher than the other three methods to reach 100%. In the kinematic mode, the percentage of BILI method calculation results in the horizontal component is 60% when the RMS is less than 15 cm, and the other three methods are lower than this value; while, in the vertical component, the BILI method exceeds 40%, and the other three methods not reaching 40%, and the percentage of BILI methods in other statistical stages is higher than the other three methods. In summary, the BILI method has higher final positioning accuracy than the other three methods in both the static and kinematic modes in the horizontal or vertical component.   6 show the percentage of convergence time. From these two pictures, whether in the 95% or 68% level, and whether it is static mode or kinematic mode, the convergence time is less than 10 min., the percentage of the BILI method calculation result is slightly lower than the other three methods, and in the subsequent convergence time less than 20 min., 30 min., etc., the percentage of BILI method calculation results has been higher than the other three methods, as shown in the green histogram. That is to say, the convergence time of the BILI interpolation method is the shortest when compared to the other three methods. Figure 7 shows the statistical results of the root mean square (RMS) error in the horizontal and vertical components in the static and kinematic mode, respectively. The RMS errors are based on the statistics of the last 15 min. of the position solution error [44]. In the static mode, the percentage of BILI method calculations exceeds 40% in the horizontal and vertical components when the RMS is less than 5 cm, and the other three methods are lower than 40% and, in the subsequent statistical stages, the percentage of the BILI method is higher than the other three methods to reach 100%. In the kinematic mode, the percentage of BILI method calculation results in the horizontal component is 60% when the RMS is less than 15 cm, and the other three methods are lower than this value; while, in the vertical component, the BILI method exceeds 40%, and the other three methods not reaching 40%, and the percentage of BILI methods in other statistical stages is higher than the other three methods. In summary, the BILI method has higher final positioning accuracy than the other three methods in both the static and kinematic modes in the horizontal or vertical component.   The average convergence time (removed the data that the convergence time is equal to 180 min., that is, the data that does not convergence within three hours) is shown in Table 3, as compared to the JUNK method with the longest average convergence time, the average convergence time of BILI method at the 95% level is improved from 3.87 to 3.14 min. (18.9%) in horizontal (H(95%)) and by from 19.44 to 15.18 min. (21.9%) in the vertical components (V(95%)) in the static mode, while, in the  The average convergence time (removed the data that the convergence time is equal to 180 min., that is, the data that does not convergence within three hours) is shown in Table 3, as compared to the JUNK method with the longest average convergence time, the average convergence time of BILI method at the 95% level is improved from 3.87 to 3.14 min.  The average convergence time (removed the data that the convergence time is equal to 180 min, that is, the data that does not convergence within three hours) is shown in Table 3, as compared to the JUNK method with the longest average convergence time, the average convergence time of BILI method at the 95% level is improved from 3.87 to 3.14 min (18.9%) in horizontal (H(95%)) and by from 19.44 to 15.18 min (21.9%) in the vertical components (V(95%)) in the static mode, while, in the kinematic mode, it increased by 15.3% in the horizontal component and 13.2% in the vertical component. At the 68% level, the average convergence time in the horizontal component (H(68%)) significantly reduced by 17.9% from 17.55 to 14.40 min and by 23.6% from 42.92 to 32.81 min in the vertical component (V(68%)) in the static mode, and in the kinematic mode, lessened by 6.8% and 2.3% in the horizontal component and vertical component, respectively. The calculation time of these 18720 tests was normalized to FLRS station (Table 4). We used a Lenovo computer, which is configured with Windows 8.1 professional computer with an i5-3470 Central Processing Unit (CPU) and 4GB of installed Random Access Memory (RAM). From Table 4, it can be concluded that the BILI interpolation method calculation time is almost the same as the Junkins weighted interpolation method, which takes the shortest time. According to all of the above analysis, the interpolation method that is adopted for GIM products in the following ionosphere-constrained model is bilinear interpolation.

Single-Frequency PPP with the Constraints
In this section, we analyse the effects of single-frequency GPS on unconstrained, troposphereconstrained, ionosphere-constrained, and their dual constrains on convergence time. Table 5 summarizes the four constraint schemes.

TIC0
Standard single-frequency PPP model TIC1 Troposphere-constrained single-frequency PPP model TIC2 Ionosphere-constrained single-frequency PPP model TIC3 Troposphere-and ionosphere-constrained single-frequency PPP model Figures 8 and 9 show the percentage statistics of the effects of different constraint methods on the convergence time of GPS. It can be concluded from these two figures that in the horizontal component, the convergence time of the TIC3 has the largest percentage in the same range, and the TIC1 has less influence on the convergence time. In the vertical component, the percentage of convergence time of the TIC3 in the same range is still the largest, but the influence of TIC1 on the convergence time is higher than the percentage of TIC2 in a certain range. On the whole, the convergence performance of TIC0 is worse in both components when not using GIM products or tropospheric products as constraints. The convergence performance is improved as compared to TIC0 after considering TIC1 or (and) TIC2 in GPS SF PPP processing.     At the same time, Table 6 shows that TIC1 performs a notably better than TIC0 in the vertical component and marginally better than TIC0 in the horizontal component. TIC2 performs notably better than TIC0 and TIC1 in the horizontal component, while, in the vertical component, TIC2 is better than TIC0 and worse than TIC1 only in 68% level static mode. TIC3 performs best among the handling schemes. Compared with TIC0, the average convergence time of TIC1 at the 95% level is improved from 27.17 to 26.07 min (4.0%) in horizontal and by from 23.21 to 18.17 min (21.7%) in the vertical components. At the 68% level, the convergence time in the horizontal component reduced by 4.3% from 42.11 to 40.29 min and by 28.8% from 39.34 to 27.98 min in the vertical component. However, in kinematic mode, the average convergence time of TIC1 is not much higher than that of TIC0, with a maximum increase of 11.7% and a minimum increase of 2.1%. TIC2 has a larger boost when compared with TIC0. In static mode, it is increased by 16.6% in the vertical component by 68% level, and the maximum increase is 88.4% in the horizontal component of 95% level; in kinematic mode, they are 2.5% and 73.9%, respectively. TIC3 has the largest boost as compared with the other two methods. In static mode, it reaches 88.8% from 27.17 to 3.04 min and the lowest is 46.5% from 39.34 to 21.05 min In the kinematic mode, the minimum increase is 5.4%, and the maximum increase is 74.5%.
Three stations (FLRS, TRO1, and ZECK station) were randomly selected from 78 IGS stations to measure RMS error in the horizontal component, vertical component, and three-dimensional. It can be seen from Figures 10 and 11 that the final positioning error reaches the sub-decimetre level in the kinematic mode, and it reaches the centimetre level in the static mode. However, the final positioning errors are substantially close in the four schemes. At the same time, Table 6 shows that TIC1 performs a notably better than TIC0 in the vertical component and marginally better than TIC0 in the horizontal component. TIC2 performs notably better than TIC0 and TIC1 in the horizontal component, while, in the vertical component, TIC2 is better than TIC0 and worse than TIC1 only in 68% level static mode. TIC3 performs best among the handling schemes. Compared with TIC0, the average convergence time of TIC1 at the 95% level is improved from 27 However, in kinematic mode, the average convergence time of TIC1 is not much higher than that of TIC0, with a maximum increase of 11.7% and a minimum increase of 2.1%. TIC2 has a larger boost when compared with TIC0. In static mode, it is increased by 16.6% in the vertical component by 68% level, and the maximum increase is 88.4% in the horizontal component of 95% level; in kinematic mode, they are 2.5% and 73.9%, respectively. TIC3 has the largest boost as compared with the other two methods. In static mode, it reaches 88.8% from 27.17 to 3.04 min. and the lowest is 46.5% from 39.34 to 21.05 min. In the kinematic mode, the minimum increase is 5.4%, and the maximum increase is 74.5%.
Three stations (FLRS, TRO1, and ZECK station) were randomly selected from 78 IGS stations to measure RMS error in the horizontal component, vertical component, and three-dimensional. It can be seen from Figures 10 and 11 that the final positioning error reaches the sub-decimetre level in the kinematic mode, and it reaches the centimetre level in the static mode. However, the final positioning errors are substantially close in the four schemes.

Conclusions
In this study, we discussed the effects of four methods, including NENE, BILI, BICU, and JUNK, for the percentage of convergence time, the average convergence time and the calculation time of single station of SF TIC2 PPP. The numerical results show that: (1) the percentage of the convergence time of the estimated parameters through the BILI method is the shortest in every condition; and, (2) in the static mode, when the RMS is less than 5 cm, only the percentage of the BILI method calculation exceeds 40%. While in the kinematic mode, when the RMS is less than 15 cm, only the BILI method is more than 60% and 40% in the horizontal and vertical component, respectively. Under the other identical thresholds of the RMS, the percentage of the convergence time of the estimated parameters through BILI method is higher than the other three methods; (3) in the static mode, the average convergence time of the estimated parameters through BILI method is reduced by at least 2.2% when compared to the other three methods and the maximum is shortened by 23.6%; while, in the kinematic mode, they are 1.9% and 16.6%, respectively; and, (4) for the average computation time of a station, the BILI method has 0.121 s more than the shortest calculation time in the static mode and 0.038 s in the kinematic mode. Therefore, we chose the BILI method for calculating the TEC at the required time when using the TIC2 or TIC3 calculation based on the GIM grid products.
In order to verify the different constraint methods, a total of 18,720 tests were performed with 78 stations. The experimental results revealed the following findings: (1) the TIC1 convergence time

Conclusions
In this study, we discussed the effects of four methods, including NENE, BILI, BICU, and JUNK, for the percentage of convergence time, the average convergence time and the calculation time of single station of SF TIC2 PPP. The numerical results show that: (1) the percentage of the convergence time of the estimated parameters through the BILI method is the shortest in every condition; and, (2) in the static mode, when the RMS is less than 5 cm, only the percentage of the BILI method calculation exceeds 40%. While in the kinematic mode, when the RMS is less than 15 cm, only the BILI method is more than 60% and 40% in the horizontal and vertical component, respectively. Under the other identical thresholds of the RMS, the percentage of the convergence time of the estimated parameters through BILI method is higher than the other three methods; (3) in the static mode, the average convergence time of the estimated parameters through BILI method is reduced by at least 2.2% when compared to the other three methods and the maximum is shortened by 23.6%; while, in the kinematic mode, they are 1.9% and 16.6%, respectively; and, (4) for the average computation time of a station, the BILI method has 0.121 s more than the shortest calculation time in the static mode and 0.038 s in the kinematic mode. Therefore, we chose the BILI method for calculating the TEC at the required time when using the TIC2 or TIC3 calculation based on the GIM grid products.
In order to verify the different constraint methods, a total of 18,720 tests were performed with 78 stations. The experimental results revealed the following findings: (1) the TIC1 convergence time percentage is slightly higher than the TIC0 in each convergence time threshold, while the TIC2 is larger and the TIC3 is the largest. That is to say, the convergence can be accelerated under the constraint conditions; (2) as compared with TIC0 method, the average convergence time of the TIC1 method is shortened by at least 4.0% and, at most, 28.9% in static mode. The TIC2 method has a larger percentage of shortening, at least up to 16.6%, the maximum up to 88.4%. However, in the vertical component, the TIC1 method is larger than the TIC2 method at the 68% level. Compared with TIC0, the TIC3 method is shortened by 46.5% at least and 88.8% at most. Additionally TIC3 method is better than TIC1 method and TIC2 method. In kinematic mode, the percentage of time that is shortened after the constraint is less than in static mode. The average convergence time of the TIC1 method is at least 2.1% shorter than the TIC0 method, and the maximum is shortened by 11.7%. The TIC2 method has a larger percentage of shortening compared to the TIC1 method, with a minimum of 2.5% and a maximum of 73.9%. The TIC3 method is superior to the TIC1 method and the TIC2 method, with a minimum of 5.4% reduction and a maximum reduction of 74.5% as compared to TIC0; and, (3) through comparative analysis, it is concluded that, the final positioning error of TIC0, TIC1, TIC2 and TIC3 methods are basically the same, no method is better. After the above analysis and conclusions, under the constraints of atmospheric error products, convergence can be accelerated, and there is no impact on the final accuracy improvement.
Author Contributions: R.W., J.G., N.Z. and Z.L. conceived and designed the experiments; R.W. and Y.Y. performed the experiments and analyzed the data; L.Z. drew pictures; R.W. and Y.W. wrote the paper; all authors reviewed the paper.