Spherical Coordinates

Spherical coordinates define a point in 3D space (may be generalized to higher dimensions) by using a radial distance from the origin and angles from reference axes. It is a higher-dimensional equivalent to 2D polar coordinates.

Conventions

There are several conventions for representing spherical coordinates. Each is discussed here.

Latitude

Latitude angles are measured from an equatorial plane passing through the origin (canonically the xy-plane). This type of measurement is used typically in geospatial contexts.

Given spherical coordinates with radius, longitude, and latitude $⟨ r; θ; ϕ ⟩$ , the corresponding point in Cartesian coordinates is given by:

$\begin{aligned} x & = r \cos ϕ \cos θ \\ y & = r \cos ϕ \sin θ \\ z & = r \sin ϕ \end{aligned}$

This convention will be used for other areas of this page unless otherwise specified.

Colatitude

Colatitude angles are measured from the z-axis. This convention is seldom used outside the realm of mathematics. Even so, it is used in conjunction with other conventions.

Given spherical coordinates with radius, longitude, and colatitude $⟨ r; θ; φ ⟩$ , the corresponding point in Cartesian coordinates is given by:

$\begin{aligned} x & = r \sin φ \cos θ \\ y & = r \sin φ \sin θ \\ z & = r \cos φ \end{aligned}$

Angle Between Vectors

Theorem. Given vectors ${\vec{v}}_{1} = ⟨ r_{1}; θ_{1}; ϕ_{1} ⟩$ and ${\vec{v}}_{2} = ⟨ r_{2}; θ_{2}; ϕ_{2} ⟩$ , the angle $ψ$ between them is:

ψ = \arccos (\cos ϕ_{1} \cos ϕ_{2} \cos (Δ θ) + \sin ϕ_{1} \sin ϕ_{2})

Where $Δ θ = | θ_{2} - θ_{1} |$ .

Proof. Note that $‖ {\vec{v}}_{1} ‖ = r_{1}$ and $‖ {\vec{v}}_{2} ‖ = r_{2}$ . The dot product is given by:

\begin{aligned} {\vec{v}}_{1} \cdot {\vec{v}}_{2} & = r_{1} r_{2} \cos ψ \\ = x_{1} x_{2} + y_{1} y_{2} + z_{1} z_{2} \end{aligned}

Where $⟨ x_{1}, y_{1}, z_{1} ⟩$ and $⟨ x_{2}, y_{2}, z_{2} ⟩$ are the Cartesian coordinates of ${\vec{v}}_{1}$ and ${\vec{v}}_{2}$ , respectively. Equating both definitions and solving for $ψ$ gives:

\begin{aligned} \cos ψ & = \frac{x_{1} x_{2} + y_{1} y_{2} + z_{1} z_{2}}{r_{1} r_{2}} \\ ψ & = \arccos (\frac{x_{1} x_{2} + y_{1} y_{2} + z_{1} z_{2}}{r_{1} r_{2}}) \end{aligned}

Substituting the above formulae for these coordinates into the alternate dot product definition gives:

\begin{aligned} ψ & = \arccos (\frac{r_{1} r_{2} \cos ϕ_{1} \cos ϕ_{2} \cos θ_{1} \cos θ_{2} + r_{1} r_{2} \cos ϕ_{1} \cos ϕ_{2} \sin θ_{1} \sin θ_{2} + r_{1} r_{2} \sin ϕ_{1} \sin ϕ_{2}}{r_{1} r_{2}}) \\ = \arccos (\cos ϕ_{1} \cos ϕ_{2} \cos θ_{1} \cos θ_{2} + \cos ϕ_{1} \cos ϕ_{2} \sin θ_{1} \sin θ_{2} + \sin ϕ_{1} \sin ϕ_{2}) \\ = \arccos (\cos ϕ_{1} \cos ϕ_{2} (\cos θ_{1} \cos θ_{2} + \sin θ_{1} \sin θ_{2}) + \sin ϕ_{1} \sin ϕ_{2}) \\ = \arccos (\cos ϕ_{1} \cos ϕ_{2} \cos (θ_{1} - θ_{2}) + \sin ϕ_{1} \sin ϕ_{2}) \end{aligned}

Note the use of the subtractive trigonometric identity for cosine. Since cosine is an even function (i.e. $\cos (θ) = \cos (- θ)$ , the order of subtraction does not matter; only the difference between $θ_{1}$ and $θ_{2}$ matters.

quod erat demonstrandum

Computational Formula

In situations where floating point rounding errors pose a challenge, the following alternate formula may be used (from wikipedia:Great-circle distance):

ψ = \arctan \frac{\sqrt{{(\cos ϕ_{2} \sin (Δ θ))}^{2} + {(\cos ϕ_{1} \sin ϕ_{2} - \sin ϕ_{1} \cos ϕ_{2} \cos (Δ θ))}^{2}}}{\sin ϕ_{1} \sin ϕ_{2} + \cos ϕ_{1} \cos ϕ_{2} \cos (Δ θ)}

Below is a Python implementation (using the NumPy library) of this formula.

import numpy as np

def angle_between(lng1, lat1, lng2, lat2):
    s1 = np.sin(lat1)
    c1 = np.cos(lat1)
    s2 = np.sin(lat2)
    c2 = np.cos(lat2)
    dlng = lng2 - lng1
    c2cdlng = c2 * np.cos(dlng)

    return np.arctan2(
        np.sqrt((c2*np.sin(dlng))**2 + (c1*s2 - s1*c2cdlng)**2),
        s1*s2 + c1*c2cdlng
    )

Spherical Coordinates

Contents

Conventions

Latitude

Colatitude

Angle Between Vectors

Computational Formula

Navigation menu

Spherical Coordinates

Conventions

Latitude

Colatitude

Angle Between Vectors

Computational Formula

Navigation menu

Search