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

OSGB36 — where it appears

The Ordnance Survey Great Britain 1936 datum, on the Airy 1830 ellipsoid, which the National Grid is defined on. Its transformation to WGS84 is about a hundred metres and is published both as seven parameters and as a shift grid, because the network's own strain will not fit seven numbers.

Named by 13 essays across 4 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 OSGB36 to WGS84, one parameter at a time. Each bar is how far the mark at 2.0° west, 54.5° north moves under one of the seven parameters with the other six set to zero. seven of the seven are not zero for OSGB36. The units hide the comparison: one arcsecond of rotation moves this point 30.8 metres and one part per million of scale moves it 6.36 metres, so OSGB36's scale term contributes 130 metres — more than one of its three translations.

The seven parameters, and what each one does

A datum transformation is published as three metres, three arcseconds and one part per million, in three different units — which hides the fact that all seven are the same size of thing. Converted to ground displacement, the scale term beats two of the translations.

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 What 4 parts per million of network strain leaves behind. Each arrow is where the best-fitting seven-parameter transformation leaves a marker, over 36 markers laid out as a grid across OSGB36's ground. The RMS residual is 1.64 metres and the worst is 3.00 metres. Seven parameters span the constant and linear parts of a displacement field; this one is quadratic, so no choice of the seven can reach it. The arrows are drawn 27× life size.

Where a fit leaves residuals

Seven parameters can carry a rigid motion and a size exactly. A triangulation network is neither, so the best possible transformation between two datums leaves metres on the table — in a pattern, not as noise — and which seven parameters come out depends on where the markers were.

datums · Datum 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 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.

wrong · Ellipsoid 53° north, 2° west: four places one pair of numbers could be. The same two numbers, read as three datums and as the pair in the other order. The centre is the WGS84 reading; the marks are where the ground is if the numbers were measured against OSGB36, ED50, NAD27 instead, at 103 m, 138 m, 195 m. Every one of those is a plausible position — near enough to look like a survey difference and far enough to be a different field. Reversing the two numbers instead moves the point 7942 kilometres, and is the mistake nobody worries about because it is usually obvious.

A coordinate without its system is not a location

The same two numbers name ground 103 metres apart on one datum, 195 on another, and 7,942 kilometres away if the pair is read in the other order. The datum errors are the dangerous ones because they are plausible — and the axis swap, which everyone calls obvious, does nothing at all along a line that crosses three continents.

applied · Dataset 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, and the height changes by 51, 44, 47 metres, which is the component nobody quotes.

The third coordinate moves too

A datum shift is quoted as a horizontal displacement because horizontal is what people look at. The transformation acts on a three-dimensional point, and its vertical component is between a quarter and a half of the horizontal one — 51 metres, for a British coordinate.

datums · Height Only one of a grid's five parameters changes the map. The same ground point written on the British National Grid and on UTM zone 31N, with the difference between the two coordinate pairs taken apart. The largest term by four orders of magnitude is where the two grids put zero — 5544 km, being a different central meridian, a different true origin and different false constants, all of which move every coordinate and no distance. The datum, which is the term everybody names, moves the ground 106 m. And the whole geometric difference between two transverse Mercators is their scale factors — 0.9996012717 against 0.9996 — which at this point is 50 cm. Four of a grid's five parameters are bookkeeping; the fifth is the map.

Two grids over the same ground

One point in the English Midlands is 406,788 east on the British grid and 167,478 east on UTM. Separating the difference properly leaves a surprise — only one of a grid's five parameters changes the map, and it is not any of the four that dominate the number.

practice · Reduction A published coordinate is a definition being quoted. Re-observing a control point perfectly, with an instrument good to 20 mm, and comparing the answer to what is published for it. The published coordinate was fixed on a datum realisation that has since been superseded, which accounts for 106 m, and the ground it marks has moved 750 mm in 30 years at 25 mm a year. Both numbers are larger than the observation and neither is an error in it: subtracting a published coordinate from an observed one measures the interval between two definitions.

A published coordinate is a result

Re-observe a control mark perfectly, with an instrument good to twenty millimetres, and the answer disagrees with the published value by a hundred metres. Neither number is wrong. The difference measures the interval between two definitions.

practice · Reduction A 1-hectare parcel across British National Grid, at 52°N. The departure of a plan's dimensions from the ground's, across the width of the zone, as a length and as an area. The area curve is the length curve doubled — an areal scale factor is the square of a linear one, and squaring a small departure doubles it — so a parcel whose sides are each 399 parts per million short on the central meridian is 797 short in area. On a 1-hectare parcel that is 8.0 square metres at the central meridian and 8.0 at 0.00° out.

An area on the grid is not an area on the ground

A grid's scale factor is a property of lengths and what a surveyor sells is an area. Squaring a departure doubles it, so a hectare drawn on the British grid at its central meridian has 10,008 square metres of ground under it — and correcting an area with the line factor instead of its square leaves half the error behind.

practice · Grid A plane fit swallows a datum shift, and keeps swallowing it. A map of OSGB36's ground drawn as though it were on WGS84, at seven sheet sizes. The upper curve is how far the drawing moves — a real 99-metre error on the ground — and the lower one is what survives the best similarity between the two, which is what a fit for the projection removes for free. The absorption runs from 99.68 per cent on a 1° sheet to 85.16 on a 60° one, and what is left is 3.9 parts per million of the map's own size even then — two microns on a sheet half a metre across.

The datum hides inside the projection's parameters

Six essays recover a projection from control points with the body it was drawn on taken as known. It is not known, and the fit cannot find it: a 99-metre datum shift is absorbed to 99.7 per cent on a two-degree sheet and to 85 per cent on a hemisphere, leaving eleven parts per million of the map behind.

wrong · Identify What reading the seven parameters under the wrong convention costs. The distance between the two published conventions' answers for the same parameter set, at the worst point of the world. The difference is 2s(r × X) exactly — twice the rotation, crossed into the position — so a parameter set with no rotations is immune and one with a large rotation is not: DHDN reaches 152.3 metres and NAD27, whose published transformation is three translations and nothing else, reaches zero. The note beside each bar is the size of that datum's rotation.

The rotation has two sign conventions

Eleven essays price a datum transformation's parameters, their fit, their residuals and their uncertainty. None asks what the numbers mean: three of the seven are published under two conventions whose rotations differ in sign, and reading one set with the other formula costs exactly twice the rotation — 55 metres on OSGB36 and 152 on DHDN.

datums · Datum

Named alongside it

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

DatumRealisationHelmert transformationNational GridToleranceVerificationWGS84ConventionEllipsoidED50EpochResidual

All concepts