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

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

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 73 to 94 of 94.
There is no such thing as the radius of the Earth at a latitude. The two principal radii of curvature of WGS84, against latitude, with the mean radius that every table prints drawn across them. The meridional radius M runs from 6335.44 km at the equator to 6399.59 at the pole; the prime-vertical radius N from 6378.14 to the same value. They differ by 42.70 km at the equator and meet only at the pole, and the single number 6,371 km lies between them at no latitude where either is right. What the numbers refer to

The radius of curvature is two numbers

Twelve essays work on the ellipsoid and every one of them takes a radius when it needs one. There are two at every point, they differ by 42.70 kilometres at the equator, and the single number every table prints is out by 5,583 parts per million on a line running north.

 a transform boundary: what the ground is doing to itself. The same 72 places as the velocity field, with the strain rate computed from the motion's own four partial derivatives rather than from its size. Each cross carries two principal rates: the long stroke is the greater extension, the barred one is shortening. The largest second invariant anywhere here is 241.0 nanostrain/yr, and the arithmetic is the same arithmetic that reads a projection's indicatrix. Drawn in Azimuthal equidistant centred on the window. What the numbers refer to

The ground has an indicatrix too

Two hundred and forty-six essays hold the Earth still while the page is measured. The ground is moving at tens of millimetres a year, and the motion is a map with four partial derivatives — so the same construction that draws Tissot's ellipse draws one for the ground, and the fastest plate in the model returns exactly nothing.

 One velocity field, three frames. a transform boundary, drawn three times. The ground is the same ground and the arrows are not the same arrows: the fastest velocity anywhere in the 3 panels is 70.6 mm/yr and in the quietest panel the same places are nearly still. A frame is a choice of which rotation to subtract, and there is no measurement that picks one. Every panel has the same strain rate at every point. What the numbers refer to

A velocity needs a frame and a strain rate does not

One place on a plate boundary moves at 5.3, 7.0, 10.2 or 54.6 millimetres a year, and towards the south-east or the north-west, depending on which rotation was subtracted first. Every one of those readings reports the same strain rate, to two parts in a hundred thousand million.

 A 900 km circular accuracy at 55° north, projected. Six thousand ground positions drawn from a circular error of 900 kilometres about one place, each projected in Mercator and plotted as a displacement from the projected place. The curve is the nominal 95 per cent ellipse, computed the standard way — the ground covariance sandwiched between the projection's own derivatives. It holds 93.83 per cent of the points, the cloud is measurably longer than it along its own long axis by 5.52 per cent, and it is not symmetric: the third moment along the page's second axis is 0.727 rather than zero. Measuring distortion

The error ellipse is not an ellipse

Rung two pushed a covariance through a projection with the same matrix sandwich that draws an indicatrix. That is a first-order operation on a map with a second derivative, so the propagated distribution is not the ellipse the sandwich draws — and a nominal 95 per cent ellipse holds 93.06 per cent on one projection and 95.63 on another, in opposite directions, from the same input.

 The same resolution, the grid moved, and a different answer. The same field at a fixed 18 × 9 division, with the grid slid by fractions of a cell. Nothing is lost — every cell is the same size as before and there are exactly as many of them — and the largest reported value moves over a range of 12.59 per cent. That is more than a whole halving of the resolution costs, which is 10.15 per cent on the same field. The scale effect has an excuse and this one has none. What a machine does with it

The answer depends on the cells it was counted in

Nine essays price the cell as a shape. The number reported out of it is priced nowhere: sliding a grid without changing its resolution moves the largest reported value by 12.6 per cent, which is more than halving the resolution costs.

 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

The previous rung found that alignment decides 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 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.

 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

Twelve rungs vary the body's shape and none varies its topology, which is what every impossibility here 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