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

Every essay — page 6

Essays 121 to 144 of 339, in the same order.
A raster warped to Lambert azimuthal equal-area and back, nearest against bilinear. The left panel is the field the raster carries — a smooth analytic function, so that the error of an interpolation is the interpolation's error and not a photograph's history. The other panels are what is left after warping into Lambert azimuthal equal-area and back to Equirectangular, shown as the difference from the original at six times the contrast. Nothing moved: the coordinates go through the maps exactly. What is lost is that a target pixel's centre does not fall on a source pixel's centre, so a value has to be invented for it. Bilinear is closer to the field — RMS 0.0022 against 0.0212 — and has given up 0.63 per cent of its variance to get there. The panels are drawn at 48 by 32 cells; the measurement is made at the same resolution. What a machine does with it

Reprojecting a raster invents values

Moving a picture from one projection to another moves no coordinate — the maps are exact both ways. What is lost is that a target cell's centre does not land on a source cell's centre, so a value has to be made up for it, and the making-up has an order of convergence: 1.00 for nearest, 1.98 for bilinear, 2.93 for a cubic, measured by refining the grid.

4 figures · Dataset Clairaut's theorem, and the term it drops. Clairaut's theorem says the flattening of a rotating body plus the flattening of the gravity on it equals five halves of the ratio of centrifugal to gravitational acceleration at the equator. Measured on WGS84: f = 3.3528e-3, f* = 5.3024e-3, and their sum is 8.65525e-3 against the theorem's 8.62447e-3. The residual is 3.08e-5, which is 2.74 times f² — second order, which is what a first-order theorem is entitled to be wrong by. What the numbers refer to

The flattening is not a free parameter

An ellipsoid is usually presented as two numbers somebody fitted. One of them is not free — Clairaut's theorem relates the shape of a rotating body to the gravity on it, and the relation holds on WGS84 with a residual of 3.1×10⁻⁵ — which is 2.74 times f², exactly what a first-order theorem is entitled to.

7 figures · Ellipsoid The conformal map onto a square face of a cube. The spherical face of a cube carried onto its flat face by a map that is conformal everywhere — the measured angular deformation over the drawn interior is 1.63e-6°, which is the arithmetic's own floor. The rings and spokes are circles and radii on the sphere, and they cross at right angles here because that is what conformal means. The map was solved for as 16 terms of a series rather than written down: the face's edge comes out straight to 0.33 per cent of its own half-width, and that residual — not the conformality — is what more terms buy. At each corner the map behaves like ζ^0.75, so the scale factor there is infinite. The families

A conformal map onto a face

The polyhedral ladder ended owing a conformal face map, on the grounds that it needs elliptic functions. It does not: a conformal map of the sphere is an analytic function of one conformal coordinate, so the map is a power series, choosing it is a least-squares fit — and its scale factor is infinite at the corners, which is the angle deficit arriving as a singularity.

7 figures · Polyhedral The map that is continuous, and the pair it pays with. The orthographic is defined and continuous at every place on the Earth — it is written in the components of the place itself, with no longitude in it to jump. What it gives up is being one to one, and it gives it up almost everywhere: 47 per cent of the sphere shares its page point with the place directly behind it. Borsuk–Ulam guarantees at least one ANTIPODAL pair among those, and here it is exactly one — the centre and the place on the far side of the world, both at the middle of the picture, found to a residual of 1.5e-14. The impossibility

Two opposite places on the same spot

A projection may be continuous everywhere or one to one everywhere, and the first two rungs price both. What neither says is that the choice is not symmetric: a map that keeps continuity does not lose injectivity somewhere arbitrary. It loses it, always and at minimum, on a pair of places directly opposite each other on the Earth.

8 figures · Topology Everywhere within 200 km of the route from London to Tokyo. The shortest route between the two places, and the set of places within 200 kilometres of it, with both edges computed on the sphere and then projected. On the ground the set holds 3.949 million square kilometres — the band 3.823 and the two end caps 0.126, which between them make one disc of the corridor's own width. Drawn in Mollweide the corridor is visibly wider at one end than the other, and the ground it stands for is not. Paths and directions

A corridor has a width the page cannot keep

Buffering a line is the most-run operation in spatial analysis and the corridor it produces is a reach set with two failures nobody separates: stroked on the page, one width covers 190 to 471 kilometres of ground along a single route; computed in closed form, the formula stops being the area at a width the route's own length fixes, and eventually claims more ground than the sphere has.

7 figures · Reach What a solved map is, as a list of numbers. A map with no formula is a list of coefficients, and this is the list. The Chebyshev map of an elongated region, 30° by 10° has its coefficients falling by a factor of 1.5e+13 from the first to the fourteenth, so a table of a dozen numbers carries the whole projection; the conformal cube face's coefficients, marked separately, fall far more slowly because the map has a singularity at each corner. How fast this line falls is exactly how portable the map is — and neither map has a name, an inverse in closed form, or a formula anybody could quote. What each projection optimises

A map with no formula

The solved projection has no name, no formula and no closed-form inverse. It is fourteen numbers — and the rate at which those numbers fall away decides whether a map can be shipped at all: geometrically for a smooth region, and like a power for one with corners.

7 figures · Condition How badly two arcs determine the flattening. The reciprocal flattening recovered by inverting two measured degree lengths, at 1.5° and 66.33°, against an error introduced into the equatorial one. The exact pair returns 298.26. Ten metres of error — 90 parts per million of a 110-kilometre arc — returns 301.5, and the relative error in the flattening is 118 times the relative error in the arc. At 1000 metres the inversion returns a negative flattening: an Earth longer through the poles than across the equator, which is the answer the Paris Observatory defended for a generation. What the numbers refer to

The figure of the Earth was measured

Two expeditions went to Lapland and Peru to measure the length of a degree of latitude, and the whole signal separating a flattened Earth from a round one is a kilometre in a hundred and eleven. Inverting two arcs amplifies their error by 118 — and a kilometre of error returns a lemon-shaped planet.

7 figures · Ellipsoid Past the five solids: what more faces buy, and what they cost. The icosahedron subdivided 2, 3, 4, 6, 8 ways, giving 20, 80, 180, 320, 720, 1280 faces. The worst angular deformation inside a face falls from 13.1° to 0.26°, with a fitted exponent of -0.947 against the face count — the reciprocal, as it must be, because a face's angular size goes as the inverse square root of the count and the gnomonic's deformation goes as the square of that. The total cut length rises from 11.6 to 96 sphere radii, fitted at 0.509. Both exponents together say the whole economics of the family in one line: halving the distortion costs √2 times the cutting, for ever. The families

More faces, less distortion, more cutting

The five regular solids are where polyhedral mapping stops, and they stop because there are only five rather than for any reason about maps. Subdivide instead, and the family runs to any number of faces with two fitted exponents: distortion falls as the reciprocal of the count and cutting rises as its square root.

8 figures · Polyhedral five planes a dataset might be stored in, scored on three operations. Each candidate measured over -10° to 30° east and 35° to 60° north: the worst areal error, the worst angular deformation, and the spread of the scale factor, which are what an area query, a shape and a distance respectively depend on. The best plane for areas is Gall–Peters, for shapes Lambert conformal conic, and for distances Lambert conformal conic — three different answers, and no fourth candidate would collapse them, because a projection exact in two of these columns has a = b = 1 everywhere and is the isometry Gauss's theorem forbids. area of a polygon costs 3.06× too large in the wrong plane; drawing a line between two points costs 194 km from the ground it claims. What a machine does with it

The operation decides the coordinate system

Five candidate planes over one region, scored on the three things a spatial operation depends on. The conformal conic wins shape and distance and is 11.7 per cent out on area; the equal-area member is exact on area and 38.9° out on shape. No candidate is exact in two columns, and no candidate ever will be, because one that was would be an isometry.

5 figures · Dataset Two rules, three populations. The taught rule keys on latitude; the replacement keys on shape. On the thirty regions the replacement was read from it scores 83 per cent against the taught rule's 63. On forty-five different regions built the same way it scores 91 — higher, not lower, so the generalisation the shortfall doubted is real. On the seven named regions this collection actually uses, both rules score 29 per cent, and following either costs a mean factor of 14.6. The third bar of each group puts the taught rule's polar clause back into the shape rule, which is what the out-of-sample failures ask for: it repairs four of them, breaks two that were right, and reaches 43 of 45. What is taught wrongly

The rule scored out of sample

A replacement rule was read off thirty regions and scored on the same thirty, and this collection recorded that as not being evidence about any other thirty. It is: on forty-five different regions the rule scores 91 per cent against the 83 it managed at home. What it cannot do is the seven regions the collection actually uses, where both it and the rule it replaced name the winner twice out of seven and cost a mean factor of 14.6.

6 figures · Audit A published coordinate is a definition being quoted. Re-observing a control point perfectly, with an instrument good to 20 mm, and comparing the answer to what is published for it. The published coordinate was fixed on a datum realisation that has since been superseded, which accounts for 106 m, and the ground it marks has moved 750 mm in 30 years at 25 mm a year. Both numbers are larger than the observation and neither is an error in it: subtracting a published coordinate from an observed one measures the interval between two definitions. Grids, and what a survey does

A published coordinate is a result

Re-observe a control mark perfectly, with an instrument good to twenty millimetres, and the answer disagrees with the published value by a hundred metres. Neither number is wrong. The difference measures the interval between two definitions.

5 figures · Reduction Every aspect of Mercator over Japan, and the line the old sweep searched. The regional distortion of Mercator over Japan for every position of the projection's pole — darker is better — with the meridian the site's one-dimensional sweep searches drawn on it. The sweep's best is a gain of 24.4× over the normal aspect; the two-parameter search finds 61.4×, which is 2.51 times better again, at a pole 45° of longitude away from anything the sweep could reach. The shortfall recorded when the sweep was written said this would happen for a region whose long axis runs diagonally, and this is the measurement of it. What each projection optimises

The aspect has three numbers, not one

The site's aspect search has swung the projection's axis through one plane for a long time, and recorded that a region whose long axis runs diagonally has its optimum somewhere that plane never reaches. Searching the whole sphere of pole positions finds 2.6 times more improvement over Japan and 3.4 over the conterminous United States.

8 figures · Choosing What the drawn line costs, projection by projection. The ground length of the page-straight route from London to Tokyo, as a percentage above the shortest route. Two of these have names. On the gnomonic, centred on the route, the drawn line IS the shortest route, at -0.0000 per cent. On Mercator it is the rhumb, matching the rhumb's own length to 1 parts per million — which is why that projection exists. On the other eight it is a curve with no name and a cost between 0.0 and 19.6 per cent. Paths and directions

The line drawn straight on the page is a route

Eleven essays draw the route on the map. Nobody has drawn the map's own proposal: the ground curve somebody follows by laying a ruler on the page. It has a name on exactly two projections and is 18.2 per cent long on the one where it is famous.

6 figures · Paths The world after a map of degree 2. The graticule of the sphere, sent through the square of the stereographic coordinate and then drawn in Mollweide. Every parallel is still a closed curve and every meridian still runs pole to pole, and the whole world has been wrapped round the sphere 2 times: the 2 meridians that used to be 180° apart now lie on top of one another. Nothing has been torn. The degree measured by integrating the area the map sweeps is 2.0004. The impossibility

How many times, not whether

Three rungs of this ladder answer yes or no and have no other kind of answer. The degree is the first quantity here that counts: a continuous map of the sphere to itself covers it a whole number of times, injectivity forces that number to ±1, and the number is recoverable three ways — by counting preimages, by integrating swept area, and from the rate at which the preimages coalesce.

6 figures · Topology The circle a reprojection is about to turn into an ellipse. A circle drawn on the Mollweide plane, and its image under the map that carries that plane to Mercator — normalised so each pair has the same area, because the whole map can be rescaled and the shape is what is being shown. The dashed circle is what an undistorted reprojection would leave. This is Tissot's construction with a plane in place of the sphere, and it is the right picture for a reprojection because both ends are pictures. Worst angular deformation over the sampled points: 81.2°. Drawn on the source plane, in Mollweide. Measuring distortion

A projection between two projections

Every distortion measured on this site so far compares a map with the sphere. The operation a machine actually performs compares a map with another map — and that map has its own two principal scales, its own areal factor and its own angular deformation, none of which is the difference of the two it was built from.

8 figures · Tissot A reach set with a cost that depends on direction. Everywhere reachable in the time it takes to cover 3000 kilometres in still air, under a steady westerly, with the still-air set drawn inside it. The anisotropic set runs from 1653 kilometres against the flow to 4261 with it — a ratio of 2.58 — while the still-air set is round to 1.023, which is the lattice's own floor and not a shape. Drawn in an azimuthal equidistant centred on the source, so every radius on the page is a ground distance and none of the shape is the projection's. Paths and directions

A reach set with a cost that depends on direction

Three rungs build reach sets out of a distance, which is symmetric and isotropic by construction. Nothing anybody travels is: in a flow at 45 per cent of a vehicle's own speed the same vehicle gets 4,261 kilometres one way and 1,653 the other — a ratio of 2.58 — while the ground it covers grows by ten per cent.

6 figures · Reach Everywhere Albers is exactly right, and the band round it. The set on which both principal scale factors are one — the only ground where a ruler on this map, at the map's own stated scale, measures the true distance in every direction. It is the parallels at 20.000° and 60.000°, drawn as a curve, with the band within 0.01 of true scale shaded round it. That band is 2.673 per cent of the sphere, and it narrows as ε as the tolerance tightens. The curve at its centre has no width at all, and no tolerance makes it have one. The impossibility

The places where a map is exactly right

Nine essays on this ladder say a map cannot be right everywhere. None asks where it IS right — and the answer is a curve, a pair of curves, or two isolated places, never a patch. Measured across fourteen projections the set's neighbourhood shrinks with an exponent of 0.48, 1.0 or 2.0, and the value the impossibility forbids is 0.

6 figures · Curvature The scale asked for and the scale the pyramid has, at 0°. A tiling scheme exists only at integer zoom levels, a factor of two apart in resolution, so a request for any scale between them is answered by the nearest rung. The ratio runs from 0.707 to 1.405 — 1/√2 to √2 — and repeats identically at every doubling, which is four times that in area. A request for 1:10,000 is served at zoom 16, which is 1:8,531; A request for 1:25,000 is served at zoom 14, which is 1:34,124; A request for 1:50,000 is served at zoom 13, which is 1:68,247. Nothing anywhere reports it, because the map that arrives is a perfectly good map of something. What a machine does with it

Zoom is a ladder

A tiling scheme exists only at integer zoom levels a factor of two apart, so a request for 1:25,000 is answered with 1:34,124 — 36 per cent coarser, and 86 per cent coarser in area. The mismatch runs from 1/√2 to √2 and repeats identically at every doubling, and nothing anywhere reports it, because the map that arrives is a perfectly good map of something.

7 figures · Screen The boundary of the region the series is the map in. The curve where the transverse coordinate reaches 2.918, which is where the terms stop shrinking. It crosses the equator 83.81° from the central meridian and closes towards the poles, because the same longitude is a smaller transverse coordinate at a higher latitude — the boundary is a curve rather than a meridian. The narrow band beside it is a 3° zone, the width national grids actually use, drawn to the same scale: the practical world sits in about a fiftieth of what the series can reach. Drawn in Mollweide. The families

Where the series stops being the map

The transverse Mercator has no closed form on an ellipsoid, so every national grid computes a truncated series. Asking how many terms it needs has an answer everywhere; asking whether more terms always help has an answer only within 83.81° of the central meridian, and the boundary is computed rather than assumed.

9 figures · Ellipsoid The same map on four differently prepared pages. Mollweide at 30°E 20°N, then the same projection with a rotation and a magnification applied to the page, then with one axis stretched by 1.6, then with a shear of 0.5. The similarity changes nothing: flexion, skewness, ω and the anisotropy are identical to every printed figure, and only the last column — the same turning measured per unit of page arc rather than per unit of ground arc — moves, by exactly the magnification. The stretch and the shear change all of them, and flexion by 35 per cent. Measuring distortion

The second derivative is not an invariant

The first-order ladder established which quantities survive a change of coordinates and which are artefacts of the parameterisation. Asked of the second order, the answer is that flexion survives a rotation and a magnification of the page exactly, and survives nothing else: a stretch of 1.6 in one axis moves it by eleven per cent, a shear of 0.5 by thirty-five, and a shear of 3 leaves a ranking of eight world maps with a rank correlation of −0.07 to the one it started with.

6 figures · Flexion How thinly a sphere can be covered by a few equal caps. The covering density of the best arrangement of n equal caps found for each n — the total area of the caps divided by the sphere's, so a value of one would be a perfect tiling with no overlap. The horizontal line is 2π/√27 = 1.2092, the thinnest covering density of the PLANE by equal discs, which this site has used for the sphere since its first atlas essay. It is wrong in both directions: at 2 caps the sphere is covered more thinly than any plane can be, because a cap may be a hemisphere, and at every count from 3 upwards more thickly — 1.5092 at 3, and 1.3377 at 14. The ringed points are the four counts whose optimum is proved: 2 at 90.00°, 4 at 70.53°, 6 at 54.74°, 12 at 37.38°. Everything else is an upper bound from a search, drawn as one, and the bound loosens as the count rises — the search reaches the proved optimum to 3.4 per cent at these counts and has no such check anywhere else. What each projection optimises

The sphere is not the plane at small counts

The site's atlas arithmetic multiplies an ideal sheet count by 2π/√27, the thinnest covering density of the plane. At the four counts whose optimal covering of the sphere is a theorem the plane's number is 21 per cent high at two caps and 9, 5 and 2 per cent low at four, six and twelve — wrong in both directions, and the direction changes with the count.

8 figures · Choosing Two parameter sets 100 metres apart, over the region they were fitted to. Each marker is drawn at a size proportional to how far the two transformations put it apart. The second set differs from the first by 100 metres of translation along the direction this network can least see, with the rotations and the scale re-fitted to absorb it — which is what a second agency's adjustment does when it chooses a different constraint. The worst disagreement anywhere in the region is 5.64 metres and the mean is 3.61. Applied at south-eastern Australia the same two sets differ by 193 metres, because the rotation that absorbed the translation here is a rotation of the whole Earth. What the numbers refer to

Two parameter sets, one transformation

Agencies publish seven-parameter datum transformations that differ by hundreds of metres in translation, and the usual reading is that one of them is better. Over the region either was fitted to they are the same transformation: a hundred metres of translation, re-absorbed by the rotations and the scale, moves a British coordinate by 5.6 metres and an Australian one by 193.

8 figures · Datum The same ranking, with the right answer removed from the library. A map drawn in Mercator, fitted by every candidate except Mercator. Something still wins: Conformal conic, by a factor of 1.64 over the runner-up, leaving 0.4 per cent of the map's width unexplained. With Mercator in the library the winner's margin is 1.2e+13. The ranking always produces a name; what tells the two situations apart is how far ahead the name is. What is taught wrongly

When the answer is not in the library

Fitting twenty candidates to a map and ranking the residuals always produces a winner, which makes it a ceremony unless it can also produce a refusal. Held out of its own library, a Mercator map is named as a conformal conic, leaving 0.4 per cent of the map's width unexplained — and the quantity that tells the two situations apart is not the residual but the margin, which is 1.6 when the truth is absent and 10¹³ when it is present.

6 figures · Identify A 1-hectare parcel across British National Grid, at 52°N. The departure of a plan's dimensions from the ground's, across the width of the zone, as a length and as an area. The area curve is the length curve doubled — an areal scale factor is the square of a linear one, and squaring a small departure doubles it — so a parcel whose sides are each 399 parts per million short on the central meridian is 797 short in area. On a 1-hectare parcel that is 8.0 square metres at the central meridian and 8.0 at 0.00° out. Grids, and what a survey does

An area on the grid is not an area on the ground

A grid's scale factor is a property of lengths and what a surveyor sells is an area. Squaring a departure doubles it, so a hectare drawn on the British grid at its central meridian has 10,008 square metres of ground under it — and correcting an area with the line factor instead of its square leaves half the error behind.

5 figures · Grid