Ellipsoid

An ellipsoid is a type of quadric surface that is a higher dimensional analogue of an ellipse. It is similar to a sphere, except that the radii along one of the axis are longer than the radii along the other axis. It is important to understand ellipsoids in geography because the shape of the Earth is an ellipsoid, with the distance around the Equator longer than the distance around the Prime Meridian. Therefore, all map projections must take into account that the Earth is an ellipsoid to help minimize distortion.

The equation of a standard axis-aligned ellipsoid body in an X-Y-Z Cartesian coordinate system is
 * $${X^2 \over a_x^2}+{Y^2 \over a_y^2}+{Z^2 \over b^2}=1\,\!$$

where ax and ay are the transverse, equatorial radii (along the X- and Y- axes) and b is the conjugate, polar radius (along the Z-axis), all of which are fixed positive real numbers determining the shape of the ellipsoid (traditionally, ax and ay are denoted as a and b, respectively——thereby defining the equator as an ellipse——and b as c; however, in most cases——such as Earth——the equator is considered spherical, with the equatorial radius being defined simply as a and the polar as b——like an ellipse——thus a,b,c notation can be needlessly confusing).

More generally, a not-necessarily-axis-aligned ellipsoid is defined by the equation
 * $$\mathbf{x}^\mathrm{T} A \mathbf{x} = 1$$

where A is a symmetric positive definite matrix and x is a vector. In that case, the eigenvectors of A define the principal directions of the ellipsoid and the inverse of the square root of the eigenvalues are the corresponding equatorial radii.

If all three radii are equal, the solid body is a sphere; if two radii are equal, the ellipsoid is a spheroid:


 * $$a_x=a_y=b:\,\!$$ Sphere;
 * $$a_x=a_y>b:\,\!$$ Oblate spheroid, or oblatum (disk-shaped);
 * $$a_x=a_ya_y>b:\,\!$$ Scalene ellipsoid ("three unequal sides").

The points (ax,0,0), (0,ay,0) and (0,0,b) lie on the surface and the line segments from the origin to these points are called the semi-principal axes. These correspond to the semi-major axis and semi-minor axis of the appropriate ellipses.

Scalene ellipsoids are frequently called "triaxial ellipsoids", the implication being that all three axes need to be specified to define the shape.

Parameterization
A biaxial ellipse can be parameterized with Cartesian coordinates as


 * $$\begin{align}

x&=a\,\cos(\beta);\\ y&=b\,\sin(\beta);\end{align}\,\!$$

Where $$\beta\,\!$$ is a point's reduced, or parametric, latitude and $${\color{white}+}\!\!\!\lambda{\color{white}'}\,\!$$ is its planetographic longitude, an ellipsoid can be parameterized either triaxially or biaxially, by making $$a\,\!$$ a variable with respect to $$\lambda\,\!$$, $$a(\lambda)\,\!$$ :


 * $$\begin{align}

X&=a_x\cos(\beta)\cos(\lambda);\,{\color{white}|}\\ Y&=a_y\cos(\beta)\sin(\lambda);\\ Z&=b\sin(\beta);\end{align}\,\!$$ or
 * $$\begin{align}

a(\lambda)&=\sqrt{(a_x\cos(\lambda))^2+(a_y\sin(\lambda))^2}^{\color{white}|};\\ x(\lambda)&=\sqrt{X^2+Y^2}^{\color{white}|}=a(\lambda)\cos(\beta);\\ y&=\qquad Z\qquad=b\sin(\beta);\end{align}\,\!$$
 * $$\begin{matrix}-\frac{\pi}{2}\leq\beta\leq+\frac{\pi}{2};

\quad-\pi\leq\lambda\leq+\pi;\!{\color{white}\big|}\end{matrix}\,\!$$
 * (Note that this parameterization is not 1-1 at the poles, where $$\scriptstyle{{\color{white}|}\beta=\pm{\frac{\pi}{2}}}{\color{white}|}\,\!$$)

Volume
The volume of an ellipsoid is given by the formula
 * $$\frac{4}{3}\pi a_xa_yb;\,\!$$

Note that this equation reduces to that of the volume of a sphere when all three elliptic radii are equal, and to that of an oblate or prolate spheroid when two of them are equal.

Surface area
The surface area of an ellipsoid is given by:
 * $$2\pi\left(a_xa_y\sin(o\!\varepsilon_x)E(o\!\varepsilon_x,m)+a_yb\cot(o\!\varepsilon_x)F(o\!\varepsilon_x,m)+b^2\right);\,\!$$

where
 * $$o\!\varepsilon=\arccos\left(\frac{b}{a}\right)\;\textrm{(oblate\;or\;scalene)\;\;or\;\;}\arccos\left(\frac{a}{b}\right)\;(\textrm{prolate}),\,\!$$

is the modular angle, or angular eccentricity; $$m=\Big(\tfrac{\sin(o\!\varepsilon_y)}{\sin(o\!\varepsilon_x)}\Big)^2\,\!$$; and $$E(o\!\varepsilon_x,m)\,\!$$, $$F(o\!\varepsilon_x,m)\,\!$$ are the incomplete elliptic integrals of the second and first kind, respectively:
 * $$\begin{align}

E(o\!\varepsilon_x,m)&=E\!\left(o\!\varepsilon_x,\Big(\tfrac{\sin(o\!\varepsilon_y)}{\sin(o\!\varepsilon_x)}\Big)^2\right)=\int_0^{o\!\varepsilon_x}\sqrt{1-\Big(\sin(o\!\varepsilon)\tfrac{\sin(o\!\varepsilon_y)}{\sin(o\!\varepsilon_x)}\Big)^2}^{\color{white}|}\mbox{d}o\!\varepsilon;\\ F(o\!\varepsilon_x,m)&=F\!\left(o\!\varepsilon_x,\Big(\tfrac{\sin(o\!\varepsilon_y)}{\sin(o\!\varepsilon_x)}\Big)^2\right)=\int_0^{o\!\varepsilon_x}\frac{\mbox{d}o\!\varepsilon}{\sqrt{1-\Big(\sin(o\!\varepsilon)\tfrac{\sin(o\!\varepsilon_y)}{\sin(o\!\varepsilon_x)}\Big)^2}},\\ &{\color{white}|}\qquad\qquad\qquad\qquad\quad\,\!= \int_0^{o\!\varepsilon_x}\frac{\sin(o\!\varepsilon_x)\mbox{d}o\!\varepsilon}{ \sqrt{\sin^2(o\!\varepsilon_x)-\Big(\sin(o\!\varepsilon)\sin(o\!\varepsilon_y)\Big)^2}};\end{align}\,\!$$ In the case of a spheroid, the above simplifies to closed-form expressions:
 * Oblate: $$2\pi\!\left(a_x^2+b^2\frac{\operatorname{arctanh}(\sin(o\!\varepsilon))}{\sin(o\!\varepsilon)}\right);\,\!$$  Prolate: $$2\pi\!\left(a_x^2+b^2\frac{o\!\varepsilon}{\tan(o\!\varepsilon)}\right);\,\!$$

An approximate formula for any ellipsoid is:
 * $$\approx 4\pi\!\left(\frac{ a_x^p a_y^p+a_x^p b^p+a_y^p b^p }{3}\right)^{1/p};\,\!$$

where p ≈ 1.6075 yields a relative error of at most 1.061% (Knud Thomsen's formula); a value of p = 8/5 = 1.6 is optimal for nearly spherical ellipsoids, with a relative error of at most 1.178% (David W. Cantrell's formula).

In the "flat" limit of $$b \ll a_x, a_y\,\!$$, the area is approximately $$2\pi a_xa_y.\,\!$$

Mass properties
The mass of an ellipsoid of uniform density is:
 * $$m = \rho V = \rho \frac{4}{3} \pi a_xa_yb;\,\!$$

where $$\rho\,\!$$ is the density.

The mass moments of inertia of an ellipsoid of uniform density are:
 * $$I_{\mathrm{xx}} = m {a_y^2+b^2 \over 5};\,\!$$
 * $$I_{\mathrm{yy}} = m {a_x^2+b^2 \over 5};\,\!$$
 * $$I_{\mathrm{zz}} = m {a_x^2+a_y^2 \over 5};\,\!$$

where $$I_{\mathrm{xx}}\,\!$$, $$I_{\mathrm{yy}}\,\!$$, and $$I_{\mathrm{zz}}\,\!$$ are the moments of inertia about the X, Y, and Z axes, respectively. Products of inertia are zero.

It can easily be shown that if ax=ay=b, then the moments of inertia reduce to those for a uniform-density sphere.

Conversely, if the mass and principle inertias of an arbitrary rigid body are known, an equivalent ellipsoid of uniform density can be constructed, with the following characteristics:
 * $$a_x = \sqrt{{5 \over 2} {I_{\mathrm{yy}}+I_{\mathrm{zz}}-I_{\mathrm{xx}} \over m}}$$
 * $$a_y = \sqrt{{5 \over 2} {I_{\mathrm{zz}}+I_{\mathrm{xx}}-I_{\mathrm{yy}} \over m}}$$
 * $$b = \sqrt{{5 \over 2} {I_{\mathrm{xx}}+I_{\mathrm{yy}}-I_{\mathrm{zz}} \over m}}$$
 * $$\rho = \frac{3}{4} {m \over \pi a_xa_yb}.\,\!$$

Rotational equilibrium
Scalene ellipsoids and cuboids rotate stably along their major or minor axes, but not along their median axis. This can be seen experimentally by throwing an eraser with some spin. In addition, moment of inertia considerations mean that rotation along the major axis is more easily perturbed than rotation along the minor axis. One practical effect of this is that scalene astronomical bodies such as generally rotate along their minor axes (as does the Earth, which is merely oblate); in addition, because of tidal locking, scalene moons in synchronous orbit such as those of Saturn orbit with their major axis aligned radially to their planet.

A relaxed ellipsoid, that is, one in hydrostatic equilibrium, has an oblateness ax − b directly proportional to its mean density and mean radius. Ellipsoids with a differentiated interior—that is, a denser core than mantle—have a lower oblateness than a homogeneous body. Over all, the ratio (ay–b)/(ax−b) is approximately 0.25, though this drops for rapidly rotating bodies.

Linear transformations
If one applies an invertible linear transformation to a sphere, one obtains an ellipsoid; it can be brought into the above standard form by a suitable rotation, a consequence of the spectral theorem. If the linear transformation is represented by a symmetric 3-by-3 matrix, then the eigenvectors of the matrix are orthogonal (due to the spectral theorem) and represent the directions of the axes of the ellipsoid: the lengths of the semiaxes are given by the eigenvalues.

The intersection of an ellipsoid with a plane is either empty, a single point, or an ellipse (including a circle).

One can also define ellipsoids in higher dimensions, as the images of spheres under invertible linear transformations. The spectral theorem can again be used to obtain a standard equation akin to the one given above.

Egg shape
The shape of a chicken egg is approximately that of half each a prolate and roughly spherical (potentially even minorly oblate) ellipsoid joined at the equator, sharing a principal axis of rotational symmetry. Although the term egg-shaped usually implies a lack of reflection symmetry across the equatorial plane, it may also refer to true prolate ellipsoids. It can also be used to describe the 2D figure that, revolved around its major axis, produces the 3D surface. See also oval.

Applications in Cartography
Although on small scales earth's true shape is practically irrelevant, it is essential to take the earth's shape into account when measuring distances between two points, as the distance across the surface of a sphere could greatly differ from that of an ellipsoid. Also, when dealing with coordinate systems, it is critical to be attentive to the ellipsoid used. The same point on earth can have drastically different locations according to different ellipsoids.