Dilution of precision (navigation) explained

Dilution of precision (DOP), or geometric dilution of precision (GDOP), is a term used in satellite navigation and geomatics engineering to specify the error propagation as a mathematical effect of navigation satellite geometry on positional measurement precision.

Introduction

The concept of dilution of precision (DOP) originated with users of the Loran-C navigation system.^[1] The idea of geometric DOP is to state how errors in the measurement will affect the final state estimation. This can be defined as:^[2]

\operatorname{GDOP} = (\Delta ( outputlocation ))/(\Delta ( measureddata ))

Conceptually you can geometrically imagine errors on a measurement resulting in the

\Delta ( measureddata ) term changing. Ideally small changes in the measured data will not result in large changes in output location. The opposite of this ideal is the situation where the solution is very sensitive to measurement errors. The interpretation of this formula is shown in the figure to the right, showing two possible scenarios with acceptable and poor GDOP.

With the wide adoption of satellite navigation systems, the term has come into much wider usage. Neglecting ionospheric ^[3] and tropospheric^[4] effects, the signal from navigation satellites has a fixed precision. Therefore, the relative satellite-receiver geometry plays a major role in determining the precision of estimated positions and times. Due to the relative geometry of any given satellite to a receiver, the precision in the pseudorange of the satellite translates to a corresponding component in each of the four dimensions of position measured by the receiver (i.e.,

x , y , z , and t ). The precision of multiple satellites in view of a receiver combine according to the relative position of the satellites to determine the level of precision in each dimension of the receiver measurement. When visible navigation satellites are close together in the sky, the geometry is said to be weak and the DOP value is high; when far apart, the geometry is strong and the DOP value is low. Consider two overlapping rings, or annuli, of different centres. If they overlap at right angles, the greatest extent of the overlap is much smaller than if they overlap in near parallel. Thus a low DOP value represents a better positional precision due to the wider angular separation between the satellites used to calculate a unit's position. Other factors that can increase the effective DOP are obstructions such as nearby mountains or buildings.

DOP can be expressed as a number of separate measurements:

HDOP: Horizontal dilution of precision

VDOP: Vertical dilution of precision

PDOP: Position (3D) dilution of precision

TDOP: Time dilution of precision

GDOP: Geometric dilution of precision

These values follow mathematically from the positions of the usable satellites. Signal receivers allow the display of these positions (skyplot) as well as the DOP values.

The term can also be applied to other location systems that employ several geographical spaced sites. It can occur in electronic-counter-counter-measures (electronic warfare) when computing the location of enemy emitters (radar jammers and radio communications devices). Using such an interferometry technique can provide certain geometric layout where there are degrees of freedom that cannot be accounted for due to inadequate configurations.

The effect of geometry of the satellites on position error is called geometric dilution of precision (GDOP) and it is roughly interpreted as ratio of position error to the range error. Imagine that a square pyramid is formed by lines joining four satellites with the receiver at the tip of the pyramid. The larger the volume of the pyramid, the better (lower) the value of GDOP; the smaller its volume, the worse (higher) the value of GDOP will be. Similarly, the greater the number of satellites, the better the value of GDOP.

Interpretation

DOP Value	Rating^[5]	Description
< 1	Ideal	Highest possible confidence level to be used for applications demanding the highest possible precision at all times.
1–2	Excellent	At this confidence level, positional measurements are considered accurate enough to meet all but the most sensitive applications.
2–5	Good	Represents a level that marks the minimum appropriate for making accurate decisions. Positional measurements could be used to make reliable in-route navigation suggestions to the user.
5–10	Moderate	Positional measurements could be used for calculations, but the fix quality could still be improved. A more open view of the sky is recommended.
10–20	Fair	Represents a low confidence level. Positional measurements should be discarded or used only to indicate a very rough estimate of the current location.
> 20	Poor	At this level, measurements should be discarded.

The DOP factors are functions of the diagonal elements of the covariance matrix of the parameters, expressed either in a global or a local geodetic frame.

Computation

As a first step in computing DOP, consider the unit vectors from the receiver to satellite

i : \begin{align} &\left( (x _{i} - x)/(R _{i}) , (y _{i} - y)/(R _{i}) , (z _{i} - z)/(R _{i}) \right), & R _{i} & = \sqrt{(x _{i} - x) ^{2} + (y _{i} - y) ^{2} + (z _{i} - z) ^{2} \end{align}}

where

x, y, z denote the position of the receiver and x _{i,} y _{i,} z _{i} denote the position of satellite i. Formulate the matrix, A, which (for 4 pseudorange measurement residual equations) is:A = \begin{bmatrix} (x _{1} - x)/(R _{1}) & (y _{1} - y)/(R _{1}) & (z _{1} - z)/(R _{1}) & 1 \\ (x _{2} - x)/(R _{2}) & (y _{2} - y)/(R _{2}) & (z _{2} - z)/(R _{2}) & 1 \\ (x _{3} - x)/(R _{3}) & (y _{3} - y)/(R _{3}) & (z _{3} - z)/(R _{3}) & 1 \\ (x _{4} - x)/(R _{4}) & (y _{4} - y)/(R _{4}) & (z _{4} - z)/(R _{4}) & 1 \end{bmatrix}

The first three elements of each row of A are the components of a unit vector from the receiver to the indicated satellite. The last element of each row refers to the partial derivative of pseudorange w.r.t. receiver's clock bias.Formulate the matrix, Q, as the covariance matrix resulting from the least-squares normal matrix:

Q = \left( A ^{T} A \right) ^{-1}

In general:

Q = \left(J ^{T} _{x} \left(J _{d} C _{d} J ^{T \right)} _{d} ^{-1} J _{x} \right) ^{-1}

where

J _{x} is the Jacobian of the sensor measurement residual equations f _{i} \left(\underline{x}, \underline{d}\right) = 0 , with respect to the unknowns, \underline{x} ; J _{d} is the Jacobian of the sensor measurement residual equations with respect to the measured quantities \underline{d} , and C _{d} is the correlation matrix for noise in the measured quantities.

For the preceding case of 4 range measurement residual equations:

\underline{x} = (x, y, z, \tau) ^{T} , \underline{d} = \left(\tau _{1} , \tau _{2} , \tau _{3} , \tau _{4} \right) ^{T} , \tau = ct , \tau _{i} = ct _{i} , R _{i} = |\tau _{i} - \tau| = \sqrt{(\tau _{i} - \tau) ^{2} }, f _{i} \left(\underline{x}, \underline{d}\right) = \sqrt{(x _{i} - x) ^{2} + (y _{i} - y) ^{2} + (z _{i} - z) ^{2} } - \sqrt, J _{x} = A , J _{d} = - I and the measurement noises for the different \tau _{i} have been assumed to be independent which makes C _{d} = I .

This formula for Q arises from applying best linear unbiased estimation to a linearized version of the sensor measurement residual equations about the current solution

\Delta\underline{x} = - Q*\left(J ^{T \left(J} _{d} C _{d} J ^{T \right)} _{d} ^{-1} f\right) , except in the case of B.L.U.E. C _{d} is a noise covariance matrix rather than the noise correlation matrix used in DOP, and the reason DOP makes this substitution is to obtain a relative error. When C _{d} is a noise covariance matrix, Q is an estimate of the matrix of covariance of noise in the unknowns due to the noise in the measured quantities. It is the estimate obtained by the first-order second moment (F.O.S.M.) uncertainty quantification technique which was state of the art in the 1980s. In order for the F.O.S.M. theory to be strictly applicable, either the input noise distributions need to be Gaussian or the measurement noise standard deviations need to be small relative to rate of change in the output near the solution. In this context, the second criteria is typically the one that is satisfied.

This (i.e. for the 4 time of arrival/range measurement residual equations) computation is in accordance with [6] where the weighting matrix,

P = \left(J _{d} C _{d} J ^{T \right)} _{d} ^{-1} happens to simplify down to the identity matrix.

Note that P only simplifies down to the identity matrix because all the sensor measurement residual equations are time of arrival (pseudo range) equations. In other cases, for example when trying to locate someone broadcasting on an international distress frequency,

P would not simplify down to the identity matrix and in that case there would be a "frequency DOP" or FDOP component either in addition to or in place of the TDOP component. (Regarding "in place of the TDOP component": Since the clocks on the legacy International Cospas-Sarsat Programme LEO satellites are much less accurate than GPS clocks, discarding their time measurements would actually increase the geolocation solution accuracy.)

The elements of

Q are designated as:Q = \begin{bmatrix} \sigma ^{2} _{x} & \sigma _{xy} & \sigma _{xz} & \sigma _{xt} \\ \sigma _{xy} & \sigma ^{2} _{y} & \sigma _{yz} & \sigma _{yt} \\ \sigma _{xz} & \sigma _{yz} & \sigma ^{2} _{z} & \sigma _{zt} \\ \sigma _{xt} & \sigma _{yt} & \sigma _{zt} & \sigma ^{2 \end{bmatrix}} _{t}

PDOP, TDOP, and GDOP are given by:^[6]

\begin{align} \operatorname{PDOP} & = \sqrt{\sigma ^{2} _{x} + \sigma ^{2} _{y} + \sigma ^{2 }\\} _{z} \operatorname{TDOP} & = \sqrt{\sigma ^{2 }\\} _{t} \operatorname{GDOP} & = \sqrt{\operatorname{PDOP} ^{2} + \operatorname{TDOP} ^{2}\\} & = \sqrt{\operatorname{tr} Q}\\ \end{align}

Notice GDOP is the square root of the trace of the

Q matrix.

The horizontal and vertical dilution of precision,

\begin{align} \operatorname{HDOP} & = \sqrt{\sigma ^{2} _{n} + \sigma ^{2 }} _{e} \\ \operatorname{VDOP} & = \sqrt{\sigma ^{2 } \end{align}} _{u} ,are both dependent on the coordinate system used. To correspond to the local east-north-up coordinate system,

EDOP^2 x x x x NDOP^2 x x x x VDOP^2 x x x x TDOP^2

and the derived dilutions:

\begin{align} \operatorname{GDOP} & = \sqrt{\operatorname{EDOP} ^{2} + \operatorname{NDOP} ^{2} + \operatorname{VDOP} ^{2} + \operatorname{TDOP} ^{2}} \\ \operatorname{HDOP} & = \sqrt{\operatorname{EDOP} ^{2} + \operatorname{NDOP} ^{2}} \\ \operatorname{PDOP} & = \sqrt{\operatorname{EDOP} ^{2} + \operatorname{NDOP} ^{2} + \operatorname{VDOP} ^{2} \end{align}}

See also

• Circular error probable
• GNSS positioning calculation

Further reading

• DOP Factors
• manually calculating GDOP
• HDOP AND GPS HORIZONTAL POSITION ERRORS
• Article on DOP and Trimble's program: Determining Local GPS Satellite Geometry Effects On Position Accuracy.
• Notes & GIF image on manually calculating GDOP: Geographer's Craft
• GPS Errors & Estimating Your Receiver's Accuracy: Sam Wormley's GPS Accuracy Web Page
• GPS Accuracy, Errors & Precision: Radio-Electronics.com

Notes and References

1. Web site: Dilution of Precision . Richard B. Langley . GPS World . May 1999 . 2011-10-12 . 2011-10-04 . https://web.archive.org/web/20111004205627/http://gauss.gge.unb.ca/papers.pdf/gpsworld.may99.pdf . live.
2. Book: Dudek . Gregory . Gregory Dudek . Jenkin . Michael . Computational Principles of Mobile Robotics . registration . . 2000 . 0-521-56876-5.
3. Web site: GNSS and Ionospheric Scintillation: How to Survive the Next Solar Maximum . Paul Kintner, Cornell University . Todd Humphreys . University of Texas-Austin . Joanna Hinks . Cornell University . . July–August 2009 . 2011-10-12 . https://web.archive.org/web/20111106011219/http://www.insidegnss.com/node/1579 . 2011-11-06 . dead.
4. Web site: GPS errors (Trimble tutorial) . 2016-02-08 . 2016-03-07 . https://web.archive.org/web/20160307124113/http://www.trimble.com/gps_tutorial/howgps-error.aspx . live .
5. Isik . Oguz Kagan . Hong . Juhyeon . Petrunin . Ivan . Tsourdos . Antonios . Integrity Analysis for GPS-Based Navigation of UAVs in Urban Environment . Robotics . 25 August 2020 . 9 . 3 . 66 . 10.3390/robotics9030066 . free .
6. https://web.archive.org/web/20141122153439/http://www.gmat.unsw.edu.au/snap/gps/gps_survey/chap1/149.htm Section 1.4.9 of Principles of Satellite Positioning