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

The thread: What survives a change of coordinates — page 5

Scale along the meridian and scale along the parallel depend on how the sphere was parameterised. The principal scales, the areal factor and the angular deformation do not, and only those are properties of the map itself. Essays 97 to 115 of 115.
How aligned a region's ellipses are, for every projection and every region. The resultant length of the doubled indicatrix orientations, over 12 projections and 9 regions. One means every ellipse in the region points the same way; zero means they are spread evenly and cancel. Three cylindrical rows are 1.00 throughout, and everything else varies down the row and across it — which is the answer to the question the number was first stated without: alignment is a property of the pair. Mercator's row is the exception that is not a measurement, because a conformal projection's indicatrix is a circle and a circle has no orientation. Measuring distortion

Whether the ellipses point the same way

Alignment was found to decide whether the average of a region's deformations is above or below the deformation of its average, and stated it at ten projections over one region. Swept over a hundred and eight pairs, the answer is that alignment belongs to the projection 28 per cent, to the region 38, and to neither 33 — and two rows of the table turn out not to be measurements at all.

 A territory, and two enclosures somebody might score instead. A stated territory at 45° north — an ellipse twenty-eight degrees long and eight wide, its long axis bearing 45° — with its latitude–longitude bounding box and the smallest cap around it, on a globe. The territory fills 40.6 per cent of the box and 28.6 per cent of the cap. Scoring the library over each names a winner: an azimuthal for the territory itself, an azimuthal for its bounding box, an azimuthal for the smallest cap around it. The cost of working from the box rather than the territory, measured as the box's winner scored on the territory against the territory's own winner, is a factor of 1.0584. What is taught wrongly

The region was a rectangle in a projection

Every verdict this collection hands out is a verdict about a region, and almost every region it has been given is a latitude–longitude box. A box discards a territory's shape, which costs a few per cent, and it records the territory's width in degrees of longitude, which costs up to a factor of 1.72 and changes which family wins on four of forty stated territories. Convert that width to arc and the box gives the right answer every time — so the comparison had a projection in it before it began.

 A ground that is not deforming, read off a Mercator sheet. Every place on this map is carried by a rigid rotation of the whole Earth at forty millimetres a year — the motion that deforms nothing, and that this collection has already shown deforms nothing. The circles are the strain rate a geodesist would report from the grid coordinates alone, up to 17.4 nanostrain/yr. None of it is on the ground. It is the projection's own scale factor changing along the displacement, which is a second derivative of the map arriving in a first-order measurement. What the numbers refer to

The strain a map adds to the ground's

A ground carried rigidly at forty millimetres a year deforms nothing, and read off a Mercator sheet at sixty degrees north it reports 10.9 nanostrain a year. Along a profile through a real boundary the invented part is 0.0 per cent of the answer where the zone is loud and 884 per cent where it is quiet.

 A compacting basin, and the tilt it produces. The stated vertical velocity field — a bowl of subsidence 120 km across, with no mass leaving — drawn as circles proportional to the rate, with the tilt of the ground surface as the arrows. The subsidence is largest at the centre and the tilt is exactly zero there, because a smooth bowl has no gradient at its own bottom. The largest tilt is 125.5 nanoradians a year, on a ring at the bowl's own scale length over root two, and it is the quantity a levelling network measures. What the numbers refer to

A vertical rate needs a height system

Four rungs of this anchor measure the two-by-two horizontal tensor, because that is what a tangent chart returns. The larger signal in a subsiding basin is vertical — 126 nanoradians a year of tilt against 0.76 nanostrain a year of horizontal strain — and it is not a measurement at all until the surface it is measured against is named, because that surface is moving too.

 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.

 Seven published numbers, twenty-eight actual ones. The correlation matrix of a seven-parameter fit to 64 common points over a region 9° across. Only the diagonal is ever published — seven standard deviations — and the twenty-one off-diagonal entries are not small: the strongest is ty against rx at 0.940. A translation and the rotation that mimics it over a small patch are very nearly the same parameter, so the fit cannot tell them apart and its errors in the two are locked together. What the numbers refer to

The parameters are not independent

Rung seven gives the seven parameters their own uncertainty and stops at seven numbers. There are twenty-eight, and the twenty-one nobody publishes are not small: a translation and the rotation that mimics it correlate at 0.94, the normal matrix has a condition number of 4 × 10¹⁶, and propagating from the diagonal alone overstates the transformation's uncertainty by up to a factor of thirty-six.

 The coast copied twelve miles out gives up the sea of every narrow bay. A stated coast 116 km long with five bays 5 km deep, 6, 10, 16, 24, 36 km wide. Dashed: the envelope of arcs, the line every point of which is twelve nautical miles from the nearest point of the coast. Solid: the coast itself moved twelve miles seaward, which dents the limit by the full 5 km over every bay. Tinted: the sea between them, 240 km² in all, 9% of a straight coast's territorial sea. Over the 6 km bay the copy lies 4.8 km inside the envelope; over the 36 km bay the two meet. The coast is stated, not surveyed. Grids, and what a survey does

Copied out to sea, a coast keeps its length and loses its bays

Before 1930 a territorial sea was drawn by following the coast in all its sinuosities, and the phrase fits two constructions. Moved twelve miles seaward, a stated coast of five bays gives up 240 square kilometres against the envelope of arcs, most of it over the narrowest bay; moved along its own normals, it ties a loop over every bay. On a fractal coast charted at finer detail the copy grows by a quarter at every level, to 480 kilometres, while the envelope settles at 133.0.

 Conformality helps and does not save it. The proportion of a stated terrain whose plan curvature changes sign when it is read off a grid in each projection. A conformal map turns every direction through the same angle, so the contour and the slope line stay perpendicular and the sign ought to survive — and it mostly does, at 0.30 per cent against 13.3. It is not zero, and the term that flips it is the gradient of the scale factor: the curvature of a curve under a conformal map is (κ − ∂ₙ log λ)/λ, and Mercator's λ has a gradient. Measuring distortion

The curvature of a field is not the curvature of its picture

Whether a place is a spur or a hollow is the sign of a second derivative, and every automatic terrain classification is built on it. Read off an equal-area grid, that sign is wrong on 13.25 per cent of a sheet — and conformality reduces it to 0.30 per cent without removing it.

 A surface whose curvature changes sign, and integrates to nothing. a wide ring of major radius 3 and minor radius 1, shaded by its Gaussian curvature. The outer half is positively curved like a sphere, at up to 0.250; the inner half is saddle-shaped and negative, down to -0.500; and the two circles between them, drawn as lines, are exactly flat. The integral of the curvature over the whole surface is zero, which is 2π times the Euler characteristic — and every impossibility in this collection rests on that number being 2 rather than 0. What the numbers refer to

On a body with a hole, north can be up everywhere

Varying a body's shape never changes its topology, which is what every impossibility about maps actually rests on. A torus has a nowhere-zero tangent field, a total curvature of zero rather than 4π, and a conformal world map with no cut and no singular point — and it still cannot be flattened, for the one reason that survives.

 Millimetres of horizontal move per kilometre of assumed height. How far the transformed latitude and longitude move when the height fed into a datum transformation changes by one kilometre, for four published transformations at five places. It runs from 7.7 to 121 millimetres per kilometre. A latitude and a longitude do not carry a height, so a two-dimensional coordinate cannot be transformed until somebody supplies one that is not in it. What the numbers refer to

The height a coordinate does not carry

Ten rungs treat a datum shift as a map from one pair of angles to another. It is not one: the transformation runs through Cartesian coordinates, so the answer depends on the height — by 7.7 to 121 millimetres per kilometre depending on the datum and the place, which puts Lhasa 432 millimetres from where a height of zero would have said.

 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.

 The points a projection cannot move. The critical points of a stated field, found twice: once on the sphere from its own gradient, and once in Lambert cylindrical's page coordinates from the page's own numbers, with nothing shared between the two searches. seven points, 3 maxima, 2 saddles and 2 minima, and the two sets agree in position to 2.3e-7 degrees and in type at every one. The page's Jacobian is invertible wherever the map is a map, so it sends a zero gradient to a zero gradient and cannot move a critical point anywhere. Measuring distortion

What the page cannot move

Six rungs measure what a page does to a field's readings and every one of them moves. The critical points do not: found on the sphere and found again in a projection's own page coordinates with nothing shared between the searches, they agree to 10⁻⁷ degrees and in type at every point, on every projection. What the page does move is their shape, by a factor bounded exactly by the indicatrix's axis ratio squared.

 One of these averages exists. Mercator's area-weighted mean areal factor, and its Kavrayskiy number, against how close the sampled band comes to the pole. The first is artanh(sin Φ)/sin Φ in closed form and has no limit: it passes 7.04 at a tenth of a degree from the pole and keeps going, gaining a fixed amount every time the remaining gap is halved. The second settles by 85° and does not move again. Both are published as summary distortion figures for the same map. Measuring distortion

A mean that does not exist can still be printed

Mercator's area-weighted mean areal factor is artanh(sin Φ)/sin Φ, and it has no limit. A sampler asked for it returns the logarithm of its own sample count plus 1.512 — measured slope 1.001 against ln n — so the number is a property of the person who computed it. The Kavrayskiy number for the same map over the same sphere settles at 0.52124 and is a number.

 Which of this collection's numbers are the sampler's. How much each of two published quantities moves when the sampler behind it goes from twenty samples a side to sixty, for eight projections over the whole sphere. The Kavrayskiy numbers — the summary means the rankings are built from — move by at most 1.15 per cent. The worst-point angular deformations move by up to 19.0 per cent, all in the same direction, because they are maxima over a sample and a maximum over a sample is a lower bound. The split is clean, and it says which numbers here need refining and which do not. Measuring distortion

Which of these numbers are the sampler's

Four rungs have shown that a sampled maximum understates, a sampled mean can be a report on the sampler, no arrangement of points is neutral, and refining until the answer settles proves nothing. So the collection re-measured itself. The means move by at most 1.1 per cent, the worst points by up to 19, and the rankings — which is what the essays actually argue with — do not move at all.

 Two special lines, and they are not the same line. Each projection's true-scale parallels, found by solving for a parallel scale of exactly one, and the parallels along which its second-order failure vanishes, found by minimising the flexion. Five of eleven put them more than five degrees apart, and Gall–Peters puts them forty-five. The pattern is tangency: a projection with one standard parallel has both lines there, and a secant one has its true-scale lines moved off the centre while the bending zero stays. Measuring distortion

The lines where the bending vanishes

Every projection has a line printed in its margin — the standard parallel, where the scale is exactly one. It has a second special line nobody prints: the one along which a geodesic is drawn straight to second order. On Gall–Peters they are forty-five degrees apart, and the rule turns out to be tangency — a tangent construction puts both lines at its point of contact and a secant one moves only the first.

 Densifying helps one and hurts the other. The standard deviation of the shoelace area, and the bias of the perimeter, against the number of vertices on the same 100-metre circle with the same 5-centimetre noise on each. The area's spread falls from 9.9 m² to 2.0 — the closed form says it goes as the square root of the vertex count's reciprocal, because a vertex's influence on the area is the vector between its two neighbours and densifying shortens it. The perimeter's bias rises from 0.07 mm to 261, a factor of 3697, because every leg contributes its own σ²/d and shorter legs contribute more. Measuring distortion

The area is unbiased and the perimeter is not

A boundary measured from noisy vertices comes out long, always, by σ²/d on every leg. The area enclosed by the same vertices comes out exactly right, because a shoelace is bilinear and the cross terms vanish. So densifying a boundary makes its area five times more precise and its perimeter three thousand times more wrong, and every compactness score computed from it falls short.

 Two maps with the same ellipses. Mollweide, and the same projection with the sheet turned 37° on the desk. Every ellipse is the same size and the same shape as its counterpart, and its angle to its own graticule is the same — the two indicatrix fields agree to 1.8e-12, which is arithmetic noise. The rotation differs by 37° at every point. So the classical description of a map's distortion is exactly blind to the difference between these two, and the difference is the whole of what a surveyor calls convergence. Measuring distortion

The fourth number the ellipse does not carry

A Jacobian has four independent entries and an indicatrix reports three. The missing one is the rigid rotation in the polar decomposition A = R·S, the indicatrix is exactly S, and R is what turns north into grid north: for the transverse Mercator it agrees with the survey formula for convergence to 5 × 10⁻¹⁰ degrees, from a different library and a different derivative.

 What a reprojection does, drawn where it does it. The angular deformation of the map that takes Sinusoidal's page to Gall–Peters's — which is what a reprojection is — shaded over the sphere and drawn on Mollweide so that equal ground areas are equal page areas. Both parents are equal-area, so the composite preserves area exactly; its mean angular deformation is 55.5° against 38.9° and 32.2° for the two maps separately. Composing two maps did not average their shape errors; it made a larger one. Measuring distortion

Two indicatrices do not make a third

Reprojecting is composing, and the composite's indicatrix is not a function of its parents'. Mollweide followed by Hammer is gentler than either map over eighty per cent of the sphere; the sinusoidal followed by Gall–Peters is worse than both everywhere. What separates them is one angle, and it is the number rung ten showed the ellipse does not carry.

 One current, two maps, and only one of them conserves it. A flow that has no sources or sinks anywhere on the sphere — it is the perpendicular of a stream function's gradient, so its divergence is zero by construction — drawn on two projections. The arrows on the right are the images of the same ground velocities as the ones on the left. On the equal-area map the drawn field is still divergence-free; on the other it is not, and the arrows drawn in the warning colour are where a reader measuring the picture would find a source or a sink that is not there. Measuring distortion

A current drawn on a page has sources

Seven rungs project a scalar field and ask what the page does to its gradient. A wind or a current is the other half of what gets mapped, and the operator that matters for it is the divergence — which is preserved by an equal-area map exactly, by no other map at all, and by a conformal map least of anybody's expectation.

All threads