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

Tissot's ellipse, and the three things it does not capture

indicatrix-map scatters it across a whole map, indicatrix-detail magnifies one, indicatrix-compare puts the same point under several projections, departure-rate and mean-indicatrix reduce the field to a number. cell-shapes is the same question asked of a finite cell rather than an infinitesimal one, direction-field draws the direction the stretch points in, and finite-angle is the error the ellipse cannot see because it lives at finite size rather than at a point.

indicatrix-map scatters it across a whole map, indicatrix-detail magnifies one, indicatrix-compare puts the same point under several projections, departure-rate and mean-indicatrix reduce the field to a number. cell-shapes is the same question asked of a finite cell rather than an infinitesimal one, direction-field draws the direction the stretch points in, and finite-angle is the error the ellipse cannot see because it lives at finite size rather than at a point.

Every one of these is the same generator answering a different question, which is why they share a file, a set of colour roles and a set of assertions. Changing it changes all 8.

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

What it draws

One picture per question the family answers. Each has a page of its own.

Tissot's indicatrix across Mercator. A small circle on the sphere, drawn where the projection puts it. The dashed circle behind each is what an undistorted map would show. On Mercator every ellipse is a circle — ω never exceeds 1.0e-6°, and the areal factor reaches 4.0.

indicatrix-map

Tissot's indicatrix across Mercator 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.

indicatrix-detail

The indicatrix at 45°, 55° on Mercator The same point at 30°, 55° under four projections. One circle on the sphere, four projections, four ellipses. A conformal projection keeps it circular and lets the area run; an equal-area projection keeps the area and lets the shape go. Nothing keeps both, and the dashed circle shows what keeping both would look like.

indicatrix-compare

The same point at 30°, 55° under four projections How large a circle each indicatrix describes, at 30°, 40°. The radius at which the image of a circle departs from the ellipse Tissot's construction predicts by 1 per cent, measured by bisection at 30°, 40° on six projections. Albers equal-area conic keeps its ellipse honest to 829 kilometres and Gnomonic to 46 — a factor of 17.8, decided by how fast each projection's own scale factor changes rather than by how large its distortion is.

departure-rate

How large a circle each indicatrix describes, at 30°, 40° Four averages of one region's indicatrices — Robinson over the whole sphere. The faint ellipses are the indicatrix at ninety sampled places; the four drawn over them are four candidate averages of exactly that set. The arithmetic mean of the tensors gives semi-axes 1.1621 and 1.0182; the log-Euclidean mean gives 1.0053 and 0.9608; the Karcher mean gives 1.0042 and 0.9619; and taking the mean of a and the mean of b separately — which is what a published average distortion usually is — gives 1.2335 and 0.8200. Their maximum angular deformations are 7.57°, 2.60°, 2.47°, 23.23°, against a mean of the pointwise deformations of 20.67°. Five numbers, one set of ellipses.

mean-indicatrix

Four averages of one region's indicatrices — Robinson over the whole sphere One 20° × 10° cell at 60°, under four projections. The same patch of the sphere, drawn by four projections and scaled to fit. The shapes differ and so do the areas: the figure under each panel is the areal inflation relative to the same cell at the equator, so 1.00 means the projection treats the two fairly.

cell-shapes

One 20° × 10° cell at 60°, under four projections 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.

direction-field

Which way Sinusoidal stretches the ground A triangle of 25.0° excess, drawn on Mercator. The three geodesic sides as curves and the three straight sides a ruler draws, with each vertex labelled by how far the ruler's angle is from the true one. The straight-sided triangle's angles sum to exactly 180° because it lies in a plane, and the real one's sum to 205.02°, so the three errors have to account for the whole 25.02° of spherical excess between them — and they do so on this conformal projection exactly as they do on any other.

finite-angle

A triangle of 25.0° excess, drawn on Mercator

Where it is called

Changing this changes every one of these figures.

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.

 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 machinery here found that without being told to look.

 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.

 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.

 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.

 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 five named here are Mercator, Stereographic, Gall–Peters, Mollweide, Winkel tripel. The impossibility

The trade-off is two lines

Conformal means the two principal scales are equal. Equal-area means their product is one. Both at once forces both to one, which is an isometry, which the curvature forbids. That is the entire argument.

 A 1,000 km bar on Mercator, drawn at 0° and read elsewhere. The same length of paper, carried up the map. Each pair of bars is the ground distance that length actually spans at that latitude — the filled bar along the parallel, the outline along the meridian. At 75° the reading along the parallel is 259 km against the 1,000 the bar claims, an error of 74 per cent. The two readings agree everywhere, because Mercator is conformal — so one number per latitude corrects any measurement taken off it. What a machine does with it

A scale bar is right in one place

The bar in the corner of a world map is a picture of a distance, and it is a true picture along one line. On Mercator it reads 500 kilometres for a thousand at 60° north — and on an equal-area map it reads 500 one way and 2,000 the other, so the projection recommended for measuring is the one on which no single correction exists.

 One 20° × 10° cell at 60°, under four projections. The same patch of the sphere, drawn by four projections and scaled to fit. The shapes differ and so do the areas: the figure under each panel is the areal inflation relative to the same cell at the equator, so 1.00 means the projection treats the two fairly. What is taught wrongly

The projection that shows true size

There is no such thing, and the phrase hides a real question. Equal-area projections preserve area and destroy shape; nothing preserves size in the sense the phrase implies; and the tools that make the point best are not maps at all.

The whole library