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

Radius of curvature — where it appears

How sharply a surface bends at a point in a stated direction, of which an ellipsoid has two different ones everywhere except at its poles. The meridian and prime-vertical radii differ by about a per cent on the Earth, and the difference is why a normal section is not a geodesic.

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

The same coordinate on four datums. One pair of numbers — 2.0° west, 54.5° north — read as a coordinate on OSGB36, ED50, NAD27 and on WGS84, and plotted where each reading puts the mark on the ground. The spread runs to 195 metres. The numbers are identical; only what they refer to differs.

What a coordinate refers to

A latitude and a longitude are not a place. They are a claim about a model, and the model has a shape, a position, a set of marked stones and a date — change any one of the four and the ground under the numbers moves by hundreds of metres.

datums · Datum Fitting a sphere to 10° of latitude at 54.5°. The Earth's Gaussian radius of curvature across the band, with two candidate spheres drawn against it. The global mean radius is 6371.0 kilometres and sits 2203 parts per million away from the ground here — 2.2 metres in every kilometre measured. The best local radius is 6385.0 kilometres and has no bias at all by construction, leaving 325 parts per million of spread that no sphere can remove, because the curvature varies across the band and a sphere's does not. The gain is a factor of 6.9.

A datum is fitted to a region

Forty countries adopted forty ellipsoids, and the usual explanation is that measurement was poor. It was not — a figure fitted to one country beats the global one over that country by a factor of seven, and what it removes is a systematic bias of two metres in every kilometre.

datums · Datum One pixel at zoom 11: 76.4 m projected, 47.6 m at 51.5°. The grid squares are pixels, at their own size. The open mark is the stored coordinate and the filled one is where it is drawn: rounding moves it 20.2 metres, against a worst case of 33.6 — half a pixel's diagonal, which is the whole of the bound. The second point sits 43 metres away, 0.58 of a pixel east and 0.70 north, so whether the two are drawn as one dot or two is decided by where the tile grid happens to fall: they merge for 12 per cent of the possible offsets, against the 12 per cent the two fractions predict.

The pixel is a place with a size

Drawing a coordinate rounds it to a pixel, which moves it by up to half a diagonal — 16.8 metres at zoom 12. Whether two points 43 metres apart appear as two dots is not a property of the data at all: they merge for 12 per cent of the positions the tile grid could take, and the closed form predicts 12.4.

applied · Screen 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 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 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.

impossibility · Curvature "Within 0.01°" on the ground, at five latitudes. The same condition drawn at five latitudes, all at one scale. A degree of latitude is a fixed distance — 1106 metres here, varying by less than one per cent from equator to pole — while a degree of longitude collapses as cos φ, from 1113 metres to 289 at 75°. So the "circle" is an ellipse of 3.86:1 there, and it encloses 26 per cent of the ground the same condition covers on the equator. Even on the equator it is not round: M is smaller than N by the flattening, so the shape is 6694 parts per million shorter north–south than east–west.

A degree is not a unit of length

"Within 0.01 degrees" is a condition anybody can write and no instrument can measure. On the ground it is an ellipse — 1,106 metres north–south and 558 east–west at 60° — and even on the equator it is not a circle, because the meridian's radius of curvature is smaller than the parallel's by the flattening.

applied · Dataset How much curvature varies from place to place, body by body. The Gaussian curvature along a meridian, each body against its own smallest value, so the curves are comparable. Written from the bodies actually drawn: Earth varies by 1.0135, Mars varies by 1.0239, Vesta varies by 2.71, Phobos varies by 4.16, from the flattest place on each to the sharpest. A surface of constant curvature is a sphere and nothing else is, so the impossibility argument this field is built on has a local version on every real body: how flat a patch is depends on where the patch is.

Curvature that varies from place to place

The impossibility this site is built on was argued on a sphere, where the curvature is one number. On the Earth it varies by 1.35 per cent, on Jupiter by 31, on Vesta by a factor of 2.7 — and on a body with three axes it varies along a parallel, which no formula in latitude can express.

impossibility · Curvature Four steps between a tape and a drawing. A slope distance of 76.895 km measured at 3.2° on ground 220 m above sea level, reduced to the British National Grid, with every step drawn on a logarithmic scale in millimetres. The slope reduction is the largest by a wide margin at 119.90 m and is the one everybody applies. The grid reduction is 30.23 m. The fourth bar is not a step at all: it is what using the levelled height where the height above the ellipsoid is wanted costs, with a geoid separation of 48.5 m — 583 mm, which is 1.9% of the grid reduction it sits beside and 29 times the tolerance the job closes to. Reading a correction's importance off its share of the chain is how it gets dropped.

What a tape measures

Four steps stand between an instrument reading and a coordinate, and their sizes are not in the order anybody expects. On a 77-kilometre line the slope reduction is 120 metres, the grid reduction 30, and a step nobody names — using the height above sea level where the height above the ellipsoid is wanted — is 583 millimetres.

practice · Reduction 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.

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.

wrong · Ellipsoid The route at height does not lie above the route on the ground. London to Tokyo, solved at the surface and again at a stated height, with the two ground tracks compared point for point. At a cruising altitude of eleven kilometres the two part by 11.8 metres; the departure is proportional to the height, at a fitted slope of 0.997. On a sphere the same measurement returns 2.7e-9 metres, because the offset of a sphere is a sphere and the two geodesics coincide exactly. The offset of an ellipsoid is not an ellipsoid, and this is what that costs.

The shortest route is not at sea level

Ten rungs route on a surface and nothing is ever flown on one. The offset of a sphere is a sphere, so at altitude the great circle is the great circle. The offset of an ellipsoid is not an ellipsoid — its radii of curvature are M + h and N + h, which belong to no ellipsoid — so the shortest route at cruising height does not lie above the shortest route on the ground.

paths · Paths The two reductions, at a grid factor of 1.0004. A measurement made on the ground has to be brought to the ellipsoid and then to the grid, and the two corrections have opposite signs. The elevation factor is R/(R+h) and always shrinks; the grid factor here is 1.0004, which stretches. Their product is the only number a surveyor can use. They cancel exactly at 2551 metres. At 1500 metres the combined factor is 165 parts per million, which is 1.65 metres on a 10 kilometre baseline.

The ground is not the grid

A tape measure on a hillside has to be brought down to the ellipsoid and then out onto the map, and the two corrections have opposite signs. On a grid whose scale factor exceeds one there is exactly one elevation where they cancel — 2,551 metres, for a factor of 1.0004.

datums · Height Which corrections a 10 mm job may leave out. Each correction inverted: the line length — or, for the last row, the patch radius — at which that correction alone reaches 10 mm, for a job 150 m above the ellipsoid on 2° of ground, 0.0° from the British National Grid's central meridian. The tightest is the slope reduction at 16 m. Three of the four rows are linear in the tolerance, so this ranking is the same at a millimetre and at a decimetre; what does change it is the site — move the job up a mountain or onto steeper ground and the order is different, which is why a specification's list is not transferable.

The tolerance decides the model

Every correction inverts to a distance — the line length at which it alone exceeds a stated tolerance. The ranking those distances produce is the same at a millimetre and at a decimetre, and it is different for two jobs at the same tolerance on different ground.

practice · Reduction Clairaut's theorem, and the term it drops. Clairaut's theorem says the flattening of a rotating body plus the flattening of the gravity on it equals five halves of the ratio of centrifugal to gravitational acceleration at the equator. Measured on WGS84: f = 3.3528e-3, f* = 5.3024e-3, and their sum is 8.65525e-3 against the theorem's 8.62447e-3. The residual is 3.08e-5, which is 2.74 times f² — second order, which is what a first-order theorem is entitled to be wrong by.

The flattening is not a free parameter

An ellipsoid is usually presented as two numbers somebody fitted. One of them is not free — Clairaut's theorem relates the shape of a rotating body to the gravity on it, and the relation holds on WGS84 with a residual of 3.1×10⁻⁵ — which is 2.74 times f², exactly what a first-order theorem is entitled to.

datums · Ellipsoid 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 The signal, and four instruments' noise. The spherical excess of an equilateral triangle at 45° north against its side, with the standard deviation of a measured excess — σ√3 — ruled for four instrument accuracies. A fifty-kilometre triangle, which is about the largest anybody routinely observed, has an excess of 5.49 seconds of arc. A theodolite reading to one second gives that excess a standard deviation of 1.73 seconds, so the measurement carries about three significant bits. Everything in this rung follows from that ratio.

How big a triangle it takes

Gauss proved that a surface-dweller can read the curvature off a triangle's angles. Doing it is another matter: a fifty-kilometre triangle has 5.49 seconds of excess, a one-second theodolite gives that excess a standard deviation of 1.73, and reading K to one per cent needs a side of 281 kilometres. Not one of the great surveys built a triangle within a factor of three of that.

impossibility · Curvature Four radii of the Earth, and one that is a range. The four constants called the mean radius of the Earth, on a scale of kilometres, with the range of the local Gaussian radius √(MN) drawn behind them. Three of the four agree to about a part per million; the rectifying radius is 3560 metres smaller, which is 559 parts per million. The Gaussian radius spans twelve times that range on its own, which is why there is no such thing as the conformal sphere.

Four radii of the Earth

Ten rungs handle the ellipsoid with an auxiliary latitude. Every one of those constructions also needs a radius, and the radius that makes each property exact is a different number: the published 6371 km is right for an area to half a part per million and wrong for a meridian distance by 559 — while the radius a conformal map needs is not a constant at all, and spans 6,739.

families · Ellipsoid

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.

datums · Ellipsoid

Named alongside it

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

EllipsoidToleranceVerificationFlatteningScale factorSpherical approximationEllipsoidal heightGaussian curvatureNational GridWGS84Closed formCombined factor

All concepts