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

Auxiliary latitude — where it appears

A latitude-like variable derived from the geodetic latitude so that some spherical formula becomes exact on the ellipsoid. The conformal one preserves angles, the authalic one preserves areas and the rectifying one preserves meridian distance, and no single substitution does more than one.

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

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.

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.

wrong · Ellipsoid The radius of equal curvature, from equator to pole. The radius of the sphere that has the same Gaussian curvature as the ellipsoid does, at each latitude. On WGS84 it runs from 6357 km at the equator to 6400 km at the pole — the polar region is the flattened part of a squashed ball and is therefore the LEAST curved — so K itself varies by 1.35%. A surface whose curvature varies cannot be laid on one whose curvature does not, so no map from the sphere to the ellipsoid is faithful either, and the least a conformal one can vary in scale is 6739 parts per million.

The curvature of the Earth is not one number

An ellipsoid's Gaussian curvature varies by 1.35% from equator to pole, so the sphere is not developable onto the ellipsoid any more than the plane is onto the sphere. The spherical approximation is a projection with an irreducible cost, and the cost is 6,739 parts per million.

impossibility · Curvature Four ways to build the same equal-area projection. The cylindrical equal-area projection's own areal scale factor, measured from its Jacobian against the metric of the body it is drawn for. On a sphere with the spherical formula it is one everywhere, which is the control. Feed the same formula the geodetic latitude a coordinate actually carries and measure against the ellipsoid it refers to, and it is 1.00674 on the equator and 0.99332 at 88° — a spread of 1.34 per cent, on a projection whose entire purpose is that there is no spread. Rescaling to make the totals agree does not repair it. The authalic northing q/2 does, exactly.

Equal-area on the wrong body

The site's headline is that Web Mercator puts geodetic latitudes into a spherical conformal formula and stops being conformal. The same sentence is true with "equal-area" in it and nobody says it: the areal factor is 1.00674 at the equator, 0.99332 at 88°, and averages to almost exactly one — so every check that adds up areas passes while every cell is wrong.

wrong · Audit A deflection of 10 arcseconds, and what it hides. The ellipsoid normal and the plumb line at one point, with the geoid tilted against the ellipsoid by 10 arcseconds — drawn 3000× steeper than life, because at true scale the two lines are indistinguishable. The relation is exact and linear: an arcsecond of deflection is the geoid rising 4.85 millimetres in a kilometre, so 10 arcseconds is 48.5 millimetres per kilometre. A star sight measures the plumb line's direction, so astronomic latitude differs from geodetic by exactly this angle — 309 metres of ground, at a point where the coordinate itself is correct.

The plumb line is not the normal

A latitude measured from the stars and a latitude that means a position on the ellipsoid are different angles, because a plumb bob hangs along gravity and gravity is not perpendicular to a mathematical surface. Ten arcseconds of difference is 309 metres of ground.

datums · Height The equator, with its ellipticity exaggerated forty thousand times. The dashed circle is the equator every projection formula on this site assumes. The solid curve is the equator satellite geodesy reports, drawn with its departure multiplied by 40,000 so that seventy metres on a six-thousand-kilometre radius can be seen at all. The long axis is at 14.9° west and the short one ninety degrees from it, and the difference between them is 70.0 metres — a real quantity, about the height of a twenty-storey building, on a body every geodetic computation treats as a surface of revolution.

The equator is not a circle either

Eleven rungs price what pretending the Earth is a sphere costs, and every one of them replaces the sphere with a surface of revolution — a body whose equator is a circle. It is not. The two equatorial radii differ by seventy metres, the two surfaces part by thirty-five, and the auxiliary latitudes every ellipsoidal formula is written in stop existing.

wrong · Ellipsoid Going out and coming back, at three orders. A point three degrees from the central meridian, taken forward into the projection and back again with the series truncated at the same order both ways, and the ground distance between where it started and where it returned. The second-order pair is out by a tenth of a metre at some latitudes; the third by half a millimetre; the fourth — the order every national grid formula in ordinary use is written to — by 0.39 mm. The dips are where one of the two series passes through a node, not where the map is better.

The inverse of the series is not the series of the inverse

Transverse Mercator on an ellipsoid has no closed form, so every national grid computes it as a truncated series. There are two of them — one out and one back — and they are separate approximations. Composing them does not give the identity: at the order every grid formula in ordinary use is written to, a point three degrees from the central meridian comes back 0.39 mm north of where it started.

families · Ellipsoid Vesta, with its own coordinate lines. Vesta drawn in an orthographic projection — every figure here is a projection, including this one — with the parametric coordinate lines its published coordinates are written in. Its three semi-axes are 286.3, 278.6 and 223.2 km, so the equator is an ellipse rather than a circle and the body has no axis of revolution to define a latitude against. The surface's own radius runs from 223.2 to 286.2 km, a ratio of 1.282, and the outward normal departs from the direction out of the centre by up to 14.10°.

A map of a body with three axes

Every projection on this site is written in a latitude, and a body with three unequal axes has none: the coordinate lines are not perpendicular, Lambert's equal-area formula spreads areas by 28 per cent on Vesta, and the equal-area map that does work has a different one on every meridian.

datums · Bodies Four radii of the Earth, and one that is a range. The four constants called the mean radius of the Earth, on a scale of kilometres, with the range of the local Gaussian radius √(MN) drawn behind them. Three of the four agree to about a part per million; the rectifying radius is 3560 metres smaller, which is 559 parts per million. The Gaussian radius spans twelve times that range on its own, which is why there is no such thing as the conformal sphere.

Four radii of the Earth

The ellipsoid is usually handled with an auxiliary latitude. Every such construction also needs a radius, and the radius that makes each property exact is a different number: the published 6371 km is right for an area to half a part per million and wrong for a meridian distance by 559 — while the radius a conformal map needs is not a constant at all, and spans 6,739.

families · Ellipsoid Four spheres are wrong about a direction before anything has moved. At 45° of latitude, the worst error over twenty-four starting directions of the azimuth each substituted sphere gives for a geodesic on WGS84, against the distance to its end. Geodetic latitude is wrong by 347″ at one kilometre, and the parametric, authalic and rectifying latitudes by 174″, 116″ and 86.8″ — a half, a third and a quarter of it — flat until a few hundred kilometres. The conformal latitude is right at the start and wrong by 54.3 mas at one kilometre, 5.43″ at 100 and 54.3″ at 1,000; the geocentric latitude starts at 620 mas and joins it. Bessel's sphere, which also corrects the longitude, is exact at every distance to the arithmetic's floor of 1.2e-7″.

A bearing on a sphere is decided by its latitude, not its radius

A sphere that stands in for the ellipsoid cannot be wrong about a direction by being the wrong size, because enlarging a sphere moves no angle. What it can be wrong by is the latitude put on it. Geodetic latitude turns every bearing at 45° by up to 347 arcseconds before anything has moved; the conformal latitude starts exact and drifts 54 arcseconds in a thousand kilometres; and only Bessel's sphere, which changes the longitude too, is exact at every distance.

wrong · Ellipsoid An angle is a difference of two directions, and the error does not cancel between them. A geodesic triangle with sides of 30 km at 45° N, its first side laid off on a bearing of 45° and the second 60° round from it. Substituting a sphere carrying the geodetic latitude turns each direction at a corner by an amount that runs as the sine of twice its azimuth, so the two directions at A are turned by different amounts and the angle between them is wrong by their difference. Here the worst single direction is out by 345.79″ and the worst angle by 521.09″ — larger, not smaller. Bessel's sphere, which also corrects the longitude, keeps every angle to 6.2e-9″.

An angle is a difference, and the difference doubles the error

Substituting a sphere for the ellipsoid turns every direction at a point, and a surveyor measures angles rather than directions — so the obvious hope is that a turn common to both directions cancels in their difference. It does not cancel, because the turn is not common: it runs as the sine of twice the azimuth, and the largest angle error is exactly twice the largest direction error at a well-shaped corner. A first-order triangulation computed on the geodetic-latitude sphere carries angles nearly nine minutes of arc wrong, and its closure check passes perfectly.

wrong · Ellipsoid Where the correction overtakes the sphere it is correcting for. The worst angle error of four substituted spheres against the size of the triangle, with Legendre's spherical excess on the same axes. The two that are wrong at first order in the flattening are flat in the side, because the error belongs to the latitude and not to the distance. The two that are right at first order have no such floor, so what is left is a term rising in proportion to the side. The excess rises as its square, so each pair meets exactly once, and the meeting decides which of the two corrections a computation is entitled to ignore: the geodetic latitude at 490 km, the rectifying latitude at 247 km, the conformal latitude at 23.8 km, the geocentric latitude at 17.6 km.

The correction is the smaller of the two corrections

Legendre's theorem subtracts a third of the spherical excess from each angle of a triangle and solves the rest as a plane figure. It is the correction every nineteenth-century computer applied and the one every textbook explains. Beside the error the choice of substituted sphere leaves in the same triangle it is small: on the thirty-kilometre triangle a first-order chain is made of, omitting it costs 166 millimetres in a computed side and using the geodetic-latitude sphere costs 43 metres. The two terms scale differently, so each sphere has a size at which the correction overtakes it — and for the two spheres a careful computer would have chosen, that size falls inside a real chain.

wrong · Ellipsoid Along a meridian a good sphere's error adds up triangle by triangle; a bad base does not. A strip of thirty-kilometre equilateral triangles running north from 45°, each side computed from the last by Legendre's theorem on a substituted sphere. The scale error of the side just computed, against how many triangles lie between it and the measured base. The conformal sphere adds 7.90 ppm a triangle and reaches 158.1 ppm after 20; the geocentric, 158.9 ppm; the authalic, −52.4 ppm; the rectifying, 0.29 ppm. Shaded: three standard deviations of the scale error one-second angles leave, 3.96 ppm a triangle, growing as the square root of the count. Dashed: a base measured 8 ppm long, which stays 8 ppm at every triangle.

A strip of triangles carries the sphere's error to the next base

One triangle cannot tell a substituted sphere from a bad base: both leave the computed side a few parts per million long. A triangulation is not one triangle but a strip, each computed side the next one's base, and along it the two come apart completely. A bad base is one factor, the same at every triangle. A sphere's error is made afresh in each triangle and the strip adds it up, and what it adds up to at the far end is exactly the change in that sphere's own meridian scale between the two ends — 158 parts per million for the conformal sphere over twenty triangles north from 45°, predicted to half a part. Along a parallel it is nothing; at a check base three triangles on it is already three times the noise of one-second angles.

wrong · Ellipsoid Adjusted to close on the second base, a sphere's drift of hundreds of ppm becomes a bow of under one. A strip of 40 thirty-kilometre equilateral triangles north from 45°, computed on the conformal sphere, between a base at each end. Unadjusted (dashed), the side error climbs to 315.0 ppm at the second base. Adjusted to close there, each triangle's share of the misclosure equal, what is left (solid) is the sphere's stretch less a straight line between the bases: 0.63 ppm at its largest, at triangle 23. Shaded: three standard deviations of what one-second angles leave in the adjusted strip, 37.6 ppm at the middle and nothing at either base.

An adjustment hides the sphere only in a strip of equal triangles

A strip computed on a substituted sphere drifts by the sphere's own stretch, 315 parts per million over forty triangles north from 45°, and an adjustment that closes it on a second base takes almost all of that away: in a strip of equal triangles what is left is the stretch less a straight line, 0.63 parts per million, a third-order bow sixty times below the noise of the old angles. A real chain is not made of equal triangles. The adjustment spreads the misclosure by how uncertain each triangle's angles are, the sphere put it in by how far north each triangle reaches, and the difference is a first-order residual that never clears the noise at any one side and that the old network's own geometry predicts exactly — so a re-observation matched against it reads the sphere at seven standard deviations.

wrong · Ellipsoid

Named alongside it

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

Geodetic latitudeEllipsoidSpherical approximationVerificationConformal latitudeFlatteningScale factorSurface normalSurvey networkToleranceTriangulationClosed form

All concepts