Concept

Optimisation — where it appears

Choosing a map by minimising a stated quantity over a stated region, rather than by picking a construction and measuring what it does. It changes the question from which construction to use into what to minimise, and the answer then depends on the region, the criterion and nothing else.

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

Angular deformation against latitude, four projections. The same quantity for mercator, gallPeters, winkelTripel, robinson, from the equator to 80°. A conformal projection sits flat on zero in the first of these plots and runs away in the second; an equal-area projection does the reverse. Nothing is flat in both.

Every projection minimises something

A projection is the solution to an optimisation problem, and naming the objective explains more than naming the family. Some objectives are exact constraints, some are least-squares fits, and one is a table of numbers a man adjusted until it looked right.

choosing · Choosing
Fitting a sphere to 10° of latitude at 54.5°. The Earth's Gaussian radius of curvature across the band, with two candidate spheres drawn against it. The global mean radius is 6371.0 kilometres and sits 2203 parts per million away from the ground here — 2.2 metres in every kilometre measured. The best local radius is 6385.0 kilometres and has no bias at all by construction, leaving 325 parts per million of spread that no sphere can remove, because the curvature varies across the band and a sphere's does not. The gain is a factor of 6.9.

A datum is fitted to a region

Forty countries adopted forty ellipsoids, and the usual explanation is that measurement was poor. It was not — a figure fitted to one country beats the global one over that country by a factor of seven, and what it removes is a systematic bias of two metres in every kilometre.

datums · Datum
The scale factor was chosen, and the choice is a V. The worst departure of the grid scale factor from unity across a region 9° wide and 9° tall, for every scale factor between 0.998108 and 1.0002, on a central meridian of 2°W. A tangent projection touches at one meridian and is too big everywhere else, at 1987 ppm. Scaling the whole grid down slides that interval until it straddles unity, and the minimum sits at 1/√k_max = 0.999008102, where the worst departure is 993 ppm — a factor of 2.00, which is the most this construction can buy and is reached exactly. The published value marked beside it was chosen the same way, for this region, in the 1930s.

The scale factor was chosen

A grid's scale factor is not a constant of nature. Britain's 0.9996012717 is the reciprocal square root of the worst distortion a tangent projection would have had over the country, it halves the worst-case error exactly, and the halving is the most that construction can ever buy.

practice · Grid
Mercator swung through every tilt, over Chile and the tropics. The regional distortion of one projection as its axis is tilted from the normal aspect at 0° to the transverse at 90°, each curve divided by its own value in the normal aspect so the two regions can share an axis. Chile is best at a tilt of 90°, a factor of 183.5 better than north-up; the tropics is best at a tilt of 0°, which is north-up. Rotating the sphere costs nothing and changes none of the projection's own properties, which makes this the cheapest improvement available and the one most often left unmade.

Fitting the aspect to the region

Choosing a projection is a choice among a few dozen named things. Choosing its aspect is a choice among a continuum, it costs nothing, it changes none of the projection's own properties, and for a long thin country it is worth a factor of 183.

choosing · Choosing
What a grid of one's own is worth. For every candidate central meridian, the best worst-case distortion achievable over the region once the scale factor has also been optimised — so each point on the curve is already the bottom of its own V. The minimum is at 2.9°W with k₀ = 0.9993009, giving 700 ppm. The dashed line is UTM zone 30, whose meridian and scale factor were chosen for no region in particular, at 1059 ppm. The national grid is better by a factor of 1.51 — which is the whole of the answer to why a country publishes a grid rather than using the zones, and it is a smaller factor than the argument is usually made to sound.

Designing a grid for one region

A country with a grid of its own beats the international zone that would otherwise carry it by a factor of 1.51. That is worth having and it is much less than the argument is usually made to sound, and the difference between the two numbers is the whole case for and against national grids.

practice · Grid
The scale factor across Europe, from the middle outwards. Each band is the range the larger principal scale factor takes on the ring at that distance from the centre of Europe: 0 is the middle, 1 the frontier. two of the three projections drawn are conformal, and every one of those reaches its maximum at the right-hand edge, because the logarithm of a conformal map's scale factor is subharmonic. None of them reaches a maximum inside the region.

Where the worst point is

The largest scale error on a conformal map of a country is always on the frontier, never inside it, whatever the country's shape and whichever conformal projection was chosen. It is a theorem rather than a tendency, and it is the reason Chebyshev's criterion works.

distortion · Tissot
London to Tokyo, round a 15° exclusion. The direct great circle, dashed, runs through the disc. The admissible shortest route leaves it along a great circle tangent to the rim, follows the rim, and leaves along another tangent — which is the closed-form answer and is checked against a shortest-path search over the rim that knows nothing about tangents. It costs 39 kilometres on 9559, which is 0.41 per cent. The rim stretch is 293 kilometres of it, and it is the only part of the route that is not a geodesic anywhere along its length. Drawn in Orthographic.

A route that must go round

Every route on this site so far has been free to go anywhere, and no real route is. The shortest path past a circular exclusion is two tangent great circles and an arc of the rim — a closed form that agrees with a shortest-path search to three metres in 9,598 kilometres — and it costs not the obstacle's size but the square of how far the obstacle reaches past the route.

paths · Paths
Every mixture of Equirectangular at 50.46° and Aitoff. Airy's criterion over the whole sphere, for every weighted average of the Winkel tripel's own pair. Each mixture is normalised to its own best constant scale first, so the comparison is about shape rather than size. The curve has its minimum at w = 0.45, where the criterion is 0.3037 against 0.4114 and 0.4290 at the two ends. That is a projection nobody constructed beating both projections somebody did, by 26 per cent.

The average of two projections

The Winkel tripel is literally the arithmetic mean of two other projections, and this site's implementation of it agrees with that mean to the last bit. Averaging beats both ingredients by 26 per cent — and it preserves conformality exactly, destroys equal-area completely, and can turn eighteen per cent of the world inside out without either distortion measure saying so.

choosing · Choosing
The azimuthal family, as one function of one variable. Every azimuthal projection is a rule for how far from the centre to draw a point at angular distance ρ, and nothing else. five named projections are plotted as their radial functions — tan ρ, 2 tan(ρ/2), ρ, 2 sin(ρ/2), sin ρ — and the two dashed curves are the solutions of the two conditions, integrated from the centre with no more than "true scale at the centre" to start from. They land on the stereographic and the Lambert azimuthal to 8e-10. The names are the answers to the equations, not descriptions of viewpoints.

The azimuthal family is one function

Five azimuthal projections are taught as five viewpoints — from the centre, from the far pole, from infinity. They are one projection with five choices of how far out to draw a point at angular distance ρ, and each named property is a differential equation in that one function: f′ = f/sin ρ integrates to the stereographic projection and f f′ = sin ρ to Lambert's, from nothing but true scale at the centre.

families · Families
The taught rule against the measurement, on 30 regions. Each cell is a region built to order — a box of the stated height and width-to-height ratio, centred at the stated latitude — with the family that actually scores best over it, and whether that is what the rule says. Every family is given its own parameters for the region: the conic its cone constant, the cylindrical its standard parallel and the choice of a normal or transverse axis, the azimuthal its centre. The rule is right on 19 of 30, and where it is wrong it is wrong in one direction: it keys on latitude, and what decides the answer is shape.

The rule of thumb, scored

Cylindrical near the equator, conic in the middle latitudes, azimuthal at the poles. It is the most repeated piece of practical advice in cartography and it has never been run against a population of regions. Run against thirty, it is right nineteen times, and every one of its failures has the same shape.

wrong · Audit
Two rules, three populations. The taught rule keys on latitude; the replacement keys on shape. On the thirty regions the replacement was read from it scores 83 per cent against the taught rule's 63. On forty-five different regions built the same way it scores 91 — higher, not lower, so the generalisation the shortfall doubted is real. On the seven named regions this collection actually uses, both rules score 29 per cent, and following either costs a mean factor of 14.6. The third bar of each group puts the taught rule's polar clause back into the shape rule, which is what the out-of-sample failures ask for: it repairs four of them, breaks two that were right, and reaches 43 of 45.

The rule scored out of sample

A replacement rule was read off thirty regions and scored on the same thirty, and this collection recorded that as not being evidence about any other thirty. It is: on forty-five different regions the rule scores 91 per cent against the 83 it managed at home. What it cannot do is the seven regions the collection actually uses, where both it and the rule it replaced name the winner twice out of seven and cost a mean factor of 14.6.

wrong · Audit
Every aspect of Mercator over Japan, and the line the old sweep searched. The regional distortion of Mercator over Japan for every position of the projection's pole — darker is better — with the meridian the site's one-dimensional sweep searches drawn on it. The sweep's best is a gain of 24.4× over the normal aspect; the two-parameter search finds 61.4×, which is 2.51 times better again, at a pole 45° of longitude away from anything the sweep could reach. The shortfall recorded when the sweep was written said this would happen for a region whose long axis runs diagonally, and this is the measurement of it.

The aspect has three numbers, not one

The site's aspect search has swung the projection's axis through one plane for a long time, and recorded that a region whose long axis runs diagonally has its optimum somewhere that plane never reaches. Searching the whole sphere of pole positions finds 2.6 times more improvement over Japan and 3.4 over the conterminous United States.

choosing · Choosing
How thinly a sphere can be covered by a few equal caps. The covering density of the best arrangement of n equal caps found for each n — the total area of the caps divided by the sphere's, so a value of one would be a perfect tiling with no overlap. The horizontal line is 2π/√27 = 1.2092, the thinnest covering density of the PLANE by equal discs, which this site has used for the sphere since its first atlas essay. It is wrong in both directions: at 2 caps the sphere is covered more thinly than any plane can be, because a cap may be a hemisphere, and at every count from 3 upwards more thickly — 1.5092 at 3, and 1.3377 at 14. The ringed points are the four counts whose optimum is proved: 2 at 90.00°, 4 at 70.53°, 6 at 54.74°, 12 at 37.38°. Everything else is an upper bound from a search, drawn as one, and the bound loosens as the count rises — the search reaches the proved optimum to 3.4 per cent at these counts and has no such check anywhere else.

The sphere is not the plane at small counts

The site's atlas arithmetic multiplies an ideal sheet count by 2π/√27, the thinnest covering density of the plane. At the four counts whose optimal covering of the sphere is a theorem the plane's number is 21 per cent high at two caps and 9, 5 and 2 per cent low at four, six and twelve — wrong in both directions, and the direction changes with the count.

choosing · Choosing
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
The scale spread of every grid Britain could have adopted. A grid is conformal by requirement, so its angular deformation is zero everywhere and the whole design problem is the spread of its one remaining number: the largest scale factor over the region divided by the smallest, in parts per million. The grid's own scale factor does not enter — multiplying every scale by a constant leaves the ratio alone, which is why it is chosen last. The bottom bar is Chebyshev's optimum, the conformal map of this region whose scale is constant on its boundary, which no map of any family can beat; the adopted grid sits 1.98 times above it. Measured over a stated box rather than a coastline, because a coastline would put the vendor's generalisation into the answer.

The best grid a country could have had

A national grid is a conformal map chosen for one region, so its whole design problem is one number: the spread of its scale factor. That number has a theoretical floor, this site can now compute it, and the adopted grid turns out to be either exactly optimal or half as good again — depending entirely on which box the country is declared to be.

practice · Grid
The shortest route, and the shortest route a vehicle can fly. A leg of 60 kilometres for a vehicle whose minimum turning radius is 5, arriving on a heading 120° off the line and required to leave on one -60° off it. The straight line is the geodesic; the curve is the shortest curvature-bounded path, which is Dubins's RSR — a turn, a straight, a turn — at 67.29 kilometres against 60. The second curve is the runner-up word, drawn to show that the choice between them is a real one rather than a formality.

The shortest route a vehicle can fly

Nine rungs find the shortest path under a metric and none asks whether the thing travelling can follow it. Bound the curvature and the route depends on two headings as well as two positions: the turning cost is a fixed 1.81 kilometres whatever the leg length, so it is 18 per cent of a short leg and 0.28 per cent of a long one, and the whole of it vanishes when the vehicle happens to be pointing the right way.

paths · Paths

The landscape the search walks on

The three-parameter aspect search was run and its answer recorded with a note admitting nothing proved it global. Mapping the objective finds 26 to 34 local minima for every projection and region tried, a downhill walk from a random start reaching the best of them 6 to 35 per cent of the time — and one seed from the coarse grid the search already uses reaching it in all four cases. The score is reproducible to two per cent across a sevenfold refinement; the pole it names moves 60 degrees.

choosing · Choosing

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

Sixty zones was a decision about one latitude

Twelve rungs price a grid, a zone, an origin and a reference, and every one of them works inside a single zone. The number of zones has never been asked about: six degrees meets its tolerance at the equator and is loose everywhere else, so a system spending the same tolerance evenly would use 51 zones at the equator and 8 at 82° — and UTM's worst error is 981 parts per million, not the 400 always quoted.

practice · Grid

What the net heuristic cannot find

Rung seven chose a net by unfolding outwards from a face and showed the choice beats guessing, and recorded that its own limit was unknown. Enumerated in full, the heuristic turns out to be exactly optimal on every solid small enough to check — rank one of 384 — and above that ceiling no sample can tell whether it still is, because eight hundred random nets never reach it.

families · Polyhedral

The shape of the valley

An aspect search returns three numbers, two searches return triples that differ by a hemisphere, and the maps they produce agree. One cause is an exact degeneracy and the rest was called a valley and left unmeasured. Sampled densely, it is neither a valley nor a basin: a connected sheet spanning 170° of pole that fractures into fourteen pieces once the threshold tightens.

choosing · Choosing

Where the valley breaks in two

The previous rung found a near-optimal set that is one connected sheet at a loose threshold and fourteen basins at a tight one, and explained the transition without testing it. The explanation is a prediction about region size and projection sharpness: swept over both, the threshold falls from 2.13 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.

choosing · Choosing

The shortest route between two coasts

Ten rungs find the shortest path between two points. Between two regions the answer is attained at a pair of boundary points, and which pair is not the pair any page makes look nearest — a cylindrical map picks a pair 2,135 kilometres from the right one across the Atlantic, and 1,716 kilometres wrong between Chile and New Zealand, where its route is 22.5 per cent long.

paths · Paths

The nearest map to an impossible request

Rung seven found that not every distortion can be asked for. It never asked what happens when one is asked for anyway — and the answer has a shape: the achievable fields are the solutions of an elliptic equation, a request is a point off that set, and the nearest point leaves a residual whose floor is the curvature rather than the size of the ask.

choosing · Condition

How wrong a flat picture has to be

The rung below proves no flat picture of four places is exact and leaves the size of the failure to a determinant nobody can read. Measured directly, the least error falls as the square of how much sphere the places span — fitted exponent 2.0088 — and the same exponent comes back from five different arrangements while the constant in front of it moves by a factor of seventeen.

impossibility · Embedding

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.

impossibility · Embedding

A family is not closed under averaging

Nine rungs treat a family as a set of maps with a parameter running through it, and this collection's own compromise projections are averages of members. An average of two members is not a member: two conics far apart average to something 31 per cent of their own separation outside the family — and averaging within a family makes the map worse every time, while averaging across two makes it 16 per cent better.

families · Families

The basins have widths as well as depths

The previous rung measured the height of the pass and recorded a shortfall: shape means widths too. Measured, the basin has three of them — 32°, 16° and 7° at Japan — it gets wider rather than narrower as the region grows, and the exponent it predicts overshoots the measured one by half again.

choosing · Choosing

Which projection a weighting can make best

Rung eight finds the seven world projections nothing beats on both counts and tells 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.

wrong · Audit

The ranking is not an order

The previous rung showed that a weighting can make almost any projection best. 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.

wrong · Audit

The first break is mostly its denominator

Three rungs 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 pass's own depth is constant to 12 per cent below twenty degrees of span, and the whole of the threshold's movement there is the denominator — the best score the region admits at all — rising with exponent 0.92.

choosing · Choosing

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 family is a function, not a list

The equal-area pseudocylindricals are the solutions of one equation in two unknown functions, so one function is free and the named members are points in a space of them. Writing that function as six numbers and searching over them finds a map at 25.03° of mean angular deformation against Eckert IV's 28.08 — and two numbers are already enough to beat every map the family has a name for.

families · Families

The pooled score abandons a region

Fifteen rungs optimise 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.

choosing · Choosing

The maps with no family are simply better

Rung eight found seven projections with no continuous symmetry and noticed they are almost exactly the set anybody would choose for a world map, then offered a conjecture with 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.

families · Families

The cheapest map that meets its areas

An earlier essay bracketed a cartogram's least cost between a construction charging eighty degrees and a bound valid only for symmetric densities, and recorded the gap as a shortfall. One request settles it: a density of contrast eighty whose least cost is exactly zero, met by a map written down in closed form, while the standard construction charges 43.8° for it.

distortion · Cartogram

The best compromise for angle is not the best for bending

Every compromise projection in the library is an average of two others, and averaging is a first-order operation — so the second derivative was never part of the bargain. Swept along five ordinary blend paths, the weight that minimises angular deformation and the weight that minimises flexion are between a quarter and a half of the axis apart, and how much a compromise buys at one order predicts nothing about the other.

distortion · Flexion

Named alongside it

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

PurposeAspectVerificationClosed formObjective functionTrade-offConformalityKavrayskiy's criterionRegional distortionToleranceDegeneracyRegion

All concepts