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

Five latitudes that are not the latitude, on WGS84

Drawn here at the parameters it defaults to, with every essay that calls it.
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.

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.

It is drawn by ellipsoid-figure with show: "latitude-ladder" — one member of a family of 12 figures that share a generator, so the drawing above is what that generator returns when it is asked for this one and given nothing else.

12 essays call it. Every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Changing this changes every one of these figures.

Why there are two latitudes, at 45°. A meridian section of an ellipsoid with the flattening exaggerated 12× — at the true value of 1/298 this outline would be indistinguishable from a circle. Geodetic latitude is the angle the surface normal makes with the equatorial plane; geocentric latitude is the angle the radius makes with it. The two lines meet the plane at different points and the angles differ by 11.55 arcminutes on the real WGS84 ellipsoid, which is about 21 km of ground distance. What is taught wrongly

The Earth is a sphere, and when it is not

Every essay before this one treated the Earth as a ball, and said so. The flattening is one part in three hundred, which is nothing for a distance, everything for a latitude, and exactly enough to make the most-used projection in the world fail the property in its own name.

 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. What the numbers refer to

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.

 What treating the Earth as a sphere costs, per journey. The ellipsoidal geodesic minus the spherical great circle, in kilometres, for five journeys. The correction is a few tenths of a per cent and it changes sign: a route running east–west at mid latitude is longer on the ellipsoid, and one running along a meridian is shorter, because an oblate body is fatter round the equator and flatter pole to pole. Both distances are computed — Vincenty's iteration against the haversine formula — and the ellipsoidal one is checked against a geodesic obtained by integrating its own differential equation. Paths and directions

Geodesics on the ellipsoid, and why they are hard

The shortest path on a flattened Earth is not a plane curve, has no closed form, and can be longer or shorter than the spherical answer depending on which way it runs. Every practical method is a series or an iteration, and the correction changes sign.

 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. What is taught wrongly

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.

 What changing the datum alone does to a coordinate. The distance on the ground between a point as read on its national datum and the same numbers read on WGS84, computed through the published seven-parameter transformation. The shifts run from 49 to 166 metres. For comparison, the scale error a UTM zone introduces at its edge is under a metre per kilometre — so the datum, which is usually left unstated, dominates the projection, which is usually argued about. What is taught wrongly

Datum shifts dwarf projection errors

The projection argument is conducted in parts per million and the datum question is answered in hundreds of metres. A coordinate whose datum is unstated is out by more than any projection choice could ever put it, and almost nobody checks.

 The radius of equal curvature, from equator to pole. The radius of the sphere that has the same Gaussian curvature as the ellipsoid does, at each latitude. On WGS84 it runs from 6357 km at the equator to 6400 km at the pole — the polar region is the flattened part of a squashed ball and is therefore the LEAST curved — so K itself varies by 1.35%. A surface whose curvature varies cannot be laid on one whose curvature does not, so no map from the sphere to the ellipsoid is faithful either, and the least a conformal one can vary in scale is 6739 parts per million. The impossibility

The curvature of the Earth is not one number

An ellipsoid's Gaussian curvature varies by 1.35% from equator to pole, so the sphere is not developable onto the ellipsoid any more than the plane is onto the sphere. The spherical approximation is a projection with an irreducible cost, and the cost is 6,739 parts per million.

 Three curves between 45°N and 55°N, 3026 km apart. The normal section observed from the first point, the normal section observed from the second, and the geodesic, each plotted as its distance to one side of the great-circle chord between the two ends. The two sections are 79.0 metres apart at their widest and the geodesic runs between them, 53.3 metres from the first. In LENGTH they differ by almost nothing — 2.40 millimetres over 3026 kilometres — so an observation that follows the wrong one measures the right distance along the wrong ground. On a sphere all three coincide exactly. What is taught wrongly

The normal section is not the geodesic

A theodolite at A sighted on B swings in one plane and the instrument at B sighted back swings in another, so the two observations trace different curves on the ground — 79 metres apart over 3,026 kilometres — and the shortest path is neither of them. The lengths differ by 2.4 millimetres, so the wrong curve measures the right distance along the wrong ground.

 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. What the numbers refer to

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.

 The boundary of the region the series is the map in. The curve where the transverse coordinate reaches 2.918, which is where the terms stop shrinking. It crosses the equator 83.81° from the central meridian and closes towards the poles, because the same longitude is a smaller transverse coordinate at a higher latitude — the boundary is a curve rather than a meridian. The narrow band beside it is a 3° zone, the width national grids actually use, drawn to the same scale: the practical world sits in about a fiftieth of what the series can reach. Drawn in Mollweide. The families

Where the series stops being the map

The transverse Mercator has no closed form on an ellipsoid, so every national grid computes a truncated series. Asking how many terms it needs has an answer everywhere; asking whether more terms always help has an answer only within 83.81° of the central meridian, and the boundary is computed rather than assumed.

 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. What is taught wrongly

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.

 Going out and coming back, at three orders. A point three degrees from the central meridian, taken forward into the projection and back again with the series truncated at the same order both ways, and the ground distance between where it started and where it returned. The second-order pair is out by a tenth of a metre at some latitudes; the third by half a millimetre; the fourth — the order every national grid formula in ordinary use is written to — by 0.39 mm. The dips are where one of the two series passes through a node, not where the map is better. The families

The inverse of the series is not the series of the inverse

Transverse Mercator on an ellipsoid has no closed form, so every national grid computes it as a truncated series. There are two of them — one out and one back — and they are separate approximations. Composing them does not give the identity: at the order every grid formula in ordinary use is written to, a point three degrees from the central meridian comes back 0.39 mm north of where it started.

 There is no such thing as the radius of the Earth at a latitude. The two principal radii of curvature of WGS84, against latitude, with the mean radius that every table prints drawn across them. The meridional radius M runs from 6335.44 km at the equator to 6399.59 at the pole; the prime-vertical radius N from 6378.14 to the same value. They differ by 42.70 km at the equator and meet only at the pole, and the single number 6,371 km lies between them at no latitude where either is right. What the numbers refer to

The radius of curvature is two numbers

Twelve essays work on the ellipsoid and every one of them takes a radius when it needs one. There are two at every point, they differ by 42.70 kilometres at the equator, and the single number every table prints is out by 5,583 parts per million on a line running north.

The whole library