Next Article in Journal
Validity Evaluation Method Based on Data Driving for On-Line Monitoring Data of Transformer under DC-Bias
Next Article in Special Issue
New Strategy for Improving the Accuracy of Aircraft Positioning Based on GPS SPP Solution
Previous Article in Journal
Smartphone App with an Accelerometer Enhances Patients’ Physical Activity Following Elective Orthopedic Surgery: A Pilot Study
Previous Article in Special Issue
Examination of Multi-Receiver GPS/EGNOS Positioning with Kalman Filtering and Validation Based on CORS Stations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Letter

Precise Method of Ambiguity Initialization for Short Baselines with L1-L5 or E5-E5a GPS/GALILEO Data

by
Mieczysław Bakuła
1,2
1
Faculty of Geoengineering, University of Warmia and Mazury, 10-719 Olsztyn, Poland
2
Institute of Navigation, Military University of Aviation, 08-521 Dęblin, Poland
Sensors 2020, 20(15), 4318; https://doi.org/10.3390/s20154318
Submission received: 13 July 2020 / Revised: 28 July 2020 / Accepted: 30 July 2020 / Published: 2 August 2020
(This article belongs to the Special Issue GNSS Sensors in Aerial Navigation)

Abstract

:
This paper presents a precise and fast method of ambiguity resolution (PREFMAR) for frequencies L1/E1 and L5/E5a of GPS/GALILEO data. The developed method is designed for precise and fast determination of ambiguities in GNSS phase observations. Ambiguities are chosen based on mathematical search functions. The fact that no variance–covariance matrix (VC matrix) with a so-called float solution is needed proves the innovativeness of the developed method. The developed method enables determination of the ambiguities for short baseline double-difference (DD) observations. The presented algorithms for the developed method enable unique and reliable calculation of the ambiguity if the actual errors of code measurements of DD observations are less than 0.38 m and the relative errors of phase observations are in the range of ±3 cm. The paper presents both mathematical derivations of the functions used in the PREFMAR and numerical calculations based on real double-difference GPS observations (L1-L5). The elaborated algorithms can be easily implemented into GNSS receivers or mobile phones. Therefore, they can be widely used in many geoscience applications, as well as in precise GPS/GALILEO navigation.

1. Introduction

GNSS positioning, with centimeter precision, using precise phase observations requires the inclusion of some unknown integer values, the so-called ambiguities, in satellite observations. Double-difference observations using two observation points and two satellites are used in relative positioning. For such an observation system, certain unknown N -values must be determined. These are constant over time if the observations are continuous and have no data gaps. In addition, if two frequencies are used, e.g., L1-L5, ambiguity needs to be determined for each frequency, i.e., N L 1 and N L 5 . The problem of determining integer ambiguities of phase measurements in GPS observations dates back to 1978 [1,2]. Ambiguities are different for each pair of satellites and for each frequency, which poses an extremely complex mathematical and physical problem. Many scientists have worked to find the best and the most effective solution to this problem since the development of GPS systems [3]. In general, when using one or two frequencies, existing calculation methods can be divided into three groups:
  • The first group of methods is based on linear combinations;
  • The second group of methods uses the original search functions in 3D coordinate area;
  • The third group is the most popular and employs the least-squares method.
The greatest step in ambiguity resolution was taken by Hatch–Melbourne–Wubbena [4,5,6] linear combinations, especially by using the so-called wide lane and narrow lane combinations. However, the first approach of ambiguity resolution was based on the use of mathematical functions [7]. Initially, the use of mathematical functions made it possible to achieve centimeter accuracy for vectors with a length of up to 10 km during static sessions lasting 2–3 h [8]. Later studies on mathematical functions improved efficiency of this approach [9,10,11]. The second method only uses the fractional part of the carrier phase measurements to minimize the objective function with the use of the least-squares estimation [12]. The third group is based on the integer least-squares estimation [13,14,15]. Methods of the first and the third group can effectively complement each other, especially in iterative computational strategies based on triple frequencies [16,17,18], in short baseline single-epoch solutions [19,20], and in longer baseline solutions [21].
Most ambiguity solutions start with the determination of estimated coordinates in the global XYZ system in order to determine the approximate ambiguity values, followed by a second and a third stage: estimation and validation. The second stage (estimation) is specific to various methods, as it indicates the most probable sets of ambiguities and is different in each of the ambiguity determination methods. The developed precise and fast method of ambiguity resolution (PREFMAR) uses mathematical functions to choose sets of ambiguities for GNSS observations transmitted on at least two frequencies. Therefore, each frequency combination has its own special properties and different efficiency in the reliable determination of unknown N -values. This research presents a method for determining the most probable ambiguities for GNSS signals transmitted on two frequencies: L1/E1 = 1575.42 MHz and L5/E5a = 1176.45 MHz. Signals of such frequencies are transmitted by GPS satellites, designated as L1 and L5, and GALILEO satellites, designated as E1 and E5a. Therefore, the L1/E1-L5/E5a frequency combination allows the integration of GPS and GALILEO for precise and reliable GNSS positioning [22,23,24].The first part of this contribution presents the mathematical background of the new method of ambiguity estimation based on new functions: Ψ ( N L 1 ) N L 1 N L 5 and Ψ ( N L 5 ) N L 5 N L 1 . The second part of the work provides a detailed numerical example of ambiguity estimation ( N L 1 , N L 5 ) based on real double-difference L1-L5 GPS data.

2. Observation of Double-Difference Equations for GPS Observations on L1 and L5 Frequencies and Their Intercorrelation

For GPS observations transmitted on two frequencies (L1 and L5), we can write the following observation equations for double-difference phase observations ( ϕ L 1 , ϕ L 5 ) and code observations ( P L 1 ,   P L 5 ) for a given measurement epoch ( t ):
λ L 1 ϕ L 1 ( t ) = ϱ ( t ) + λ L 1 N L 1 + λ L 1 ε ϕ L 1 ( t )
P L 1 ( t ) = ϱ ( t ) + ε P L 1 ( t )
λ L 5 ϕ L 5 ( t ) = ϱ ( t ) + λ L 5 N L 5 + λ L 5 ε ϕ L 5 ( t )
P L 5 ( t ) = ϱ ( t ) + ε P L 5 ( t )
where ϕ L 1 , ϕ L 5 are observations of double-difference (DD) phase measurements for L1 and L5 frequencies (in cycles), P L 1 , P L 5 are observations of DD in code measurements (m), ε ϕ L 1 , ε ϕ L 5   are errors of phase DD observations, λ L 1 = 0.190293672798365 (m), and λ L 5 = 0.254828048790854   (m).
There are three unknowns in the above equations: the value of ϱ ( t ) and the sought integral values N L 1 and N L 5 . By separately examining the frequencies L1 and L5 (1-2) and (3-4), we can calculate the ambiguities N L 1 and N L 5 for a single measurement epoch, which are known in the literature as the so-called geometry-free (GF) integer solutions [25]:
N L 1 = [ ϕ L 1 ( t ) P L 1 ( t ) λ L 1 ] roundoff
N L 5 = [ ϕ L 5 ( t ) P L 5 ( t ) λ L 5 ] roundoff
Due to the transmission of signals at two frequencies at the same time, code and phase observations are strongly correlated. Similarly, as in Equation (6), Equation (5) may only be dependent on code measurements at frequency L5, as follows:
N L 1 = [ ϕ L 1 ( t ) P L 5 ( t ) λ L 1 ] roundoff
Instead of value P L 5 ( t ) , P L 1 ( t ) could be used, but practice shows that code observations   P L 5 are much more accurate than P L 1 observations (see Figure 1). Code observations P L 5 are also much more accurate than P L 1 in smartphones [26]. The values of GF ambiguity for N L 1 using code P L 1 (based on Equation (5)) and P L 5 (based on Equation (7)) have been presented in Figure 1. The time series of single-epoch DD geometry-free float solutions presented in Figure 1 have been calculated using real GPS observations, made with Javad GNSS receivers, for a 10-min session with a measurement interval of 1 s.
.
GPS observations, expressed in formulas (6) and (7), are strongly correlated as they depend on the same code observation [27], e.g., observation P L 5 . According to Teunissen’s research, strong correlation can also be obtained when the code P L 5 is replaced by a code P L 1 or by 2 1 ( P L 1 + P L 5 ) . An example of correlated ambiguities for a pair of satellites (G01-G06) and observation ϕ L 1 , ϕ L 5 , P L 5 is shown in Figure 2, both for sessions with a length of five measurement epochs and for sessions of 600 observation epochs.
Figure 2 shows correlated ambiguities N L 1 and N L 5 , using only the P L 5 code. Please note that the float ambiguities are arranged precisely along a straight line, both for 5-s and 10-min sessions.

3. Ambiguity Regression Line Equations

If the correlated float ambiguities are arranged along a straight line, the regression equation for the ambiguity N L 5 in relation to ambiguity N L 1 in the ambiguity system N L 1 N L 5 may be written as
N L 5 = a L 1 , L 5 N L 1 + b L 1 , L 5
or
N ˜ L 5 = a L 1 , L 5 N ˜ L 1 + b L 1 , L 5
where
a L 1 , L 5 = f L 5 f L 1 = λ L 1 λ L 5 = 115 154
As ambiguities in the form of real numbers ( N ˜ L 1 , N ˜ L 5 ) are located along a straight line expressed by (9), the value of b L 1 , L 5 can be calculated based on approximate float values ( N ˜ L 1 , 0 , N ˜ L 5 , 0 ):
b L 1 , L 5 = N ˜ L 5 , 0 115 154 N ˜ L 1 , 0
Thus, taking (11) into account, Equation (8) can be written as follows:
N L 5 = 115 154 ( N L 1 N ˜ L 1 , 0 ) + N ˜ L 5 , 0
If both integer ( N L 5 ) and float ( N ˜ L 5 ) ambiguities lie on the same straight line, the N ˜ L 5 can also be written as
N ˜ L 5 = 115 154 ( N L 1 N ˜ L 1 , 0 ) + N ˜ L 5 , 0
Thus, Equations (12) and (13) represent regression line equations for both real ambiguities (13) and the integer ambiguities (12) and are dependent on approximate values of N ˜ L 1 , 0 and N ˜ L 5 , 0 .

4. Ambiguity Functions for L1-L5 GPS Data

Let us, therefore, define a certain function, Ψ ( N L 1 ) N L 1 N L 5 , as follows:
Ψ ( N L 1 ) N L 1 N L 5 = λ L 5 ( N ˜ L 5 [ N ˜ L 5 ] roundoff )
Equation (14), with Equation (13) taken into account, for the system N L 1 N L 5   can be written as follows:
Ψ ( N L 1 ) N L 1 N L 5 = λ L 5 ( N ˜ L 5 [ N ˜ L 5 ] roundoff ) = λ L 5 ( 115 154 ( N L 1 N ˜ L 1 , 0 ) + N ˜ L 5 , 0 [ 115 154 ( N L 1 N ˜ L 1 , 0 ) + N ˜ L 5 , 0 ] roundoff )
where N ˜ L 1 , 0 and N ˜ L 5 , 0 represent a float solution or are calculated as follows:
N ˜ L 1 , 0 = ϕ L 1 ( t ) P L 5 ( t ) λ L 1
N ˜ L 5 , 0 = ϕ L 5 ( t ) P L 5 ( t ) λ L 5
Similarly, for system N L 5 N L 1 and function Ψ ( N L 5 ) N L 5 N L 1 we obtain
Ψ ( N L 5 ) N L 5 N L 1 = λ L 1 ( N ˜ L 1 [ N ˜ L 1 ] roundoff ) = λ L 1 ( 154 115 ( N L 5 N ˜ L 5 , 0 ) + N ˜ L 1 , 0 [ 154 115 ( N L 5 N ˜ L 5 , 0 ) + N ˜ L 1 , 0 ] roundoff )
Behaviors of the functions | Ψ ( N L 1 ) N L 1 N L 5 | and | Ψ ( N L 5 ) N L 5 N L 1 | have been presented in Figure 3 and Figure 4, where their minima, and thus their periodic character (repeatability), can be clearly seen. The horizontal axes represent integer candidates, but the vertical ones represent values of the ambiguity functions.
Based on an analysis of values of the function Ψ ( N L 1 ) N L 1 N L 5 , in the ( 0.127; + 0.126 m), interval, it can be observed that the values of the function Ψ ( N L 1 ) N L 1 N L 5 for f L 1 and f L 5 frequencies of GPS observations repeat precisely every 154 N L 1 cycles. Values of the function Ψ ( N L 5 ) N L 5 N L 1 are in the interval ( 0.094; + 0.094 m) and repeat precisely every 115 N L 5   cycles. The functions Ψ ( N L 1 ) N L 1 N L 5 and Ψ ( N L 5 ) N L 5 N L 1 have the same wavelength ( λ Ψ = 29.3052256109482 m) and frequency ( f Ψ = 10.23 MHz) and are equivalent in the process of determining ambiguities N L 1 and N L 5 , for values
| ε L 1 , L 5 | 0.5 λ L 1 0.0951468363991824   m
where ε L 1 , L 5 represents both relative errors in phase observations, i.e.,
ε L 1 , L 5 ( t ) = λ L 5 ε ϕ L 5 ( t ) λ L 1 ε ϕ L 1 ( t )
and the value ε L 1 , L 5   meets the relation expressed by Equation (14), i.e.,
ε L 1 , L 5 = Ψ ( N L 1 ) N L 1 N L 5 = λ L 5 ( N ˜ L 5 [ N ˜ L 5 ] roundoff )
Selection of the most likely sets of ambiguities ( N L 1 i , N L 5 i ) in the PREFMAR depends on approximate values of the float solution ( N ˜ L 1 , 0 ,   N ˜ L 5 , 0 ) and relative error values in L1 and L2 carrier phase observations.

5. Ambiguity Search Space in the PREFMAR

Generally, we can assume that the sought integer ambiguities N L 1 and N L 5 lie exactly on the regression line, which additionally needs to pass through an unknown and sought point with coordinates N ^ L 1   and N ^ L 5 . Unfortunately, this is only when (I) ε ϕ L 1 = 0 and ε ϕ L 5 = 0 , and therefore when the errors in phase observations are equal to zero, or (II) when the relative errors in phase observations are equal to zero, i.e., ε L 1 , L 5 = 0 . The first situation corresponds to ideal phase observations, whereas the second one is in the situation when λ L 1 ε ϕ L 1 = λ L 5 ε ϕ L 5 . Detailed analysis of the functions | Ψ ( N L 1 ) N L 1 N L 5 | and | Ψ ( N L 5 ) N L 5 N L 1 | shows that the minima of the functions | Ψ ( N L 1 ) N L 1 N L 5 | and | Ψ ( N L 5 ) N L 5 N L 1 | indicate the sought ambiguity values ( N ^ L 1   and N ^ L 5 ), but only if the errors in phase observations are equal to or close to zero. However, the question that remains is how big the errors in DD phase observations can be for the ambiguities to be determined by the minima of these functions. In the case of phase measurements L1-L5, this depends on the smallest value of the function Ψ ( N L 1 ) N L 1 N L 5 or Ψ ( N L 5 ) N L 5 N L 1 , and these are the values Δ Ψ ( N L 1 ) N L 1 N L 5 and Δ Ψ ( N L 5 ) N L 5 N L 1 , which we calculate as
Δ Ψ ( N L 1 ) N L 1 N L 5 = Δ Ψ ( N L 5 ) N L 5 N L 1 = 1 154 λ L 5 = 1 115 λ L 1 = 0.0016   m
The formula (22), therefore, defines the smallest unit of function Ψ ( N L 1 ) N L 1 N L 5 in the N L 1 N L 5 system and the smallest unit of function Ψ ( N L 5 ) N L 5 N L 1 in the N L 5 N L 1 system. Thus, if the relative errors in DD phase observations are less than half of this unit (22), i.e., less than 0.0008 m, then the minima of the functions | Ψ ( N L 1 ) N L 1 N L 5 | and | Ψ ( N L 5 ) N L 5 N L 1 | determine, on the horizontal axis, the sought ambiguities, with N L 1 using the minimum of the | Ψ ( N L 1 ) N L 1 N L 5 | function and N L 5 using the minimum of the | Ψ ( N L 5 ) N L 5 N L 1 | function. In both cases, approximate values N ˜ L 1 , 0   and N ˜ L 5 , 0   must be at a distance of less than 14.653 m from the true values sought ( N ^ L 1 and N ^ L 5 ).
If, however, the absolute values Δ Ψ ( N L 1 ) N L 1 N L 5 or Δ Ψ ( N L 5 ) N L 5 N L 1 are different from zero and larger than 0.0008 m, then ambiguity values will be located at different points than the minima of the | Ψ ( N L 1 ) N L 1 N L 5 | and | Ψ ( N L 5 ) N L 5 N L 1 | functions, and their selection will strongly depend on the real values ε L 1 , L 5 and on the approximate N ˜ L 1 , 0   and N ˜ L 5 , 0   values.
When analyzing the graph of the function Ψ ( N L 1 ) N L 1 N L 5 in the system N L 1 N L 5   and the graph of the function   Ψ ( N L 5 ) N L 5 N L 1 in the system N L 5 N L 1 , we should start the search for the first pair of ambiguities [ N L 1 I , N L 5 I ] in the area of errors of code measurements below 2 λ L 1 float [ N ˜ L 1 , 0 , N ˜ L 5 , 0 ], i.e., | ( N L 1 I N ˜ L 1 , 0 ) | < 2 λ L 1 ,   and for relative errors of phase observations   ε L 1 , L 5 < 32   m m ; + 32   m m > ; we perform the search for the second pair [ N L 1 I I , N L 5 I I ] in the area   4 λ L 1 , i.e., | ( N L 1 I I N ˜ L 1 , 0 ) | < 4 λ L 1 ; while the third pair [ N L 1 I I I , N L 5 I I I ] of ambiguities must be at a distance of 6 λ L 1 from the value N ˜ L 1 , 0 and N ˜ L 5 , 0 .
Let us, therefore, assume that the values of the function Ψ ( N L 1 ) N L 1 N L 5 lie within the range ±32 mm, then the function Ψ ( N L 1 ) N L 1 N L 5 behaves as shown in Figure 5. Similarly, for the function   Ψ ( N L 5 ) N L 5 N L 1 , we assume a range of values of ±32 mm (Figure 6). Please note that for relative errors in DD phase observations with values up to ± 3 cm, the ambiguity N L 1 changes exactly by four cycles, whereas ambiguity N L 5 changes by three cycles. For comparison, the search areas for   Ψ ( N L 1 ) N L 1 N L 5 and Ψ ( N L 5 ) N L 5 N L 1 which lie within ±62 mm have been presented in Figure 7 and Figure 8.
For relative errors of carrier phases (20) in the range < 32 ; + 32   mm > , these will be the following proposals for N L 1 (Figure 5) and for N L 5 (Figure 6):
N L 1 { 0 ; 4 ; 8 ; 12 ; 16 ; 20 ; 24 ; 28 ; 32 ; 36 ; 39 ; 43 ; 47 ; 51 ; 55 ; 59 ; 63 ; 67 ; 71 ; 75 ; 79 ; 83 ; 87 ; 91 ; 95 ; 99 ;   103 ; 107 ; 111 ; 115 ; 118 ; 122 ; 126 ; 130 ; 134 ; 138 ; 142 ; 146 ; 150 ; 154 } , N L 5 { 0 ; 3 ; 6 ; 9 ; 12 ; 15 ; 18 ; 21 ; 24 ; 27 ; 29 ; 32 ; 35 ; 38 ; 41 ; 44 ; 47 ; 50 ; 53 ; 56 ; 59 ; 62 ; 65 ; 68 ; 71 ; 74 ;   77 ; 80 ; 83 ; 86 ; 88 ; 91 ; 94 ; 97 ; 100 ; 103 ; 106 ; 109 ; 112 ; 115 }
Figure 7 and Figure 8 show a template of possible ambiguities for relative errors in DD phase observations up to ± 62 mm, which seems to be sufficient from a practical point of view even for baselines of several kilometers or much more and for noisy GPS/GALILEO data in mobile phones.

6. Calculation Scheme for N Measurement Epochs Using the PREFMAR

For continuous GNSS observations, and using more than one measurement epoch, general formulas for the PREFMAR are used, which can be summarized in the following main points:
  • Calculation of correlated approximate values of N ˜ L 1 , 0 and N ˜ L 5 , 0 .
    Determination of float value based on the global solution N ˜ L 1 , 0   and N ˜ L 5 , 0 , or using only DD observations for any pair of satellites (for n epochs), i.e.,
    N ˜ L 1 , 0 = n 1 i = 1 n ( ϕ L 1 , i P L 5 λ L 1 )
    N ˜ L 5 , 0 = n 1 i = 1 n ( ϕ L 5 , i P L 5 λ L 5 )
  • Regression line estimation.
    N ˜ L 5 = 115 154 ( N L 1 N ˜ L 1 , 0 ) + N ˜ L 5 , 0
    N ˜ L 1 = 154 115 ( N L 5 N ˜ L 5 , 0 ) + N ˜ L 1 , 0
  • Determination of search functions Ψ ( N L 1 ) N L 1 N L 5 (m) or Ψ ( N L 5 ) N L 5 N L 1 (m).
    Ψ ( N L 1 ) N L 1 N L 2 = = λ L 5 ( 115 154 ( N L 1 N ˜ L 1 , 0 ) + N ˜ L 5 , 0 [ 115 154 ( N L 1 N ˜ L 1 , 0 ) + N ˜ L 5 , 0 ] roundoff )
    Ψ ( N L 5 ) N L 5 N L 1 = = λ L 1 ( 154 115 ( N L 5 N ˜ L 5 , 0 ) + N ˜ L 1 , 0 [ 154 115 ( N L 5 N ˜ L 5 , 0 ) + N ˜ L 1 , 0 ] roundoff )
  • Selection of the most likely sets of ambiguities, depending on approximate values ( N ˜ L 1 , 0 , N ˜ L 5 , 0 ) and using function   Ψ ( N L 1 ) N L 1 N L 5 or function Ψ ( N L 5 ) N L 5 N L 1 .
For short baselines, and assuming that ε L 1 , L 5 < 32 ; + 32   mm > , the four most likely sets of ambiguities ( N L 1 i , N L 5 i ) will be at a distance of up to 1.55 m from the float value and will be selected as follows:
Solution   no   I : = { N L 1 I N L 5 I ,   where   | ( N L 1 I N ˜ L 1 , 0 ) | < 2 λ L 1 ,   ( < 0.381   m )
Solution   no   II : = { N L 1 I I N L 5 I I ,   where   | ( N L 1 I I N ˜ L 1 , 0 ) | < 4 λ L 1 ,   ( < 0.761   m )
Solution   no   III : = { N L 1 I I I N L 5 I I I ,   where   | ( N L 1 I I I N ˜ L 1 , 0 ) | < 6 λ L 1 ,   ( < 1.142   m )
Solution   no   IV : = { N L 1 I V N L 5 I V ,   where   | ( N L 1 I V N ˜ L 1 , 0 ) | < 8 λ L 1 , ( < 1.552   m )
However, if we assume that the relative errors ε L 1 , L 5 < 62 ; + 62   mm > , then solutions I and II remain the same as for ε L 1 , L 5 < 32 ; + 32   mm > , whereas solutions III and IV are at a distance of 3 λ L 1 and 4 λ L 1 from the float value, respectively. Limiting the search range to 4 λ L 1 requires an accuracy of the float position of about 0.76 m if the search range is set for ε L 1 , L 5 < 62 ; + 62   mm > . However, given the number of satellites currently available, obtaining an accuracy of the float solution better than 0.76 m should not be a problem, which allows a quick indication of the most likely ambiguities for validation, even for GNSS measurements under difficult observational conditions or for mobile phones. However, based on preliminary tests, for high-quality GNSS receivers, the search area of ambiguity resolution N L 1 and N L 5 for relative errors ε L 1 , L 5 < 32 ; + 32   mm > should be sufficient for short baseline RTK positioning.

7. Numerical Example Using the PREFMAR and Float Solution

The formulas presented for the developed method were applied to real GPS data, for a very short baseline, based on DD observations for the G01–G06 satellite pair, using five observation epochs. The DD observations for these satellites are shown in Table 1.
Sets of the most probable ambiguities ( N L 1 , N L 5 ) can be determined both with the use of function Ψ ( N L 1 ) N L 1 N L 5 and function Ψ ( N L 5 ) N L 5 N L 1 . Let us thus present the necessary calculations in the form of Table 2 for the function Ψ ( N L 1 ) N L 1 N L 5 and using 5 DD observations for which the float values are respectively equal to N ˜ L 1 , 0 = 5.141 and N ˜ L 5 , 0 = 8.349 . Then, for subsequent integer values N L 1 , located around the float solution N ˜ L 1 , 0 , within the range of, e.g., ± 8 N L 1 cycles (i.e., at a distance of 1.5 m) from approximate values of N ˜ L 1 , 0 = 5.141, we calculate elements of the column N ˜ L 5 , i using the formula
N ˜ L 5 , i = 115 154 N L 1 , i + 8.349 115 156 5.141
and the values of the Ψ ( N L 1 ) N L 1 N L 5 function are calculated with the formula
Ψ ( N L 1 ) N L 1 N L 5 = λ L 5 ( 115 154 N L 1 , i + 8.349 115 154 5.141 [ 115 154 N L 1 , i + 8.349 115 154 5.141 ] roundoff )
Similarly, for the function Ψ ( N L 5 ) N L 5 N L 1 , the calculations are shown in Table 3. The values of the function Ψ ( N L 5 ) N L 5 N L 1 are calculated from the following formula:
Ψ ( N L 5 ) N L 5 N L 1 = λ L 1 ( 154 115 N L 5 + 5.141 154 115 8.349 [ 154 115 N L 5 + 5.141 154 115 8.349 ] roundoff )
Identical sets of the most probable ambiguities have been obtained based on Table 2 and Table 3, which proves the reliability of the mathematical functions Ψ ( N L 1 ) N L 1 N L 5 and Ψ ( N L 5 ) N L 5 N L 1 used in the PREFMAR. The lowest values of the search functions were obtained for [ N L 1 = 2 ;   N L 5 = 6 ] ; the second proposal was [ N L 1 = 6 ;   N L 5 = 9 ] , and the third proposal was [ N L 1 = 10 ;   N L 5 = 12 ] . Similar results are obtained using the data from Table 1, for any of the individual DD epochs. The information contained in Table 2 and Table 3 can also be presented in graphical form, as shown in Figure 9 and Figure 10.
However, as already mentioned before, taking into account the accuracy of the code measurements, the most probable pair of ambiguities must be selected in a distance from a float solution below 2 λ L 1 (i.e.,   N L 1 I = 6 ;   N L 5 I = 9 ); the search for the second pair is performed in the area 4 λ L 1 from float solution (i.e.,   N L 1 I I = 2 ;   N L 5 I I = 6 ), while the third pair of ambiguities must be at a distance 6 λ L 1 (i.e.,   N L 1 I I I = 10 ;   N L 5 I I I = 12 ). However, for the relative errors up to ±62 mm, the first four most probable pair of ambiguities will be as follows: [ N L 1 I = 6 ;   N L 5 I = 9 ], [ N L 1 I I = 2 ;   N L 5 I I = 6 ], [ N L 1 I I I = 5 ;   N L 5 I I I = 8 ], and [ N L 1 I V = 9 ;   N L 5 I V = 11 ].

8. Discussion

Generally, the search area in the PREFMAR is represented by a parallelogram. The orientation of this parallelogram in a given system of ambiguities is only dependent on the frequency of the GNSS signals transmitted. For system N L 1 N L 5 , the short sides of this parallelogram are parallel to the vertical lines of the system N L 1 N L 5 , i.e., to the N L 5 axis (Figure 9). The longer side of this parallelogram is slanted with respect to the N L 5 axis under an angle of α L 1 , L 5 = t a n 1 ( 115 / 154 ) = 36.750655 o, and the geometrical center of the search area is located at point ( N ˜ L 1 , 0 , N ˜ L 5 , 0 ) . Similarly, for the system N L 5 N L 1 (Figure 10), α L 5 , L 1 = t a n 1 ( 154 / 115 ) = 53.249344 o.
Please note that the presented functions fulfil the relation Ψ ( N L 1 ) N L 1 N L 5 = Ψ ( N L 5 ) N L 5 N L 1 because the values of the functions represent relative errors of phase observations between L1-L5, which can be additionally written in a simplified manner:
Ψ ( N L 1 ) N L 1 N L 5 = [ λ L 5 ϕ L 5 ( t ) λ L 5 N L 5 ] [ λ L 1 ϕ L 1 ( t ) λ L 1 N L 1 ]
Ψ ( N L 5 ) N L 5 N L 1 = [ λ L 1 ϕ L 1 ( t ) λ L 1 N L 1 ] [ λ L 5 ϕ L 5 ( t ) λ L 5 N L 5 ]
The formulas presented above also allow indicating integer ambiguities by substituting appropriate sets of ambiguities N L 1 and N L 5 , near the float ambiguity, e.g., for the case presented above using 5 DD epochs, we have ϕ L 1 ( t ) = 42.689 and ϕ L 5 ( t ) = 36.388 thus for [ N L 1 = 6 ;   N L 5 = 9 ] ; based on formulas (36) and (37) we obtain:
Ψ ( N L 1 ) N L 1 N L 5 = [ λ L 5 36.388 λ L 5 9 ] [ λ L 1 42.689 λ L 1 6 ] = 0.002   m
Ψ ( N L 5 ) N L 5 N L 1 = [ λ L 1 42.689 λ L 1 6 ] [ λ L 5 36.388 λ L 5 9 ] = 0.002   m
which is equivalent to the values of the functions Ψ ( N L 1 ) N L 1 N L 5 and Ψ ( N L 5 ) N L 5 N L 1 obtained based on formulas (16) and (19) and presented in Table 2 and Table 3, respectively.
Presented equations of ambiguity functions can be used for different double or triple frequencies; however, they have specific properties in every case. Therefore, we are going to present the most popular combinations of GNSS frequencies in the near future to show the power of the PREFMAR.
It should be emphasized that the PREFMAR gives results based on GNSS measurements from only two available GNSS satellites (see Section 7); therefore, we can use this method effectively in challenging observational satellite conditions.
Furthermore, to obtain the most accurate float solutions, we can effectively use the Kalman filter in differential [28] or relative code GNSS positioning.

9. Summary and Conclusions

This research presents a precise and fast method of ambiguity resolution (PREFMAR) for indicating the most probable ambiguities for GNSS observations transmitted on two frequencies, i.e., L1/E1 = 1575.42 MHz and L5/E5a = 1176.45 MHz. Signals with such frequencies are transmitted by satellites of the American GPS system and the European GALILEO system. Therefore, the combination of these frequencies allows the precise integration of the various GNSS positioning systems. The described method is used to determine ambiguities using only one observation epoch immediately. This PREFMAR is intended for both short and long baselines. Its efficiency mostly depends on the values of relative errors in DD phase observations and on the precision of code measurements. Ambiguity is chosen based on mathematical functions Ψ ( N L 1 ) N L 1 N L 5 or Ψ ( N L 5 ) N L 5 N L 1 , using a correlation between the value of a precisely defined ambiguity and dependent only on the frequency of satellite signals. The PREFMAR allows the ambiguities for single measurement epochs to be determined without using the VC matrix from the float solution and also allows the most probable ambiguities to be indicated even for individual DD observations. The presented mathematical functions can also be used for precise and immediate re-initialization of the ambiguities. Additionally, the contribution presents interpretations of the derived functions Ψ ( N L 1 ) N L 1 N L 5 and Ψ ( N L 5 ) N L 5 N L 1 in relation to the relative errors in phase observations. As a result, a simplified calculation method is included for values of the functions Ψ ( N L 1 ) N L 1 N L 5   and   Ψ ( N L 5 ) N L 5 N L 1 , which are equivalent to the mathematically detailed versions, but only in the range 0.5 λ L 5 for function Ψ ( N L 1 ) N L 1 N L 5 and in the range 0.5 λ L 1 for function Ψ ( N L 5 ) N L 5 N L 1 .

Funding

This paper was funded by the University of Warmia and Mazury.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hofmann-Wellenhof, B.; Lichtenegger, H.; Collins, J. GNSS—Global Navigation Satellite Systems: GPS, GLONASS, Galileo, and More; Springer: New York, NY, USA, 2008. [Google Scholar]
  2. Leick, A.; Rapoport, L.; Tatarnikov, D. GPS Satellite Surveying, 3rd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015. [Google Scholar]
  3. Kim, D.; Langley, R.B. GPS ambiguity resolution and validation: Methodologies, trends and issues. In Proceedings of the 7th GNSS Workshop International Symposium on GPS/GNSS, Seoul, Korea, 30 November–2 December 2000; pp. 213–221. [Google Scholar]
  4. Hatch, R.R. The synergism of GPS code and carrier measurements. In Proceedings of the Third International Geodetic Symposium on Satellite Doppler Positioning, Las Cruces, NM, USA, 8–12 February 1982; Volume 2, pp. 1213–1231. [Google Scholar]
  5. Melbourne, W.G. The case for ranging in GPS based geodetic systems. In Proceedings of the 1st International Symposium on Precise Positioning with the GPS, Rockville, MD, USA, 15–19 April 1985; pp. 373–386. [Google Scholar]
  6. Wubbena, G. Software developments for geodetic positioning with GPS using TI4100 code and carrier measurements. In Proceedings of the 1st International Symposium on Precise Positioning with the GPS, Rockville, MD, USA, 15–19 April 1985; pp. 403–412. [Google Scholar]
  7. Counselman, C.C.; Gourevitch, S.A. Miniature interferometer terminals for earth surveying: Ambiguity and multipath with the global positioning system. IEEE Trans. Geosci. Remote Sens. 1981, 19, 244–252. [Google Scholar] [CrossRef]
  8. Counselman, C.C.; Abbot, R.I.; Gourevitch, S.A.; King, R.W.; Paradis, A.R. Centimeter-level relative positioning with GPS. J. Surv. Eng. 1983, 109, 81–89. [Google Scholar] [CrossRef]
  9. Remondi, B.W. Pseudo-kinematic GPS Results Using the Ambiguity Function Method. Navigation 1991, 38, 17–36. [Google Scholar] [CrossRef]
  10. Han, S.; Rizos, C. Improving the computational efficiency of the ambiguity function algorithm. J. Geod. 1996, 70, 330–341. [Google Scholar] [CrossRef]
  11. Baselga, S. Global optimization applied to GPS positioning by ambiguity functions. Meas. Sci. Technol. 2010, 21, 125102. [Google Scholar] [CrossRef]
  12. Cellmer, S. Search procedure for improving modified ambiguity function approach. Surv. Rev. 2013, 45, 380–385. [Google Scholar] [CrossRef]
  13. Frei, E.; Beutler, G. Rapid static positioning based on the fast ambiguity resolution approach “FARA”: Theory and first results. Manus. Geod. 1990, 15, 325–356. [Google Scholar]
  14. Teunissen, P.J.G. The least-squares ambiguity decorrelation adjustment: A method for fast GPS integer ambiguity estimation. J. Geod. 1995, 70, 65–82. [Google Scholar] [CrossRef]
  15. Chang, X.W.; Yang, X.; Zhou, T. MLAMBDA: A modified LAMBDA method for integer least-squares estimation. J. Geod. 2005, 79, 552–565. [Google Scholar] [CrossRef] [Green Version]
  16. Hatch, R.; Jung, J.; Enge, P.; Pervan, B. Civilian GPS: The benefits of three frequencies. GPS Solut. 2000, 3, 1–9. [Google Scholar] [CrossRef]
  17. Wang, X.; Liu, W.; Sun, G. An Improved Geometry-free Three Carrier Ambiguity Resolution Method for the BeiDou Navigation Satellite System. J. Navig. 2016, 69, 1393–1408. [Google Scholar] [CrossRef]
  18. Wang, S.; Deng, J.; Lu, X.; Song, Z.; Xu, Y. A new GNSS single-epoch ambiguity resolution method based on triple-frequency signals. ISPRS Int. J. Geo Inf. 2017, 6, 46. [Google Scholar] [CrossRef] [Green Version]
  19. Deng, C.; Tang, W.; Liu, J.; Shi, C. Reliable single-epoch ambiguity resolution for short baselines using combined GPS/BeiDou system. GPS Solut. 2014, 18, 375–386. [Google Scholar] [CrossRef]
  20. Wu, Z.; Bian, S.; Ji, B.; Xiang, C.; Jiang, D. Short baseline GPS multi-frequency single-epoch precise positioning: Utilizing a new carrier phase combination method. J. Geod. 2016, 20, 373–384. [Google Scholar] [CrossRef]
  21. Xu, L.; Liu, H.; Shu, B.; Zheng, F.; Zhang, M.; Qian, C.; Duan, Y. GLONASS real-time wide-lane ambiguity resolution with an enhanced geometry-based model for medium-range baselines. Adv. Space Res. 2018, 62, 2467–2479. [Google Scholar] [CrossRef]
  22. Tiberus, C.; Pany, T.; Eissfeller, B.; Joosen, P.; Verhagen, S. 0.99999999 confidence ambiguity resolution with GPS and Galileo. GPS Solut. 2002, 6, 96–99. [Google Scholar] [CrossRef]
  23. Cellmer, S.; Paziewski, J.; Wielgosz, P. Fast and precise positioning using MAFA method and new GPS and Galileo signals. Acta Geodyn. Geomater. 2013, 10, 393–400. [Google Scholar] [CrossRef] [Green Version]
  24. Bolla, P.; Borre, K. Performance analysis of dual-frequency receiver using combinations of GPS L1, L5, and L2 civil signals. J. Geod. 2019, 93, 437–447. [Google Scholar] [CrossRef]
  25. Misra, P.; Enge, P. Global Positioning System. Signals, Measurements, and Performance; Ganga-Jamuna Press: Nanded, India, 2006. [Google Scholar]
  26. Uradziński, M.; Bakuła, M. Assessment of static positioning accuracy using low-cost smartphone GPS devices for geodetic survey points’ determination and monitoring. Appl. Sci. 2020, 10, 5308. [Google Scholar]
  27. Teunissen, P.J.G. An analytical study of ambiguity decorrelation using dual frequency code and carrier phase. J. Geod. 1996, 70, 515–528. [Google Scholar] [CrossRef]
  28. Ciećko, A.; Bakuła, M.; Grunwald, G.; Ćwiklak, J. Examination of Multi-Receiver GPS/EGNOS Positioning with Kalman Filtering and Validation Based on CORS Stations. Sensors 2020, 20, 2732. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Time series of single-epoch double-difference (DD) geometry-free float solutions N L 1 with the use of P L 1 and P L 5 .
Figure 1. Time series of single-epoch double-difference (DD) geometry-free float solutions N L 1 with the use of P L 1 and P L 5 .
Sensors 20 04318 g001
Figure 2. Scatter plot of correlated single-epoch geometry-free ambiguities N L 1 and N L 5 .
Figure 2. Scatter plot of correlated single-epoch geometry-free ambiguities N L 1 and N L 5 .
Sensors 20 04318 g002
Figure 3. Behavior of a full period of the | Ψ ( N L 1 ) N L 1 N L 5 | function in the N L 1 N L 5 system, in the N L 1 < 0 ; 154 > interval, for GNSS observations and frequencies L1/E1 and L5/E5a.
Figure 3. Behavior of a full period of the | Ψ ( N L 1 ) N L 1 N L 5 | function in the N L 1 N L 5 system, in the N L 1 < 0 ; 154 > interval, for GNSS observations and frequencies L1/E1 and L5/E5a.
Sensors 20 04318 g003
Figure 4. Behavior of a full period of the function | Ψ ( N L 5 ) N L 5 N L 1 | in the N L 5 N L 1 system, in the N L 5 < 0 ; 115 > interval, for GNSS observations and frequencies L1/E1 and L5/E5a.
Figure 4. Behavior of a full period of the function | Ψ ( N L 5 ) N L 5 N L 1 | in the N L 5 N L 1 system, in the N L 5 < 0 ; 115 > interval, for GNSS observations and frequencies L1/E1 and L5/E5a.
Sensors 20 04318 g004
Figure 5. Behavior of the function Ψ ( N L 1 ) N L 1 N L 5 (m), for values from interval < 32 ; + 32   mm > .
Figure 5. Behavior of the function Ψ ( N L 1 ) N L 1 N L 5 (m), for values from interval < 32 ; + 32   mm > .
Sensors 20 04318 g005
Figure 6. Behavior of the function Ψ ( N L 5 ) N L 5 N L 1 (m), for values from the interval < 32 ; + 32   mm > .
Figure 6. Behavior of the function Ψ ( N L 5 ) N L 5 N L 1 (m), for values from the interval < 32 ; + 32   mm > .
Sensors 20 04318 g006
Figure 7. Behavior of the function Ψ ( N L 1 ) N L 1 N L 5 (m), for values from interval < 62 ; + 62   mm > .
Figure 7. Behavior of the function Ψ ( N L 1 ) N L 1 N L 5 (m), for values from interval < 62 ; + 62   mm > .
Sensors 20 04318 g007
Figure 8. Behavior of the function Ψ ( N L 5 ) N L 5 N L 1 (m), for values from interval < 62 ; + 62   mm > .
Figure 8. Behavior of the function Ψ ( N L 5 ) N L 5 N L 1 (m), for values from interval < 62 ; + 62   mm > .
Sensors 20 04318 g008
Figure 9. Search area of ambiguities in N L 1 N L 5 system, for N ˜ L 1 , 0 = 5.141 and N ˜ L 5 , 0 = 8.349 , with the use of the function | Ψ ( N L 1 ) N L 1 N L 5 | , for | ε L 1 , L 5 | < 5.5   cm .
Figure 9. Search area of ambiguities in N L 1 N L 5 system, for N ˜ L 1 , 0 = 5.141 and N ˜ L 5 , 0 = 8.349 , with the use of the function | Ψ ( N L 1 ) N L 1 N L 5 | , for | ε L 1 , L 5 | < 5.5   cm .
Sensors 20 04318 g009
Figure 10. Search area of ambiguities in N L 5 N L 1 system, for N ˜ L 5 , 0 = 8.349 and N ˜ L 1 , 0 = 5.141 , with the use of the function | Ψ ( N L 5 ) N L 5 N L 1 | , for | ε L 5 , L 1 | < 5.5   cm .
Figure 10. Search area of ambiguities in N L 5 N L 1 system, for N ˜ L 5 , 0 = 8.349 and N ˜ L 1 , 0 = 5.141 , with the use of the function | Ψ ( N L 5 ) N L 5 N L 1 | , for | ε L 5 , L 1 | < 5.5   cm .
Sensors 20 04318 g010
Table 1. DD observations for the G01–G06 satellite pair.
Table 1. DD observations for the G01–G06 satellite pair.
EpochC1 (m)L1 (c)P5 (m)L5 (c)NL1 (c)NL5 (c)
17.43742.6866.96736.3906.0749.050
27.09342.7257.07336.4075.5568.651
35.40242.6547.13236.3655.1758.378
47.39342.6997.33336.3944.1647.618
53.77142.6817.22136.3844.7348.047
Average6.21942.6897.14536.3885.1418.349
Table 2. The set of integer candidates based on the function Ψ ( N L 1 ) N L 1 N L 2 , for N ˜ L 1 , 0 = 5.141 and N ˜ L 5 , 0 = 8.349 , where N L 1 , i < 3 ; 13 > .
Table 2. The set of integer candidates based on the function Ψ ( N L 1 ) N L 1 N L 2 , for N ˜ L 1 , 0 = 5.141 and N ˜ L 5 , 0 = 8.349 , where N L 1 , i < 3 ; 13 > .
N L 1 N ˜ L 5 Ψ ( N L 1 ) N L 1 N L 5   ( m ) Solution No. N L 5 = [ N ˜ L 5 ] r o u n d o f f
−32.2700.069 2
−23.0160.004IV3
−13.763−0.060 4
04.510−0.125 5
15.2570.065 5
26.0030.001II6
36.750−0.064 7
47.4970.127 7
58.2440.062 8
68.990−0.002I9
79.737−0.067 10
810.4840.123 10
911.2310.059 11
1011.977−0.006III12
1112.724−0.070 13
1213.4710.120 13
1314.2180.055 14
Table 3. Set of integer candidates based on the function Ψ ( N L 5 ) N L 5 N L 1 , for N ˜ L 5 , 0 = 8.349 and N ˜ L 1 , 0 = 5.141 , where N L 5 < 2 ; 14 > .
Table 3. Set of integer candidates based on the function Ψ ( N L 5 ) N L 5 N L 1 , for N ˜ L 5 , 0 = 8.349 and N ˜ L 1 , 0 = 5.141 , where N L 5 < 2 ; 14 > .
N L 5 N ˜ L 1 Ψ ( N L 5 ) N L 5 N L 1   ( m ) Solution No. N L 1 = [ N ˜ L 1 ] r o u n d o f f
2−3.361−0.069 −3
3−2.022−0.004IV−2
4−0.6830.060 −1
50.656−0.065 1
61.995−0.001II2
73.3350.064 3
84.674−0.062 5
96.0130.002I6
107.3520.067 7
118.691−0.059 9
1210.0300.006III10
1311.3690.070 11
1412.708−0.055 13

Share and Cite

MDPI and ACS Style

Bakuła, M. Precise Method of Ambiguity Initialization for Short Baselines with L1-L5 or E5-E5a GPS/GALILEO Data. Sensors 2020, 20, 4318. https://doi.org/10.3390/s20154318

AMA Style

Bakuła M. Precise Method of Ambiguity Initialization for Short Baselines with L1-L5 or E5-E5a GPS/GALILEO Data. Sensors. 2020; 20(15):4318. https://doi.org/10.3390/s20154318

Chicago/Turabian Style

Bakuła, Mieczysław. 2020. "Precise Method of Ambiguity Initialization for Short Baselines with L1-L5 or E5-E5a GPS/GALILEO Data" Sensors 20, no. 15: 4318. https://doi.org/10.3390/s20154318

APA Style

Bakuła, M. (2020). Precise Method of Ambiguity Initialization for Short Baselines with L1-L5 or E5-E5a GPS/GALILEO Data. Sensors, 20(15), 4318. https://doi.org/10.3390/s20154318

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