Meridian arc

Last updated

In geodesy and navigation, a meridian arc is the curve between two points on the Earth's surface having the same longitude. The term may refer either to a segment of the meridian, or to its length.

Contents

The purpose of measuring meridian arcs is to determine a figure of the Earth. One or more measurements of meridian arcs can be used to infer the shape of the reference ellipsoid that best approximates the geoid in the region of the measurements. Measurements of meridian arcs at several latitudes along many meridians around the world can be combined in order to approximate a geocentric ellipsoid intended to fit the entire world.

The earliest determinations of the size of a spherical Earth required a single arc. Accurate survey work beginning in the 19th century required several arc measurements in the region the survey was to be conducted, leading to a proliferation of reference ellipsoids around the world. The latest determinations use astro-geodetic measurements and the methods of satellite geodesy to determine reference ellipsoids, especially the geocentric ellipsoids now used for global coordinate systems such as WGS 84 (see numerical expressions).

History of measurement

Early estimations of Earth's size are recorded from Greece in the 4th century BC, and from scholars at the caliph's House of Wisdom in Baghdad in the 9th century. The first realistic value was calculated by Alexandrian scientist Eratosthenes about 240 BC. He estimated that the meridian has a length of 252,000 stadia, with an error on the real value between -2.4% and +0.8% (assuming a value for the stadion between 155 and 160 metres). [1] Eratosthenes described his technique in a book entitled On the measure of the Earth, which has not been preserved. A similar method was used by Posidonius about 150 years later, and slightly better results were calculated in 827 by the arc measurement method, [2] attributed to the Caliph Al-Ma'mun.[ citation needed ]

Ellipsoidal Earth

Early literature uses the term oblate spheroid to describe a sphere "squashed at the poles". Modern literature uses the term ellipsoid of revolution in place of spheroid, although the qualifying words "of revolution" are usually dropped. An ellipsoid that is not an ellipsoid of revolution is called a triaxial ellipsoid. Spheroid and ellipsoid are used interchangeably in this article, with oblate implied if not stated.

17th and 18th centuries

Although it had been known since classical antiquity that the Earth was spherical, by the 17th century, evidence was accumulating that it was not a perfect sphere. In 1672, Jean Richer found the first evidence that gravity was not constant over the Earth (as it would be if the Earth were a sphere); he took a pendulum clock to Cayenne, French Guiana and found that it lost 2+12 minutes per day compared to its rate at Paris. [3] [4] This indicated the acceleration of gravity was less at Cayenne than at Paris. Pendulum gravimeters began to be taken on voyages to remote parts of the world, and it was slowly discovered that gravity increases smoothly with increasing latitude, gravitational acceleration being about 0.5% greater at the geographical poles than at the Equator.

In 1687, Isaac Newton had published in the Principia as a proof that the Earth was an oblate spheroid of flattening equal to 1/230. [5] This was disputed by some, but not all, French scientists. A meridian arc of Jean Picard was extended to a longer arc by Giovanni Domenico Cassini and his son Jacques Cassini over the period 1684–1718. [6] The arc was measured with at least three latitude determinations, so they were able to deduce mean curvatures for the northern and southern halves of the arc, allowing a determination of the overall shape. The results indicated that the Earth was a prolate spheroid (with an equatorial radius less than the polar radius). To resolve the issue, the French Academy of Sciences (1735) undertook expeditions to Peru (Bouguer, Louis Godin, de La Condamine, Antonio de Ulloa, Jorge Juan) and to Lapland (Maupertuis, Clairaut, Camus, Le Monnier, Abbe Outhier, Anders Celsius). The resulting measurements at equatorial and polar latitudes confirmed that the Earth was best modelled by an oblate spheroid, supporting Newton. [6] However, by 1743, Clairaut's theorem had completely supplanted Newton's approach.

By the end of the century, Jean Baptiste Joseph Delambre had remeasured and extended the French arc from Dunkirk to the Mediterranean Sea (the meridian arc of Delambre and Méchain). It was divided into five parts by four intermediate determinations of latitude. By combining the measurements together with those for the arc of Peru, ellipsoid shape parameters were determined and the distance between the Equator and pole along the Paris Meridian was calculated as 5130762  toises as specified by the standard toise bar in Paris. Defining this distance as exactly 10000000 m led to the construction of a new standard metre bar as 0.5130762 toises. [6] :22

19th century

In the 19th century, many astronomers and geodesists were engaged in detailed studies of the Earth's curvature along different meridian arcs. The analyses resulted in a great many model ellipsoids such as Plessis 1817, Airy 1830, Bessel 1841, Everest 1830, and Clarke 1866. [7] A comprehensive list of ellipsoids is given under Earth ellipsoid.

The nautical mile

Historically a nautical mile was defined as the length of one minute of arc along a meridian of a spherical earth. An ellipsoid model leads to a variation of the nautical mile with latitude. This was resolved by defining the nautical mile to be exactly 1,852 metres. However, for all practical purposes, distances are measured from the latitude scale of charts. As the Royal Yachting Association says in its manual for day skippers: "1 (minute) of Latitude = 1 sea mile", followed by "For most practical purposes distance is measured from the latitude scale, assuming that one minute of latitude equals one nautical mile". [8]

Calculation

On a sphere, the meridian arc length is simply the circular arc length. On an ellipsoid of revolution, for short meridian arcs, their length can be approximated using the Earth's meridional radius of curvature and the circular arc formulation. For longer arcs, the length follows from the subtraction of two meridian distances, the distance from the equator to a point at a latitude φ. This is an important problem in the theory of map projections, particularly the transverse Mercator projection.

The main ellipsoidal parameters are, a, b, f, but in theoretical work it is useful to define extra parameters, particularly the eccentricity, e, and the third flattening n. Only two of these parameters are independent and there are many relations between them:

Definition

The meridian radius of curvature can be shown to be equal to: [9] [10]

The arc length of an infinitesimal element of the meridian is dm = M(φ) (with φ in radians). Therefore, the meridian distance from the equator to latitude φ is

The distance formula is simpler when written in terms of the parametric latitude,

where tan β = (1 − f)tan φ and e2 = e2/1 − e2.

Even though latitude is normally confined to the range [−π/2,π/2], all the formulae given here apply to measuring distance around the complete meridian ellipse (including the anti-meridian). Thus the ranges of φ, β, and the rectifying latitude μ, are unrestricted.

Relation to elliptic integrals

The above integral is related to a special case of an incomplete elliptic integral of the third kind. In the notation of the online NIST handbook [11] (Section 19.2(ii)),

It may also be written in terms of incomplete elliptic integrals of the second kind (See the NIST handbook Section 19.6(iv)),

The calculation (to arbitrary precision) of the elliptic integrals and approximations are also discussed in the NIST handbook. These functions are also implemented in computer algebra programs such as Mathematica [12] and Maxima. [13]

Series expansions

The above integral may be expressed as an infinite truncated series by expanding the integrand in a Taylor series, performing the resulting integrals term by term, and expressing the result as a trigonometric series. In 1755, Leonhard Euler derived an expansion in the third eccentricity squared. [14]

Expansions in the eccentricity (e)

Delambre in 1799 [15] derived a widely used expansion on e2,

where

Richard Rapp gives a detailed derivation of this result. [16]

Expansions in the third flattening (n)

Series with considerably faster convergence can be obtained by expanding in terms of the third flattening n instead of the eccentricity. They are related by

In 1837, Friedrich Bessel obtained one such series, [17] which was put into a simpler form by Helmert, [18] [19]

with

Because n changes sign when a and b are interchanged, and because the initial factor 1/2(a + b) is constant under this interchange, half the terms in the expansions of H2k vanish.

The series can be expressed with either a or b as the initial factor by writing, for example,

and expanding the result as a series in n. Even though this results in more slowly converging series, such series are used in the specification for the transverse Mercator projection by the National Geospatial-Intelligence Agency [20] and the Ordnance Survey of Great Britain. [21]

Series in terms of the parametric latitude

In 1825, Bessel [22] derived an expansion of the meridian distance in terms of the parametric latitude β in connection with his work on geodesics,

with

Because this series provides an expansion for the elliptic integral of the second kind, it can be used to write the arc length in terms of the geodetic latitude as

Generalized series

The above series, to eighth order in eccentricity or fourth order in third flattening, provide millimetre accuracy. With the aid of symbolic algebra systems, they can easily be extended to sixth order in the third flattening which provides full double precision accuracy for terrestrial applications.

Delambre [15] and Bessel [22] both wrote their series in a form that allows them to be generalized to arbitrary order. The coefficients in Bessel's series can expressed particularly simply

where

and k!! is the double factorial, extended to negative values via the recursion relation: (−1)!! = 1 and (−3)!! = −1.

The coefficients in Helmert's series can similarly be expressed generally by

This result was conjectured by Friedrich Helmert [23] and proved by Kazushige Kawase. [24]

The extra factor (1 − 2k)(1 + 2k) originates from the additional expansion of appearing in the above formula and results in poorer convergence of the series in terms of φ compared to the one in β.

Numerical expressions

The trigonometric series given above can be conveniently evaluated using Clenshaw summation. This method avoids the calculation of most of the trigonometric functions and allows the series to be summed rapidly and accurately. The technique can also be used to evaluate the difference m(φ1) − m(φ2) while maintaining high relative accuracy.

Substituting the values for the semi-major axis and eccentricity of the WGS84 ellipsoid gives

where φ(°) = φ/ is φ expressed in degrees (and similarly for β(°)).

On the ellipsoid the exact distance between parallels at φ1 and φ2 is m(φ1) − m(φ2). For WGS84 an approximate expression for the distance Δm between the two parallels at ±0.5° from the circle at latitude φ is given by

Quarter meridian

A quarter meridian or Earth quadrant. Longitudinaler Erdquadrant.svg
A quarter meridian or Earth quadrant.

The distance from the equator to the pole, the quarter meridian (analogous to the quarter-circle), also known as the Earth quadrant, is

It was part of the historical definition of the metre and of the nautical mile, and used in the definition of the hebdomometre.

The quarter meridian can be expressed in terms of the complete elliptic integral of the second kind,

where are the first and second eccentricities.

The quarter meridian is also given by the following generalized series:

(For the formula of c0, see section #Generalized series above.) This result was first obtained by James Ivory. [25]

The numerical expression for the quarter meridian on the WGS84 ellipsoid is

The polar Earth's circumference is simply four times quarter meridian:

The perimeter of a meridian ellipse can also be rewritten in the form of a rectifying circle perimeter, Cp = 2πMr. Therefore, the rectifying Earth radius is:

It can be evaluated as 6367449.146 m.

The inverse meridian problem for the ellipsoid

In some problems, we need to be able to solve the inverse problem: given m, determine φ. This may be solved by Newton's method, iterating

until convergence. A suitable starting guess is given by φ0 = μ where

is the rectifying latitude. Note that it there is no need to differentiate the series for m(φ), since the formula for the meridian radius of curvature M(φ) can be used instead.

Alternatively, Helmert's series for the meridian distance can be reverted to give [26] [27]

where

Similarly, Bessel's series for m in terms of β can be reverted to give [28]

where

Adrien-Marie Legendre showed that the distance along a geodesic on a spheroid is the same as the distance along the perimeter of an ellipse. [29] For this reason, the expression for m in terms of β and its inverse given above play a key role in the solution of the geodesic problem with m replaced by s, the distance along the geodesic, and β replaced by σ, the arc length on the auxiliary sphere. [22] [30] The requisite series extended to sixth order are given by Charles Karney, [31] Eqs. (17) & (21), with ε playing the role of n and τ playing the role of μ.

See also

Related Research Articles

<span class="mw-page-title-main">Latitude</span> Geographic coordinate specifying north–south position

In geography, latitude is a coordinate that specifies the north–south position of a point on the surface of the Earth or another celestial body. Latitude is given as an angle that ranges from –90° at the south pole to 90° at the north pole, with 0° at the Equator. Lines of constant latitude, or parallels, run east–west as circles parallel to the equator. Latitude and longitude are used together as a coordinate pair to specify a location on the surface of the Earth.

<span class="mw-page-title-main">Mercator projection</span> Cylindrical conformal map projection

The Mercator projection is a cylindrical map projection presented by Flemish geographer and cartographer Gerardus Mercator in 1569. It became the standard map projection for navigation because it is unique in representing north as up and south as down everywhere while preserving local directions and shapes. The map is thereby conformal. As a side effect, the Mercator projection inflates the size of objects away from the equator. This inflation is very small near the equator but accelerates with increasing latitude to become infinite at the poles. As a result, landmasses such as Greenland, Antarctica, Canada and Russia appear far larger than they actually are relative to landmasses near the equator, such as Central Africa.

<span class="mw-page-title-main">Pauli matrices</span> Matrices important in quantum mechanics and the study of spin

In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices that are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

<span class="mw-page-title-main">Euler's totient function</span> Number of integers coprime to and not exceeding n

In number theory, Euler's totient function counts the positive integers up to a given integer n that are relatively prime to n. It is written using the Greek letter phi as or , and may also be called Euler's phi function. In other words, it is the number of integers k in the range 1 ≤ kn for which the greatest common divisor gcd(n, k) is equal to 1. The integers k of this form are sometimes referred to as totatives of n.

<span class="mw-page-title-main">Ellipsoid</span> Quadric surface that looks like a deformed sphere

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

<span class="mw-page-title-main">Fermat's spiral</span> Spiral that surrounds equal area per turn

A Fermat's spiral or parabolic spiral is a plane curve with the property that the area between any two consecutive full turns around the spiral is invariant. As a result, the distance between turns grows in inverse proportion to their distance from the spiral center, contrasting with the Archimedean spiral and the logarithmic spiral. Fermat spirals are named after Pierre de Fermat.

<span class="mw-page-title-main">Rhumb line</span> Arc crossing all meridians of longitude at the same angle

In navigation, a rhumb line, rhumb, or loxodrome is an arc crossing all meridians of longitude at the same angle, that is, a path with constant bearing as measured relative to true north.

<span class="mw-page-title-main">Cardioid</span> Type of plane curve

In geometry, a cardioid is a plane curve traced by a point on the perimeter of a circle that is rolling around a fixed circle of the same radius. It can also be defined as an epicycloid having a single cusp. It is also a type of sinusoidal spiral, and an inverse curve of the parabola with the focus as the center of inversion. A cardioid can also be defined as the set of points of reflections of a fixed point on a circle through all tangents to the circle.

<span class="mw-page-title-main">Stable distribution</span> Distribution of variables which satisfies a stability property under linear combinations

In probability theory, a distribution is said to be stable if a linear combination of two independent random variables with this distribution has the same distribution, up to location and scale parameters. A random variable is said to be stable if its distribution is stable. The stable distribution family is also sometimes referred to as the Lévy alpha-stable distribution, after Paul Lévy, the first mathematician to have studied it.

<span class="mw-page-title-main">Universal Transverse Mercator coordinate system</span> Map projection system

The Universal Transverse Mercator (UTM) is a map projection system for assigning coordinates to locations on the surface of the Earth. Like the traditional method of latitude and longitude, it is a horizontal position representation, which means it ignores altitude and treats the earth surface as a perfect ellipsoid. However, it differs from global latitude/longitude in that it divides earth into 60 zones and projects each to the plane as a basis for its coordinates. Specifying a location means specifying the zone and the x, y coordinate in that plane. The projection from spheroid to a UTM zone is some parameterization of the transverse Mercator projection. The parameters vary by nation or region or mapping system.

<span class="mw-page-title-main">Tissot's indicatrix</span> Characterization of distortion in map protections

In cartography, a Tissot's indicatrix is a mathematical contrivance presented by French mathematician Nicolas Auguste Tissot in 1859 and 1871 in order to characterize local distortions due to map projection. It is the geometry that results from projecting a circle of infinitesimal radius from a curved geometric model, such as a globe, onto a map. Tissot proved that the resulting diagram is an ellipse whose axes indicate the two principal directions along which scale is maximal and minimal at that point on the map.

<span class="mw-page-title-main">Great-circle navigation</span> Flight or sailing route along the shortest path between two points on a globes surface

Great-circle navigation or orthodromic navigation is the practice of navigating a vessel along a great circle. Such routes yield the shortest distance between two points on the globe.

<span class="mw-page-title-main">Geographical distance</span> Distance measured along the surface of the Earth

Geographical distance or geodetic distance is the distance measured along the surface of the Earth, or the shortest arch length.

Transverse Mercator projection has many implementations. Louis Krüger in 1912 developed one of his two implementations that expressed as a power series in the longitude difference from the central meridian. These series were recalculated by Lee in 1946, by Redfearn in 1948, and by Thomas in 1952. They are often referred to as the Redfearn series, or the Thomas series. This implementation is of great importance since it is widely used in the U.S. State Plane Coordinate System, in national and also international mapping systems, including the Universal Transverse Mercator coordinate system (UTM). They are also incorporated into the Geotrans coordinate converter made available by the United States National Geospatial-Intelligence Agency. When paired with a suitable geodetic datum, the series deliver high accuracy in zones less than a few degrees in east-west extent.

In geodesy and geophysics, theoretical gravity or normal gravity is an approximation of the true gravity on Earth's surface by means of a mathematical model representing Earth. The most common model of a smoothed Earth is a rotating Earth ellipsoid of revolution.

In mathematics, infinite compositions of analytic functions (ICAF) offer alternative formulations of analytic continued fractions, series, products and other infinite expansions, and the theory evolving from such compositions may shed light on the convergence/divergence of these expansions. Some functions can actually be expanded directly as infinite compositions. In addition, it is possible to use ICAF to evaluate solutions of fixed point equations involving infinite expansions. Complex dynamics offers another venue for iteration of systems of functions rather than a single function. For infinite compositions of a single function see Iterated function. For compositions of a finite number of functions, useful in fractal theory, see Iterated function system.

<span class="mw-page-title-main">Geodesics on an ellipsoid</span> Shortest paths on a bounded deformed sphere-like quadric surface

The study of geodesics on an ellipsoid arose in connection with geodesy specifically with the solution of triangulation networks. The figure of the Earth is well approximated by an oblate ellipsoid, a slightly flattened sphere. A geodesic is the shortest path between two points on a curved surface, analogous to a straight line on a plane surface. The solution of a triangulation network on an ellipsoid is therefore a set of exercises in spheroidal trigonometry.

<span class="mw-page-title-main">Earth section paths</span> Plane curved by the intersection of an earth ellipsoid and a plane

Earth section paths are plane curves defined by the intersection of an earth ellipsoid and a plane. Common examples include the great ellipse and normal sections. Earth section paths are useful as approximate solutions for geodetic problems, the direct and inverse calculation of geographic distances. The rigorous solution of geodetic problems involves skew curves known as geodesics.

References

  1. Russo, Lucio (2004). The Forgotten Revolution . Berlin: Springer. p.  273-277.
  2. Torge, W.; Müller, J. (2012). Geodesy. De Gruyter Textbook. De Gruyter. p. 5. ISBN   978-3-11-025000-8 . Retrieved 2021-05-02.
  3. Poynting, John Henry; Joseph John Thompson (1907). A Textbook of Physics, 4th Ed. London: Charles Griffin & Co. p.  20.
  4. Victor F., Lenzen; Robert P. Multauf (1964). "Paper 44: Development of gravity pendulums in the 19th century". United States National Museum Bulletin 240: Contributions from the Museum of History and Technology reprinted in Bulletin of the Smithsonian Institution. Washington: Smithsonian Institution Press. p. 307. Retrieved 2009-01-28.
  5. Isaac Newton: Principia, Book III, Proposition XIX, Problem III, translated into English by Andrew Motte. A searchable modern translation is available at 17centurymaths. Search the following pdf file for 'spheroid'.
  6. 1 2 3 Clarke, Alexander Ross (1880). Geodesy. Oxford: Clarendon Press. OCLC   2484948.. Freely available online at Archive.org and Forgotten Books ( ISBN   9781440088650). In addition the book has been reprinted by Nabu Press ( ISBN   978-1286804131), the first chapter covers the history of early surveys.
  7. Clarke, Alexander Ross; James, Henry (1866). Comparisons of the standards of length of England, France, Belgium, Prussia, Russia, India, Australia, made at the Ordnance survey office, Southampton. London: G.E. Eyre and W. Spottiswoode for H.M. Stationery Office. pp. 281–87. OCLC   906501. Appendix on Figure of the Earth.
  8. Hopkinson, Sara (2012). RYA day skipper handbook - sail. Hamble: The Royal Yachting Association. p. 76. ISBN   9781-9051-04949.
  9. Rapp, R, (1991): Geometric Geodesy, Part I, §3.5.1, pp. 28–32.
  10. Osborne, Peter (2013), The Mercator Projections, doi:10.5281/zenodo.35392 Section 5.6. This reference includes the derivation of curvature formulae from first principles and a proof of Meusnier's theorem. (Supplements: Maxima files and Latex code and figures)
  11. F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, editors, 2010, NIST Handbook of Mathematical Functions (Cambridge University Press).
  12. Mathematica guide: Elliptic Integrals
  13. Maxima, 2009, A computer algebra system, version 5.20.1.
  14. Euler, L. (1755). "Élémens de la trigonométrie sphéroïdique tirés de la méthode des plus grands et plus petits" [Elements of spheroidal trigonometry taken from the method of maxima and minima]. Mémoires de l'Académie Royale des Sciences de Berlin 1753 (in French). 9: 258–293. Figures.
  15. 1 2 Delambre, J. B. J. (1799): Méthodes Analytiques pour la Détermination d'un Arc du Méridien; précédées d'un mémoire sur le même sujet par A. M. Legendre, De L'Imprimerie de Crapelet, Paris, 72–73
  16. Rapp, R, (1991), §3.6, pp. 36–40.
  17. Bessel, F. W. (1837). "Bestimmung der Axen des elliptischen Rotationssphäroids, welches den vorhandenen Messungen von Meridianbögen der Erde am meisten entspricht" [Estimation of the axes of the ellipsoid through measurements of the meridian arc]. Astronomische Nachrichten (in German). 14 (333): 333–346. Bibcode:1837AN.....14..333B. doi:10.1002/asna.18370142301.
  18. Helmert, F. R. (1880): Die mathematischen und physikalischen Theorieen der höheren Geodäsie, Einleitung und 1 Teil, Druck und Verlag von B. G. Teubner, Leipzig, § 1.7, pp. 44–48. English translation (by the Aeronautical Chart and Information Center, St. Louis) available at doi : 10.5281/zenodo.32050
  19. Krüger, L. (1912): Konforme Abbildung des Erdellipsoids in der Ebene . Royal Prussian Geodetic Institute, New Series 52, page 12
  20. J. W. Hager, J.F. Behensky, and B.W. Drew, 1989. Defense Mapping Agency Technical Report TM 8358.2. The universal grids: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS)
  21. A guide to coordinate systems in Great Britain, Ordnance Survey of Great Britain.
  22. 1 2 3 Bessel, F. W. (2010). "The calculation of longitude and latitude from geodesic measurements (1825)". Astron. Nachr. 331 (8): 852–861. arXiv: 0908.1824 . Bibcode:2010AN....331..852K. doi:10.1002/asna.201011352. S2CID   118760590. English translation of Astron. Nachr. 4, 241–254 (1825), §5.
  23. Helmert (1880), §1.11
  24. Kawase, K. (2011): A General Formula for Calculating Meridian Arc Length and its Application to Coordinate Conversion in the Gauss-Krüger Projection, Bulletin of the Geospatial Information Authority of Japan, 59, 1–13
  25. Ivory, J. (1798). "A new series for the rectification of the ellipsis". Transactions of the Royal Society of Edinburgh . 4 (2): 177–190. doi:10.1017/s0080456800030817. S2CID   251572677.
  26. Helmert (1880), §1.10
  27. Adams, Oscar S (1921). Latitude Developments Connected With Geodesy and Cartography. US Coast and Geodetic Survey Special Publication No. 67. p. 127.
  28. Helmert (1880), §5.6
  29. Legendre, A. M. (1811). Exercices de Calcul Intégral sur Divers Ordres de Transcendantes et sur les Quadratures [Exercises in Integral Calculus] (in French). Paris: Courcier. p.  180. OCLC   312469983.
  30. Helmert (1880), Chap. 5
  31. Karney, C. F. F. (2013). "Algorithms for geodesics". Journal of Geodesy. 87 (1): 43–55. arXiv: 1109.4448 . Bibcode:2013JGeod..87...43K. doi:10.1007/s00190-012-0578-z. S2CID   119310141. Open Access logo PLoS transparent.svg Addenda.