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

Ellipsoid — the ladder

15 distinct arguments against one idea, from the one that introduces it to the one that assumes the rest.
  1. 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.

    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.

    rung 1 · wrong
  2. Transverse Mercator. The graticule of the Transverse Mercator projection at 30° of longitude and 15° of latitude. the only projection on which a constant compass bearing is a straight line. It is conformal.

    Transverse Mercator and the series that computes it

    The projection most of the world's survey data lives in has no closed form. What every national grid actually computes is a truncated power series in the flattening, and how far it can be trusted is an engineering parameter rather than a property of the projection.

    rung 2 · families
  3. 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.

    rung 2 · wrong
  4. The sixty zones, each six degrees wide. Every zone is a separate transverse Mercator projection about its own central meridian, so the world is covered by sixty maps rather than one. Zone 31 is picked out, running from 0° to 6° with its axis on 3°. Coordinates do not carry across a zone boundary — a point on either side of one has two entirely different eastings, and nothing in the numbers says which zone they belong to.

    UTM and the zone system

    Sixty separate maps of the world, each six degrees wide, each with a scale factor of 0.9996 chosen so the projection is wrong everywhere and less wrong at the edges. Every constant in the definition is a measured trade rather than a convention.

    rung 3 · families
  5. 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.

    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.

    rung 3 · wrong
  6. Grid north against true north at 52°, across a 6° zone. seven meridians drawn in the grid's own coordinates. The grid's north is straight up the page everywhere by construction; the meridians are not, and the angle between them is the convergence. It is exactly zero on the central meridian, reaches 2.36° — 142 minutes of arc — at 3° out, and changes sign across the middle. The manual's γ = Δλ sin φ is 3.0″ short of it at the zone edge.

    Grid north is not north

    A grid has one north, parallel everywhere on the sheet by construction. The Earth has a different one at every point. The angle between them reaches two degrees at the edge of a UTM zone, and a straight line on the grid is not a straight line on the ground either.

    rung 4 · families
  7. 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.

    rung 4 · wrong
  8. 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.

    rung 5 · datums
  9. 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.

    rung 5 · datums
  10. 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.

    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.

    rung 6 · families
  11. 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 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.

    rung 7 · families
  12. 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.

    rung 8 · wrong
  13. 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.

    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.

    rung 9 · datums
  14. Two expansions of one integral. The worst error along the whole meridian, against the number of sine terms kept, for the series in the third flattening and the classical series in e². Both are checked against a Simpson's rule on the defining integral, which shares no algebra with either. At one and two terms they are the same number to four digits and the e² series is fractionally ahead; from the third term the n series pulls away, and at four it is 1175 times more accurate — 7.6e-8 metres against 9.0e-5.

    Which small quantity the series is in

    Every ellipsoidal formula in this collection is a truncated series in the third flattening, inherited from Krüger in 1912 and justified nowhere. Measured against the classical expansion in e², at four sine terms it is 1,175 times more accurate — and at one and two terms it is fractionally worse, which is not what the folklore implies.

    rung 10 · wrong
  15. 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.

    rung 11 · families

All ladders