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

Indicatrix — where it appears

The ellipse a projection maps an infinitesimal circle to, whose semi-axes are the two principal scale factors. Everything a distortion measure reports is a property of it, and it exists only where the projection has a derivative — which excludes a polyhedral net's corners.

Named by 7 essays across 2 fields — each of them below, with the objects they name alongside it.

The indicatrix is a limit, and here it is being taken. The axis ratio of the image of a circle of radius ρ, measured from the picture rather than computed from the derivatives, against the ratio the derivatives give. It converges, and it converges as ρ² — each halving of the radius quarters the departure. At the smallest radius drawn the two agree to 7.6e-4 relative. This is what "the indicatrix is a limit" means when it is working: there is a number for the ratio to settle on, and it settles.

The indicatrix at a point that has none

Ten essays on this ladder take the derivative for granted. Rung three proved the indicatrix is a limit and every figure since has assumed the limit exists. There are points on the maps this collection draws where it does not, and they are not exotic: the horizon of a gnomonic map, and every corner of every polyhedral net.

distortion · Tissot Tissot's ellipse and the one that governs gradients, at 20°E 48°N. The solid ellipse is the image of a small circle — Tissot's indicatrix, semi-axes a and b. The dashed one is the image of a unit gradient, whose semi-axes are 1/b and 1/a because a gradient transforms by the inverse transpose of the Jacobian rather than by the Jacobian. Its long axis therefore lies where the indicatrix's short one does. On Mercator both are circles, so a gradient's direction survives; on Lambert cylindrical the two ellipses are the same shape turned through a right angle, so the worst direction for a gradient is the best direction for a shape.

A slope is not a shape

Every map in this collection has carried geometry. An applied map far more often carries a field — elevation, pressure, a density — and the first thing anybody does with one is differentiate it. A gradient is a covector, it transforms by the inverse transpose of the Jacobian, and the ellipse that governs it is the indicatrix turned inside out.

distortion · Gradient A geodesic circle on Mollweide, and the two models of where it goes. A circle of geodesic radius ρ about 20°E 45°N, projected. The dashed outline is where the indicatrix says it goes — the ellipse a Tissot figure draws — and the thin solid one adds the quadratic term, which is the flexion and skewness this ladder measures. At 16° the indicatrix is out by 13.5 per cent of the figure's own size and the second-order model by 2.20 — a factor of 6.16, which is what one more derivative buys.

The size at which the second derivative arrives

Six essays have measured a projection's second derivative at points, over regions and under transformations, and none of them says at what size it stops being a curiosity. The answer needs a figure with an extent rather than a point, and it is smaller than anybody drawing a national map would guess: the indicatrix alone places a shape to one part in a thousand out to thirteen kilometres.

distortion · Flexion 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.

An average of ellipses is not an ellipse

Rung three integrates distortion over a region and finds the weighting is somebody's opinion. It never asked what was being averaged. An indicatrix is a positive-definite matrix, matrices form a cone rather than a vector space, and the arithmetic average of an equal-area map's indicatrices comes back inflating area by up to 11 per cent.

distortion · Tissot The same tolerance, applied in two orders, at 65°. The faint line is the region's boundary as built, 1025 vertices across 400 km of ground. Both pipelines were given the same tolerance of 2000 m on the ground. Simplifying in degrees and then projecting keeps 311 of them; projecting into Mercator and then simplifying keeps 129; doing it on the ground itself, which no pipeline does, keeps 129. The two drawn lines separate by 1883 m, which is 94 per cent of the tolerance that was supposed to bound the whole operation.

Simplification does not commute with the projection

A pipeline either simplifies the geometry and then projects it, or projects it and then simplifies. Both orders are in use, neither is recorded, and given the same tolerance in ground metres they keep different vertices — 129 of them on the ground, 367 in degree space at 80°, and 459 on an equal-area page.

applied · Generalise A 900 km circular accuracy at 55° north, projected. Six thousand ground positions drawn from a circular error of 900 kilometres about one place, each projected in Mercator and plotted as a displacement from the projected place. The curve is the nominal 95 per cent ellipse, computed the standard way — the ground covariance sandwiched between the projection's own derivatives. It holds 93.83 per cent of the points, the cloud is measurably longer than it along its own long axis by 5.52 per cent, and it is not symmetric: the third moment along the page's second axis is 0.727 rather than zero.

The error ellipse is not an ellipse

Rung two pushed a covariance through a projection with the same matrix sandwich that draws an indicatrix. That is a first-order operation on a map with a second derivative, so the propagated distribution is not the ellipse the sandwich draws — and a nominal 95 per cent ellipse holds 93.06 per cent on one projection and 95.63 on another, in opposite directions, from the same input.

distortion · Precision How aligned a region's ellipses are, for every projection and every region. The resultant length of the doubled indicatrix orientations, over 12 projections and 9 regions. One means every ellipse in the region points the same way; zero means they are spread evenly and cancel. Three cylindrical rows are 1.00 throughout, and everything else varies down the row and across it — which is the answer to the question the number was first stated without: alignment is a property of the pair. Mercator's row is the exception that is not a measurement, because a conformal projection's indicatrix is a circle and a circle has no orientation.

Whether the ellipses point the same way

The previous rung found that alignment decides whether the average of a region's deformations is above or below the deformation of its average, and stated it at ten projections over one region. Swept over a hundred and eight pairs, the answer is that alignment belongs to the projection 28 per cent, to the region 38, and to neither 33 — and two rows of the table turn out not to be measurements at all.

distortion · Tissot

Named alongside it

The objects these essays reach for when they reach for this one.

ConformalityAnisotropyEqual-areaPrincipal directionAngular deformationJacobianLimitProjectionScale factorSecond-orderSkewnessTaylor expansion

All concepts