Gps Dilution Of Precision⚠ unverified
Aerospace / Navigation · Compute the geometric dilution of precision (GDOP)
Parameters
| Input | Symbol | Unit | Default | Description |
|---|---|---|---|---|
| hdop | hdop | — | 1.0 | Horizontal dilution of precision, dimensionless |
| vdop | vdop | — | 1.0 | Vertical dilution of precision, dimensionless |
| Output | Symbol | Unit | Description |
|---|---|---|---|
| result | GDOP | — | Geometric dilution of precision, dimensionless |
The science & history
Understanding the Parameters
-
DOP is a geometry amplifier, not an error itself — it is a pure number (typically 1–20) that multiplies the per-satellite ranging error (the "user-equivalent range error", a few metres) to estimate the position error: $\sigma_{\text{position}} \approx DOP \times \sigma_{\text{range}}$. DOP of 1 is ideal; under 2 is excellent; above 6 is marginal; above 20 the fix is nearly useless — all from geometry alone.
-
Why geometry matters — a GPS fix trilaterates from the intersection of ranging spheres. If the satellites are spread widely across the sky, the spheres intersect at sharp, well-defined angles and the position is crisp. If they are clustered, the spheres intersect at grazing angles and a small ranging error smears the intersection over a large volume — a large DOP. It is the same reason a survey network with widely-spaced reference points is more accurate than one with points all in a line.
-
The DOP family — different projections of the same geometry:
- HDOP — horizontal (lat/lon) position;
-
VDOP — vertical (altitude), usually worse than HDOP because all satellites are above the horizon (no ranging from below), weakening vertical geometry;
-
PDOP — 3-D position, $\sqrt{HDOP^2+VDOP^2}$ (what this calculator computes);
- TDOP — receiver-clock (time) estimate;
- GDOP — the grand total, $\sqrt{PDOP^2+TDOP^2}$, geometry's effect on the full position-and-time solution.
-
VDOP is the weak axis — because usable satellites all sit above the horizon, vertical geometry is inherently one-sided, so GPS altitude is typically 1.5–3× less accurate than horizontal position. This is why aviation historically distrusted GPS altitude and relied on barometric altitude for vertical guidance until augmentation systems (WAAS/SBAS) tightened VDOP-limited accuracy.
-
DOP changes minute to minute — as satellites rise and set, the geometry — and hence DOP — continuously changes. Receivers report the current DOP so users (and RAIM integrity monitors) know how much to trust the fix; a good receiver picks the satellite subset that minimises DOP.
Derivation (Approaching a Proof)
DOP falls out of the least-squares GPS solution. Linearising the pseudorange equations about an estimated position gives an observation matrix $\mathbf{G}$ (each row a unit line-of-sight vector to a satellite, plus a 1 for the clock term). The covariance of the position/time estimate is
$$\operatorname{cov} = \sigma_{\text{range}}^2\,(\mathbf{G}^\top\mathbf{G})^{-1},$$
so the geometry alone determines the error shape through $(\mathbf{G}^\top\mathbf{G})^{-1}$. The diagonal elements of this matrix are the squared DOPs for each coordinate:
$$HDOP^2 = q_{xx}+q_{yy}, \quad VDOP^2 = q_{zz}, \quad TDOP^2 = q_{tt},$$
where $q_{ii}$ are the diagonal entries of $(\mathbf{G}^\top\mathbf{G})^{-1}$. The composite DOPs are then just sub-sums of these variances (because variances of independent coordinate errors add):
$$PDOP = \sqrt{HDOP^2 + VDOP^2} = \sqrt{q_{xx}+q_{yy}+q_{zz}},$$ $$GDOP = \sqrt{PDOP^2 + TDOP^2} = \sqrt{q_{xx}+q_{yy}+q_{zz}+q_{tt}}.$$
The calculator computes the first line — PDOP — the root-sum-square of horizontal and vertical DOP. To reach true GDOP one would add the $TDOP^2$ (clock) term under the root, which the registry formula omits. The root-sum-square structure is simply the Pythagorean addition of independent error variances. $\qquad\blacksquare$
Dimensional check. $HDOP$ and $VDOP$ are dimensionless geometry factors, so $\sqrt{HDOP^2+VDOP^2}$ is dimensionless ✓ — DOP is a pure multiplier.
History and Development
-
Roots in surveying. Dilution of precision predates GPS: geodesists long knew that the configuration of reference points controls how measurement errors propagate into computed positions. GPS inherited the concept wholesale, with satellites playing the role of reference points.
-
The GPS error budget. As the Global Positioning System matured (declared operational 1995), the position error was cleanly factored into a ranging part (clocks, ephemeris, ionosphere, troposphere, multipath — the "UERE") and a geometric part (DOP). This factorisation, $\sigma_{\text{pos}} = DOP \times UERE$, is the organising principle of GPS accuracy analysis and receiver design.
-
Constellation design and RAIM. Minimising DOP drove the 24-satellite GPS constellation layout (satellites spread so that good geometry is available almost everywhere, almost always) and underlies RAIM (receiver autonomous integrity monitoring), which uses redundant satellites and their geometry to detect faults — essential for aviation use of GPS.
Related Concepts: Barometric Altitude, Great-Circle Distance, Cross Track Error, Magnetic Heading, Circular Orbit Velocity
Notes: This is PDOP, not GDOP — $\sqrt{HDOP^2+VDOP^2}$ is the position DOP; true GDOP = $\sqrt{PDOP^2+TDOP^2}$ adds the clock term the registry omits (misnamed). DOP is a dimensionless geometry amplifier: $\sigma_{\text{pos}}\approx DOP\times$ ranging error. Small = satellites well spread; large = clustered. VDOP > HDOP typically (satellites only above the horizon), so GPS altitude is the weak axis. From $\operatorname{cov}=\sigma^2(\mathbf G^\top\mathbf G)^{-1}$. Defaults ⇒ $\sqrt2$.