Jump to content

GNSS enhancement: Difference between revisions

From Wikipedia, the free encyclopedia
Content deleted Content added
No edit summary
m Clean up spacing around commas and other punctuation fixes, replaced: ,j → , j (18), ,k → , k (18), ,m → , m (8), ,n → , n (8)
 
(27 intermediate revisions by 18 users not shown)
Line 1: Line 1:
{{short description|GPS accuracy improvement techniques}}
{{Expert needed|1=Geography|reason=possibly duplicate references to the same phase stuff, esp. between RTK and carrier-phase tracking|date=March 2020}}
'''GNSS enhancement''' refers to techniques used to improve the accuracy of positioning information provided by the [[Global Positioning System]] or other [[global navigation satellite systems]] in general, a network of satellites used for navigation.
'''GNSS enhancement''' refers to techniques used to improve the accuracy of positioning information provided by the [[Global Positioning System]] or other [[global navigation satellite systems]] in general, a network of satellites used for navigation.
Enhancement methods of improving accuracy rely on external information being integrated into the calculation process. There are many such systems in place and they are generally named or described based on how the GPS sensor receives the information. Some systems transmit additional information about sources of error (such as clock drift, ephemeris, or [[ionospheric delay]]), others provide direct measurements of how much the signal was off in the past, while a third group provides additional navigational or vehicle information to be integrated into the calculation process.
{{duplication|dupe=Global Positioning Systems|date=February 2016}}
Enhancement methods of improving accuracy rely on external information being integrated into the calculation process. There are many such systems in place and they are generally named or described based on how the GPS sensor receives the information. Some systems transmit additional information about sources of error (such as clock drift, ephemeris, or [[ionospheric delay]]), others provide direct measurements of how much the signal was off in the past, while a third group provide additional navigational or vehicle information to be integrated in the calculation process.

Examples of augmentation systems include the [[Wide Area Augmentation System]], [[Differential GPS]], [[Inertial Navigation System]]s and [[Assisted GPS]].


== Background ==
== Background ==
The [[Global Positioning System]] (GPS) is the American satellite-based system for positioning and navigation. Receivers on or near the earth's surface can determine their locations based on signals received from any four or more of the satellites in the network.
The [[Global Positioning System]] (GPS) is the American satellite-based system for positioning and navigation. Receivers on or near the Earth's surface can determine their locations based on signals received from any four or more of the satellites in the network.


All GPS satellites broadcast on the same two frequencies, known as L1 (1575.42 MHz) and L2 (1227.60 MHz). The network uses [[code division multiple access]] (CDMA) to allow separate messages from the individual satellites to be distinguished. Two distinct CDMA encodings are used: the coarse/acquisition (C/A) code, which is accessible by the general public, and the precise (P) code, that is encrypted so that only the U.S. military can access it. The messages sent from each satellite contain information ranging from the satellite health, the satellite's orbital path, the clock state of the satellite, and the configuration of the entire satellite network.
All GPS satellites broadcast on the same two frequencies, known as L1 (1575.42 MHz) and L2 (1227.60 MHz). The network uses [[code-division multiple access]] (CDMA) to allow separate messages from the individual satellites to be distinguished. Two distinct CDMA encodings are used: the coarse/acquisition (C/A) code, which is accessible by the general public, and the precise (P) code, which is encrypted, so that only the U.S. military can access it. The messages sent from each satellite contain information ranging from the satellite's health, the satellite's orbital path, the clock state of the satellite, and the configuration of the entire satellite network.


== Precise monitoring ==
== Precise monitoring ==
The accuracy of a calculation can also be improved through precise monitoring and measuring of the existing GPS signals in additional or alternate ways.
The accuracy of a calculation can also be improved through precise monitoring and measuring of the existing GPS signals in additional or alternate ways.


After ''Selective Availability'' was turned off by the U.S. government, the largest error in GPS was usually the unpredictable delay through the ionosphere. The spacecraft broadcast ionospheric model parameters, but they are necessarily imperfect. This is one reason the GPS spacecraft transmit on at least two frequencies, L1 and L2. [[Error analysis for the Global Positioning System#Atmospheric effects|Ionospheric delay]] is a well-defined function of frequency and the [[total electron content]] (TEC) along the path, so measuring the arrival time difference between the frequencies determines TEC and thus the precise ionospheric delay at each frequency.
After ''Selective Availability'' was turned off by the U.S. government, the largest error in GPS was usually the unpredictable delay through the ionosphere. The spacecraft broadcast ionospheric model parameters, but they are necessarily imperfect. This is one reason the GPS spacecraft transmit on at least two frequencies, L1 and L2. [[Error analysis for the Global Positioning System#Atmospheric effects|Ionospheric delay]] is a well-defined function of frequency and the [[total electron content]] (TEC) along the path, so measuring the arrival-time difference between the frequencies determines TEC and thus the precise ionospheric delay at each frequency.


Receivers with decryption keys can decode the P(Y)-code transmitted on both L1 and L2. However, these keys are reserved for the military and authorized agencies and are not available to the public. Without keys, it is still possible to use a ''codeless'' technique to compare the P(Y) codes on L1 and L2 to gain much of the same error information. However, this technique is slow, so it is currently limited to specialized surveying equipment. In the future, additional civilian codes are expected to be transmitted on the L2 and L5 frequencies (see [[GPS modernization]]). Then all users will be able to perform dual-frequency measurements and directly compute ionospheric delay errors.
Receivers with decryption keys can decode the P(Y)-code transmitted on both L1 and L2. However, these keys are reserved for the military and authorized agencies and are unavailable to the public. Without keys, it is still possible to use a ''codeless'' technique to compare the P(Y) codes on L1 and L2 to gain much of the same error information. However, this technique is slow and is currently limited to specialized surveying equipment. In the future, additional civilian codes are expected to be transmitted on the L2 and L5 frequencies (see [[GPS modernization]]). Then all users can perform dual-frequency measurements and directly compute ionospheric-delay errors.


A second form of precise monitoring is called ''Carrier-Phase Enhancement'' (CPGPS). <ref name="ubiquitous" /> The error, which this corrects, arises because the pulse transition of the [[Pseudorandom number generator|PRN]] is not instantaneous, and thus the [[cross-correlation|correlation]] (satellite-receiver sequence matching) operation is imperfect. The CPGPS approach utilizes the L1 carrier wave, which has a [[frequency|period]] of
A second form of precise monitoring is called ''carrier-phase enhancement'' (CPGPS).<ref name="ubiquitous" /> The error, which this corrects, arises because the pulse transition of the [[Pseudorandom number generator|PRN]] is not instantaneous, and thus the [[cross-correlation|correlation]] (satellite–receiver sequence matching) operation is imperfect. The CPGPS approach utilizes the L1 carrier wave, which has a [[frequency|period]] of


:<math> \frac{1 \ \mathrm{s}}{1575.42 \times 10^6} = 0.63475 \ \mathrm{ns} \approx 1 \ \mathrm{ns} \ </math>
:<math> \frac{1~\text{s}}{1575.42 \times 10^6} = 0.63475~\text{ns} \approx 1~\text{ns}, </math>


which is about one-thousandth of the C/A Gold code bit period of
which is about one-thousandth of the C/A Gold code bit period of


:<math> \frac{1 \ \mathrm{s}}{1023 \times 10^3} = 977.5 \ \mathrm{ns} \ \approx 1000 \ \mathrm{ns} \ </math>
:<math> \frac{1~\text{s}}{1023 \times 10^3} = 977.5~\text{ns} \approx 1000~\text{ns}, </math>

to act as an additional [[clock signal]] and resolve the uncertainty. The phase difference error in the normal GPS amounts to between 2 and 3 meters (6 to 10&nbsp;ft) of ambiguity. CPGPS working to within 1% of perfect transition reduces this error to 3 centimeters (1&nbsp;inch) of ambiguity. By eliminating this source of error, CPGPS coupled with [[Differential GPS|DGPS]] normally realizes between 20 and 30 centimeters (8 to 12&nbsp;inches) of absolute accuracy.


to act as an additional [[clock signal]] and resolve the uncertainty. The phase-difference error in the normal GPS amounts to between 2 and 3 meters (6 to 10&nbsp;ft) of ambiguity. CPGPS working to within 1% of perfect transition reduces this error to 3 centimeters (1&nbsp;inch) of ambiguity. By eliminating this source of error, CPGPS coupled with [[Differential GPS|DGPS]] normally realizes between 20 and 30 centimeters (8 to 12&nbsp;inches) of absolute accuracy.
== Real-time Kinematic Positioning ==
{{main article|Real Time Kinematic}}
''Real-time Kinematic Positioning'' (RKP) is another approach for a precise GPS-based positioning system. In this approach, determination of range signal can be resolved to a precision of less than 10 [[centimeters]] (4&nbsp;in). This is done by resolving the number of cycles in which the signal is transmitted and received by the receiver. This can be accomplished by using a combination of differential GPS (DGPS) correction data, transmitting GPS signal phase information and ambiguity resolution techniques via statistical tests—possibly with processing in real-time.


== Real-time kinematic positioning ==
== Timekeeping {{anchor|GPS time|GPS time and date}} ==
{{main article|Real-time kinematic}}
While most clocks are synchronized to [[Coordinated Universal Time]] (UTC), the [[atomic clock]]s on the satellites are set to ''GPS time''. The difference is that GPS time is not corrected to match the rotation of the Earth, so it does not contain [[leap second]]s or other corrections which are periodically added to UTC. GPS time was set to match [[Coordinated Universal Time]] (UTC) in 1980, but has since diverged. The lack of corrections means that GPS time remains at a constant offset with [[International Atomic Time]] (TAI) (TAI - GPS = 19 seconds). Periodic corrections are performed on the on-board clocks to correct relativistic effects and keep them synchronized with ground clocks.
''Real-time kinematic positioning'' (RTK) is another approach for a precise GPS-based positioning system. In this approach, the determination of the range signal can be resolved to a precision of less than 10 [[centimeters]] (4&nbsp;in). This is done by resolving the number of cycles in which the signal is transmitted and received by the receiver. This can be accomplished by using a combination of differential GPS (DGPS) correction data, transmitting GPS signal phase information, and ambiguity resolution techniques via statistical tests, possibly with processing in real time.


== Carrier-phase tracking (surveying) ==
The GPS navigation message includes the difference between GPS time and UTC, which as of 2016 is 17 seconds due to leap seconds added since 1980. Receivers subtract this offset from GPS time to calculate UTC and specific timezone values. The GPS-UTC offset, a signed eight-bit field,<ref>[http://www.gps.gov/technical/icwg/IS-GPS-200G.pdf Global Positioning Systems Directorate Systems Engineering & Integration Interface Specification p. 114]</ref> can accommodate up to 128 negative leaps or 127 positive leap seconds. At the current rate of Earth's rotation (with one positive leap second introduced approximately every 18 months), the offset field is sufficient for approximately one more century. However, as the Earth is not only rotating slower than when UTC was adopted in 1972, but also [[Coordinated Universal Time#Future|slowing down]] at a somewhat unpredictable rate, leap seconds will begin to occur more frequently. So, the GPS offset field may [[Integer overflow|overflow]] within the 21st century.
Utilizing the navigation message to measure pseudo-range has been discussed. Another method that is used in GPS surveying applications is carrier-phase tracking. The period of the carrier frequency times the speed of light gives the wavelength, which is about 0.19 meters for the L1 carrier. With a 1% of wavelength accuracy in detecting the leading edge, this component of pseudo-range error might be as low as 2 millimeters. This compares to 3 meters for the C/A code and 0.3 meters for the P code.


However, this 2-millimeter accuracy requires measuring the total phase, the total number of wavelengths, plus the fractional wavelength. This requires specially equipped receivers. This method has many applications in the field of surveying.
As opposed to the year, month, and day format of the [[Gregorian calendar]], the GPS date is expressed as a week number and a seconds-into-week number. The week number is transmitted as a ten-[[bit]] field in the C/A and P(Y) navigation messages, and so it becomes zero again every 1,024 weeks (19.6 years). GPS week zero started at 00:00:00 UTC (00:00:19 TAI) on January 6, 1980, and the week number became zero again for the first time at 23:59:47 UTC on August 21, 1999 (00:00:19 TAI on August 22, 1999). To determine the current Gregorian date, a GPS receiver must be provided with the approximate date (to within 3,584 days) to correctly translate the GPS date signal. To address this concern the modernized GPS navigation message uses a 13-bit field, which only repeats every 8,192 weeks (157 years), thus lasting until the year 2137 (157 years after GPS week zero).


We now describe a method that could potentially be used to estimate the position of receiver 2 given the position of receiver 1 using triple differencing followed by numerical root finding and a mathematical technique called [[least squares]]. A detailed discussion of the errors is omitted in order to avoid detracting from the description of the methodology. In this description, differences are taken in the order of differencing between satellites, differencing between receivers, and differencing between epochs. This should not be construed to mean that this is the only order that can be used. Indeed, other orders of taking differences are equally valid.
== Carrier phase tracking (surveying) ==
Utilizing the navigation message to measure pseudorange has been discussed. Another method that is used in GPS surveying applications is carrier phase tracking. The period of the carrier frequency times the speed of light gives the wavelength, which is about 0.19 meters for the L1 carrier. With a 1% of wavelength accuracy in detecting the leading edge, this component of pseudorange error might be as low as 2 millimeters. This compares to 3 meters for the C/A code and 0.3 meters for the P code.


The satellite carrier's total phase can be measured with ambiguity regarding the number of cycles. Let <math>\phi_{i,j,k}</math> denote the phase of the carrier of satellite <math>j</math> measured by receiver <math>i</math> at time <math>t_k</math>. Also, we define three functions: <math>\Delta^r, \Delta^s, \Delta^t</math>, which perform differences between receivers, satellites, and time points, respectively. Each function has a linear combination of variables with three subscripts as its argument. These three functions are defined below:
However, this 2 millimeter accuracy requires measuring the total phase, that is the total number of wavelengths plus the fractional wavelength. This requires specially equipped receivers. This method has many applications in the field of surveying.


:<math>\Delta^r(\phi_{i,j,k}) = \phi_{i+1,j,k} - \phi_{i,j,k},</math>
We now describe a method which could potentially be used to estimate the position of receiver 2 given the position of receiver 1 using triple differencing followed by numerical root finding, and a mathematical technique called [[least squares]]. A detailed discussion of the errors is omitted in order to avoid detracting from the description of the methodology. In this description differences are taken in the order of differencing between satellites, differencing between receivers, and differencing between epochs. This should not be construed to mean that this is the only order which can be used. Indeed, other orders of taking differences are equally valid.
:<math>\Delta^s(\phi_{i,j,k}) = \phi_{i,j+1,k} - \phi_{i,j,k},</math>
:<math>\Delta^t(\phi_{i,j,k}) = \phi_{i,j,k+1} - \phi_{i,j,k}.</math>


Also if <math>\phi_{i,j,k}</math> and <math>\phi_{l,m,n}</math> are valid arguments for the three functions, and <math>a</math> and <math>b</math> are constants, then <math>(a \phi_{i,j,k} + b \phi_{l,m,n})</math> is a valid argument with values defined as
The satellite carrier total phase can be measured with ambiguity as to the number of cycles. Let <math>\ \phi(r_i, s_j, t_k) </math> denote the phase of the carrier of satellite j measured by receiver i at time <math>\ \ t_k </math>. This notation has been chosen so as to make it clear what the subscripts i, j, and k mean. In view of the fact that the receiver, satellite, and time come in alphabetical order as arguments of <math>\ \phi </math> and to strike a balance between readability and conciseness, let <math>\ \phi_{i,j,k} = \phi(r_i, s_j, t_k) </math> so as to have a concise abbreviation. Also we define three functions, :<math>\ \Delta^r, \Delta^s, \Delta^t </math> which perform differences between receivers, satellites, and time points respectively. Each of these functions has a linear combination of variables with three subscripts as its argument. These three functions are defined below. If <math>\ \alpha_{i,j,k} </math> is a function of the three integer arguments, i, j, and k then it is a valid argument for the functions, :<math>\ \Delta^r, \Delta^s, \Delta^t </math>, with the values defined as


:<math>\ \Delta^r(\alpha_{i,j,k}) = \alpha_{i+1,j,k} - \alpha_{i,j,k} </math> ,
:<math>\Delta^r(a \phi_{i,j,k} + b \phi_{l,m,n}) = a \Delta^r(\phi_{i,j,k}) + b \Delta^r(\phi_{l,m,n}),</math>
:<math>\ \Delta^s(\alpha_{i,j,k}) = \alpha_{i,j+1,k} - \alpha_{i,j,k} </math> , and
:<math>\Delta^s(a \phi_{i,j,k} + b \phi_{l,m,n}) = a \Delta^s(\phi_{i,j,k}) + b \Delta^s(\phi_{l,m,n}),</math>
:<math>\ \Delta^t(\alpha_{i,j,k}) = \alpha_{i,j,k+1} - \alpha_{i,j,k} </math> .
:<math>\Delta^t(a \phi_{i,j,k} + b \phi_{l,m,n}) = a \Delta^t(\phi_{i,j,k}) + b \Delta^t(\phi_{l,m,n}).</math>


Receiver-clock errors can be approximately eliminated by differencing the phases measured from satellite 1 with that from satellite 2 at the same epoch.<ref>[http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/633.htm Between-satellite differencing] {{webarchive|url=https://web.archive.org/web/20110306051358/http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/633.htm |date=2011-03-06 }}</ref> This difference is designated as <math>\Delta^s(\phi_{1,1,1}) = \phi_{1,2,1} - \phi_{1,1,1}</math>.
Also if <math>\ \alpha_{i,j,k}\ and\ \beta_{l,m,n} </math> are valid arguments for the three functions and a and b are constants then
<math>\ ( a\ \alpha_{i,j,k} + b\ \beta_{l,m,n} ) </math> is a valid argument with values defined as


Double differencing can be performed by taking the differences of the satellite difference observed by receiver 1 with that observed by receiver 2.<ref>[http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/635.htm Double differencing] {{webarchive|url=https://web.archive.org/web/20110306051809/http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/635.htm |date=2011-03-06 }}</ref> The satellite-clock errors will be approximately eliminated by this between receiver differencing. This double difference is
:<math>\ \Delta^r(a\ \alpha_{i,j,k} + b\ \beta_{l,m,n}) = a \ \Delta^r(\alpha_{i,j,k}) + b \ \Delta^r(\beta_{l,m,n})</math> ,
:<math>\ \Delta^s(a\ \alpha_{i,j,k} + b\ \beta_{l,m,n} )= a \ \Delta^s(\alpha_{i,j,k}) + b \ \Delta^s(\beta_{l,m,n})</math> , and
:<math>\Delta^r(\Delta^s(\phi_{1,1,1})) = \Delta^r(\phi_{1,2,1} - \phi_{1,1,1})
:<math>\ \Delta^t(a\ \alpha_{i,j,k} + b\ \beta_{l,m,n} )= a \ \Delta^t(\alpha_{i,j,k}) + b \ \Delta^t(\beta_{l,m,n})</math> .
= \Delta^r(\phi_{1,2,1}) - \Delta^r(\phi_{1,1,1})
= (\phi_{2,2,1} - \phi_{1,2,1}) - (\phi_{2,1,1} - \phi_{1,1,1}).</math>


Triple differencing can be performed by taking the difference of double differencing performed at time <math>t_2</math> with that performed at time <math>t_1</math>.<ref>[http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/636.htm Triple differencing] {{webarchive|url=https://web.archive.org/web/20110306051924/http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/636.htm |date=2011-03-06 }}</ref> This will eliminate the ambiguity associated with the integral number of wavelengths in carrier phase, provided this ambiguity does not change with time. Thus the triple difference result has eliminated all or practically all clock bias errors and integer ambiguity. Also, errors associated with atmospheric delay and satellite ephemeris have been significantly reduced. This triple difference is
Receiver clock errors can be approximately eliminated by differencing the phases measured from satellite 1 with that from satellite 2 at the same epoch.<ref>[http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/633.htm Between-satellite differencing]</ref> This difference is designated as <math>\ \Delta^s(\phi_{1,1,1}) = \phi_{1,2,1} - \phi_{1,1,1}</math>
:<math>\Delta^t(\Delta^r(\Delta^s(\phi_{1,1,1}))).</math>


Triple-difference results can be used to estimate unknown variables. For example, if the position of receiver 1 is known, but the position of receiver 2 is unknown, it may be possible to estimate the position of receiver 2 using numerical root finding and [[least squares]]. Triple difference results for three independent time pairs quite possibly will be sufficient to solve for the three components of the position of receiver 2. This may require the use of a numerical procedure such as one of those found in the chapter on root finding and nonlinear sets of equations in Numerical Recipes.<ref name="NR">Press (1986), p. 959.</ref> To use such a numerical method, an initial approximation of the position of receiver 2 is required. This initial value could probably be provided by a position approximation based on the navigation message and the intersection of sphere surfaces. Although multidimensional numerical root finding can have problems, this disadvantage may be overcome with this good initial estimate. This procedure using three-time pairs and a fairly good initial value followed by iteration will result in one observed triple-difference result for the receiver 2 positions. Greater accuracy may be obtained by processing triple-difference results for additional sets of three independent time pairs. This will result in an overdetermined system with multiple solutions. To get estimates for an overdetermined system, least squares can be used. The least-squares procedure determines the position of receiver 2 that best fits the observed triple-difference results for receiver 2 positions under the criterion of minimizing the sum of the squares.
Double differencing can be performed by taking the differences of the between satellite difference observed by receiver 1 with that observed by receiver 2.<ref>[http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/635.htm Double differencing]</ref> The satellite clock errors will be approximately eliminated by this between receiver differencing. This double difference is:
:<math>\begin{align}
\Delta^r(\Delta^s(\phi_{1,1,1}))\,&=\,\Delta^r(\phi_{1,2,1} - \phi_{1,1,1})
&=\,\Delta^r(\phi_{1,2,1}) - \Delta^r(\phi_{1,1,1})
&=\,(\phi_{2,2,1} - \phi_{1,2,1}) - (\phi_{2,1,1} - \phi_{1,1,1})
\end{align}</math>


==GNSS augmentation==
Triple differencing can be performed by taking the difference of double differencing performed at time <math>\ \ t_2 </math> with that performed at time <math>\ \ t_1 </math>.<ref>[http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap6/636.htm Triple differencing]</ref> This will eliminate the ambiguity associated with the integral number of wavelengths in carrier phase provided this ambiguity does not change with time. Thus the triple difference result has eliminated all or practically all clock bias errors and the integer ambiguity. Also errors associated with atmospheric delay and satellite ephemeris have been significantly reduced. This triple difference is:
{{excerpt|GNSS augmentation}}
:<math>\ \Delta^t(\Delta^r(\Delta^s(\phi_{1,1,1}))) </math>


==Other enhancements==
Triple difference results can be used to estimate unknown variables. For example, if the position of receiver 1 is known but the position of receiver 2 unknown, it may be possible to estimate the position of receiver 2 using numerical root finding and [[least squares]]. Triple difference results for three independent time pairs quite possibly will be sufficient to solve for the three components of position of receiver 2. This may require the use of a numerical procedure such as one of those found in the chapter on root finding and nonlinear sets of equations in Numerical Recipes.<ref name="NR">Press (1986), p. 959.</ref> To use such a numerical method, an initial approximation of the position of receiver 2 is required. This initial value could probably be provided by a position approximation based on the navigation message and the intersection of sphere surfaces. Although multidimensional numerical root finding can have problems, this disadvantage may be overcome with this good initial estimate. This procedure using three time pairs and a fairly good initial value followed by iteration will result in one observed triple difference result for receiver 2 position. Greater accuracy may be obtained by processing triple difference results for additional sets of three independent time pairs. This will result in an over determined system with multiple solutions. To get estimates for an over determined system, least squares can be used. The least squares procedure determines the position of receiver 2 which best fits the observed triple difference results for receiver 2 positions under the criterion of minimizing the sum of the squares.
Other examples of GNSS enhancements include [[Inertial Navigation System]]s and [[Assisted GPS]].


== See also ==
== See also ==
Line 78: Line 73:


== Notes ==
== Notes ==
{{reflist}}
{{reflist|refs=
<ref name="ubiquitous">{{cite book |last=Mannings |first=Robin |date=2008 |title=Ubiquitous Positioning |url=https://books.google.com/books?id=kei6zO5iBsMC|publisher=Artech House |page=102 |isbn=978-1596931046}}</ref>
}}


== References ==
== References ==
* {{cite book |author1=Press, Flannery |author2=Tekolsky, Vetterling |last-author-amp=yes | year=1986 | title=Numerical Recipes, The Art of Scientific Computing | publisher=Cambridge University Press}}
* {{cite book |author1=Press, Flannery |author2=Tekolsky, Vetterling |name-list-style=amp | year=1986 | title=Numerical Recipes, The Art of Scientific Computing | publisher=Cambridge University Press}}
* <ref name="ubiquitous">{{cite book |last=Mannings |first=Robin |date=2008 |title=Ubiquitous Positioning |url=https://books.google.es/books/about/Ubiquitous_Positioning.html?id=kei6zO5iBsMC|publisher=Artech House |page=102 |isbn=1596931043}}</ref>


== External links ==
== External links ==
{{commons|Global Positioning System}}
{{commons|Global Positioning System}}
* [http://www.gps.gov/ GPS.gov]—General public education website created by the U.S. Government
* [http://www.gps.gov/ GPS.gov]—General public education website created by the U.S. Government
* U.S. Army Corps of Engineers manual: [http://www.usace.army.mil/publications/eng-manuals/em1110-1-1003/toc.htm NAVSTAR HTML] and [http://www.usace.army.mil/publications/eng-manuals/em1110-1-1003/entire.pdf PDF (22.6 MB, 328 pages)]
* U.S. Army Corps of Engineers manual: [https://web.archive.org/web/20080822132227/http://www.usace.army.mil/publications/eng-manuals/em1110-1-1003/toc.htm NAVSTAR HTML] and [https://web.archive.org/web/20080625111519/http://www.usace.army.mil/publications/eng-manuals/em1110-1-1003/entire.pdf PDF (22.6 MB, 328 pages)]
* [http://pnt.gov/public/docs/2008/spsps2008.pdf GPS SPS Performance Standard]—The official Standard Positioning Service specification (2008 version).
* [https://web.archive.org/web/20091229023237/http://pnt.gov/public/docs/2008/spsps2008.pdf GPS SPS Performance Standard]—The official Standard Positioning Service specification (2008 version).
* [http://gps.afspc.af.mil/gpsoc/documents/GPS_Signal_Spec.pdf GPS SPS Performance Standard]—The official Standard Positioning Service specification (2001 version).
* [http://gps.afspc.af.mil/gpsoc/documents/GPS_Signal_Spec.pdf GPS SPS Performance Standard]{{Dead link|date=February 2022 |bot=InternetArchiveBot |fix-attempted=yes }}—The official Standard Positioning Service specification (2001 version).
* [http://gps.afspc.af.mil/gpsoc/documents/PPS_PS_Signed_Final_23_Feb_07.pdf GPS PPS Performance Standard]—The official Precise Positioning Service specification.
* [http://gps.afspc.af.mil/gpsoc/documents/PPS_PS_Signed_Final_23_Feb_07.pdf GPS PPS Performance Standard] {{Webarchive|url=https://web.archive.org/web/20091224204720/http://gps.afspc.af.mil/gpsoc/documents/PPS_PS_Signed_Final_23_Feb_07.pdf |date=2009-12-24 }}—The official Precise Positioning Service specification.


{{Satellite navigation systems}}
{{Satellite navigation systems}}

Latest revision as of 13:32, 1 October 2024

GNSS enhancement refers to techniques used to improve the accuracy of positioning information provided by the Global Positioning System or other global navigation satellite systems in general, a network of satellites used for navigation. Enhancement methods of improving accuracy rely on external information being integrated into the calculation process. There are many such systems in place and they are generally named or described based on how the GPS sensor receives the information. Some systems transmit additional information about sources of error (such as clock drift, ephemeris, or ionospheric delay), others provide direct measurements of how much the signal was off in the past, while a third group provides additional navigational or vehicle information to be integrated into the calculation process.

Background

[edit]

The Global Positioning System (GPS) is the American satellite-based system for positioning and navigation. Receivers on or near the Earth's surface can determine their locations based on signals received from any four or more of the satellites in the network.

All GPS satellites broadcast on the same two frequencies, known as L1 (1575.42 MHz) and L2 (1227.60 MHz). The network uses code-division multiple access (CDMA) to allow separate messages from the individual satellites to be distinguished. Two distinct CDMA encodings are used: the coarse/acquisition (C/A) code, which is accessible by the general public, and the precise (P) code, which is encrypted, so that only the U.S. military can access it. The messages sent from each satellite contain information ranging from the satellite's health, the satellite's orbital path, the clock state of the satellite, and the configuration of the entire satellite network.

Precise monitoring

[edit]

The accuracy of a calculation can also be improved through precise monitoring and measuring of the existing GPS signals in additional or alternate ways.

After Selective Availability was turned off by the U.S. government, the largest error in GPS was usually the unpredictable delay through the ionosphere. The spacecraft broadcast ionospheric model parameters, but they are necessarily imperfect. This is one reason the GPS spacecraft transmit on at least two frequencies, L1 and L2. Ionospheric delay is a well-defined function of frequency and the total electron content (TEC) along the path, so measuring the arrival-time difference between the frequencies determines TEC and thus the precise ionospheric delay at each frequency.

Receivers with decryption keys can decode the P(Y)-code transmitted on both L1 and L2. However, these keys are reserved for the military and authorized agencies and are unavailable to the public. Without keys, it is still possible to use a codeless technique to compare the P(Y) codes on L1 and L2 to gain much of the same error information. However, this technique is slow and is currently limited to specialized surveying equipment. In the future, additional civilian codes are expected to be transmitted on the L2 and L5 frequencies (see GPS modernization). Then all users can perform dual-frequency measurements and directly compute ionospheric-delay errors.

A second form of precise monitoring is called carrier-phase enhancement (CPGPS).[1] The error, which this corrects, arises because the pulse transition of the PRN is not instantaneous, and thus the correlation (satellite–receiver sequence matching) operation is imperfect. The CPGPS approach utilizes the L1 carrier wave, which has a period of

which is about one-thousandth of the C/A Gold code bit period of

to act as an additional clock signal and resolve the uncertainty. The phase-difference error in the normal GPS amounts to between 2 and 3 meters (6 to 10 ft) of ambiguity. CPGPS working to within 1% of perfect transition reduces this error to 3 centimeters (1 inch) of ambiguity. By eliminating this source of error, CPGPS coupled with DGPS normally realizes between 20 and 30 centimeters (8 to 12 inches) of absolute accuracy.

Real-time kinematic positioning

[edit]

Real-time kinematic positioning (RTK) is another approach for a precise GPS-based positioning system. In this approach, the determination of the range signal can be resolved to a precision of less than 10 centimeters (4 in). This is done by resolving the number of cycles in which the signal is transmitted and received by the receiver. This can be accomplished by using a combination of differential GPS (DGPS) correction data, transmitting GPS signal phase information, and ambiguity resolution techniques via statistical tests, possibly with processing in real time.

Carrier-phase tracking (surveying)

[edit]

Utilizing the navigation message to measure pseudo-range has been discussed. Another method that is used in GPS surveying applications is carrier-phase tracking. The period of the carrier frequency times the speed of light gives the wavelength, which is about 0.19 meters for the L1 carrier. With a 1% of wavelength accuracy in detecting the leading edge, this component of pseudo-range error might be as low as 2 millimeters. This compares to 3 meters for the C/A code and 0.3 meters for the P code.

However, this 2-millimeter accuracy requires measuring the total phase, the total number of wavelengths, plus the fractional wavelength. This requires specially equipped receivers. This method has many applications in the field of surveying.

We now describe a method that could potentially be used to estimate the position of receiver 2 given the position of receiver 1 using triple differencing followed by numerical root finding and a mathematical technique called least squares. A detailed discussion of the errors is omitted in order to avoid detracting from the description of the methodology. In this description, differences are taken in the order of differencing between satellites, differencing between receivers, and differencing between epochs. This should not be construed to mean that this is the only order that can be used. Indeed, other orders of taking differences are equally valid.

The satellite carrier's total phase can be measured with ambiguity regarding the number of cycles. Let denote the phase of the carrier of satellite measured by receiver at time . Also, we define three functions: , which perform differences between receivers, satellites, and time points, respectively. Each function has a linear combination of variables with three subscripts as its argument. These three functions are defined below:

Also if and are valid arguments for the three functions, and and are constants, then is a valid argument with values defined as

Receiver-clock errors can be approximately eliminated by differencing the phases measured from satellite 1 with that from satellite 2 at the same epoch.[2] This difference is designated as .

Double differencing can be performed by taking the differences of the satellite difference observed by receiver 1 with that observed by receiver 2.[3] The satellite-clock errors will be approximately eliminated by this between receiver differencing. This double difference is

Triple differencing can be performed by taking the difference of double differencing performed at time with that performed at time .[4] This will eliminate the ambiguity associated with the integral number of wavelengths in carrier phase, provided this ambiguity does not change with time. Thus the triple difference result has eliminated all or practically all clock bias errors and integer ambiguity. Also, errors associated with atmospheric delay and satellite ephemeris have been significantly reduced. This triple difference is

Triple-difference results can be used to estimate unknown variables. For example, if the position of receiver 1 is known, but the position of receiver 2 is unknown, it may be possible to estimate the position of receiver 2 using numerical root finding and least squares. Triple difference results for three independent time pairs quite possibly will be sufficient to solve for the three components of the position of receiver 2. This may require the use of a numerical procedure such as one of those found in the chapter on root finding and nonlinear sets of equations in Numerical Recipes.[5] To use such a numerical method, an initial approximation of the position of receiver 2 is required. This initial value could probably be provided by a position approximation based on the navigation message and the intersection of sphere surfaces. Although multidimensional numerical root finding can have problems, this disadvantage may be overcome with this good initial estimate. This procedure using three-time pairs and a fairly good initial value followed by iteration will result in one observed triple-difference result for the receiver 2 positions. Greater accuracy may be obtained by processing triple-difference results for additional sets of three independent time pairs. This will result in an overdetermined system with multiple solutions. To get estimates for an overdetermined system, least squares can be used. The least-squares procedure determines the position of receiver 2 that best fits the observed triple-difference results for receiver 2 positions under the criterion of minimizing the sum of the squares.

GNSS augmentation

[edit]
Augmentation of a global navigation satellite system (GNSS) is a method of improving the navigation system's attributes, such as precision, reliability, and availability, through the integration of external information into the calculation process. There are many such systems in place, and they are generally named or described based on how the GNSS sensor receives the external information. Some systems transmit additional information about sources of error (such as clock drift, ephemeris, or ionospheric delay), others provide direct measurements of how much the signal was off in the past, while a third group provides additional vehicle information to be integrated in the calculation process.

Other enhancements

[edit]

Other examples of GNSS enhancements include Inertial Navigation Systems and Assisted GPS.

See also

[edit]

Notes

[edit]
  1. ^ Mannings, Robin (2008). Ubiquitous Positioning. Artech House. p. 102. ISBN 978-1596931046.
  2. ^ Between-satellite differencing Archived 2011-03-06 at the Wayback Machine
  3. ^ Double differencing Archived 2011-03-06 at the Wayback Machine
  4. ^ Triple differencing Archived 2011-03-06 at the Wayback Machine
  5. ^ Press (1986), p. 959.

References

[edit]
  • Press, Flannery & Tekolsky, Vetterling (1986). Numerical Recipes, The Art of Scientific Computing. Cambridge University Press.
[edit]