If you wish to contribute or participate in the discussions about articles you are invited to contact the Editor

Precise modelling terms for PPP: Difference between revisions

From Navipedia
Jump to navigation Jump to search
No edit summary
Line 9: Line 9:
The [[Precise Point Positioning|Precise Point Positioning (PPP)]] technique allows achieving centimetre level accuracy for static positioning and decimetre level for kinematic, typically, after the best part of one hour. This high accuracy requires an accurate measurements modelling, where all model terms must be taken into account (up to the centimetre level or better).<ref name="Kouba">Kouba, Jan (et al.), Precise Point Positioning, Chapter 25, Handbook of Global Navigation Satellite Systems, 2017</ref><ref name="Hernandez">M. Hernandez-Pajares, J.M. Juan, J. Sanz, R. Orus, A. Garcia-Rigo, J. Feltens, A. Komjathy, S.C. Schaer, A. Krankowski: The IGS VTEC maps: A reliable source of ionospheric information since 1998, J. Geodesy 83(3/4), 263–275 (2009)</ref><ref name="King">M.A. King, Z. Altamimi, J. Boehm, M. Bos, R. Dach, P. Elosegui, F. Fund, M. Hernández-Pajares, D. Lavallee, P.J. Mendes Cerveira, N. Penna, R.E.M. Riva, P. Steigenberger, T. van Dam, L. Vittuari, S. Williams, P. Willis: Improved constraints on models of glacial isostatic adjustment: A review of the contribution of ground-based geodetic observations, Surveys Geophys. 31(5), 465–507 (2010)</ref>
The [[Precise Point Positioning|Precise Point Positioning (PPP)]] technique allows achieving centimetre level accuracy for static positioning and decimetre level for kinematic, typically, after the best part of one hour. This high accuracy requires an accurate measurements modelling, where all model terms must be taken into account (up to the centimetre level or better).<ref name="Kouba">Kouba, Jan (et al.), Precise Point Positioning, Chapter 25, Handbook of Global Navigation Satellite Systems, 2017</ref><ref name="Hernandez">M. Hernandez-Pajares, J.M. Juan, J. Sanz, R. Orus, A. Garcia-Rigo, J. Feltens, A. Komjathy, S.C. Schaer, A. Krankowski: The IGS VTEC maps: A reliable source of ionospheric information since 1998, J. Geodesy 83(3/4), 263–275 (2009)</ref><ref name="King">M.A. King, Z. Altamimi, J. Boehm, M. Bos, R. Dach, P. Elosegui, F. Fund, M. Hernández-Pajares, D. Lavallee, P.J. Mendes Cerveira, N. Penna, R.E.M. Riva, P. Steigenberger, T. van Dam, L. Vittuari, S. Williams, P. Willis: Improved constraints on models of glacial isostatic adjustment: A review of the contribution of ground-based geodetic observations, Surveys Geophys. 31(5), 465–507 (2010)</ref>


This modelling involves the following terms, among those considered in the [[Code Based Positioning (SPS)|SPS]]:
This modelling involves the following terms, among those considered in [[Code Based Positioning (SPS)|Code Based Positioning]]:




*''' Precise satellite orbits and clocks:'''
*''' Precise satellite orbits and clocks:'''
: ''The precise orbits and clocks files'' (see [[Precise GNSS Satellite Coordinates Computation]]) must be used instead of the broadcast ones used in conventional code-based pseudorange positioning.
: ''The precise orbits and clocks files'' (see [[Precise GNSS Satellite Coordinates Computation|Precise GNSS Satellite Coordinates Computation]]) must be used instead of the broadcast ones used in conventional [[Code Based Positioning (SPS)|Code Based Positioning]].


:The satellite coordinates between epochs can be computed by polynomial interpolation.  
:The satellite coordinates between epochs can be computed by polynomial interpolation.  
Line 40: Line 40:




: ''Relativistic effects''(see <ref name=Ashby">Ashby, N., 2003. Relativity in the Global Positioning System</ref> for a very good review of relativistic effects on GPS): The gravitational path range correction (3) (see [[Relativistic Path Range Effect]]), can be added among the satellite clock correction due to the orbit eccentricity (4) considered in the [[Code Based Positioning (SPS)|SPS]] (see [[Relativistic Clock Correction]]).
: ''Relativistic effects''(see <ref name=Ashby">Ashby, N., 2003. Relativity in the Global Positioning System</ref> for a very good review of relativistic effects on GPS): The gravitational path range correction (3) (see [[Relativistic Path Range Effect]]), can be added among the satellite clock correction due to the orbit eccentricity (4) considered in [[Code Based Positioning (SPS)|Code Based Positioning]] (see [[Relativistic Clock Correction]]).


::<math>
::<math>
Line 83: Line 83:


* '''Antenna biases and orientation:'''
* '''Antenna biases and orientation:'''
: The ''satellite and receiver antenna phase centres'' can be found in the IGS ANTEX files, after the GPS week <math>1400</math> (see [[Antenna Phase Centre]]).
: The ''satellite and receiver antenna phase centres'' can be found in the IGS ANTEX files, after the GPS week 1400 (see [[Antenna Phase Centre]]).





Revision as of 13:30, 15 June 2020


FundamentalsFundamentals
Title Precise modelling terms for PPP
Author(s) J. Sanz Subirana, J.M. Juan Zornoza and M. Hernández-Pajares, Technical University of Catalonia, Spain.
Level Advanced
Year of Publication 2011
Logo gAGE.png

The Precise Point Positioning (PPP) technique allows achieving centimetre level accuracy for static positioning and decimetre level for kinematic, typically, after the best part of one hour. This high accuracy requires an accurate measurements modelling, where all model terms must be taken into account (up to the centimetre level or better).[1][2][3]

This modelling involves the following terms, among those considered in Code Based Positioning:


  • Precise satellite orbits and clocks:
The precise orbits and clocks files (see Precise GNSS Satellite Coordinates Computation) must be used instead of the broadcast ones used in conventional Code Based Positioning.
The satellite coordinates between epochs can be computed by polynomial interpolation.
As an example of polynomial interpolation, the Lagrange method is presented as follows. Given a table of values [math]\displaystyle{ (x_i,y_i), i=1,…,n, }[/math] the interpolated value [math]\displaystyle{ y≃P_n(x) }[/math] at a given x can be computed as (see Precise GNSS Satellite Coordinates Computation):
[math]\displaystyle{ \begin{array}{lll} P_n(x)&=& \sum_{i=1}^n{y_i \frac{\prod_{j\neq i}{(x-x_j)}}{\prod_{j\neq i}{(x_i-x_j)}}}=\\ &=& y_1 \frac{x-x_2}{(x_1-x_2)} \cdots \frac{x-x_n}{(x_1-x_n)}+\cdots\\ & &+y_i \frac{x-x_1}{(x_i-x_1)} \cdots \frac{x-x_{i-1}}{(x_i-x_{i-1})}\frac{x-x_{i+1}}{(x_i-x_{i+1})}\cdots \frac{x-x_n}{(x_i-x_n)}+\cdots\\ & & +y_n \frac{x-x_1}{(x_n-x_1)} \cdots\frac{x-x_{n-1}}{(x_n-x_{n-1})} \end{array} \qquad \mbox{(1)} }[/math]


Notice that, the orbits are referred to the satellite mass centre, and thence, the Satellite Antenna Phase Centre offset must be applied to compute the [math]\displaystyle{ {\boldsymbol \Delta}_{APC} }[/math] vector offset (see Satellite Antenna Phase Centre).
[math]\displaystyle{ {\mathbf r}_{sat_{_{APC}}}={\mathbf r}_{sat_{_{MC}}}+{\mathbf R}\cdot {\boldsymbol \Delta}_{_{APC}} \qquad \mbox{(2)} }[/math]


The satellite clocks should not be interpolated and thence, only epochs having clocks available must be used (Precise GNSS Satellite Coordinates Computation).


Relativistic effects(see [4] for a very good review of relativistic effects on GPS): The gravitational path range correction (3) (see Relativistic Path Range Effect), can be added among the satellite clock correction due to the orbit eccentricity (4) considered in Code Based Positioning (see Relativistic Clock Correction).
[math]\displaystyle{ \Delta \rho_{_{rel}}=\frac{2\,\mu}{c^2}ln\frac{r^{sat}+r_{rcv}+r_{rcv}^{sat}}{r^{sat}+r_{rcv}-r_{rcv}^{sat}} \qquad \mbox{(3)} }[/math]
[math]\displaystyle{ \Delta_{rel}=- 2\, \frac{\mathbf{r} \cdot \mathbf{v}}{c^2} \qquad \mbox{(4)} }[/math]


  • Atmospheric effects:
The ionospheric refraction and DCBs are removed using the ionosphere free combination of measurements (see details in Ionosphere-free Combination for Dual Frequency Receivers).


The tropospheric refraction can be modelled by (5) (see Tropospheric Delay )
[math]\displaystyle{ T(E)=T_{z,dry}\cdot M_{dry}(E)+T_{z,wet}\cdot M_{wet}(E) \qquad \mbox{(5)} }[/math]


where the mapping of Niell is used. The dry and wet tropospheric delays are given by
[math]\displaystyle{ \begin{array}{l} T_{z,dry}= a\; e^{-b\; H}\\ T_{z,wet}= T_{z_0,wet}+\Delta T_{z,wet} \qquad \mbox{(6)} \end{array} }[/math]


The deviation of the zenith tropospheric delay [math]\displaystyle{ \Delta T_{z,wet} }[/math] regarding to the nominal [math]\displaystyle{ T_{z_0,wet} }[/math] must be estimated in the Kalman filter, together with the coordinates, clock and carrier phase biases.


  • Antenna biases and orientation:
The satellite and receiver antenna phase centres can be found in the IGS ANTEX files, after the GPS week 1400 (see Antenna Phase Centre).


The carrier phase wind-up effect due to the satellite motion is given by (7) (see Carrier Phase Wind-up Effect)
[math]\displaystyle{ \Delta \phi=2N\pi+sign(\zeta)\cdot \arccos \left( \frac{\vec{D}^{\prime}\cdot\vec{D}}{\|\vec{D}^{\prime}\|\cdot\|\vec{D}\|}\right) \qquad \mbox{(7)} }[/math]


The satellites under eclipse conditions should be removed from the computation due to the largest orbit error. The eclipse condition is given by (8) (see Satellite Eclipses)
[math]\displaystyle{ \cos \phi = \frac{{\mathbf r}_{sat}\cdot {\mathbf r}_{sun}}{|{\mathbf r}_{sat}\cdot {\mathbf r}_{sun}|} \lt 0 \qquad \mbox{and} \qquad r_{sat} \sqrt{1-\cos^2 \phi}\, \lt a_e \qquad \mbox{(8)} }[/math]


  • Earth deformation effects:
Solid tides can be modelled by equations (9) and (10) (see Solid Tides )
[math]\displaystyle{ \Delta {\mathbf r}= \sum_{j=2}^{3}{\frac{G\,M_j\,R_e^4}{G\,M_\oplus\,R_j^3}} \left \{h_2 \,\hat{\mathbf r} \left ( \frac{3}{2} (\hat{\mathbf R}_j \cdot \hat{\mathbf r})^2 -\frac{1}{2}\right) + 3\,l_2\,(\hat{\mathbf R}_j \cdot \hat{\mathbf r}) \left [\hat{\mathbf R}_j-(\hat{\mathbf R}_j \cdot \hat{\mathbf r})\,\hat{\mathbf r} \right ]\right \} \qquad \mbox{(9)} }[/math]
[math]\displaystyle{ \Delta {\mathbf r}= \sum_{j=2}^{3}{\frac{G\,M_j\,R_e^5}{G\,M_\oplus\,R_j^4}} \left \{h_3 \,\hat{\mathbf r} \left ( \frac{5}{2} (\hat{\mathbf R}_j \cdot \hat{\mathbf r})^3 -\frac{3}{2} (\hat{\mathbf R}_j \cdot \hat{\mathbf r})\right) + l_3\,\left (\frac{15}{2} (\hat{\mathbf R}_j \cdot \hat{\mathbf r})^2 - \frac{3}{2} \right ) \left [\hat{\mathbf R}_j-(\hat{\mathbf R}_j \cdot \hat{\mathbf r})\,\hat{\mathbf r} \right ]\right \} \qquad \mbox{(10)} }[/math]


Ocean loading and Pole Tides are second order effects and can be neglected for PPP accuracies at the centimetre level (see comments in Ocean loading and Pole Tide).


Related Articles

References

  1. ^ Kouba, Jan (et al.), Precise Point Positioning, Chapter 25, Handbook of Global Navigation Satellite Systems, 2017
  2. ^ M. Hernandez-Pajares, J.M. Juan, J. Sanz, R. Orus, A. Garcia-Rigo, J. Feltens, A. Komjathy, S.C. Schaer, A. Krankowski: The IGS VTEC maps: A reliable source of ionospheric information since 1998, J. Geodesy 83(3/4), 263–275 (2009)
  3. ^ M.A. King, Z. Altamimi, J. Boehm, M. Bos, R. Dach, P. Elosegui, F. Fund, M. Hernández-Pajares, D. Lavallee, P.J. Mendes Cerveira, N. Penna, R.E.M. Riva, P. Steigenberger, T. van Dam, L. Vittuari, S. Williams, P. Willis: Improved constraints on models of glacial isostatic adjustment: A review of the contribution of ground-based geodetic observations, Surveys Geophys. 31(5), 465–507 (2010)
  4. ^ Ashby, N., 2003. Relativity in the Global Positioning System