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

Distortion over a region, and the map fitted to one

region-ranking is the ranking a region induces; chebyshev-cap is the one case with a proved optimum; sheet-ladder is what a sheet series costs; boundary-extremes and corridor-scale are the same quantity read along a boundary and along a corridor; and solved-map and family-solve are the map that comes back when the region is handed to a solver rather than to a taxonomy.

region-ranking is the ranking a region induces; chebyshev-cap is the one case with a proved optimum; sheet-ladder is what a sheet series costs; boundary-extremes and corridor-scale are the same quantity read along a boundary and along a corridor; and solved-map and family-solve are the map that comes back when the region is handed to a solver rather than to a taxonomy.

Every one of these is the same generator answering a different question, which is why they share a file, a set of colour roles and a set of assertions. Changing it changes all 7.

36 essays call it. Every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

What it draws

One picture per question the family answers. Each has a page of its own.

The same six projections, ranked over Europe and Chile. Each column orders the projections by Kavrayskiy's regional criterion — the root-mean-square departure of the two principal scales from unity, integrated over the region with the area element. The lines cross, which is the point: Lambert conformal conic leads over Europe and comes fifth over Chile. A table of projections ordered by distortion is a table about somebody's region.

region-ranking

The same six projections, ranked over Europe and Chile The least distortion possible over a 30° region. Scale factor along a radius of the cap, each projection normalised to unit scale at the centre. Chebyshev's criterion names the projection whose scale is constant on the boundary as the conformal map of least scale variation, and for a cap that is the stereographic projection centred on it — reaching exactly sec²(ρ/2) = 1.0718 at the rim, marked. Every other conformal projection drawn here rises past that line before it gets there. This is the only figure on the site showing an optimum rather than a comparison.

chebyshev-cap

The least distortion possible over a 30° region Sheets for a tolerance, from Chebyshev's bound and a covering. For each stated tolerance on the scale error, the cap radius at which the best possible conformal projection just meets it — sec²(ρ/2) − 1 = tolerance, which is Chebyshev's bound and has no fitting in it — and then the number of such caps needed to cover the sphere at the packing density a real arrangement achieves. One part in a thousand costs 1210 sheets of 403 kilometres radius. The slope is -0.989: a factor of ten in what the job will accept is a factor of ten in the atlas.

sheet-ladder

Sheets for a tolerance, from Chebyshev's bound and a covering The minimum-distortion conformal map of an elongated region, 30° by 10°. The conformal projection of an elongated region, 30° by 10° whose scale is constant on the boundary, which is Chebyshev's criterion, obtained by fitting eight terms of a series rather than by choosing a named projection. The scale factor runs from 0.98640 to 1.00000, a spread of 1.01379; on the boundary itself the largest departure from constancy is 2.13e-7 in the log, which is what the fit achieved and not what it was told. Each dot is an interior sample shaded by its own departure from the boundary's scale.

solved-map

The minimum-distortion conformal map of an elongated region, 30° by 10° Solving the conformal condition, and finding Mercator. The cylindrical family is x = λ and y = Y(φ), one unknown function. Requiring the map to be conformal makes the condition a first-order differential equation in Y — Y′ = sec φ — and this curve is the Runge–Kutta integration of it, drawn against Mercator's own closed form as a dashed line. They agree to 1.4e-12 of the ordinate, so the named projection is not quoted anywhere in the calculation: it is what the condition returns.

family-solve

Solving the conformal condition, and finding Mercator 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.

boundary-extremes

The scale factor across Europe, from the middle outwards The scale along the Madrid–Tokyo corridor. The scale factor of four projections along the great circle from Madrid to Tokyo, each normalised to its own average over the route so the comparison is of variation rather than of size. The oblique Mercator whose own equator is laid along the corridor holds the scale to 0 parts per million; Mercator varies by 118.7% over the same line. A corridor is a curve, not a region, and the projection an area criterion picks is not the one a curve wants.

corridor-scale

The scale along the Madrid–Tokyo corridor

Where it is called

Changing this changes every one of these figures.

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. Grids, and what a survey does

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.

 The least distortion possible over a 30° region. Scale factor along a radius of the cap, each projection normalised to unit scale at the centre. Chebyshev's criterion names the projection whose scale is constant on the boundary as the conformal map of least scale variation, and for a cap that is the stereographic projection centred on it — reaching exactly sec²(ρ/2) = 1.0718 at the rim, marked. Every other conformal projection drawn here rises past that line before it gets there. This is the only figure on the site showing an optimum rather than a comparison. The impossibility

Total curvature and the scale rule

The impossibility has a size. A region covering a fraction of the sphere carries a fixed amount of curvature that any flat map must absorb, and for a circular region the least distortion any conformal projection can achieve is a closed form nobody can beat.

 The least distortion possible over a 30° region. Scale factor along a radius of the cap, each projection normalised to unit scale at the centre. Chebyshev's criterion names the projection whose scale is constant on the boundary as the conformal map of least scale variation, and for a cap that is the stereographic projection centred on it — reaching exactly sec²(ρ/2) = 1.0718 at the rim, marked. Every other conformal projection drawn here rises past that line before it gets there. This is the only figure on the site showing an optimum rather than a comparison. What each projection optimises

Chebyshev's criterion

The only optimality theorem in the subject. Among all conformal projections of a region, the one with least scale variation is the one whose scale factor is constant on the boundary — a criterion with a proof, a unique answer, and a test the site can run.

 The angles of a triangle, and the curvature inside it. A triangle on the sphere with sides that are great-circle arcs. Its three angles sum to 212.26°, overshooting the flat 180° by 32.26°. Integrating the curvature over the interior gives 0.56306 against an excess of 0.56306 — the same number by two routes that share no arithmetic. That is Gauss–Bonnet, and it is how a flatlander measures the curvature of a world it cannot step outside. The impossibility

Measuring curvature from inside

A flatlander with a ruler can find out that its world is round, by drawing a triangle and adding the angles. Gauss–Bonnet turns that into an exact statement, and the total curvature of a closed surface turns out to be a number obtained by counting.

 The least distortion possible over a 20° region. Scale factor along a radius of the cap, each projection normalised to unit scale at the centre. Chebyshev's criterion names the projection whose scale is constant on the boundary as the conformal map of least scale variation, and for a cap that is the stereographic projection centred on it — reaching exactly sec²(ρ/2) = 1.0311 at the rim, marked. Every other conformal projection drawn here rises past that line before it gets there. This is the only figure on the site showing an optimum rather than a comparison. Measuring distortion

Scale distortion is the third failure

A map's stated scale is its scale at one place. Everywhere else it differs, by a factor that no conformal projection can avoid and that neither of the two usual distortion measures reports. It is the failure everybody uses and nobody counts.

 The same six projections, ranked over Europe and Chile. Each column orders the projections by Kavrayskiy's regional criterion — the root-mean-square departure of the two principal scales from unity, integrated over the region with the area element. The lines cross, which is the point: Lambert conformal conic leads over Europe and comes fifth over Chile. A table of projections ordered by distortion is a table about somebody's region. Measuring distortion

Distortion over a region

An indicatrix describes a point and every practical question is about a country. Going from one to the other means integrating, integrating means choosing a weighting, and the weighting is the step that turns a measurement into somebody's opinion.

 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. Grids, and what a survey does

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.

 How large a patch can be treated as flat, at 10 parts per million. The smallest scale distortion any map of a circular patch can have, against the radius of the patch, on logarithmic axes. The line is straight with a slope of 2.00: the error grows as the SQUARE of the size, so a patch ten times wider is a hundred times worse. A tolerance of 10 ppm is reached at a radius of 40.3 km — 81 km across — and that is the number behind the boundary between plane surveying and geodesy. The impossibility

How small is flat enough

A builder works in plane coordinates and a national mapping agency does not, and the line between them is not a convention. The unavoidable error of treating a patch of the Earth as flat grows as the square of its size, and the size at any stated tolerance is a number.

The whole library