Wednesday, January 14, 2026

Changing between geographic and geocentric latitude


If the earth had been an ideal sphere, “down” can be the path to the middle of the earth, wherever you stand. However as a result of our planet is a bit flattened on the poles, a line perpendicular to the floor and a line to the middle of the earth should not the identical. They’re practically the identical as a result of the earth is sort of a sphere, however not precisely, until you’re on the equator or at one of many poles. Typically the distinction issues and typically it doesn’t.

From a given level on the earth’s floor, draw two strains: one straight down (i.e. perpendicular to the floor) and one straight to the middle of the earth. The angle φ that the previous makes with the equatorial aircraft is geographic latitude. The angle θ that the latter makes with the equatorial aircraft is geocentric latitude.

For illustration we’ll draw an ellipse that’s much more eccentric than a polar cross-section of the earth.

At first it will not be clear why geographic latitude is outlined the best way it’s; geocentric latitude is conceptually easier. However geographic latitude is less complicated to measure: a plumb bob will present you which of them path is straight down.

There could also be some slight variation between the path of a plumb bob and a perpendicular to the earth’s floor resulting from variations in floor gravity. Nevertheless, the deviations resulting from gravity are a pair orders of magnitude smaller than the variations between geographic and geocentric latitude.

Conversion formulation

The conversion between the 2 latitudes is as follows.

begin{align*} theta &= text{atan2}((1 - e^2)sinvarphi, cosvarphi)  varphi &= text{atan2}(sintheta, (1 - e^2)costheta) end{align*}

Right here e is eccentricity. The equations above work for any elliipsoid, however for earth particularly e² = 0.00669438.

The perform atan2(y, x) returns an angle in the identical quadrant as the purpose (x, y) whose tangent is y/x. [1]

As a fast sanity verify on the equations, word that when eccentricity e is zero, i.e. within the case of a circle, φ = θ. Additionally, if φ = 0 then θ = φ for all eccentricity values.

Subsequent we give a proof of the equations above.

Proof

We are able to parameterize an ellipse with semi-major axis a and semi-minor axis b by

(x(t), y(t)) = (a cos t, b sin t)

The slope at some extent (x(t), y(t)) is the ratio

frac{y^prime(t)}{x^prime(t)} = frac{b cos t}{-a sin t}

and so the slope of a line perpendicular to the tangent, i.e tan φ, is

tan varphi = frac{a sin t}{b cos t} = frac{a}{b} tan t

Now

tan theta = frac{b sin t}{a cos t} = frac{b}{a} tan t

and so

begin{align*} tan varphi &= frac{a}{b} tan t  &= frac{a}{b} left( frac{a}{b} tan theta right)  &= frac{a^2}{b^2} tan theta  &= frac{1}{1 - e^2} tan theta end{align*}

the place e² = 1 − b²/a² is the eccentricity of the ellipse. Due to this fact

(1 - e^2) tan varphi = tan theta

and the equations on the prime of the put up comply with.

Distinction

For the earth’s form, e² = 0.006694 per WGS84. For small eccentricities, the distinction between geographic and geocentric latitude is roughly symmetric round 45°.

However for bigger values of eccentricity the asymmetry turns into extra pronounced.

Associated posts

[1] There are a pair problems with programming language implementations of atan2. Some name the perform arctan2 and a few reverse the order of the arguments. Extra on that right here.

Related Articles

Latest Articles