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

The thread: What survives a change of coordinates — page 1

Scale along the meridian and scale along the parallel depend on how the sphere was parameterised. The principal scales, the areal factor and the angular deformation do not, and only those are properties of the map itself. Essays 1 to 24 of 110.
The indicatrix at 45°, 55° on Mercator. A circle on the sphere becomes this ellipse on the map. Its semi-axes are the two principal scale factors, a = 1.743 and b = 1.743; their product is the areal factor 3.040; and the maximum angular deformation is 0.00°. h and k are shown too, and depend on the coordinates rather than on the map. Measuring distortion

Tissot's indicatrix

A tiny circle on the sphere becomes an ellipse on the map, and the ellipse's two axes are the whole story. Almost every published indicatrix is drawn without them, which discards the content and keeps the decoration.

 A false origin moves every number and no geometry. British National Grid, drawn twice at the same scale. On the left the coordinates are measured from the projection's own origin, where the central meridian meets the true origin's parallel, and 52% of the country takes a negative easting or northing — the worst reaching -323 km. On the right the authority's published false origin of 400 km east and -100 km north has been applied, and none of them does. Every distance computed from the two sets of numbers agrees to the last bit a double has left after carrying six figures, and every bearing agrees exactly; the only thing that changed is that no coordinate carries a sign. The margins say the origin was fitted to the land rather than placed arbitrarily below it: 77 km spare in the west and 50 m in the south, on a grid 988 km tall. Grids, and what a survey does

A grid has an origin that is not there

Every published national grid measures from a point in the sea. The British one sits 400 kilometres west and 100 kilometres south of where the projection's own zero is, and moving it changes every coordinate in the country and not one distance or bearing.

 The indicatrix at 30°, 40° on Gall–Peters. A circle on the sphere becomes this ellipse on the map. Its semi-axes are the two principal scale factors, a = 1.083 and b = 0.923; their product is the areal factor 1.000; and the maximum angular deformation is 9.16°. h and k are shown too, and depend on the coordinates rather than on the map. Measuring distortion

What survives a change of coordinates

The scale along the meridian is a property of the map and the grid together. The principal scale factors are a property of the map alone. Only the second kind describes the projection, and the two are routinely quoted as though they were the same thing.

 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 "1:136,495" means at each latitude, at zoom 12. A screen map at zoom 12 prints one scale for the whole world. The curve is how much larger in scale the map really is, measured from the projection's own derivatives rather than from a formula: at 60° it is 1.98 times, so the map labelled 1:136,495 is a 1:68,765 map. The hollow marks are sec φ, the textbook answer. They do not sit on the curve — the worst gap is 9949 parts per million at 85° — because Web Mercator puts a geodetic latitude into a spherical formula, and the same spherical Mercator measured the same way reproduces sec φ exactly. What a machine does with it

The scale of a screen map is not one number

A zoom level prints one scale for the whole world, and the map is at that scale along exactly one line. At 60° north the picture labelled 1:136,495 is a 1:68,765 map — and the factor is not quite sec φ either, because the projection puts a geodetic latitude into a spherical formula.

 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. Grids, and what a survey does

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.

 London to Tokyo, seen four ways. The same two routes on four projections. The gnomonic projection renders every great circle as an exactly straight line, which is what it is for; Mercator renders every rhumb line straight instead. Neither path changed — only the map did. The families

The aspect is a free choice

A projection's distortion pattern is fixed relative to its own axis, and where that axis points is entirely up to the cartographer. Rotating it is the cheapest available improvement and it is the one most often left unmade.

 Angular deformation against latitude, four projections. The same quantity for Mercator, Gall–Peters, Mollweide, Winkel tripel, from the equator to 80°. A conformal projection sits flat on zero in the first of these plots and runs away in the second; an equal-area projection does the reverse. Nothing is flat in both. Measuring distortion

The two ways a map is wrong

Angles and areas fail independently. A projection can be perfect about one and catastrophic about the other, and a single number summarising "how distorted" a map is has already thrown away the distinction that matters.

 London to Tokyo on Orthographic. Two routes. The great circle is 9559 km and is the shortest path on the sphere. The rhumb line holds a single compass bearing the whole way and is 11296 km — 1737 km further, or 18.2 per cent. On Orthographic the rhumb line departs from straight by 2.9e-1 of its own length. Paths and directions

The great-circle vertex

One quantity is constant along a shortest path on a sphere, and it fixes the highest latitude that path will reach before the journey starts. That number is why polar routes exist, and it can be read off the departure bearing without tracing the route at all.

 One pixel at zoom 11: 76.4 m projected, 47.6 m at 51.5°. The grid squares are pixels, at their own size. The open mark is the stored coordinate and the filled one is where it is drawn: rounding moves it 20.2 metres, against a worst case of 33.6 — half a pixel's diagonal, which is the whole of the bound. The second point sits 43 metres away, 0.58 of a pixel east and 0.70 north, so whether the two are drawn as one dot or two is decided by where the tile grid happens to fall: they merge for 12 per cent of the possible offsets, against the 12 per cent the two fractions predict. What a machine does with it

The pixel is a place with a size

Drawing a coordinate rounds it to a pixel, which moves it by up to half a diagonal — 16.8 metres at zoom 12. Whether two points 43 metres apart appear as two dots is not a property of the data at all: they merge for 12 per cent of the positions the tile grid could take, and the closed form predicts 12.4.

 The angles of a triangle, and the curvature inside it. A triangle on the sphere with sides that are great-circle arcs. Its three angles sum to 212.26°, overshooting the flat 180° by 32.26°. Integrating the curvature over the interior gives 0.56306 against an excess of 0.56306 — the same number by two routes that share no arithmetic. That is Gauss–Bonnet, and it is how a flatlander measures the curvature of a world it cannot step outside. The impossibility

Measuring curvature from inside

A flatlander with a ruler can find out that its world is round, by drawing a triangle and adding the angles. Gauss–Bonnet turns that into an exact statement, and the total curvature of a closed surface turns out to be a number obtained by counting.

 The least distortion possible over a 20° region. Scale factor along a radius of the cap, each projection normalised to unit scale at the centre. Chebyshev's criterion names the projection whose scale is constant on the boundary as the conformal map of least scale variation, and for a cap that is the stereographic projection centred on it — reaching exactly sec²(ρ/2) = 1.0311 at the rim, marked. Every other conformal projection drawn here rises past that line before it gets there. It shows an optimum rather than a comparison. Measuring distortion

Scale distortion is the third failure

A map's stated scale is its scale at one place. Everywhere else it differs, by a factor that no conformal projection can avoid and that neither of the two usual distortion measures reports. It is the failure everybody uses and nobody counts.

 Two feet, and the one that is not a rounding. The US survey foot is 1200/3937 m and the international foot is exactly 0.3048 m, so the first is longer by 2.0000 parts per million. Read a coordinate in the wrong one and the error is proportional to the coordinate, not to any distance in the job: at 2.5 million feet from the false origin it reaches 1.52 m, and it crosses a survey's own closing tolerance of twenty millimetres at 33 thousand feet. Across a 300-foot site the same mistake changes every dimension by 0.18 mm — below every tolerance there is, which is why a plan in the wrong foot passes every check and is in the wrong place. Grids, and what a survey does

The units are part of the coordinate

Two feet were in use in the United States until 2022, differing by exactly two parts per million. A plan drawn in the wrong one has every dimension right to a fifth of a millimetre and sits three-quarters of a metre from where it belongs, which is why it passes every check anybody runs.

 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. What is taught wrongly

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.

 The same six projections, ranked over Europe and Chile. Each column orders the projections by Kavrayskiy's regional criterion — the root-mean-square departure of the two principal scales from unity, integrated over the region with the area element. The lines cross, which is the point: Lambert conformal conic leads over Europe and comes fifth over Chile. A table of projections ordered by distortion is a table about somebody's region. Measuring distortion

Distortion over a region

An indicatrix describes a point and every practical question is about a country. Going from one to the other means integrating, integrating means choosing a weighting, and the weighting is the step that turns a measurement into somebody's opinion.

 Both zones are right, and one of each is not. A 20 km baseline straddling the boundary between UTM zones 31 and 32 at 52°N, computed from grid coordinates three ways, on logarithmic bars. Taking both ends in zone 31 and taking both in zone 32 give answers 0 nanometres apart — the resolution a double has left after carrying a six-figure easting rather than a disagreement — so a job may use either and the overlap belt every zone publishes exists to let it. Taking each end from the zone it nominally belongs to gives 392 km, because the two eastings are measured from meridians six degrees apart and subtracting them measures nothing at all. The same ground point is 706 km east in one zone and 294 km east in the other, and both coordinates are correct. Grids, and what a survey does

Where two zones meet

A twenty-kilometre baseline across a zone boundary is computed correctly in either zone — the two answers differ by nanometres — and taking one end from each returns 392 kilometres. Both zones are right and the seam between them is not a small error but a category one.

 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. What a machine does with it

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.

 Which way Sinusoidal stretches the ground. The major axis of Tissot's indicatrix at 180 points, drawn as a stroke rather than an ellipse because the argument is the orientation. Longer strokes mark more elongated indicatrices. The direction departs from the graticule by up to 44.4°, so the scale "along the meridian" is not the scale along the direction that is actually stretched most. Measuring distortion

Distortion has a direction

Tissot's indicatrix is an ellipse, an ellipse has an orientation, and the orientation is never reported. On most projections it is not along the meridian, on a conformal one it does not exist at all, and both facts are computable from the same four derivatives as everything else.

 The area of one 20° × 10° cell at 50–60° north, seven ways. The cell has an exact area — R²Δλ(sin φ₂ − sin φ₁), 1,416,580 square kilometres — so every other row is a measurement of the method rather than of the ground. The equal-area projection returns it to 1.000000 and the spherical polygon formula to 1.000000, which is three routes agreeing — and the same cell integrated on the ELLIPSOID comes out 0.45 per cent away from all three, because the sphere is a model. Taking the shoelace in Mercator gives 3.06 times too much, and treating degrees as a length gives 1.75 times — about sec φ at the cell's middle, which is where that error comes from. What a machine does with it

Computing an area needs a surface

A shoelace over a ring of coordinates returns a number whatever the coordinates are. For one twenty-by-ten-degree cell it returns 1.75 times the true area in degrees, 3.06 in a conformal plane, and exactly the closed form in an equal-area one — and the closed form itself is 0.45 per cent out, because the sphere is a model too.

 A closed traverse cannot see a scale error. The same closed figure with two different errors in it. On the left every leg is 400 parts per million too long, which is roughly what forgetting the grid scale factor costs — and the figure closes to 2.3e-13 m, which is the last bit of a double rather than a measurement. Scaling every leg of a closed figure by the same factor produces a similar figure and a similar figure is still closed, so the check every specification leans on is blind to it. Every dimension in that figure is wrong: the perimeter is out by 1.17 m. On the right one angle is 20 seconds out — a far smaller disturbance in its own units — and the figure fails to close by 0.062 m. Closure tests the shape and says nothing about the size. Grids, and what a survey does

A traverse must close

The check every survey specification leans on cannot see the error this whole field is about. Scale every leg of a closed figure by four hundred parts per million and it closes to the last bit of a double, while every dimension in it is wrong.

 Three surfaces with the same curvature, one of which is a sphere. Every one of these is a surface of revolution whose Gaussian curvature is 1 at every point, built by solving r″ + r = 0 for the meridian rather than by writing a shape down. The spindle closes to a point with an angle deficit, the sphere closes smoothly, and the bulge does not close at all — it ends in two circular edges. A surveyor confined to a patch of any of them, measuring angles and distances, cannot tell which one it is. The impossibility

Two surfaces with the same curvature

The ladder's base says curvature is the obstruction to a faithful map. It has never asked whether curvature is the whole obstruction — and it is, locally: there is a whole family of surfaces with the unit sphere's curvature at every point, none of them a sphere, and a geodesic circle drawn on one agrees with the same circle on another to one part in 10¹⁴.

 A great circle and the ruled line, on Mercator. The shorter of the two routes is the curved one. The great circle between the two marked points is drawn against the straight line a ruler would give between them on Mercator; over 70° of arc the curve departs from the ruled line by 11.9 per cent of the chord. The flexion at the midpoint, which is the rate the image turns per radian of arc, is 0.788; the angular deformation there is 0.00°. Measuring distortion

Tissot stops at the first derivative

Every quantity this site has measured is read off one derivative of the projection. A map can be conformal at a point — the indicatrix a circle, the angular deformation zero to eleven figures — and still bend every geodesic through it, at a rate of 0.839 radians of turning per radian of arc.

 Why a levelled height is not a distance. The correction between raw levelling and orthometric height, for lines at 200, 500, 1000, 2000 metres above the geoid running north from 50°. It is not instrument error: level surfaces converge towards the pole, so a run that stays on one of them gains height relative to another. A line 2000 metres up reaches 647 millimetres over 400 kilometres. The dashed line is 10 millimetres, which is about what a first-order levelling network closes to over that distance — so this is not a refinement, it is the larger of the two numbers. What the numbers refer to

A levelled height is not a distance

Level surfaces converge towards the poles by five metres in a thousand, so a chain of perfectly executed levelling observations does not sum to a height difference. The correction over four hundred kilometres of northing is larger than the network's own closure.

 Six projections, at 40° north. Flexion (solid) and skewness (light) against direction of travel, drawn about a zero circle at one point of each projection. The three conformal ones have curves of exactly equal size, offset by exactly 90°, because on a conformal map both quantities are components of one vector — the gradient of the log of the scale. The others have neither property: ratios run from 0.72 to 2.72. Measuring distortion

Bending and stretching are one failure

The ladder was written expecting flexion and skewness to be two independent ways for a map to be wrong. On a conformal projection they are not independent at all: the two extremes are the same size to four figures and sit exactly ninety degrees apart, because both are components of a single vector.

All threads