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

What 4 parts per million of network strain leaves behind

Drawn here at the parameters it defaults to, with every essay that calls it.
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.

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.

It is drawn by datum-figure with show: "fit-residuals" — one member of a family of 11 figures that share a generator, so the drawing above is what that generator returns when it is asked for this one and given nothing else.

8 essays call it. Every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Changing this changes every one of these figures.

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 the numbers refer to

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.

 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. What the numbers refer to

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.

 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. What the numbers refer to

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.

 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. What the numbers refer to

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.

 Two conventions, and neither of them finds the blunder. A blunder of 350 mm added to leg 3 of a closed traverse, producing a misclosure of 350 mm — one part in 8344, which most specifications would accept. Bowditch's rule shares the misclosure out in proportion to leg length; the Transit rule shares it in proportion to each leg's component along the axis being corrected. Both close the figure exactly, so both are valid; they disagree with each other by up to 14 mm; and neither puts more than 104 mm of correction on the leg that is actually wrong. A rule for distributing a misclosure is a convention for producing consistent numbers, and it is not an instrument for finding errors. Grids, and what a survey does

The two ways to spread a misclosure

Bowditch's rule and the Transit rule take the same closed figure and the same misclosure and disagree about which legs were wrong. Both close it exactly, neither puts the correction on the leg that actually carries the blunder, and no measurement can settle which is right.

 Two parameter sets 100 metres apart, over the region they were fitted to. Each marker is drawn at a size proportional to how far the two transformations put it apart. The second set differs from the first by 100 metres of translation along the direction this network can least see, with the rotations and the scale re-fitted to absorb it — which is what a second agency's adjustment does when it chooses a different constraint. The worst disagreement anywhere in the region is 5.64 metres and the mean is 3.61. Applied at south-eastern Australia the same two sets differ by 193 metres, because the rotation that absorbed the translation here is a rotation of the whole Earth. What the numbers refer to

Two parameter sets, one transformation

Agencies publish seven-parameter datum transformations that differ by hundreds of metres in translation, and the usual reading is that one of them is better. Over the region either was fitted to they are the same transformation: a hundred metres of translation, re-absorbed by the rotations and the scale, moves a British coordinate by 5.6 metres and an Australian one by 193.

 A transformation is more certain in some places than others. The horizontal position uncertainty a Helmert transformation carries, from the stated widths of its own seven parameters, along three meridians. It runs from 40 mm to 68 mm — a factor of 1.71 — and it falls towards the poles, because the rotation terms act on the distance from the Earth's axis. A single figure quoted for "the accuracy of the transformation" is the value at some latitude nobody wrote down. What the numbers refer to

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°.

 What a seven-parameter fit's residual is made of. A published transformation accuracy is the root-mean-square residual at the common points, and here it is 1.51 metres. Almost all of it — 1.51 — is the network's own distortion, which a rigid motion and a scale cannot follow and which is present at every point of the country whether it was used in the fit or not. The transformation's own error, measured as the disagreement between the fitted parameters and the true ones over a clean grid, is 0.091 metres: 6 per cent of the quoted figure. The last bar is the control — the same fit with the distortion switched off, at 2.9e-4 metres. What the numbers refer to

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.

The whole library