Next Article in Journal
Matchup Characteristics of Sea Surface Salinity Using a High-Resolution Ocean Model
Next Article in Special Issue
A New Approach for the Development of Grid Models Calculating Tropospheric Key Parameters over China
Previous Article in Journal
Satellite-Derived Estimation of Grassland Aboveground Biomass in the Three-River Headwaters Region of China during 1982–2018
Previous Article in Special Issue
A New Method to Determine the Optimal Thin Layer Ionospheric Height and Its Application in the Polar Regions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modified Interpolation Method of Multi-Reference Station Tropospheric Delay Considering the Influence of Height Difference

1
State Key Laboratory of Geodesy and Earth’s Dynamics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430077, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(15), 2994; https://doi.org/10.3390/rs13152994
Submission received: 9 June 2021 / Revised: 22 July 2021 / Accepted: 26 July 2021 / Published: 29 July 2021
(This article belongs to the Special Issue BDS/GNSS for Earth Observation)

Abstract

:
In network real-time kinematic (NRTK) positioning, atmospheric delay information is critical for generating virtual observations at a virtual reference station (VRS). The traditional linear interpolation method (LIM) is widely used to obtain the atmospheric delay correction. However, even though the conventional LIM is robust in the horizontal direction of the atmospheric error, it ignores the influence of the vertical direction, especially for the tropospheric error. If the height difference between the reference stations and the rover is large and, subsequently, tropospheric error and height are strongly correlated, the performance of the traditional method is degraded for tropospheric delay interpolation at the VRS. Therefore, considering the height difference between the reference stations and the rover, a modified linear interpolation method (MLIM) is proposed to be applied to a conventional single Delaunay triangulated network (DTN). The systematic error of the double-differenced (DD) tropospheric delay in the vertical direction is corrected first. The LIM method is then applied to interpolate the DD tropospheric delay at the VRS. In order to verify the performance of the proposed method, we used two datasets from the American NOAA continuously operating reference stations (CORS) network with significant height differences for experiments and analysis. Results show that the DD tropospheric delay interpolation accuracy obtained by the modified method is improved by 84.1% and 69.6% on average in the two experiments compared to the conventional method. This improvement is significant, especially for low elevation satellites. In rover positioning analysis, the traditional LIM has a noticeable systematic deviation in the up component. Compared to the conventional method, the positioning accuracy of the MLIM method is improved in the horizontal and vertical directions, especially in the up component. The accuracy of the up component is reduced from tens of centimeters to a few centimeters and demonstrates better positioning stability.

1. Introduction

The traditional single-baseline real-time kinematic (RTK) positioning technology has been developed over the last several decades. However, due to the influence of atmospheric errors, orbital errors, and other distance-related errors, the distance between the rover and the reference station can only be about 10 km for centimeter-level positioning. The emergence of network real-time kinematic (NRTK), which is based on multiple reference stations, has expanded the scope of services for precise real-time positioning [1,2,3]. Virtual reference station (VRS) technology has been widely applied in NRTK data processing. When the ambiguities between reference stations are fixed in the network, distance-dependent errors, such as the double-differenced (DD) ionospheric and tropospheric delay of baselines, can easily be calculated. Using the estimated error among the reference stations, we can obtain the corresponding error between the reference stations and the VRS near the rover through the interpolation method. Finally, the generated virtual observations and the rover’s actual observations can be used for short baseline relative positioning, and the centimeter-level coordinates of the rover can be quickly obtained [4,5].
After extracting the atmospheric error from the baseline of the reference stations correctly, selecting the appropriate algorithm to generate the virtual observations at the VRS becomes a crucial issue. Therefore, researchers have proposed several interpolation models. Wanninger originally proposed a linear interpolation method (LIM) that requires at least three reference stations around the rover. The dual-frequency phase observations and the known coordinates of the reference station are applied to derive the regional DD ionospheric model [6]. The LIM can also build the model of other distance-dependent errors, such as tropospheric delay [7]. Han and Rizos introduced a linear combination model (LCM) of single-differenced (SD) observations, mainly used to eliminate orbital errors and that can be applied to model atmospheric errors [8]. The distance-based linear interpolation method (DIM), proposed by Gao et al., interpolates the rover’s ionospheric delay based on the distance between the reference stations and the rover [9]. Wübbena applied a suitable trend surface, called the low-order surface model (LSM), to simulate the distance-dependent error trend in the network. The coefficients of the LSM model can be obtained from the observation data of the reference stations by the least-squares adjustment [10]. The least-squares collocation method can also be applied to interpolate the atmospheric delay [11,12,13]. A comprehensive comparison of the advantages and disadvantages of these interpolation methods was undertaken by Dai et al. [14], Fotopoulos et al. [15], Wu et al. [16], and Al-Shaery et al. [17]. All results showed that the effects of these interpolation methods are similar. It is difficult to say which one performs the best. However, The LSM requires at least four stations to resolve the model coefficients. If the network is formed by Delaunay triangulation in VRS mode, it is difficult to use LSM in the smallest triangular unit. The rover station is generally surrounded by three reference stations in the traditional Delaunay triangulated network (DTN), so the LIM method is easier to implement for a triangulation model.
With the increase of the baseline length, the correlation of tropospheric errors among the stations decreases, which leads to large residual errors in DD observations. Significant residual errors will affect the ambiguity resolution (AR) and positioning quality. Since the tropospheric delay is highly correlated with elevation, its correlation coefficient exceeds 0.9, which will seriously affect the estimation of the station positions [18,19,20]. Landau used the experimental data of the German SAPOS network and found that when an obvious height difference exists between the rover and reference stations, the tropospheric delay systematic error caused by the height deviation can reach 6.8 cm [21]. If the interpolation algorithms mentioned above are applied in the smallest triangular unit directly, they only consider the horizontal distribution of the tropospheric errors. The tropospheric delay interpolated at the VRS is strongly constrained in the horizontal plane that is formed by the reference stations. The tropospheric delay in the vertical direction may be significantly inaccurate. Wu et al. proposed that the multi-baseline tropospheric delay interpolation method has better accuracy than the traditional LIM in the star network [22]. In addition, based on the BP artificial neural network method, Qiu et al. used a neural network to better correct the tropospheric delay and build the spatial tropospheric delay error model [23].
The above-mentioned interpolation methods either did not consider the systematic error of the tropospheric delay caused by the height difference between the user station and the reference station or required using multiple baselines and complex modeling methods, which are not convenient in practice. In our study, we proposed an improved linear interpolation method that was simple to implement and that was effective. Taking the height difference between the reference stations and the rover into account, the prior model of the troposphere was used to correct the estimated DD tropospheric delay on the network baseline first. After that, the LIM was used to interpolate the tropospheric delay between the VRS and the master reference station. Finally, the interpolated value was applied to generate virtual observations and positioning.
The main objective of this research is to improve the accuracy of the tropospheric delay correction using the interpolation algorithm in the NRTK of the Delaunay triangulated network (DTN) model in order to obtain more stable and more accurate positioning results through the modified method, especially in the vertical direction. The structure of this paper is as follows: In Section 2, the calculation method of the DD tropospheric delay of the network baseline is introduced, and the traditional LIM and MLIM methods are presented in detail. In Section 3, using the experimental data from the triangular units, the accuracy of the tropospheric delay correction when applying both the conventional and modified methods is analyzed with different satellite elevation angles. The difference in the positioning results from the two methods is compared. The results of the analyses are then discussed in Section 4. Finally, the conclusions are drawn, and further research directions are suggested in Section 4.

2. Materials and Methods

Here, we first describe the VRS’s tropospheric delay acquisition process, which includes the calculation of the DD tropospheric delay among the reference stations and the traditional modeling method of the tropospheric delay at the VRS. The modified method is then developed in Section 2.2, using the prior tropospheric model to correct the tropospheric elevation system error due to the height difference between the reference and the rover stations and combining the LIM interpolation algorithm with the calculated DD tropospheric delay to obtain the DD tropospheric delay between the VRS and the master station.

2.1. Calculation and Modeling Process of the Tropospheric Delay

In order to accurately obtain the tropospheric delay of the network baselines, we established the DD observation equations to obtain the wide-lane (WL) ambiguities, which can be expressed as follows:
{ L r b , W L i j = ρ r b i j + T r b i j + f 1 f 2 ( I r b , 1 j I r b , 1 i ) + λ W L ( N r b , W L j N r b , W L i ) + ε r b , W L i j L r b , 1 i j = ρ r b i j + T r b i j ( I r b , 1 j I r b , 1 i ) + λ 1 ( N r b , 1 j N r b , 1 i ) + ε r b , 1 i j P r b , 1 i j = ρ r b i j + T r b i j + ( I r b , 1 j I r b , 1 i ) + e r b , 1 i j P r b , 2 i j = ρ r b i j + T r b i j + f 1 2 f 2 2 ( I r b , 1 j I r b , 1 i ) + e r b , 2 i j I r b , 1 i j ˜ = I r b , 1 j I r b , 1 i
where the subscripts r and b identify the reference stations, and the superscripts i and j denote the reference and rover satellites, respectively. L r b , W L i j , L r b , 1 i j , P r b , 1 i j , and P r b , 2 i j represent the DD carrier and the pseudorange observations of the WL, and L1, P1, and P2 are in meters. ρ r b i j is the DD geometric distance from the satellite to receiver. T r b i j is the DD slant tropospheric delay of each satellite pair. I r b , 1 i and I r b , 1 j are the SD ionospheric delays on L1. N r b , W L i , N r b , W L j , N r b , 1 i , and N r b , 1 j are the SD WL and L1 ambiguities. f 1 and f 2 are the frequencies of L1 and L2, and λ W L and λ 1 are the wavelengths of WL and L1. ε r b , W L i j , ε r b , 1 i j , e r b , 1 i j , and e r b , 2 i j denote all other errors, including the noise of the DD carrier and the pseudorange. I r b , 1 i j ˜ is DD ionospheric pseudo-observation.
Using WL ambiguity obtained by Equation (1), the L1 ambiguity was then estimated using the following observation equation:
L r b , I F i j = ρ r b i j + T r b i j + C f 1 + f 2 ( N r b , 1 j N r b , 1 i ) + C f 2 f 1 2 + f 2 2 N r b , W L i j + ε r b , I F i j P r b , I F i j = ρ r b i j + T r b i j + e r b , I F i j
where IF denotes the ionosphere-free combination. N r b , W L i j is DD WL ambiguities. The other symbols are the same as in Equation (1).
The prerequisite for obtaining the DD tropospheric delay is to resolve the ambiguity correctly. Generally, the baseline difference of the reference stations is tens of kilometers or even more than 100 km, and the AR is difficult. Therefore, we used Equation (1) to estimate the SD WL ambiguity, and the SD L1 ambiguity was then estimated by applying Equation (2). The ionospheric delay is one of the major factors restricting the AR among the reference stations. This effect is handled with a temporal and spatial correlation weighting scheme for the ionospheric pseudo-observations in Equation (1) [13,24]. The DD ionospheric pseudo-observations can be obtained from the broadcast Klobuchar ionospheric model [25] or the IGS global ionosphere models (GIM). We selected the GIM model because it has better accuracy than the broadcast Klobuchar model [26]. On the other hand, the UNB3m model combined with Niell mapping functions (NMF) was used to compute the hydrostatic tropospheric corrections [27,28]. The wet tropospheric delay was estimated as the relative tropospheric zenith delay (RZTD) parameter [29].
The Kalman filter was adopted to estimate the float SD WL, L1 ambiguity, the SD ionospheric delay, and the RZTD. The DD ambiguity was obtained from the SD ambiguity through the transition matrix, and the LAMBDA method was applied to obtain the integer of the DD WL and L1 ambiguities [30]. The L2 DD ambiguity could then be directly derived from the linear relationship between the WL and L1 DD ambiguities.
After we obtained the L1 and L2 DD integer ambiguities, the tropospheric slant delay corresponding to each satellite pair between r and b was calculated as follows:
T ^ r b i j = f 1 2 f 1 2 f 2 2 ( L r b , 2 i j λ 1 N ˜ r b , 1 i j ) f 2 2 f 1 2 f 2 2 ( L r b , 2 i j λ 2 N ˜ r b , 2 i j ) ρ r b i j
where T ^ r b i j represents the calculated DD tropospheric delay. N ˜ r b , 1 i j and N ˜ r b , 2 i j denote fixed DD ambiguities on L1 and L2 in units of cycles. Observation noise, multipath effects, and other errors are ignored here.
Generally, at least three reference stations around the rover are used to build a DD tropospheric interpolation model, and the reference station closest to the rover is selected as the master reference station, while the other two are auxiliary reference stations. The interpolation model for applying the LIM can be expressed as
T ^ = B X ab
T ^ = T ^ M , 1 i j T ^ M , 2 i j T ^ M , n 1 i j B = Δ X M , 1 Δ Y M , 1 Δ X M , 2 Δ Y M , 2 Δ X M , n 1 Δ Y M , n 1   X a b = a b
where i, j denote the reference and rover satellite, subscript M is the master reference station, and 1, 2, …, n − 1 represent the auxiliary reference stations. Δ X and Δ Y are the difference of the horizontal coordinates between the master and auxiliary reference stations. T ^ is the calculated DD tropospheric delay, and B is the coefficient matrix. When the DD tropospheric delay is determined by (3), the coefficients a and b in (5) can be estimated by least-squares adjustment:
X a b = B T B 1 B T T ^
The DD tropospheric delay between the VRS and the master reference station can then be obtained as follows:
T M , V ij = Δ X M , V Δ Y M , V X ab = Δ X M , V Δ Y M , V B T B 1 B T T ^
where T M , V ij is the interpolated DD tropospheric delay between the master reference station and the VRS for the satellite pair i, j. Δ X M , V and Δ Y M , V are the difference in the horizontal coordinates between the VRS and the master reference station.

2.2. The Modified Tropospheric Interpolation Algorithm

In GNSS data processing, the ionosphere is assumed to be a single-layer model above the earth, and its height is generally 350 km. For a small area, the point where the satellite signal passes through the single layer will form an approximate plane. Therefore, when ionospheric changes are not active, using conventional LIM to interpolate the ionospheric delay will achieve a sufficient accuracy in the NRTK. Compared to the ionospheric delay, the tropospheric temporal and spatial distribution characteristics are significantly different. The tropospheric delay is affected by both horizontal and height factors. It can be seen from Section 2.1 that the traditional interpolation method only uses horizontal coordinates to interpolate tropospheric errors, but ignores the vertical element. This method constrains the interpolation of the tropospheric delay to the height level constituted by reference stations. When the height of the rover is much higher or lower than the selected reference stations, it causes the positioning performance to decrease, especially in the up component. Therefore, we proposed that the tropospheric system error caused by the height difference between the reference stations and rovers should be corrected first. The LIM was then used for the DD tropospheric delay modeling.
Assuming that the reference stations M, A, and B form a triangular unit and assuming that station M is selected as the master reference station, A and B are the two auxiliary reference stations. The coordinates of the VRS are the result of the standard point positioning (SPP) of the rover. Suppose the height of the VRS is much higher or lower than the reference stations. In that case, the tropospheric delay between the master station and the auxiliary reference stations calculated by Equation (3) needs to be corrected in advance before applying the LIM. Take the baseline between A and M as an example. The calculated DD tropospheric delay is
T M A i j ( Δ h M A ) = T A i j ( h A ) T M i j ( h M )
where T M A i j ( Δ h M A ) is the calculated DD tropospheric delay of the reference and rover satellite, i and j are on the baseline MA with the height difference Δ h M A . T A i j ( h A ) , and T M i j ( h M ) are the SD tropospheric delays of the reference and rover satellites i, j at stations A and M, with heights of h A and h M , respectively.
The SD tropospheric delay correction equations for station A and M are as follows:
T ¯ A i j ( h A ) = T A i j ( h A ) + T C o r A
T ¯ M i j ( h M ) = T M i j ( h M ) + T C o r M
where T ¯ A i j ( h A ) and T ¯ M i j ( h M ) are the corrected SD tropospheric delays for stations A and M. The specific derivations of the correction terms T C o r A and T C o r M are as follows:
T C o r A = T A i j ( h V ) T A i j ( h A ) = Z T D A ( h V ) × ( M F A j ( h V ) M F A i ( h V ) ) Z T D A ( h A ) × ( M F A j ( h A ) M F A i ( h A ) )
T C o r M = T M i j ( h V ) T M i j ( h M ) = Z T D M ( h V ) × ( M F M j ( h V ) M F M i ( h V ) ) Z T D M ( h M ) × ( M F M j ( h M ) M F M i ( h M ) )
where T C o r A and T C o r M are the correction items, and h A and h V are the height of station A and the VRS, respectively. Z T D A ( h A ) and Z T D A ( h V ) represent the ZTD at station A with height h A and h V , respectively, which can be calculated using the prior model. M F A j ( h A ) , M F A i ( h A ) , M F A j ( h V ) and M F A j ( h V ) denote the tropospheric mapping functions corresponding to the satellites i, j at station A with the height h A and h V , respectively.
Although the ZTD calculated by the prior model has model errors, the difference of the ZTD obtained by the prior model at the same horizontal position but at different heights can effectively eliminate the tropospheric system error caused by the height difference. Therefore, the T C o r A and T C o r M obtained by the prior model still have a high correction accuracy. We substituted Equations (11) and (12) into Equations (9) and (10), respectively, and replaced T A i j ( h A ) and T M i j ( h M ) in Equation (8) with T ¯ A i j ( h A ) and T ¯ M i j ( h M ) in Equations (9) and (10). Equation (8) then became
T ¯ M A i j ( Δ h V V ) = T ¯ A i j ( h A ) T ¯ M i j ( h M ) = T A i j ( h V ) T M i j ( h V )
where T ¯ M A i j ( Δ h V V ) is the corrected DD tropospheric delay of the reference and rover satellites i, j on the baseline MA.
The mean error of the ZTD calculated by the UNB series model in North America is about 2 cm [31], which represents good accuracy. The NMF and other mapping functions have similar accuracy. Taking the experimental data from North America into account, the tropospheric systematic deviation correction was calculated using the UNB3m model with the NMF in this paper. After correcting the systematic error of the DD tropospheric delay in advance, and then using the LIM for interpolation, a DD tropospheric delay with no height difference deviation on the VRS baseline and the master station could be obtained.

3. Results

In this chapter, we first describe the two experimental datasets used to verify the modified interpolation algorithm. After that, in the first experiment, we chose different periods and satellites in various triangular units in Dataset One to analyze the tropospheric delay interpolation effect of the traditional and modified methods. In the second experiment, several satellites in the same period in the same triangle unit in Dataset Two were selected to analyze the interpolation effects. Finally, the positioning performance after adopting the traditional and proposed methods was analyzed.

3.1. Experimental Data

The experimental data were from the US CORS network. We chose 21 stations as reference stations and another ten stations as simulated rovers for experimental analysis. These data were divided into two experimental datasets.
Dataset One consisted of six independent triangular units, and each unit contained a simulated rover. Dataset Two only had one triangle unit, but it included four simulated rover stations. The plane distribution and baseline length of Dataset One and Dataset Two are shown in Figure 1 and Figure 2. The reference station closest to each rover was selected as the master reference station; they are connected by the green dashed line in the figure. Figure 3 and Figure 4 show the station height distribution of the two datasets. The observation data on day of year (DOY) 180 in 2020 were selected. The sampling interval was set to 30 s, and the cut-off angle was set to 15°.
The coordinates of the reference stations and the simulated rovers were obtained from the US CORS website and were used as known coordinates in the experiments. The Kp index, proposed by Bartels, is an indicator of geomagnetic activity and is updated every three hours [32]. We checked the Kp value of DOY 180 in 2020; the maximum was 2 and the minimum was 0, which meant that the ionosphere was quiet on that day. Therefore, the DD ionospheric delay obtained by the LIM had high accuracy, and it did not affect the tropospheric delay interpolation. Multipath effects, antenna phase center correction, and other errors were ignored.
The slant DD tropospheric delay between the auxiliary and master stations is calculated using Equation (3). The traditional LIM and modified LIM are then applied to obtain the slant DD tropospheric delay between the rover and master stations. The true value of ZTD between the rover and the master reference stations could be obtained using the latest Canadian Spatial Reference System precise point positioning (CSRS-PPP) service, provided by Natural Resources Canada [33]. Since the estimation accuracy could reach 0–2 cm [34,35], the slant DD tropospheric delay of each satellite pair calculated in combination with the mapping function could be used as the reference true value. The whole experiment was conducted in post-processing mode.

3.2. Analysis of Tropospheric Delay Interpolation Results with Dataset One

The triangular units of Dataset One were distributed in different areas, and the height distribution of the reference stations and rovers also had apparent differences. They were used to verify the effectiveness of the proposed method compared to the traditional LIM. The DD tropospheric delay was obtained using the CSRS-PPP as the true value. The values of the root mean square (RMS) of the tropospheric delay obtained by the two methods were analyzed and compared.
In each triangle unit, we selected different periods, different elevation angle changes, and three continuously observed satellites to analyze the effect of tropospheric interpolation. The results were obtained using both the traditional LIM and the modified LIM, which are denoted here as TLIM and MLIM. A comparison with the true value is shown in Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10.
We selected three satellites with different elevation angle range changes on each rover for the analysis. The trend of elevation angle changes for these selected satellites was also different. Overall, from Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10, we can see that as the elevation angles of each satellite increased, the DD tropospheric delay became smaller, and vice versa. These selected satellites can be divided into three types: low-elevation, medium-elevation, and high-elevation angle satellites.
For high elevation angle satellites, we took G10, G18, and G14 in Figure 5, Figure 6 and Figure 7 as examples. The DD tropospheric delay obtained using TLIM at the VRS had a slight offset from the true value. However, the DD tropospheric delay corrected using MLIM almost coincided with the true value. For G14, which had the largest elevation angle increase, the results of the proposed method were the most consistent with the true value.
Similarly, we took G18, G20, and G14 in Figure 5, Figure 7 and Figure 8 as medium elevation angle satellites for analysis. These three satellites showed a downward trend in the range of medium elevation angles. The offset between the DD tropospheric delay obtained by the TLIM and true value was obvious. As the elevation angles of these satellites decreased, these offsets became larger. The changing trend of the TLIM value and the true value was not the same. When MLIM was adopted, the offsets between the interpolation results of the MLIM and true value were mostly eliminated, and a better correction effect was obtained.
From Figure 7, Figure 9 and Figure 10, the low elevation angle satellites G08, G16, and G11 were analyzed. Their elevation angles were relatively low, and they showed almost no change trend. Therefore, their tropospheric delay change trends according to the TLIM, MLIM, and true values were not noticeable. However, their tropospheric delay values were considerable. This is because the lower the elevation, the longer the satellite signal propagation path becomes, and the delay of signal propagation will increase the range of the troposphere. The most significant offsets between TLIM and the true value occurred with this type of elevation satellite. Although the DD tropospheric delay obtained by the MLIM still had a small difference from the true value, it effectively eliminated this offset compared to the TLIM.
It can be seen that whether the DD tropospheric delay interpolated by the proposed method is from high elevation or low elevation satellites, the offset with the true value is smaller. For all of the satellites, the trend of the MLIM results was consistent with the true value.
The statistical analysis of the results of the TLIM and MLIM for each rover is shown in Figure 11. Overall, the RMS errors of the MLIM are smaller than those of the TLIM. This indicates that the accuracy of the interpolated tropospheric delay improved after applying the proposed method. The RMS errors of the MLIM on each continuously observed satellite showed a significant improvement. The maximum increases were achieved on the G32, G29, G08, G11, G16, and G11 satellites at stations P478, EWPP, P603, P225, WINT, and P121, respectively, with an improvement of 90.1%, 92.5%, 90.0%, 93.2%, 90.6%, and 78.5%.
To analyze all tropospheric interpolation accuracies for all of the satellites with different elevation angles during the selected period for each rover, Figure 12 shows the TLIM and MLIM modeling error curves for all of the satellites with different elevation angles. Overall, from Figure 12, we can see the accuracy of the tropospheric modeling using TLIM increased with increasing satellite elevation. However, the accuracy of the MLIM was always better than that of the TLIM at any elevation interval. For low elevation satellites, the error of the TLIM was obvious. For example, the error forP603 reached 0.2 m, but the error from the MLIM was still minimal. The mean error values of the TLIM were 0.017 m, −0.027 m, 0.068 m, 0.029 m, 0.045 m, and 0.010 m, and the corresponding standard deviations (STD) were 0.014 m, 0.023 m, 0.048 m, 0.023 m, 0.031 m, and 0.006 m, respectively. The accuracy of the MLIM was significantly improved, especially for low elevation satellites. The mean errors of each rover were 0.002 m, 1.65 × 10−5 m, 0.007 m, 0.005 m, 0.005 m, and 0.003 m, and the corresponding STDs were 0.002 m, 0.004 m, 0.006 m, 0.004 m, 0.004 m, and 0.002 m, respectively. The discontinuous part in Figure 12 indicates that there were no visible or ambiguity fixed satellites within that elevation angle range.
The mean RMS errors of the DD tropospheric delay obtained by the two methods for all satellites are listed in Table 1. The statistical analysis showed that the RMS of the TLIM reached tens of centimeters, while the RMS of the interpolation using the proposed method was only several millimeters. The improvement reached 87.5%, 87.9%, 89.5%, 85.4%, 87.9%, and 66.7%, respectively. This indicated that the modified algorithm performed better than the traditional algorithm, especially for satellites with lower elevation angles.

3.3. Analysis of the Tropospheric Delay Interpolation Results with Dataset Two

Dataset Two only had one triangle unit, but it contained four rovers with different elevations. We also compared and analyzed the performance of the TLIM and MLIM at these stations.
We selected three identical satellites, G10, G18, and G27, at these four rovers to analyze the effect tropospheric interpolation during the same period. The results are illustrated in Figure 13.
As shown in Figure 13, the interpolated DD tropospheric delay by the MLIM was highly consistent with the true value for G10 at the highest elevation. For G18, which had an elevation angle on a downward trend, although there were still some deviations compared to the true value, the MLIM results were also closer to the true value. Analyzing the satellite G27, although its elevation angle was increasing, it had the lowest elevation angle among the three satellites. Its TLIM results showed that the effect of directly interpolating the tropospheric delay was not satisfactory. However, after adopting the MLIM, the value of the DD tropospheric delay interpolation was consistent with the true value. Therefore, the modified interpolation method had a better correction effect, even for low-elevation satellites.
The RMS of the DD tropospheric delays obtained using the TLIM and MLIM for each rover during the selected period are shown in Figure 14. The maximum improvement values of the MLIM at the four rovers compared to the TLIM were 93.0%, 83.4%, 78.3%, and 79.3%, respectively.
Further, troposphere modeling errors for all of the satellites over the selected observation period are shown in Figure 15. The mean errors of the TLIM were 0.011 m, 0.027 m, 0.029 m, and 0.021 m, and the corresponding STDs were 0.010 m, 0.027 m, 0.029 m, and 0.020 m, respectively. The accuracy of the MLIM was dramatically improved, with mean errors of −0.002 m, −0.006 m, 0.009 m, and 0.006 m, and the corresponding STDs were 0.007 m, 0.007 m, 0.008 m, and 0.006 m, respectively.
The mean RMS of the DD tropospheric delay using the traditional and the modified methods for all of satellites in Dataset Two is shown in Table 2. It can be seen that the mean RMS error of the TLIM of each station, P532, P539, P602, and MASW, were 0.012 m, 0.033 m, 0.035 m, and 0.026 m, respectively. However, the mean RMS error of the MLIM were 0.005 m, 0.008 m, 0.010 m, and 0.007 m. These values represent improvements of 58.3%, 75.6%, 71.4%, and 73.1%, respectively, which are consistent with the requirements of the accuracy of DD tropospheric delay interpolation that need to be met in order to generate VRS observations in the NRTK.

3.4. Comparison of the Positioning Results of Dataset One and Dataset Two

According to the proposed and traditional tropospheric interpolation delay results, the positioning performance of these rovers using the VRS was investigated in this section. The DD tropospheric delay was obtained through both the TLIM and MLIM. We then used these interpolated atmospheric delays to generate virtual carrier phase and pseudorange observations, which were applied to the final rover positioning. The RMS errors of the positioning results were analyzed in this paper. We expected that the proposed method would show a better positioning performance, especially in the up component and under significant height differences.

3.4.1. Dataset One

The positioning error of the TLIM and MLIM is shown in Figure 16, using station P225 as an example. It can be seen that the positioning error of the MLIM showed a significant improvement in the up component. However, there was little difference between the TLIM and MLIM in the horizontal direction. Only the results of the fixed solutions were counted here, and the MLIM had more fixed solutions than the TLIM. The positioning results of each rover in Dataset One are given in Table 3. In terms of the RMS, the accuracy in the three directions for all of the rovers was improved after using the MLIM. There were slight improvements in the horizontal direction, but remarkable improvements were achieved in the up component. All the RMS errors of the up component obtained using the MLIM were within 5 cm. The most obvious improvement occurred at P603, in which the accuracy of the RMS, changing from 24.95 cm to 4.24 cm, improved by 83%. It can be seen from Figure 16 that the positioning result became more stable after the proposed method was applied.

3.4.2. Dataset Two

Similarly, taking MASW as an example, it can be seen from Figure 17 that the accuracy improvement in the up direction was obvious. Table 4 shows the statistical positioning results of all of the rovers in Dataset Two. Analyzing the RMS results of the four rover stations, the accuracy of the TLIM was within 2 cm in the east and north directions. The accuracy of the up component was about 17 cm. After applying the MLIM, the positioning accuracy of all rovers was within 1 cm in the east and north directions. The accuracy in the up direction was within 5 cm, reaching the general accuracy level of the NRTK positioning. Figure 17 and Table 4 showed that in the case of the height distribution of Dataset Two, the positioning results using the modified interpolation methods were reliable and stable.
The results of these experiments demonstrate that the accuracy of the traditional LIM DD tropospheric delay at the VRS deteriorates when a relatively significant height difference exists between the rover and the reference stations. The results also confirm that as the satellite elevation angle drops, the delay of the satellite signal passing through the troposphere increases, and the effect of interpolating the tropospheric delay will worsen.
It can be seen from these results that the more significant the height difference between the rover and the reference station, the worse the tropospheric accuracy obtained by using TLIM is, and the better the correction effect by using MLIM is. The average RMS error of the conventional method was 0.035 m. The MLIM’s average RMS was only 0.006 m.
According to the results when using VRS observations for rover positioning, the modified method showed a slight improvement in the horizontal direction compared to the traditional method. However, the accuracy of the proposed method was significantly improved in the vertical direction. Therefore, when building or selecting reference station sites, one should thoroughly consider whether the height distribution of the reference station is approximately consistent. In addition, the height difference between the potential users and the reference stations should be as small as possible. Otherwise, a method similar to the MLIM should be used to solve the impact of the height difference in the NRTK.
It should be noted that this experiment only verified that a considerable height difference between the reference station and the rover would cause systematic errors when modeling the DD tropospheric delay in the NRTK. The heights of the reference stations and the rovers were quite different: Dataset One was from −29 m to about 1900 m, and Dataset Two was from 34 m to 713 m. The maximum height difference between the rover and reference stations of each triangle unit is from 245 m to 1158 m in Dataset One, and the maximum height difference is from 499 m to 679 m in Dataset Two. In addition, when the LIM is used to interpolate the atmospheric delay, the spatial distance and height difference between the rover and each reference station should be treated as the weighting factor for interpolating the DD tropospheric delay in each triangle unit. This perspective should be studied through further experiments in the future.

4. Conclusions

This study confirmed that a height difference between reference stations and rovers leads to a decrease of interpolation DD tropospheric delay accuracy. Therefore, based on the tropospheric delay correction model and traditional LIM, we proposed the MLIM correction method for the triangular VRS mode. This method corrects the systematic error of the tropospheric delay caused by height differences in advance. It then applies the LIM to interpolate the corrected value to obtain the DD tropospheric delay between the rover and the master stations. We compared the accuracy of TLIM and MLIM by modeling the DD tropospheric delay by using two datasets of triangular units in VRS mode. Finally, the positioning performance of the two methods was evaluated.
In the first experiment, we selected six independent triangular units in different areas. The accuracy of the DD tropospheric delay modeling, which used traditional LIM, was inferior. This phenomenon is evident for satellites with low elevation angles. Compared to traditional LIM, the accuracy of interpolated DD tropospheric delay was significantly improved after using the modified method, especially for the satellites with lower elevation angles. The second experiment considered only one triangular unit, but it included rover stations with different elevation distributions. From the perspective of the tropospheric interpolation results, as the height of each rover increased, the accuracy of the tropospheric interpolation using the conventional method deteriorated. After adopting the MLIM, the RMS error was significantly improved.
In terms of positioning performance, the up component accuracy obtained by the TLIM showed a noticeable systematic deviation, and the worst results in the two datasets could still reach the decimeter level. However, this situation improved significantly after using the MLIM, and the positioning results of all of stations in the up direction reached just a few centimeters. In the horizontal direction, although the accuracy of the proposed algorithm was similar to that of the traditional algorithm, there was still a slight improvement. In addition, the positioning result of the proposed method is more stable than that of the conventional method.
In summary, the proposed method significantly improved the tropospheric delay interpolation and the positioning performance in the NRTK triangulation network. It should be noted that this paper only analyzed the influence of the DD tropospheric delay interpolation in the NRTK where the height distribution of the reference stations and the rover stations were quite different, and further research should be devoted to investigating how the height difference of the reference stations and the spatial distance between the rover and reference stations can affect the DD tropospheric interpolation accuracy in the NRTK.

Author Contributions

Y.P. provided the initial idea and wrote the manuscript; Y.P. and M.S. designed and performed the research; M.S. helped with the writing and revising, and Y.Y. and M.S. partially financed the research. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Research Program (No. 2016YFB0501900) and China Natural Science Funds (No. 41604033 and 41874093). The third author is supported by the Wang Kuancheng Pioneer Talents Project Lu Jiaxi International Team Program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. The GPS observation data from the NOAA CORS network are available at https://www.ngs.noaa.gov/UFCORS/ (accessed on 26 April 2021). The GPS broadcast ephemeris data are available at https://cddis.nasa.gov/archive/gnss/data/daily/ (accessed on 26 April 2021). The GIM products from IGS can be obtained at https://cddis.nasa.gov/archive/gnss/products/ionex/ (accessed on 26 April 2021).

Acknowledgments

The authors gratefully acknowledged the Canadian Spatial Reference System Precise Point Positioning (CSRS-PPP) service, the NOAA CORS network for providing the GNSS data, and the IGS for providing the GIM products.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wanninger, L. Virtual reference stations (VRS). GPS Solut. 2003, 7, 143–144. [Google Scholar] [CrossRef]
  2. Bakuła, M.; Przestrzelski, P.; Kaźmierczak, R. Reliable technology of centimeter GPS/GLONASS surveying in forest environments. IEEE Trans. Geosci. Remote Sens. 2014, 53, 1029–1038. [Google Scholar] [CrossRef]
  3. Paziewski, J.; Sieradzki, R.; Baryla, R. Multi-GNSS high-rate RTK, PPP and novel direct phase observation processing method: Application to precise dynamic displacement detection. Meas. Sci. Technol. 2018, 29, 035002. [Google Scholar] [CrossRef]
  4. Hu, G.R.; Khoo, H.S.; Goh, P.C.; Law, C.L. Development and assessment of GPS virtual reference stations for RTK positioning. J. Geod. 2003, 77, 292–302. [Google Scholar] [CrossRef]
  5. Chen, X.; Han, S.; Rizos, C.; Goh, P.C. Improving real time positioning efficiency using the Singapore integrated multiple reference station network (SIMRSN). In Proceedings of the 13th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS 2000), Salt Lake City, UT, USA, 19–22 September 2000. [Google Scholar]
  6. Wanninger, L. Improved ambiguity resolution by regional differential modelling of the ionosphere. In Proceedings of the 8th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 1995), Palm Springs, CA, USA, 12–15 September 1995. [Google Scholar]
  7. Wanninger, L. The performance of virtual reference stations in active geodetic GPS-networks under solar maximum conditions. In Proceedings of the 12th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 1999), Nashville, TN, USA, 14–17 September 1999. [Google Scholar]
  8. Han, S.; Rizos, C. GPS network design and error mitigation for real-time continuous array monitoring systems. In Proceedings of the 9th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 1996), Kansas City, MO, USA, 17–20 September 1996. [Google Scholar]
  9. Gao, Y.; Li, Z.; McLellan, J. Carrier phase based regional area differential GPS for decimeter-level positioning and navigation. In Proceedings of the 10th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 1997), Kansas City, MO, USA, 16–19 September 1997. [Google Scholar]
  10. Wübbena, G.; Bagge, A.; Seeber, G.; Böder, V.; Hankemeier, P. Reducing distance dependent errors for real-time precise DGPS applications by establishing reference station networks. In Proceedings of the 9th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 1996), Kansas City, MO, USA, 17–20 September 1996. [Google Scholar]
  11. Odijk, D.; Marel, H.V.D.; Song, I. Precise GPS Positioning by Applying Ionospheric Corrections from an Active Control Network. GPS Solut. 2000, 3, 49–57. [Google Scholar] [CrossRef]
  12. Marel, H.V. Virtual GPS reference stations in the Netherlands. In Proceedings of the 11th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 1998), Nashville, TN, USA, 15–18 September 1998. [Google Scholar]
  13. Odijk, D. Fast precise GPS positioning in the presence of ionospheric delays. Publ. Geod. 2002, 52. [Google Scholar]
  14. Dai, L.; Han, S.; Wang, J.; Rizos, C. Comparison of interpolation algorithms in network-based GPS techniques. NAVIGATION J. Inst. Navig. 2003, 50, 277–293. [Google Scholar] [CrossRef]
  15. Fotopoulos, G.; Cannon, M. An overview of multi-reference station methods for cm-level positioning. GPS Solut. 2001, 4, 1–10. [Google Scholar] [CrossRef]
  16. Wu, S. Performance of regional atmospheric error models for NRTK in GPSnet and the implementation of a NRTK system. RMIT Univ. Aust. 2009. [Google Scholar]
  17. Al-Shaery, A.; Lim, S.; Rizos, C. Investigation of different interpolation models used in Network-RTK for the virtual reference station technique. J. Glob. Position. Syst. 2011, 10, 136–148. [Google Scholar] [CrossRef] [Green Version]
  18. Huang, L.; Liu, L.; Yao, C. A zenith tropospheric delay correction model based on the regional CORS network. Geod. Geodyn. 2012, 3, 53–62. [Google Scholar]
  19. Wielgosz, P.; Cellmer, S.; Rzepecka, Z.; Paziewski, J.; Grejner-Brzezinska, D.A. Troposphere modeling for precise GPS rapid static positioning in mountainous areas. Meas. Sci. Technol. 2011, 22, 045101. [Google Scholar] [CrossRef]
  20. Yin, H.; Huang, D.; Xiong, Y. Regional tropospheric delay modeling based on GPS reference station network. In VI Hotine-Marussi Symposium on Theoretical and Computational Geodesy; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  21. Landau, H.; Vollath, U.; Chen, X. Virtual reference stations versus broadcast solutions in network RTK–advantages and limitations. In Proceedings of the GNSS, Graz, Austria, 22–24 April 2003. [Google Scholar]
  22. Wu, B.; Gao, C.; Pan, S.; Deng, J.; Gao, W. Regional modeling of atmosphere delay in network rtk based on multiple reference station and precision analysis. In Proceedings of the China Satellite Navigation Conference (CSNC) 2015 Proceedings, Xi'an, China, 13–15 May 2015; Springer: Berlin/Heidelberg, Germany, 2015; Volume II. [Google Scholar]
  23. Lei, Q.; Lei, L.; Zemin, W. An tropospheric delay model for GPS NET RTK. In Proceedings of the 2010 Second International Conference on Information Technology and Computer Science, Kiev, Ukraine, 24–25 July 2010. [Google Scholar]
  24. Hu, G.; Abbey, D.A.; Castleden, N.; Featherstone, W.E.; Earls, C.; Ovstedal, O.; Weihing, D. An approach for instantaneous ambiguity resolution for medium-to long-range multiple reference station networks. GPS Solut. 2005, 9, 1–11. [Google Scholar] [CrossRef]
  25. Klobuchar, J.A. Ionospheric effects on GPS. Glob. Position. Syst. Theory Appl. 1991, 1, 517–546. [Google Scholar]
  26. Vstedal, O. Absolute Positioning with Single-Frequency GPS Receivers. GPS Solut. 2002, 5, 33–44. [Google Scholar] [CrossRef]
  27. Leandro, R.F.; Langley, R.B.; Santos, M.C. UNB3m_pack: A neutral atmosphere delay package for radiometric space techniques. GPS Solut. 2008, 12, 65–70. [Google Scholar] [CrossRef]
  28. Niell, A.E. Global Mapping Functions for the Atmosphere Delay at Radio Wavelengths. J. Geophys. Res. Solid Earth 1996, 101, 3227–3246. [Google Scholar] [CrossRef]
  29. Zhang, J.; Lachapelle, G. Precise estimation of residual tropospheric delays using a regional GPS network for real-time kinematic applications. J. Geod. 2001, 75, 255–266. [Google Scholar] [CrossRef]
  30. Teunissen, P.J.G. The least-square ambiguity decorrelation adjustment: A method for fast GPS ambiguity estimation. J. Geod. 1995, 70, 1–2. [Google Scholar] [CrossRef]
  31. Yao, Y.; Xu, X.; Xu, C.; Peng, W.; Wan, Y. Establishment of a Real-Time Local Tropospheric Fusion Model. Remote Sens. 2019, 11, 1321. [Google Scholar] [CrossRef] [Green Version]
  32. Bartels, J. The technique of scaling indices K and Q of geomagnetic activity. Ann. Intern. Geophys. 1957, 4, 215–226. [Google Scholar]
  33. Mireault, Y.; Tétreault, P.; Lahaye, F.; Héroux, P.; Kouba, J. Online precise point positioning. GPS World 2008, 19, 59–64. [Google Scholar]
  34. Guo, Q. Precision comparison and analysis of four online free PPP services in static positioning and tropospheric delay estimation. GPS Solut. 2015, 19, 537–544. [Google Scholar] [CrossRef]
  35. Xiao, G.; Liu, G.; Ou, J.; Liu, G.; Wang, S.; Guo, A. MG-APP: An open-source software for multi-GNSS precise point positioning and application analysis. GPS Solut. 2020, 24, 1–13. [Google Scholar] [CrossRef]
Figure 1. Geographical distribution of the selected stations in Dataset One.
Figure 1. Geographical distribution of the selected stations in Dataset One.
Remotesensing 13 02994 g001
Figure 2. Geographical distribution of the selected stations in Dataset Two.
Figure 2. Geographical distribution of the selected stations in Dataset Two.
Remotesensing 13 02994 g002
Figure 3. The elevation distribution of the stations in Dataset One.
Figure 3. The elevation distribution of the stations in Dataset One.
Remotesensing 13 02994 g003
Figure 4. The elevation distribution of the stations in Dataset Two.
Figure 4. The elevation distribution of the stations in Dataset Two.
Remotesensing 13 02994 g004
Figure 5. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G10, G18, and G32 at the rover station P478.
Figure 5. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G10, G18, and G32 at the rover station P478.
Remotesensing 13 02994 g005
Figure 6. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G15, G18, and G21 at the rover station EWPP.
Figure 6. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G15, G18, and G21 at the rover station EWPP.
Remotesensing 13 02994 g006
Figure 7. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G08, G14, and G20 at the rover station P603.
Figure 7. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G08, G14, and G20 at the rover station P603.
Remotesensing 13 02994 g007
Figure 8. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G01, G14, and G31 at the rover station P225.
Figure 8. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G01, G14, and G31 at the rover station P225.
Remotesensing 13 02994 g008
Figure 9. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G16, G22, and G26 at the rover station WINT.
Figure 9. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G16, G22, and G26 at the rover station WINT.
Remotesensing 13 02994 g009
Figure 10. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G07, G11, and G13 at the rover station P121.
Figure 10. The DD tropospheric delay obtained by different methods (a) and elevation angles (b) for satellites G07, G11, and G13 at the rover station P121.
Remotesensing 13 02994 g010
Figure 11. RMS errors of the DD tropospheric delay using modified and traditional methods for all non-reference satellites at all rover stations. The percentages on each bar are the improvements of the MILM with respect to the TLIM.
Figure 11. RMS errors of the DD tropospheric delay using modified and traditional methods for all non-reference satellites at all rover stations. The percentages on each bar are the improvements of the MILM with respect to the TLIM.
Remotesensing 13 02994 g011
Figure 12. The DD tropospheric delay interpolation errors obtained by the TLIM and MLIM (left and right) with different satellite elevation angles at different rover stations: (a) P478, (b) EWPP, (c) P603, (d) P225, (e) WINT, and (f) P121.
Figure 12. The DD tropospheric delay interpolation errors obtained by the TLIM and MLIM (left and right) with different satellite elevation angles at different rover stations: (a) P478, (b) EWPP, (c) P603, (d) P225, (e) WINT, and (f) P121.
Remotesensing 13 02994 g012
Figure 13. Elevation angles (upper row) for the reference satellites G10, G18, and G27. The DD tropospheric delay obtained using the TLIM and the PPP method (three lower rows) at different rover stations: (a) P532, (b) P539, (c) P602, and (d) MASW.
Figure 13. Elevation angles (upper row) for the reference satellites G10, G18, and G27. The DD tropospheric delay obtained using the TLIM and the PPP method (three lower rows) at different rover stations: (a) P532, (b) P539, (c) P602, and (d) MASW.
Remotesensing 13 02994 g013
Figure 14. RMS errors of the DD tropospheric delay using modified and traditional methods for all non-reference satellites at all rover stations. The percentages on each bar are the improvements of the MILM with respect to the TLIM.
Figure 14. RMS errors of the DD tropospheric delay using modified and traditional methods for all non-reference satellites at all rover stations. The percentages on each bar are the improvements of the MILM with respect to the TLIM.
Remotesensing 13 02994 g014
Figure 15. DD tropospheric delay interpolation errors obtained by the TLIM and MLIM (left and right) with different satellite elevation angles at different rover stations (from top to bottom): P532, P539, P602, and MASW.
Figure 15. DD tropospheric delay interpolation errors obtained by the TLIM and MLIM (left and right) with different satellite elevation angles at different rover stations (from top to bottom): P532, P539, P602, and MASW.
Remotesensing 13 02994 g015
Figure 16. The TLIM (a) and MLIM (b) positioning errors with respect to the epoch during the selected data period at rover P225.
Figure 16. The TLIM (a) and MLIM (b) positioning errors with respect to the epoch during the selected data period at rover P225.
Remotesensing 13 02994 g016
Figure 17. TLIM (a) and MLIM (b) positioning errors with respect to the epoch during the selected data period at rover MASW.
Figure 17. TLIM (a) and MLIM (b) positioning errors with respect to the epoch during the selected data period at rover MASW.
Remotesensing 13 02994 g017
Table 1. The mean RMS and improvement percentage of the interpolated DD tropospheric delay using the TLIM and MLIM at each rover station in Dataset One.
Table 1. The mean RMS and improvement percentage of the interpolated DD tropospheric delay using the TLIM and MLIM at each rover station in Dataset One.
Station ID P478EWPPP603P225WINTP121
TLIM RMS (m)0.0240.0330.0760.0410.0580.012
MLIM RMS (m)0.0030.0040.0080.0060.0070.004
Improvement (%)87.587.989.585.487.966.7
Table 2. The mean RMS errors and improvement percentage of the interpolated DD tropospheric delay using the TLIM and MLIM at each rover station in Dataset Two.
Table 2. The mean RMS errors and improvement percentage of the interpolated DD tropospheric delay using the TLIM and MLIM at each rover station in Dataset Two.
Station ID P532P539P602MASW
TLIM RMS (m)0.0120.0330.0350.026
MLIM RMS (m)0.0050.0080.0100.007
Improvement (%)58.375.671.473.1
Table 3. Positioning RMS statistical results of each rover station in Dataset One.
Table 3. Positioning RMS statistical results of each rover station in Dataset One.
Accuracy (cm)
Station ID ENU
P478RMSTLIM0.740.9311.59
MLIM0.630.763.42
EWPPRMSTLIM0.911.0512.91
MLIM0.790.892.67
P603RMSTLIM1.351.7624.95
MLIM0.941.054.24
P225RMSTLIM0.851.0812.34
MLIM0.750.791.96
WINTRMSTLIM0.961.2113.36
MLIM0.870.724.32
P121RMSTLIM0.761.074.42
MLIM0.781.072.65
Table 4. Positioning RMS statistical results of each rover station in Dataset Two.
Table 4. Positioning RMS statistical results of each rover station in Dataset Two.
Accuracy (cm)
Station ID ENU
P532RMSTLIM0.740.746.91
MLIM0.770.641.44
P539RMSTLIM0.931.1316.05
MLIM0.830.843.28
P602RMSTLIM0.771.0317.22
MLIM0.670.774.01
MASWRMSTLIM0.800.8514.04
MLIM0.670.762.48
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pu, Y.; Song, M.; Yuan, Y. Modified Interpolation Method of Multi-Reference Station Tropospheric Delay Considering the Influence of Height Difference. Remote Sens. 2021, 13, 2994. https://doi.org/10.3390/rs13152994

AMA Style

Pu Y, Song M, Yuan Y. Modified Interpolation Method of Multi-Reference Station Tropospheric Delay Considering the Influence of Height Difference. Remote Sensing. 2021; 13(15):2994. https://doi.org/10.3390/rs13152994

Chicago/Turabian Style

Pu, Yakun, Min Song, and Yunbin Yuan. 2021. "Modified Interpolation Method of Multi-Reference Station Tropospheric Delay Considering the Influence of Height Difference" Remote Sensing 13, no. 15: 2994. https://doi.org/10.3390/rs13152994

APA Style

Pu, Y., Song, M., & Yuan, Y. (2021). Modified Interpolation Method of Multi-Reference Station Tropospheric Delay Considering the Influence of Height Difference. Remote Sensing, 13(15), 2994. https://doi.org/10.3390/rs13152994

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