What each projection optimises

How many sheets an atlas needs

A tolerance on the scale error inverts, through Chebyshev's bound, into a sheet radius — and a covering problem turns the radius into a count. One part in a thousand costs 1,210 sheets of 403 kilometres radius, the count goes as the reciprocal of the tolerance exactly, and the projection multiplies it by anything from one to fifty-six.

Every argument about which projection a map series should use assumes the series exists. Before that decision there is an earlier one that nobody frames as a projection question at all: how many sheets.

It is a projection question, it has an answer in closed form, and the answer says that the sheet count is settled by the tolerance long before anybody chooses a projection — and then multiplied by the projection anyway, by a factor of up to fifty-six.

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.
Fig. 1 For each tolerance on the scale error, the largest cap the best possible conformal projection can cover within it, and the number of such caps needed to cover the sphere. The slope is −0.989: a factor of ten in what the job will accept is a factor of ten in the atlas. The ringed point is the International Map of the World, whose 4° × 6° sheets correspond to a tolerance of one part in 1,560, and whose predicted count of 1,888 sits against a published extent of about 2,500.

The first half is a theorem

Chebyshev’s criterion says that among all conformal projections of a region, the one with least scale variation is the one whose scale factor is constant on the boundary. For a spherical cap that projection has a name — the stereographic projection centred on the cap — and its scale spread has a closed form:

kmaxkmin  =  sec2 ⁣(ρ2)\frac{k_{\max}}{k_{\min}} \;=\; \sec^2\!\left(\frac{\rho}{2}\right)

where ρ\rho is the cap’s angular radius. That is a bound on every conformal projection of the cap, not a property of one of them, which is what Chebyshev’s criterion is about and what makes this arithmetic a lower bound rather than an estimate.

Inverting it costs one line. A stated tolerance τ\tau on the scale error gives

ρ  =  2arccos11+τ\rho \;=\; 2\arccos\frac{1}{\sqrt{1+\tau}}

and for small τ\tau that is ρ2τ\rho \approx 2\sqrt{\tau}: the radius goes as the square root of the tolerance, so the cap’s area goes as the tolerance itself.

The second half is a covering problem

A sphere of area 4π4\pi covered by caps of area 2π(1cosρ)2\pi(1-\cos\rho) needs at least 4π/2π(1cosρ)4\pi / 2\pi(1-\cos\rho) of them, and that number is achievable by nothing: discs overlap, and a disc is not a hexagon.

The thinnest covering of the plane by equal discs is the hexagonal arrangement, at a density of 2π/27=1.20922\pi/\sqrt{27} = 1.2092 — every point covered 1.21 times over on average. The sphere’s covering numbers are not exactly the plane’s at any finite count, so both figures are carried: the ideal, which is a bound, and the hexagonal, which is what an arrangement achieves.

Covering, not packing. The two are different problems with confusingly similar constants — the hexagonal packing density is 0.9069, which is how much of the plane discs fill without overlapping — and using the packing number here would be asking for sheets that leave gaps between them.

tolerance sheet radius ideal count hexagonal count
1 in 10,000 127 km 10,001 12,093
1 in 3,300 221 km 3,334 4,032
1 in 1,000 403 km 1,001 1,210
1 in 330 697 km 334 404
1 in 100 1,270 km 101 122
1 in 10 3,903 km 11 13

The exponent is the whole of the first result. The count goes as τ1\tau^{-1}, fitted at −0.989 across three decades, because the cap’s area goes as τ\tau and the count goes as one over the area. Ten times the precision is ten times the atlas, and no cleverness anywhere in cartography changes that.

Against a real series

The International Map of the World, agreed in 1913 and never finished, laid the world out in sheets of 4° of latitude by 6° of longitude at 1:1,000,000. A 4° × 6° cell at mid-latitude fits inside a cap of about 2.9° radius, which corresponds to a tolerance of one part in 1,560.

The arithmetic above predicts 1,888 sheets at that tolerance. The published series ran to roughly 2,500, with the excess accounted for by the two things this calculation deliberately leaves out — sheets are rectangles rather than discs, and the scheme merges cells above 60° instead of shrinking them.

Landing within a third of a real answer from two closed forms and no data is the check worth having. It says the tolerance-to-count chain is measuring the thing the sheet layout was actually constrained by, and not some idealisation of it.

The projection multiplies it, and the multiplier is not small

The essay was written expecting to show that the projection choice is a rounding error beside the tolerance. That is false, and the machinery refused it.

What the choice of conformal projection costs, over a 10° cap. Each conformal projection's scale excess over a cap of 10° radius, divided by the least any conformal projection can achieve — which Chebyshev's theorem says is the stereographic projection centred on the cap. The penalty is the factor by which the sheet count is multiplied at a fixed tolerance. Mercator costs 56 times the sheets, and the same projection pointed 20° away from the region costs 5.8. The tolerance sets the order of magnitude; this sets the multiplier, and it is the one a cartographer decides.
Fig. 2 Each conformal projection’s scale excess over a 10° cap, divided by the least any conformal projection can achieve — which is the number of sheets it costs at a fixed tolerance. The stereographic projection centred on the cap is the optimum by Chebyshev’s theorem, so it is one by definition. The same projection pointed 20° away costs 5.8; the Lambert conformal conic with its standard parallels elsewhere costs 4.3; Mercator costs 55.9.

Fifty-six times the sheets at the same tolerance is not a detail. At one part in a thousand it is the difference between 159 sheets over the region and 8,838.

So the count has two factors and they multiply:

N    1.209×4π2π(1cosρ(τ))  ×  ϵchosenϵoptimalN \;\approx\; \frac{1.209 \times 4\pi}{2\pi\left(1-\cos\rho(\tau)\right)} \;\times\; \frac{\epsilon_{\text{chosen}}}{\epsilon_{\text{optimal}}}

The first is arithmetic on the tolerance and nobody can improve it. The second is entirely the cartographer’s, it is available for free, and it is the one every argument about projections is actually about.

Why the second factor is so large

A conformal projection’s scale spread over a region is set by how far the region is from where the projection is exact. The stereographic projection centred on the cap is exact at its centre and the spread grows as ρ2\rho^2; Mercator over the same cap at 45° north is exact at the equator, forty-five degrees away, and the spread over the cap is dominated by the variation of secφ\sec\varphi across it rather than by the cap’s own size.

That is the same mechanism where the worst point is measures: on a conformal map the extremes are always on the boundary, so the spread is a statement about the boundary’s distance from the point of exactness.

The least distortion possible over a 10° 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.0077 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.
Fig. 3 The optimum and its rivals over a 10° cap at 45° north. The theorem’s projection achieves the bound to four decimal places, which is the check that the bound is being computed rather than quoted, and the rivals are above it by the factors the previous figure ranks. The one worth watching is the same projection centred 20° away: nothing about it is wrong except where it is pointed.

The same arithmetic for a country

A national series covers a region rather than a sphere, so the covering runs over an area rather than over 4π4\pi, and the numbers become recognisable.

Britain is about 244,000 square kilometres. At one part in a million — a tolerance a survey control network works to — a sheet may be 12.7 kilometres in radius and the country needs 579 of them. At one part in a hundred thousand a sheet may be 40 kilometres and 58 will do. At one part in ten thousand, six.

tolerance sheet radius sheets for Britain
1 in 1,000,000 12.7 km 579
1 in 100,000 40.3 km 58
1 in 10,000 127 km 6
1 in 1,000 403 km 1

The last row is the one that explains national grids. A country the size of Britain fits inside a single sheet at one part in a thousand, so a national mapping agency that wanted only that would need no zone system, no seams, no sheet index and no edge-matching — which is exactly what a national grid is, and is why designing a grid for one region finds the gain over an international zone smaller than the argument sounds. The gain in distortion is 1.51; the gain in administration is the whole of the difference between one sheet and sixty.

Conversely, the 579-sheet row is why survey control is not published as a map at all. At the tolerance a control network needs, no map covers a country, and the answer is a coordinate list plus a projection formula rather than a series of sheets.

What the choice of conformal projection costs, over a 5° cap. Each conformal projection's scale excess over a cap of 5° radius, divided by the least any conformal projection can achieve — which Chebyshev's theorem says is the stereographic projection centred on the cap. The penalty is the factor by which the sheet count is multiplied at a fixed tolerance. Mercator costs 101 times the sheets, and the same projection pointed 20° away from the region costs 11.5. The tolerance sets the order of magnitude; this sets the multiplier, and it is the one a cartographer decides.
Fig. 4 The same penalties over a 5° cap, which is a quarter of the area and a quarter of the tolerance. The ordering is unchanged and the factors are close to the 10° case, because the penalty is a ratio of two scale excesses and both scale together — so the projection multiplier is very nearly independent of the sheet size, and multiplies the tolerance’s count cleanly.

The tolerance is not a technical parameter

The uncomfortable part of this arithmetic is what it does to the word “accuracy”. A series specified at one part in ten thousand is not a hundred times better than one at one part in a hundred; it is a hundred times larger, and the reader of a single sheet cannot tell which they are holding.

The trade is stated cleanly by the two exponents. Sheet radius goes as τ1/2\tau^{1/2} and sheet count as τ1\tau^{-1}, so halving the sheet radius quadruples the atlas. A national series at 1:25,000 with a 10-kilometre sheet is not a decision about paper size; it is a decision that the country will be covered by some thousands of sheets rather than some tens.

Which is why the sensible question, and the one this arithmetic makes askable, is the inverse: given how many sheets a programme can afford to print and maintain, what tolerance does that buy? That is a division rather than a search, and it is the number that should be in the specification.

The bound is for conformal projections only, and that is a restriction

Chebyshev’s theorem is a statement about the conformal projections of a region. It says nothing about the others, and the others are not obviously worse: an equal-area projection of a small cap has an areal factor of exactly one everywhere, which is a better result than any conformal projection manages about area.

What it does not have is a single scale factor. A conformal projection’s error can be stated as one number per point, so a tolerance on it is a tolerance on everything a sheet is used for; an equal-area projection has two different scales at each point and a tolerance has to name which. That is the distinction the two ways a map is wrong makes, and it is why sheet specifications are written in terms of scale error and therefore in terms of conformal projections.

So the count above is the count for a series whose sheets are meant to be measured on with a ruler in any direction. A series meant for area — a cadastral or a land-cover product — is a different optimisation with a different bound, and this site does not have the equal-area analogue of Chebyshev’s theorem because there does not appear to be one.

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.
Fig. 5 Where each projection’s extreme scale factors fall over a region. On the two conformal members every extreme is on the frontier, which is the maximum principle and is what makes a boundary condition — Chebyshev’s criterion — the right thing to impose. On the equal-area member the extremes are not constrained that way, because the quantity being extremised is no longer the modulus of a holomorphic derivative.

What this shares with the grid essays

Designing a grid for one region asks a version of this question with the answer fixed at one sheet: given a country, how much distortion does a single projection impose, and is a national grid worth having over an international zone. The answer there was a factor of 1.51, which is real and much smaller than the argument is usually made to sound.

This is the same trade with the sheet count as the free variable rather than the distortion. A zone system is a covering: UTM and the zone system is sixty caps of a particular shape, chosen so that the worst scale error inside one is 1,000 parts per million, and sixty is what that tolerance costs at that shape.

The sixty zones, each six degrees wide. Every zone is a separate transverse Mercator projection about its own central meridian, so the world is covered by sixty maps rather than one. Zone 31 is picked out, running from 0° to 6° with its axis on 3°. Coordinates do not carry across a zone boundary — a point on either side of one has two entirely different eastings, and nothing in the numbers says which zone they belong to.
Fig. 6 UTM’s sixty zones, which are a covering of the world by sheets whose tolerance was fixed first. Each is six degrees wide with a scale factor of 0.9996 chosen so that the error at the edge and the error at the centre are equal in magnitude, which is the same balancing Chebyshev’s criterion performs in closed form for a cap.

The two schemes differ in shape — zones are gores and sheets are discs — and the arithmetic differs only in the covering density that shape implies. Gores tile exactly, with density 1, and pay for it by being long: the worst point in a zone is at its corner, far from the central meridian, so the tolerance a gore achieves is worse than a disc of the same area would.

The flat-Earth version of the same question

There is a coarser calculation a builder makes that has the same shape and a different exponent, and putting them side by side says what the sheet count is actually a statement about.

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.
Fig. 7 The error of treating a patch of the Earth as flat, against the patch’s size. That error grows as the square of the radius, so a stated tolerance buys a radius as its square root — the same square root the sheet radius obeys, because both are second-order departures from a first-order description.

The flat-Earth radius and the sheet radius obey the same law because they are the same quantity: how far a locally correct description survives. The difference is what is being described. A plane is correct at a point and wrong quadratically; a conformal projection is correct at a point and wrong quadratically; and the sheet count is the flat-Earth question asked over a whole sphere with the answers counted rather than quoted.

That also gives the honest reading of “the best possible projection”: it does not remove the second-order term, it moves the point where the term vanishes to the middle of the sheet. Chebyshev’s projection is the one that spreads the residual most evenly, and evenness is what makes a single tolerance describable at all.

What is deliberately not computed here

The sheets’ shape. A real series uses rectangles in some grid, because sheets have to be printed, stored, indexed and joined, and a disc does none of those things. The covering density used here is the one a disc arrangement achieves and a rectangular scheme’s is different; the calculation is a statement about the projection’s limit rather than about the cartographer’s layout.

The overlap. Real sheets overlap deliberately so that a feature near an edge appears whole somewhere, which adds to the count by a factor that is a policy rather than a geometry.

Anything about content. The number of sheets a survey can afford is a question about money and staff, and the only thing this arithmetic contributes to it is the exponent — which is exactly the part that is not negotiable.

Interruption. A world map may be cut into lobes instead of into sheets, which reduces the distortion by a measurable factor and pays in continuity rather than in count. Giving up continuity measures that trade, and the two are alternatives: a lobed map is a covering whose pieces are joined along some edges and not others, and the sheet count above is the limit of cutting everywhere.

What a screen map does with the same arithmetic

A tiling scheme is a covering too, and it answers the question the other way round: it fixes the sheets first and lets the tolerance fall where it may.

A tile at zoom zz covers 1/4z1/4^z of a square world, so the count is a power of four and the sheet radius halves at every level. Substituting into the law above, the scale tolerance a tile achieves improves by a factor of four per level — which is why a screen map is a pyramid of tiles never states a tolerance: it has a different one at every level, and the reader chooses the level.

That is a genuinely different design and it is worth naming the trade. A paper series picks a tolerance and pays a sheet count. A pyramid picks a sheet count per level and inherits whatever tolerance the level gives, which is why the top of the pyramid is a map nobody could measure anything on and the bottom is one nobody could print.

The law is linear, and that is the cheap direction

The bound has a shape worth reading on its own, because it says which of an atlas’s two demands is the expensive one.

The count runs as 1/τ1/\tau: one thousand and one sheets at one part in a thousand, ten thousand and one at one part in ten thousand, one hundred and one at one part in a hundred. Halving the tolerance doubles the sheet count and does nothing worse than that.

Now set that beside the other demand a series can make. A sheet of fixed paper size at scale ss covers ground area proportional to 1/s21/s^2, so doubling the scale — 1:50,000 to 1:25,000 — quadruples the number of sheets. That is the familiar arithmetic of a national series, and it is the reason a large-scale series runs to thousands of sheets while a road atlas runs to dozens.

So the two demands have different exponents, and the accuracy demand has the smaller one. Asking for twice the fidelity costs half as much as asking for twice the detail, measured in sheets, and the gap widens at every further doubling: four steps of tolerance is a factor of 16, four steps of scale is a factor of 256.

That inverts the intuition the words carry. Twice as accurate sounds like the severe request and twice as large a scale sounds like a routine one, and in sheet count it is the other way round.

The reason is that the two are demands on different things. Scale is a demand on how much of the sheet the ground occupies, which is an area, so it scales as a square. Tolerance is a demand on the scale factor’s spread over a cap, and the spread goes as the square of the cap’s radius — so inverting it gives a radius proportional to τ\sqrt{\tau}, an area proportional to τ\tau, and a count proportional to 1/τ1/\tau. The square that appears in one law is the same square that cancels in the other.

None of this makes 1,001 sheets a small number. It does mean that a series which has already accepted a sheet count for its scale has room to be much more accurate than it usually is, at a price it has already shown itself willing to pay for something else.

Where this ladder goes

The choosing field has now measured what a projection minimises, what a compromise buys, what a region does to a ranking, what a line does to it, what an aspect is worth, what interruption costs, and — here — what a tolerance costs before any of those are decided.

What remains is the case where the region is not a cap and not a country but a set of regions with different requirements, which is the atlas problem proper: not how many sheets, but which sheets, at which scales, with which projections, chosen together. That is an optimisation with a discrete part, and this site does not have it.

What it does have now is the floor: no arrangement of any conformal projections covers the world at one part in a thousand in fewer than 1,001 sheets, and the number is a theorem rather than an estimate.

Named alongside this one

Essays reaching for the same objects. Nobody chose these; they are what the concept index makes visible.

What links here

Every essay whose body links to this one.

The objects this essay names

Each one links to every other essay that touches it.

Chebyshev's boundClosed formCoveringLower boundOptimal conformalProjection selectionPurposeRegional distortionRepresentative fractionSheet layoutSpherical capTolerance