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.

“How distorted is this map” sounds like a question with a number for an answer. It is not, and the reason is that there are at least two independent things going wrong and no principled way to add them.

Angular deformation against latitude, four projectionsThe 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.20°40°60°80°MercatorGall–PetersMollweideWinkellatitudeangular deformationalong a meridian
Fig. 1 Angular deformation against latitude for four projections. Mercator lies flat on zero — it is exactly right about angles. The equal-area projections rise steeply. Read this figure alone and Mercator is the clear winner.
How much each projection inflates a cell, by latitudeFive 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.equator23°45°60°70°MercatorEquirectangularMillerGall–Peterslatitude of the cell1× means treated fairlyrelative to the equator
Fig. 2 The same four projections measured the other way. Now the equal-area ones are flat and Mercator runs away to fifteenfold. The ranking has reversed completely, and neither figure is wrong.

Two measurements, two rankings, opposite conclusions. That is not a paradox; it is what independence means.

What the two quantities are

From the indicatrix come the principal scale factors aa and bb at each point — the largest and smallest amounts by which the map stretches.

Angular deformation ω\omega depends on their ratio:

sinω2=aba+b\sin\frac{\omega}{2} = \frac{a-b}{a+b}

It is zero when a=ba = b, whatever their common value. A projection can double every distance at a point and still have zero angular deformation there, because doubling everything preserves every angle.

Areal error depends on their product, abab, and is zero when the product is one. A projection can have a=4a = 4 and b=1/4b = 1/4 and be perfectly equal-area at that point while destroying every shape.

Ratio and product are independent — knowing one says nothing about the other — and that is the whole of the independence. There is no hidden relationship, no conservation law linking them, and no way to trade a known amount of one for a known amount of the other.

What is genuinely forced

Only one thing: both cannot be zero at once over a region, because that would make the map an isometry.

Everything else is free. A projection can be:

  • zero angular error and unbounded areal error (Mercator)
  • zero areal error and large angular error (Gall–Peters)
  • moderate amounts of both (Winkel tripel, Robinson)
  • large amounts of both (equirectangular, which measures 108° and eightfold)

That last case is worth dwelling on because it disposes of an intuition. The trade-off does not mean that being bad at one buys being good at the other — a projection can be bad at both, and the plate carrée is. Being good at something requires choosing it.

Every projection in the library, measured against both propertiesMaximum 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.10⁻⁶10⁻⁴10⁻²110010⁻¹¹10⁻⁸10⁻⁵10⁻²1010000maximum angular deformation / degreesmaximum areal errornothing here, everMercatorStereographicWeb MercatorGall–PetersMollweideEquirectangulartolerances: 0.0001° and 0.00000121 projections measured
Fig. 3 The two quantities plotted against each other for every projection in the library. Points near the left edge are conformal; points near the bottom are equal-area; the shaded corner is empty by theorem. Everything else is scattered, because the two coordinates are independent.

That scatter is the essay’s argument in one picture. If the quantities were related, the points would lie on a curve. They do not.

Why a single number fails

Attempts to summarise a projection’s distortion in one figure exist and are used, and their limitations are structural rather than fixable.

The obvious approach is to combine ω\omega and the areal error into a weighted sum, integrate it over the map, and compare. Airy’s criterion and its descendants do this. Two problems follow immediately:

The weighting is a choice. How many degrees of angular deformation equal a doubling of area? There is no principled exchange rate, because they measure different things. Any weighting encodes a purpose, and once a purpose has been chosen the comparison is nearly settled anyway.

The region and the measure are choices too. Integrating over the whole sphere weights the empty ocean equally with everything else. Weighting by land area, or by population, gives different rankings — and each is defensible for a different map.

So a single number is a purpose in disguise. That is not a reason never to use one; it is a reason to state the weighting, which almost nobody does.

Where each failure shows

Because the two are independent, they are visible in different ways and matter to different users.

Angular deformation shows locally. A shape looks wrong — a country looks stretched, a circular lake looks elliptical, streets meet at wrong angles. It is noticeable at any scale and is what makes a map look distorted.

Areal error shows globally. No individual shape looks wrong; the relative sizes are wrong. This is much harder to see and much easier to be misled by, because there is nothing locally odd to notice. The Greenland effect is invisible in any single region and enormous across the map.

That asymmetry has a consequence worth naming: areal distortion is the more dangerous of the two, precisely because it does not look like anything. A reader can see that Gall–Peters has stretched Africa. A reader cannot see that Mercator has enlarged Greenland fifteenfold, because Greenland’s own shape is fine.

One 20° × 10° cell at 60°, under four projectionsThe 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.Mercator5.62× areaGall–Peters1.00× areaMollweide1.00× areaEquirectangular2.36× areathe cell at 60°scaled to fit, areas as stated
Fig. 4 The same 20°-by-10° patch of sphere at 60–70° north, under four projections, scaled to fit. The shapes differ visibly. The areal inflation is printed underneath because the picture cannot carry it — each panel has been scaled, which is exactly the operation that hides the areal error.

A third quantity, occasionally

Some treatments add a third failure: distance distortion, whether lengths are preserved.

It is not independent of the other two — it is bounded by them, since aa and bb are the extreme scale factors — but it is a different question, because a projection can preserve distances from a particular point or along particular lines without preserving them generally.

Azimuthal equidistant is the standard example: distances measured from its centre are exact, in every direction, and distances between any other pair of points are not. That is a genuine and useful property and it does not appear in either of this essay’s two measurements, which is a reminder that the two-quantity picture is a simplification chosen because those two are the ones that can hold everywhere at once.

Where the two failures are worst

The two distortions do not merely differ in kind; they differ in where on a map they concentrate, and the pattern is worth knowing.

Angular deformation for a cylindrical equal-area projection is worst at the equator and at the poles, and zero at the standard parallel. For a pseudocylindrical projection it is worst at the outer edges, where the meridians are most oblique.

Areal error for a conformal projection is worst wherever the scale factor is largest — the poles for Mercator, the antipode for stereographic — and grows without bound.

So the two failures are usually in different places, which is one more reason a single summary number is misleading: it averages a defect at the edges with a defect in the middle.

Tissot's indicatrix across MollweideA small circle on the sphere, drawn where the projection puts it. The dashed circle behind each is what an undistorted map would show. On Mollweide ω reaches 75°, and the areal factor reaches 1.0.dashed: an undistorted circledrawn in Mollweide
Fig. 5 Mollweide’s indicatrices. All the same area, and the shape distortion is mild in the centre and severe at the outer edges — a spatial pattern quite different from Mercator’s, where the shape is perfect everywhere and the areas fail toward the poles.

Two failures, two audiences

A useful way to decide which matters: ask who is reading the map and what they will do with it.

Someone recognising a shape — finding a country, matching a coastline, identifying a region — is affected by angular deformation and essentially unaffected by areal error. A map with perfect shapes and wildly wrong areas serves them well.

Someone comparing quantities — reading a choropleth, judging relative extent, estimating how much of something there is — is affected by areal error and can tolerate a good deal of shape distortion.

Someone measuring — a surveyor, a navigator, an engineer — needs the property their instrument assumes, which is usually conformality because instruments measure angles.

Those three requirements are genuinely independent, and no map serves all of them. Which is why the choice needs a purpose and why the two-quantity framing is the useful one: it maps onto three real user needs rather than onto one abstract notion of accuracy.

The correlation that does exist

One qualified exception to the independence, because it is worth being precise.

Over a fixed family of projections with a tuning parameter, the two errors often do trade against each other, because the family has already fixed most of the freedom. Sliding a standard parallel changes both, and in opposite directions over part of the range.

But that is a fact about the family rather than about the quantities. Across the whole space of projections there is no relationship, which is what the scatter in the audit plot shows and what makes a projection that is bad at both possible.

How Robinson distorts, by latitudeAngular deformation and areal scale plotted against latitude along the meridian at 0°. On Robinson the angular deformation reaches 115.3° and the areal factor reaches 2.7.-60°-30°30°60°angular deformation, to 115°areal factor, to 2.7×latitudetwo independent distortionsalong the 0° meridian
Fig. 6 Robinson measured along a meridian. Both curves are moderate and neither is flat, which is what a compromise looks like when both quantities are plotted together — and neither curve predicts the other.

The third failure, named properly

Distance distortion is the quantity most often added to the list, and it behaves differently enough to be worth separating.

It is not a global property. No projection preserves all distances — that would be the isometry the theorem forbids — so every equidistance claim is qualified: distances from the centre, or along the meridians, or along a standard parallel.

That makes it a family of local guarantees rather than a single property, and it does not appear in the two-quantity framing because the two quantities are the ones that can hold everywhere at once.

It is nonetheless the property a great many maps actually want. A map answering “how far is it from here” needs distances from one point, and the azimuthal equidistant projection provides exactly that.

Azimuthal equidistantThe graticule of the Azimuthal equidistant projection at 30° of longitude and 15° of latitude. distances from the centre are true — from the centre, and from nowhere else. It is neither conformal nor equal-area.neither conformal nor equal-areadrawn in Azimuthal equidistant
Fig. 7 The azimuthal equidistant projection. Distances from the centre are exact in every direction; distances between any two other points are not. That is a genuine guarantee of a different shape from either quantity this essay measures.

What to report instead of one number

The practical recommendation, since the essay is against scalar summaries.

Report both quantities, at stated locations or as a range over a stated region. “Angular deformation up to 12° and areal error up to 1.4 over the mapped area” says everything a single index would and does not hide the trade.

That is what this site’s captions do, and it is why they are longer than the usual ones. A caption reading “moderate distortion” is a summary somebody computed with a weighting they did not state.

Angular deformation against latitude, four projectionsThe same quantity for albers, lambertConformalConic, winkelTripel, miller, 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.20°40°60°80°AlbersLambertWinkelMillerlatitudeangular deformationalong a meridian
Fig. 8 Four projections from different families and different property classes, measured on one axis. Reading this plot alone would rank them; reading it with the areal plot would rank them differently; and reporting both is the only honest summary.

A closing note on vocabulary. “Distortion” as a mass noun — more of it, less of it — is the habit this essay is arguing against, and it is deeply embedded in how the subject is discussed. There is no such quantity. There are two of them, they are independent, and any sentence using the singular has already chosen a weighting without saying so.

The independence also explains why the subject has so many projections. If distortion were one quantity there would be a best projection and the field would have found it. Because there are two, there is a frontier rather than an optimum, and every point on that frontier is somebody’s right answer.

The two-quantity framing is also what makes the site’s audit possible. Both are computable at every point from the same four derivatives, both have a well-defined zero, and both can be asserted against. A single blended index would have neither property — there is no value of it that means “correct”, so there would be nothing to assert.

A last note on how the two quantities are reported in practice. Software and documentation typically give the scale factor and stop, which is one number describing a conformal projection adequately and describing anything else not at all. Where a projection is not conformal there is no single scale factor, and a tool reporting one has picked a direction on the reader’s behalf.

The independence is finally what makes this subject a design discipline rather than an optimisation. If distortion were scalar there would be a best map and the work would be finding it. Because there are two quantities and no exchange rate between them, there is a frontier of non-dominated choices, and picking a point on it requires knowing what the map is for — which is a judgement, informed by measurement and not replaced by it.

What was computed here

Both quantities come from the same four partial derivatives, so a systematic error in the differentiation would move both together — which is a reason to check the derivatives independently rather than to trust their agreement. The analytic scale factors are implemented for the four projections that have simple ones and agree with the numerical route to better than 3×10113\times10^{-11}.

Areal errors here are computed two ways. The pointwise version comes from the Jacobian determinant. The cell version comes from comparing a projected polygon’s planar area against the closed form R2Δλ(sinφ2sinφ1)R^2\,\Delta\lambda\,(\sin\varphi_2 - \sin\varphi_1) for a latitude–longitude cell, which is exact. The cell formula is verified by summing a full covering of the sphere and requiring 4π4\pi, to 7×10157\times10^{-15}.

Using closed-form cells rather than a coastline dataset is deliberate. A polygon’s area depends on its simplification level, so a Greenland-against-Africa figure taken from a shapefile is partly a measurement of the vendor’s generalisation. The cells have no such ambiguity.

What the pictures cannot show

The independence, directly. Every figure here shows one quantity or the other, and the claim that they are unrelated is carried by the scatter in the audit plot and by the reversal between the first two figures. A picture of two independent quantities is two pictures.

The cell figures also hide the very thing they measure: each panel is scaled to fit, and scaling is what destroys the areal comparison. The numbers are printed because there is no honest way to draw four differently-sized shapes at true relative size and still see any of them.

Who found it, and when

Tissot’s 1881 analysis separates the two explicitly, and the distinction has been standard in the technical literature ever since.

It is much less standard outside it. Popular accounts routinely speak of “distortion” as a single quantity and rank projections along one axis, which is where the intuition that Mercator is simply worse than Gall–Peters comes from — a ranking that is correct on area, reversed on angle, and meaningless without saying which.

Attempts at scalar summaries began with Airy in 1861 and continue; Kavrayskiy, Jordan, Klingach and others all proposed measures, and no two agree on the ranking of the compromise projections. The disagreement is not a failure of any of them. It is what happens when several people choose different weightings for quantities that have no natural exchange rate.

Where this goes next

The tool both measurements come from is Tissot’s indicatrix. The one relationship between them that is forced is the trade-off is two lines. And for what a projection does with the freedom the independence leaves, every projection minimises something.