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

How Mercator distorts, by latitude

Drawn here at the parameters it defaults to, with every essay that calls it.
How Mercator distorts, by latitude. Angular deformation and areal scale plotted against latitude along the meridian at 0°. On Mercator the angular deformation reaches 0.0° and the areal factor reaches 91.5.

Angular deformation and areal scale plotted against latitude along the meridian at 0°. On Mercator the angular deformation reaches 0.0° and the areal factor reaches 91.5.

It is drawn by distortion-figure with show: "distortion-profile" — one member of a family of 7 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.

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

 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.

 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, gallPeters, mollweide, winkelTripel, 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.

 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.

 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 six named here are Mercator, Web Mercator, Gall–Peters, Mollweide, Equirectangular, Stereographic. Measuring distortion

Measuring instead of naming

A projection is called conformal because that is its name. Running the definition as a computation over several hundred points takes about twenty lines, catches a projection the whole internet uses, and is almost never done.

 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. This is the only figure on the site showing 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.

 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.

 How Equirectangular distorts, by latitude. Angular deformation and areal scale plotted against latitude along the meridian at 0°. On Equirectangular the angular deformation reaches 114.2° and the areal factor reaches 11.5. What is taught wrongly

The plate carrée, the projection nobody chooses

Plotting latitude against longitude on ordinary axes is a projection. It preserves nothing, its angular deformation reaches 108° and its areal error eightfold, and it is probably the most widely produced map in the world because it is what happens when nobody decides anything.

 The image of a circle on Mercator, against its own indicatrix. A circle of three radii on the ground at 30°, 40°, projected exactly — the solid curve — against the ellipse the indicatrix predicts for it, dashed. The filled dot is the image of the circle's centre and the hollow one is the centre of area of what was actually drawn, which is not the same point. The departure runs from 3.14 per cent at 4° to 15.80 per cent at 16°, so it grows in proportion to the radius rather than to its square: halving the circle halves the relative error and does not quarter it. Measuring distortion

The indicatrix is a limit

Tissot's ellipse describes an infinitesimal circle, and every published one is drawn finite. The error of the description falls as the radius rather than as its square — halving the circle halves the lie — and on the Robinson projection at a table entry it does not fall at all: sixty metres and six kilometres are both 1.17 per cent wrong.

 The azimuthal family, as one function of one variable. Every azimuthal projection is a rule for how far from the centre to draw a point at angular distance ρ, and nothing else. five named projections are plotted as their radial functions — tan ρ, 2 tan(ρ/2), ρ, 2 sin(ρ/2), sin ρ — and the two dashed curves are the solutions of the two conditions, integrated from the centre with no more than "true scale at the centre" to start from. They land on the stereographic and the Lambert azimuthal to 8e-10. The names are the answers to the equations, not descriptions of viewpoints. The families

The azimuthal family is one function

Five azimuthal projections are taught as five viewpoints — from the centre, from the far pole, from infinity. They are one projection with five choices of how far out to draw a point at angular distance ρ, and each named property is a differential equation in that one function: f′ = f/sin ρ integrates to the stereographic projection and f f′ = sin ρ to Lambert's, from nothing but true scale at the centre.

 Conformality helps and does not save it. The proportion of a stated terrain whose plan curvature changes sign when it is read off a grid in each projection. A conformal map turns every direction through the same angle, so the contour and the slope line stay perpendicular and the sign ought to survive — and it mostly does, at 0.30 per cent against 13.3. It is not zero, and the term that flips it is the gradient of the scale factor: the curvature of a curve under a conformal map is (κ − ∂ₙ log λ)/λ, and Mercator's λ has a gradient. Measuring distortion

The curvature of a field is not the curvature of its picture

Whether a place is a spur or a hollow is the sign of a second derivative, and every automatic terrain classification is built on it. Read off an equal-area grid, that sign is wrong on 13.25 per cent of a sheet — and conformality reduces it to 0.30 per cent without removing it.

The whole library