The Geometric Algorithm of Inverse and Direct Problems with an Area Solution for the Great Elliptic Arcs

Total Page:16

File Type:pdf, Size:1020Kb

The Geometric Algorithm of Inverse and Direct Problems with an Area Solution for the Great Elliptic Arcs Journal of Marine Science and Technology, Vol. 23, No. 4, pp. 481-490 (2015) 481 DOI: 10.6119/JMST-014-1120-1 THE GEOMETRIC ALGORITHM OF INVERSE AND DIRECT PROBLEMS WITH AN AREA SOLUTION FOR THE GREAT ELLIPTIC ARCS Wei-Kuo Tseng, Jiunn-Liang Guo, and Chung-Ping Liu Key words: great ellipse, edge, geography type, computer algebra The result of any computation, e.g. the length of a path or system. the intersection of polygons, depends on the definition of the edges that connect these points. On a planar map, the edge between two points is obviously the line segment that connects ABSTRACT them, but on an ellipsoidal earth model the choice is not ob- The definition of the connecting edge(s) between two ver- vious, and it varies between different software products. While tices in a geography type is the great elliptic minor arc in differences in accuracy and performance are to be expected, it Microsoft SQL Server’s Geography Type. Similar to the way is a sad state of affairs when different software packages dis- an edge is defined by Microsoft SQL Server, an edge in agree on the theoretical results of their computations. The Geodyssey’s Hipparchus library is also defined as a great paper (Kallay, 2007) presents the definition of edges in Mi- circle arc on a reference sphere. Using Hipparchus for their crosoft’s SQL Server’s Geography Type, proposing it as an computations, IBM’s DB2 Geodetic Extender and Informix industry standard. It stands to reason that a round earth edge Geodetic Datablade share this definition. Compact formulae should satisfy the following requirements: are given for the great elliptic sailing on a spheroid providing solutions to both the forward and inverse problems with ex- 1. Locally, an edge should be experienced as straight. ceptional accuracy, and latitude in terms of longitude. The 2. A pair of points should define a unique edge between them. solution incorporates a closed form for the azimuth and the 3. An edge should admit a differentiable parameterization, derivation of the algorithm is presented and illustrated. In which assigns a point on the edge to every real number addition, the area of polygon bounded by the elliptic arcs is between 0 and 1. treated. This paper also shows that a computer algebra system is a powerful tool to solve mathematical derivations in navi- The geodesic is the curve on the surface of an ellipsoid de- gation, geodesy, and cartography. fining the shortest distance between two points. Kallay (2007) points out that geodesic curves score poorly on the require- ment 2, 3, and even on 1 they are not the obvious choice: I. INTRODUCTION Lines, polygonal paths and polygons are widely used in the 1. While the geodesic curve is the shortest path that is con- description of geospatial data, and they are usually defined in fined to the surface, most human activities take us beyond terms of their endpoints and vertices. The definition of the the surface, for example, surveyors measure along straight connecting edge between two vertices is the shorter great el- lines of sight and airplanes fly miles above the surface of liptic arc in Microsoft SQL Server’s Geography Type (Kallay, the globe. 2007; Microsoft, 2013). Similar to the way an edge is defined 2. There are numerous (not necessarily antipodal) pairs of by Microsoft SQL Server, an edge in Geodyssey’s Hipparchus points on an ellipsoid between which there are more than library is also defined as a great circle arc on a reference sphere one short geodesic (Rapp, 1991). (Geodyssey, 2013). Using Hipparchus for their computations, 3. Computing points along geodesic curves is notoriously IBM’s DB2 Geodetic Extender and Informix Geodetic Datab- difficult and expensive. An exact differentiable parame- lade share this definition (IBM, 2013). terization is not known, and approximate ones are also dif- ficult and expensive to compute. Paper submitted 03/05/14; revised 08/07/14; accepted 11/20/14. Author for Classical surveying suggests the definition of an edge as a correspondence: Wei-Kuo Tseng (e-mail: [email protected]). 1Department of Merchant Marine, National Taiwan Ocean University, Kee- normal section (Rapp, 1991), which is a plane curve created lung, Taiwan, R.O.C. by intersecting a plane containing the normal to the spheroid 482 Journal of Marine Science and Technology, Vol. 23, No. 4 (2015) However, the existed formulae need cumbersome algorithms Great Ellipse and their accuracies are not very high. In addition, the mathe- matical derivations in those literatures are a bit tedious, and abstruse, hardly suited for coding (Bowring, 1984; Pallikaris Great Circle and Latsas, 2009). The direct solution was also not completely a provided in those articles (Earle, 2011). For these reasons, in b P2 P Auxiliary Spere this paper we revisit the solution for the great elliptic arc and P1 Spheroid provide a more straightforward and compact mathematical φ a derivation of the spherical trigonometric solutions. This paper λ a a also gives a general inverse and direct solution attaining any accuracy requirement for the calculation of the great ellipse R(φ) ()cosφ cosλ, cosφ sinλ, sinφ sailing. In the mathematical derivation, we consider the direct and Fig. 1. A great ellipse on a spheroid. inverse scenarios to produce solutions determining the great ellipse from one point and its azimuth or between two points. with the surface of the spheroid. Alas, this definition is rarely Furthermore, the interpolation for latitude in terms of longi- unique. The surveyor’s plane at the other endpoint may define tude between end points of a great ellipse on the spheroid has a different normal section. not yet been found in the literature. As a consequence of these The paper (Kallay, 2007) evaluates the definition of the observations, the complete solution to the great ellipse pre- great elliptic arc against the above stated requirements: sented here will include a method to determine latitude for any specified longitude along the ellipse. Because the calculation 1. Edges are experienced as straight or approximately straight of the area of polygon bounded by the geodesics needs cum- in several senses: As a great elliptic arc, an edge is planar. bersome algorithms (Sjöberg, 2006), the alternative calcula- The angular deviation of its plane from a surveyor’s planes tion of the area of polygon bounded by the great elliptic arc is at either endpoint is no more than 12’. This translates to also provided here. The accuracies attained can satisfy the about 2.8 cm for an edge whose length is 10 km. An edge is requirement of ECDIS and GIS environments. Finally, we approximately the shortest path between its endpoints – the give the full formulae of spherical trigonometry that can easily length of no edge exceeds the geodesic distance by more be coded in a programming language so that readers should than 0.02%. In the space of direction, the edge is a line comprehensively grasp the meaning of the geometry. segment, and so are its gnomonic projections. 2. Every pair of non-antipodal points defines a unique edge. II. PARAMETERS OF THE GREAT ELLIPSE 3. The parameterization is simple and differentiable. The parameterization as a quadratic rational Bezier curve may Using geodetic and geocentric latitudes, a point P on the be slightly more expensive to set up but very efficient for surface of a spheroid such as Earth can be represented as a generating multiple points along the edge. vector function of longitude , geodetic latitude , or geo- centric latitude . The great elliptic arc on spheroid has been investigated in (Bowring, 1984; Walwyn, 1990; Williams, 1998; Earle, 2000, P(,) xyz 2008; Kally, 2007; Tseng and Lee 2010), but is rarely men- 2 tioned elsewhere. The great elliptic arc between two points P1 Ne()coscos, cossin, (1 )sin (1) and P2 on a spheroid, centered at O, is the minor arc of the ellipse of the intersection between the spheroid and the plane OP1P2 (Fig. 1). If the two points are antipodal, the collinear and points P1, O, and P2 do not determine a unique plane, in such a case it would be reasonable to choose the route passing through PxyzR( , ) ( ) cos cos , cos sin , sin the two poles of the spheroid. The azimuth at the point P1 is (2) the angle that the tangent at P1 to great ellipse P1P2 makes with 221/2 the meridian through P1, and is measured from the clockwise where e is the eccentricity, and Nae( ) /(1 sin ) is direction northerly. The azimuth at arbitrary points on the the radius of curvature of the prime vertical, and great ellipse would be similarly defined (Bowring, 1984). Some approximate formulae, the great elliptic equation and 1/2 Rae( ) (1222 ) /(1 e cos ) . (3) great circle equations have been provided in a number of pa- pers (Earle, 2000; Pallikaris and Latsas, 2009; Tseng and Lee 2007a, 2007b, 2010, 2012, 2013) that have studied this prob- The following equation can transform Cartesian coordi- lem of the great ellipse sailing and achieved remarkable results. nates of a point on the spheroid to geodetic coordinates. W.-K. Tseng et al.: Geometric Algorithm of Inverse and Direct Problems with an Area Solution for the Great Elliptic Arcs 483 ' ' TTNNN=−⋅+−⋅cos(ϕφ ) sin( ϕφ ) TTN cos( ) sin( ) (8) N N NP α N ϕ − φ T p TN v The normal to spheroid at point P and the unit velocity γ vector at point P are orthogonal, so the inner product of the T P E two vectors equals to 0. Substitute Eq. (8) into Eq. (7) to obtain (Bowring, 1984; Earle, 2008): O ψ σ λ atan2TTVE , cos( ) TT VN (9) λe E Therefore the following relations exist: Fig.
Recommended publications
  • So, What Is Actually the Distance from the Equator to the Pole? – Overview of the Meridian Distance Approximations
    the International Journal Volume 7 on Marine Navigation Number 2 http://www.transnav.eu and Safety of Sea Transportation June 2013 DOI: 10.12716/1001.07.02.14 So, What is Actually the Distance from the Equator to the Pole? – Overview of the Meridian Distance Approximations A. Weintrit Gdynia Maritime University, Gdynia, Poland ABSTRACT: In the paper the author presents overview of the meridian distance approximations. He would like to find the answer for the question what is actually the distance from the equator to the pole ‐ the polar distance. In spite of appearances this is not such a simple question. The problem of determining the polar distance is a great opportunity to demonstrate the multitude of possible solutions in common use. At the beginning of the paper the author discusses some approximations and a few exact expressions (infinite sums) to calculate perimeter and quadrant of an ellipse, he presents convenient measurement units of the distance on the surface of the Earth, existing methods for the solution of the great circle and great elliptic sailing, and in the end he analyses and compares geodetic formulas for the meridian arc length. 1 INTRODUCTION navigational receivers and navigational systems (ECDIS and ECS [Weintrit, 2009]) suggest the Unfortunately, from the early days of the necessity of a thorough examination, modification, development of the basic navigational software built verification and unification of the issue of sailing into satellite navigational receivers and later into calculations for navigational systems and receivers. electronic chart systems, it has been noted that for the The problem of determining the distance from the sake of simplicity and a number of other, often equator to the pole is a great opportunity to incomprehensible reasons, this navigational software demonstrate the multitude of possible solutions in is often based on the simple methods of limited common use.
    [Show full text]
  • From Lacaille to Gill and the Start of the Arc of the 30 Meridian
    From LaCaille to Gill and the start of the Arc of the 30th Meridian Jim R. SMITH, United Kingdom Key words: LaCaille. Maclear. Gill. Meridian Arcs. 30th Arc. SUMMARY If one looks at a map of European arcs of meridian and parallel at the beginning of the 20th century there will be seen to be a plethora them. Turning to Africa there was the complete opposite. Other than the arcs of Eratosthenes (c 300 BC) and the short arcs by LaCaille (1752) and Maclear (1841-48) the continent was empty. It was in 1879 that David Gill had the idea for a Cape to Cairo meridian arc but 1954 before it was completed. This presentation summaries the work of LaCaille and Maclear and then concentrates on a description of the 30th Meridian Arc in East Africa. The only other comparable arcs at the time were that through the centre of India by Lambton and Everest observed between 1800 and 1843; the Struve arc of 1816 to 1852 and the various arcs in France during the 18th century. The usefulness of such arcs was highlighted in 2005 with the inscription by UNESCO of the Struve Geodetic Arc on the World Heritage Monument list. A practical extension to the Struve Arc Monument is the 30th Meridian Arc since there is a connection between the two. At an ICA Symposium in Cape Town, 2003, Lindsay Braun gave a presentation that detailed the political machinations of the history of the 30th Arc; here it is hoped to fill in other aspects of this work including some facts and figures since that is what surveyors thrive upon.
    [Show full text]
  • Geodesy Methods Over Time
    Geodesy methods over time GISC3325 - Lecture 4 Astronomic Latitude/Longitude • Given knowledge of the positions of an astronomic body e.g. Sun or Polaris and our time we can determine our location in terms of astronomic latitude and longitude. US Meridian Triangulation This map shows the first project undertaken by the founding Superintendent of the Survey of the Coast Ferdinand Hassler. Triangulation • Method of indirect measurement. • Angles measured at all nodes. • Scaled provided by one or more precisely measured base lines. • First attributed to Gemma Frisius in the 16th century in the Netherlands. Early surveying instruments Left is a Quadrant for angle measurements, below is how baseline lengths were measured. A non-spherical Earth • Willebrod Snell van Royen (Snellius) did the first triangulation project for the purpose of determining the radius of the earth from measurement of a meridian arc. • Snellius was also credited with the law of refraction and incidence in optics. • He also devised the solution of the resection problem. At point P observations are made to known points A, B and C. We solve for P. Jean Picard’s Meridian Arc • Measured meridian arc through Paris between Malvoisine and Amiens using triangulation network. • First to use a telescope with cross hairs as part of the quadrant. • Value obtained used by Newton to verify his law of gravitation. Ellipsoid Earth Model • On an expedition J.D. Cassini discovered that a one-second pendulum regulated at Paris needed to be shortened to regain a one-second oscillation. • Pendulum measurements are effected by gravity! Newton • Newton used measurements of Picard and Halley and the laws of gravitation to postulate a rotational ellipsoid as the equilibrium figure for a homogeneous, fluid, rotating Earth.
    [Show full text]
  • Great Elliptic Arc Distance1
    GREAT ELLIPTIC ARC DISTANCE1 R. E. Deakin School of Mathematical & Geospatial Sciences, RMIT University, GPO Box 2476V, MELBOURNE VIC 3001, AUSTRALIA email: [email protected] January 2012 ABSTRACT These notes provide a detailed derivation of series formula for (i) Meridian distance M as a function of latitude φ and the ellipsoid constant eccentricity-squared e2, and φ and the ellipsoid constant third flattening n; (ii) Rectifying latitude µ as functions of φ,e2 and φ,n; and (iii) latitude φ as functions of µ,e2 and µ,n These series can then be used in solving the direct and inverse problems on the ellipsoid (using great elliptic arcs) and are easily obtained using the Computer Algebra System Maxima. In addition, a detailed derivation of the equation for the great elliptic arc on an ellipsoid is provided as well as defining the azimuth and the vertex of a great elliptic. And to assist in the solution of the direct and inverse problems the auxiliary sphere is introduced and equations developed. INTRODUCTION In geodesy, the great elliptic arc between P1 and P2 on the ellipsoid is the curve created by intersecting the ellipsoid with the plane containing P1 , P2 and O (the centre of the ellipsoid) and these planes are great elliptic planes or sections. Figure 1 shows P on the great elliptic arc between P1 and P2 . P is the geocentric latitude of P and P is the longitude of P. There are an infinite number of planes that cut the surface of the ellipsoid and contain the chord PP12 but only one of these will contain the centre O.
    [Show full text]
  • The Meridian Arc Measurement in Peru 1735 – 1745
    The Meridian Arc Measurement in Peru 1735 – 1745 Jim R. SMITH, United Kingdom Key words: Peru. Meridian. Arc. Triangulation. ABSTRACT: In the early 18th century the earth was recognised as having some ellipsoidal shape rather than a true sphere. Experts differed as to whether the ellipsoid was flattened at the Poles or the Equator. The French Academy of Sciences decided to settle the argument once and for all by sending one expedition to Lapland- as near to the Pole as possible; and another to Peru- as near to the Equator as possible. The result supported the view held by Newton in England rather than that of the Cassinis in Paris. CONTACT Jim R. Smith, Secretary to International Institution for History of Surveying & Measurement 24 Woodbury Ave, Petersfield Hants GU32 2EE UNITED KINGDOM Tel. & fax + 44 1730 262 619 E-mail: [email protected] Website: http://www.ddl.org/figtree/hsm/index.htm HS4 Surveying and Mapping the Americas – In the Andes of South America 1/12 Jim R. Smith The Meridian Arc Measurement in Peru 1735-1745 FIG XXII International Congress Washington, D.C. USA, April 19-26 2002 THE MERIDIAN ARC MEASUREMENT IN PERU 1735 – 1745 Jim R SMITH, United Kingdom 1. BACKGROUND The story might be said to begin just after the mid 17th century when Jean Richer was sent to Cayenne, S. America, to carry out a range of scientific experiments that included the determination of the length of a seconds pendulum. He returned to Paris convinced that in Cayenne the pendulum needed to be 11 lines (2.8 mm) shorter there than in Paris to keep the same time.
    [Show full text]
  • The History of Geodesy Told Through Maps
    The History of Geodesy Told through Maps Prof. Dr. Rahmi Nurhan Çelik & Prof. Dr. Erol KÖKTÜRK 16 th May 2015 Sofia Missionaries in 5000 years With all due respect... 3rd FIG Young Surveyors European Meeting 1 SUMMARIZED CHRONOLOGY 3000 BC : While settling, people were needed who understand geometries for building villages and dividing lands into parts. It is known that Egyptian, Assyrian, Babylonian were realized such surveying techniques. 1700 BC : After floating of Nile river, land surveying were realized to set back to lost fields’ boundaries. (32 cm wide and 5.36 m long first text book “Papyrus Rhind” explain the geometric shapes like circle, triangle, trapezoids, etc. 550+ BC : Thereafter Greeks took important role in surveying. Names in that period are well known by almost everybody in the world. Pythagoras (570–495 BC), Plato (428– 348 BC), Aristotle (384-322 BC), Eratosthenes (275–194 BC), Ptolemy (83–161 BC) 500 BC : Pythagoras thought and proposed that earth is not like a disk, it is round as a sphere 450 BC : Herodotus (484-425 BC), make a World map 350 BC : Aristotle prove Pythagoras’s thesis. 230 BC : Eratosthenes, made a survey in Egypt using sun’s angle of elevation in Alexandria and Syene (now Aswan) in order to calculate Earth circumferences. As a result of that survey he calculated the Earth circumferences about 46.000 km Moreover he also make the map of known World, c. 194 BC. 3rd FIG Young Surveyors European Meeting 2 150 : Ptolemy (AD 90-168) argued that the earth was the center of the universe.
    [Show full text]
  • B10.Pdf (774.2Kb)
    Page 1 Coordinates Series B, No. 10 The History of Cartography in a Nutshell Persistent URL for citation: http://purl.oclc.org/ coordinates/b10.pdf Vladimiro Valerio Vladimiro Valerio (e-mail: [email protected]) is full Professor of Projective Geometry at the Date of Publication: 06/03/08 Department of History of Architecture at the University IUAV of Venice. His publications include "Cartography in the Kingdom of Naples During the Early Modern Period" in vol. 3 of The History of Cartography (University of Chicago Press). Abstract This is a very short history of cartography. Notes and links to images are included at the end. Keywords: maps, cartography, history, portolan charts, map projections, Claudius Ptolemy, Gerardus Mercator, Fra Mauro, Peutinger Table, C.F. Cassini de Thury Editor’s Note: About five years ago Professor Valerio was asked to prepare a short article on the history of cartography for a multimedia presentation by the Istituto e Museo di Storia della Scienza of Florence. He was astounded to learn that his article could only be thirteen lines long, but he nonetheless complied. The following is a translation of that article. Most of us are familiar with much longer accounts of this subject, such as the multi-volume History of Cartography being published by the University of Chicago Press, or the single volume edited by James R. Akerman and Robert W. Karrow, Jr., Maps: Finding Our Place in the World (Chicago: University of Page 2 Chicago Press, 2007). But, as the Reader’s Digest has taught us, there are few limits to condensation. Prof.
    [Show full text]
  • Meridian Distance.Pdf
    MERIDIAN DISTANCE R.E.Deakin School of Mathematical & Geospatial Sciences, RMIT University, GPO Box 2476V, MELBOURNE, Australia email: [email protected] March 2006 ABSTRACT These notes provide a detailed derivation of the equations for computing meridian distance on an ellipsoid of revolution given ellipsoid parameters and latitude. Computation of meridian distance is a requirement for geodetic calculations on the ellipsoid, in particular, computations to do with conversion of geodetic coordinates φλ, (latitude, longitude) to Universal Transverse Mercator (UTM) projection coordinates E,N (East,North) that are used in Australia for survey coordination. The "opposite" transformation, E,N to φλ, requires latitude given meridian distance and these note show how series equations for meridian distances are "reversed" to give series equations for latitude. The derivation of equations follow methods set out in Lehrbuch Der Geodäsie (Baeschlin, 1948), Handbuch der Vermessungskunde (Jordan/Eggert/Kneissl, 1958), Geometric Geodesy (Rapp, 1982) and Geodesy and Map Projections (Lauf, 1983) and some of the formula derived are used in the Geocentric Datum of Australia Technical Manual (ICSM, 2002) – an on-line reference manual available from Geoscience Australia. An understanding of the methods introduced in the following pages, in particular the solution of elliptic integrals by series expansion and reversion of a series, will give the student an insight into other geodetic calculations. MATLAB functions for computing (i) meridian distance given
    [Show full text]
  • On the Ancient Determination of the Meridian Arc Length by Eratosthenes of Kyrene
    On the Ancient Determination of the Meridian Arc Length by Eratosthenes of Kyrene Dieter LELGEMANN, Germany Key words: SUMMARY For topographic/geographic purposes the size of the Earth is of utmost importance in modern time as well as for Eratosthenes facing the task to draw a map of the "oikumene". His determination of the radius of the Earth resulted in 252000 stadia. If he has used a stadion definition of 1 stadion = 158,7m = 300 Egyptian royal cubits = 600 Gudea units (length of the yardstick at the statue of Gudea (2300 b c) in the Louvre/Paris), he has already observed the meridian arc length to 252000 ⋅ 0,1587 = 40000km. How could Eratosthenes obtain in ancient times already such an accurate result? Ptolemaios describes in his "Geographike hyphegesis" a method used by the "elder" to determine the size of the Earth; this ancient method to measure the meridian arc length between the latitude circles of two cities (e.g. Alexandria/Syene, Syene/Meroe) is based on a traversing technique, as will be shown. Geographical latitudes could be measured using a "Skiotheron" (shadow seizer). An according to ancient information reconstructed instrument will be shown and explained; the accuracy of sun observations with such a kind of instrument is comparable with those of a modern sextant. A recovery of the two systems of ancient geographical stadia (Alexandrian and Greek) is presented. It is presently used for a rectification of the digitalised maps given in Ptolemy's "Geographike hyphegesis". The stadion definition Eratosthenes has used (1 meridian degree = 700 stadia) was applied also in northern and western Europe and in Asia east of the Tigris river; using it as a scale factor we got very good results for the rectification.
    [Show full text]
  • The International Association of Geodesy 1862 to 1922: from a Regional Project to an International Organization W
    Journal of Geodesy (2005) 78: 558–568 DOI 10.1007/s00190-004-0423-0 The International Association of Geodesy 1862 to 1922: from a regional project to an international organization W. Torge Institut fu¨ r Erdmessung, Universita¨ t Hannover, Schneiderberg 50, 30167 Hannover, Germany; e-mail: [email protected]; Tel: +49-511-762-2794; Fax: +49-511-762-4006 Received: 1 December 2003 / Accepted: 6 October 2004 / Published Online: 25 March 2005 Abstract. Geodesy, by definition, requires international mitment of two scientists from neutral countries, the collaboration on a global scale. An organized coopera- International Latitude Service continued to observe tion started in 1862, and has become today’s Interna- polar motion and to deliver the data to the Berlin tional Association of Geodesy (IAG). The roots of Central Bureau for evaluation. After the First World modern geodesy in the 18th century, with arc measure- War, geodesy became one of the founding members of ments in several parts of the world, and national geodetic the International Union for Geodesy and Geophysics surveys in France and Great Britain, are explained. The (IUGG), and formed one of its Sections (respectively manifold local enterprises in central Europe, which Associations). It has been officially named the Interna- happened in the first half of the 19th century, are tional Association of Geodesy (IAG) since 1932. described in some detail as they prepare the foundation for the following regional project. Simultaneously, Gauss, Bessel and others developed a more sophisticated Key words: Arc measurements – Baeyer – Helmert – definition of the Earth’s figure, which includes the effect Figure of the Earth – History of geodesy – of the gravity field.
    [Show full text]
  • A General Formula for Calculating Meridian Arc Length and Its Application to Coordinate Conversion in the Gauss-Krüger Projection 1
    A General Formula for Calculating Meridian Arc Length and its Application to Coordinate Conversion in the Gauss-Krüger Projection 1 A General Formula for Calculating Meridian Arc Length and its Application to Coordinate Conversion in the Gauss-Krüger Projection Kazushige KAWASE Abstract The meridian arc length from the equator to arbitrary latitude is of utmost importance in map projection, particularly in the Gauss-Krüger projection. In Japan, the previously used formula for the meridian arc length was a power series with respect to the first eccentricity squared of the earth ellipsoid, despite the fact that a more concise expansion using the third flattening of the earth ellipsoid has been derived. One of the reasons that this more concise formula has rarely been recognized in Japan is that its derivation is difficult to understand. This paper shows an explicit derivation of a general formula in the form of a power series with respect to the third flattening of the earth ellipsoid. Since the derived formula is suitable for implementation as a computer program, it has been applied to the calculation of coordinate conversion in the Gauss-Krüger projection for trial. 1. Introduction cannot be expressed explicitly using a combination of As is well known in geodesy, the meridian arc elementary functions. As a rule, we evaluate this elliptic length S on the earth ellipsoid from the equator to integral as an approximate expression by first expanding the geographic latitude is given by the integrand in a binomial series with respect to e2 , regarding e as a small quantity, then readjusting with a a1 e2 trigonometric function, and finally carrying out termwise S d , (1) 2 2 3 2 0 1 e sin integration.
    [Show full text]
  • Article Is Part of the Spe- Senschaft, Kunst Und Technik, 7, 397–424, 1913
    IUGG: from different spheres to a common globe Hist. Geo Space Sci., 10, 151–161, 2019 https://doi.org/10.5194/hgss-10-151-2019 © Author(s) 2019. This work is distributed under the Creative Commons Attribution 4.0 License. The International Association of Geodesy: from an ideal sphere to an irregular body subjected to global change Hermann Drewes1 and József Ádám2 1Technische Universität München, Munich 80333, Germany 2Budapest University of Technology and Economics, Budapest, P.O. Box 91, 1521, Hungary Correspondence: Hermann Drewes ([email protected]) Received: 11 November 2018 – Revised: 21 December 2018 – Accepted: 4 January 2019 – Published: 16 April 2019 Abstract. The history of geodesy can be traced back to Thales of Miletus ( ∼ 600 BC), who developed the concept of geometry, i.e. the measurement of the Earth. Eratosthenes (276–195 BC) recognized the Earth as a sphere and determined its radius. In the 18th century, Isaac Newton postulated an ellipsoidal figure due to the Earth’s rotation, and the French Academy of Sciences organized two expeditions to Lapland and the Viceroyalty of Peru to determine the different curvatures of the Earth at the pole and the Equator. The Prussian General Johann Jacob Baeyer (1794–1885) initiated the international arc measurement to observe the irregular figure of the Earth given by an equipotential surface of the gravity field. This led to the foundation of the International Geodetic Association, which was transferred in 1919 to the Section of Geodesy of the International Union of Geodesy and Geophysics. This paper presents the activities from 1919 to 2019, characterized by a continuous broadening from geometric to gravimetric observations, from exclusive solid Earth parameters to atmospheric and hydrospheric effects, and from static to dynamic models.
    [Show full text]