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

Raster grid — where it appears

A regular array of cells covering a region, each holding one value for a place with a size. Reprojecting one requires reading it at positions its cells do not sit on, which invents values however carefully it is done.

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

What grid spacing a tolerance buys. The worst bilinear interpolation error across OSGB36's ground, against the spacing of the table it is interpolated from, both axes logarithmic. The smooth shift falls by exactly four for every halving, which is what second-order interpolation does. Adding a 30 centimetre ripple 2 degrees across — three parts in a thousand of the 99 metre shift it rides on — moves the spacing needed for 20 millimetres from 1° to 0.125°, which is 64 times as many nodes.

When a formula is not enough

National mapping agencies distribute datum shifts as tables of numbers on a grid rather than as parameters. What decides the spacing that table needs is not the size of the shift — a hundred-metre shift tabulates coarsely, and a thirty-centimetre ripple riding on it does not.

datums · Datum How Equirectangular distorts, by latitude. Angular deformation and areal scale plotted against latitude along the meridian at 0°. On Equirectangular the angular deformation reaches 114.2° and the areal factor reaches 11.5.

The plate carrée, the projection nobody chooses

Plotting latitude against longitude on ordinary axes is a projection. It preserves nothing, its angular deformation reaches 108° and its areal error eightfold, and it is probably the most widely produced map in the world because it is what happens when nobody decides anything.

wrong · Audit One field, one round trip between two cell schemes. Left: a stated field binned into an equal-angle grid of 36 by 18 cells. Right: the same field after being rebinned into an equal-area grid of 30 by 15 offset by six degrees of longitude, and rebinned back. Every step is exact area-weighted averaging, the total is preserved to 2 × 10⁻¹⁶, and the root-mean-square difference between the two pictures is 0.144 on a field whose own standard deviation is 0.370.

The same data on two grids

Five essays have addressed, queried and ordered cells within one scheme and nobody has moved a number between two. Doing it exactly — area-weighted, both directions — preserves the total to 2 × 10⁻¹⁶ and loses 39 per cent of the field's own standard deviation in a single round trip; six round trips leave 23 per cent of its variance. The quantity that would reveal the damage is the one that never moves.

applied · Cells

Named alongside it

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

InterpolationAddressAggregationCellConservationData visualisationDatumDefaultDiscrete global gridEqual-areaEquirectangularGrid shift

All concepts