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

Metric — where it appears

The rule that assigns a length to a path on a surface, and the thing a distance is measured against before any projection is chosen. The metric is a property of the body and the projection is a property of the page, and every distortion quantity here is the ratio between them.

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

Why there are two latitudes, at 45°. A meridian section of an ellipsoid with the flattening exaggerated 12× — at the true value of 1/298 this outline would be indistinguishable from a circle. Geodetic latitude is the angle the surface normal makes with the equatorial plane; geocentric latitude is the angle the radius makes with it. The two lines meet the plane at different points and the angles differ by 11.55 arcminutes on the real WGS84 ellipsoid, which is about 21 km of ground distance.

The Earth is a sphere, and when it is not

Every essay before this one treated the Earth as a ball, and said so. The flattening is one part in three hundred, which is nothing for a distance, everything for a latitude, and exactly enough to make the most-used projection in the world fail the property in its own name.

wrong · Ellipsoid Which of eight places is nearest, decided on the ground and decided on the page. Every cell of the window shaded by which of eight places across Europe is nearest on the ground, with the cells the page's own answer would hand to a different site drawn in the failure colour. Plate carrée misassigns 13.67 per cent of the window's ground area — 3,058 thousand square kilometres. The misassigned cells are not scattered: they lie in bands along the boundaries, which is what a systematic error looks like and what a sampled test of a few query points is least likely to find.

Nearest of many is a partition

One reach question with one source is a disc. With several sources it is a division of the whole surface, every place belonging to whichever source is nearest — and computing that division in the plane the data is stored in hands away between 0.75 and 22.16 per cent of the ground, in unbroken strips up to 1,591 kilometres across.

paths · Reach 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 is a question about the metric

Asked in the plane the data is stored in, a nearest-neighbour query returns a different site for 5 of 400 queries over Europe on Web Mercator — and 103 of 400 on the plate carrée, with the wrong answer up to 309 kilometres further away. Every disagreement is a contest closer than the region's own scale spread, and over a city there are none.

applied · Dataset five 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 Lambert conformal conic, and for distances Lambert conformal conic — 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.

The operation decides the coordinate system

Five candidate planes over one region, scored on the three things a spatial operation depends on. The conformal conic wins shape and distance and is 11.7 per cent out on area; the equal-area member is exact on area and 38.9° out on shape. No candidate is exact in two columns, and no candidate ever will be, because one that was would be an isometry.

applied · Dataset The circle a reprojection is about to turn into an ellipse. A circle drawn on the Mollweide plane, and its image under the map that carries that plane to Mercator — normalised so each pair has the same area, because the whole map can be rescaled and the shape is what is being shown. The dashed circle is what an undistorted reprojection would leave. This is Tissot's construction with a plane in place of the sphere, and it is the right picture for a reprojection because both ends are pictures. Worst angular deformation over the sampled points: 81.2°. Drawn on the source plane, in Mollweide.

A projection between two projections

Every distortion measured on this site so far compares a map with the sphere. The operation a machine actually performs compares a map with another map — and that map has its own two principal scales, its own areal factor and its own angular deformation, none of which is the difference of the two it was built from.

distortion · Tissot The set that can be reached is not the set that can reach. Two sets under a steady westerly: everywhere reachable from the marked place in a stated time, and everywhere from which the marked place can be reached in the same time. They have the same area to a fraction of a per cent — reversing a uniform flow is a reflection — and they overlap on only 32.3 per cent of their union. 67.7 per cent of the ground in one of them is not in the other.

The set that can be reached is not the set that can reach

The moment a cost stops being symmetric, two questions that read alike stop having the same answer. Under a flow at 45 per cent of a vehicle's own speed the set reachable from a place and the set from which the place is reachable have the same area to five significant figures and share 32 per cent of their union — so 68 per cent of the ground in one of them is not in the other.

paths · Reach 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.

A centroid belongs to a plane

Every renderer labels a region at its centroid, and every centroid is a shoelace over coordinates as stored — which is a statement about the plane they are in. Six planes put the middle of one 20° × 20° region up to 273 kilometres apart, the equal-area member is 66 kilometres out, and the disagreement falls as the square of the region's size.

applied · Dataset The query, the box it becomes, and the 14 per cent it loses. A query for everything within 400 kilometres of a place at 55° north, drawn in the stored coordinates of Web Mercator. The ring is the true answer's edge; the rectangle is the box the index is given, sized from the scale factor at the query point. The arcs outside the rectangle are true neighbours the index never returns — 13.6 per cent of the rim, reaching 5.0 per cent of the box's own half-width beyond it — and nothing downstream can tell they are missing.

The query a fast path actually answers

Fifteen rungs price what a stored coordinate means and none asks what it is searched with. No index answers "within two hundred kilometres of here"; an index answers "inside this rectangle of stored coordinates", and the rectangle is built by somebody's arithmetic — which at 55° north silently drops a tenth of the true answer on a conformal projection and three fifths of it on an equal-area one.

applied · Dataset

Named alongside it

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

VerificationConformalityPurposeEqual-areaConventionCoordinate reference systemEquidistantGeodesicNearest-neighbourPartitionRegional distortionScale factor

All concepts