Next Article in Journal
Detection of Change Points in Pseudo-Invariant Calibration Sites Time Series Using Multi-Sensor Satellite Imagery
Next Article in Special Issue
An Improved Multipath Mitigation Method and Its Application in Real-Time Bridge Deformation Monitoring
Previous Article in Journal
Deep Learning for Detection of Visible Land Boundaries from UAV Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Algorithm for Real-Time Cycle Slip Detection and Repair for Low Elevation GPS Undifferenced Data in Different Environments

1
College of Geology Engineering and Geomatics, Chang’an University, Xi’an 710054, China
2
Shaanxi Satellite Application Technology Center of Natural Resources, Research Institute of Investigation and Monitoring, Xi’an 710119, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(11), 2078; https://doi.org/10.3390/rs13112078
Submission received: 11 May 2021 / Revised: 22 May 2021 / Accepted: 23 May 2021 / Published: 25 May 2021

Abstract

:
Real-time cycle slip detection and repair is one of the key issues in global positioning system (GPS) high precision data processing and application. In particular, when GPS stations are in special environments, such as strong ionospheric disturbance, sea, and high-voltage transmission line interference, cycle slip detection and repair in low elevation GPS observation data are more complicated than those in normal environments. For low elevation GPS undifferenced carrier phase data in different environments, a combined cycle slip detection algorithm is proposed. This method uses the first-order Gauss–Markov stochastic process to model the pseudorange multipath in the wide-lane phase minus narrow-lane pseudorange observation equation, and establishes the state equation of the wide-lane ambiguity with the pseudorange multipath as a parameter, and it uses the Kalman filter for real-time estimation and detects cycle slips based on statistical hypothesis testing with a predicted residual sequence. Meanwhile, considering there are certain correlations among low elevation, observation epoch interval, and ionospheric delay error, a second-order difference geometry-free combination cycle slip test is constructed that takes into account the elevation. By combining the two methods, real-time cycle slip detection for GPS low elevation satellite undifferenced data is achieved. A cycle slip repair method based on spatial search and objective function minimization criterion is further proposed to determine the correct solution of the cycle slips after they are detected. The whole algorithm is experimentally verified using the static and kinematic measured data of low elevation satellites under four different environments: normal condition, high-voltage transmission lines, dynamic condition in the sea, and ionospheric disturbances. The experimental results show that the algorithm can detect and repair cycle slips accurately for low elevation GPS undifferenced data, the difference between the float solution and the true value for the cycle slip does not exceed 0.5 cycle, and the differences obey the normal distribution overall. At the same time, the wide-lane ambiguity and second-order difference GF combination sequence calculated by the algorithm is smoother, which give further evidence that the algorithm for cycle slip detection and repair is feasible and effective, and has the advantage of being immune to the special observation environments.

Graphical Abstract

1. Introduction

High-quality global positioning system (GPS) observation data play a key role in obtaining high precision GPS positioning. However, when a satellite signal is blocked by an obstacle and fails to reach a receiver, or the satellite signal experiences temporary loss of lock due to external interference and harsh environment in which the receiver is located, an integer number of cycles in the phase observable jumps suddenly. This phenomenon is referred to as cycle slip [1,2,3,4]. Cycle slips will lead to destruction of continuous primary carrier phase observation values, which will seriously affect high precision gps data processing results. Therefore, accurately detecting and repairing cycle slips is an important data pre-processing task in high precision GPS positioning and applications. Basically, there are two approaches for processing cycle slips: one is to detect and repair; and the other is to use cycle slips as an unknown ambiguity parameter, which is estimated together with other unknown parameters in subsequent data processing [5,6]. Obviously, the second approach of dealing with cycle slips is not satisfactory, because it makes the data processing more complicated [5,6]. For frequent cycle slips, unpredictable situations may arise. Therefore, many cycle slip detection and repair algorithms have been developed based on the first approach [7,8]. Commonly used cycle slip detection methods include the polynomial fitting method [9,10,11], TurboEdit method [12], Bayesian theory [13], and the combination method of pseudorange and carrier phase [14,15,16]. These methods all have shortcomings. For example, the polynomial fitting method is mostly suitable for cycle slip detection of static GPS data. The cycle slip detection of Bayesian theory is effective only for high sampling GPS data without ionospheric disturbance. In the TurboEdit and combination method of pseudorange and carrier phase, small cycle slips are often missed owing to the observation noise with a large pseudorange [17,18,19]. In particular, in special observation environments like sea and high-voltage transmission lines, these methods more likely cause greater multipath errors in low elevation GPS data, and thus affect wide-lane phase minus narrow-lane pseudorange combination failures when detecting cycle slip of two cycles or less. Moreover, in a strong ionospheric disturbance environment, cycle slip detection and repair in low elevation GPS observation data are more complicated than in normal environments, and common cycle slip detection methods will also cause cycle slip detection failure or inaccuracy [20,21,22]. Nevertheless, low elevation undifferenced GPS observation data usually contain a significant amount of useful information (such as tropospheric delay and water vapor information). Therefore, how to detect and repair the cycle slip of GPS phase data in different environments is an important problem to be solved.
To tackle the above problems, a cycle slip detection combination algorithm is proposed in this study, based on the observation equation of wide-lane phase minus narrow-lane pseudorange; the correlation between the pseudorange multipath and certain time intervals is considered in the algorithm; and the first-order Gauss–Markov stochastic model is used for processing. A state equation with wide-lane ambiguity and pseudorange multipath as parameters is established and the Kalman filter is used to conduct real-time estimation. The separation of wide-lane ambiguity and pseudorange multipath is achieved and cycle slip detection is performed by constructing a statistical hypothesis testing method containing prediction residual sequences. At the same time, on the basis of the geometry-free observation equation, the correlation between low elevation, observation epoch interval, and ionospheric delay error is considered, and a second-order difference geometry-free cycle slip test quantity taking into account the elevation factor is constructed to minimize the impact of ionospheric residuals generated by ionospheric fluctuations for cycle slip detection. Real-time cycle slip detection is achieved by combining the above two methods. According to the integer number nature of cycle slips, the method of spatial search and objective function minimization criterion are used to determine the correct solution of cycle slips for repairing. The algorithm is verified using the measured static and kinematic data of each low elevation satellite under normal environment, high-voltage transmission lines, sea environment, and ionospheric disturbances in the experiment. The experimental results show that the algorithm can accurately detect and repair cycle slips in low elevation GPS satellite phase data, which is feasible in real-time processing of GPS undifferenced data cycle slips.

2. Materials and Methods

In this section, the pseudorange multipath error and cycle slip detection performance of wide-lane phase minus narrow-lane pseudorange combination are firstly analyzed. Then, the influence of ionospheric residuals on the cycle slip detection method of first-order difference geometry-free combination is discussed. On this basis, the combination algorithm of cycle slip real-time detection and repair is proposed, and the principle and steps of the whole algorithm are described in detail.

2.1. Wide-Lane Phase Minus Narrow-Lane Pseudorange (WL-NL) Combination with Pseudorange Multipath

For single epoch GPS satellites, the undifferenced observation equation of dual-frequency code pseudorange and carrier phase can be derived as follows [23,24,25]:
P 1 = ρ + c d t u d t s + I e f 1 2 + T + M P 1 + η 1
P 2 = ρ + c d t u d t s + I e f 2 2 + T + M P 2 + η 2
λ 1 φ 1 = ρ + c d t u d t s λ 1 N 1 I e f 1 2 + T + m p 1 + ε 1
λ 2 φ 2 = ρ + c d t u d t s λ 2 N 2 I e f 2 2 + T + m p 2 + ε 2
In Equations (1)–(4), ρ is the geometric distance between the satellite and receiver antenna; d t u   and   d t s are the receiver clock error and satellite clock error, respectively; I e is the ionospheric delay error; T is the tropospheric delay error; P 1   and   P 2 are the pseudorange measurements; M P 1   and   M P 2 are the pseudorange multipath errors; η 1   and   η 2 are the pseudo-range noise errors; φ 1   and   φ 2 are the carrier phase measurements; N 1   and   N 2 are the ambiguities; m p 1   and   m p 2 are the carrier phase multipath errors; ε 1   and   ε 2 are the phase noise errors; c is the speed of light; λ 1   and   λ 2 are the wavelengths; and f 1   and   f 2 are the frequencies of GPS signals. By subtracting Equation (2) from Equation (1), the pseudorange combination equation can be obtained that eliminates the clock error, tropospheric delay, and geometric distance between the satellite and the Earth:
P 1 P 2 = I e f 1 2 f 2 2 f 1 2 f 2 2 + M P 1 M P 2 + η 1 η 2
From Equation (5), the ionospheric delay error can be calculated:
I e = f 1 2 f 2 2 f 1 2 f 2 2 P 1 P 2 M P 1 M P 2 η 1 η 2
At the same time, dividing two sides of Equations (3) and (4) by λ i and then subtracting the equations yields the wide-lane phase combination equation:
φ 1 φ 2 + N 1 N 2 = f 1 f 2 c ρ + c d t u d t s + T + f 1 f 2 c f 1 f 2 I e + f 1 m p 1 f 2 m p 2 c + f 1 ε 1 f 2 ε 2 c
Substituting Equation (6) into Equation (1) to solve the term ( ρ + c d t u d t s + T yields
ρ + c d t u d t s + T = f 1 2 P 1 f 2 2 P 2 f 1 2 f 2 2 f 1 2 M P 1 f 2 2 M P 2 f 1 2 f 2 2 f 1 2 η 1 f 2 2 η 2 f 1 2 f 2 2
Substituting Equations (8) and (6) into the terms ρ + c d t u d t s + T and I e of Equation (7) yields [26]
λ w φ 1 φ 2 + λ w N w = f 1 P 1 + f 2 P 2 f 1 + f 2 f 1 M P 1 + f 2 M P 2 f 1 + f 2   f 1 η 1 + f 2 η 2 f 1 + f 2 + f 1 m p 1 f 2 m p 2     f 1 f 2 + f 1 ε 1 f 2 ε 2     f 1 f 2
In Equation (9), N w = N 1 N 2 is the wide-lane ambiguity and λ w = c f 1 f 2 is the wide-lane carrier phase wavelength. Because the observation noise of the carrier phase is 2 mm, the maximum multipath error does not exceed 5 cm, so the item f 1 m p 1 f 2 m p 2     f 1 f 2 + f 1 ε 1 f 2 ε 2 f 1 f 2 can be ignored. Meanwhile, letting M P = f 1 M P 1 + f 2 M P 2 f 1 + f 2 and η = f 1 η 1 + f 2 η 2 f 1 + f 2 , the relationship between the combined pseudorange multipath error M P and wide-lane ambiguity N w can be obtained as
λ w φ 1 φ 2 + λ w N w = f 1 P 1 + f 2 P 2 f 1 + f 2 M P η
It can be seen from Equation (10) that, when using it for cycle slip detection, the accuracy of cycle slip detection mainly comes from the pseudorange multipath error M P and measurement noise η ; in an ideal situation, the mean value of the multipath error and measurement noise of pseudorange measurement should be close to 0, and the pseudorange measurement error M e r r o r introduced by M P and η can be approximately calculated by
M L W = N M W + M P + η = f 1 P 1 f 2 P 2 f 1 + f 2 λ w φ 1 φ 2
M e r r o r = M L W M e a n M L W = M P + η
In Equations (11) and (12), M L W is the measurement error of the pseudorange combination including the wide-lane ambiguity, and N M W = λ w N w , N M W should be a constant when there is no cycle slip. M e a n M L W is the mean of M L W within the observation arc.
To intuitively analyze the magnitude of the pseudorange measurement error ( M e r r o r ), the measured GPS data from the BJFS reference station are used as an example. The data were observed on 1 January 2014 and the sampling rate was 30 s. The code observation types are C/A and P2. The dual frequency observations corresponding to satellites G26 and G08 are selected and calculated according to Equation (12). The relationship between the pseudorange measurement error sequence of G26 and G08 and their elevation is shown in Figure 1.
It can be seen from Figure 1 that, even when the satellite is at a high elevation (above 30°), its pseudorange measurement error can reach ±0.4 m; when the satellite is at a low elevation (5–20°), the effect of its pseudorange measurement error can reach ±1.0 m; therefore, when the wide-lane phase minus the narrow-lane pseudorange combination is used to calculate the ambiguity N M W , the pseudorange measurement error must be dealt with, otherwise the accuracy of ambiguity will be reduced.

2.2. Inaccurate Detection of Cycle Slip by Conventional MW Combination

The analysis in Section 2.1 indicates that the pseudorange multipath in the wide-lane phase minus narrow-lane pseudorange combination (i.e., MW combination) is the main factor affecting the accuracy of the wide-lane ambiguity calculation; when the conventional MW combination is used to detect cycle slips, the mean recursive algorithm is used to obtain the mean of the wide-lane ambiguity N ¯ M W i and the root mean square σ i for each epoch, and the cycle slip testing is performed by the following equation [12,18,26]:
N M W i N ¯ M W i 1 4 σ i
Because the accuracy of the P code pseudorange is approximately 0.3 m, and the MW combination is also affected by the pseudorange measurement error, this method is not sensitive to small cycle slips (two cycles or below), and thus inaccurate cycle slip detection may occur. The above mentioned G08 satellite is used as an example again; at the 50th and 250th epochs of the L1 carrier phase observation, cycle slips of injecting one and two cycles are simulated. The conventional MW combination cycle slip detection results are shown in Figure 2. It can be seen from the results in Figure 2 that, for low elevation satellites, the cycle slips identified by conventional MW combination for two cycles or less are invalid.

2.3. Cycle Slip Detection with Fused Kalman Filter and Hypothesis Testing

As GPS receivers cannot usually obtain P 1 code observation, in general, C/A code is used instead of P 1 code to perform the MW combination in Equation (10). The observation accuracy of C/A code is lower than that of P 1 code, which is greatly affected by pseudorange multipath, so the accuracy of the calculated wide-lane ambiguity value is reduced. Therefore, the pseudorange multipath error M P in Equation (11) must be separated from the wide-lane ambiguity N M W . As the pseudorange multipath error is not constant and correlated with a certain time interval, the first-order Gauss–Markov process can be used to estimate the multipath error [27,28]:
M P ^ t = 1 τ M P t
where M P ^ t is the change rate of multipath error and τ is the correlation time. Meanwhile, when the carrier phase has no cycle slip, the ambiguity is equal, and a state equation with the wide-lane ambiguity and pseudorange multipath error as the state parameters can be established:
M P ^ t 0 = 1 / τ 0 0 0 M P t N M W
Using Equation (11) as the observation equation of the discrete system, combined with Equation (14), the following Kalman filter model can be established [29,30,31]:
X k = ϕ k , k 1 X k 1 + ω k ω k 0 , Q ω k L k = A k X k + v k v k 0 , R v k
where X k   and   X k 1 are the state parameter vectors, ϕ k , k 1 is the state transition matrix and ϕ k , k 1 = e x p Δ t τ 0 0 1 , L k is the observation vector and L k = f 1 P 1 f 2 P 2 f 1 + f 2 λ w φ 1 φ 2 , and A k is the coefficient matrix and A k = 1 1 . ω k   and   Q ω k are the state noise and its variance matrix, respectively, Q k = q τ 2 1 e x p 2 Δ t τ 0 0 10 15 , and v k   and   R v k are the observation noise and its variance matrix, respectively [32,33]. The Kalman filter recursive equation is used to achieve the separation of wide-lane ambiguity and pseudorange multipath error by epoch, and the best estimation of wide lane ambiguity is obtained.
Meanwhile, assuming that there is no cycle slip in the carrier phase observation of the first k 1 epochs and cycle slip occurs at the k th epoch, the Kalman filter model in Equation (17) becomes
X k = ϕ k , k 1 X k 1 + ω k L k = A k X k + C k + v k
The prediction residual of the k th epoch is then
V k = L k A k X ¯ k , k 1 = A k X k + C k + v k A k ϕ k , k 1 X k 1 , k 1 = A k ϕ k , k 1 X k 1 + ω k + C k + v k A k ϕ k , k 1 X k 1 , k 1 = A k ϕ k , k 1 X k 1 X k 1 , k 1 + A k ω k + C k + v k
At this point, E V k = C k ,   Q V k = A k Q ¯ k , k 1 A k T + R v k . It can be seen from Equation (18) that, when a cycle slip occurs, the prediction residual sequence V k will shift and no longer obey the zero-mean normal distribution [34,35]. Thus, cycle slip detection is performed by constructing a statistical hypothesis testing containing the predicted residual sequence, with the null hypothesis being that the cycle slip does not occur and the alternative hypothesis being that the cycle slip occurs, namely,
H 0 V :   cycle   slip   does   not   occur ,         V k ~ N 0 , Q V k H 1 V :   cycle   slip   occurs ,         V k ~ N C k , Q V k
The testing is performed based on the following criteria:
If V k < n Q V k , then H 0 V is accepted.
If V k n Q V k , then H 0 V is rejected and H 1 V is accepted.
Here, n is a weighting coefficient, which can be three (99.7% confidence level) or four (99.9% confidence level). Cycle slip detection is performed up to the last epoch using this criterion.

2.4. Detecting Cycle Slips with Second-Order Difference Geometry-Free Combination Considering Elevation

Based on the MW combination and using the fused method of the Kalman filter and hypothesis testing, cycle slip detection can be performed effectively. However, when the equal cycle slip occurs at frequencies L1 and L2, detection blind spots may appear. Therefore, a difference geometry-free (GF) combination taking into account the elevation factor is proposed here to correct the blind spots of the cycle slip detection method described in Section 2.3. By subtracting Equation (4) from Equation (3), the GF combination observation equation can be obtained to eliminate the clock error and tropospheric delay error [36,37]:
L G F = λ 1 φ 1 λ 2 φ 2 = λ 2 N 2 λ 1 N 1 + 1 f 2 2 1 f 1 2 I e + m p 1 m p 1 ) + ( ε 1 ε 2
Here, the multipath m p 1 m p 1 and observation noise ε 1 ε 2 on frequencies L1 and L2 are approximately equal, so the m p 1 m p 1 + ε 1 ε 2 term on the right-hand side of the equation can be ignored, and the corresponding GF combination only contains the effects of the ionosphere and the ambiguity. When a cycle slip occurs in the carrier phase observation of the i th epoch, differencing the above equation between adjacent epochs yields the first-order difference GF combination [18,25] containing the ionospheric residual I e r e s = I e i I e i 1 and cycle slip information λ 2 Δ N 2 λ 1 Δ N 1 :
Δ L G F = L G F i L G F i 1 = λ 2 Δ N 2 λ 1 Δ N 1 + 1 f 2 2 1 f 1 2 I e r e s
Under normal circumstances, the Δ L G F sequence appears as a rapidly changing curve due to cycle slip. Meanwhile, it can be seen from Equation (21) that the cycle slip information is affected by the ionospheric residual I e r e s , and especially when GPS satellites are in an active period of the ionosphere, the ionospheric residual term will inevitably cause disturbances in the Δ L G F sequence of Equation (21), resulting in inaccurate cycle slip detection [20,21]. Therefore, to minimize the effect of the ionospheric residual term, considering the correlation between Δ L G F and the elevation (especially the low elevation), and the epoch time interval, the second-order difference GF combination cycle slip detection quantity is constructed that takes into account the elevation, as expressed by the following equation:
Δ L G F E l e = λ 2 Δ N 2 λ 1 Δ N 1 s i n e l e = L G F i 2 L G F i 1 + L G F i 2 1 f 2 2 1 f 1 2 I e i 2 I e i 1 + I e i 2 s i n e l e
The root mean square of the first i epochs is recursively calculated by Eqution (23):
σ Δ L G F E l e 2 i = i 2 i 1 σ Δ L G F E l e 2 i 1 + 1 i Δ L G F E l e i 2
Then, the cycle slip detection threshold condition is expressed by the following equation:
Δ L G F E l e i 5 σ Δ L G F E l e i , cycle   slip   occurs Δ L G F E l e i < 5 σ Δ L G F E l e i , cycle   slip   does   not   occurs
The cycle slip detection of low elevation GPS carrier phases in different environments is performed based on the combination of the above two methods.

2.5. Cycle Slip Repair Based on Spatial Search and Objective Function Minimization Criterion

The methods proposed in Section 2.3 and Section 2.4 can achieve real-time detection of cycle slips. By analyzing Equations (17) and (18), it can be seen that the value of the cycle slips can be expressed by residuals. Equation (22) shows that the cycle slip value can be expressed by the difference between the epochs of L G F . However, when the two frequencies have equal cycle slips, or the cycle slips of the two frequencies are proportional to the wavelength, a detection blind zone that is not sensitive to cycle slips manifests when the two methods are used alone. Therefore, in actual application, the two methods need to be combined effectively. Together, they make up for each others’ shortcomings, and the occurring frequency and value of cycle slips can be determined by Equation (25).
c f 1 f 2 c f 1 f 2 λ 1 λ 2 Δ N 1 Δ N 2 = C k   Δ L G F
Cycle slips are expressed in integer numbers, while the solution given by Equation (22) is the float solution of cycle slips [38]; to avoid the effect of rounding errors caused by direct rounding, a method using spatial search and objective function minimization criteria is proposed. To obtain the integer solution of cycle slip, the specific steps are as follows:
  • Round the float value calculated by Equation (25) as the initial value of the cycle slip.
  • The effects of low elevation and ionospheric disturbance period are fully considered, the initial value of cycle slip is taken as the center, ±5 cycles is the search range, and one cycle is the search step, in order to form a cycle slip candidate combination.
  • Use the cycle slip candidate combination constructed in Step 2 to repair the carrier phase observation in real time, then re-check it through the cycle slip detection methods proposed in this study, and store the cycle slip combination that the two cycle slip detection methods meet the condition of no cycle slip simultaneously.
  • Use the cycle slip combination that meets the detection condition to repair the carrier phase observation, and substitute the repaired carrier phase into the objective function, so that the set of cycle slip combinations that makes the objective function in Equation (26) meet the minimization criterion is the correct cycle slip solution.
    L k L k 1 2 + L G F k L G F k 1 2 = m i n

3. Results

In this section, the GPS data in four different environments are firstly introduced, then the pseudorange multipath error of low elevation satellite is analyzed. Based on the experimental schemes of small cycle slips, large cycle slips and special cycle slips were added into the simulation at different epochs of GPS data, and the cycle slip detection and repair method proposed in this study is verified and analyzed.

3.1. Measured GPS Data in Four Different Environments

To analyze and verify the feasibility of the cycle slip detection and repair methods proposed in this paper, the following four different environments and different sampling rates of GPS static and kinematic measured data were used for testing.

3.1.1. Normal Environment

The GPS data of the P090 reference station measured on 3 January 2010 were selected for testing. The sampling rate of the data was 15 s. The satellite G09 is a low elevation satellite. The elevation and pseudorange measurement error sequence are shown in Figure 3, from which it can be seen that, even if the IGS reference station is in a good location environment, the influence of low elevation on the pseudorange measurement error is still large.

3.1.2. High-Voltage Transmission Lines Environment

The measured data of the GPS receiver at station TB06, which is located directly under high-voltage transmission lines and is closest to the transmission tower (see Figure 4), were selected for testing. The data were taken on 2 June 2016 with a sampling rate of 1 s. G03 is a low elevation satellite; its elevation and pseudorange measurement error sequence are shown in Figure 5.

3.1.3. Dynamic Environment in the Sea Area

The measured kinematic data of GPS receiver installed on a buoy in a sea region of Indonesia were selected for testing. The GPS trajectory on the buoy is shown in Figure 6. The data were taken on 3 March 2010 with a sampling rate of 5 s. G10 is a low elevation satellite; the corresponding elevation and pseudorange measurement error sequence are shown in Figure 7.

3.1.4. Ionospheric Disturbance Environment

Cai et al. [18] show that the solar flare intensity on 15 April 2001 was 14.4, which was the highest flare intensity in the last maximum solar activity cycle from the year 2000 to 2002. The ionosphere was thus in a period of strong disturbance. The actual measurement data of IGS reference station AMC2 in this environment were used for testing, and its sampling rate was 30 s. G26 is a low elevation satellite; the corresponding elevation and pseudorange measurement error sequence are shown in Figure 8.
It can be seen from Figure 3, Figure 5, Figure 7 and Figure 8 that the pseudorange measurement error of each low elevation satellite varies greatly in the different special environmental conditions, which will inevitably affect the accuracy of conventional MW and GF combination cycle slip detection.

3.2. Validation of Cycle Slip Detection and Repair Results

To test the feasibility of the cycle slip detection and repair algorithm developed in this study, large cycle slips, small cycle slips, and special cycle slip pairs were added into the simulation at different epochs of the low elevation measured data sequence under the above four different environmental conditions. The cycle slip pairs are shown in Table 1.
The corresponding data of each low elevation in the above four different environmental conditions are used to test the cycle slip detection and repair method proposed in this paper. The cycle slip detection results are shown in Figure 9, Figure 10, Figure 11 and Figure 12.
It can be seen from Figure 9, Figure 10, Figure 11 and Figure 12 that the cycle slip combination detection algorithm proposed in this paper accurately detects the epoch position in which the cycle slips occur, and is not affected by the sampling rate of the observation data. Thus, this algorithm is feasible for cycle slip detection of low elevation dual-frequency undifferenced GPS phase data in different special environments, and it is easy to implement.
At the same time, to evaluate the cycle slip repair method based on the spatial search and objective function minimization criterion proposed in this paper, Table 2 shows the cycle slip repair results of each low elevation under the aforementioned four different environments. The statistics of the difference between the float solution and its respective true value are shown in Figure 13.
It can be seen from Table 2 and Figure 13 that, when the cycle slip repair method proposed in this paper is used to obtain the initial float solution, the difference between the float solution and the true value for the cycle slip does not exceed 0.5 cycles, and the differences obey the normal distribution overall. This indicates that the float solution has high accuracy, and the cycle slip results obtained through spatial search and objective function minimization criterion are accurate and reliable, which further demonstrates the effectiveness of the cycle slip repair method.

4. Discussion

According to the experimental results in Section 3.2, for different sampling rates of low elevation GPS phase data, it is necessary to effectively combine the two cycle slip detection methods proposed in this paper to achieve accurate cycle slip detection and repair in four different environments. Only at the 500th epoch in Figure 11, the cycle slip detection method with the fused Kalman filter and hypothesis testing based on MW combination cannot detect the small cycle slip of one cycle. The reason is that the acquisition environment of the dynamic data is located in the sea, which is one of the main error sources of multipath error, and as there is no P1 code observation in the data, the observation noise is large when the C/A code is used to replace the P code for MW combination. However, the cycle slip is accurately detected by the second-order differential GF method considering the elevation, and the correct cycle slip value is finally found. Therefore, for small cycle slips such as (1,0), if the observation data contains P code, the performance of the combined cycle slip detection algorithm proposed in this paper will be perfect.
To further evaluate and analyze the performance of the cycle slip detection algorithm, satellite G10 was used under the dynamic environment in the sea area, and satellite G26 in the ionospheric disturbance environment as examples to give the ambiguity sequence after separating the pseudorange multipath based on Kalman filter, and the second-order difference GF combination sequence taking into account the satellite elevation; the corresponding results are shown in Figure 14 and Figure 15, respectively.
It can be seen from Figure 14 and Figure 15 that the wide-lane ambiguity in the MW combination estimated by the Kalman filter is smoother than the wide-lane ambiguity sequence determined by the conventional MW combination, which indicated that the Kalman filter method used in this paper can effectively separate the pseudorange multipath error, improve the reliability of the wide-lane ambiguity estimation, and improve the cycle slip detection capability of the MW combination. At the same time, the second-order difference GF combination method that considers the satellite elevation is more stable than the conventional first-order difference GF combination. In particular, in the case of ionospheric disturbance, the second-order difference GF combination that takes into account the low satellite elevation is better in reducing the effect of ionospheric error on the cycle slip testing quantities, and further improves the GF combination cycle slip detection capability under special environments.
Analysis of experimental results further shows that, on the basis of MW combination, the wide-lane ambiguity can be separated from the pseudorange multipath through the Kalman filter, thereby improving the accuracy of float solutions of wide-lane ambiguity and providing a reference data processing method for subsequent ambiguity fixes. At the same time, the second-order difference GF combination taking into account the low elevation can effectively reduce the ionospheric residual error, and has the advantage of being immune to the strong ionospheric disturbance observation environment.

5. Conclusions

In this study, for the cycle slip detection problem in GPS undifferenced carrier phase data at low elevation in different environments, on the basis of MW combination, cycle slip detection is carried out by constructing cycle slip test quantities, fused Kalman filter, hypothesis testing, and a second-order difference GF combination is proposed that takes into account the elevation factor. The experimental test results of data with different sampling rates at various low elevation under four special environments show that the cycle slip combination detection algorithm proposed in this paper can detect the epoch positions of small and large cycle slips accurately in real time, and it is feasible for cycle slip detection and easy to implement.
After the cycle slips are accurately detected, the experimental results of cycle slip repair show that the proposed method based on spatial search and objective function minimization criterion can achieve accurate results for the cycle slip integer solution. This method is effective for repairing cycle slips of GPS low elevation satellite data.

Author Contributions

N.L. designed the experiments and wrote the main manuscript. Q.Z. and S.Z. reviewed the paper. X.W. edited the paper. All components of this research were carried out under N.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by China Postdoctoral Science Foundation Funded Project (No.: 2016M592733) and the National Key Research and Development Program of China (No.: 2020YFC1512000).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets analyzed in this study are managed by the College of Geology Engineering and Geomatics, Chang’an University and can be available on request from the corresponding author.

Acknowledgments

The authors would like to thank IGS for providing the relevant GPS observations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bisnath, S.; Kim, D.; Langley, R. A new approach to an old problem: Carrier phase cycle slips. GPS World 2001, 13, 42–49. [Google Scholar]
  2. Lee, H.K.; Wang, J.L.; Rizos, C. Effective cycle slip detection and identification for high precision GPS/INS integrated systems. J. Navig. 2003, 56, 475–486. [Google Scholar] [CrossRef] [Green Version]
  3. Karaim, M.; Karamat, T.; Noureldin, A.; Ei-Shafie, A. GPS cycle slip detection and correction at measurement level. Br. J. Appl. Sci. Technol. 2014, 4, 4239–4251. [Google Scholar]
  4. Wang, Z.; Wang, X.; Ji, S.; Chen, W. Optimal cycle slip detection and correction with reliability. Surv. Rev. 2016, 48, 233–239. [Google Scholar] [CrossRef]
  5. Liu, Z.Z. A new automated cycle slip detection and repair method for a single dual-frequency GPS receiver. J. Geod. 2011, 85, 171–183. [Google Scholar]
  6. Wei, H.; Li, J.C.; Zhang, S.J.; Xu, X.Y. Cycle slip detection and repair for dual-frequency LEO satellite GPS carrier phase observations with orbit dynamic model information. Remote Sens. 2019, 11, 1273. [Google Scholar] [CrossRef] [Green Version]
  7. Zhen, D. MATLAB software for GPS cycle-slip processing. GPS Solut. 2012, 16, 267–272. [Google Scholar]
  8. Liu, Z.Z. A new approach for cycle slip detection and fix using single GPS receiver’s single satellite dual frequency data containing arbitrarily large pseudorange errors. J. Glob. Pos. Syst. 2018, 16, 1–21. [Google Scholar] [CrossRef] [Green Version]
  9. Lichtenegger, H.; Hofmann-Wellenhof, B. GPS-data preprocessing for cycle-slip detection. In Global Positioning System: An Overview; Springer: Berlin/Heidelberg, Germany, 1990; pp. 57–68. [Google Scholar]
  10. Ju, B.; Gu, D.F.; Chang, X.; Herring, T.; Duan, X.J.; Wang, Z.M. Enhanced cycle slip detection method for dual-frequency BeiDou GEO carrier phase observations. GPS Solut. 2017, 21, 1227–1238. [Google Scholar] [CrossRef]
  11. Li, B.F.; Liu, T.X.; Nie, L.W.; Qin, Y.N. Single-frequency GNSS cycle slip estimation with positional polynomial constraint. J. Geod. 2019, 93, 1781–1803. [Google Scholar] [CrossRef]
  12. Blewitt, G. An Automatic Editing Algorithm for GPS data. Geophys. Res. Lett. 1990, 17, 199–202. [Google Scholar] [CrossRef] [Green Version]
  13. Lacy, M.; Reguzzoni, M.; Sanso, F.; Venuti, G. The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem. J. Geod. 2008, 82, 527–542. [Google Scholar] [CrossRef]
  14. Melbourne, W.G. The case for ranging in GPS based geodetic systems. In Proceedings of the 1st International Symposium on Precise Positioning with the Global Positioning System, Rockville, ML, USA, 15–19 April 1985; pp. 373–386. [Google Scholar]
  15. Hatch, R. The synergism of GPS code and carrier measurements. In Proceedings of the Third International Symposium on Satellite Doppler Positioning at Physical Sciences Laboratory of New Mexico State University, Las Cruces, NM, USA, 8–12 February 1982; pp. 1213–1231. [Google Scholar]
  16. Gao, Y.; Li, Z.F. Cycle slip detection and ambiguity resolution algorithms for dual-frequency GPS data processing. Mar. Geod. 1999, 22, 169–181. [Google Scholar]
  17. Ji, S.Y.; Chen, W.; Weng, D.J.; Wang, Z.J.; Ding, X.L. A study on cycle slip detection and correction in case of ionospheric scintillation. Adv. Space Res. 2013, 51, 742–753. [Google Scholar] [CrossRef]
  18. Cai, C.S.; Liu, Z.Z.; Xia, P.F.; Dai, W.J. Cycle slip detection and repair for undifferenced GPS observations under high ionospheric activity. GPS Solut. 2013, 17, 247–260. [Google Scholar] [CrossRef]
  19. Chen, D.Z.; Ye, S.R.; Zhou, W.; Liu, Y.Y.; Jiang, P.; Tang, W.M.; Yuan, B.; Zhao, L.W. A double-differenced cycle slip detection and repair method for GNSS CORS network. GPS Solut. 2016, 20, 439–450. [Google Scholar] [CrossRef]
  20. Ji, S.Y.; Wang, Z.J.; Chen, W.; Weng, D.j.; Xu, Y.; Fan, S.J.; Huang, B.H.; Sun, G.Y.; Wang, H.Q.; He, W.Y. Geometry-free and non-geometry-free testing quantities for cycle slip detection and correction in case of strong atmospheric variations with static observations. Surv. Rev. 2014, 335, 104–111. [Google Scholar] [CrossRef]
  21. Tang, L.; Zheng, K.; Li, X.X. Analysis of geometry-free residuals in case of traveling ionosphere disturbances and their impact cycle slip detection. GPS Solut. 2017, 21, 1221–1226. [Google Scholar] [CrossRef]
  22. Luo, X.M.; Liu, Z.Z.; Lou, Y.D.; Gu, S.F.; Chen, B.Y. A study of multi-GNSS ionospheric scintillation and cycle-slip over Hong Kong region for moderate solar flux conditions. Adv. Space Res. 2017, 60, 1039–1053. [Google Scholar] [CrossRef]
  23. Miao, Y.; Sun, Z.W.; Wu, S.N. Error analysis and cycle-slip detection research on satellite-borne GPS observation. J. Aerosp. Eng. 2011, 24, 95–101. [Google Scholar] [CrossRef]
  24. Li, P.; Jiang, X.Y.; Zhang, X.H.; Ge, M.R.; Schuh, H. Kalman-filter-based undifferenced cycle slip estimation in real-time precise point positioning. GPS Solut. 2019, 23, 99. [Google Scholar] [CrossRef]
  25. Ye, S.R.; Liu, Y.Y.; Song, W.W.; Lou, Y.D.; Yi, W.T.; Zhang, R.; Jiang, P.; Xiang, Y. A cycle slip fixing method with GPS + GLONASS observations in real-time kinematic PPP. GPS Solut. 2016, 20, 101–110. [Google Scholar] [CrossRef]
  26. Bisnath, S.B. Efficient automated cycle-slip correction of dual-frequency kinematic GPS data. In Proceedings of the Institute of Navigation GPS 2000, Salt Lake City, UT, USA, 19–22 September 2000; pp. 145–154. [Google Scholar]
  27. Yang, Y.; Hatch, R.R.; Sharpe, R.T. GPS multipath mitigation in measurement domain and its applications for high accuracy navigation. In Proceedings of the International Technical Meeting of the Satellite Division of the Institute of Navigation, Long Beach, CA, USA, 21–24 September 2004. [Google Scholar]
  28. Farrell, J.; Givargis, T. Differential GPS reference station algorithm—Design and analysis. IEEE Trans. Contr. Syst. Technol. 2000, 8, 519–531. [Google Scholar] [CrossRef]
  29. Bastos, L.; Landau, H. Fixing cycle slips in dual-frequency kinematic GPS-applications using Kalman filtering. Manuscr. Geod. 1988, 13, 249–256. [Google Scholar]
  30. Lin, S.G.; Yu, F.C. Cycle slips detection algorithm for low cost single frequency GPS RTK positioning. Surv. Rev. 2013, 330, 206–214. [Google Scholar] [CrossRef]
  31. Mertikas, S.P.; Rizos, C. On-line detection of abrupt changes in the carrier-phase measurements of GPS. J. Geod. 1997, 71, 469–482. [Google Scholar] [CrossRef]
  32. Yu, X.; Gao, J. Kinematic Precise Point Positioning Using Multi-Constellation Global Navigation Satellite System (GNSS) Observations. ISPRS Int. J. GeoInf. 2017, 6, 1–15. [Google Scholar] [CrossRef] [Green Version]
  33. Ansari, K. Real-Time Positioning Based on Kalman Filter and Implication of Singular Spectrum Analysis. IEEE Geosci. Remote Sens. Lett. 2021, 18, 58–61. [Google Scholar] [CrossRef]
  34. Fujita, S.; Saito, S.; Yoshihara, T. Cycle slip detection and correction methods with time-differenced model for single frequency GNSS applications. Trans. Inst. Syst. Control. Inf. Eng. 2013, 26, 8–15. [Google Scholar] [CrossRef] [Green Version]
  35. Teunissen, P.J.G. Quality control in integrated navigation systems. IEEE Aerosp. Electron. Syst. Mag. 2002, 5, 35–41. [Google Scholar] [CrossRef]
  36. Feng, W.; Zhao, Y.H.; Zhou, L.T.; Huang, D.F.; Hassan, A. Fast cycle slip determination for high-rate multi-GNSS RTK using modified geometry-free phase combination. GPS Solut. 2020, 24, 42. [Google Scholar] [CrossRef]
  37. Qian, C.; Liu, H.; Zhang, M.; Shu, B.; Xu, L.W.; Zhang, R.F. A geometry-based cycle slip detection and repair method with time-differenced carrier phase (TDCP) for a single frequency Global Position System (GPS) + BeiDou navigation satellite system (BDS) receiver. Sensors 2016, 16, 2064. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Zhang, W.H.; Wang, J.L. A real-time cycle slip repair method using the multi-epoch geometry-based model. GPS Solut. 2021, 25, 60–71. [Google Scholar] [CrossRef]
Figure 1. (a) Elevation and pseudorange measurement error sequence of satellite G26; (b) elevation and pseudorange measurement error sequence of satellite G08.
Figure 1. (a) Elevation and pseudorange measurement error sequence of satellite G26; (b) elevation and pseudorange measurement error sequence of satellite G08.
Remotesensing 13 02078 g001
Figure 2. Results of conventional MW combination detection of small cycle slips.
Figure 2. Results of conventional MW combination detection of small cycle slips.
Remotesensing 13 02078 g002
Figure 3. Elevation and pseudorange measurement error sequence of satellite G09.
Figure 3. Elevation and pseudorange measurement error sequence of satellite G09.
Remotesensing 13 02078 g003
Figure 4. GPS station TB06 under high-voltage transmission line environment (the latitude and longitude of station TB06 are 34°08′01″ and 107°54′38″, respectively).
Figure 4. GPS station TB06 under high-voltage transmission line environment (the latitude and longitude of station TB06 are 34°08′01″ and 107°54′38″, respectively).
Remotesensing 13 02078 g004
Figure 5. Elevation and pseudorange measurement error sequence of satellite G03.
Figure 5. Elevation and pseudorange measurement error sequence of satellite G03.
Remotesensing 13 02078 g005
Figure 6. GPS trajectory on buoy.
Figure 6. GPS trajectory on buoy.
Remotesensing 13 02078 g006
Figure 7. Elevation and pseudorange measurement error sequence of satellite G10.
Figure 7. Elevation and pseudorange measurement error sequence of satellite G10.
Remotesensing 13 02078 g007
Figure 8. Elevation and pseudorange measurement error sequence of satellite G26.
Figure 8. Elevation and pseudorange measurement error sequence of satellite G26.
Remotesensing 13 02078 g008
Figure 9. Cycle slip detection results of satellite G09 at different epochs under normal environment.
Figure 9. Cycle slip detection results of satellite G09 at different epochs under normal environment.
Remotesensing 13 02078 g009aRemotesensing 13 02078 g009bRemotesensing 13 02078 g009c
Figure 10. Cycle slip detection results of satellite G03 at different epochs under high-voltage transmission line environment.
Figure 10. Cycle slip detection results of satellite G03 at different epochs under high-voltage transmission line environment.
Remotesensing 13 02078 g010aRemotesensing 13 02078 g010bRemotesensing 13 02078 g010c
Figure 11. Cycle slip detection results of satellite G10 at different epochs under dynamic sea environment.
Figure 11. Cycle slip detection results of satellite G10 at different epochs under dynamic sea environment.
Remotesensing 13 02078 g011aRemotesensing 13 02078 g011bRemotesensing 13 02078 g011c
Figure 12. Cycle slip detection results of satellite G26 at different epochs in ionospheric disturbance environment.
Figure 12. Cycle slip detection results of satellite G26 at different epochs in ionospheric disturbance environment.
Remotesensing 13 02078 g012aRemotesensing 13 02078 g012b
Figure 13. Statistics of difference between float solution and true value for cycle slip.
Figure 13. Statistics of difference between float solution and true value for cycle slip.
Remotesensing 13 02078 g013
Figure 14. Wide-lane ambiguity and second-order difference GF combination sequence of satellite G10 under dynamic sea environment.
Figure 14. Wide-lane ambiguity and second-order difference GF combination sequence of satellite G10 under dynamic sea environment.
Remotesensing 13 02078 g014
Figure 15. Wide-lane ambiguity and second-order difference GF combination sequence of satellite G26 in ionospheric disturbance environment.
Figure 15. Wide-lane ambiguity and second-order difference GF combination sequence of satellite G26 in ionospheric disturbance environment.
Remotesensing 13 02078 g015
Table 1. Cycle slip pairs in carrier phase data of each low elevation satellite.
Table 1. Cycle slip pairs in carrier phase data of each low elevation satellite.
Environment and Satellite NumberSampling RateEpochCycle Slip Pairs
Normal, G0915 s50th(1, 1)
100th(0, 2)
150th(0, 1)
200th(9, 7)
250th(−10, 10)
300th(50, −50)
350th(77, 60)
High-voltage power lines, G031 s100th(1, 1)
150th(0, 2)
200th(9, 7)
250th(−5, 5)
300th(1, 0)
350th(77, 60)
400th(−5, −4)
450th(10, −10)
Dynamic sea environment, G105 s100th(1, 1)
200th(0, 2)
300th(9, 7)
400th(−5, 5)
500th(1, 0)
600th(77, 60)
700th(−4, −5)
Ionospheric disturbance, G2630 s50th(1, 1)
100th(0, 2)
150th(5, 4)
200th(9, 7)
250th(−77, −60)
Table 2. Cycle slip repair results.
Table 2. Cycle slip repair results.
Environment and Satellite NumberEpochFloat SolutionObjective Function ValueCycle Slip Search Results
Normal, G0950th(1.038, 1.038)3.719 × 10−3(1, 1)
100th(−0.021, 1.979)6.897 × 10−2(0, 2)
150th(0.035, 1.035)4.656 × 10−2(0, 1)
200th(9.022, 7.022)5.780 × 10−3(9, 7)
250th(−9.743, 10.257)6.176 × 10−2(−10, 10)
300th(49.439, −50.561)9.803 × 10−2(50, −50)
350th(77.334, 60.334)9.281 × 10−2(77, 60)
High-voltage power lines, G03100th(1.024, 1.024)1.022 × 10−4(1, 1)
150th(−0.020, 1.979)6.975 × 10−2(0, 2)
200th(9.066, 7.066)6.292 × 10−2(9, 7)
250th(−5.042, 4.958)4.904 × 10−1(−5, 5)
300th(0.978, −0.022)3.560 × 10−2(1, 0)
350th(77.020, 60.020)3.144 × 10−2(77, 60)
400th(−4.996, −3.996)1.784 × 10−2(−5, −4)
450th(10.013, −9.987)2.059 × 10−2(10, −10)
Dynamic sea environment, G10100th(1.046, 1.046)9.862 × 10−4(1, 1)
200th(−0.004, 1.996)5.427 × 10−3(0, 2)
300th(8.992, 6.992)3.087 × 10−2(9, 7)
400th(−4.965, 5.035)7.906 × 10−3(−5, 5)
500th(1.063, 0.063)3.188 × 10−3(1, 0)
600th(77.084, 60.084)2.363 × 10−2(77, 60)
700th(−4.051, −5.051)6.672 × 10−2(−4, −5)
Ionospheric disturbance, G2650th(1.017, 1.017)1.358 × 10−4(1, 1)
100th(−0.114, 1.886)6.133 × 10−2(0, 2)
150th(4.995, 3.995)2.241 × 10−1(5, 4)
200th(9.080, 7.080)1.196 × 10−1(9, 7)
250th(−77.006, −60.006)2.114 × 10−1(−77, −60)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, N.; Zhang, Q.; Zhang, S.; Wu, X. Algorithm for Real-Time Cycle Slip Detection and Repair for Low Elevation GPS Undifferenced Data in Different Environments. Remote Sens. 2021, 13, 2078. https://doi.org/10.3390/rs13112078

AMA Style

Liu N, Zhang Q, Zhang S, Wu X. Algorithm for Real-Time Cycle Slip Detection and Repair for Low Elevation GPS Undifferenced Data in Different Environments. Remote Sensing. 2021; 13(11):2078. https://doi.org/10.3390/rs13112078

Chicago/Turabian Style

Liu, Ning, Qin Zhang, Shuangcheng Zhang, and Xiaoli Wu. 2021. "Algorithm for Real-Time Cycle Slip Detection and Repair for Low Elevation GPS Undifferenced Data in Different Environments" Remote Sensing 13, no. 11: 2078. https://doi.org/10.3390/rs13112078

APA Style

Liu, N., Zhang, Q., Zhang, S., & Wu, X. (2021). Algorithm for Real-Time Cycle Slip Detection and Repair for Low Elevation GPS Undifferenced Data in Different Environments. Remote Sensing, 13(11), 2078. https://doi.org/10.3390/rs13112078

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