Concept

Datum — where it appears

A model of the Earth's shape, placed and oriented, against which a latitude and longitude are numbers rather than positions. Two datums differ by a position, an orientation and a scale, which is seven numbers — each of them the output of a fit and each with a width the tables do not carry.

Named by 33 essays across 6 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
Every projection in the library, measured against both properties. Maximum angular deformation across the bottom, maximum areal error up the side, both on logarithmic scales and both measured over several hundred sample points rather than taken from the projection's description. Conformal projections lie on the left edge, equal-area ones along the bottom, and the corner where both vanish is empty because a projection there would be an isometry of the sphere onto the plane. All 26 are plotted; the four named here are Mercator, Mercator (ellipsoidal), Web Mercator, Gall–Peters.

Web Mercator is not conformal

It carries almost every map on the internet, it is named after the projection whose entire purpose is preserving angles, and it does not preserve angles. The machinery here found that without being told to look.

wrong · Audit
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
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.

What a grid is made of

A projected coordinate declares five things — a datum, a projection, its parameters, a unit and an axis order — and four of them are invisible in the numbers. Every one has a documented case of being guessed wrong, and the cheapest of those cost a hundred metres.

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

wrong · Ellipsoid
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
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.

families · Ellipsoid
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
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.

wrong · Ellipsoid
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
Where the ground is going, six plates. The velocity of the ground under a coordinate, computed as the cross product of each plate's rotation vector with the position, and drawn along the band where that rotation is fastest rather than over an outline this site does not hold. The quickest arrow is 75 millimetres a year on the Pacific plate, which is 1.89 metres in twenty-five years. drawn in Robinson — the arrows' directions are the projection's as much as the plates'.

The epoch is part of the coordinate

A latitude and a longitude with no date on them are incomplete, because the ground they refer to is moving at between ten and seventy millimetres a year. Australia crosses an ordinary survey tolerance in under two years and moved 1.8 metres between two versions of its own datum.

datums · Datum
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
Three surfaces, and the two heights between them. The ellipsoid, the geoid and the ground, with the two heights a coordinate can carry. Ellipsoidal height h is what a satellite fix returns and is measured from a surface defined by four constants. Orthometric height H is what a level and a staff measure and is referred to the geoid — the equipotential surface that best fits mean sea level. They differ by the separation N, drawn here as 45 metres because that is a stated input rather than a computed one: a geoid model is a data product with a truncation degree in it, and this site computes rather than downloads. The arithmetic h = H + N is exact whatever N is.

Height above what?

A satellite reports a height above a mathematical surface. A level and a staff report a height above the surface water settles on. The two disagree by tens of metres, both are correct, and only one of them decides which way a pipe drains.

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

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.

datums · Height

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

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 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 residual has more than one explanation

The method names a projection by fitting every candidate to a set of control points and taking the smallest residual. It has never been asked what else a small residual could be. A map drawn in the right projection from coordinates on the wrong datum leaves a residual of one part in a million — indistinguishable from noise, at every region size, because a similarity fit absorbs a datum shift almost exactly.

wrong · Identify

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.

wrong · Ellipsoid

A coordinate is the output of a solve

Six essays measure a tape, close a traverse, spread a misclosure and reduce a chain. The coordinate that comes out of the far end is the solution of a least-squares problem, and the problem has a decision in it that is not a measurement: what to hold fixed. Change it and every coordinate moves by centimetres while not one residual moves at all.

practice · Reduction

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

The difference of two coordinates

Three essays give a single coordinate a width. Every practical use of one is a difference of two — a distance, a bearing, a movement, an area — and the width of a difference is not the two widths combined, because the errors are not independent. Far from its datum a one-leg baseline is six times more certain than the positions it joins.

distortion · Precision

A rotation is not absorbed the way a shift is

The previous rung expected a datum's rotations to be the part a plane fit could not swallow. They are the part it swallows best — 93 times better than the translations over a hemisphere, and 24 times better per metre moved — and the reason is that the rotation which matters is a change of longitude, which is a symmetry of the map.

wrong · Identify

The seven parameters have their own uncertainty

Nine essays on this ladder print a datum transformation as seven exact numbers. Every published set is the output of a least-squares fit and arrives with standard errors as much a part of the result as the parameters — and pushing those widths through to the ground gives an ellipse, not a number, that is 68 mm across at the equator and 43 mm at 70°.

datums · Datum

A meridian boundary moves when its datum does

The 141st meridian is the one boundary description in this collection with no geometric ambiguity in it: every reading of it is the same curve, exactly. It has a different one. A longitude refers to a datum, the 1825 convention named none, and the line of longitude exactly 141° west sits 129.0 metres apart on NAD27 and WGS84 — a ribbon 1,039 kilometres long and 134 square kilometres in area.

practice · Boundary

What another common point buys

Rung three finds that a seven-parameter datum fit leaves a pattern rather than noise. Six per cent of the residual it reports is the transformation's own error and the other ninety-four is distortion no seven parameters can follow — so adding common points improves a term that was already small and cannot touch the one that is quoted.

datums · Datum

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

A grid stops fitting the ground it was laid on

A hundred-kilometre baseline across a plate boundary changes by 815 millimetres in fifty years, and the same baseline on the fastest plate in the model changes by nothing at all in the same fifty. The interval before an ordinary first-order specification is broken is six years in one place and never in the other.

datums · Strain

Every country's zero is a different surface

Ten rungs measure a height against a geoid and treat the geoid as one object. No national datum is on it: each is pinned to the mean sea level at one tide gauge, the sea surface stands up to two metres from the geoid, and the eight European zeros measured here spread over 462 millimetres — a step no levelling can remove and nobody's error.

datums · Height

The parameters are not independent

Rung seven gives the seven parameters their own uncertainty and stops at seven numbers. There are twenty-eight, and the twenty-one nobody publishes are not small: a translation and the rotation that mimics it correlate at 0.94, the normal matrix has a condition number of 4 × 10¹⁶, and propagating from the diagonal alone overstates the transformation's uncertainty by up to a factor of thirty-six.

datums · Datum

The height a coordinate does not carry

Ten rungs treat a datum shift as a map from one pair of angles to another. It is not one: the transformation runs through Cartesian coordinates, so the answer depends on the height — by 7.7 to 121 millimetres per kilometre depending on the datum and the place, which puts Lhasa 432 millimetres from where a height of zero would have said.

datums · Datum

Named alongside it

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

EllipsoidRealisationHelmert transformationOSGB36ToleranceVerificationConventionWGS84National GridReference frameResidualFlattening

All concepts