Positioning | |
Practical multivariate statistical multipath detection methods
|
|||||||||||
|
|||||||||||
Phase multipath is one of the most crucial error sources in centimetre or millimetre level GNSS high precision positioning. Short-delay multipath is still especially difficult to detect or mitigate by the state-of-the-art hardwarebased techniques. Therefore, processing algorithm-based multipath mitigation methods are crucial for the further improvement of positioning accuracy, either integrated with other techniques or in a stand-alone mode. The effectiveness of some of these is, however, limited by the degrees of freedom in currently available solutions, i.e. insufficient satellites and signals. This problem is similar to the un-robustness and unreliability of some outlier detection techniques used in RAIM and other integrity algorithms in the current GPS system. Fig. 1 Simulated GPS three-frequency (red: L1, blue: L2, green: L5) multipath error in PRN02 in the LCPC dataset (relative permittivity =3.9) GPS modernization is being undertaken. The GPS Block IIR-M satellites are already transmitting an unencrypted civil signal (L2C) on L2 frequency and Block III satellites will transmit a new civil signal (L1C) on L1. Moreover, the signal power of L2 will be increased. This will make tracking of L2 much easier and more reliable and will increase the use of L2 in high precision kinematic applications. An additional signal, the so-called L5, will be available on GPS Block IIF satellites scheduled for launch beginning in late 2009. Both the modernised L2 and the new L5 civil signals allow coherent tracking of code and phase and so avoid the losses that occur when tracking the current P(Y) code in L2. This had led to the extensive current interest, e.g. (Hatch et al., 2000) in investigating the potential of three-frequency data for a wide range of applications. On the other hand, the European GNSS, named Galileo, is being developed to provide four carrier frequencies and its Full Operational Capability (FOC) is scheduled to be in 2013. Galileo signals are expected to be available to users in four categories: Open Service (OS), Safety-of-Life (SoL) service, Commercial Service (CS), and Public Regulated Service (PRS). This paper investigates the application of a sophisticated multiple outlier detection technique, which the author refers to as the cocktail multiple outlier detection algorithm, to multipath contaminated measurements in the three-frequency GPS and Galileo systems (only OS Galileo signals are used in this investigation, i.e., L1, E5a, and E5b) and a combined multiple-frequency GPS and Galileo system. It is tested with simulated data as real GPS L5 (only available in one Block IIR-M satellite now) and Galileo data are not yet available. The results, possibilities and weaknesses of the method are analysed. Note that the results of this investigation are not affected by the movement of the roving receiver because single-epoch data processing is used. Moreover, ambiguity resolution is not the main research interest in this paper, which assumes that ambiguities are fixed before the proposed algorithms are applied therefore no ambiguities are simulated in the carrier-phase data. Cocktail multiple outlier detection algorithmFor high precision GPS positioning, the observation equations used in parameter estimation are usually linearized as: Description of testingDescription of simulated test datasetsThe first two simulated datasets are referenced to the setup of a real multipath experiment carried out at the Laboratoire Central des Ponts et Chaussées (LCPC) near Nantes in France during May 2002 [6]. The real data was used in the validation of the multipath model used in the simulation of multipath data for the investigations of this paper, the validation results show that the model can generate very realistic phase multipath error (Lau and Cross, 2007). In the experiment, two Leica Geosystems System 530 receivers attached to lightweight AT502 antennas were used and a 5m by 2.5m steel panel was constructed and placed about 5m to one of the receiving antennas in order to create a sufficiently large multipath signal. The antenna-reflector geometry is shown in Fig. 2. The geometry of the first two simulated datasets are the same but the reflectors with the assumed relative permittivities of 10 (e.g. flint glass) and 20 are used in multipath simulation. Therefore, the test datasets are denoted as LCPC-10 and LCPC-20 according to the relative permittivity used. The baseline length is about 86 m. The information used in multipath simulation is summarised in Table 1. The sky plot of the satellites-reflector-antenna geometry is shown in Fig. 3. Moreover, details of the satellites whose signals are contaminated by multipath are shown in Table 2. Another virtual test site is at the IGS global tracking station LBCH in Long Beach, United States. In the LBCH test dataset, a concrete wall (reflector) is set 5m to the north of the antenna and a reference station 100m to its south. A relative permittivity of 7 (see Table 1) was used in the simulation of multipath in order to create strong multipath. This value is based on Stavrou and Saunders (2003), which found that the real part of the complex permittivity of concrete varied from 6.2 to 7 (for signals in the range of 1 to 95.9 GHz). Table 2 summarises the multipath simulation. Also a kinematic data set, denoted as K-HK7-300, was simulated based on a 720m trajectory along a real street in Hong Kong. The roving antenna was assumed to set about 29cm (measured to L1 phase centre) above a 1m by 1m steel carry platform that travel at 1ms-1. The buildings were assumed to be made of concrete with a relative permittivity of 7. The reference station was assumed to be about 500 m to the north of the road. Table 3 summarises the simulation – note that some satellites are blocked for some of the time. In the all of the test datasets, Leica System 530 receivers and AT502 antennas were assumed to have been used in data collection. Moreover, normal distributed random measurement noise with the standard deviation of 1 mm was generated in each phase measurement. It should be pointed out that it was necessary to assume the use of a particular receiver and antenna because factors such as receiver correlator spacing and antenna gain pattern have an impact on the simulated multipath. The author does not believe the choices the author has made affect the overall conclusions of this research. Testing methodology Five scenarios with different GNSSs or combinations of frequencies have been tested and results compared with the known positions in all cases: • Scenario 1: the current single frequency GPS data; • Scenario 2: the modernised dual-frequency GPS data, • Scenario 3: the future threefrequency GPS data, • Scenario 4: the future OS threefrequency Galileo data, and
• Scenario 5: the future OS Galileo Note that the reason for Scenario 1 is that |
Results of the cocktail multiple outlier detection (CMOD) algorithm for multipath error detectionSummaries of the positioning errors after rejection of measurements with detected ultipath by CMOD in the test scenarios 3 to 5 for LCPC-10, LCPC-20, LBCH-7 and K-HK7-300 datasets are shown in Tables 4, 6, 8, and 10 respectively. The percentage improvements (or deteriorations) when compared with single-frequency GPS and the same scenarios using MOD are also shown in the table. The results of scenarios 1 to 5 without using MOD and CMOD (i.e., standard least squares) are also shown in the table for comparison. Tables 5, 7 and 9 show the successful rates for detection of multipath errors in scenarios 3, 4 and 5 for LCPC-10, LCPC-20, and LBCH-7 datasets respectively. In the results of K-HK7-300 dataset, CMOD shows deteriorations in 3D position accuracy in scenarios 3 to 5 (see the right-most column in Table 10) when compared with the current singlefrequency GPS system or the multiplefrequency least squares only solutions. This is because the multipath errors on GPS L2 and L5 and/or Galileo E5a and E5b frequencies are highly correlated (their frequencies are very similar) when reflections occur at the very close carrying platform. The introduction of another highly correlated multipath error clearly further drags the estimated position away from the true position. Table 11 shows the coherence of phase multipath errors against the differential path delays and Table 12 shows the correlation functions between frequencies and differential path delays. The orange highlighted values in Table 12 show that the high correlations (correlation functions are greater than 0.95) of some frequency pairs occur when the differential path delays (also antenna-reflector distance) are very short (about less than or equal to one metre). The performance of CMOD is worse than that of MOD in this dataset. There are two potential reasons. Firstly, the relatively small multipath errors may not be detected by MOD but they may be detected by CMOD, which should have a positive impact. However, since there are many multipath contaminated measurements at each epoch in this dataset as shown in Table 3, many measurements are rejected, which then weakens the satellite geometry. Secondly, the highly correlated low-frequency multipath errors from the very close reflector (the carrying platform) may lead to false detection of measurements without multipath error. In addition to the problem of correlation of multipath errors, the number of multipath contaminated measurements in this dataset is much more than the number of good measurements. This in itself leads to difficulty in the detection of multipath errors. ConclusionA practical cocktail multiple outlier detection algorithm, called CMOD, is proposed to tackle the undetected outlier problem in classical multiple outlier detection method (MOD) when phase multipath errors in two or more frequencies of a satellite are in-phase. Tests with static test datasets showed that using CMOD with GPS or Galileo three-frequency data may not improve positioning accuracy when compared with MOD results. The author believes that this is due to the fact that there is insufficient redundancy to take advantage of the multiple outlier detection process used in CMOD. The result is the rejection of too many measurements/ satellites, which then weakens the satellite geometry. However, using CMOD with combined GPS and Galileo multiplefrequency data (scenario 5) shows a substantial increase in correct detection of multipath errors and significant reduction in false detection. Only one false detection occurred in only one of the static test datasets that were tested. Table 4 RMS positioning errors in millimetres and percentage improvement of using CMOD for the test scenarios in LCPC-10 dataset Table 5 Approximate percentages of correct and wrong detections of using CMOD for scenarios 3 to 5 in the LCPC-10 dataset Table 6 RMS positioning errors in millimetres and percentage improvement of using CMOD for the test scenarios in LCPC-20 dataset Table 7 Approx. percentages of correct and wrong detections of using CMOD for scenarios 3 to 5 in the LCPC-20 dataset Table 8 RMS positioning errors in millimetres and percentage improvement of using CMOD for the test scenarios in LBCH-7 Table 9 Approx. percentages of correct and wrong detections of using CMOD for scenarios 3 to 5 in the LBCH-7 dataset Table 10 RMS positioning errors in millimetres and percentage improvement of the test scenarios in K-HK7-300 dataset However CMOD shows a significant deterioration when compared with MOD and the current single-frequency GPS system when dealing with our kinematic data set. The reason has been identified as being due to a very close reflector, which led to a large number of highly correlated multipath contaminated observations (almost all signals are affected by multipath). Actually the author believe that any very close reflector (less than one metre) has the potential to lessen the advantage of using multiple-frequency GNSS data, this is because the frequencies of GPS L2 and L5, and Galileo E5a and E5b are extremely close and hence the magnitude and phase of carrier-phase multipath errors are extremely close. Note that the results of this paper are not affected by the movement (static/kinematic) of the roving receiver because single-epoch data processing is used. Table 11 Phase multipath errors of the GPS and Galileo frequencies for various differential path delay (DPD) Table 12 Correlation functions of phase multipath errors among different GPS and Galileo frequencies in Table 11 against the differential path delay (DPD) Finally the author remarks that, when combining GPS and Galileo, if the number of measurements contaminated by multipath is small compared to the total number of measurements at any epoch, and if there are no very close reflectors (or at least no reflectors that lead to small additional path lengths), then the author believes that the performance of CMOD will be better than that of MOD. Moreover in this case CMOD will be an extremely effective way to reduce the impact of multipath by 3-12% when comparing with MOD and by 65-73% when comparing with the standard least squares solution using the current reliable GPS L1 data. This in turns leads to the recommendation to avoid, as far as possible, locating GNSS antennas close to reflecting surfaces even for the future multiple-frequency GNSS. References[1] Hatch R, Jung J, Enge P, Pervan B (2000) Civilian GPS: The benefits of three frequencies. GPS Solutions, 3(4):1-9 [2] Cross PA (1994) Advanced Least Squares Applied to Position-Fixing. Working Paper No. 6, Department of Land Surveying, University of East London. ISSN 0260-9142 [3] Baarda,W (1968) A Testing Procedure for Use in Geodetic Networks. Publication on Geodesy, New Series, 2(5), Netherlands Geodetic Commission [4] Pope AJ (1976) The Statistics of Residuals and The Detection of Outliers. NOAA Technical Report NOS 65 NGS 1 [5] Lau, L. (2005) Phase Multipath Modelling and Mitigation in Multiple Frequency GPS and Galileo Positioning. PhD Dissertation, University of London. [6] Bétaille D, Maenpa J, Euler HJ, Cross PA (2003) New Approach to GPS Phase Multipath Mitigation. In: Proceedings of ION National Technical Meeting 2003 [7] Lau L, Cross P (2007) Development and Testing of a New Rigorous Ray-Tracing Approach to GNSS Carrier-Phase Multipath Modelling. J. Geodesy, 81(11):713-732. DOI 10.1007/s00190-007-0139-z [8] Stavrou S, Saunders S (2003) Review of Constitutive Parameters of Building Materials. In: Proceedings of the 12th International Conference on Antennas and Propagation, Exeter, 31st March – 3rd April 2003, pp 211-215 |
||||||||||||||||||
|
Pages: 1 2
Leave your response!