HomeSite MapMollweideHow Projections are CreatedEquidistant Cylindrical, Winkel I/IIMap Projections

Two Aspects for Two Arbitrary Azimuthal Projections

General Polar Azimuthal Projections

The mathematical development of the azimuthal orthographic projection is purely geometric.  Even though some azimuthal projections do not follow such a perspective process, all can be reduced to a general pattern, given the tangency point T:
Sphere to map
Outlined projection plane tangent to sphere at point T
theta = theta' due to the azimuthal property
rho = f (r)
x = rho cos theta
y = rho sin theta
For the North polar aspect, theta = lambda.

The Azimuthal Equidistant Projection

For the azimuthal equidistant, an important projection for navigational applications, distance rho from the center of the map is directly proportional to radial distance from the tangent point.  In the North polar aspect:

rho = (pi / 2 - phi)R

The austral aspect is just as easy, with

rho = (pi / 2 + phi)R and theta = -lambda

Lambert's Azimuthal Equal-area Projection

In the only projection both azimuthal and equal-area, created by Lambert and suitable for world maps, distance from the center of the map is progressively reduced in order to keep areal equivalence.  Formulas follow from basic integral calculus.
Area element on sphere, given colatitude Phi = pi / 2 - phi:
ds = 2piR sin Phi R dPhi = 2piR² sin Phi dPhi
Area element on globe
Corresponding element on map: ds = 2pirdr Area element on map
For a given Phi1, we want rho: Radial distance for Lambert map

A similar sign change applies if the south polar aspect is intended.

Azimuthal equidistant and equal-area polar maps
Combined azimuthal equidistant (top) and equal-area (bottom) polar maps
The combined map shows both projections virtually identical at latitudes above 70°N.  Beyond that, parallels get closer and closer together in Lambert's half, while remaining equally spaced in the azimuthal equidistant portion.  The resulting areal difference is clearly visible in Antarctica.

General Equatorial Aspect for Azimuthal Maps

Calculating other aspects for azimuthal maps is possible applying coordinate transformations and rotations in space.  However, the important equatorial aspect can be obtained in a more direct way, using two properties of triangles on a spherical surface.

Spherical triangle Given A, B, C angles on a spherical triangle's vertices, and alpha, beta, gamma the corresponding angles between edges connecting triangle vertices and center O of sphere,
Law of sines: sin A / alpha = sin B / beta = sin C / gamma
Law of cosines: cos gamma = cos alpha cos beta + sin alpha sin beta sin  C

Spherical triangle on Earth On the equatorial aspect, tangent point T lies on the intersection of Equator and an arbitrary central meridian.  The projected point P marks a shaded triangle, whose vertices define central angle alpha, latitude phi and longitude lambda
Formulas for equatorial aspect

The Equatorial Azimuthal Equidistant

Just substituting results for cos theta and sin theta,

rho = r = alphaR
x = alphaR cos phi sin lambda / sin alpha
y = alphaR sin phi / sin alpha

If lambda = phi = 0, sin alpha = 0, x = y = 0.

The Equatorial Azimuthal Equal-Area

In equations for Lambert's equal-area azimuthal projection, substitute alpha for Phi1:

x = rho cos theta = rho cos phi / sin lambda sin alpha
y = rho sin theta = rho sin phi / sin alpha

sin (a + b) = sin a cos b + cos a sin b, therefore sin 2a = 2 sin a cos a
cos (a + b) = cos a cos b - sin a sin b; cos 2a = cos2a - sin2a = 2 cos2a - 1, therefore (cos 2a + 1) / 2 = cos2a

Formulas for equatorial aspect of Lambert's projection

Azimuthal equidistant and equal-area polar maps
Combined azimuthal equidistant (top) and equal-area (bottom) equatorial maps
Again, both projections are very similar near the tangent point: the northern and southern portions of Africa join almost seamlessly.

HomeSite MapMollweideHow Projections are CreatedEquidistant Cylindrical, Winkel I/II  www.progonos.com/furuti    January 14, 2004
Copyright © 2001 Carlos A. Furuti