Ellipsoid

An ellipsoid is a surface that may be obtained from a sphere by deforming it by means of directional scalings, or more generally, of an affine transformation.

Examples of ellipsoids with equation
sphere (top, a=b=c=4),
spheroid (bottom left, a=b=5, c=3),
tri-axial ellipsoid (bottom right, a=4.5, b=6, c=3)

An ellipsoid is a quadric surface;  that is, a surface that may be defined as the zero set of a polynomial of degree two in three variables. Among quadric surfaces, an ellipsoid is characterized by either of the two following properties. Every planar cross section is either an ellipse, or is empty, or is reduced to a single point (this explains the name, meaning "ellipse like"). It is bounded, which means that it may be enclosed in a sufficiently large sphere.

An ellipsoid has three pairwise perpendicular axes of symmetry which intersect at a center of symmetry, called the center of the ellipsoid. The line segments that are delimited on the axes of symmetry by the ellipsoid are called the principal axes, or simply axes of the ellipsoid. If the three axes have different lengths, the ellipsoid is said to be tri-axial or rarely scalene, and the axes are uniquely defined.

If two of the axes have the same length, then the ellipsoid is an ellipsoid of revolution, also called a spheroid. In this case, the ellipsoid is invariant under a rotation around the third axis, and there are thus infinitely many ways of choosing the two perpendicular axes of the same length. If the third axis is shorter, the ellipsoid is an oblate spheroid; if it is longer, it is a prolate spheroid. If the three axes have the same length, the ellipsoid is a sphere.

Standard equation

Using a Cartesian coordinate system in which the origin is the center of the ellipsoid and the coordinate axes are axes of the ellipsoid, the implicit equation of the ellipsoid has the standard form

where a, b, c are positive real numbers.

The points (a, 0, 0), (0, b, 0) and (0, 0, c) lie on the surface. The line segments from the origin to these points are called the principal semi-axes of the ellipsoid, because a, b, c are half the length of the principal axes. They correspond to the semi-major axis and semi-minor axis of an ellipse.

If one has an oblate spheroid; if one has a prolate spheroid; if one has a sphere.

Parameterization

The ellipsoid may be parameterized in several ways, which are simpler to express when the ellipsoid axes coincide with coordinate axes. A common choice is

where

These parameters may be interpreted as spherical coordinates, where is the polar angle and is the azimuth angle of the point (x, y, z) of the ellipsoid.[1]

Measuring from the center rather than a pole,

where

is the reduced latitude, parametric latitude, or eccentric anomaly and is azimuth or longitude.

Measuring angles directly to the surface of the ellipsoid, not to the circumscribed sphere,

where

would be geocentric latitude on the Earth, and is azimuth or longitude. These are true spherical coordinates with the origin at the center of the ellipsoid.

For geodesy, geodetic latitude, the angle between the vertical and the equatorial plane, is most commonly used. Geodetic latitude is not defined for a general ellipsoid because it depends upon longitude.

Volume and surface area

Volume

The volume bounded by the ellipsoid is

Alternatively expressed, where A, B and C are the lengths of the principal axes (A=2a, B=2b and C=2c):

.

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.

The volume of an ellipsoid is the volume of a circumscribed elliptic cylinder, and the volume of the circumscribed box.

The volumes of the inscribed and circumscribed boxes are respectively:

Surface area

The surface area of a general (tri-axial) ellipsoid is[2][3]

where

and where F(φ,k) and E(φ,k) are incomplete elliptic integrals of the first and second kind respectively.

The surface area of an ellipsoid of revolution (or spheroid) may be expressed in terms of elementary functions:

which, as follows from basic trigonometric identities, are equivalent expressions (i.e. the formula for can be used to calculate the surface area of a prolate ellipsoid and vice versa). In both cases e may again be identified as the eccentricity of the ellipse formed by the cross section through the symmetry axis. (See ellipse). Derivations of these results may be found in standard sources, for example Mathworld.[4]

Approximate formula

Here p ≈ 1.6075 yields a relative error of at most 1.061%;[5] a value of p = 8/5 = 1.6 is optimal for nearly spherical ellipsoids, with a relative error of at most 1.178%.

In the "flat" limit of c much smaller than a, b, the area is approximately 2πab, equivalent to p ≈ 1.5850.

Plane sections

Properties

Plane section of an ellipsoid

The intersection of a plane and a sphere is a circle (or is reduced to a single point, or is empty). Any ellipsoid is the image of the unit sphere under some affine transformation, and any plane is the image of some other plane under the same transformation. So, because affine transformations map circles to ellipses, the intersection of a plane with an ellipsoid is an ellipse or a single point, or is empty.[6] Obviously, spheroids contain circles. This is also true, but less obvious, for triaxial ellipsoids (see Circular section).

Determining the ellipse of a plane section

Plane section of an ellipsoid (See example)

Given: Ellipsoid and the plane with equation which have an ellipse in common.
Wanted: Three vectors (center) and (conjugate vectors), such that the ellipse can be represented by the parametric equation

(See ellipse).
Plane section of the unit sphere (See example)

Solution: The scaling transforms the ellipsoid onto the unit sphere and the given plane onto the plane with equation . Let be the Hesse normal form of the new plane and its unit normal vector. Hence
is the center of the intersection circle and its radius (See diagram).
In case of let be (The plane is horizontal !)
In case of let be
In any case the vectors are orthogonal, parallel to the intersection plane and have length (radius of the circle). Hence the intersection circle can be described by the parametric equation
The reverse scaling (See above) transforms the unit sphere back to the ellipsoid and the vectors are mapped onto vectors , which were wanted for the parametric representation of the intersection ellipse.
How to find the vertices and semi-axes of the ellipse is described in ellipse.

Example: The diagrams show an ellipsoid with the semi-axes which is cut by the plane

Pins-and-string construction

Pins-and-string construction of an ellipse
length of the string (red)
Pins-and-string construction of an ellipsoid, blue: focal conics
Determination of the semi axis of the ellipsoid

The pins-and-string construction of an ellipsoid is a transfer of the idea constructing an ellipse using two pins and a string (see diagram).

A pins-and-string construction of an ellipsoid of revolution is given by the pins-and-string construction of the rotated ellipse.

The construction of points of a 3-axial ellipsoid is more complicated. First ideas are due to the Scottish physicist J. C. Maxwell (1868) [7]. Main investigations and the extension to quadrics was done by the German mathematician O. Staude in 1882, 1886 and 1898 [8] [9] [10]. The description of the pins-and-string construction of ellipsoids and hyperboloids is containd in the book Geometry and the imagination written by D. Hilbert & S. Vossen[11], too.

Steps of the construction
(1) Choose an ellipse and a hyperbola, which are a pair of focal conics:
Ellipse: and
Hyperbola:
with the vertices and foci of the ellipse
and a string (in diagram red) of length .
(2) Pin one end of the string to vertex and the other to focus . The string is kept tight at a point with positive y- and z-coordinates, such that the string runs from to behind the upper part of the hyperbola (see diagram) and is free to slide on the hyperbola. The part of the string from to runs and slides in front of the ellipse. The string runs through that point of the hyperbola, for which the distance over any hyperbola point is minimal. The analog statement on the second part of the string and the ellipse has to be true, too.
(3) Then: is a point of the ellipsoid with equation
and
(4) The remaining points of the ellipsoid can be constructed by suitable changes of the string at the focal conics.

Equations for the semi axes of the generated ellipsoid can be derived by special choices for point ::

The lower part of the diagram shows: are the foci of the ellipse in the x-y-plane, too. Hence, it is confocal to the given ellipse and the length of the string is . Solving for yields: . Further more: .
From the upper diagram one gets: are the foci of the ellipse (of the ellipsoid) in the x-z-plane and the equation .

Converse

If conversely a 3-axial ellipsoid is given by its equation then from the equations in step (3) one can derive the parameters for a pins-and-string construction.

Limit case: ellipsoid of revolution
In case of one gets , which means: the focal ellipse degenerates to a line segment and the focal hyperbola collapses to two infinite line segments on the x-axis. The ellipsoid is rotational symmetric with the x-axis as axis of rotation and .

In general position

As quadric

More generally, an arbitrarily oriented ellipsoid, centered at v, is defined by the solutions x to the equation

where A is a positive definite matrix and x, v are vectors.

The eigenvectors of A define the principal axes of the ellipsoid and the eigenvalues of A are the reciprocals of the squares of the semi-axes: , and .[12] An invertible linear transformation applied to a sphere produces an ellipsoid, which can be brought into the above standard form by a suitable rotation, a consequence of the polar decomposition (also, see 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 semi-axes are computed from the eigenvalues. The singular value decomposition and polar decomposition are matrix decompositions closely related to these geometric observations.

Parametric representation

ellipsoid as an affine image of the unit sphere

The key to a parametric representation of an ellipsoid in general position is the alternative definition:

  • An ellipsoid is an affine image of the unit sphere.

An affine transformation can be represented by a translation with a vector and a regular 3×3-matrix :

,

where are the column vectors of matrix .

A parametric representation of an ellipsoid in general position can be obtained by the parametric representation of a unit sphere (see above) and an affine transformation:

  • .

If the vectors form an orthogonal system, the points with vectors are the vertices of the ellipsoid and are the semi principal axes.

A surface normal vector at point is

For any ellipsoid there exists an implicit representation . If for simplicity the center of the ellipsoid is the origin, i.e. , the following equation describes the ellipsoid above:[13]

Applications

The ellipsoidal shape finds many practical applications:

Geodesy
  • Earth ellipsoid, a mathematical figure approximating the shape of the Earth.
  • Reference ellipsoid, a mathematical figure approximating the shape of planetary bodies in general.
Mechanics
  • Poinsot's ellipsoid, a geometrical method for visualizing the torque-free motion of a rotating rigid body.
  • Lamé's stress ellipsoid, an alternative to Mohr's circle for the graphical representation of the stress state at a point.
  • Manipulability ellipsoid, used to describe a robot's freedom of motion.
Crystallography
  • Index ellipsoid, a diagram of an ellipsoid that depicts the orientation and relative magnitude of refractive indices in a crystal.
  • Thermal ellipsoid, ellipsoids used in crystallography to indicate the magnitudes and directions of the thermal vibration of atoms in crystal structures.
Lighting
  • Ellipsoidal reflector floodlight
  • Ellipsoidal reflector spotlight
Medicine
  • Measurements obtained from MRI imaging of the prostate can be used to determine the volume of the gland using the approximation L × W × H × 0.52 (where 0.52 is an approximation for π/6) [14]

Dynamical properties

The mass of an ellipsoid of uniform density ρ is:

The moments of inertia of an ellipsoid of uniform density are:

For these moments of inertia reduce to those for a sphere of uniform density.

Artist's conception of Haumea, a Jacobi-ellipsoid dwarf planet, with its two moons

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.[15]

One practical effect of this is that scalene astronomical bodies such as Haumea generally rotate along their minor axes (as does Earth, which is merely oblate); in addition, because of tidal locking, moons in synchronous orbit such as Mimas orbit with their major axis aligned radially to their planet.

A spinning body of homogeneous self-gravitating fluid will assume the form of either a Maclaurin spheroid (oblate spheroid) or Jacobi ellipsoid (scalene ellipsoid) when in hydrostatic equilibrium, and for moderate rates of rotation. At faster rotations, non-ellipsoidal piriform or oviform shapes can be expected, but these are not stable.

Fluid dynamics

The ellipsoid is the most general shape for which it has been possible to calculate the creeping flow of fluid around the solid shape. The calculations include the force required to translate through a fluid and to rotate within it. Applications include determining the size and shape of large molecules, the sinking rate of small particles, and the swimming abilities of microorganisms.[16]

In probability and statistics

The elliptical distributions, which generalize the multivariate normal distribution and are used in finance, can be defined in terms of their density functions. When they exist, the density functions f have the structure:

where is a scale factor, is an -dimensional random row vector with median vector (which is also the mean vector if the latter exists), is a positive definite matrix which is proportional to the covariance matrix if the latter exists, and is a function mapping from the non-negative reals to the non-negative reals giving a finite area under the curve.[17] The multivariate normal distribution is the special case in which for quadratic form .

Thus the density function is a scalar-to-scalar transformation of a quadric expression. Moreover, the equation for any iso-density surface states that the quadric expression equals some constant specific to that value of the density, and the iso-density surface is an ellipsoid.

In higher dimensions

The volume of a higher-dimensional ellipsoid (a hyperellipsoid) can be calculated using the dimensional constant given for the volume of a hypersphere. One can also define hyperellipsoids 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.

See also

  • Ellipsoid method
  • Ellipsoidal coordinates
  • Elliptical distribution, in statistics
  • Flattening, also called ellipticity and oblateness, is a measure of the compression of a circle or sphere along a diameter to form an ellipse or an ellipsoid of revolution (spheroid), respectively.
  • Focaloid, a shell bounded by two concentric, confocal ellipsoids
  • Geodesics on an ellipsoid
  • Geodetic datum, the gravitational Earth modeled by a best fitted elipsoid
  • Homoeoid, a shell bounded by two concentric, similar ellipsoids
  • List of surfaces

Notes

  1. Kreyszig (1972, pp. 455–456)
  2. F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, editors, 2010, NIST Handbook of Mathematical Functions (Cambridge University Press), available on line at "Archived copy". Archived from the original on 2012-12-02. Retrieved 2012-01-08.CS1 maint: archived copy as title (link) (see next reference).
  3. NIST (National Institute of Standards and Technology) at http://www.nist.gov Archived 2015-06-17 at the Wayback Machine
  4. W., Weisstein, Eric. "Prolate Spheroid". mathworld.wolfram.com. Archived from the original on 3 August 2017. Retrieved 25 March 2018.
  5. Final answers Archived 2011-09-30 at the Wayback Machine by Gerard P. Michon (2004-05-13). See Thomsen's formulas and Cantrell's comments.
  6. Albert, Abraham Adrian (2016) [1949], Solid Analytic Geometry, Dover, p. 117, ISBN 978-0-486-81026-3
  7. W. Böhm: Die FadenKonstruktion der Flächen zweiter Ordnung, Mathemat. Nachrichten 13, 1955, S. 151
  8. Staude, O.: Ueber Fadenconstructionen des Ellipsoides. Math. Ann. 20, 147–184 (1882)
  9. Staude, O.: Ueber neue Focaleigenschaften der Flächen 2. Grades. Math. Ann. 27, 253–271 (1886).
  10. Staude, O.: Die algebraischen Grundlagen der Focaleigenschaften der Flächen 2. Ordnung Math. Ann. 50, 398 - 428 (1898).
  11. D. Hilbert & S Cohn-Vossen: Geometry and the imagination, Chelsea New York, 1952, ISBN 0-8284-1087-9, p. 20 .
  12. "Archived copy" (PDF). Archived (PDF) from the original on 2013-06-26. Retrieved 2013-10-12.CS1 maint: archived copy as title (link) pp. 17–18.
  13. Computerunterstützte Darstellende und Konstruktive Geometrie. Archived 2013-11-10 at the Wayback Machine Uni Darmstadt (PDF; 3,4 MB), S. 88.
  14. Bezinque, Adam; et al. (2018). "Determination of Prostate Volume: A Comparison of Contemporary Methods". Academic Radiology. 25 (12): 1582–1587. doi:10.1016/j.acra.2018.03.014. PMID 29609953.
  15. Goldstein, H G (1980). Classical Mechanics, (2nd edition) Chapter 5.
  16. Dusenbery, David B. (2009).Living at Micro Scale, Harvard University Press, Cambridge, Massachusetts ISBN 978-0-674-03116-6.
  17. Frahm, G., Junker, M., & Szimayer, A. (2003). Elliptical copulas: applicability and limitations. Statistics & Probability Letters, 63(3), 275–286.

References

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.