Next Article in Journal
Dynamic Remote Sensing Prediction for Wheat Fusarium Head Blight by Combining Host and Habitat Conditions
Previous Article in Journal
Detection of Irrigated and Rainfed Crops in Temperate Areas Using Sentinel-1 and Sentinel-2 Time Series
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method for Mitigating the GPS Multipath Effect Based on a Multi-Point Hemispherical Grid Model

1
GNSS Research Center, Wuhan University, 129 Luoyu Road, Wuhan 430079, China
2
Wuhan Panda Space & Time Technology Co., Ltd., 95 Luoyu Road, Wuhan 430079, China
3
PLA the Rocket Force Command College, 145 Erqi Road, Wuhan 430012, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(18), 3045; https://doi.org/10.3390/rs12183045
Submission received: 19 July 2020 / Revised: 8 September 2020 / Accepted: 10 September 2020 / Published: 18 September 2020

Abstract

:
The multipath effect is a crucial error source caused by the environment around the station and cannot be eliminated or mitigated by differential algorithms. Theoretically, the maximum value for the carrier phase is a quarter the wavelength, i.e., about 4.8 cm for the GPS L1 signal. Considering the increasing demands of high-precision applications, the multipath error has become a major factor affecting the accuracy and reliability of GPS millimeter-level data processing. This paper proposes a multi-point hemispherical grid model (MHGM) to mitigate the multipath effect. In this method, the hemisphere centered on each station is divided into a grid, and the multipath error at the station is estimated based on the parameterization of the grid points. The double-differenced (DD) observed-minus-calculated (OMC) values on some previous days are treated as the observation values to model the present multipath error. Contrary to the present methods which rely much on the platform of data collection and processing, MHGM can be potentially applied to GPS data processing with the existing hardware and software. Experiments in high-multipath and low-multipath environments are designed by mounting a baffle or not. The experimental results show that MHGM is effective in mitigating the multipath effect. When using data from the previous day, an average improvement of about 63.3% in the RMS of DD OMC can be made compared with that without correction, and this is basically consistent with the sidereal filtering (SF) method which is 63.0%. Furthermore, the effectiveness of the above two methods is better than that of the empirical site model (ESM). The kinematic positioning results are also basically consistent with the statistical results of the RMS values of DD OMC. Historical data from more than one day can more explicitly and effectively model the MHGM. Furthermore, compared with the SF method, the MHGM can be used not only to mitigate the multipath error, but also to orientate the sources of the multipath error around the station, and give guidance in the physical elimination of these sources.

1. Introduction

When the direct signals, reflected signals and diffracted signals of global positioning system (GPS) satellites reach the receiver antenna simultaneously, the interference will lead to an important error called the multipath error [1]. For carrier phase observations, the maximum error can be as large as 1/4 the wavelength [2]. It is thus a major factor affecting the accuracy and reliability in GPS data processing [3,4].
Existing methods adopted to mitigate the influence of the multipath error fall into the following categories: (1) selection of an appropriate observation environment and avoidance of specular environmental reflection or radiation sources in the area around the station; (2) adoption of technical means to improve the robustness of reception equipment to the multipath errors. For instance, the narrow correlator uses signal delay to improve the accuracy of the code-tracking loop [5]. Multipath estimation delay lock loop (MEDLL) uses the maximum likelihood criterion to estimate the correlation coefficients of direct and reflected signals to reduce the multipath error [6]. Multipath elimination technology (MET) adopts a new tracking loop to improve upon the narrow correlator technology [7]. The above techniques all carry out some improvements in the correlator to mitigate pseudorange multipath errors. The influence of the reflected signal on the receiver’s antenna can also be mitigated through the use of choke rings [8,9], or the reflected signal can be prevented from reaching the lower part of the antenna through the mounting of antenna ground planes [10]. All these methods acting on the antenna can be applicable for both pseudorange and carrier multipath mitigation; or (3) designing of algorithms in GPS data processing to mitigate the multipath error when the observation environment and receiver equipment cannot be guaranteed. Numerous studies in the literature have made a contribution to this class of approach. A brief introduction of their contributions is as follows, and these methods are mainly used to reduce the carrier phase multipath errors.
Due to the generating mechanism of the multipath effect, signals observed from satellites at low elevation angles are susceptible to more serious multipath errors than those from satellites at higher elevation angles. Thus, the influence of the multipath error on the positioning results can be reduced by using a weighting scheme based on the satellite elevation angle [11,12]. Similarly, signal-to-noise ratio (SNR) information in observed signals can be used to separate reflected or refracted signals from directly transmitted satellite signals [13], or adopted as the basis of a weighting scheme [14]. In addition, many scholars have attempted to use filtering methods to mitigate the multipath error. These attempts have included wavelet analysis [15,16], Kalman filtering [17], use of a Vondrak filter with cross-validation [18], adaptive filtering [19,20], use of a finite impulse response (FIR) filter [21], and so on. These methods can reduce the impact of multipath signals to some extent, but interference from other random noises in the signals limits their ability to eliminate multipath errors. Sidereal filtering (SF) is a widely used method for mitigating the multipath error. This method relies upon the observation that when the spatial relationship between the antenna of the station and the surrounding environment remains unchanged, multipath effects at the station are mainly dependent on the change in the direction of propagation of the satellite signal. The multipath error can thus be identified using repetition of the spatial positional relationship between the satellite and station every sidereal day, allowing correction of observations on subsequent days, and thereby mitigating its effects [22,23,24,25]. However, as the understanding of the generating mechanism of multipath errors has gradually deepened, researchers have discovered that not all GPS satellites have orbital repeat periods of one sidereal day; there are also differences between satellites [26,27]. Choi found through further research that it is more reasonable to use the average orbit period as a basis for filtering: a period that was about 9 s in advance relative to the generally assumed sidereal repeat period for his data. This observation formed the basis for his proposed filtering method—modified sidereal filtering (MSF) [28]. Larson et al. used the time between a given satellite’s arrivals at the same topocentric position as its orbital repeat period, and proposed the improved aspect repeat time adjustment (ARTA) filtering method based on this [29,30]. Wang et al. summarized and compared three different methods for estimating the GPS multipath repeat time [31]. All of the methods described above conduct SF in the coordinate domain, that is, the multipath error is extracted from the coordinate sequence for correction of positioning results on subsequent sidereal days. Corrections based on these sidereal filtering methods can also be applied to the observation domain. Alber (2000), during research into the moisture content of the atmosphere, proposed a method to isolate the multipath effect by transforming double-differenced residuals [32]. Ragheb (2007) compared SF methods based on the coordinate domain with those based on the observation domain, concluding that the former was more computationally efficient, but the latter was more accurate [33]. Zhong (2010), using SF on a short baseline of a single-differenced observation, was able to obtain more accurate positioning results than those possible with corrections applied in the coordinate domain [34].
As with modeling the multipath error in the coordinate or observation domain, based on the repetitive characteristics of satellite orbits, modeling the error in the spatial domain can also be an effective solution for mitigation of its influence. Cohen and Parkinson (1991) used observed data to set up a lookup table mapping the multipath environment around the receiver antenna [35]. Fuhrmann (2015) pointed out that it is unreasonable that stacking cells have a fixed azimuthal resolution, so an advanced method that makes use of congruent cells which have a similar shape and size was proposed [36]. Further, in its realization, the precise point positioning (PPP) residuals are used to generate the multipath stacking maps, which may add satellite clock errors and orbital errors in the final model. Moore (2014) used an empirical site model (ESM) to mitigate unmodeled site-specific errors [37]. Zero-differenced (ZD) carrier phase residuals are used to generate a site-specific correction map in this method. However, the technique to transform double-differenced residuals into ZD residuals may be not suitable due to the assumption of “zero mean” [38]. Dong (2016) proposed a multipath error correction method based on a multipath hemispherical map [39]. This method averages the residual of observed values corresponding to the fixed solution of single-difference ambiguity between stations of the satellite when a common-view satellite passes through a grid on the hemisphere between stations. Then, the low frequency multipath effect is retained and can be modeled while the noise effect is weakened. It solves the disadvantages due to the assumption of “zero mean”, but its available distance is limited by a common receiver clock with multiple receivers.
In this work, we first give a brief presentation of the characteristics of the multipath effect in the space domain. Based on that, a novel method called the multi-point hemispherical grid model (MHGM) for mitigating the influence of the multipath error is introduced. The new method divides the hemisphere centered on each station into a grid, and estimates the multipath error at the station based on the parameterization of the grid points. Conventional double-differenced observations are used for modeling, so it can be potentially applied to GPS data processing with the existing hardware and software. Afterward, experiments in different multipath environments are designed, and observation residuals and positioning results before and after correction using MHGM are analyzed and compared with the results of the SF method and ESM. Finally, we draw some conclusions and give some suggestions of potential applications of MHGM in the future.

2. The Multi-Point Hemispherical Grid Model

In order to model the carrier phase multipath in the spatial domain, it is necessary to understand the characteristic of the multipath error. A model for the effects of the carrier phase multipath is built. Assuming that the receiver antenna is in a single planar horizontal reflector which is infinitely large and that its transmissions are only affected by a single reflected signal, the error M caused by the multipath effect on the carrier phase observation can be described by the following formula [40]:
M = λ 2 π t a n 1 α s i n ( 4 π H λ s i n ε ) 1 + α c o s ( 4 π H λ s i n ε )
where H is the antenna height, ε is the angle of incidence of the reflected signal, α is the reflection coefficient of the reflected signal, and λ is the wavelength of the satellite signal.
It can be seen from (1) that for a satellite signal of a certain wavelength, the multipath error depends upon the reflection coefficient, the angle of incidence of the reflected signal, and the spatial relationship between the antenna and the reflector. If the spatial relationship between the antenna and the surrounding environment remains unchanged during the observation process, changes in the propagation path of the reflected signal will mainly be due to changes in the direction of propagation of the satellite signal; its angle of incidence and reflection coefficient will change accordingly, ultimately affecting the multipath error present in the observed signal.
Based on the characteristics of the multipath effect in the space domain described above, a hemisphere can be defined centered upon the antenna’s position, allowing the position of points to be described by a pair of angular coordinates (i.e., elevation and azimuth angles). Assuming an insignificant change in the environment around the station, and a given signal frequency, for a point in the hemispherical coordinate system based on the antenna, the multipath effect is related only to the elevation and azimuth angles of the satellite, and is not dependent on the signal’s time of observation, or the satellite’s pseudo-random noise (PRN) code.
This paper thus proposes the multi-point hemispherical grid model. This method aims to isolate multipath information, using inter-station double-differenced (DD) observed-minus-calculated (OMC), allowing the establishment of a multi-point hemispherical grid model consistent across stations. The lower hardware requirements of the novel approach allow it to be potentially applied in existing GPS network data processing systems. In particular, the use of DD OMC to model the multipath error does not require specialized receiver equipment, and the distance between stations is not limited by hardware considerations.
A grid is constructed for each station for the MHGM. For the elevation angle, the minimum value for grid points is set to   E 0 , and the maximum value is set to   E 1 . The azimuth angle of grid points is set to 0° to 360°, the range of possible values for it. A grid point is also set at the zenith of the coordinate system. Grid points are separated by an interval of d e for the elevation angle and d a for the azimuth angle. The number of parameters to be estimated in the model is determined by the value of d e and   d a . The more parameters to be estimated, the more memory resources of the computer are consumed and the longer computation time is taken, but the description of the multipath error at the station is more detailed. The above factors should be taken into account to determine the value of d e and d a in modeling. An example of how the grid could be divided for a set of stations is given below ( E 0   = 0°, E 1 = 60°, d e = 30°, d a = 30°); the corresponding model parameters must be estimated for each of the black grid points shown.
The Positioning and Navigation Data Analyst (PANDA) software was used for attaining DD OMC between satellites and stations [41]. PANDA works under the ZD processing mode and the zero-differenced observed-minus-calculated (OMC) between each station and each visible satellite can be attained directly [42]. After that, the ZD ambiguities are mapped to a maximum set of independent DD ambiguities. Finally, the DD ambiguities could be solved and fixed sequentially [43]. For each set of fixed DD ambiguities processed by the PANDA software, the four groups of ZD OMC in the corresponding period could be obtained, so that the DD OMC of the carrier phase observation corresponding to the DD ambiguity is calculated according to the following equation as
D D   O M C = ( Z D m j Z D m k ) ( Z D n j Z D n k ) Δ N
The Δ N in Equation (2) stands for the DD ambiguities. Normal equations corresponding to the grid points to be evaluated are established for the DD OMC s between the stations ( m , n ) and the satellites ( j , k ), corresponding to an arbitrary set of fixed-ambiguity solutions. Figure 1 shows that the multipath error between the station m and the satellite j is assumed to involve four points Q m 1 j , Q m 2 j , Q m 3 j and Q m 4 j (where the elevation angle of satellite j is between E 0 and E 1 ), or three points Q m 1 j , Q m 2 j and Q m 3 j (where the elevation angle of satellite j exceeds E 1 ), for this grid model.
If the spatial relationship between the pierce point of the propagation of the signal of satellite j and the hemispherical grid point parameters is as shown in the diagram on the upper right side of Figure 1, and it is assumed that the elevation angles and azimuth angles corresponding to the four grid point parameters Q m 1 j , Q m 2 j , Q m 3 j and Q m 4 j are respectively ( e 1 , a 1 ), ( e 1 , a 2 ), ( e 2 , a 2 ) and ( e 2 , a 1 ), then estimates of the value of model parameters at V m j can be calculated by bilinear interpolation, as follows:
V m j = w m 1 j Q m 1 j + w m 2 j Q m 2 j + w m 3 j Q m 3 j + w m 4 j Q m 4 j
where w m 1 j = ( e 2 e ) ( a 2 a ) ( e 2 e 1 ) ( a 2 a 1 ) , w m 2 j = ( e 2 e ) ( a a 1 ) ( e 2 e 1 ) ( a 2 a 1 ) , w m 3 j = ( e e 1 ) ( a a 1 ) ( e 2 e 1 ) ( a 2 a 1 ) , w m 4 j = ( e e 1 ) ( a 2 a ) ( e 2 e 1 ) ( a 2 a 1 ) .
If the spatial relationship between the pierce point of the propagation of the signal of satellite j , and the hemispherical grid point parameters is as shown in the diagram on the lower right side of Figure 1, and assuming that corresponding elevation angles and azimuth angles of the three grid point parameters Q m 1 j , Q m 2 j and Q m 3 j are respectively ( e 1 , a 1 ), ( e 2 , a 2 ) and ( e 2 , a 3 ), then estimates of the value of model parameters at V m j can be calculated by plane interpolation, as follows:
V m j = ( 1 u v ) Q m 1 j + u Q m 2 j + v Q m 3 j
where u   = ( a a 1 ) ( e 2 e 1 ) ( e e 1 ) ( a 3 a 1 ) ( e 2 e 1 ) ( a 2 a 3 ) , v = ( a 2 a 1 ) ( e e 1 ) ( e 2 e 1 ) ( a a 1 ) ( e 2 e 1 ) ( a 2 a 3 ) .
Similarly, expressions for the values V m k , V n j and V n k of model parameters, between station m and satellite k , station n and satellite j and station n and satellite k , can be established. Then, the DD OMC s also can be expressed as follows:
s = ( V m j V m k ) ( V n j V n k )
Observation (5) is used to derive normal equations corresponding to the relevant grid point parameters, and the fused solution, for several different pairs of satellites and ground stations and for data observed during different periods, is generated through superposition of these normal equations. The accuracy and reliability of the model’s solution can also be further enhanced by imposing constraints on the correlation among grid point parameters, and imposing a requirement that the absolute values of grid point parameters must be less than 1/4 the wavelength. The specific process of the model parameter estimation is displayed in Figure 2.

3. Experiments and Discussion

3.1. Experimental Design

Experimental verification was carried out on the top floor of the academic experiment building at Wuhan University during days 226–233 in 2018. Three GPS stations, A, B and C, were set up, with a metal baffle mounted to the west of station A to create a simulated high-multipath environment. Stations B and C were positioned in normal observation environments. The hardware specifications of stations A, B and C used in the experiment are shown in Table 1. The antenna used in each station was a common survey antenna without any anti-multipath techniques. The distances between the three stations were 2.7, 4.7 and 6.9 m, respectively, and their spatial relationship with the surrounding observation environments is shown in Figure 3.
Using the accurately known coordinates of the three stations, multipath error modeling was applied to data from L1 observations on days 226–233 in 2018. The effectiveness of multipath error correction for the observations from the 233rd day of 2018 was analyzed. The models and parameters used in the data processing of the PANDA software are listed in Table 2. The maximum DD multipath error should not exceed 1/2 the wavelength when the multipath errors reach the theoretical maximum and the signs are opposite at both ends of the baseline. For this reason, the DD OMC larger than 1/2 the wavelength will be eliminated in the data preprocessing before the process of modeling.
To demonstrate the influence of the duration of observation data available on the effectiveness of multipath error modeling, MHGMs were separately estimated based on 1–7 days of data during the experiment. To keep a balance of the computational efficiency and the level of details of the model, the model parameters were set as follows: E 0 = 5°, E 1 = 85°, d e = 2° and d a = 2°. In the process of modeling the multipath error, there are three kinds of observation equations.
(1)
The observation (3) constructed by DD OMC observations, where the weights of this kind of observation equation are determined by the accuracy of DD OMC observations, whose accuracy can achieve a millimeter or better [44].
(2)
The observation equation constructed by imposing a requirement on the absolute values of the grid point parameters, where the weights of this kind of observation equation are determined by 1/4 of the L1 wavelength [2].
(3)
The observation equation constructed by imposing constraints on the correlation among grid point parameters, where the weights of this kind of observation equation are determined by an empirical variation value between grid point parameters, and this value is set to 1.0 cm/degree based on the results of multiple experiments, and it is found that setting the value to 5 mm/degree ~5 cm/degree has little effect on the final model.
For each day, the same DD ambiguities fixed results are used to establish MHGM, ESM and SF models. Among them, the MHGM also has been modeled using multi-day historical data. Therefore, this paper provides the results with different kinds of multipath error correction and the result without correction for comparison. The specific solution strategies adopted in the experiment (i.e., combinations of methods used and datasets to which they were applied) are shown in Table 3.

3.2. Analysis of Results

PANDA works under the ZD processing mode and the zero-differenced observed-minus-calculated between each station and each visible satellite can be attained directly. Corresponding to the strategies shown in Table 3, Figure 4 shows the ZD OMC for each satellite at the three stations A, B and C during a period of fixed ambiguity, in the form of a sub-satellite point track. With no correction for the multipath error (N), the results for all three stations, A (N), B (N) and C (N), include relatively large ZD OMCs. This is particularly noticeable for paths via the west side of station A where the baffle was mounted. So, the multipath effect is particularly obvious on the left side of diagram A (N) with many groups of large observation residuals. For the SF method of multipath error correction, the satellites’ orbit repeat period advance must be taken into account. In this experiment, the SF in the observation domain is adopted, in which the average orbit repeat period of the whole GPS constellation is set to 23 h 55 m 55 s (86,155 s), with a daily advance of 245 s [28]. The results of A (SF), B (SF) and C (SF) show that the multipath effect has been effectively mitigated most of the time, but the residuals towards the end of the satellite trajectories are larger and it is a sequence that cannot be corrected effectively using the SF model. Using the MHGM proposed in this paper and the ESM, even only using the data of the previous day for modeling, as shown in A (M1/E1), B (M1/E1) and C (M1/E1), these uncorrected satellite trajectories are less obvious (i.e., the phenomenon of larger residuals towards the end of the satellite observation periods does not exist). The modeling effectiveness of MHGM can be further improved by using data observed during the prior two or more days. This is because the orbit repeat period of GPS satellites is about one day, and the multipath error in GPS observation data can be effectively modeled using data observed over the most recent two days. In order to verify the validity of MHGMs constructed by using multi-day GPS data, the M7 strategy is also shown in the following conclusions. Figure 4, abbreviated for reasons of space, shows only the ZD OMC at the three stations for the strategies of N, SF, ESM, M1 and M7.
To quantify the effects of the MHGM proposed in this paper in comparison with those of the SF and ESM methods, and to evaluate differences in the effectiveness of multipath error modeling based on data observed over 1~7 days, Table 4 shows the mean root mean square (RMS) of DD OMC during fixed-ambiguity periods on day 233 for all 10 strategies.
The HMP in Table 4 stands for the statistical results of baseline A–C which is in a high-multipath environment. If data observed on the previous day only are used for multipath error modeling, the effect of the ESM method is very obvious. The DD RMS improvement of strategy ESM compared with strategy N can reach 68.3%, but it is worse than that of SF. For strategy M1, it is comparable to the improvement level of strategy SF, and the improvement increases from 72.6% to 75.4%. For multipath error modeling based on data observed over the previous two days or more, the improvements offered by M2–7 average about 75.6%. The improvement between M2 and M7 is not very obvious, that is because MHGM can make full use of the correlation of the multipath error in the neighboring directions, and then the direction without valid modeling data also can be appropriately corrected by extrapolation or interpolation.
The LMP in Table 4 stands for the statistical results of baseline B–C which is in a low-multipath environment. Strategy M1 is slightly worse than that of strategy SF, but the results of multiple days can further improve the mitigation effect of MHGM. Figure 5 provides a graphical display of the distribution of the RMS values of the observation residuals from the double-difference ambiguity fixed arcs from the 233rd day, derived using the different processing strategies.
As can be seen in Figure 6, the low-frequency patterns are basically eliminated after the correction using different models in the double-differenced residual series of G01 and G03 on station A and station B. Compared with SF, the methods based on the spatial domain can also mitigate the low-frequency part of the multipath effects, and the series are slightly more stable with model correction using MHGM. Power spectrum analysis indicates that the low-frequency multipath effects can be reduced with these multipath corrections. In the frequency band higher than 0.1 Hz, however, the power of the MHGM is similar to the ESM and SF.
Precise positioning is a more intuitive form to assess the effect of multipath error correction. For this reason, the influence of multipath error correction on the positioning results using SF, ESM and MHGM methods is processed and shown in Table 5. Station C is fixed in the data processing, and the three-dimensional coordinates of static stations A and B are solved according to the kinematic relative positioning mode. For the new method, the strategies M1 and M7 are used to model the multipath error.
When the multipath error is uncorrected, the positioning accuracy of station A is worse than that of station B because a metal baffle is mounted to station A to simulate the high-multipath environment. The horizontal positioning accuracy of stations A and B is 5.48 and 3.35 mm, respectively, and the vertical positioning accuracy is 9.90 and 5.15 mm, respectively. The correction effect achieved by the MHGM can be basically consistent with the SF method only based on data from the previous day, and the improvement of the three-dimensional positioning accuracy is about 59.7% and 60.9%. Further, the positioning results using strategy ESM are worse than M1 and SF. When seven days of observations are used for modeling the multipath error using MHGM, the positioning results of M7 can improve 10.2% in three directions compared with those of M1. The kinematic positioning results are basically consistent with the statistical results of the RMS values of DD OMC. In order to further demonstrate the influence of five processing strategies, N, SF, ESM, M1 and M7, on the positioning results at different times, the real-time kinematic positioning error sequences and the corresponding RMSs of station A (HMP) in the north, east and up directions are given in Figure 7.
Compared with the SF method, the MHGM proposed in this paper not only enhances the improvements obtained through multipath error correction, but also effectively provides a graphical display of multipath error interference around the stations. Figure 8 shows the hemispherical modeling results from strategy M1 at the stations. In this, “overlook”, “west” and “east” are visualizations of the effect from above, looking to the west side and looking to the east side, respectively. The top view shows higher modeled values in large areas of the left side area of A (M1). Further, in the west side view, an obvious and continuous abnormal area at the bottom is visible, the location of which coincides with the position where the metal baffles were mounted: this reflects the influence of these baffles on the multipath error at station A. This characteristic of the MHGM implies it can not only be used to mitigate the influence of the multipath error at stations, but also to orientate the sources of the multipath error around the station, providing guidance as to the measures which could be useful for physically eliminating the influence of these error sources.
In order to show the changes in the model results over multiple days intuitively, the modeling results of the three stations based on daily GPS observations from 1 day to 7 days are also attached to the Supplementary Materials (Supplemental-Gif-1.mp4). Due to space considerations, only the modeling results of MHGM using seven daily observations are shown below. It can be seen that the modeling results of MHGM from strategy M7 are generally similar to those of strategy M1, but the results of strategy M7 have more details in the top view of Figure 9. As previously mentioned, the orbit repeat period of GPS satellites is about one day, so multipath error modeling based on one-day observations will be feasible, and modeling base on multi-day observations can further make up for the lack of data that may exist in one-day observations.

4. Conclusions

This paper proposes a novel method for the mitigation of multipath error at stations based on a multi-point hemispherical grid model. The method divides the hemisphere centered on each station into meshes, and estimates the multipath error at the station based on the parameterization of the grid points. During the process of grid point parameter estimation, the multi-point hemispherical grid model for each station is obtained, taking the DD OMC corresponding to the solutions of fixed ambiguity for different satellite pairs and for different periods as input information.
The experimental results show that for the multipath error modeled by data from only the previous day, the correction performance of the MHGM can be generally consistent with the SF method, and better than ESM. When the previous two or more days’ observation data are used for modeling the multipath error, the MHGM can be described in more detail. The average improvement of the mean RMS of DD OMC increases from 63.3% (M1) to 66.3% (M7), in comparison to the uncorrected results. Further, the kinematic positioning results also coincide with the statistical RMS values of DD OMC. The results also show that the multipath error modeling from one-day observations will be feasible and modeling from multi-day observations can further make up the lack of data that may exist in one-day observations.
Compared with the SF method, the proposed MHGM can not only mitigate the multipath error effectively, but also provide a graphical distribution of multipath error interference around the stations, and thus give guidance in physically eliminating the sources of the multipath error in the future.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-4292/12/18/3045/s1, Supplemental-Gif-1.mp4.

Author Contributions

Y.W., X.Z. and W.T. provided the initial idea and designed the research; X.Z. and Y.W. processed data, and wrote the manuscript; Y.L. helped to accomplish some tests; C.D., Y.Z. and J.F. helped with the writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (Grant No. 2018YFB0505201), National Natural Science Foundation of China (Grant Nos. 41874037, 41804028), Fundamental Research Funds for the Central Universities (Grant No. 2042018gf0001) and the Wuhan Science and Technology Project (Grant No. 2018010401011271).

Acknowledgments

The authors would like to acknowledge the IGS for providing the broadcast ephemeris and precise products.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Roberts, G.W.; Meng, X.; Psimoulis, P.; Brown, C.J. Time series analysis of rapid GNSS measurements for quasi-static and dynamic bridge monitoring. In Geodetic Time Series Analysis in Earth Sciences; Montillet, J.P., Bos, M., Eds.; Springer Geophysics: Cham, Switzerland, 2020; pp. 345–417. [Google Scholar]
  2. Hofmann-Wellenhof, B.; Lichtenegger, H.; Wasle, E. GNSS–Global Navigation Satellite Systems: GPS, Glonass, Galileo, and More; Springer: Vienna, Austria, 2008. [Google Scholar]
  3. Xie, L.; Cui, X.; Zhao, S.; Lu, M. Mitigating multipath bias using a dual-polarization antenna: Theoretical performance, algorithm design, and simulation. Sensors 2017, 17, 359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Chang, G.; Chen, C.; Yang, Y.; Xu, T. Tikhonov regularization based modeling and sidereal filtering mitigation of GNSS multipath errors. Remote Sens. 2018, 10, 1801. [Google Scholar] [CrossRef] [Green Version]
  5. Dierendonck, A.J.v.; Fenton, P.; Ford, T. Theory and performance of narrow correlator spacing in a GPS receiver. Navigation 1992, 39, 265–283. [Google Scholar] [CrossRef]
  6. Nee, R.D.J.v. The multipath estimating delay lock loop. In Proceedings of the IEEE 2nd International Symposium on Spread Spectrum Techniques and Applications, Yokohama, Japan, 29 November–2 December 1992; pp. 39–42. [Google Scholar]
  7. Townsend, B.R.; Fenton, P. A practical approach to the reduction of pseudorange multipath errors in a L1 GPS receiver. In Proceedings of the 7th International Technical Meeting of the Satellite Division of the Institute of Navigation, Salt Lake City, UT, USA, 20–23 September 1994. [Google Scholar]
  8. Filippov, V.; Tatarnicov, D.; Ashjaee, J.; Astakhov, A.; Sutiagin, I. The first dual-depth dual-frequency choke ring. 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; pp. 1035–1040. [Google Scholar]
  9. Tranquilla, J.M.; Carr, J.P.; Al-Rizzo, H.M. Analysis of a choke ring ground plane for multipath control in global positioning system (GPS) applications. IEEE Trans. Antennas Propag. 1994, 42, 905–911. [Google Scholar] [CrossRef]
  10. Tatarnikov, D.; Filippov, V.; Soutiaguine, I.; Astahov, A.; Stepanenko, A.; Shamatulsky, P. Multipath mitigation by conventional antennas with ground planes and passive vertical structures. GPS Solut. 2005, 9, 194–201. [Google Scholar] [CrossRef]
  11. Rothacher, M.; Springer, T.; Schaer, S.; Beutler, G. Processing strategies for regional GPS networks. In Advances in Positioning and Reference Frames; Springer: Heidelberg, Germany, 1998; pp. 93–100. [Google Scholar]
  12. Meindl, M.; Schaer, S.; Hugentobler, U.; Beutler, G. Tropospheric gradient estimation at CODE: Results from global solutions. J. Meteorol. Soc. Japan. Ser. II 2004, 82, 331–338. [Google Scholar] [CrossRef] [Green Version]
  13. Comp, C.; Axelrad, P. Adaptive SNR based carrier phase multipath mitigation technique. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 264–276. [Google Scholar] [CrossRef]
  14. Lau, L.; Mok, E. Improvement of GPS relative positioning accuracy by using SNR. J. Surv. Eng. 1999, 125, 185–202. [Google Scholar] [CrossRef]
  15. Souza, E.M.; Monico, J.F.G. Wavelet shrinkage: High frequency multipath reduction from GPS relative positioning. GPS Solut. 2004, 8, 152–159. [Google Scholar] [CrossRef]
  16. Satirapod, C.; Rizos, C. Multipath mitigation by wavelet analysis for GPS base station applications. Surv. Rev. 2005, 38, 2–10. [Google Scholar] [CrossRef]
  17. Ince, C.D.; Sahin, M. Real-time deformation monitoring with GPS & Kalman filter. Earth Planets Space 2000, 52, 837–840. [Google Scholar]
  18. Zheng, D.; Zhong, P.; Ding, X.; Chen, W. Filtering GPS time-series using a Vondrak filter and cross-validation. J. Geod. 2005, 79, 363–369. [Google Scholar] [CrossRef]
  19. Ge, L.; Han, S.; Rizos, C. Multipath mitigation of continuous GPS measurements using an adaptive filter. GPS Solut. 2000, 4, 19–30. [Google Scholar] [CrossRef]
  20. Yu, S.; Guo, F.; Zhang, X.; Liu, W.; Li, X.; Wu, R. A new method for GNSS multipath mitigation with an adaptive frequency domain filter. Sensors 2018, 18, 2514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Han, S.; Rizos, C. Multipath effects on GPS in mine environments. In Proceedings of the 10th International Congress of the International Society for Mine Surveying, Fremantle, Australia, 2–6 November 1997; pp. 447–457. [Google Scholar]
  22. Genrich, J.; Bock, Y. Rapid resolution of crustal motion at short ranges with the global positioning system. J. Geophys. Res. 1992, 97, 3261–3269. [Google Scholar] [CrossRef]
  23. Bock, Y.; Nikolaidis, R.M.; Jonge, P.J.; Michael, B. Instantaneous geodetic positioning at medium distances with the global positioning system. J. Geophys. Res. 2000, 105, 28223–28254. [Google Scholar] [CrossRef]
  24. Bock, Y.; Prawirodirdjo, L.; Melbourne, T.I. Detection of arbitrarily large dynamic ground motions with a dense high-rate GPS network. Geophys. Res. Lett. 2004, 31, L06604. [Google Scholar] [CrossRef] [Green Version]
  25. Nikolaidis, R.M.; Bock, Y.; Jonge, P.J.D.; Shearer, P.; Agnew, D.C.; Domselaar, M.V. Seismic wave observations with the global positioning system. J. Geophys. Res. 2001, 106, 21897–21916. [Google Scholar] [CrossRef] [Green Version]
  26. Seeber, G.; Menge, F.; Völksen, C.; Wübbena, G.; Schmitz, M. Precise GPS Positioning Improvements by Reducing Antenna and Site Dependent Effects; Springer: Heidelberg, Germany, 1998; Volume 118. [Google Scholar]
  27. Ding, X.; Chen, Y.; Zhu, J.; Huang, D. Surface deformation detection using GPS. In Proceedings of the 12th International Technical Meeting of the Satellite Division of The Institute of Navigation, Nashville, TN, USA, 14–17 September 1999; pp. 53–62. [Google Scholar]
  28. Choi, K.; Bilich, A.; Larson, K.M.; Axelrad, P. Modified sidereal filtering: Implications for high-rate GPS positioning. Geophys. Res. Lett. 2004, 31, 178–198. [Google Scholar] [CrossRef] [Green Version]
  29. Agnew, D.C.; Larson, K.M. Finding the repeat times of the GPS constellation. GPS Solut. 2007, 11, 71–76. [Google Scholar] [CrossRef]
  30. Larson, K.; Bilich, A.; Axelrad, P. Improving the precision of high-rate GPS. J. Geophys. Res. 2007, 112, B05422. [Google Scholar] [CrossRef]
  31. Wang, M.; Wang, J.; Dong, D.; Li, H.; Han, L.; Chen, W. Comparison of three methods for estimating GPS multipath repeat time. Remote Sens. 2018, 10, 6. [Google Scholar] [CrossRef] [Green Version]
  32. Alber, C.; Ware, R.; Rocken, C.; Braun, J. Obtaining single path phase delays from GPS double differences. Geophys. Res. Lett. 2000, 27, 2661–2664. [Google Scholar] [CrossRef] [Green Version]
  33. Ragheb, A.E.; Clarke, P.J.; Edwards, S.J. GPS sidereal filtering: Coordinate- and carrier-phase-level strategies. J. Geod. 2007, 81, 325–335. [Google Scholar] [CrossRef]
  34. Zhong, P.; Ding, X.; Yuan, L.; Xu, Y.; Kwok, K.; Chen, Y. Sidereal filtering based on single differences for mitigating GPS multipath effects on short baselines. J. Geod. 2010, 84, 145–158. [Google Scholar] [CrossRef]
  35. Cohen, C.E.; Parkinson, B.W. Mitigating multipath error in GPS-based attitude determination. In Proceedings of the Annual Rocky Mountain Guidance and Control Conference, Keystone, CO, USA, 2–6 February 1991; pp. 53–68. [Google Scholar]
  36. Fuhrmann, T.; Luo, X.; Knoepfler, A.; Mayer, M. Generating statistically robust multipath stacking maps using congruent cells. GPS Solut. 2015, 19, 83–92. [Google Scholar] [CrossRef]
  37. Moore, M.; Watson, C.; King, M.; McClusky, M.; Tregoning, P. Empirical modelling of site-specific errors in continuous GPS data. J. Geod. 2014, 88, 887–900. [Google Scholar] [CrossRef]
  38. Shi, Q.; Dai, W.; Zeng, F.; Kuang, C. The BDS Multipath Hemispherical Map Based on Double Difference Residuals and Its Application Analysis. In China Satellite Navigation Conference (CSNC) 2016 Proceedings: Volume I; Lecture Notes in Electrical Engineering; Sun, J., Liu, J., Fan, S., Wang, F., Eds.; Springer: Singapore, 2016; Volume 388. [Google Scholar]
  39. Dong, D.; Wang, M.; Chen, W.; Zeng, Z.; Song, L.; Zhang, Q.; Cai, M.; Cheng, Y.; Lv, J. Mitigation of multipath effect in GNSS short baseline positioning by the multipath hemispherical map. J. Geod. 2016, 90, 255–262. [Google Scholar] [CrossRef]
  40. Elósegui, P.; Davis, J.L.; Jaldehag, R.T.K.; Johansson, J.M.; Niell, A.E.; Shapiro, I.I. Geodesy using the global positioning system: The effects of signal scattering on estimates of site position. J. Geophys. Res. 1995, 100, 9921–9934. [Google Scholar] [CrossRef]
  41. Liu, J.; Ge, M. PANDA software and its preliminary result of positioning and orbit determination. Wuhan Univ. J. Nat. Sci. 2003, 8, 603–609. [Google Scholar]
  42. Shi, C.; Zhao, Q.; Geng, J.; Lou, Y.; Ge, M.; Liu, J. Recent development of PANDA software in GNSS data processing. In International Conference on Earth Observation Data Processing and Analysis (ICEODPA); International Society for Optics and Photonics: Bellingham, WA, USA, 2008; p. 7285. [Google Scholar]
  43. Ge, M.; Gendt, G.; Dick, G.; Zhang, F.P. Improving carrier-phase ambiguity resolution in global GPS network solutions. J. Geod. 2005, 79, 103–110. [Google Scholar] [CrossRef]
  44. Bona, P. Accuracy of GPS phase and code observations in practice. Acta Geod. Geoph. Hung. 2000, 35, 433–451. [Google Scholar]
Figure 1. Distribution of hemispherical grid points.
Figure 1. Distribution of hemispherical grid points.
Remotesensing 12 03045 g001
Figure 2. Process of model parameter estimation.
Figure 2. Process of model parameter estimation.
Remotesensing 12 03045 g002
Figure 3. Spatial relationship between stations and surrounding observation environment ((A) is for station A, (B) is for station B, (C) is for station C).
Figure 3. Spatial relationship between stations and surrounding observation environment ((A) is for station A, (B) is for station B, (C) is for station C).
Remotesensing 12 03045 g003
Figure 4. Satellite zero-differenced (ZD) observed-minus-calculated (OMC) at stations (A), (B) and (C).
Figure 4. Satellite zero-differenced (ZD) observed-minus-calculated (OMC) at stations (A), (B) and (C).
Remotesensing 12 03045 g004
Figure 5. Distribution of RMS DD OMC for a fixed-ambiguity period under different processing strategies.
Figure 5. Distribution of RMS DD OMC for a fixed-ambiguity period under different processing strategies.
Remotesensing 12 03045 g005
Figure 6. The double-differenced residual series and power spectral density in the period of fixed ambiguity, 233th day, 2018.
Figure 6. The double-differenced residual series and power spectral density in the period of fixed ambiguity, 233th day, 2018.
Remotesensing 12 03045 g006
Figure 7. Real-time kinematic positioning error sequences and the corresponding RMSs of station A.
Figure 7. Real-time kinematic positioning error sequences and the corresponding RMSs of station A.
Remotesensing 12 03045 g007
Figure 8. Hemispherical modeling results from strategy M1 ((A) is for the MHGM from strategy M1 of station A, (B) is for the MHGM from strategy M1 of station B, (C) is for the MHGM from strategy M1 of station C).
Figure 8. Hemispherical modeling results from strategy M1 ((A) is for the MHGM from strategy M1 of station A, (B) is for the MHGM from strategy M1 of station B, (C) is for the MHGM from strategy M1 of station C).
Remotesensing 12 03045 g008
Figure 9. Hemispherical modeling results from strategy M7 ((A) is for the MHGM from strategy M7 of station A, (B) is for the MHGM from strategy M7 of station B, (C) is for the MHGM from strategy M7 of station C).
Figure 9. Hemispherical modeling results from strategy M7 ((A) is for the MHGM from strategy M7 of station A, (B) is for the MHGM from strategy M7 of station B, (C) is for the MHGM from strategy M7 of station C).
Remotesensing 12 03045 g009
Table 1. The hardware specifications of stations A, B and C.
Table 1. The hardware specifications of stations A, B and C.
StationABC
ReceiverSeptentrio PolaRx5SPANDA PD318Septentrio PolaRx5S
AntennaHarxon HX-CSX607AHarxon HX-CSX607AHarxon HX-CSX607A
with baffle or notYesNoNo
Table 2. Model and parameters for data processing.
Table 2. Model and parameters for data processing.
ParameterModelConstraint
ObservationL1 + P10.02 cycle, 1.0 m
Observation weightSatellite elevation/
Cutoff elevation7 degree/
Phase center patternigs14.atx/
Tropospheric delaySaastamoinen + GMFpriori covariance 0.300 m
Satellite clockBroadcast + Processpriori covariance 3000 m
Receiver clockRange estimating + White noisepriori covariance 9000 m
EOPFixed to IERS/
Satellite orbitBroadcastpriori covariance 10 m
Station displacementSolid earth, pole tide, ocean loading/
Table 3. Solution strategies for multipath error analysis.
Table 3. Solution strategies for multipath error analysis.
StrategyModeling MethodDays of the Year 2018
Modeling DataTest Data
N//233
SFSF232233
ESMESM232233
M1MHGM232233
M2MHGM231~232233
M3MHGM230~232233
M4MHGM229~232233
M5MHGM228~232233
M6MHGM227~232233
M7MHGM226~232233
Table 4. Mean RMS values of double-differenced (DD) OMC for a fixed-ambiguity period under different processing strategies.
Table 4. Mean RMS values of double-differenced (DD) OMC for a fixed-ambiguity period under different processing strategies.
Strategy Mean RMS of DD OMC/cm Improvement Compared with Strategy N
HMPLMPHMPLMP
N1.055 0.629 //
SF0.289 0.294 72.6%53.4%
ESM0.3350.34568.3%45.2%
M10.260 0.307 75.4%51.2%
M20.257 0.281 75.6%55.3%
M30.257 0.282 75.6%55.2%
M40.258 0.273 75.5%56.7%
M50.257 0.270 75.7%57.0%
M60.257 0.270 75.7%57.0%
M70.257 0.271 75.6%57.0%
Table 5. RMS of kinematic positioning results.
Table 5. RMS of kinematic positioning results.
Strategy RMS (mm)
A (HMP)B (LMP)
HorizontalVerticalHorizontalVertical
N5.48 9.90 3.35 5.15
SF1.543.421.522.31
ESM2.114.031.953.07
M11.522.711.732.77
M71.50 2.46 1.54 2.37

Share and Cite

MDPI and ACS Style

Wang, Y.; Zou, X.; Deng, C.; Tang, W.; Li, Y.; Zhang, Y.; Feng, J. A Novel Method for Mitigating the GPS Multipath Effect Based on a Multi-Point Hemispherical Grid Model. Remote Sens. 2020, 12, 3045. https://doi.org/10.3390/rs12183045

AMA Style

Wang Y, Zou X, Deng C, Tang W, Li Y, Zhang Y, Feng J. A Novel Method for Mitigating the GPS Multipath Effect Based on a Multi-Point Hemispherical Grid Model. Remote Sensing. 2020; 12(18):3045. https://doi.org/10.3390/rs12183045

Chicago/Turabian Style

Wang, Yawei, Xuan Zou, Chenlong Deng, Weiming Tang, Yangyang Li, Yongfeng Zhang, and Jin Feng. 2020. "A Novel Method for Mitigating the GPS Multipath Effect Based on a Multi-Point Hemispherical Grid Model" Remote Sensing 12, no. 18: 3045. https://doi.org/10.3390/rs12183045

APA Style

Wang, Y., Zou, X., Deng, C., Tang, W., Li, Y., Zhang, Y., & Feng, J. (2020). A Novel Method for Mitigating the GPS Multipath Effect Based on a Multi-Point Hemispherical Grid Model. Remote Sensing, 12(18), 3045. https://doi.org/10.3390/rs12183045

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