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

Positioning Error

From Navipedia
Revision as of 13:09, 23 January 2012 by Jaume.Sanz (talk | contribs)
Jump to navigation Jump to search


FundamentalsFundamentals
Title Positioning Error
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 formal, predicted and measured accuracy parameters are discussed in this article. The first one is a measure of the uncertainty of the estimates, according to the statistical characterisation of the errors and the linear model used for the position estimate.

The predicted accuracy, provides the expected position accuracy based on a simple statistical description of the measurement errors. Its computation does not require the measurements, but only its standard deviation [math]\displaystyle{ \sigma }[/math], and the approximated satellites and user location coordinates. Thence, it can be computed at any point and at any time with the almanac, without needing the measurements.

The measured accuracy is the measure of the true error and must be assessed with the actual measurements.

Formal Accuracy

Denoting as [math]\displaystyle{ P_{xx} }[/math], [math]\displaystyle{ P_{yy} }[/math], [math]\displaystyle{ P_{zz} }[/math],[math]\displaystyle{ P_{tt} }[/math] the diagonal elements of matrix [math]\displaystyle{ {\mathbf P}_{_{\mathbf \Delta X}} }[/math] of (see Best Linear Unbiased Minimum-Variance Estimator (BLUE))

[math]\displaystyle{ {\mathbf P}=({\mathbf G}^T\,{\mathbf R}^{-1}\,{\mathbf G})^{-1} \qquad \mbox{(1)} }[/math]


we call formal errors of the estimated components of vector [math]\displaystyle{ {\mathbf \Delta X}=(dx,dy,dz,dt) }[/math] to the standard deviations:

[math]\displaystyle{ \sigma_X=\sqrt{P_{xx}} \qquad \sigma_Y=\sqrt{P_{yy}} \qquad \sigma_Z=\sqrt{P_{zz}} \qquad \sigma_T=\sqrt{P_{tt}} \qquad \mbox{(2)} }[/math]


The previous expressions provides a characterisation of the quality of the coordinates and clock estimates (linked to the measurement error model assumed). Thence, they are not the actual error, but only a measure of the uncertainty of the error estimates.


Previous expressions give the errors in the ECEF XYZ coordinates. Nevertheless, it is usually more meaningful to a user to think in terms of horizontal and vertical position error, or East, North and Up (ENU) coordinates.


Let [math]\displaystyle{ {\mathbf R} }[/math] the transformation matrix (3) of ENU coordinates to XYZ, (i.e., which columns are the unit orthogonal vectors [math]\displaystyle{ \{ \hat{\mathbf e}, \hat{\mathbf n}, \hat{\mathbf u}\} }[/math] as expressed in the XYZ coordinates system [footnotes 1] at a point of longitude [math]\displaystyle{ \lambda }[/math] and latitude [math]\displaystyle{ \varphi }[/math]) (see Transformations between ECEF and ENU coordinates):

[math]\displaystyle{ {\mathbf R}= \left [ \begin{array}{ccc} - \sin \lambda & - \sin \varphi \cos \lambda & \cos \varphi \cos \lambda \\ \cos \lambda & -\sin \varphi \sin \lambda & \cos \varphi \sin \lambda \\ 0 & \cos \varphi & \sin \varphi \\ \end{array} \right ] \qquad \mbox{(3)} }[/math]


thence, [math]\displaystyle{ {\mathbf P}_{enu}={\mathbf R}^T\,{\mathbf P}_{xyz}\,{\mathbf R} }[/math] is given, where [math]\displaystyle{ {\mathbf P}_{xyz} }[/math] is the sub-matrix of [math]\displaystyle{ {\mathbf P}_{_{\mathbf \Delta X}} }[/math] containing solely the geometric components [footnotes 2]. From this covariance matrix, they can be defined:

[math]\displaystyle{ \sigma_{East}=\sqrt{P_{ee}} \qquad \sigma_{North}=\sqrt{P_{nn}} \qquad \sigma_{Vertical}=\sqrt{P_{uu}} \qquad \mbox{(4)} }[/math]


For the horizontal error, there are usually also defined:

[math]\displaystyle{ \sigma_{Horizontal}=\sqrt{P_{ee}+ P_{nn}} \qquad \mbox{(5)} }[/math]


It is usual to express the Horizontal error in two main directions where the covariance matrix [math]\displaystyle{ {\mathbf P}_{enu} }[/math] is diagonal. They define two orthogonal axis, associated to the minor and major axis of the error ellipse. The major axis is given by:

[math]\displaystyle{ \sigma_{H_{major}}=\sqrt{\frac{P_{ee}+ P_{nn}}{2}+\sqrt{\left ( \frac{P_{ee}-P_{nn}}{2} \right )^2+P_{en}^2}} \qquad \mbox{(6)} }[/math]


Predicted Accuracy: Dilution of Precision

Assuming the simple statistical model:

[math]\displaystyle{ \begin{array}{l} {\mathbf m}=E[{\boldsymbol \varepsilon}]=0\;\;\;; {\mathbf R}=E[{\boldsymbol \varepsilon} \, {\boldsymbol \varepsilon}^T]=\sigma^2\,{\mathbf I} \end{array} \qquad \mbox{(7)} }[/math]


the covariance matrix of [math]\displaystyle{ {\mathbf{\Delta X}} }[/math] estimate is given by (see Parameters adjustment)

[math]\displaystyle{ {\mathbf P}_{_{\mathbf{\Delta X}}}=\sigma^2\,({\mathbf G}^T\,{\mathbf G})^{-1}\;;\;\;\; \mbox{MSE}_{_{\mathbf \Delta X}}=n\,\sigma^2 \qquad \mbox{(8)} }[/math]


as:
[math]\displaystyle{ {\mathbf P}_{_{\mathbf{\Delta X}}}=\sigma^2\,({\mathbf G}^T\,{\mathbf G})^{-1} \qquad \mbox{(9)} }[/math]


This expression depends upon two factors: 1) the measurements (prefit-residuals) variance [math]\displaystyle{ \sigma^2 }[/math], and 2) the geometry matrix ([math]\displaystyle{ {\mathbf G} }[/math]), which is linked only to the receiver-satellite rays geometry, (see Code Based Positioning (SPP)).


The Root Mean Square (RMS) is given by:

[math]\displaystyle{ \mbox{RMS}_{_{\mathbf \Delta X}}= \sqrt{trace ({\mathbf P}_{_{\mathbf \Delta X}})}=\sigma \sqrt{trace \left [({\mathbf G}^T\,{\mathbf G})^{-1}\right ]} \qquad \mbox{(10)} }[/math]


what means that, the trace of matrix [math]\displaystyle{ ({\mathbf G}^T\,{\mathbf G})^{-1} }[/math] is a scale factor on [math]\displaystyle{ \sigma }[/math] for the RMS[math]\displaystyle{ _{_{\mathbf \Delta X}} }[/math].


Notice that, as the matrix [math]\displaystyle{ {\mathbf G} }[/math] does not depend on the measurements, but only on the geometry, it can be computed from the almanac (because not accurate satellite positions are needed), i.e., does not require receiver measurements.


On the basis of this simple approach, the following Dilution of Precision (DOP) parameters are defined:

[math]\displaystyle{ {\mathbf Q}_{_{\mathbf \Delta X}} \equiv({\mathbf G}^T\,{\mathbf G})^{-1}=\left[ \begin{array}{cccc} q_{xx}&q_{xy}&q_{xz}&q_{xt} \\ q_{xy}&q_{yy}&q_{yz}&q_{yt} \\ q_{xz}&q_{yz}&q_{zz}&q_{zt} \\ q_{xt}&q_{yt}&q_{zt}&q_{tt} \\ \end{array} \right] \qquad \mbox{(11)} }[/math]


  • Geometric Dilution of Precision:
[math]\displaystyle{ GDOP=\sqrt{q_{xx}+q_{yy}+q_{zz}+q_{tt}} }[/math]
  • Position Dilution of Precision:
[math]\displaystyle{ PDOP=\sqrt{q_{xx}+q_{yy}+q_{zz}} }[/math]
  • Time Dilution of Precision:
[math]\displaystyle{ TDOP=\sqrt{q_{tt}} }[/math]


Using (3) the sub-matrix matrix [math]\displaystyle{ {\mathbf Q}_{xyz} }[/math] of [math]\displaystyle{ {\mathbf Q}_{_{\mathbf \Delta X}} }[/math] can be transformed to ENU coordinates as [math]\displaystyle{ {\mathbf Q}_{enu}={\mathbf R}^T\,{\mathbf Q}_{xyz}\,{\mathbf R} }[/math], in order to define:


  • Horizontal Dilution of Precision:
[math]\displaystyle{ HDOP=\sqrt{q_{ee}+q_{nn}} }[/math]


  • Vertical Dilution of Precision:
[math]\displaystyle{ VDOP=\sqrt{q_{uu}} }[/math]


Thence, estimations of the expected accuracy are given by:

  • GDOP, [math]\displaystyle{ \sigma_{0} }[/math] ... geometric precision in position and time
  • PDOP, [math]\displaystyle{ \sigma_{0} }[/math] ... precision in position
  • TDOP , [math]\displaystyle{ \sigma_{0} }[/math] ... precision in time
  • HDOP , [math]\displaystyle{ \sigma_{0} }[/math] ... precision in horizontal positioning
  • VDOP , [math]\displaystyle{ \sigma_{0} }[/math] ... precision in vertical positioning


where, basically, DOP represents an approximate ratio factor between the precision in the measurements ([math]\displaystyle{ \sigma_0 }[/math]) and in positioning. This ratio is computed only from the satellites-receiver geometry.


Measured Accuracy

Let [math]\displaystyle{ \Delta E_i }[/math], [math]\displaystyle{ \Delta N_i }[/math] and [math]\displaystyle{ \Delta U_i }[/math] the errors in the East, North and Vertical components of the [math]\displaystyle{ i }[/math]-th position estimate sample. The RMS vertical, horizontal (2-D) and 3-D errors are defined as [footnotes 3]:

[math]\displaystyle{ \mbox{rms vertical error}=\sqrt{\frac{1}{n}\sum_{i=1}^n{\Delta U_i^2}} \qquad \mbox{(12)} }[/math]
[math]\displaystyle{ \mbox{2-D rms horizontal error}=\sqrt{\frac{1}{n}\sum_{i=1}^n{\left ( \Delta E_i^2+\Delta N_i^2 \right )^2}} \qquad \mbox{(13)} }[/math]
[math]\displaystyle{ \mbox{3-D rms error}=\sqrt{\frac{1}{n}\sum_{i=1}^n{\left ( \Delta E_i^2+\Delta N_i^2 +\Delta U_i^2\right )^2}} \qquad \mbox{(14)} }[/math]


Other measures of the quality of the position estimates are:

  • 50th or 95th percentiles of horizontal, vertical and 3-D errors.
  • CEP: Circular Error Probable, as the 50th percentile of horizontal error.
  • SEP: Spherical Error Probable, as the 50th percentile of 3-D error.
  • 2drms: calculated as twice the 2-D rms horizontal error given by (13).


Assuming that the position estimates follows a multivariate normal distribution centred at the true position and the errors in north, east and up are uncorrelated, with [math]\displaystyle{ \sigma_E=\sigma_N }[/math] and [math]\displaystyle{ \sigma_U=2(\sigma_e^2+\sigma_N^2)^{1/2} }[/math], thence, the following relations can be considered ( [Misra and Enge, 2001] [1]):


rms vertical error [math]\displaystyle{ \qquad \simeq 0.5 \times }[/math] vertical error (95%)
[math]\displaystyle{ \simeq }[/math] 2drms
0.9 [math]\displaystyle{ \times }[/math] 3-D rms error
2-D rms error [math]\displaystyle{ \qquad =0.5 \times }[/math] 2drms
[math]\displaystyle{ \simeq 0.6 \times }[/math] horizontal error (95%)
1.2 [math]\displaystyle{ \times }[/math] CEP
3-D rms error [math]\displaystyle{ \qquad =2.2 \times }[/math] 2-D rms error
1.2 [math]\displaystyle{ \times }[/math] horizontal error (95%)
1.3 [math]\displaystyle{ \times }[/math] SEP


A wide discussion on such accuracy equivalences can be found in [Diggelen, 2007] [2].


Notes

  1. ^ Note [math]\displaystyle{ (\lambda,\varphi) }[/math] are the ellipsoidal coordinates and, thence, the vector [math]\displaystyle{ {\mathbf u} }[/math] is orthogonal to the tangent plane to the ellipsoid, which is defined by [math]\displaystyle{ (\hat{\mathbf e}, \hat{\mathbf n}) }[/math]. If [math]\displaystyle{ (\lambda,\varphi) }[/math] are taken as the spherical latitude and longitude, thence, the vector [math]\displaystyle{ {\mathbf u} }[/math] is in the radial direction and [math]\displaystyle{ (\hat{\mathbf e}, \hat{\mathbf n}) }[/math] defines the tangent plane to the sphere.
  2. ^ The matrix [math]\displaystyle{ {\mathbf P}_{enu} }[/math] can be obtained directly in the ENU coordinates, using the geometry matrix [math]\displaystyle{ {\mathbf G} }[/math] computed in the ENU coordinates. In this local system the rows of matrix [math]\displaystyle{ {\mathbf G} }[/math] (see Code Based Positioning (SPP)) are [math]\displaystyle{ [\cos el^i \sin az^i, \cos el^i \cos az^i, \sin el^i, 1] }[/math], where [math]\displaystyle{ el^i }[/math] and [math]\displaystyle{ az^i }[/math] are the elevation and azimuth angle of the satellite [math]\displaystyle{ i }[/math] observed from the receiver position, (see Transformations between ECEF and ENU coordinates).
  3. ^ Notice that the [math]\displaystyle{ rms }[/math] equals to the standard derivation [math]\displaystyle{ s }[/math] only when the mean position error is zero. Indeed, let [math]\displaystyle{ m=\frac{1}{n}\sum_{i=1}^n{x_i} \mbox{ ; }rms=\sqrt{\frac{1}{n}\sum_{i=1}^n{x_i^2}} \mbox{ ; } s=\sqrt{\frac{1}{n}\sum_{i=1}^n{(x_i-m)^2}} \qquad \mbox{ thence } \qquad rms^2=s^2+m^2 }[/math]

References

  1. ^ [Misra and Enge, 2001] Misra, P. and Enge, P., 2001. Global Positioning System. Signals, Measurements and Performance.. Ganga-Jamuna Press, Massachusetts,US.
  2. ^ [Diggelen, 2007] Diggelen, F. v., 2007. GNSS Accuracy: Lies, Damn Lies, and Statistics. GPS World. 18(1), pp. 26{32.