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

An operation on stored geometry, drawn

The subject is an operation rather than a place, so most modes show two answers to one question and the distance between them: area-methods, nearest-query, kernel-order, centroid-planes, inside-band, segment-departure, departure-scaling, warp-loss and reprojection each run one operation two ways. coordinate-readings, swap-band, ground-ellipse and seam-failure are the readings and the failures that come of performing it in the wrong system, and operation-table is the whole set at once.

The subject is an operation rather than a place, so most modes show two answers to one question and the distance between them: area-methods, nearest-query, kernel-order, centroid-planes, inside-band, segment-departure, departure-scaling, warp-loss and reprojection each run one operation two ways. coordinate-readings, swap-band, ground-ellipse and seam-failure are the readings and the failures that come of performing it in the wrong system, and operation-table is the whole set at once.

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 14.

27 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.

53° north, 2° west: four places one pair of numbers could be. The same two numbers, read as three datums and as the pair in the other order. The centre is the WGS84 reading; the marks are where the ground is if the numbers were measured against OSGB36, ED50, NAD27 instead, at 103 m, 138 m, 195 m. Every one of those is a plausible position — near enough to look like a survey difference and far enough to be a different field. Reversing the two numbers instead moves the point 7942 kilometres, and is the mistake nobody worries about because it is usually obvious.

coordinate-readings

53° north, 2° west: four places one pair of numbers could be Where longitude and latitude can be exchanged without noticing, to 1 km. An equirectangular frame with the line λ = φ drawn on it. Reversing the two numbers of a coordinate is the most common coordinate mistake there is and is usually called obvious; it is obvious everywhere except along this line, where it does exactly nothing. The marked points are 0°E 60°N, 40°E 40°N, 70°W 40°N, displaced by 8390 km, 0 km, 14872 km. The band inside which the mistake stays under 1 kilometre is 0.013° wide at 10° and 0.016° at 60° — narrow, and not zero.

swap-band

Where longitude and latitude can be exchanged without noticing, to 1 km "Within 0.01°" on the ground, at five latitudes. The same condition drawn at five latitudes, all at one scale. A degree of latitude is a fixed distance — 1106 metres here, varying by less than one per cent from equator to pole — while a degree of longitude collapses as cos φ, from 1113 metres to 289 at 75°. So the "circle" is an ellipse of 3.86:1 there, and it encloses 26 per cent of the ground the same condition covers on the equator. Even on the equator it is not round: M is smaller than N by the flattening, so the shape is 6694 parts per million shorter north–south than east–west.

ground-ellipse

"Within 0.01°" on the ground, at five latitudes The area of one 20° × 10° cell at 50–60° north, seven ways. The cell has an exact area — R²Δλ(sin φ₂ − sin φ₁), 1,416,580 square kilometres — so every other row is a measurement of the method rather than of the ground. The equal-area projection returns it to 1.000000 and the spherical polygon formula to 1.000000, which is three routes agreeing — and the same cell integrated on the ELLIPSOID comes out 0.45 per cent away from all three, because the sphere is a model. Taking the shoelace in Mercator gives 3.06 times too much, and treating degrees as a length gives 1.75 times — about sec φ at the cell's middle, which is where that error comes from.

area-methods

The area of one 20° × 10° cell at 50–60° north, seven ways A straight line stored in Web Mercator, and where it really goes. Two points 5570 kilometres apart, joined by a straight segment in the plane the file's coordinates are in. Drawn on that plane it is a straight line and looks like the route; on the ground it is the curve marked measured, and the geodesic between the same two endpoints is the other one. The profile below is the distance between them along the line, reaching 718 kilometres at 49 per cent of the way across. Nothing here is an error in the data: both endpoints are exact, and the whole of the discrepancy is the word straight.

segment-departure

A straight line stored in Web Mercator, and where it really goes The departure against length, on Web Mercator at 45°. six east–west segments from 50 to 1600 kilometres, each stored as two points and drawn straight, with the worst distance from the geodesic measured for each. On logarithmic axes the points lie on a line of slope 2.001 — the departure of a chord from an arc goes as κL²/8, so two is the prediction and the fit is the test of it. A segment of 50 km is 49 metres out; sixteen times the length is 1029 times the error rather than sixteen.

departure-scaling

The departure against length, on Web Mercator at 45° A 20° shape across the antimeridian, in the space where the numbers live. Longitude runs across the page from −180° to 180°, which is where the failure is: the shape is one rectangle on the ground and two pieces in the numbering, and every operation that treats longitude as a real number sees the two. The bounding box comes out 359° wide instead of 20°, the planar area comes out 17 times too large because the shoelace encloses the complement, and the midpoint of a segment from one edge to the other lands 20015 kilometres away — the antipode of where it belongs. The true area, from the closed form, is 4,920,667 square kilometres.

seam-failure

A 20° shape across the antimeridian, in the space where the numbers live 5 of 400 nearest-neighbour queries change answer in the plane. 40 sites and 400 queries over -20° to 40° east and 35° to 70° north, each query answered twice — once by geodesic distance and once by straight-line distance in the stored plane. They agree 98.8 per cent of the time, which is why the operation survives, and the 5 that differ are marked. The mechanism is not that the plane is wrong by a lot but that its scale factor varies: over this region it spans 139 per cent, and every disagreement is a contest closer than that — the worst margin measured is 6.4 per cent. The furthest a wrong answer is from the right one is 24 kilometres.

nearest-query

5 of 400 nearest-neighbour queries change answer in the plane A raster warped to Lambert azimuthal equal-area and back, nearest against bilinear. The left panel is the field the raster carries — a smooth analytic function, so that the error of an interpolation is the interpolation's error and not a photograph's history. The other panels are what is left after warping into Lambert azimuthal equal-area and back to Equirectangular, shown as the difference from the original at six times the contrast. Nothing moved: the coordinates go through the maps exactly. What is lost is that a target pixel's centre does not fall on a source pixel's centre, so a value has to be invented for it. Bilinear is closer to the field — RMS 0.0022 against 0.0212 — and has given up 0.63 per cent of its variance to get there. The panels are drawn at 48 by 32 cells; the measurement is made at the same resolution.

warp-loss

A raster warped to Lambert azimuthal equal-area and back, nearest against bilinear The order of three resampling kernels, measured. Each kernel run over five rasters from 48×32 to 160×110 cells, warped into Lambert azimuthal equal-area and back, with the RMS departure from the field measured on the grid it started on. On logarithmic axes each kernel is a straight line and its slope is its order of convergence: nearest 1.00, bilinear 1.98, cubic 2.93. That is the difference between them stated as a number rather than as an adjective — halving the cell size halves nearest's error, quarters bilinear's, and divides the cubic's by eight.

kernel-order

The order of three resampling kernels, measured four planes a dataset might be stored in, scored on three operations. Each candidate measured over -10° to 30° east and 35° to 60° north: the worst areal error, the worst angular deformation, and the spread of the scale factor, which are what an area query, a shape and a distance respectively depend on. The best plane for areas is Gall–Peters, for shapes Web Mercator, and for distances Lambert azimuthal equal-area — three different answers, and no fourth candidate would collapse them, because a projection exact in two of these columns has a = b = 1 everywhere and is the isometry Gauss's theorem forbids. area of a polygon costs 3.06× too large in the wrong plane; drawing a line between two points costs 194 km from the ground it claims.

operation-table

four planes a dataset might be stored in, scored on three operations Where the middle of a 20° × 20° region is, in five planes. The region is drawn in longitude and latitude — which is itself a projection, and one of the ones being compared. Each filled mark is the shoelace centroid computed in one projected plane and inverted back to the ground; the hollow mark is the centre of area on the sphere, by integration. They spread over 273 kilometres. The equal-area member is 66 kilometres out, because a centroid is a first moment and preserving area says nothing about where the area sits.

centroid-planes

Where the middle of a 20° × 20° region is, in five planes Two readings of the same triangle, and the band between them. The same three vertices, joined two ways: with straight lines in the plane the coordinates are stored in, and along the ground. Every dot is a point the two readings disagree about — inside on one and outside on the other. The band covers 29.5 per cent of the polygon, which is 4241 thousand square kilometres, and it is not an error in either reading: the file does not say which one it means. Densifying the stored boundary removes it, which is the only fix there is. Drawn in Equirectangular.

inside-band

Two readings of the same triangle, and the band between them What reprojecting from Mollweide to Mercator does, drawn on the Mollweide plane. Every dot is a point of the source picture, sized by the angular deformation the reprojection itself applies there — not the deformation of either projection against the sphere, but of the map that carries one plane to the other. It reaches 116.5° here. This is the quantity a reprojection actually costs a shape, and it is computed from the two Jacobians alone: the sphere's own metric cancels exactly, so the transformation does not depend on the shape of the Earth at all. Drawn on the source plane, in Mollweide.

reprojection

What reprojecting from Mollweide to Mercator does, drawn on the Mollweide plane

Where it is called

Changing this changes every one of these figures.

What 4 parts per million of network strain leaves behind. Each arrow is where the best-fitting seven-parameter transformation leaves a marker, over 36 markers laid out as a grid across OSGB36's ground. The RMS residual is 1.64 metres and the worst is 3.00 metres. Seven parameters span the constant and linear parts of a displacement field; this one is quadratic, so no choice of the seven can reach it. The arrows are drawn 27× life size. What the numbers refer to

Where a fit leaves residuals

Seven parameters can carry a rigid motion and a size exactly. A triangulation network is neither, so the best possible transformation between two datums leaves metres on the table — in a pattern, not as noise — and which seven parameters come out depends on where the markers were.

 One pixel at zoom 11: 76.4 m projected, 47.6 m at 51.5°. The grid squares are pixels, at their own size. The open mark is the stored coordinate and the filled one is where it is drawn: rounding moves it 20.2 metres, against a worst case of 33.6 — half a pixel's diagonal, which is the whole of the bound. The second point sits 43 metres away, 0.58 of a pixel east and 0.70 north, so whether the two are drawn as one dot or two is decided by where the tile grid happens to fall: they merge for 12 per cent of the possible offsets, against the 12 per cent the two fractions predict. What a machine does with it

The pixel is a place with a size

Drawing a coordinate rounds it to a pixel, which moves it by up to half a diagonal — 16.8 metres at zoom 12. Whether two points 43 metres apart appear as two dots is not a property of the data at all: they merge for 12 per cent of the positions the tile grid could take, and the closed form predicts 12.4.

 Five latitudes that are not the latitude, on WGS84. Each curve is the amount by which one auxiliary latitude falls below the geodetic latitude a coordinate actually means, in arcminutes. All five vanish at the equator and at the poles and peak near 45°, where the geocentric latitude is 11.55 arcminutes below the geodetic one — about 21.4 km on the ground. The curves never cross, which is forced by the algebra rather than by this ellipsoid's particular flattening. What is taught wrongly

Geodetic against geocentric latitude

Two angles, both called latitude, differing by eleven and a half arcminutes at their worst. One is what every coordinate means and the other is what every spherical formula assumes, and the gap between them is twenty-one kilometres on the ground.

 London to Tokyo in three straight legs. The great circle, and the route a plan of three constant-heading legs actually follows between waypoints on it. The two touch at the waypoints and part between them by up to 456 km, and the flown route is 253 km longer than the direct one. The headings are 50°, 109°, 149°. Two routes. The great circle is 9559 km and is the shortest path on the sphere. The rhumb line holds a single compass bearing the whole way and is 11296 km — 1737 km further, or 18.2 per cent. On Orthographic the rhumb line departs from straight by 2.4e-1 of its own length. Paths and directions

Flying a curve in straight legs

Nobody steers a great circle, because a great circle requires the heading to change continuously. What is actually flown is a handful of constant-heading legs between waypoints on it, and the gap between plan and curve falls as the square of the number of legs.

 53° north, 2° west: four places one pair of numbers could be. The same two numbers, read as three datums and as the pair in the other order. The centre is the WGS84 reading; the marks are where the ground is if the numbers were measured against OSGB36, ED50, NAD27 instead, at 103 m, 138 m, 195 m. Every one of those is a plausible position — near enough to look like a survey difference and far enough to be a different field. Reversing the two numbers instead moves the point 7942 kilometres, and is the mistake nobody worries about because it is usually obvious. What a machine does with it

A coordinate without its system is not a location

The same two numbers name ground 103 metres apart on one datum, 195 on another, and 7,942 kilometres away if the pair is read in the other order. The datum errors are the dangerous ones because they are plausible — and the axis swap, which everyone calls obvious, does nothing at all along a line that crosses three continents.

 "Within 0.01°" on the ground, at five latitudes. The same condition drawn at five latitudes, all at one scale. A degree of latitude is a fixed distance — 1106 metres here, varying by less than one per cent from equator to pole — while a degree of longitude collapses as cos φ, from 1113 metres to 289 at 75°. So the "circle" is an ellipse of 3.86:1 there, and it encloses 26 per cent of the ground the same condition covers on the equator. Even on the equator it is not round: M is smaller than N by the flattening, so the shape is 6694 parts per million shorter north–south than east–west. What a machine does with it

A degree is not a unit of length

"Within 0.01 degrees" is a condition anybody can write and no instrument can measure. On the ground it is an ellipse — 1,106 metres north–south and 558 east–west at 60° — and even on the equator it is not a circle, because the meridian's radius of curvature is smaller than the parallel's by the flattening.

 The area of one 20° × 10° cell at 50–60° north, seven ways. The cell has an exact area — R²Δλ(sin φ₂ − sin φ₁), 1,416,580 square kilometres — so every other row is a measurement of the method rather than of the ground. The equal-area projection returns it to 1.000000 and the spherical polygon formula to 1.000000, which is three routes agreeing — and the same cell integrated on the ELLIPSOID comes out 0.45 per cent away from all three, because the sphere is a model. Taking the shoelace in Mercator gives 3.06 times too much, and treating degrees as a length gives 1.75 times — about sec φ at the cell's middle, which is where that error comes from. What a machine does with it

Computing an area needs a surface

A shoelace over a ring of coordinates returns a number whatever the coordinates are. For one twenty-by-ten-degree cell it returns 1.75 times the true area in degrees, 3.06 in a conformal plane, and exactly the closed form in an equal-area one — and the closed form itself is 0.45 per cent out, because the sphere is a model too.

 Three curves between 45°N and 55°N, 3026 km apart. The normal section observed from the first point, the normal section observed from the second, and the geodesic, each plotted as its distance to one side of the great-circle chord between the two ends. The two sections are 79.0 metres apart at their widest and the geodesic runs between them, 53.3 metres from the first. In LENGTH they differ by almost nothing — 2.40 millimetres over 3026 kilometres — so an observation that follows the wrong one measures the right distance along the wrong ground. On a sphere all three coincide exactly. What is taught wrongly

The normal section is not the geodesic

A theodolite at A sighted on B swings in one plane and the instrument at B sighted back swings in another, so the two observations trace different curves on the ground — 79 metres apart over 3,026 kilometres — and the shortest path is neither of them. The lengths differ by 2.4 millimetres, so the wrong curve measures the right distance along the wrong ground.

The whole library