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

The thread: Purpose before property — page 5

Asking which projection is best is asking an incomplete question. Every one of them minimises something, and the only useful comparison is between a projection and the job it was chosen for. Essays 97 to 120 of 153.
Believing a group takes away the redundancy that would test it. The thirteen degrees of freedom in this network, divided between the two groups, against the weight ratio the adjustment was told. At a ratio of a quarter the edges carry 12.1 of them and at sixteen they carry 0.07. A group's own variance can only be estimated from its own share, so a surveyor who is confident about an instrument has taken away the arithmetic that would have caught the confidence. The two curves sum to thirteen at every ratio, which is the identity that makes this a redistribution rather than a loss. Grids, and what a survey does

The weights are a guess the solve believes

Choosing what to hold fixed is a decision inside the least-squares problem that no residual can see. There is a second, made more often and thought about less. Every observation enters with a weight nobody measured, the weights move the coordinates by a factor of 1.8, and the standard check on them can be made to pass by a scaling that moves nothing at all.

 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. Measuring distortion

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.

 Where Robinson's valley breaks, against how large the region is. The threshold at which the set of near-optimal aspects stops being one connected piece, for square regions of growing size at 38° north. It falls from 2.10 at 6° to 0.27 at 30°, a factor of 7.9. The earlier measurement was made at one region and quoted "about twice the optimum"; that value belongs to a small region, and the prediction that it should fall as the region grows is what this tests. What each projection optimises

Where the valley breaks in two

A near-optimal set is one connected sheet at a loose threshold and thirteen basins at a tight one, and the transition was explained without being tested. The explanation is a prediction about region size and projection sharpness: swept over both, the threshold falls from 2.10 to 0.27 as a region grows from 6° to 30°, three projections lie on nearly one curve, and a fourth declines to join for a reason worth having.

 Hilbert order on one face, as a curve. The order in which Hilbert numbering visits the 64 cells of one cube face at level 3. The line never leaves a cell without entering one that shares an edge with it — that is what makes it a space-filling curve, and it is why two cells with nearby identifiers are usually near each other on the ground. What a machine does with it

The address is a curve through the sphere

A database does not fetch a set of cells, it reads ranges of identifiers — so the cost of a query is how many runs its cells form, not how many cells it needs. Hilbert order wins that measurement and loses the one usually quoted for it: its neighbouring cells are further apart in identifier than row-major's, on average and at worst.

 The same kernels, across an edge. Each kernel measured over the same five rasters as the smooth measurement, on a field that carries a step across a tilted line instead of a smooth function. The solid lines are the edge and the faint ones the smooth field. The orders on the smooth field are nearest 1.00, bilinear 1.98, cubic 2.93; across the edge they are nearest 0.78, bilinear 0.58, cubic 0.60 — within a factor of 1.33 of one another, and the fastest of them belongs to nearest-neighbour, which does no interpolating at all. The ranking a smooth field establishes does not survive a discontinuity, and a real raster is mostly edges. What a machine does with it

An edge has no order of convergence

On a smooth field the three resampling kernels converge at orders 1, 2 and 3 and the choice is obvious. Across a discontinuity they converge at 0.78, 0.58 and 0.60 — within a factor of 1.4 of each other, in an order that puts nearest-neighbour first, and a real raster is mostly edges.

 How large a blunder has to be before the test notices. A blunder of increasing size put into the least-checked observation of a braced quadrilateral, with the standardised residual it produces. The horizontal line is the critical value the test uses, and the vertical one is the minimal detectable bias — δ₀σ/√r, which is 87 millimetres for this observation and is computed from the network's DESIGN, before any observation is made. Below it nothing is flagged; above it everything is. The observation's redundancy number is 0.144, so it is checked by a seventh of an observation and hides six-sevenths of whatever is wrong with it. Grids, and what a survey does

The blunder the network cannot see

A least-squares adjustment has no concept of a mistake. The smallest blunder its test will find in the least-checked leg of a braced quadrilateral is 87 millimetres, and by the time it fires a station has moved by nearly ten times the accuracy the same adjustment reports for it.

 The plan, the crossing that was flown, and the one a right forecast would have given. A crossing of NaN km planned against a jet forecast to sit on the fortieth parallel, flown through a jet that is actually 6 degrees further south, and re-planned six times on the way. Committing to the plan costs 167.81 hours. Re-planning costs 165.84. A vehicle that had known where the jet was would have taken 163.21. So the wrong forecast is worth 4.60 hours and re-planning recovers 1.97 of them. Paths and directions

A crossing is a chain of decisions

A route is handed back as a curve, and nothing anybody flies is decided once. Re-planning is worth exactly nothing when the forecast turns out right — a sub-path of an optimal path is optimal, and the chain is the plan. When it turns out wrong it recovers 19 per cent of the cost at two decisions and 83 at eighteen, and never all of it.

 Every projection against a picture that is not a map. The worst relative distance error over London, New York, Tokyo, Sydney, for each projection in the library and for the best flat arrangement of the same distances. Each azimuthal member is centred on the set's own centroid, which is the fair comparison. The free picture reaches 0.63% and the best projection, Azimuthal equidistant, reaches 8.28% — a ratio of 13.17. The free picture cannot lose, because every projection's own layout was handed to the search as a starting point; what the figure measures is how much the freedom is worth, and it is worth different amounts at different sizes. The impossibility

The best flat picture is not a map

A set of dots whose separations are as nearly right as separations can be made beats every projection in this collection — by a factor of thirteen on four world cities, and by eight per cent on sixteen. The collapse between those two numbers is not about cartography. It is that a picture of n places has 2n − 3 free numbers and n(n − 1)/2 distances to spend them on.

 Four layouts, the same five stations, the same ten distances. Every panel has five stations, all ten distances between them, the same instrument precision and the same three degrees of freedom. The ellipses are the error ellipses of the adjusted coordinates, drawn at one common exaggeration, and they are computed from the geometry and the weights alone — no observation value enters any of them. The worst semi-axis runs from 11.1 millimetres to 232, a factor of 20.9, and the difference is entirely where the marks were put. Grids, and what a survey does

The network's answer is decided before it is measured

An adjustment is usually judged by what it does with observations. Every quantity a specification is written about — the error ellipses, the redundancy numbers, the smallest detectable blunder — is a function of the geometry and the weights alone, and does not contain an observed value anywhere. Four layouts of five stations with the same ten distances differ by a factor of 20.9 in their worst coordinate.

 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. Measuring distortion

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.

 One pair of numbers, a hundred and twenty places. The easting 412,000 and northing 5,678,000, interpreted in each of the sixty zones and in both hemispheres. The northern candidates are a ring of sixty at 51.247° north, spaced exactly six degrees apart; the southern ones are a second ring at -39.043°, which is not the mirror of the first because the southern convention subtracts the northing from ten million. Every one of the hundred and twenty is a perfectly valid reading of the same two numbers. Grids, and what a survey does

One pair of numbers, a hundred and twenty places

A UTM coordinate is two numbers and a zone. Drop the zone and the numbers are still valid in each of the sixty; drop the hemisphere too and the pair names a hundred and twenty places. They form two rings at one latitude each, spaced exactly six degrees apart, and every one of them has the same grid convergence and the same scale factor — so no further geometric measurement can choose between them.

 Four of the seven on the front can never be first. Every library projection over the whole sphere, with both errors normalised to the table's own range. The seven filled circles are the Pareto front — nothing beats them on both counts. The line through three of them is the lower convex hull, and a weighted sum of the two errors is a straight line in this space, so only a hull vertex can ever come first. Four projections — Web Mercator, Miller cylindrical, Equirectangular, Winkel tripel — sit in the dents, undominated and unchoosable. What is taught wrongly

Which projection a weighting can make best

The seven world projections nothing beats on both counts tell a reader with a preference that one of them is their answer. Four of the seven are not: they sit in dents of the front, undominated and unreachable, and no weighting of angle against area can ever put them first. Robinson can be first, on 2.8 per cent of the weight range.

 The shortest crossing that never goes more than 745 km from an airfield. New York to London, drawn on a conformal conic, with a circle of 745 km round each of twenty airfields on both sides of the North Atlantic. The dashed line is the great circle, 5,540 km. The solid line is the shortest route that stays inside the circles, 5,698 km — 158 km or 2.85 per cent longer — and it bends only where two circles cross: between Narsarsuaq and Keflavík; Keflavík and Vágar. Nowhere else can it bend, because anywhere else a straight cut stays inside a circle. Paths and directions

The rule that keeps a route near land is pinned at both ends

A crossing that may never be more than a stated distance from an airfield is refused below one radius and free above another, and each wall is set by a single feature of the ocean. Between them the price falls continuously, with corners and no jumps — and at the lower wall it falls as the square root of the allowance, so the first kilometre is worth twelve of route.

 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.

 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.

 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.

 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.

 The same deflection, from a compact mass and from a broad root. Each curve is the deflection of the vertical across a mass buried 8 km down, with the mass solved so that all of them peak at 10″. The narrow one is the buried sphere the earlier measurement used; the broad ones are crustal roots 40 and 160 km wide, modelled as that mass spread along a line. They agree where it matters most and disagree everywhere else: the signal is 30 km wide for the compact body and 188 km for the widest, which is the difference the plumb line feels as it descends. What the numbers refer to

A mountain is not a buried sphere

The plumb line's drift was measured over a compact buried body and grows as the 0.69 power of the column's height — an exponent that is a statement about how quickly a buried sphere's field weakens with distance rather than about mountains. Spread the same mass into a crustal root and the exponent climbs to 0.84, while the proportionality to the deflection survives exactly.

 Two regions, three answers. Britain and New Zealand, 166° apart, with the pole of the best oblique conic under each of three objectives. Pooling the samples and taking an area-weighted score puts the pole in one place; refusing to let either region be worse than the other puts it somewhere else. The regions are drawn on Mollweide so that equal ground areas are equal page areas. What each projection optimises

The pooled score abandons a region

Every criterion in this subject optimises for one region. An atlas is several, and pooling their samples into one area-weighted score is what everybody does — which on Britain and New Zealand serves Britain 1.2 times worse than it could be served alone and New Zealand 125 times worse. The worst-case objective makes them equal at 33 and 59, and the cost of sharing rises with separation from 1.4 to 59.

 The test the conjecture asked for, and it refutes the conjecture. How much better the best asymmetric projection is than the best symmetric one, under weightings of four different symmetries, with every symmetric map allowed to re-aim its axis at twelve candidate poles. The conjecture was that the seven earn their place by PLACING distortion where a symmetric map cannot, so their advantage should collapse under a criterion with no place preference. It does the opposite: the advantage is largest at 1.343 under the uniform weighting and smallest at 1.144 under a band, with the fully asymmetric concentration at 1.204 in between. The winner is named on each row and the map it beat is Equirectangular throughout. The seven are simply better maps. The families

The maps with no family are simply better

Seven projections have no continuous symmetry, and they are almost exactly the set anybody would choose for a world map. That coincidence came with a conjecture and a test attached: their advantage should collapse under a criterion that does not care where anything is. Run, it does the opposite — 1.343 times under a uniform weighting and 1.204 under a concentration. The conjecture is refuted.

 A meander grows, and its boundary moves with it. Three loops of a river whose valley repeats every 5 km, drawn when its bends are gentle — 30° at most from the valley's direction — and again at the moment its necks close to the channel's own 150 m width, at 120.3°. Under accretion the boundary is the deep-water line at every moment, so it moves with the river, and the shaded ground is the ground that has changed side of it between the two moments: 83.5 km², 27.8 km² a loop. The river is a stated sine-generated curve, not a surveyed one. Grids, and what a survey does

A river boundary goes where the river goes, or stays where it was

A boundary that follows a river moves when the river moves, and the law of every such boundary turns on one word: whether the river moved gradually or suddenly. Over one meander's life, gradual movement changes the country of about twenty-eight square kilometres of land a loop. A single flood that cuts the loop off strands twenty-six — nearly the same land, arriving as one parcel on the wrong side of the river.

All threads