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

The thread: Purpose before property — page 6

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 121 to 135 of 135.
Twenty-four versions of one shape, and not one of them gains area. The same closed boundary rotated twenty-four times and simplified at the same tolerance. If the area error were noise the values would straddle zero and their mean would fall towards it; they do not. Every one is negative, the mean is -0.4644 per cent, and the mean is 71 standard errors from zero. A bias of that size cannot be removed by averaging over more boundaries, which is the only defence anybody has against a rounding error. What a machine does with it

A thousand features are wrong in the same direction

The area a simplification costs is unpredictable in sign for one feature. Over a population it is not: twenty-four presentations of one shape all lose area, the mean is seventy standard errors below zero, and no amount of aggregation removes it.

 How long a network has before it breaks its own tolerance. The years until the worst baseline in a 400-kilometre network exceeds five millimetres plus one part per million of its length, which is an ordinary first-order specification. A boundary zone breaks it in 16 years; a plate interior takes 1819; a rigid plate never does, at any speed, because a rigid body keeps every distance it has. Bars are clipped at five thousand years. What the numbers refer to

A grid stops fitting the ground it was laid on

A hundred-kilometre baseline across a plate boundary changes by 815 millimetres in fifty years, and the same baseline on the fastest plate in the model changes by nothing at all in the same fifty. The interval before an ordinary first-order specification is broken is six years in one place and never in the other.

 One of these three lines does not slope. Three lengths at each zoom, for a world coordinate near the antimeridian. The pixel halves with every level, as it must. The spacing between representable double-precision values is far below it and halves with it. The spacing between representable SINGLE-precision values does not move at all — 2.0 metres at every zoom, because it is a property of the size of the number and the world does not get smaller. The two cross at zoom 17, and past it a vertex snaps to a lattice coarser than the pixels it is drawn into. What a machine does with it

The renderer runs out of numbers before the zoom does

Eleven essays price the pyramid in exact arithmetic. The pipeline that draws it carries single precision, where a world coordinate near the antimeridian quantises to two metres — at every zoom, because the number does not get smaller when the pixel does. At zoom 22 that is fifty-four pixels.

 The query, the box it becomes, and the 14 per cent it loses. A query for everything within 400 kilometres of a place at 55° north, drawn in the stored coordinates of Web Mercator. The ring is the true answer's edge; the rectangle is the box the index is given, sized from the scale factor at the query point. The arcs outside the rectangle are true neighbours the index never returns — 13.6 per cent of the rim, reaching 5.0 per cent of the box's own half-width beyond it — and nothing downstream can tell they are missing. What a machine does with it

The query a fast path actually answers

Fifteen rungs price what a stored coordinate means and none asks what it is searched with. No index answers "within two hundred kilometres of here"; an index answers "inside this rectangle of stored coordinates", and the rectangle is built by somebody's arithmetic — which at 55° north silently drops a tenth of the true answer on a conformal projection and three fifths of it on an equal-area one.

 The sea surface, against the surface heights are measured from. The stated mean dynamic topography — how far the sea stands above the geoid — drawn on Mollweide, with the eight tide gauges marked. It ranges over 2.09 metres, it is smooth, and it is what every national datum's zero is sitting on. Nothing here is a model's output: it is a stated closed form, chosen to have the observed sign structure and the observed size. What the numbers refer to

Every country's zero is a different surface

Ten rungs measure a height against a geoid and treat the geoid as one object. No national datum is on it: each is pinned to the mean sea level at one tide gauge, the sea surface stands up to two metres from the geoid, and the eight European zeros measured here spread over 462 millimetres — a step no levelling can remove and nobody's error.

 A ruler on a cartogram is not measuring anything. Twelve pairs of places, each measured on the ground and on the page of the triangular, x first cartogram of four cities, with the page scaled so that the median pair reads exactly right — the most generous calibration available. The line is that calibration. The worst pair, Tokyo to Quito, reads 21598 km for a ground distance of 14439 km, an error of 49.6%. The open marks are the same pairs on the equal-area map the cartogram was built from. Measuring distortion

Everything else on the page pays for the areas

A cartogram gets one quantity exactly right and every other reading a page supports is collateral. A ruler on it is out by 35 per cent after the most generous calibration available, and ten of sixty triples of places change which one is in the middle.

 One rule, applied once, to one dataset. Two features 3 metres apart on the ground, and a snapping tolerance of 5 map units on Web Mercator. The tolerance's reach is the curve; the pair's separation is the line. They merge everywhere below about 52.8° and stay separate above it — the same features, the same rule, the same run. A pipeline that cleans a global dataset in one pass produces a topology that changes at a latitude nobody chose and nothing records. What a machine does with it

A tolerance in map units is not a tolerance

A snapping tolerance is a number, and the number is in whatever units the file is in. Five map units on Web Mercator is 4.97 metres of ground at the equator and 0.87 at eighty degrees — so a rule that merges two features three metres apart merges them everywhere below 52.8° north and refuses everywhere above it, in one pass, over one dataset, with nothing recording where the boundary is.

 The density every real cartogram is handed. Two islands in an ocean that is not lightly populated but empty: the density is exactly zero over 83.7 per cent of the sphere. Four rungs of this anchor assume a positive density, because the construction divides by it — and the conditional cumulative of a column with no mass in it is zero over zero. Drawn on the equal-area base, so a cell's ink is a density and its area is ground. Measuring distortion

A density that asks for no room at all

Five rungs assume the density is positive everywhere, because the construction divides by it. Every cartogram anybody draws has an ocean, and an ocean is not sparsely populated but empty — 83.7 per cent of the sphere, exactly zero, and the construction returns nothing at all for a third of the probes.

 The same tiling, turned. The cube's eight vertices, ringed, and eight cities, with the aspect of the cell each city falls in written beside it. The two panels are the same tiling: the same cells, the same areas, the same shapes, in the same numbers. Only where they sit has changed, and the eight cities get cells whose areas spread by 1.26 in one and 1.39 in the other. Drawn on Mollweide, in which equal ground areas are equal page areas. What a machine does with it

The orientation is a policy

A polyhedral cell system has three free angles nobody scores. They cannot improve it: rotating the solid rotates every cell rigidly, so the distribution of cell areas is identical for every orientation there is. What they decide is who stands on the bad cells — and the eight cities measured here get a spread of cell area of 1.00 under the best turn and 1.50 under the worst.

 Waiting longer buys almost nothing under the spectrum that is actually there. The scatter of block means against the block length, for the three spectra. White noise falls with a slope of -0.484 — the inverse square root everybody assumes. Flicker falls with a slope of -0.111, so a hundredfold longer average buys a factor of 1.67 rather than ten. A random walk is very nearly flat. Measuring distortion

The error that does not average down

Five rungs give a coordinate a width, and all five assume the observations it was averaged from are independent. They are not. Under the noise spectrum every published analysis of a position time series reports, the standard error falls as the eleventh root rather than the square root — so a factor of ten costs a hundred observations if the noise is white and a thousand million if it is not.

 The same field, the same statistic, two pages. A single northern concentration — v = 10 + 90 exp(−d²/(25°)²), d the distance from 20°E 55°N — drawn twice, with the shading rule and the data identical. The ground's own area-weighted mean is 14.009. A reader weighting by the area actually on the page takes 17.490 off Mercator, which is +24.85%, and 14.009 off Gall–Peters, which is +0.00%. Neither map has misdrawn a single cell. Measuring distortion

A choropleth is read by area

Every cartography course states the rule — use an equal-area projection for a thematic map — and states it as advice. It is a theorem, and it has a residual: the error a page puts into a reading is exactly the covariance of the value with the areal factor, which is 24.85 per cent for a northern concentration read off Mercator and 0.00 per cent for the same field read off a map that spreads area by 7.7 to one.

 The same totals, drawn as symbols, on two pages. A band field, largest in the mid latitudes carried by each region and drawn as a circle whose area is proportional to the total. Every symbol is right: a total is a total wherever it is drawn, and the two maps carry identical numbers. What differs is the ground under each symbol. On Mercator the symbol at 60° sits on a region drawn 4.0 times larger than it would be on an equal-area sheet, so the density a reader forms — symbol against region — is out by that factor. Measuring distortion

A symbol has a size on the page and an area on the ground

A proportional symbol is right as a total wherever it is drawn — a count is a count. Read against the region beneath it, which is how a reader forms a density, it is wrong by exactly the reciprocal of the areal factor: 0.083 at 73° north on Mercator, a factor of twelve, with the correction available as one multiplication that no atlas makes.

 One field, one classifier, two sets of breaks. A field that varies with latitude put into 5 classes by area-weighted quantiles, on Mercator. The upper panel weights each region by its ground area and the lower by the area it occupies on this page, which is what a classifier handed projected geometry does. 60 of 150 regions land in a different class, marked in the lower panel. The data has not changed and neither has the number of classes. Measuring distortion

The class breaks were computed on the page

The three rungs below price what a reader does with a finished map. A classifier is software, it runs on the geometry it has, and the geometry it has is projected: a five-class quantile classification of one stated field puts half of the three hundred and eighty-four regions in a different colour on Mercator, and 87.5 per cent of them at nine classes.

 Six regions, before and after. Six circular regions of unequal size on the equal-area rectangle, and the same six after the cartogram of four cities has been solved. Each drawn area is exactly its base area times the density it was asked for. The four that grow stay recognisably round; the two that shrink are drawn out into shapes that no longer resemble what they were, and nothing in the construction chose to treat them differently. Measuring distortion

A cartogram keeps the shapes it inflates

Seven rungs build cartograms and none reads one back. Reading means recognising a region and dividing its drawn area by its base area, and the construction is against the reader twice: the correlation between how much a region grows and how much of its shape it keeps is −0.996, and the base area a reader has to divide by is the map the cartogram replaced.

 Two places, one written coordinate. The lattice a coordinate written to 5 decimal places lives on, at 45° of latitude, where a cell is 1.11 metres north to south. Two real places a fifth of a cell apart are drawn as open marks and the single value they both round to as a filled one. The distance between them, computed from what was written, is exactly zero — an error a noise model cannot produce, since independent noise never puts two different points in the same place. Measuring distortion

Rounding is not noise

Eight rungs treat a coordinate's error as noise that averages down. A published coordinate has a second error that does not: it is deterministic, it is shared between every point in the same cell, and at five centimetres apart ninety-three per cent of pairs come out as the same place — which no amount of independent noise can produce.

All threads