Concept

Convergence — where it appears

The approach of a sequence of approximations to a limit, and the rate at which the approach happens. It is the quantity this collection refines a grid or shortens a step to measure, because a rate is evidence about a method where a single error is not.

Named by 27 essays across 9 fields — each of them below, with the objects they name alongside it.

Transverse Mercator. The graticule of the Transverse Mercator projection at 30° of longitude and 15° of latitude. the only projection on which a constant compass bearing is a straight line. It is conformal.

Transverse Mercator and the series that computes it

The projection most of the world's survey data lives in has no closed form. What every national grid actually computes is a truncated power series in the flattening, and how far it can be trusted is an engineering parameter rather than a property of the projection.

families · Ellipsoid
Where the inverse problem stops converging, 3° around the antipode of London. Every target in a 6° square centred on the point diametrically opposite London, shaded by how many iterations Vincenty's inverse formula needed to find the geodesic to it. 3 of 625 — 0% — never converged at all, and the worst success took 113 steps against a handful anywhere else on Earth. The failure is not a defect in the formula: near the antipode the shortest path is nearly ambiguous, and an iteration looking for one answer is being asked which of many.

The route with no shortest path

Between a point and the point diametrically opposite there are infinitely many shortest routes and no shortest route, and the standard formula for the distance between two places stops converging in a neighbourhood of it. The failure is a property of the question rather than a defect in the answer.

paths · Paths
A direction carried once round the sphere. A vector transported round a closed loop on the sphere, kept as parallel to itself as the surface allows at every step — the component that leaves the tangent plane is removed and nothing else is done to it. The heavy arrows are its direction at the start and at the finish, drawn from the same point; the light ones are its direction along the way. It comes back turned through 282.0 degrees, which is 0.783 of a revolution, and nothing in the transport turned it. Drawn in an orthographic projection of the embedding, which is itself a map and has its own distortion.

A direction carried round a loop

Carry a bearing round a circuit, keeping it as parallel to itself as the surface allows, and it comes back turned. Round a parallel at 45° the turn is 4.443 radians and the cap enclosed is 1.840, and they sum to exactly one revolution — so the turning is not the curvature, and on a cone it is all of one and none of the other.

impossibility · Curvature
Grid north against true north at 52°, across a 6° zone. seven meridians drawn in the grid's own coordinates. The grid's north is straight up the page everywhere by construction; the meridians are not, and the angle between them is the convergence. It is exactly zero on the central meridian, reaches 2.36° — 142 minutes of arc — at 3° out, and changes sign across the middle. The manual's γ = Δλ sin φ is 3.0″ short of it at the zone edge.

Grid north is not north

A grid has one north, parallel everywhere on the sheet by construction. The Earth has a different one at every point. The angle between them reaches two degrees at the edge of a UTM zone, and a straight line on the grid is not a straight line on the ground either.

families · Ellipsoid
Where two equal geodesics meet, from 15° north. A geodesic leaving at azimuth α and its mirror image at −α have the same length wherever they meet, and by symmetry they meet on the antipodal meridian. On a sphere they all meet at one point, the antipode, to 1.4e-7° — the flat line. On the ellipsoid the meeting latitude moves with the azimuth, from -15.563° to -15.009°, so the set of points with two shortest routes is an arc 61 km long, and the distance to it varies by 31 km along its own length.

Where the shortest route stops being the only one

On a sphere there is exactly one point with no shortest route from a given place: the antipode. On the ellipsoid the Earth actually is, that point is an arc — sixty-six kilometres of the antipodal meridian for a point on the equator, half a kilometre for one at 85°, and every point of it reachable by two different geodesics of exactly equal length.

paths · Paths
One length, three corrugations. A straight segment shortened by a factor of 0.9, then wiggled across its own direction until its length is back to what it was. The wiggle's amplitude is what buys the length and the number of wiggles is free, so all three curves have exactly the target length while the third stays 16 times closer to the shortened segment than the first. That is the mechanism: the family converges to a map that is not isometric, while every member of it is.

Impossible in two derivatives, possible in one

The impossibility this whole collection rests on computes a second derivative, so it is a statement about maps that have two. Take one away and it is false: a corrugation restores an exact length while converging to the map that does not, and iterating it gives a flattening whose derivative converges and whose curvature runs to half a million.

impossibility · Curvature
The order the corrections go in, and what it costs. Setting out a 51.8 km line on the UTM zone 31N means turning a grid bearing into an azimuth to observe, and there are two corrections: the meridian convergence, 0.8669° here, and the arc-to-chord correction, 6.242″. Each bar carries the lateral offset it produces at the far end of the line, because a bearing error is a number nobody can picture and a sideways miss is the thing that misses. The exact arc-to-chord and the classical formula every manual gives differ by 0.0044″, which is 1.10 mm at the far end — small, real, and the reason the corrections have an order rather than a sum.

Setting out runs the chain backwards

Turning a design coordinate into something to observe on the ground means undoing the reduction chain, and undoing a chain reverses the order as well as the operations. Two of the corrections do not commute, and getting them the wrong way round misses by a millimetre at fifty kilometres.

practice · Reduction
How large a map has to be before a wrong projection stops fitting it. Three rivals fitted to a Mercator graticule centred at 45° north, over regions from 1° to 70° of half-extent. The residual is the shape difference alone, with the best scale, rotation and offset removed, and both axes are logarithmic. The horizontal rule is a fifth of a per cent of the map's width — about the width of a drawn line on a printed sheet — and where a curve is below it, no measurement of that map can tell the two projections apart, however carefully it is made.

Two projections that cannot be told apart

Over a small enough region every projection is the same picture, so the question is how small. The answer is not one number: two conformal projections need twelve degrees of extent before their graticules can be separated, and two projections with different anisotropy separate below one.

wrong · Identify
The boundary of the region the series is the map in. The curve where the transverse coordinate reaches 2.918, which is where the terms stop shrinking. It crosses the equator 83.81° from the central meridian and closes towards the poles, because the same longitude is a smaller transverse coordinate at a higher latitude — the boundary is a curve rather than a meridian. The narrow band beside it is a 3° zone, the width national grids actually use, drawn to the same scale: the practical world sits in about a fiftieth of what the series can reach. Drawn in Mollweide.

Where the series stops being the map

The transverse Mercator has no closed form on an ellipsoid, so every national grid computes a truncated series. Asking how many terms it needs has an answer everywhere; asking whether more terms always help has an answer only within 83.81° of the central meridian, and the boundary is computed rather than assumed.

families · Ellipsoid
Where the projection's pole should go, and the answer the third rotation moves it to. Every dot is a pole position the search tried, sized by the best score it can reach there when the third rotation is also free — the score is the distortion of Robinson over Europe by Kavrayskiy's criterion, so smaller is better and the large green dots are the good regions. The circled mark is the two-parameter optimum and the square is the three-parameter one: they are 74 pixels apart on this map, which is a different aspect rather than a refinement of the same one. The search costs 8 times the evaluations of the two-parameter one. Drawn in Mollweide.

The third parameter, run

An aspect has three numbers and this site has been searching two of them, with a note admitting it. Searching all three is worth up to 2.1 times — and the obvious way to do it, starting from the two-parameter answer and letting the third move, finds a fraction of that or nothing at all.

choosing · Choosing
Going out and coming back, at three orders. A point three degrees from the central meridian, taken forward into the projection and back again with the series truncated at the same order both ways, and the ground distance between where it started and where it returned. The second-order pair is out by a tenth of a metre at some latitudes; the third by half a millimetre; the fourth — the order every national grid formula in ordinary use is written to — by 0.39 mm. The dips are where one of the two series passes through a node, not where the map is better.

The inverse of the series is not the series of the inverse

Transverse Mercator on an ellipsoid has no closed form, so every national grid computes it as a truncated series. There are two of them — one out and one back — and they are separate approximations. Composing them does not give the identity: at the order every grid formula in ordinary use is written to, a point three degrees from the central meridian comes back 0.39 mm north of where it started.

families · Ellipsoid
The score does not settle at any resolution. The compactness of one stated boundary — a circle with cosine ripples at eight geometrically spaced wavenumbers, so it has structure at every scale — read at sixteen vertices up to two thousand and forty-eight. The ground score falls from 0.980 to 0.834, and it keeps falling: the boundary's length grows without bound as it is resolved while the area it encloses converges, so the quotient has no limit. The four page curves sit within a fraction of a per cent of the ground curve and of each other, which is the comparison this rung exists to make.

The score is not stable at any scale

One boundary, read at eight resolutions from sixteen points to two thousand and forty-eight: the compactness score falls from 0.980 to 0.834 and is still falling. Changing the projection instead moves it by 0.69 per cent. The two decisions are made by the same person on the same afternoon and only one of them is ever reported.

paths · Reach
The indicatrix is a limit, and here it is being taken. The axis ratio of the image of a circle of radius ρ, measured from the picture rather than computed from the derivatives, against the ratio the derivatives give. It converges, and it converges as ρ² — each halving of the radius quarters the departure. At the smallest radius drawn the two agree to 7.6e-4 relative. This is what "the indicatrix is a limit" means when it is working: there is a number for the ratio to settle on, and it settles.

The indicatrix at a point that has none

Ten essays on this ladder take the derivative for granted. Rung three proved the indicatrix is a limit and every figure since has assumed the limit exists. There are points on the maps this collection draws where it does not, and they are not exotic: the horizon of a gnomonic map, and every corner of every polyhedral net.

distortion · Tissot
The parameters are not reproducible and the map is. The three-parameter aspect search run at three grid resolutions and compared with the finest, twice over. Compared on the numbers it returns, the answers are 65° of pole apart. Compared on what they do to the region — the root-mean-square difference in angular deformation at every sample — they are 0.29° apart, against a map whose own deformation over that region averages about a degree. The disagreement recorded as a shortfall is a disagreement about coordinates for one map.

Report the map, not the parameters

The previous rung found the aspect search returning the same score to 2.3 per cent from poles sixty degrees of latitude apart, and recorded that as a shortfall: the answer was not reproducible. The shortfall assumed the disagreeing triples make disagreeing maps. They do not — the three answers agree on the distortion field to a quarter of the deformation the map already has.

choosing · Choosing
Equal area or steady shape, and not both. Four cell schemes plotted by how much their cells vary in area and how far from square the worst of them is. The bottom-left corner is the scheme that has both, and it is empty: the equal-area cube holds area to 1.003 and has the most elongated cells, the tangent-warped cube has the tightest shapes and lets area vary by 1.20, and the lon/lat scheme is off the scale on both. Neither axis can be driven to one while the other stays there.

A cell system trades area for shape

A grid can hold every cell to exactly the same area or hold every cell nearly square, and the measurement says it cannot do both: the equal-area cube's areas agree to a part in a thousand and its worst cell is 1.29 times as long as it is wide, while the tangent-warped cube holds shape to 1.19 and lets area vary by 20 per cent.

applied · Cells
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.

A crossing is a chain of decisions

Eleven rungs hand back a curve and stop, 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.

paths · Paths

Further north on the grid is not further north

Thirteen rungs price a grid's origin, scale, units and zones, and every one treats a coordinate as a position. It is also an ordering, and the two disagree: a parallel climbs 4,379 metres of northing on its way to the edge of a UTM zone, so two places 4.4 kilometres apart in latitude can be listed in the wrong order — and the share of pairs it happens to is the convergence, in a different unit.

practice · Grid

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.

practice · Grid

A condition imposed at points is not a condition

Nine rungs state a condition and solve it, and every solve imposes the 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.

choosing · Condition

A line has a length only at a scale

Every measurement on this site so far has been of a curve given by a formula, sampled as finely as the picture needed. 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.

applied · Generalise

The nodes were evenly spaced

The previous rung showed that 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.

choosing · Condition

A long window and a square one

Rung nine held the patch's area and found that where the cut goes still changes the map by a third, with the curvature it encloses predicting the change at r = 0.969. It recorded that it had held area and curvature and not shape. Sweeping the shape at constant area separates two things that had looked like one: the map is curvature and the cost is shape.

datums · Bodies

Two routes to one scale

A national series is cascaded — the million is derived from the quarter-million, which was derived from the fifty — and the folklore is that the errors accumulate. They do not: Douglas–Peucker and Visvalingam both cascade to the same line the direct route produces, bit for bit, because both output a sublevel set of a per-vertex number. What breaks it is anything else in the chain, and a moving average puts the product nine tolerances away.

applied · Generalise

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.

distortion · Sampling

A refinement that stops moving

Doubling the sample and watching the answer settle is how every quadrature in every field is checked. On the Robinson projection the doubling ladder — 4, 8, 16, 32, 64, 128 — converges beautifully, with its increments halving at every step, on a limit that is wrong by a factor of twenty-seven. Whether a grid finds the answer is decided by whether n is a multiple of four.

distortion · Sampling

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.

distortion · Sampling

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.

distortion · Tissot

Named alongside it

The objects these essays reach for when they reach for this one.

VerificationToleranceClosed formConformalityTransverse MercatorUTMAnisotropyEstimatorNational GridSamplingZoneDegeneracy

All concepts