Bodies — the ladder · Cartographic Projection

Two latitudes for one point on Mars. A cross-section of Mars with its flattening exaggerated 6× so the two angles can be told apart. The line to the centre defines the planetocentric latitude and the outward normal defines the planetographic one; at 45° they differ by 0.338° on the real body, which is 20 kilometres along the surface. Both numbers are published for Mars, and a coordinate that does not say which it is is ambiguous by that much.

A coordinate on another body

Mars publishes two latitudes for every point and they differ by up to 0.338°, which is twenty kilometres of ground — almost exactly the same distance as the Earth's own 0.192°, because Mars is smaller by nearly the same factor. The libraries that compute either have known both numbers since this collection's early essays and had never been pointed anywhere but here.

rung 1 · datums

The error of a spherical formula is twice the flattening. The worst angular deformation of the spherical Mercator formulae applied to each body's own latitudes, against that body's flattening, on logarithmic axes. The points are measured by differencing the projection; the line is 2f radians, which is not fitted. The two agree to a third of a per cent for Mercury, the Earth and Mars, and depart by 3.3 and 5.1 per cent for Jupiter and Saturn, whose flattenings are too large for the first term to be the whole story. The site's headline number — 0.3848° on Earth — is an instance of this law rather than a fact about Web Mercator.

The same projection on a different body

A projection's distortion does not depend on the body's size at all — only on its flattening — so Web Mercator's 0.3848° of angular deformation is not a fact about Web Mercator. It is 2f radians, and the same mistake made on Mars measures 0.6765°, on Jupiter 7.68°, and on a body scaled ten times larger than Mars the identical number to twelve figures.

rung 2 · datums

On a triaxial body, latitude depends on longitude. Walk round each body at a constant planetocentric latitude of 45° and watch the direction of the surface normal, which is what the planetographic latitude is. On Mars it does not move at all — that is what having an axis of revolution means. On Vesta it swings by 1.40° and on Phobos it swings by 6.40°. A body without an axis has no latitude that is a function of position alone, and every coordinate on it is a convention with a body-fixed frame attached.

A body that is not an ellipsoid

Vesta's three axes are 286.3, 278.6 and 223.2 kilometres, all different, so it has no axis of revolution — and on such a body the planetographic latitude of a point at 45° planetocentric swings by 1.4° as one walks round it in longitude, and by 6.4° on Phobos. Latitude stops being a function of position.

rung 3 · datums

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.

rung 4 · datums

Shape predicted from field, against shape as published. Clairaut's theorem gives a body's flattening from two numbers of its gravity field: J₂, which is how its mass is arranged, and m = ω²a³/GM, which is how fast it spins. For a body in hydrostatic equilibrium the prediction is the shape, and the diagonal is where such a body sits. Earth is on it to 0.05 per cent — 12 metres at the pole, out of twenty-one kilometres of flattening. Mars is 12.5 per cent off it, which is 2.23 kilometres, and the excess is Tharsis: a body carrying a continent-sized volcanic load is not a fluid figure, so its ellipsoid is not one of its own level surfaces, and its zero of height has to be chosen rather than found.

A body with no sea level

On Earth the zero of height is found rather than chosen — water settles onto the equipotential surface by itself. Nowhere else has one, and the difference is measurable: Clairaut's theorem predicts the Earth's flattening from its own gravity field to twelve metres at the pole and misses Mars's by 2.2 kilometres.

rung 5 · datums

A conformal map of Vesta, in the coordinates that make one possible. The body's own isothermal coordinate, used as the map — which is Mercator's construction carried out on a body that has no axis to be Mercator about. The lines are Jacobi's ellipsoidal coordinates, which are the surface's lines of curvature, at 30° and 20° spacing; the two quadratures that turn them into an isothermal pair are one-dimensional. The measured angular deformation away from the marked points is 1.1e-6°, which is this site's noise floor, and the areal factor spans a factor of 162 — conformal, and emphatically not equal-area. The four marks are the umbilics, where the coordinate collapses and the map has nothing to say.

A conformal map of a body with three axes

This site built an equal-area map of a triaxial body and wrote down what it could not do: the conformal one, which needs an isothermal coordinate that a surface with no axis of revolution was said not to have. It has one, Jacobi found it in 1839 for a different reason, and it takes two one-dimensional integrals.

rung 6 · datums

A conformal map solved rather than written down. The same patch of parameters on two bodies, mapped to the plane by solving the discrete Cauchy–Riemann equations — one complex equation per triangle, 1568 triangles, least squares, conjugate gradients, and no formula for either surface. Left: a sphere, where the answer is known in closed form and is not used. Right: a body with a bump on it, which has no isothermal coordinate and therefore no closed form at all. The parameter lines cross at right angles in both, to a median of 0.60° and 1.25° of angular deformation.

A conformal map of a body that is not a quadric

Jacobi's ellipsoidal coordinates give a triaxial body a conformal map by two quadratures, and this collection wrote down what that argument uses: the surface has to be a quadric. A real body is not. Solving the discrete Cauchy–Riemann equations instead — one complex equation per triangle, two thousand triangles, conjugate gradients — gives a conformal map of a bumped body to a median of 1.10°, converging at first order in the mesh, with the areal factor spreading by 3.07 and refusing to converge at all.

rung 7 · datums

The same patch, pinned six ways. Two vertices have to be held or the conformal energy has a similarity's worth of null space. Which two turns out to decide two of the three numbers reported. The median angular deformation is the same to 8 per cent across all six — that is the map. The areal spread runs from 2.68 to 7.44, so the 3.07 reported for this body was a statement about its corners. And pinning two adjacent vertices, which fixes the similarity through a very short lever, ruins the worst point without touching the median: a badly conditioned constraint pays for its scale in one corner.

The map depends on where it was cut

The previous rung solved the discrete conformal equations on a triangulated body and reported an areal spread of 3.07, then recorded that the number might belong to the patch, the boundary and the two pinned vertices rather than to the surface. It belongs to the pins: hold a different pair and it runs from 2.68 to 7.44, while the typical angular deformation does not move at all.

rung 8 · datums

Cuts of the same size in different places

Where a body is cut decides how well it can be mapped, by a factor of ten — established with five windows of five different sizes, so *where* and *how much* were confounded and the factor could have been entirely about extent. Held to the same surface area to a quarter of a per cent, the answer survives at a factor of 1.32, and what predicts it is the curvature the window encloses.

rung 9 · datums

The solve's cost is a U in the shape; the curvature is not. Solid: how many conjugate-gradient steps the conformal solve needs, against the window's aspect ratio, at constant surface area. It has a minimum at the square window — 174 steps — and rises to 265 and 275 at the two extremes, which are elongated by the same factor in opposite directions. Dashed: the total Gaussian curvature the window encloses, on its own scale, which falls from 0.311 to -0.012 across the same sweep and is least at one end of it. The cost has its minimum where the window is square and the curvature has its minimum somewhere else, so whatever is making the solve expensive is not what is making the map spread.

A long window and a square one

Rung nine held the patch's area and found that where the cut goes still changes the map by a third, with the curvature it encloses predicting the change at r = 0.969. It recorded that it had held area and curvature and not shape. Sweeping the shape at constant area separates two things that had looked like one: the map is curvature and the cost is shape.

rung 10 · datums

An error in a rotation rate is a longitude error that never stops growing. Where the prime meridian of each body has got to, if its published rotation rate is wrong by one unit in its last published decimal. Every line is straight through the origin, because the error is the rate error multiplied by the elapsed time and nothing else — there is no date after which it settles. Jupiter reaches 2279 metres at the equator after a century; Mars, whose rate is published to twelve decimals rather than seven, reaches 0.0011 metres over the same interval.

A longitude that drifts with the rotation rate

Ten essays here map bodies whose shape is the problem. A longitude is not about shape: it is a landmark plus an extrapolation over however many days have passed, and an error in the last published decimal of a rotation rate is a coordinate error that grows without bound in time.

rung 11 · datums

A section through a body with a neck, and the rays that leave it twice. An equatorial section of a stated contact binary — two lobes of radius 1 centred at ±1.5, joined by a neck of radius 0.35, with the origin in the neck. The lines are rays from the origin and the marks are where each one crosses the surface. two of the 25 drawn cross more than once, so along those directions there is no such thing as "the" radius, and a longitude and a latitude do not name a place.

A ray from the centre hits the surface twice

Eleven rungs map bodies that are lumpy, triaxial and turning at a drifting rate, and every one assumes the surface is star-shaped about the centre — which is what makes a longitude and a latitude a coordinate at all. A contact binary is not: on a stated body with a neck a third of a lobe wide, 10.9 per cent of the sky has no single radius, and the shape model everybody publishes fills the neck in and adds 1.67 per cent of the volume.

rung 12 · datums

A surface whose curvature changes sign, and integrates to nothing. a wide ring of major radius 3 and minor radius 1, shaded by its Gaussian curvature. The outer half is positively curved like a sphere, at up to 0.250; the inner half is saddle-shaped and negative, down to -0.500; and the two circles between them, drawn as lines, are exactly flat. The integral of the curvature over the whole surface is zero, which is 2π times the Euler characteristic — and every impossibility in this collection rests on that number being 2 rather than 0.

On a body with a hole, north can be up everywhere

Twelve rungs vary the body's shape and none varies its topology, which is what every impossibility here actually rests on. A torus has a nowhere-zero tangent field, a total curvature of zero rather than 4π, and a conformal world map with no cut and no singular point — and it still cannot be flattened, for the one reason that survives.

rung 13 · datums