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

Gradient — where it appears

The rate and direction in which a field changes, whose horizontal part turns a geoid into a deflection and a potential into an attraction. The horizontal part of the gravity potential's gradient is the deflection of the vertical, and its integral down a plumb line is how far the line bends.

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

Six projections, at 40° north. Flexion (solid) and skewness (light) against direction of travel, drawn about a zero circle at one point of each projection. The three conformal ones have curves of exactly equal size, offset by exactly 90°, because on a conformal map both quantities are components of one vector — the gradient of the log of the scale. The others have neither property: ratios run from 0.72 to 2.72.

Bending and stretching are one failure

The ladder was written expecting flexion and skewness to be two independent ways for a map to be wrong. On a conformal projection they are not independent at all: the two extremes are the same size to four figures and sit exactly ninety degrees apart, because both are components of a single vector.

distortion · Flexion Ground north, on a map that has the pole on it. Every arrow points along its own meridian, towards the pole, and the pole is at the centre. Walking once anticlockwise round any loop enclosing it turns the arrow once anticlockwise as well: the index is 1, counted as a winding number with no distance anywhere in the calculation. A field like this cannot be combed flat. There is no way to choose a page direction for north at every point of the neighbourhood without the choice tearing somewhere, and the somewhere is the point in the middle.

North cannot be up everywhere

Ground north is a field of arrows on the sphere, and a field of arrows on a sphere must vanish somewhere. The failure is not measured, it is counted: the indices of the zeros sum to two, obtained here as a winding number in seven different charts with no distance anywhere in the calculation, and it is the same two that Gauss–Bonnet gets by integrating curvature.

impossibility · Topology 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 The same six contours on Mercator and Lambert cylindrical. The value of a harmonic sum with a summit and a basin in the northern mid-latitudes travels with the point, so the set of points at a stated level is the same set on every map and each contour is exactly right on both panels. Everything a reader measures from them is not: the spacing between neighbouring contours, their lengths, and the area between two of them all change from one panel to the other, and the two panels are the same field.

The contour is right and the reading is wrong

There is exactly one thing about a field that no projection can get wrong: which points share a value. The contour lines on any two maps of the same field are the same set of points. Every quantity a reader takes off them — the spacing, the length, the area between two of them, the hypsometric curve — is not, and the two kinds of map get different ones wrong.

distortion · Gradient A mass, the geoid it raises, and the plumb lines that lean towards it. A sphere of 5 km radius buried 8 km down, denser than its surroundings by 500 kg per cubic metre — a salt dome, or an ore body. Its mass is 261.8 × 10¹² kg, and outside itself its field is a point mass's exactly, so everything above is a closed form. The geoid rises 22 centimetres over it, drawn at 20,000× the true slope; the plumb line leans by at most 2.21 arcseconds, and it does so 5.7 km to the side rather than above the body, because the deflection is the geoid's SLOPE and a slope is zero at a summit. That is the number this site has been able to relate and unable to compute since its practice phase.

A deflection is the slope of a mass

The site has computed the exact relation between a geoid slope and a deflection of the vertical — one arcsecond is 4.85 millimetres in a kilometre — and computed no deflection anywhere, because a deflection is the gradient of a geoid and the geoid is a citation. State a mass instead, and every quantity is a closed form: a salt dome five kilometres across bends the plumb line by 2.21 arcseconds.

datums · Height Steepest descent on Lambert cylindrical, computed on the ground and on the page. Fifteen routes, each started at the same place twice. The solid line follows the true direction of steepest descent on the sphere; the dashed line follows the direction read off the page at every step, which is what an analysis of a projected grid does. Both take the same length of step on the ground, so the only difference between them is the direction. They part by up to 1517.8 kilometres, against a bearing error of 45.6°. The faint lines are contours of the field, which is seven caps at stated centres and widths.

Water runs downhill on the ground, not on the page

A drainage network is the set of steepest-descent trajectories of a field, so it is built entirely out of directions. A conformal map preserves those directions exactly and therefore preserves the whole network; an equal-area map does not, and sends a route up to 892 kilometres away from where the water actually goes.

distortion · Gradient A finer grid makes a measured slope worse. The error in the direction of steepest ascent, against the spacing the field was sampled at, for four noise levels. With exact values the curve falls at a fitted slope of 2.00 — second order, which is what a central difference is. Add noise and the same curve turns over: a finite difference divides the noise by the spacing, so halving the grid doubles the noise in the slope while quartering an error that was already negligible. The minimum is where the two meet, and it is not at the fine end.

The slope of a field that was measured

Three rungs differentiate a formula, which is what makes the projection the only thing under test. A real field is a grid of numbers with an error on each of them, and differencing such a thing divides the noise by the spacing — so a finer grid gives a worse slope, there is a best spacing, and it is the cube root of the noise.

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

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.

distortion · Gradient The points a projection cannot move. The critical points of a stated field, found twice: once on the sphere from its own gradient, and once in Lambert cylindrical's page coordinates from the page's own numbers, with nothing shared between the two searches. seven points, 3 maxima, 2 saddles and 2 minima, and the two sets agree in position to 2.3e-7 degrees and in type at every one. The page's Jacobian is invertible wherever the map is a map, so it sends a zero gradient to a zero gradient and cannot move a critical point anywhere.

What the page cannot move

Six rungs measure what a page does to a field's readings and every one of them moves. The critical points do not: found on the sphere and found again in a projection's own page coordinates with nothing shared between the searches, they agree to 10⁻⁷ degrees and in type at every point, on every projection. What the page does move is their shape, by a factor bounded exactly by the indicatrix's axis ratio squared.

distortion · Gradient One current, two maps, and only one of them conserves it. A flow that has no sources or sinks anywhere on the sphere — it is the perpendicular of a stream function's gradient, so its divergence is zero by construction — drawn on two projections. The arrows on the right are the images of the same ground velocities as the ones on the left. On the equal-area map the drawn field is still divergence-free; on the other it is not, and the arrows drawn in the warning colour are where a reader measuring the picture would find a source or a sink that is not there.

A current drawn on a page has sources

Seven rungs project a scalar field and ask what the page does to its gradient. A wind or a current is the other half of what gets mapped, and the operator that matters for it is the divergence — which is preserved by an equal-area map exactly, by no other map at all, and by a conformal map least of anybody's expectation.

distortion · Gradient

Named alongside it

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

ConformalityEqual-areaInvariantTest fieldThematic mappingDualityPrincipal directionScale factorAnisotropyAreal factorClosed formCritical point

All concepts