Cartographic Projection
every map is wrong, and by exactly how much
Concept

Geodetic latitude — where it appears

The angle between the equatorial plane and the normal to the ellipsoid at a point, which is the latitude every map and every receiver reports. It differs from the geocentric latitude by up to about eleven arcminutes, at 45 degrees, on the Earth.

Named by 12 essays across 2 fields — each of them below, with the objects they name alongside it.

Molodensky against the exact route, NAD27. The distance on the ground between where each shortcut puts the transformed point and where the exact Cartesian route puts it, at 0 metres of ellipsoidal height, on a logarithmic scale. The full formulae stay within 0.8 centimetres across every latitude drawn; the abridged form, which replaces two ellipsoid-difference coefficients with one combined term, is worst at 40 centimetres — a factor of 48. The abridged curve dips near 45°, where the combined coefficient happens to equal the pair it replaces.

Molodensky's shortcut

The exact way to shift a datum needs an iteration nobody wanted to run on a 1950s machine, so a direct formula was derived instead. It lands within centimetres — and the abridged version everybody quotes is a different formula, worse by a factor of fifty.

datums · Datum Five latitudes that are not the latitude, on WGS84. Each curve is the amount by which one auxiliary latitude falls below the geodetic latitude a coordinate actually means, in arcminutes. All five vanish at the equator and at the poles and peak near 45°, where the geocentric latitude is 11.55 arcminutes below the geodetic one — about 21.4 km on the ground. The curves never cross, which is forced by the algebra rather than by this ellipsoid's particular flattening.

Geodetic against geocentric latitude

Two angles, both called latitude, differing by eleven and a half arcminutes at their worst. One is what every coordinate means and the other is what every spherical formula assumes, and the gap between them is twenty-one kilometres on the ground.

wrong · Ellipsoid Where the ground is going, six plates. The velocity of the ground under a coordinate, computed as the cross product of each plate's rotation vector with the position, and drawn along the band where that rotation is fastest rather than over an outline this site does not hold. The quickest arrow is 75 millimetres a year on the Pacific plate, which is 1.89 metres in twenty-five years. drawn in Robinson — the arrows' directions are the projection's as much as the plates'.

The epoch is part of the coordinate

A latitude and a longitude with no date on them are incomplete, because the ground they refer to is moving at between ten and seventy millimetres a year. Australia crosses an ordinary survey tolerance in under two years and moved 1.8 metres between two versions of its own datum.

datums · Datum Four ways to build the same equal-area projection. The cylindrical equal-area projection's own areal scale factor, measured from its Jacobian against the metric of the body it is drawn for. On a sphere with the spherical formula it is one everywhere, which is the control. Feed the same formula the geodetic latitude a coordinate actually carries and measure against the ellipsoid it refers to, and it is 1.00674 on the equator and 0.99332 at 88° — a spread of 1.34 per cent, on a projection whose entire purpose is that there is no spread. Rescaling to make the totals agree does not repair it. The authalic northing q/2 does, exactly.

Equal-area on the wrong body

The site's headline is that Web Mercator puts geodetic latitudes into a spherical conformal formula and stops being conformal. The same sentence is true with "equal-area" in it and nobody says it: the areal factor is 1.00674 at the equator, 0.99332 at 88°, and averages to almost exactly one — so every check that adds up areas passes while every cell is wrong.

wrong · Audit A deflection of 10 arcseconds, and what it hides. The ellipsoid normal and the plumb line at one point, with the geoid tilted against the ellipsoid by 10 arcseconds — drawn 3000× steeper than life, because at true scale the two lines are indistinguishable. The relation is exact and linear: an arcsecond of deflection is the geoid rising 4.85 millimetres in a kilometre, so 10 arcseconds is 48.5 millimetres per kilometre. A star sight measures the plumb line's direction, so astronomic latitude differs from geodetic by exactly this angle — 309 metres of ground, at a point where the coordinate itself is correct.

The plumb line is not the normal

A latitude measured from the stars and a latitude that means a position on the ellipsoid are different angles, because a plumb bob hangs along gravity and gravity is not perpendicular to a mathematical surface. Ten arcseconds of difference is 309 metres of ground.

datums · Height How badly two arcs determine the flattening. The reciprocal flattening recovered by inverting two measured degree lengths, at 1.5° and 66.33°, against an error introduced into the equatorial one. The exact pair returns 298.26. Ten metres of error — 90 parts per million of a 110-kilometre arc — returns 301.5, and the relative error in the flattening is 118 times the relative error in the arc. At 1000 metres the inversion returns a negative flattening: an Earth longer through the poles than across the equator, which is the answer the Paris Observatory defended for a generation.

The figure of the Earth was measured

Two expeditions went to Lapland and Peru to measure the length of a degree of latitude, and the whole signal separating a flattened Earth from a round one is a kilometre in a hundred and eleven. Inverting two arcs amplifies their error by 118 — and a kilometre of error returns a lemon-shaped planet.

datums · Ellipsoid Two latitudes for one point on Mars. A cross-section of Mars with its flattening exaggerated 6× so the two angles can be told apart. The line to the centre defines the planetocentric latitude and the outward normal defines the planetographic one; at 45° they differ by 0.338° on the real body, which is 20 kilometres along the surface. Both numbers are published for Mars, and a coordinate that does not say which it is is ambiguous by that much.

A coordinate on another body

Mars publishes two latitudes for every point and they differ by up to 0.338°, which is twenty kilometres of ground — almost exactly the same distance as the Earth's own 0.192°, because Mars is smaller by nearly the same factor. The libraries that compute either have known both numbers since this collection's early essays and had never been pointed anywhere but here.

datums · Bodies The equator, with its ellipticity exaggerated forty thousand times. The dashed circle is the equator every projection formula on this site assumes. The solid curve is the equator satellite geodesy reports, drawn with its departure multiplied by 40,000 so that seventy metres on a six-thousand-kilometre radius can be seen at all. The long axis is at 14.9° west and the short one ninety degrees from it, and the difference between them is 70.0 metres — a real quantity, about the height of a twenty-storey building, on a body every geodetic computation treats as a surface of revolution.

The equator is not a circle either

Eleven rungs price what pretending the Earth is a sphere costs, and every one of them replaces the sphere with a surface of revolution — a body whose equator is a circle. It is not. The two equatorial radii differ by seventy metres, the two surfaces part by thirty-five, and the auxiliary latitudes every ellipsoidal formula is written in stop existing.

wrong · Ellipsoid Four spheres are wrong about a direction before anything has moved. At 45° of latitude, the worst error over twenty-four starting directions of the azimuth each substituted sphere gives for a geodesic on WGS84, against the distance to its end. Geodetic latitude is wrong by 347″ at one kilometre, and the parametric, authalic and rectifying latitudes by 174″, 116″ and 86.8″ — a half, a third and a quarter of it — flat until a few hundred kilometres. The conformal latitude is right at the start and wrong by 54.3 mas at one kilometre, 5.43″ at 100 and 54.3″ at 1,000; the geocentric latitude starts at 620 mas and joins it. Bessel's sphere, which also corrects the longitude, is exact at every distance to the arithmetic's floor of 1.2e-7″.

A bearing on a sphere is decided by its latitude, not its radius

A sphere that stands in for the ellipsoid cannot be wrong about a direction by being the wrong size, because enlarging a sphere moves no angle. What it can be wrong by is the latitude put on it. Geodetic latitude turns every bearing at 45° by up to 347 arcseconds before anything has moved; the conformal latitude starts exact and drifts 54 arcseconds in a thousand kilometres; and only Bessel's sphere, which changes the longitude too, is exact at every distance.

wrong · Ellipsoid An angle is a difference of two directions, and the error does not cancel between them. A geodesic triangle with sides of 30 km at 45° N, its first side laid off on a bearing of 45° and the second 60° round from it. Substituting a sphere carrying the geodetic latitude turns each direction at a corner by an amount that runs as the sine of twice its azimuth, so the two directions at A are turned by different amounts and the angle between them is wrong by their difference. Here the worst single direction is out by 345.79″ and the worst angle by 521.09″ — larger, not smaller. Bessel's sphere, which also corrects the longitude, keeps every angle to 6.2e-9″.

An angle is a difference, and the difference doubles the error

Substituting a sphere for the ellipsoid turns every direction at a point, and a surveyor measures angles rather than directions — so the obvious hope is that a turn common to both directions cancels in their difference. It does not cancel, because the turn is not common: it runs as the sine of twice the azimuth, and the largest angle error is exactly twice the largest direction error at a well-shaped corner. A first-order triangulation computed on the geodetic-latitude sphere carries angles nearly nine minutes of arc wrong, and its closure check passes perfectly.

wrong · Ellipsoid Where the correction overtakes the sphere it is correcting for. The worst angle error of four substituted spheres against the size of the triangle, with Legendre's spherical excess on the same axes. The two that are wrong at first order in the flattening are flat in the side, because the error belongs to the latitude and not to the distance. The two that are right at first order have no such floor, so what is left is a term rising in proportion to the side. The excess rises as its square, so each pair meets exactly once, and the meeting decides which of the two corrections a computation is entitled to ignore: the geodetic latitude at 490 km, the rectifying latitude at 247 km, the conformal latitude at 23.8 km, the geocentric latitude at 17.6 km.

The correction is the smaller of the two corrections

Legendre's theorem subtracts a third of the spherical excess from each angle of a triangle and solves the rest as a plane figure. It is the correction every nineteenth-century computer applied and the one every textbook explains. Beside the error the choice of substituted sphere leaves in the same triangle it is small: on the thirty-kilometre triangle a first-order chain is made of, omitting it costs 166 millimetres in a computed side and using the geodetic-latitude sphere costs 43 metres. The two terms scale differently, so each sphere has a size at which the correction overtakes it — and for the two spheres a careful computer would have chosen, that size falls inside a real chain.

wrong · Ellipsoid Along a meridian a good sphere's error adds up triangle by triangle; a bad base does not. A strip of thirty-kilometre equilateral triangles running north from 45°, each side computed from the last by Legendre's theorem on a substituted sphere. The scale error of the side just computed, against how many triangles lie between it and the measured base. The conformal sphere adds 7.90 ppm a triangle and reaches 158.1 ppm after 20; the geocentric, 158.9 ppm; the authalic, −52.4 ppm; the rectifying, 0.29 ppm. Shaded: three standard deviations of the scale error one-second angles leave, 3.96 ppm a triangle, growing as the square root of the count. Dashed: a base measured 8 ppm long, which stays 8 ppm at every triangle.

A strip of triangles carries the sphere's error to the next base

One triangle cannot tell a substituted sphere from a bad base: both leave the computed side a few parts per million long. A triangulation is not one triangle but a strip, each computed side the next one's base, and along it the two come apart completely. A bad base is one factor, the same at every triangle. A sphere's error is made afresh in each triangle and the strip adds it up, and what it adds up to at the far end is exactly the change in that sphere's own meridian scale between the two ends — 158 parts per million for the conformal sphere over twenty triangles north from 45°, predicted to half a part. Along a parallel it is nothing; at a check base three triangles on it is already three times the noise of one-second angles.

wrong · Ellipsoid

Named alongside it

The objects these essays reach for when they reach for this one.

Auxiliary latitudeVerificationEllipsoidToleranceDatumFlatteningConformal latitudeSpherical approximationGeodesicSpherical excessSurface normalSurvey network

All concepts