Space-oblique Mercator projection

From Wikipedia, the free encyclopedia
Jump to: navigation, search
Space-oblique Mercator projection.

Space-oblique Mercator projection is a map projection.

History[edit]

The space-oblique Mercator projection (SOM) was developed by John P. Snyder, Alden Partridge Colvocoresses and John L. Junkins in 1976. Snyder had an interest in maps, originating back to his childhood and he regularly attended cartography conferences while on vacation. When the United States Geological Survey (USGS) needed to develop a system for reducing the amount of distortion caused when satellite pictures of the ellipsoidal Earth were printed on a flat page, they appealed for help at one such conference. Snyder worked on the problem armed with his newly purchased pocket calculator and devised the mathematical formulas needed to solve the problem. He submitted these to the USGS at no charge, starting off a new career at USGS. His formulas were used to produce maps from Landsat 4 images launched in the summer of 1978.

Projection description[edit]

The space-oblique Mercator projection provides continual conformal mapping of the swath[disambiguation needed] sensed by a satellite. Scale is true along the ground track, varying 0.01 percent within the normal sensing range of the satellite. Conformality is correct within a few parts per million for the sensing range. Distortion is essentially constant along lines of constant distance parallel to the ground track. SOM is the only projection presented that takes the rotation of Earth into account.

Equations[edit]

The forward equations for the Space Oblique Mercator projection for the sphere are as follows:


\begin{align}
\frac{x}{R} &= \int_{0}^{\lambda'} \frac{H-S^2}{\left(1+S^2\right)^{1/2}}d\lambda' - \frac{S}{\left(1+S^2\right)^{1/2}}\ln\tan\left(\frac{\pi}{4}+\frac{\phi'}{2}\right) \\
\frac{y}{R} &= \left(H+1\right) \int_{0}^{\lambda'} \frac{S}{\left(1+S^2\right)^{1/2}}d\lambda' + \frac{1}{\left(1+S^2\right)^{1/2}}\ln\tan\left(\frac{\pi}{4}+\frac{\phi'}{2}\right) \\
S &= \left(P_{2}/P_{1}\right) \sin i \cos \lambda' \\
H &= 1 - \left(P_{2}/P_{1}\right) \cos i \\
\tan\lambda' &= \cos i \tan \lambda_{t} + \sin i \tan \phi / \cos \lambda_{t} \\
\sin\phi' &= \cos i \sin \phi - \sin i \cos \phi \sin \lambda_{t} \\
\lambda_{t} &= \lambda + \left(P_{2}/P_{1}\right) \lambda'. \\
\phi &= \text{geodetic (or geographic) latitude.} \\
\lambda &= \text{geodetic (or geographic) longitude.} \\
P_{2} &= \text{time required for revolution of satellite.} \\
P_{1} &= \text{length of Earth rotation.} \\
i &= \text{angle of inclination.} \\
R &= \text{radius of Earth.} \\
x,y &= \text{rectangular map coordinates.}
\end{align}

References[edit]