Angular eccentricity

From wiki.gis.com
Jump to: navigation, search

In the study of ellipses and related geometry, various parameters in the distortion of a circle into an ellipse are identified and employed: Aspect ratio, flattening and eccentricity. All of these parameters are ultimately trigonometric functions of the ellipse's modular angle, or angular eccentricity.

Basic trigonometric functions

With the basic right triangle, the two sides adjoining the 90° angle (here, "b" and "c") are the triangle's "legs" and the third, longest, opposite side ("a") is the "hypotenuse".

Trig Fnctns.PNG

The doubled and squared half-angle functions (or "versed", meaning "turned", here, through 90°)—versine, vercosine, coversine and covercosine—have fallen into general obscurity, with function designation and abbreviation becoming ambiguous and even interchangeable (e.g., coversine is termed in some references as "vercosine" and "ver(C)" is also denoted as "vers(C)" and "versin(C)", while "cov(C)" is sometimes denoted as either "cvs(C)" or more commonly, "coversin(C)"). Also, \begin{matrix}{}^{}\sin^2\left(\frac{C}{2}\right)\end{matrix}\,\! is separately identified as "haversine" and \begin{matrix}{}^{}\sin^2\left(\frac{90^\circ-C}{2}\right)\end{matrix}\,\! as "hacoversine" (extending that terminology, \begin{matrix}{}^{}\tan^2\left(\frac{C}{2}\right)\end{matrix}\,\! can be regarded as "havertangent").

Linear Eccentricity

The parameters of an ellipse involve the same components and behave the same way as any right triangle, with one major exception: Physically speaking, there is no hypotenuse, only two "legs"—the semi-major, or transverse, and semi-minor, or conjugate, axes, or (as applied to a sphere or ellipsoid) the equatorial (or transverse) and polar (or conjugate) radii, a\,\! and b\,\!. Instead, an equivalent right triangle is created and defined, where a\,\! is the hypotenuse, b\,\! is the leg adjoining a\,\! at angle o\!\varepsilon\,\! and the complementary, imaginary "leg" is the half-focal separation, or linear eccentricity, E\,\!:  E=\sqrt{a^2-b^2}.\,\!

This "imaginary leg" equals the distance from the center of the ellipse to the focus:

El'iptic TriAngle.PNG

Notation

The generally accepted symbol for angular eccentricity is the Greek letter alpha, \alpha\,\![1].
However, not only is it more universally recognized as the denotation for navigational azimuth, but \alpha\,\! is also nearly identical, glyphically, to italic a, a\,\!, which is used to denote the major radius of an ellipse which, in turn, is used in the different calculations of \alpha\,\!, itself: E.g., \alpha=\arccos\left(\frac{b}{a}\right)\,\!.
Besides italic a, it is also glyphically a look-alike to the ancient Greek double character ("diphthong") omicron iota, o\iota\,\!, joined together ("ligaturized"), o\!\iota\,\!.
Over time, an alternative to the Latin form of o\iota\,\!, oi, came into use, the diphthong "ethel" (or "oethel"), oe, ligaturized as o\!e\,\!. In even more recent times, like its cousin ae, in most cases oe has simplified to just "e" (e.g., "oeconomy" → "economy"),[2] which is the Latin form of the Greek letter (var)epsilon, ε, the symbol for eccentricity.
Given this background, an alternative notation for angular eccentricity is used here, based on the graphenic derivation of e from o\!\iota\,\!: \alpha\approx{o\!\iota}\to{o\iota}\to{oe}\to{o\!e}\to{e}\to\varepsilon\to{o\!\varepsilon}.\,\!.
Where o\!e\,\! is ethel, o\!\varepsilon\,\! is Greek ethel, or "grethel".
The LaTex source code for this ligature is "o\!\varepsilon" (or, if——like "e" for ε——just "œ" is used, "o\!e").
There is also a more conceptual glyphical basis for o\!\varepsilon\,\!, as "o" can be regarded as a circle, with eccentricity, ε, pressing into it, making it elliptical: \mathrm{o}\;\!\varepsilon\to{\mathit{o}\varepsilon}\to{o\!\varepsilon}\to{\mathfrak{o}\!\varepsilon}\to{\mathfrak{o}\!\;\!\!\!\;\!\varepsilon}.\,\!.

Elliptic parameters

Like any angle, o\!\varepsilon\,\! can be found via the inverse of any trigonometric function it is the argument of:

\begin{matrix}{}_{\color{white}.}\\\;o\!\varepsilon\!\!\!&=&\!\!\!\arcsin\left(\frac{E}{a}\right)=\boldsymbol{\arccos\left(\frac{b}{a}\right)}=\arctan\left(\frac{E}{b}\right),\\&=&\!\!\!\!\boldsymbol{2\arctan\left(\!\sqrt{\frac{a-b}{a+b}}\;\right)}=2\arctan\left(\frac{E}{a+b}\right).\quad\\&&\!\!\!\!{}_{\mathrm{(As\;\mathit{a}\;and\;\mathit{b}\;are\;the\;known\;sides,\;the\;cosine\;and}}\qquad\\&&\!\!\!\!\!\!{}^{\mathrm{''\!havertangent''\;are\;the\;prime\;functions\;to\;inverse.)}}\end{matrix}\,\!

There are three primary parameters used in defining and constructing an elliptic figure: Circularity, eccentricity and flattening.

Circularity

The most concrete, tangible characteristic of an ellipse is its circularity, or diametric ratio, d. Also known as the aspect ratio, the circularity is the ratio of the conjugate to transverse radii, which equals the ellipse's angular eccentricity cosine:
d\;\;=\;\;\frac{b}{a}\;\;=\;\;\cos(o\!\varepsilon).\,\!
It is this defining measurement that is visually discernible. For example, if the circularity of an ellipse is .5, then the (central) vertical diameter is one-half that of the horizontal, if .1, then one-tenth, if .01, then one-hundredth, etc. The extremes and middle valued ellipses work out to the following:
OE Interval Ellps.PNG

Eccentricity

The eccentricity is actually a trio of factors: The primary, or first, eccentricity, ε (for basic text notation, "ε" is commonly expressed as "e"), is o\!\varepsilon\,\!'s sine, the second eccentricity, ε', is its tangent, and the third, ε" (also denoted in its squared form as m), is (in terms of function identity) ambiguous:

\begin{matrix}
\;\varepsilon^2&=&\frac{a^2-b^2}{a^2}&=&\frac{\sin^2(o\!\varepsilon)}{1}&=&\sin^2(o\!\varepsilon);\\\\
\varepsilon'^2&=&\frac{a^2-b^2}{b^2}&=&\frac{\sin^2(o\!\varepsilon)}{1-\sin^2(o\!\varepsilon)}&=&\tan^2(o\!\varepsilon);\\\\
\varepsilon''^2&=&\frac{a^2-b^2}{a^2+b^2}&=&\frac{\sin^2(o\!\varepsilon)}{2-\sin^2(o\!\varepsilon)}.\\
{}^{}\end{matrix}\,\!
Since they are mostly used in that form anyway, the eccentricities are usually found and kept in their squared form.
The primary eccentricity could be regarded as the complementary circularity, as it is the ratio of the linear eccentricity to the semi-major axis:  \begin{matrix}\varepsilon=\frac{\sqrt{a^2-b^2}}{a}=\frac{E}{a}.\end{matrix}\,\!

Flattening

The flattening, or ellipticity, in contrast, is self-explanatory, as it defines the degree of "squashing", from no flattening (a perfect circle) to complete flattening (a straight line). Just as the eccentricity is based on o\!\varepsilon\,\!'s sine, the flattening is based on its versine. Also like the eccentricity, there are three forms of flattening—the primary, or first, flattening, f, which is o\!\varepsilon\,\!'s versine, a second, f' (rarely used or even acknowledged), and a third, f'' (more commonly denoted as n), which is o\!\varepsilon\,\!'s "havertangent"[3]. The symmetry of the flattenings is similar in nature to the eccentricities, with the differences being that the radii aren't squared and the sines are doubled with half-angles:
\begin{matrix}
{}_{}\\\;f&=&\frac{a-b}{a}&=&\frac{2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}{1}&=&1-\cos(o\!\varepsilon)&=&\mbox{ver}(o\!\varepsilon);\\\\
f'&=&\frac{a-b}{b}&=&\frac{2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}{1-2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}&=&\sec(o\!\varepsilon)-1&=&\frac{\mbox{ver}(o\!\varepsilon)}{1-\mbox{ver}(o\!\varepsilon)};\\\\
f''&=&\frac{a-b}{a+b}&=&\frac{2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}{2-2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}&=&\tan^2\left(\frac{o\!\varepsilon}{2}\right)&=&\frac{\mbox{ver}(o\!\varepsilon)}{\mbox{vrc}(o\!\varepsilon)}.\\
{}^{}\end{matrix}\,\!
While the circularity would seem to be the ideal parameter to find an unknown axis (usually b), it is usually the inverse (primary) flattening that is provided:
\begin{matrix}{}_{}\\\mathrm{E.g.,\;\;}a=6378,\;\frac{1}{f}=300\!:\;b=a(1-\frac{1}{300})=a\cos(o\!\varepsilon)=6356.74\;\\{}^{}\end{matrix}.\,\!

Oblate vs. prolate

The basic object of elliptic geometry is the circle. If the two dimensional circle is expanded into a three dimensional solid, it becomes a sphere. Likewise, if one expands a two dimensional ellipse into a three dimensional solid, it becomes an ellipsoid. If the ellipse is rotated about its polar axis, it is known as an ellipsoid of revolution, specifically an oblate spheroid, or oblatum, where a > b—like an ellipse. If it is rotated about its equatorial axis, it is a prolate spheroid, or prolatum.

OblateSpheroid.PNG ← Oblate;   Prolate → ProlateSpheroid.png

In dealing with elliptic parameters, the main, oblate-prolate contrasting characteristic is the radius catagorization.

  • In all situations, the equatorial radius, a\,\!, is the transverse radius, while the polar, b\,\!, is the conjugate radius (putting it in proper perspective, however, the conjugate (b\,\! radius) is considered the base reference, while the transverse——meaning perpendicular or "crosswise"——(a\,\! radius) is the contrast).
  • The smaller axis/radius is always the "minor", while the larger axis/radius is always the "major" (thus, for an oblatum, the equator possesses the transverse-major axis/radius and the poles the conjugate-minor, while the equator of a prolatum defines the transverse-minor axis/radius and the poles the conjugate-major).

Due to their rotation, most of the planets (including Earth) and their satellites are (even if minimally) oblate spheroids. As such, planetodetic formulation utilizes the oblate format, which follows standard elliptic parameterization.

Prolate valuations

In all cases of eccentricity and flattening, the numerator is always the major minus the minor radii, while the denominators of the eccentricities and third flattening are also fixed as to the major-minor status, but the first and second flattenings of prolate valuations are fixed to the conjugate-transverse status [4](thus their o\!\varepsilon\,\! function assignments are reversed):

o\!\varepsilon=\arccos\left(\frac{a}{b}\right)=2\arctan\left(\!\sqrt{\frac{b-a}{b+a}}\;\right);\,\!
\begin{matrix}
\;\varepsilon^2&=&\frac{b^2-a^2}{b^2}&=&\frac{\sin^2(o\!\varepsilon)}{1}&=&\sin^2(o\!\varepsilon);\\\\
\varepsilon'^2&=&\frac{b^2-a^2}{a^2}&=&\frac{\sin^2(o\!\varepsilon)}{1-\sin^2(o\!\varepsilon)}&=&\tan^2(o\!\varepsilon);\\\\
\varepsilon''^2&=&\frac{b^2-a^2}{b^2+a^2}&=&\frac{\sin^2(o\!\varepsilon)}{2-\sin^2(o\!\varepsilon)}.\\
{}^{}\end{matrix}\,\!
\begin{matrix}
{}_{}\\\;f&=&\frac{b-a}{a}&=&\frac{2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}{1-2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}&=&\sec(o\!\varepsilon)-1&=&\frac{\mbox{ver}(o\!\varepsilon)}{1-\mbox{ver}(o\!\varepsilon)};\\\\
f'&=&\frac{b-a}{b}&=&\frac{2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}{1}&=&1-\cos(o\!\varepsilon)&=&\mbox{ver}(o\!\varepsilon);\\\\
f''&=&\frac{b-a}{b+a}&=&\frac{2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}{2-2\sin^2\left(\frac{o\!\varepsilon}{2}\right)}&=&\tan^2\left(\frac{o\!\varepsilon}{2}\right)&=&\frac{\mbox{ver}(o\!\varepsilon)}{\mbox{vrc}(o\!\varepsilon)}.\\
{}^{}\end{matrix}\,\!

Applications

For the most part, elliptic formularies ignore the angular eccentricity for the more familiar and notationally concise ε2, ε' 2, and f. However, these parameters don't provide the clear and obvious transformational relationships and structure. Consider the basic elliptic integrand at point P:

\begin{matrix}{}_{}\\
\sqrt{1-\sin^2(P)\varepsilon^2}\!&=&\!\!\sqrt{1-(\sin(P)\sin(o\!\varepsilon))^2},\qquad\qquad\qquad\qquad\qquad\qquad\\
&&\!\!\!\!\!\!{}_{=\;\sqrt{1-(1-\cos^2(P))\sin^2(o\!\varepsilon)},}\qquad\qquad\qquad\qquad\qquad\\\\
&=&\!\!\!\!\!\sqrt{\cos^2(o\!\varepsilon)+(\cos(P)\sin(o\!\varepsilon))^2},\qquad\qquad\qquad\qquad\\
&&{}_{\;=\;\cos(o\!\varepsilon)\sqrt{(\cos(P)\tan(o\!\varepsilon))^2+1},}\qquad\qquad\qquad\qquad\qquad\\\\
&=&\!\!\!\!\!\!\!\!\!\!\!\!\sqrt{\cos^2(o\!\varepsilon)+\frac{1}{2}(1+\cos(2P))\sin^2(o\!\varepsilon)},\qquad\qquad\\
&=&\!\!\!\!\!\!\sqrt{\cos^2(o\!\varepsilon)+\frac{1}{2}\sin^2(o\!\varepsilon)+\frac{1}{2}cos(2P)\sin^2(o\!\varepsilon)},\qquad\\
&&\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!{}_{=\sqrt{1-\frac{1}{2}\varepsilon^2+\frac{1}{2}cos(2P)\varepsilon^2},}\qquad\qquad\qquad\qquad\qquad\\\\
&=&\!\!\!\sqrt{\cos^4\left(\frac{o\!\varepsilon}{2}\right)+\sin^4\left(\frac{o\!\varepsilon}{2}\right)+2cos(2P)\sin^2\left(\frac{o\!\varepsilon}{2}\right)\cos^2\left(\frac{o\!\varepsilon}{2}\right)},\\\\
&=&\!\!\!\!\!\!\cos^2\left(\frac{o\!\varepsilon}{2}\right)\sqrt{1+2\cos(2P)\tan^2\left(\frac{o\!\varepsilon}{2}\right)+\tan^4\left(\frac{o\!\varepsilon}{2}\right)},\qquad\\\\
&=&\!\!\!\!\sqrt{\frac{1}{(1+f')^2}+\frac{1}{2}\cos(2P)\varepsilon^2+\frac{1}{4}f^2}\;=\;\frac{\sqrt{1+2\cos(2P)f'+f'^2}}{1+f'}.
\\{}^{}\end{matrix}\,\!

While one may consider such ability to convert as just gratuitously frivolous, there is at least one valid reason, as the Binomial series expansion (which planetodetic formularies frequently use) for {}^{\sqrt{1+2\cos(2P)\tan^2\left(\frac{o\!\varepsilon}{2}\right)+\tan^4\left(\frac{o\!\varepsilon}{2}\right)}}\,\! converges a lot quicker than the one for {}^{\sqrt{1-(\sin(P)\sin(o\!\varepsilon))^2}}\,\! which, in turn, converges quicker than {}_{\sqrt{(\cos(P)\tan(o\!\varepsilon))^2+1}}\,\!'s (which—in line with basic, series expansion theory—doesn't even converge when {}^{(\cos(P)\tan(o\!\varepsilon))^2}\,\! ≥ 1). Furthermore, as {}_{\cos^{\frac{1}{2}}(o\!\varepsilon)\approx\;\cos^2\left(\frac{o\!\varepsilon}{2}\right)\approx\;\cos^8\left(\frac{o\!\varepsilon}{4}\right)\approx\;\cos^{\!\!\frac{4^{{}^{q}}}{2}}\left(\frac{o\!\varepsilon}{2^q}\right)}\,\!, there are likely other, even more efficient, series expansions possible (if not even efficiently practical approximations to a general transcendental elliptic integral).
Another example is the equation for authalic surface area:

\begin{matrix}{}_{}\\{}^{}\mathrm{Oblate}\!\!\!&=&\!\!\!\!\!\!\!2\pi\left(a^2+\frac{b^2}{\varepsilon}\ln\left(\frac{\sqrt{1-\varepsilon^2}}{1-\varepsilon}\right)\right),\qquad\qquad\qquad\qquad\qquad\qquad\qquad\\\\
&=&\!\!\!2\pi\left(a^2+b^2\frac{\ln\left(\sec(o\!\varepsilon)(1+\sin(o\!\varepsilon)) \right)}{\sin(o\!\varepsilon)}\right)
=2\pi\left(a^2+b^2\frac{\operatorname{arcosh}(\sec(o\!\varepsilon))}{\sin(o\!\varepsilon)}\right);\\{}^{}\end{matrix}\,\!
\begin{matrix}{}_{}\\\mathrm{Prolate}\!\!\!&=&\!\!\!2\pi\left[ab\frac{\arcsin(\varepsilon)}{\varepsilon}+b^2\right]=2\pi\left[a^2\frac{o\!\varepsilon}{\sin(o\!\varepsilon)\cos(o\!\varepsilon)}+b^2\right],\\\\
&=&\!\!\!\!\!\!\!\!2\pi\left[a^2\frac{2o\!\varepsilon}{\sin(2o\!\varepsilon)}+b^2\right]=2\pi\left[\frac{a^2}{\operatorname{sin\!c}(2o\!\varepsilon)}+b^2\right].\qquad\\{}^{}\end{matrix}\,\!

While one certainly can use ε to define and express this type of equation, using o\!\varepsilon\,\! frequently provides a more illustrative—if not even its definitively mathematical—origin.

Notes

  1. Rapp, Richard H., Geometric Geodesy, Part I, , (April 1991), Dept. of Geodetic Science and Surveying, Ohio State Univ., Columbus, Ohio, sec.3.1, pp.12–16.
  2. (Wikipedia) American and British English spelling differences: Simplification of ae and oe
  3. König, R. and Weise, K. H. (1951): Mathematische Grundlagen der höheren Geodäsie und Kartographie, Band 1, Das Erdsphäroid und seine konformen Abbildungen, Springer-Verlag, Berlin/Göttingen/Heidelberg, Sec I.1
  4. MathWorld: Flattening.

See also