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

The thread: Computed, not quoted — page 10

Scale factors, areal ratios, angular deformations, geodesic distances and the areas of countries are all computed from the projection's own formulae. None is a number recalled from a table. Essays 217 to 240 of 383.
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.

 Sixteen cities, each left out in turn. The least worst-case distance error any flat picture of these sixteen cities can have is 53.3%; the best projection in the library, Azimuthal equidistant, manages 57.6%. Each bar is that least error with one city left out. Leaving out Sydney brings it to 25.5%; leaving out Cape Town, the next most costly, to 46.7%; and leaving out any of five — New York, Tokyo, Anchorage, Singapore, Los Angeles — changes it not at all. The note on a row counts the pairs holding the picture at its worst that the city belongs to. The impossibility

The error belongs to a few of the places

The least error of any flat picture of sixteen cities is 53.3 per cent, and it is not spread over the sixteen. Six pairs hold it there and Sydney is in three of them: leave Sydney out and it falls to 25.5, while leaving out any of five other cities changes it not at all. Choosing which places to draw is worth more than any choice of how to draw them.

 One unit of a vector tile, in metres of ground. A vector tile's coordinates are integers on a lattice 4096 units across the tile, and the tile halves at every level, so one unit is a distance that halves too: 5.48 m at z10 and 0.086 m at z16, at 55°. It is also a different distance at every latitude, by cos φ, because the tile is in Web Mercator — the same factor that makes a grid metre a different quantity of ground at every latitude, arriving in the file format rather than in the projection. What a machine does with it

A vector tile has an integer grid

Six essays on this ladder treat a vector tile as the thing a raster tile is not: geometry, resolution-free, styled at draw time. Its coordinates are integers on a lattice 4,096 units across a tile, the tile halves at every level, and at 55° north one unit is 88 metres at zoom 6 and 21 millimetres at zoom 18.

 The same patch, pinned six ways. Two vertices have to be held or the conformal energy has a similarity's worth of null space. Which two turns out to decide two of the three numbers reported. The median angular deformation is the same to 8 per cent across all six — that is the map. The areal spread runs from 2.68 to 7.44, so the 3.07 reported for this body was a statement about its corners. And pinning two adjacent vertices, which fixes the similarity through a very short lever, ruins the worst point without touching the median: a badly conditioned constraint pays for its scale in one corner. What the numbers refer to

The map depends on where it was cut

Solving the discrete conformal equations on a triangulated body gave an areal spread of 3.07, with the number possibly belonging to the patch, the boundary and the two pinned vertices rather than to the surface. It belongs to the pins: hold a different pair and it runs from 2.68 to 7.44, while the typical angular deformation does not move at all.

 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 projections that run in a circle. Mercator beats Sinusoidal beats Eckert IV beats Mercator, each on a majority of the same seven criteria over the whole sphere. Every margin is four to three, the narrowest a majority of seven can be, and the criteria that decide each edge are different ones. There is no way to place these three in an order that agrees with all three comparisons, and the obstruction is not a measurement error: every number is exact to the precision the sampler reaches. What is taught wrongly

The ranking is not an order

A weighting can be made to put almost any projection first. Remove the weights entirely, let each of the seven criteria vote once, and the answer is worse: over the whole sphere Mercator beats the sinusoidal, the sinusoidal beats Eckert IV, and Eckert IV beats Mercator — four such circles, every margin four to three, with a Condorcet winner sitting above them all.

 The line a commission can actually run. A boundary described as a parallel of latitude and marked by monuments 220 kilometres apart, with the offset exaggerated 700 times so that it can be seen at all. A commission cannot run a parallel: it can set a monument, sight a straight line to the next and clear the trees between, and a straight line between two points of equal latitude is a geodesic, which passes POLEWARD of the parallel everywhere between them. So the marked line lies north of the described one, by 1176.5 metres at the middle of each of its 9 chords, and encloses 1611.2 square kilometres that the words put on the other side. At the 20-kilometre spacing used throughout, the same offset is 9.03 metres. Grids, and what a survey does

The line a commission can actually run

A boundary commission cannot run a parallel of latitude. It can sight a straight line between monuments, and a straight line between two points of equal latitude passes poleward of the parallel — by s² tan φ / 8R, which at a mile of spacing is fifty-eight millimetres and at a hundred kilometres is two hundred and twenty-six metres. The described line and the marked line are different curves, and the marked one governs.

 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.

 In a wind every circle moves upwind, and the route moves with them. New York to London with the sixty-minute allowance at a stated one-engine speed of 745 km/h. In still air every airfield's circle has a radius of 745 km, drawn faint and dashed. In a northerly of 150 km/h the places from which each airfield can be reached in sixty minutes form a circle of the same radius moved 150 km upwind, drawn solid. The dashed route is the shortest the still-air rule allows, 5,698 km; from one point of it the quickest diversion in this wind takes 68.7 minutes. The solid route stays inside the moved circles, 5,765 km, and its longest diversion takes 60.8 minutes — the sixty the moved circles stand for, to within what the Earth's curvature does to a wind that keeps one bearing. Paths and directions

A diversion allowance in a wind is the same circle, moved upwind

A diversion rule written in minutes is turned into a distance by assuming still air. In a uniform wind the places from which an airfield can be reached in sixty minutes are not a stretched region but the still-air circle moved upwind by the wind's run. A 150 km/h northerly makes the still-air route across the North Atlantic need a diversion of 68.7 minutes, a southerly makes it 54.0 — and the least radius at which the ocean can be crossed barely moves.

 How many projections the map could be in. The number of candidates whose residual sits below the measurement noise, against the size of the region, at four noise levels. At one per cent of the map's width — a hand-digitised graticule — a four-degree region admits ten of the twenty candidates and a forty-degree one admits exactly one. Every curve falls, none of them crosses another, and all four end at one: identification works, and what it needs is extent rather than precision. What is taught wrongly

The answer is a set

Fitting a projection to control points produces a best fit — one projection, ranked first, with a margin. A best fit without a spread is not a measurement, and the spread is free: a control point has a width, and every candidate whose residual is inside that width has not been ruled out. At one per cent noise a four-degree region admits ten of twenty candidates and a forty-degree one admits exactly one.

 The threshold, and the two things it is a ratio of. The near-optimal set's fracture threshold on Robinson, against the size of the region, with the two quantities it is a ratio of drawn beside it. The threshold falls with fitted slope -1.341. The best score a region admits at all rises with slope 0.923 — a bigger region is harder to map — and the absolute score of the pass falls with slope -0.418. The first is the sum of the other two by construction, and the arithmetic says which of them is doing the work: the denominator carries 69 per cent of it. What each projection optimises

The first break is mostly its denominator

Three earlier measurements have fitted the near-optimal set's fracture threshold against region size and read the answer as a statement about the landscape. It is a ratio, and separating it takes one multiplication: the best score the region admits at all rises with exponent 0.92 and carries 65 to 70 per cent of the threshold's movement on three projections, while the pass's own depth falls slowly and steadily, with exponent −0.42 to −0.49.

 The drift is a straight line in the deflection. The horizontal distance between where a plumb line hangs at the top of a column of rock and where it hangs at the bottom, against the deflection of the vertical the mass produces at the surface. Four heights of column. Every line is straight through the origin: 12.54 mm of drift per arcsecond of deflection over a 4,000 m line, to two parts in ten thousand across a fortyfold range of deflection. Which is what makes the number transferable — the 47 mm the first calculation started from was a statement about one buried sphere, and this is a statement about any mass that produces the same deflection. What the numbers refer to

How far the plumb line bends

An earlier measurement dropped a plumb line down a four-kilometre column of rock beside one buried mass and found it arrived 47 millimetres from the point below the summit. That is a number about that mass. Parameterising by the deflection of the vertical instead — the quantity surveyors actually measure — gives 12.54 mm per arcsecond, exactly linear across a fortyfold range.

 The piece two schemes share, clipped rather than assumed. A cell of a gnomonic cube, whose four edges are great-circle arcs because a straight line on a gnomonic face is one, against a cell of a longitude–latitude grid, whose north and south edges are parallels and are not. Their overlap is neither a rectangle nor a spherical polygon of any standard kind, and it is 5965687 km² of the cube cell's 5965687 km² — 100.0 per cent. Computing it needs the arc of one boundary intersected with the plane of the other, which is three equations and two roots, and it is exact. What a machine does with it

Cells that are rectangles in no coordinate

What moving a field between two cell schemes costs was measured between two schemes whose cells are longitude–latitude rectangles — which is what made every overlap a rectangle with a closed-form area. The schemes anybody actually argues about have cells that are rectangles in no coordinate, and their overlaps have to be clipped.

 What a geoid model leaves out, against the degree it stops at. The RMS of everything above the model's highest degree, from Kaula's rule — the statement that the normalised coefficients at degree n are about 10⁻⁵/n². The line is R × 10⁻⁵ ÷ n, so a model to degree 360 omits 17.7 centimetres and one to 2190 omits 2.9. Every orthometric height derived from such a model carries that as an error, and it is not quoted with the height. What the numbers refer to

The geoid model stops at a degree

Eleven essays treat the geoid as a surface that exists. Every geoid anybody uses is a series truncated at a degree, so every orthometric height derived from one carries an omission error nobody quotes with the height — eighteen centimetres at degree 360 — and the same truncation removes two thirds of the slope, which does not converge at all.

 One meridian, two datums, and the ground between them. "the meridian line of the 141st degree of west longitude" — Anglo-Russian Convention, 1825. The line whose longitude is exactly 141° west, drawn twice: once on nad27 and once on WGS84, with the east–west separation exaggerated 3,000 times. The two are 129.0 metres apart, and remarkably constant — the shift changes by five centimetres over nine degrees of latitude — so the strip between them is a ribbon 129 metres wide and 1039 kilometres long, which is 134.0 square kilometres. The sentence has not changed. The surface the number refers to has. Grids, and what a survey does

A meridian boundary moves when its datum does

The 141st meridian is the one boundary description in this collection with no geometric ambiguity in it: every reading of it is the same curve, exactly. It has a different one. A longitude refers to a datum, the 1825 convention named none, and the line of longitude exactly 141° west sits 129.0 metres apart on NAD27 and WGS84 — a ribbon 1,039 kilometres long and 134 square kilometres in area.

 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.

 Seven cuts of the same size, in different places on one body. The body's own colatitude and longitude, with the seven windows drawn on it. Each covers the same surface area to 0.24 per cent — the longitude extent is divided by sin θ and a scale is then solved per window — so their sizes are held and only their positions differ. Each is shaded by the areal spread of the conformal map solved on it, listed beside the grid and running from 1.03 to 1.36. The darker windows are the worse ones, and they are the ones over the body's lobes. What the numbers refer to

Cuts of the same size in different places

Where a body is cut decides how well it can be mapped, by a factor of ten — established with five windows of five different sizes, so *where* and *how much* were confounded and the factor could have been entirely about extent. Held to the same surface area to a quarter of a per cent, the answer survives at a factor of 1.32, and what predicts it is the curvature the window encloses.

 A curve built to have dimension 1.2619. The generator replaces every segment with four of equal length at headings 0, +60.0°, −60.0° and 0. Closing the displacement fixes the length ratio at 0.33335, and four copies at that ratio give a dimension of exactly log 4 / log(2 + 2 cos θ) = 1.2619. Nothing here is measured yet: this is the construction the measurement will be checked against. Drawn at depth 5, which is 1024 segments, with the second-level shape shown faint beneath it. What a machine does with it

A line has a length only at a scale

A curve given by a formula can be sampled as finely as the picture needs. A map is not that: the geometry that reaches the page has been through an algorithm whose job is to throw most of it away. The first thing that goes is the idea that the line had a length.

 One of these settles. The largest departure of the fitted map's boundary scale from constant — the quantity the solver cannot see — against the number of collocation nodes, for the two placements. The clustered fit reaches 7.030e-5 at forty-eight nodes and returns exactly that at every count above it. The evenly spaced fit does not settle at all: it wanders by a factor of 1.43 across the same range, going up as often as down. Refining an evenly collocated fit is not convergence, and the finding that more samples improve the report and not the map is this seen from one side. What each projection optimises

The nodes were evenly spaced

Refining an evenly collocated fit improves the solver's report and leaves the map alone. Moving the same number of nodes to the Chebyshev positions — crowded toward the corners, where a conformal map of a polygon is singular — makes the fit settle: 7.030 × 10⁻⁵ at forty-eight nodes and exactly that at every count above it, against an even fit that wanders by 43 per cent and never converges at all.

 One quantity, one region, and the exponent left free. The scale departure of six projections over the world, aggregated as a p-norm, against p on a logarithmic axis. At p = 1 the best is Eckert IV; at p = 64 it is Winkel tripel. Nothing about the maps changed between the two ends of the axis — only how much of the region a bad point is allowed to spoil. Drawn in no projection: the axes are an exponent and a score. What is taught wrongly

The average was a choice of norm

Ten projections, one region, one measured quantity, and the only free decision left is how to turn a field into a number. Over the world's scale departure the ordering at the mean and the ordering at the worst case have a rank correlation of −0.04, all ten maps change position, and the exponent that produced each answer is stated nowhere.

 One equidistance line, computed five ways. The line equidistant from two facing coasts — Jan Mayen and Greenland — traced by bisection along a fan of parallels, with the two distances measured on the ellipsoid, on the sphere, and with a ruler on three different pages. The basepoints are the same four in every case and the rule is the same words in every case. The sphere sits 20 metres from the ellipsoid; with a ruler on Lambert's cylindrical sits 13.2 kilometres from it, with 324 km² of seabed in between. A delimitation is a sentence about distances, and a distance is a statement about a surface. Grids, and what a survey does

An equidistance line belongs to a surface

A maritime boundary is very often defined as the line equidistant from two coasts — a description with no coordinate in it and no curve to choose between. It has a third ambiguity: equidistant measured how. On the ellipsoid, on the sphere, and with a ruler on three different charts, the same four basepoints give lines up to 39.6 kilometres apart and 2,816 square kilometres of seabed between them.

 Three named members of the family, and one the family has no name for. The equal-area pseudocylindricals are not a list of maps. They are the solutions of one equation — C(φ)·Y′(φ) = cos φ — in two unknown functions, so one function is free and the named members are points in a space of them. Writing Y as six numbers and searching over them finds a map at 25.03° of mean angular deformation against Eckert IV's 28.08, Mollweide's 31.83 and the sinusoidal's 38.86 — 12.2 per cent better than the best map the family has a name for, and every one of the four is equal-area to the same precision. The families

A family is a function, not a list

The equal-area pseudocylindricals are the solutions of one equation in two unknown functions, so one function is free and the named members are points in a space of them. Writing that function as six numbers and searching over them finds a map at 25.03° of mean angular deformation against Eckert IV's 28.08 — and two numbers are already enough to beat every map the family has a name for.

 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.

All threads