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

The thread: A theorem, not a limitation — page 3

No map is faithful, and this is not a shortcoming awaiting a cleverer cartographer. Gaussian curvature is intrinsic, a sphere has some and a plane has none, and that settles it. Essays 49 to 72 of 106.
The least error a flat picture can have, against how much sphere it spans. The same configuration of places, shrunk about its own centroid so that every bearing is kept and only the span changes, with the least worst-case relative error of the best flat picture at each size. Both axes are logarithmic. The fitted slope over the rows below ninety degrees is 2.0246: the error falls as the SQUARE of the diameter. That is Gauss's theorem arriving as a number for a finite set — curvature is a second derivative, so its first effect on a distance is quadratic in the separation — and it is why a county fits on a sheet and a hemisphere does not. The impossibility

How wrong a flat picture has to be

A determinant proves no flat picture of four places is exact and leaves the size of the failure where nobody can read it. Measured directly, the least error falls as the square of how much sphere the places span — fitted exponent 2.0088 — and the same exponent comes back from five different arrangements while the constant in front of it moves by a factor of seventeen.

 The pass between the basins, measured rather than bracketed. The height of the lowest path from one near-optimal basin to the other, for regions of growing size. It is found by sorting the score surface and joining cells in order, so it is the minimax path's own height rather than the level at which a bisection stops finding two pieces. The fitted exponent is -1.34 at a coefficient of determination of 0.986. The earlier level-set measurement on the same projection gave -1.57 and the argument predicts −1, so measuring the height directly moves the answer towards the prediction, though not onto it. What each projection optimises

The height of the pass between two basins

One thing was left owed: the fracture threshold falls faster than the argument predicts, and the difference was supposed to be the level set's blindness to the height of the pass. Measured directly, the pass between the two deepest basins falls as the −1.34 power at R-squared 0.986 — nearer the predicted −1 than the level set's −1.57 on the same projection, and not onto it. It is not the pass that breaks the set first, and at five of six sizes the second basin is the first one's mirror twin.

 The corner against where along the seam the feature crosses. A feature crossing the seam of a cube at right angles, drawn on each face's own map and unfolded, with the crossing moved along the edge. Every curve starts at zero, because the face's mirror through the edge's own midpoint reverses the along-edge direction and forces the shear term in the Jacobian to be odd. Away from it the gnomonic — the map measured earlier as having no corner at all — reaches 44.4°, three times the equal-area map's and far past the conformal one's. The edge's half-length is 35.3°, so the right-hand end is still well inside it. The families

The corner is not at the midpoint

An earlier measurement took the corner a feature gets crossing a polyhedral seam, and found none at all under the gnomonic face map. It crossed at the edge's own midpoint every time — the one point on the edge where a face's own mirror symmetry forces the corner to vanish. Two fifths of the way to the vertex the gnomonic gives 20.15°, the equal-area map 7.28° and the conformal map under a degree, which reverses the ordering entirely.

 A conformal map of Vesta, in the coordinates that make one possible. The body's own isothermal coordinate, used as the map — which is Mercator's construction carried out on a body that has no axis to be Mercator about. The lines are Jacobi's ellipsoidal coordinates, which are the surface's lines of curvature, at 30° and 20° spacing; the two quadratures that turn them into an isothermal pair are one-dimensional. The measured angular deformation away from the marked points is 1.1e-6°, which is this site's noise floor, and the areal factor spans a factor of 162 — conformal, and emphatically not equal-area. The four marks are the umbilics, where the coordinate collapses and the map has nothing to say. What the numbers refer to

A conformal map of a body with three axes

This site built an equal-area map of a triaxial body and wrote down what it could not do: the conformal one, which needs an isothermal coordinate that a surface with no axis of revolution was said not to have. It has one, Jacobi found it in 1839 for a different reason, and it takes two one-dimensional integrals.

 Everywhere reachable in 1 h 15 min on roads 5 km apart. The ground reachable in 1 h 15 min from a road crossing, with roads every 5 km in two directions at 80 km/h and the ground between them at 20 km/h, computed on a one-kilometre lattice. It reaches 100 km along a road and 67 km along the diagonal, and covers 18,873 km². The dashed diamond is the shape the set converges to as the roads get denser, whose area is 20,000 km². The dotted ring is what a range ring at road speed claims — 31,416 km², of which 39.9% cannot be reached — and the small dotted circle is what the ground alone would give. Paths and directions

The reach set takes the shape of the roads

On open ground a reach set is a disc. On a grid of roads it converges to the convex hull of a diamond and a disc, and once the roads are √2 times faster than the ground between them the ground's speed stops mattering at all. A range ring drawn at road speed then claims 57 per cent more ground than anybody can reach.

 Two members, their average, and the family's best answer to it. Two equal-area conics at cone constants 0.25 and 0.85, the average of the two, and the member of the family nearest that average — at 0.540, which is not the parameter midpoint. The average is not a conic at all: its parallels are still arcs but they are arcs of circles about different centres, so no single cone constant reproduces it and the residual is 0.1722. The families

A family is not closed under averaging

A family is a set of maps with a parameter running through it, and the compromise projections here are averages of members. An average of two members is not a member: two conics far apart average to something 31 per cent of their own separation outside the family — and averaging within a family makes the map worse every time, while averaging across two makes it 16 per cent better.

 The share of its distances a flat picture can hold. A picture of n places on a plane has 2n coordinates and is unchanged by two translations and a rotation, so 2n − 3 numbers in it are genuinely free; each distance held exactly is one equation. So at most 2n − 3 of the n(n−1)/2 distances can be right, whatever the places are and however hard the map maker tries. The share falls as 4/n: five of six for four places, thirteen of twenty-eight for eight, and 197 of 4,950 — 4.0% — for a hundred. The dashed curve is 4/n, which the count approaches from below and never crosses. The impossibility

Five distances of six, and never more

A flat picture of n places has 2n − 3 free numbers and n(n−1)/2 distances, so it can be exactly right about 2n − 3 of them and no more — five of six for four places, twenty-nine of a hundred and twenty for sixteen, and 3.98 per cent for a hundred. The bound is reached by a construction with no fitting in it, and an azimuthal equidistant map spends only n − 1.

 What 24 bearings draw beside a river crossable at eight bridges. Everywhere within a 10 km walk of the centre, with a river 3 km to the north crossable only at bridges every 2.5 km. Short of the river the reachable ground is the disc; beyond it, each bridge reaches a half-disc of whatever walk is left, and those half-discs are 30.3% of all the reachable ground. The solid outline is the polygon a fan of 24 bearings draws: every bearing that meets the river stops at the bank, so the polygon holds none of the ground beyond it and is short by 31.6% in total. The river is drawn wider than its 80 m. Paths and directions

A drawn reach set stops at the river

A polygon drawn through a fan of bearings round a reach set is short by a chord deficit that falls as the square of the count — on ground with nothing in the way. Put a river three kilometres off, crossable at bridges, and the polygon misses 30 per cent of the reachable ground at every count from forty-eight bearings to a thousand, and a convergence check passes on it.

 A conformal map solved rather than written down. The same patch of parameters on two bodies, mapped to the plane by solving the discrete Cauchy–Riemann equations — one complex equation per triangle, 1568 triangles, least squares, conjugate gradients, and no formula for either surface. Left: a sphere, where the answer is known in closed form and is not used. Right: a body with a bump on it, which has no isothermal coordinate and therefore no closed form at all. The parameter lines cross at right angles in both, to a median of 0.60° and 1.25° of angular deformation. What the numbers refer to

A conformal map of a body that is not a quadric

Jacobi's ellipsoidal coordinates give a triaxial body a conformal map by two quadratures, and this collection wrote down what that argument uses: the surface has to be a quadric. A real body is not. Solving the discrete Cauchy–Riemann equations instead — one complex equation per triangle, two thousand triangles, conjugate gradients — gives a conformal map of a bumped body to a median of 1.10°, converging at first order in the mesh, with the areal factor spreading by 3.07 and refusing to converge at all.

 One of these is a corner. The corner reported at the exact conformal seam, and the corner reported at a cube's gnomonic seam fifteen degrees along its edge, each measured over a shrinking arc. The gnomonic's is 21.4572° at every span from twenty-four degrees down to three — the same number to four decimals — because it is a corner. The conformal one halves whenever the span does, fitted exponent 0.990, because a tangent read from a finite chord of a curved image departs from the true tangent in proportion to the chord. It is not a corner; it is the instrument. The families

The exact map says the seam is smooth

The conformal seam's corner could only be bounded at about two degrees, because the series it was fitted with holds its boundary condition to three parts in a thousand. An exact map exists — the stereographic projection composed with ∫dt/√(1 − t⁴) — and it settles it: the corner falls in exact proportion to the arc the tangent is read over, fitted exponent 0.99, while the gnomonic's 21.4572° is the same to four decimals at every span.

 The page's curvature is not a free parameter. The same trilateration, carried out on a sphere of stated radius instead of on a plane, with the error of the distances the construction decides rather than holds. At the radius the distances were measured on it is 5.1e-15 — exact, by construction, which is the refusal this figure exists to make. Two per cent either side of it the worst decided distance is already out by around 10%. Below 90 per cent of the Earth's radius the construction cannot be completed at all: two circles that must cross do not, and the page is simply too small to hold the places. The flat sheet is the right-hand limit, at 152%. The impossibility

The escape is not a dimension

Four places that will not lie in a plane surely lie in a space — and they do not. The double-centred matrix of their great-circle distances has a negative eigenvalue, so no Euclidean space of any dimension holds them, and buying a third dimension improves the picture by nothing whatever, to fifteen decimal places. What does work is a page with curvature, and the curvature is pinned to within two per cent.

 The basin has three widths, and they differ by a factor of 4.6. Two sections through the near-optimal basin of the Robinson aspect over Japan, drawn at one scale. Each ellipse is the set of aspects whose score is twice the optimum's, from the objective's own second derivative at the optimum: 32.1°, 16.0°, 6.9° along the three principal directions. A single number for "the width of the basin" is the cube root of their product, 15.3°, and it is not any of them. What each projection optimises

The basins have widths as well as depths

Measuring the height of the pass left one thing owed: shape means widths too. Measured, the basin has three of them — 32°, 16° and 7° at Japan — it gets wider rather than narrower as the region grows, and the exponent it predicts overshoots the measured one by half again.

 The piece count rises and falls. The number of connected pieces of the near-optimal aspect set for Robinson over Japan, swept finely through the threshold rather than sampled once below it. It is one piece at a wide threshold, reaches 13 at 1.010, and returns to one as the set shrinks onto the single best aspect. The set first disconnects at 1.941, which is above the peak: the pieces keep multiplying after the first break. This is the sweep the earlier measurement could not afford and it costs one grid, because every threshold reads the same 4992 evaluations. What each projection optimises

The threshold is not a percolation

The near-optimal aspect set was found breaking into twelve pieces rather than two, the transition was called a percolation, and the exponent went unmeasured. Swept finely, the piece count rises from one to thirteen and falls back to one — and refining the grid by a factor of fifteen moves the peak from ten to nineteen, at the rate the surface's own minima grow, while an uncorrelated field on the same lattice grows by a factor of twelve.

 The quickest way out of the ice bends at the edge. A ship twelve kilometres inside a straight ice edge, making 6 km/h in the ice and 25 km/h in open water, bound for a point ten kilometres out into the water and twenty along. The dashed line is the straight route, 3.244 hours. The solid one is the quickest, 2.849 hours — 24 min sooner — and it reaches the edge only 2.55 km along, leaving the slow ice as directly as it can. At the edge it meets the perpendicular at 12.0° in the ice and leaves at 60.2° in the water, and the sines of those angles stand in the ratio of the speeds to within 4e-17. Paths and directions

A crossing bends by a law only a conformal chart can show

A ship leaving pack ice for open water takes its quickest route by bending at the ice edge, so that the sines of its angles to the perpendicular stand in the ratio of its two speeds. The law has nothing in it but angles, so on Mercator it can be checked with a protractor — and on a plate carrée at sixty degrees north the same crossing reads as a ship making ten kilometres an hour in the ice rather than six.

 Where the condition holds, and where it was asked to. The boundary scale of a fit collocated at 20 points, drawn all the way round the boundary. The marked points are the ones the condition was imposed at, and the curve passes very near zero at every one of them; between them it does not. The largest departure on the samples is 3.92e-5 and the largest anywhere is 3.27e-4, and the second is the one the map has. What each projection optimises

A condition imposed at points is not a condition

Every projection written as a condition is solved by imposing that condition at a finite set of samples because that is what a linear system is. With barely more equations than unknowns the residual the solver reports is 8.3 times too good — and refining the collocation twentyfold does not improve the map at all, it only makes the report honest.

 Three face maps on a cube, over a shrinking span. The corner a feature gets crossing the seam of a cube, measured by reading the tangent over an arc and then shrinking the arc by a factor of sixteen. A chord differs from a tangent in proportion to the arc, so a perfectly smooth join reports a corner that HALVES when the span halves — a slope of one on these axes. None of the three lines has a slope of one. The fitted slopes are 0.000, 0.000, -0.018, which is a flat line in each case, and a flat line is a real corner. The three differ in size and not in kind: 20.1513°, 7.2772°, 0.7687°. The families

The span ladder, run on all five

A recorded shortfall said the exact map's span ladder was one loop away from settling the seam on the five Platonic solids. The loop was run and the answer is the other one: the series conformal face map has a corner that does not shrink with the measurement span on any of them — 4.10° on the tetrahedron, 0.77 on the cube, 0.004 on the icosahedron — and it belongs to the truncation rather than to conformality.

 A radius of 37.7 km makes a flat picture of four cities exact. Each row is one pair of London, New York, Tokyo, Sydney. The bar is the range of distances between two places of radius 37.7 km — the true centre-to-centre distance, plus or minus 75.4 km — and the solid dot is where the best picture for places of that size draws the pair. Every dot is inside its bar, and the pairs at the bars' ends are the ones that fix the radius. The hollow dots are the picture that is best about ratios instead: its worst pair is 103.3 km off, so it would need places of 51.6 km. The best picture of points is not the best picture of places. The impossibility

A place with a size can be drawn to scale

Points on a sphere have distances no flat picture holds. Places are not points: give every place a radius and the distance between two of them becomes a range, and some flat picture is right about every range once the radius passes a threshold — 38 kilometres for London, New York, Tokyo and Sydney, six metres for five towns in Britain, growing as the cube of the set.

 Walking the crossing all the way to a cube's vertex. The corner a feature gets crossing the seam, at crossings that approach the vertex geometrically — the last is 3.5e-4 degrees from it. Every face map is singular at a vertex, so the expectation is that the corner runs away. It does not. The gnomonic settles at 53.1295° and the equal-area map at 12.9656°, and both are finite. The dashed line is the solid's angle deficit, 90.0°, which is what the surface loses at that point and is not what either map's corner reaches. The families

The corner that is the curvature

Every face map is singular at a vertex, so a corner that grows as the crossing walks towards one should run away. It does not: the gnomonic settles at 53.1295° on a cube and the equal-area map at 12.9656, both finite, both reached like the first power of the remaining gap. What does not settle is the deficit beside them — 90 degrees, fixed by Descartes before any projection is chosen.

 One line cannot point both ways between New York and Tokyo. The great-circle bearing from New York to Tokyo is 333.0°, and from Tokyo back to New York 25.1°. A flat picture draws one straight line between them, whose two directions differ by exactly 180°; the two bearings differ by 180° and another 127.9°. The dashed direction is the best a straight line can do for both ends at once, and it is 63.9° from the truth at each — which is a floor on the bearing error of every flat picture that contains these two cities. The impossibility

A flat picture has one direction between two places, and the Earth has two

The bearing from New York to Tokyo and the bearing back are not reverses of each other: they miss by 128 degrees. A flat picture draws one line between the two cities, so no picture containing both can have a bearing error below 64 degrees. On any set of places up to about fifty degrees across, the best flat picture sits exactly on that floor — and away from the poles, so does Mercator.

 In a jet the sixty-minute sets stop moving together. The North Atlantic airfields with a southerly jet of 200 km/h whose core, the solid line, runs north past the southern tip of Greenland; dotted lines mark where it has slowed to 150, 100 and 50 km/h. Faint dashed circles are the still-air sets of 745 km. Solid outlines are the places within sixty minutes of each airfield in the jet, computed as quickest flights through it: each is a circle moved by the wind it meets, and airfields on opposite sides of a gap meet different winds. The least allowance at which New York is joined to London moves from 624.7 km of still-air radius to 587.7 (−36.9 km), pinned between Goose Bay–Narsarsuaq, whose sets first touch at the marked point. A uniform southerly of the same speed moves it +3.0 km. Paths and directions

A jet moves the floor that a uniform wind cannot

A uniform wind leaves the least radius at which the North Atlantic can be crossed almost where it was, because it blows the same way on both halves of every gap between airfields. A jet does not. A southerly of 200 km/h past the southern tip of Greenland lowers that floor by 36.9 kilometres, the same jet 400 kilometres further west raises it by 23.3, and one integral of the wind along each gap predicts every such move to within a few kilometres.

 50° S to 40° S, drawn conformally without a cut and with one. The band of latitude from 50° S to 40° S, drawn two ways that both keep every angle. On the left it is not cut: the only such map with least scale variation is the stereographic from the south pole, the ring closes up, and its scale varies by a factor of 1.0750 across the band. On the right it is cut along one meridian and opened into a conformal conic with cone constant 0.7080: the ring spans only 254.9° and leaves a gap of 105.1°, and its scale varies by 1.0038. Not cutting costs a factor of 1.071 in scale variation. The impossibility

A ring can be drawn whole, and only one way

A band of latitude has a hole in it, and a region with a hole has a number no angle-keeping map can change. Draw the band without cutting it and that number forces the map: the stereographic from one pole is the least-varying picture there is. One cut along a meridian frees the cone constant, and the scale variation the cut buys back runs from nothing on a polar cap to a factor of 2.13 between the tropics.

 Five bearings of four cities, held exactly. A flat picture of London, New York, Tokyo, Sydney holding five one-way bearings exactly, each drawn as an arrow from the city it is measured at: London to Tokyo, Sydney to London, New York to Tokyo, New York to Sydney, Tokyo to Sydney. Five is the most any flat picture of four places can hold, and of the 192 ways of choosing five pairs and a direction on each, this is one of eight that can be drawn with every arrow pointing forwards; all eight leave London–New York, dashed, unheld. The price is paid by the other seven bearings: the worst, London to New York, is 165.0° from the truth. The impossibility

Five bearings of twelve, and only eight ways to draw them

A flat picture of four places can hold five of their twelve one-way bearings exactly, the same count as distances, and the five are found by solving equations rather than searching. The equations cannot tell a line from its reverse. Of the 192 ways to choose five, eight can be drawn with every arrow pointing forwards, all eight leave London and New York unheld, and on sixteen cities no choice made at random can be drawn at all.

 The solve's cost is a U in the shape; the curvature is not. Solid: how many conjugate-gradient steps the conformal solve needs, against the window's aspect ratio, at constant surface area. It has a minimum at the square window — 174 steps — and rises to 265 and 275 at the two extremes, which are elongated by the same factor in opposite directions. Dashed: the total Gaussian curvature the window encloses, on its own scale, which falls from 0.311 to -0.012 across the same sweep and is least at one end of it. The cost has its minimum where the window is square and the curvature has its minimum somewhere else, so whatever is making the solve expensive is not what is making the map spread. What the numbers refer to

A long window and a square one

Rung nine held the patch's area and found that where the cut goes still changes the map by a third, with the curvature it encloses predicting the change at r = 0.969. It recorded that it had held area and curvature and not shape. Sweeping the shape at constant area separates two things that had looked like one: the map is curvature and the cost is shape.

 Four lines of position, two unknowns, and four gaps that do not close. Four bodies at azimuths of 20°, 110°, 200°, 290° from the ship, each sight given an independent error of about 2 km along its azimuth. No point lies on all four lines. The least-squares point — the hollow dot — is the one whose squared distances to the four lines are smallest, and the dashed perpendiculars are those distances: 0.31, 2.83, 0.32, 2.83 km. The ship, the solid dot, is 3.21 km away. Three lines would have left a triangle and no residual worth the name; four leave two degrees of freedom, and the four gaps are the first quantity here that can be read as a measurement of the sights themselves. Paths and directions

The residual reports the error the fix was immune to

A fourth sight gives a four-line fit something three lines never had: a leftover. The leftover is a real measurement of the sights, and it is useless exactly where it is needed. One number computed from the four azimuths alone — before any star is observed — says how much of a common error reaches the residual, and it is 1 when the bodies stand at the quarters of the compass, where such an error moves the fix by nothing at all, and 0.004 when they are bunched within sixty degrees, where it moves the fix by two kilometres.

All threads