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

The thread: Computed, not quoted — page 1

Scale factors, areal ratios, angular deformations, geodesic distances and the areas of countries are all computed from the projection's own formulae. None is a number recalled from a table. Essays 1 to 24 of 333.
London to Tokyo on Mercator. 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 Mercator the rhumb line departs from straight by 5.0e-9 of its own length. Paths and directions

The shortest route is not straight

The shortest path between two points on a sphere is an arc of a great circle, and on almost every map it is a curve. The straight line on a Mercator chart is a different route entirely, and on some journeys it is twenty-eight per cent longer.

 The tile pyramid, four levels down to quadkey 120. The world as a square, quartered three times. Each level's tile is exactly half the width of its parent, so a tile is four tiles at the next level and never needs resampling to serve one — the property the whole scheme rests on, and one that holds only because the projected world is square. Level 3 has 64 tiles at 19567.9 metres per pixel, which at 51.5° north is 12181.3 metres of ground per pixel rather than the number the scheme publishes. What a machine does with it

A screen map is a pyramid of tiles

The scheme every slippy map runs on is a coordinate system with three integers and one projection, and almost all of it is forced. A square world is what makes the quadtree work, the levels are exact powers of two, and the published resolution — 156,543 metres per pixel at zoom zero — is a distance on the ground at exactly one latitude.

 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.

 Six surfaces and their Gaussian curvature. Curvature computed at the centre of each surface, by both available routes — from the way the surface sits in space, and from distances measured inside it alone. The two agree, which is Gauss's theorem. A surface with K = 0 can be unrolled flat without stretching: here the plane, the cylinder and the cone do, and the sphere, the torus and the pseudosphere do not. The impossibility

No map is faithful

Not "no map yet", and not "no map at page size". Gaussian curvature can be computed from inside a surface, a sphere has some and a plane has none, and that closes the question permanently.

 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.

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

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 test that found it was not looking for it.

 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.

 A cylinder unrolls exactly; a sphere does not. Both pictures show the same grid. On the left it is wrapped round a cylinder of radius 1, on the right it is laid flat, and every distance in the grid is the same in both — the circumference is 2π and so is the width of the rectangle, checked to 10⁻⁹. This is possible because a cylinder has zero Gaussian curvature. No corresponding picture exists for a sphere. The impossibility

What can be unrolled

A cylinder is obviously curved and is intrinsically flat, so it lays out on a table with nothing stretched. A sphere is not. The distinction is exactly zero Gaussian curvature, and it is why paper tubes exist and paper globes do not.

 Angular deformation around three standard parallels. Three equal-area cylindrical projections differing only in where they are exact — standard parallels at equator, 30°, 45°. Each has zero angular deformation at its own standard parallel and grows away from it in both directions. Choosing a standard parallel is choosing which latitudes to treat well, and there is no choice that treats them all well. The families

What a standard parallel buys

A standard parallel is a line where the projection is exact. Choosing one does not reduce the distortion — it decides where the distortion is zero and lets everything grow away from it.

 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.

 New York to Madrid on Mercator. Two routes. The great circle is 5768 km and is the shortest path on the sphere. The rhumb line holds a single compass bearing the whole way and is 5939 km — 171 km further, or 3.0 per cent. On Mercator the rhumb line departs from straight by 9.6e-16 of its own length. Paths and directions

Why Mercator exists

A ship can hold a compass bearing and cannot easily hold a great circle. Mercator is the answer to one question — what must a map do so that a constant bearing is a straight line — and it answers it exactly.

 How much each projection inflates a cell, by latitude. Five patches of the sphere, each 20° by 10°, and the factor by which each projection enlarges them relative to the equatorial one. An equal-area projection sits flat on 1. Mercator reaches 15.4× at 70°, which is the mechanism behind every complaint about the size of Greenland. What is taught wrongly

Mercator against Peters

The most-argued question in cartography, conducted almost entirely without anyone measuring anything. Both projections are exactly what they claim, each destroys what the other keeps, and the numbers are computable in either direction.

 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.

 6 projections of the same sphere. The same graticule under Robinson, Winkel tripel, Mollweide, Mercator, Gall–Peters, Eckert IV. Every one of them is a faithful drawing of the same object and no two agree, because each has chosen a different thing to preserve. What each projection optimises

Compromise projections

A projection that preserves nothing exactly can distort everything less than one that preserves something exactly. For a general-purpose world map that is the right trade, and it is why the two most widely used ones today have no exact property at all.

 A scale factor is a property of a point, and a distance is not. The grid scale factor sampled along a 362 km line on the British National Grid at a bearing of 90°, with the three rules in use drawn as the constants they replace it by. In a transverse Mercator k depends on the easting almost alone and goes as its square, so along an east–west line the curve is very nearly a parabola. The endpoint mean is the chord across it and lands 268.6 ppm high, which is 97.3 m over this line. Simpson's rule integrates a parabola exactly and lands 0.004 ppm out — 1.3 mm. Neither number is about the line's length. Grids, and what a survey does

The scale factor of a line

A scale factor is a property of a point and a distance is not measured at a point. The three rules in use for averaging one along a line differ by a hundred metres on a 362-kilometre baseline, and which of them is adequate is decided by the line's direction rather than its length.

 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. Paths and directions

The gnomonic companion

One projection turns every great circle into a straight line, and it is the only one that does. It shows less than half the sphere, distorts enormously, and was indispensable for three centuries because of that single exact property.

 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.

 The scale factor was chosen, and the choice is a V. The worst departure of the grid scale factor from unity across a region 9° wide and 9° tall, for every scale factor between 0.998108 and 1.0002, on a central meridian of 2°W. A tangent projection touches at one meridian and is too big everywhere else, at 1987 ppm. Scaling the whole grid down slides that interval until it straddles unity, and the minimum sits at 1/√k_max = 0.999008102, where the worst departure is 993 ppm — a factor of 2.00, which is the most this construction can buy and is reached exactly. The published value marked beside it was chosen the same way, for this region, in the 1930s. Grids, and what a survey does

The scale factor was chosen

A grid's scale factor is not a constant of nature. Britain's 0.9996012717 is the reciprocal square root of the worst distortion a tangent projection would have had over the country, it halves the worst-case error exactly, and the halving is the most that construction can ever buy.

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

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.

 The cut at 85.0511°, and the alternatives that are not square. Mercator's northing runs to infinity at the pole, so a tiling has to stop somewhere, and the latitude is not a rounding: 85.0511° is where the northing equals half the world's width, which is the only cut that makes the projected world a square. The bars are what other cuts would give — at 89° the world is 1.51 times as tall as it is wide, and a single square root tile cannot cover it. The price is 0.373 per cent of the Earth's surface, 1,901,487 square kilometres in two caps, computed from 2πR²(1 − sin φ) rather than estimated. What a machine does with it

The square costs the poles

The cut at 85.0511287798° is the most-quoted number in web mapping and is almost never derived. It is where Mercator's northing equals half the world's width — the condition for a square — and it drops 1,901,487 square kilometres. A two-tile root would have reached 89.786° and dropped 3,558.

 Transverse Mercator. The graticule of the Transverse Mercator projection at 30° of longitude and 15° of latitude. the only projection on which a constant compass bearing is a straight line. It is conformal. The families

Transverse Mercator and the series that computes it

The projection most of the world's survey data lives in has no closed form. What every national grid actually computes is a truncated power series in the flattening, and how far it can be trusted is an engineering parameter rather than a property of the projection.

All threads