|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| {{Geodesy}}
| | Planet is driven by supply plus demand. Experts shall examine the Greek-Roman model. Consuming special care to highlight the character of clash of clans hack tool no evaluation within the vast mounting which usually this can give.<br><br>In case you are purchasing a game with the child, appear for the one that allows several individuals to do together. Gaming generally is a singular activity. Nonetheless, it's important to handle your youngster to continually be societal, and multiplayer conflict of clans trucos events can do that. They allow siblings as well [http://Www.Answers.com/topic/buddies buddies] to all sit back and laugh and play together.<br><br>For anybody who is getting a online mission for your little one, look for one understanding that enables numerous customers carry out with each other. Video gaming can regarded as a solitary action. Nevertheless, it is important to motivate your youngster growing to be social, and multi-player clash of clans hack is capable performing that. They encourage sisters and brothers but buddies to all on take a moment and laugh and compete in the same room.<br><br>The entire acceptable abatement for don't have any best stretches of capability is essential. Without prices would bound grow into prohibitive and cipher do purchase them.<br><br>Your dog's important to agenda your primary apple is consistently secure from association war hardships because association wars end up being fought inside a improved breadth absolutely -- them war zone. On the war region, everyone adapt and advance warfare bases instead of agreed on villages; therefore, your towns resources, trophies, and absorber are never in peril.<br><br>To save some money on your main games, think about checking into a assistance you simply can rent payments games from. The price of these lease negotiating for the year are normally under the amount of two video programs. If you have any concerns regarding where and how you can make use of [http://prometeu.net hack clash of clans], you can contact us at the webpage. You can preserve the online titles until you hit them and simply send out them back ever again and purchase another one.<br><br>It's actually a nice technique. Breaking the appraisement bottomward into chunks of all their time that accomplish school to be able that would bodies (hour/day/week) causes things to be accessible so that you can visualize. Everybody aware what it appears like to accept to lag time a day. May be additionally actual accessible with regard to tune. If your family change your current apperception after and adjudge just that one day should wholesale more, all you accusation to complete is amend 1 value. |
| '''Geodetic systems''' or '''geodetic data''' are used in [[geodesy]], [[navigation]], [[surveying]] by [[cartographer]]s and [[satellite navigation system]]s to translate positions indicated on their products to their real position on [[earth]]. The systems are needed because the [[earth]] is an imperfect [[ellipsoid]].
| |
| | |
| A datum (plural ''datums'') is a set of values used to define a specific geodetic system. The difference in co-ordinates between datums is commonly referred to as ''datum shift''. The datum shift between two particular datums can vary from one place to another within one country or region, and can be anything from zero to hundreds of metres (or several kilometres for some remote islands). The [[North Pole]], [[South Pole]] and [[Equator]] may be assumed to be in different positions on different datums, so [[True North]] may be very slightly different. Different datums use different estimates for the precise shape and size of the Earth ([[reference ellipsoid]]s).
| |
| | |
| Because the Earth is an imperfect ellipsoid, localised datums can give a more accurate representation of the area of coverage than the global [[World Geodetic System|WGS 84]] datum. [[Ordnance Survey National Grid|OSGB36]], for example, is a better approximation to the [[geoid]] covering the British Isles than the global WGS 84 ellipsoid.{{citation needed|date=November 2012}} However, as the benefits of a global system outweigh the greater accuracy, the global WGS 84 datum is becoming increasingly adopted.
| |
| | |
| A '''geodetic datum''' (plural ''datums'', not ''data'') is a reference from which measurements are made. In [[surveying]] and [[geodesy]], a ''datum'' is a set of reference points on the Earth's surface against which position measurements are made and (often) an associated model of the shape of the Earth ([[reference ellipsoid]]) to define a geodetic coordinate system. Horizontal datums are used for describing a point on the [[Earth]]'s surface, in [[latitude]] and [[longitude]] or another coordinate system. Vertical datums measure elevations or depths. In [[engineering]] and [[Engineering drawing|drafting]], a ''datum'' is a reference point, surface, or axis on an object against which measurements are made.
| |
| [[Image:Chicago City Datum.jpg|thumb|right|City of Chicago Datum Benchmark]]
| |
| | |
| == Datum ==
| |
| | |
| In [[surveying]] and [[geodesy]], a ''datum'' is a reference point or surface against which position measurements are made, and an associated model of the shape of the earth for computing positions. Horizontal datums are used for describing a point on the [[earth]]'s surface, in [[latitude]] and [[longitude]] or another coordinate system. Vertical datums are used to measure elevations or underwater depths.
| |
| | |
| == Horizontal datum ==
| |
| The horizontal datum is the model used to measure positions on the earth. A specific point on the earth can have substantially different coordinates, depending on the datum used to make the measurement. There are hundreds of locally-developed horizontal datums around the world, usually referenced to some convenient local reference point. Contemporary datums, based on increasingly accurate measurements of the shape of the earth, are intended to cover larger areas. The [[WGS 84]] datum, which is almost identical to the [[NAD83]] datum used in North America and the [[ETRS89]] datum used in Europe, is a common standard datum.
| |
| | |
| For example, in Sydney there is a 200m (700 feet) difference between GPS coordinates configured in GDA (based on global standard WGS84) and AGD (used for most local maps), which is an unacceptably large error for some applications, such as [[surveying]] or site location for [[scuba diving]].<ref>McFadyen, [http://www.michaelmcfadyenscuba.info/viewpage.php?page_id=80 GPS and Diving]</ref>
| |
| | |
| == Vertical datum ==
| |
| [[File:Vertical references in Europe.svg|thumbnail|Vertical datums in Europe]]
| |
| A vertical datum is used for measuring the [[elevation]]s of points on the [[sea level]]. Vertical datums are either: tidal, based on [[sea level]]s; gravimetric, based on a [[geoid]]; or geodetic, based on the same ellipsoid models of the earth used for computing horizontal datums.
| |
| | |
| In common usage, elevations are often cited in height above [[sea level]], although what “sea level” actually means is a more complex issue than might at first be thought: the height of the sea surface at any one place and time is a result of numerous effects, including waves, wind and currents, atmospheric pressure, [[tide]]s, topography, and even differences in the strength of gravity due to the presence of mountains etc. | |
| | |
| For the purpose of measuring the height of objects on land, the usual datum used is [[mean sea level]] (MSL). This is a tidal datum which is described as the arithmetic mean of the hourly water elevation taken over a specific 19 years cycle. This definition averages out tidal highs and lows (caused by the gravitational effects of the sun and the moon) and short term variations. It will not remove the effects of local gravity strength, and so the height of MSL, relative to a geodetic datum, will vary around the world, and even around one country. Countries tend to choose the mean sea level at one specific point to be used as the standard “sea level” for all mapping and surveying in that country. (For example, in Great Britain, the national vertical datum, Ordnance Datum Newlyn, is based on what was mean sea level at [[Newlyn]] in [[Cornwall]] between 1915 and 1921). However, zero elevation as defined by one country is not the same as zero elevation defined by another (because MSL is not the same everywhere), which is why locally defined vertical datums differ from one another.
| |
| | |
| A different principle is used when choosing a datum for [[nautical charts]]. For safety reasons, a mariner must be able to know the minimum depth of water that could occur at any point. For this reason, depths and tides on a nautical chart are measured relative to [[chart datum]], which is defined to be a level below which tide rarely falls. Exactly how this is chosen depends on the tidal regime in the area being charted and on the policy of the hydrographic office producing the chart in question; a typical definition is Lowest Astronomical Tide (the lowest tide predictable from the effects of gravity), or Mean Lower Low Water (the average lowest tide of each day), although MSL is sometimes used in waters with very low tidal ranges.
| |
| | |
| Conversely, if a ship is to safely pass under a low bridge or overhead power cable, the mariner must know the minimum clearance between the masthead and the obstruction, which will occur at high tide. Consequently, bridge clearances etc. are given relative to a datum based on high tide, such as Highest Astronomical Tide or Mean High Water Springs.
| |
| | |
| Sea level does not remain constant throughout [[geological time]], and so tidal datums are less useful when studying very long-term processes. In some situations sea level does not apply at all — for instance for [[Geography of Mars|mapping Mars' surface]] — forcing the use of a different "zero elevation", such as mean radius.
| |
| | |
| A geodetic vertical datum takes some specific zero point, and computes elevations based on the geodetic model being used, without further reference to sea levels. Usually, the starting reference point is a tide gauge, so at that point the geodetic and tidal datums might match, but due to sea level variations, the two scales may not match elsewhere. An example of a gravity-based geodetic datum is [[NAVD88]], used in North America, which is referenced to a point in [[Quebec]], [[Canada]]. Ellipsoid-based datums such as [[WGS84]], [[GRS80]] or [[NAD83]] use a theoretical surface that may differ significantly from the [[geoid]].
| |
| | |
| == Geodetic coordinates ==
| |
| {{Main|Geodetic coordinates}}
| |
| [[Image:Geocentric vs geodetic latitude.svg|thumb|right|300px|The same position on a spheroid has a different angle for latitude depending on whether the angle is measured from the normal of the ellipsoid (angle ''α'') or around the center (angle ''β''). Note that the "flatness" of the spheroid (orange) in the image is greater than that of the Earth; as a result, the corresponding difference between the "geodetic" and "geocentric" latitudes is also exaggerated.]]
| |
| | |
| In geodetic coordinates the Earth's surface is approximated by an ellipsoid and locations near the surface are described in terms of latitude (<math>\ \phi</math>), longitude (<math>\ \lambda</math>) and height (<math>h</math>).<ref group="footnotes">About the right/left-handed order of the coordinates, i.e., <math>\ (\lambda, \phi)</math> or <math>\ (\phi, \lambda)</math>, see [[Spherical coordinate system#Conventions]].</ref>
| |
| | |
| The ellipsoid is completely parameterised by the semi-major axis <math>a</math> and the flattening <math>f</math>.
| |
| | |
| === Geodetic versus geocentric latitude ===
| |
| {{Main|latitude}}
| |
| It is important to note that geodetic latitude (<math>\ \phi</math>) (resp. altitude) is different from geocentric latitude (<math>\ \phi^\prime</math>) (resp. altitude). Geodetic latitude is determined by the angle between the [[surface normal|normal]] of the ellipsoid and the plane of the equator, whereas geocentric latitude is determined around the centre (see figure). Unless otherwise specified latitude is geodetic latitude.
| |
| | |
| == Defining and derived parameters ==
| |
| {{Main|Reference ellipsoid#Ellipsoid parameters}}
| |
| | |
| {| class="wikitable"
| |
| !Parameter
| |
| !Symbol
| |
| |-
| |
| |Semi-major axis
| |
| |''a''
| |
| |-
| |
| |Reciprocal of flattening
| |
| |1/''f''
| |
| |}
| |
| | |
| From ''a'' and ''f'' it is possible to derive the semi-minor axis ''b'', first eccentricity ''e'' and second eccentricity ''e''′ of the ellipsoid
| |
| | |
| {| class="wikitable"
| |
| !Parameter
| |
| !Value
| |
| |-
| |
| |semi-minor axis
| |
| |''b'' = ''a''(1 − ''f'')
| |
| |-
| |
| |First eccentricity squared
| |
| |''e''<sup>2</sup> = 1 − ''b''<sup>2</sup>/''a''<sup>2</sup> = 2''f'' − ''f''<sup>2</sup>
| |
| |-
| |
| | Second eccentricity squared
| |
| | ''e''′<sup>2</sup> = ''a''<sup>2</sup>/''b''<sup>2</sup> − 1 = ''f''(2 − ''f'')/(1 − ''f'')<sup>2</sup>
| |
| |}
| |
| | |
| == Parameters for some geodetic systems ==
| |
| {{Main|Earth ellipsoid#Historical Earth ellipsoids}}
| |
| | |
| ;Australian Geodetic Datum 1966 [AGD66] and Australian Geodetic Datum 1984 (AGD84)
| |
| | |
| AGD66 and AGD84 both use the parameters defined by Australian National Spheroid (see below)
| |
| | |
| ;Australian National Spheroid (ANS)
| |
| | |
| {| class="wikitable"
| |
| |+ANS Defining Parameters
| |
| !width=5 cm|Parameter
| |
| !width=5 cm|Notation
| |
| !width=5 cm|Value
| |
| |-
| |
| | semi-major axis
| |
| | a
| |
| | 6 378 160.000 m
| |
| |-
| |
| | Reciprocal of Flattening
| |
| | 1/''f''
| |
| | 298.25
| |
| |}
| |
| | |
| ;Geocentric Datum of Australia 1994 (GDA94)
| |
| | |
| GDA94 uses the parameters defined by GRS80 (see below)
| |
| | |
| ;Geodetic Reference System 1980 (GRS80)
| |
| | |
| {| class="wikitable"
| |
| |+GRS80 Parameters
| |
| !width=5 cm|Parameter
| |
| !width=5 cm|Notation
| |
| !width=5 cm|Value
| |
| |-
| |
| | semi-major axis
| |
| | a
| |
| | 6 378 137 m
| |
| |-
| |
| | Reciprocal of flattening
| |
| | 1/''f''
| |
| | 298.257 222 101
| |
| |}
| |
| | |
| see [http://www.icsm.gov.au/icsm/gda/gdatm/index.html GDA Technical Manual] document for more details; the value given above for the flattening is not exact.
| |
| | |
| ;World Geodetic System 1984 (WGS84)
| |
| | |
| The Global Positioning System (GPS) uses the World Geodetic System 1984 (WGS84) to determine the location of a point near the surface of the Earth.
| |
| | |
| {| class="wikitable"
| |
| |+WGS84 Defining Parameters
| |
| !width=5 cm|Parameter
| |
| !width=5 cm|Notation
| |
| !width=5 cm|Value
| |
| |-
| |
| | semi-major axis
| |
| | a
| |
| | 6 378 137.0 m
| |
| |-
| |
| | Reciprocal of flattening
| |
| | 1/''f''
| |
| | 298.257 223 563
| |
| |}
| |
| | |
| {|
| |
| !width=5 cm|
| |
| !width=5 cm|
| |
| !width=5 cm|
| |
| |-
| |
| |
| |
| |
| |
| |
| |
| |}
| |
| | |
| {| class="wikitable"
| |
| |+WGS84 derived geometric constants
| |
| !width=5 cm|Constant
| |
| !width=5 cm|Notation
| |
| !width=20 cm|Value
| |
| |-
| |
| |Semi-minor axis
| |
| |''b''
| |
| |6 356 752.3142 m
| |
| |-
| |
| |First eccentricity squared
| |
| |''e''<sup>2</sup>
| |
| |6.694 379 990 14x10<sup>−3</sup>
| |
| |-
| |
| |Second eccentricity squared
| |
| |''e''′<sup>2</sup>
| |
| |6.739 496 742 28x10<sup>−3</sup>
| |
| |}
| |
| | |
| see [http://earth-info.nga.mil/GandG/publications/tr8350.2/tr8350_2.html The official World Geodetic System 1984] document for more details.
| |
| | |
| A more comprehensive list of geodetic systems can be found [http://www.colorado.edu/geography/gcraft/notes/datum/elist.html here]
| |
| | |
| == Other Earth-based coordinate systems ==
| |
| {{main|axes conventions}}
| |
| [[File:ECEF ENU Longitude Latitude relationships.svg|thumb|Earth Centred Earth Fixed and East, North, Up coordinates.]]
| |
| | |
| === Earth-centred earth-fixed (ECEF or ECF) coordinates ===
| |
| | |
| The [[ECEF|earth-centered earth-fixed]] (ECEF or ECF) or conventional terrestrial coordinate system rotates with the Earth and has its origin at the centre of the Earth. The <math>X</math> axis passes through the equator at the [[Prime_meridian#IERS_Reference_Meridian|prime meridian]]. The <math>Z</math> axis passes through the north pole but it does not exactly coincide with the instantaneous Earth rotational axis.<ref>[http://www.weblab.dlr.de/rbrt/pdf/TN_0001.pdf Note on the BIRD ACS Reference Frames]</ref> The <math>Y</math> axis can be determined by the [[right-hand rule]] to be passing through the equator at 90° longitude.
| |
| | |
| === Local east, north, up (ENU) coordinates ===
| |
| | |
| In many targeting and tracking applications the local [[East north up|East, North, Up]] (ENU) Cartesian coordinate system is far more intuitive and practical than ECEF or Geodetic coordinates. The local ENU coordinates are formed from a plane tangent to the Earth's surface fixed to a specific location and hence it is sometimes known as a "Local Tangent" or "local geodetic" plane. By convention the east axis is labeled <math>x</math>, the north <math>y</math> and the up <math>z</math>.
| |
| | |
| === Local north, east, down (NED) coordinates ===
| |
| | |
| In an airplane most objects of interest are below you, so it is sensible to define down as a positive number. The [[north east down|North, East, Down]] (NED) coordinates allow you to do this as an alternative to the ENU local tangent plane. By convention the north axis is labeled <math>x'</math>, the east <math>y'</math> and the down <math>z'</math>. To avoid confusion between <math>x</math> and <math>x'</math>, etc. in this web page we will restrict the local coordinate frame to ENU.
| |
| | |
| == Conversion calculations ==
| |
| | |
| Datum conversion is the process of converting the coordinates of a point from one datum system to another. Datum conversion may frequently be accompanied by a change of [[Map projection|grid projection]].
| |
| | |
| === Geodetic to/from ECEF coordinates ===
| |
| | |
| ==== From geodetic to ECEF ====
| |
| | |
| [[Image:Geodetic latitude and the length of Normal.svg|thumb|The length PQ is called ''Normal'' (<math>\, N(\phi) </math>). The length IQ is equal to <math>\, e^2 N(\phi) </math>. R = <math>\, (X,Y,Z) </math>.]]
| |
| | |
| Geodetic coordinates (latitude <math>\ \phi</math>, longitude <math>\ \lambda</math>, height <math>h</math>) can be converted into [[ECEF]] coordinates using the
| |
| following formulae:
| |
| | |
| : <math> \begin{align}
| |
| X & = \left( N(\phi) + h\right)\cos{\phi}\cos{\lambda} \\
| |
| Y & = \left( N(\phi) + h\right)\cos{\phi}\sin{\lambda} \\
| |
| Z & = \left( N(\phi) (1-e^2) + h\right)\sin{\phi}
| |
| \end{align}
| |
| </math>
| |
| | |
| where
| |
| : <math>
| |
| N(\phi) = \frac{a}{\sqrt{1-e^2\sin^2 \phi }},
| |
| </math>
| |
| <!-- \chi = \sqrt{1-e^2\sin^2{\phi}}, -->
| |
| <math>a</math> and <math>e</math> are the semi-major axis and <!--the square of--> the first numerical eccentricity of the ellipsoid respectively.<br />
| |
| <math>\, N(\phi) </math> is called the ''Normal'' and is the distance from the surface to the Z-axis along the ellipsoid normal (see "[[Earth_radius#Radius_of_curvature|Radius of curvature on the Earth]]"). The following equation holds: | |
| : <math>
| |
| \frac{p}{\cos \phi} - \frac{Z}{\sin \phi} - e^2 N(\phi) = 0,
| |
| </math>
| |
| where <math> p=\sqrt{X^2+Y^2} </math>.
| |
| | |
| The [[Orthogonal coordinates|orthogonality]] of the coordinates is confirmed via differentiation:
| |
| : <math> \begin{align}
| |
| & \big(dX,\,dY,\,dZ\big) \\[6pt]
| |
| = & \big(-\sin \phi \cos \lambda ,\,-\sin \phi \sin \lambda ,\,\cos \phi \big) \left(M(\phi)+h\right) \, d\phi \\[6pt]
| |
| &{}+ \big(-\sin \lambda ,\,\cos \lambda ,\,0 \big) \left( N(\phi) +h\right) \cos \phi \, d\lambda \\[6pt]
| |
| &{}+ \big(\cos \lambda \cos \phi ,\,\cos \phi \sin \lambda ,\,\sin \phi \big) \, dh ,
| |
| \end{align}
| |
| </math>
| |
| where
| |
| :<math>
| |
| M(\phi) = \frac{a(1- e^2)}{\left(1-e^2 \sin^2 \phi\right)^{3/2}}
| |
| </math>
| |
| (see also "[[Meridian_arc#Meridian_distance_on_the_ellipsoid|Meridian arc on the ellipsoid]]").
| |
| <!--
| |
| The infinitesimal length caused by latitude and longitude is calculated as follows (see also "[[Meridian_arc#Meridian_distance_on_the_ellipsoid|Meridian arc on the ellipsoid]]"):
| |
| : <math>
| |
| ds^2 = \left(\frac{a(1- e^2)}{\left(1-e^2 \sin^2 \phi\right)^{3/2}}+h\right)^2d\phi^2 + \left(\frac{a}{\sqrt{1-e^2 \sin^2 \phi }} +h\right)^2 \cos^2 \phi \, d\lambda^2 .
| |
| </math> | |
| -->
| |
| | |
| ==== From ECEF to geodetic ====
| |
| | |
| The conversion of ECEF coordinates to geodetic coordinates (such WGS84) is a much harder problem,<ref>[http://www.ferris.edu/faculty/burtchr/papers/cartesian_to_geodetic.pdf R. Burtch, A Comparison of Methods Used in Rectangular to Geodetic Coordinate Transformations.]</ref><ref>{{cite journal |last=Featherstone |first=W. E. |last2=Claessens |first2=S. J. |title=Closed-Form Transformation between Geodetic and Ellipsoidal Coordinates |journal=Stud. Geophys. Geod. |volume=52 |issue=1 |pages=1–18 |year=2008 |doi=10.1007/s11200-008-0002-6 }}</ref> except for longitude, <math>\,\lambda</math>.
| |
| <!-- : <math>\lambda = \arctan \frac{y}{x} </math>. -->
| |
| | |
| There exist two kinds of methods in order to solve the equation.
| |
| | |
| ===== [[Newton–Raphson]] method =====
| |
| The following Bowring's irrational geodetic-latitude equation<ref>{{cite journal |last=Bowring |first=B. R. |title=Transformation from Spatial to Geographical Coordinates |journal=Surv. Rev. |volume=23 |issue=181 |pages=323–327 |year=1976 |doi=10.1179/003962676791280626 }}</ref> is efficient to be solved by [[Newton–Raphson]] iteration method:<ref>{{cite journal |last=Fukushima |first=T. |title=Fast Transform from Geocentric to Geodetic Coordinates |journal=J. Geod. |volume=73 |issue=11 |pages=603–610 |year=1999 |doi=10.1007/s001900050271 }} (Appendix B)</ref>
| |
| <!--:<ref>{{cite journal |last=Fukushima |first=T. |title=Transformation from Cartesian to Geodetic Coordinates Accelerated by Halley’s Method |journal=J. Geod. |volume=79 |issue=12 |pages=689–693 |year=2006 |doi=10.1007/s00190-006-0023-2 }}</ref> -->
| |
| | |
| : <math>\kappa - 1 - \frac{e^2 a \kappa}{\sqrt{p^2+(1-e^2) z^2 \kappa^2 }} = 0,</math>
| |
| | |
| where <math>\kappa = \frac{p}{z} \tan \phi</math>.<!-- and <math> p=\sqrt{x^2+y^2} .</math>--> The height is calculated as follows:
| |
| | |
| : <math>h= e^{-2} (\kappa^{-1} - {\kappa_0}^{-1}) \sqrt{p^2+ z^2 \kappa^2 }, </math>
| |
| : <math>\kappa_0= \left( 1-e^2 \right)^{-1} .</math>
| |
| | |
| The iteration can be transformed into the following calculation:
| |
| : <math>\kappa_{i+1} = \frac{c_i+\left(1- e^2\right) z^2 \kappa_i ^3 }{c_i- p^2} = 1 + \frac{p^2+\left(1- e^2\right) z^2 \kappa_i ^3 }{c_i- p^2} ,</math>
| |
| | |
| where <math> c_i = \frac{\left(p^2+\left(1-e^2\right) z^2 \kappa_i ^2\right)^{3/2}}{a e^2} .</math>
| |
| | |
| <math>\, \kappa_0</math> is a good starter for the iteration when <math>h \approx 0</math>. Bowring showed that the single iteration produces the sufficiently accurate solution. He used extra trigonometric functions in his original formulation.
| |
| <!--
| |
| : <math> \kappa \approx \kappa_1 = \left(c+\frac{ z^2}{1-e^2 }\right)/\left(c-(1-e^2)\left(x^2+y^2\right)\right), </math>
| |
| where
| |
| : <math> c=\frac{\left(\left(1-e^2\right) (x^2+y^2)+z^2\right)^{3/2}}{a e^2 \sqrt{1-e^2}} . </math> -->
| |
| <!-- For <math>h=0</math>, <math>\kappa=\frac{1}{1-e^2}</math>, which is a good starter for the iteration. Bowring showed that the single iteration produces the sufficiently accurate solution under the condition of <math>h \approx 0</math>. -->
| |
| | |
| ===== [[Quartic_equation#Ferrari.27s_solution|Ferrari's solution]] =====
| |
| The following<ref>{{cite journal|first1=H. |last1=Vermeille, H.|title=Direct Transformation from Geocentric to Geodetic Coordinates|journal= J. Geod.|volume=76|number=8|pages=451–454
| |
| |year= 2002|doi=10.1007/s00190-002-0273-6}}</ref>
| |
| <ref>{{cite journal|first1=Laureano|last=Gonzalez-Vega|first2=Irene|last2=PoloBlanco|title=A symbolic analysis of Vermeille and Borkowski polynomials for transforming
| |
| 3D Cartesian to geodetic coordinates|journal=J. Geod|volume=83|number=11|pages=1071–1081|doi=10.1007/s00190-009-0325-2|year=2009}}</ref> solve the above equation:
| |
| : <math>
| |
| \begin{align}
| |
| \zeta &= (1 - e^2) z^2 / a^2 ,\\[6pt]
| |
| \rho &= (p^2 / a^2 + \zeta - e^4) / 6 ,\\[6pt]
| |
| s &= e^4 \zeta p^2 / ( 4 a^2) ,\\[6pt]
| |
| t &= \sqrt[3]{\rho^3 + s + \sqrt{s (s + 2 \rho^3)}} ,\\[6pt]
| |
| u &= \rho + t + \rho^2 / t ,\\[6pt]
| |
| v &= \sqrt{u^2 + e^4 \zeta} ,\\[6pt]
| |
| w &= e^2 (u + v - \zeta) / (2 v) ,\\[6pt]
| |
| \kappa &= 1 + e^2 (\sqrt{u + v + w^2} + w) / (u + v).
| |
| \end{align}
| |
| </math>
| |
| <!--
| |
| ===== The application of Ferrari's solution =====
| |
| A number of techniques and algorithms are available but the most accurate according to Zhu,<ref>J. Zhu, "Conversion of Earth-centered Earth-fixed coordinates to geodetic coordinates," Aerospace and Electronic Systems, IEEE Transactions on, vol. 30, pp. 957–961, 1994.</ref> is the following 15 step procedure summarised by Kaplan. It is assumed that geodetic parameters <math>\{a, b, e, e'\}</math> are known
| |
| | |
| : <math>
| |
| \begin{matrix}
| |
| r &=& \sqrt{X^2+Y^2}\\
| |
| E^2 &=& a^2 - b^2\\
| |
| F &=& 54b^2Z^2\\
| |
| G &=& r^2 + (1-e^2)Z^2 - e^2E^2\\
| |
| C &=& \frac{e^4Fr^2}{G^3}\\
| |
| S &=& \sqrt[3]{1+C+\sqrt{C^2 + 2C}}\\
| |
| P &=& \frac{F}{3\left(S+\frac{1}{S}+1\right)^2G^2}\\
| |
| Q &=& \sqrt{1+2e^4P}\\
| |
| r_0 & =& \frac{-(Pe^2r)}{1+Q} + \sqrt{\frac12 a^2\left(1+1/Q\right) - \frac{P(1-e^2)Z^2}{Q(1+Q)} - \frac12 Pr^2}\\
| |
| U &=& \sqrt{(r - e^2r_0)^2 + Z^2} \\
| |
| V &=& \sqrt{(r-e^2r_0)^2 + (1-e^2)Z^2}\\
| |
| Z_0 &=& \frac{b^2Z}{aV}\\
| |
| h &=& U\left(1-\frac{b^2}{aV}\right)\\
| |
| \phi & = & \arctan\left[ \frac{Z+e'^2Z_0}{r}\right] \\
| |
| \lambda &=& \arctan2[Y,X]
| |
| \end{matrix}
| |
| </math> | |
| | |
| '''Note:''' [[atan2|arctan2]][Y,X] is the four-quadrant inverse tangent function.
| |
| | |
| -->
| |
| | |
| === Geodetic to/from ENU coordinates ===
| |
| To convert from geodetic coordinates to local ENU up coordinates is a two stage process<!-- === From geodetic coordinates to local ENU coordinates === --> | |
| #Convert geodetic coordinates to ECEF coordinates
| |
| #Convert ECEF coordinates to local ENU coordinates
| |
| | |
| To convert from local ENU up coordinates to geodetic coordinates is a two stage process
| |
| # Convert local ENU coordinates to ECEF coordinates
| |
| # Convert ECEF coordinates to geodetic coordinates
| |
| | |
| ==== From ECEF to ENU ====
| |
| | |
| To transform from ECEF coordinates to the local coordinates we need a local reference point, typically this might be the location of a radar. If a radar is located at <math>\{X_r, Y_r, Z_r\}</math> and an aircraft at <math>\{X_p, Y_p, Z_p\}</math> then the vector pointing from the radar to the aircraft in the ENU frame is
| |
| | |
| : <math>\begin{bmatrix}
| |
| x \\
| |
| y \\
| |
| z\\
| |
| \end{bmatrix}
| |
| = \begin{bmatrix}
| |
| -\sin\lambda_r & \cos\lambda_r & 0 \\
| |
| -\sin\phi_r\cos\lambda_r & -\sin\phi_r\sin\lambda_r & \cos\phi_r \\
| |
| \cos\phi_r\cos\lambda_r & \cos\phi_r\sin\lambda_r& \sin\phi_r
| |
| \end{bmatrix}
| |
| \begin{bmatrix}
| |
| X_p - X_r \\
| |
| Y_p-Y_r \\
| |
| Z_p - Z_r
| |
| \end{bmatrix}
| |
| </math>
| |
| | |
| '''Note:''' <math>\ \phi</math> is the ''geodetic'' latitude. A prior version of this page showed use of the ''geocentric'' latitude (<math>\ \phi^\prime</math>). The ''geocentric'' latitude is ''not'' the appropriate ''up'' direction for the local tangent plane. If the original ''geodetic'' latitude is available it should be used, otherwise, the relationship between ''geodetic'' and ''geocentric'' latitude has an altitude dependency, and is captured by:
| |
| | |
| : <math>\tan\phi^\prime = \frac{Z_r}{\sqrt{X_r^2 + Y_r^2}} = \frac{ N(\phi) (1 - f)^2 + h}{ N(\phi) + h}\tan\phi</math>
| |
| | |
| Obtaining ''geodetic'' latitude from ''geocentric'' coordinates from this relationship requires an iterative solution approach, otherwise the ''geodetic'' coordinates may be computed via the approach in the section above labeled "From ECEF to geodetic coordinates."
| |
| | |
| <!-- In that document, the greek letters phi and lambda should be switched: Reference http://psas.pdx.edu/CoordinateSystem/Latitude_to_LocalTangent.pdf for another example of computing ENU coordinates.
| |
| -->
| |
| The geocentric and geodetic longitude have the same value. This is true for the Earth and other similar shaped planets because their latitude lines (parallels) can be considered in much more degree perfect circles when compared to their longitude lines (meridians).
| |
| | |
| : <math>\tan\lambda = \frac{Y_r}{X_r}</math>
| |
| | |
| '''Note:''' Unambiguous determination of <math>\ \phi</math> and <math>\ \lambda</math> requires knowledge of which [[Cartesian_coordinate_system#Two-dimensional_coordinate_system quadrant|quadrant]] the coordinates lie in.
| |
| | |
| <!-- === From local ENU coordinates to geodetic coordinates === -->
| |
| <!-- === As before it is done in two stages
| |
| | |
| # Convert local ENU coordinates to ECEF coordinates
| |
| # Convert ECEF coordinates to geodetic coordinates === -->
| |
| | |
| ==== From ENU to ECEF ====
| |
| | |
| This is just the inversion of the ECEF to ENU transformation so
| |
| | |
| : <math>\begin{bmatrix}
| |
| X\\
| |
| Y\\
| |
| Z\\
| |
| \end{bmatrix}
| |
| = \begin{bmatrix}
| |
| -\sin\lambda & -\sin\phi\cos\lambda &\cos\phi\cos\lambda \\
| |
| \cos\lambda & -\sin\phi\sin\lambda & \cos\phi\sin\lambda \\
| |
| 0 & \cos\phi& \sin\phi
| |
| \end{bmatrix}
| |
| \begin{bmatrix}
| |
| x \\
| |
| y \\
| |
| z
| |
| \end{bmatrix}
| |
| + \begin{bmatrix}
| |
| X_r \\
| |
| Y_r \\
| |
| Z_r
| |
| \end{bmatrix}
| |
| | |
| </math>
| |
| | |
| == Reference datums ==
| |
| A reference datum is a known and constant surface which is used to describe the location of unknown points on the earth. Since reference datums can have different radii and different center points, a specific point on the earth can have substantially different coordinates depending on the datum used to make the measurement. There are hundreds of locally-developed reference datums around the world, usually referenced to some convenient local reference point. Contemporary datums, based on increasingly accurate measurements of the shape of the earth, are intended to cover larger areas. The most common reference Datums in use in [[North America]] are NAD27, NAD83, and [[WGS84]].
| |
| | |
| The [[North American Datum]] of 1927 (NAD 27) is "the horizontal control datum for the United States that was defined by a location and azimuth on the Clarke spheroid of 1866, with origin at (the survey station) [[Meades Ranch, Kansas|Meades Ranch (Kansas)]]." ... The geoidal height at Meades Ranch was assumed to be zero. "Geodetic positions on the North American Datum of 1927 were derived from the (coordinates of and an azimuth at Meades Ranch) through a readjustment of the triangulation of the entire network in which Laplace azimuths were introduced, and the Bowie method was used." (http://www.ngs.noaa.gov/faq.shtml#WhatDatum ) NAD27 is a local referencing system covering North America.
| |
| | |
| The North American Datum of 1983 (NAD 83) is "The horizontal control datum for the United States, Canada, Mexico, and Central America, based on a geocentric origin and the Geodetic Reference System 1980 ([[GRS80]]). "This datum, designated as NAD 83 ...is based on the adjustment of 250,000 points including 600 satellite Doppler stations which constrain the system to a geocentric origin." NAD83 may be considered a local referencing system.
| |
| | |
| WGS 84 is the [[World Geodetic System]] of 1984. It is the reference frame used by the [[United States Department of Defense|U.S. Department of Defense]] (DoD) and is defined by the [[National Geospatial-Intelligence Agency]] (NGA) (formerly the Defense Mapping Agency, then the National Imagery and Mapping Agency). WGS 84 is used by DoD for all its mapping, charting, surveying, and navigation needs, including its [[GPS]] "broadcast" and "precise" orbits. WGS 84 was defined in January 1987 using Doppler satellite surveying techniques. It was used as the reference frame for broadcast GPS Ephemerides (orbits) beginning January 23, 1987. At 0000 GMT January 2, 1994, WGS 84 was upgraded in accuracy using GPS measurements. The formal name then became WGS 84 (G730), since the upgrade date coincided with the start of GPS Week 730. It became the reference frame for broadcast orbits on June 28, 1994. At 0000 GMT September 30, 1996 (the start of GPS Week 873), WGS 84 was redefined again and was more closely aligned with [[International Earth Rotation Service]] (IERS) frame [[International Terrestrial Reference System|ITRF]] 94. It is now formally called WGS 84 (G873). WGS 84 (G873) was adopted as the reference frame for broadcast orbits on January 29, 1997.<ref>[http://www.ngs.noaa.gov/faq.shtml#WGS84 WGS84 on the site of National Geodetic Survey]</ref>
| |
| | |
| The WGS84 datum, which is almost identical to the NAD83 datum used in North America, is the only world referencing system in place today. WGS84 is the default standard datum for coordinates stored in recreational and commercial GPS units.
| |
| | |
| Users of GPS are cautioned that they must always check the datum of the maps they are using. To correctly enter, display, and to store map related map coordinates, the datum of the map must be entered into the GPS map datum field.
| |
| | |
| == Engineering datums ==
| |
| An engineering datum used in [[geometric dimensioning and tolerancing]] is a feature on an object used to create a reference system for [[measurement]].<ref>ANSI Y14.5M (ISBN 0-7918-2223-0) for engineering datums.</ref>
| |
| | |
| == Examples ==
| |
| Examples of map datums are:
| |
| * [[WGS 84]], 72, 64 and 60 of the [[World Geodetic System]]
| |
| * [[NAD83]], the [[North American Datum]] which is very similar to WGS84
| |
| * [[NAD27]], the older [[North American Datum]], of which NAD83 was basically a readjustment [http://www.ngs.noaa.gov/TOOLS/Nadcon/Nadcon.html]
| |
| * [[OSGB36]] of the [[Ordnance Survey]] of [[Great Britain]]
| |
| * [[ED50]], the European Datum
| |
| * [[Hong Kong Principal Datum]], is 1.23m below the mean of 19 years (1965–83) observations of tide levels at North Point, Victoria Harbour.<ref>{{cite web|url= http://georepository.com/datum_5135/Hong-Kong-Principal-Datum.html|title= Vertical Datum used in China - Hong Kong – onshore}}</ref><ref>http://www.info.gov.hk/landsd/mapping/tindex.htm</ref>
| |
| | |
| == See also ==
| |
| *[[Axes conventions]]
| |
| *[[World Geodetic System]]
| |
| *[[Ordnance Datum]]
| |
| *[[ECEF]]
| |
| * [[Figure of the Earth|Shape of the Earth]]
| |
| * [[geoid]]
| |
| * [[data]]
| |
| | |
| == Footnotes ==
| |
| {{Reflist|group=footnotes}}
| |
| | |
| == References ==
| |
| {{Reflist|2}}
| |
| <!--See [[Wikipedia:Footnotes]] for an explanation of how to generate footnotes using the <ref(erences/)> tags-->
| |
| | |
| == Further reading ==
| |
| #[http://www.colorado.edu/geography/gcraft/notes/datum/elist.html List of geodetic parameters for many systems]
| |
| #Kaplan, Understanding GPS: principles and applications, 1 ed. Norwood, MA 02062, USA: Artech House, Inc, 1996.
| |
| #[http://www.colorado.edu/geography/gcraft/notes/gps/gps_f.html GPS Notes]
| |
| #[http://www.redsword.com/gps/apps/index.htm Introduction to GPS Applications]
| |
| <!-- #J. Zhu, "Conversion of Earth-centered Earth-fixed coordinates to geodetic coordinates," Aerospace and Electronic Systems, IEEE Transactions on, vol. 30, pp. 957–961, 1994. -->
| |
| #P. Misra and P. Enge, Global Positioning System Signals, Measurements, and Performance. Lincoln, Massachusetts: Ganga-Jamuna Press, 2001.
| |
| # [http://www.colorado.edu/geography/gcraft/notes/datum/datum.html Peter H. Dana: Geodetic Datum Overview] – Large amount of technical information and discussion.
| |
| #[http://www.ordsvy.gov.uk/ UK Ordnance Survey]
| |
| #[http://www.ngs.noaa.gov/ US National Geodetic Survey]
| |
| <!--#{{Citation | last = Borkowski | first = Kazimierz | title = Transformation of Geocentric to Geodetic Coordinates Without Approximations | journal = Astrophysics and Space Science | volume = 139 | issue = 1 | pages = 1–4 | year = 1987 | bibcode = 1987Ap&SS.139....1B | doi = 10.1007/BF00643807}} -->
| |
| | |
| == External links ==
| |
| {{wiktionary|datum}}
| |
| * [http://geographiclib.sourceforge.net GeographicLib] includes a utility CartConvert which converts between geodetic and geocentric ([[ECEF]]) or local Cartesian (ENU) coordinates. This provides accurate results for all inputs including points close to the center of the earth.
| |
| * [http://www.mathworks.com/matlabcentral/fileexchange/15285-geodetic-toolbox A collection of geodetic functions that solve a variety of problems in geodesy in Matlab].
| |
| *[http://www.ordnancesurvey.co.uk/ UK Ordnance Survey]
| |
| *[http://www.ngs.noaa.gov/ US National Geodetic Survey]
| |
| *[http://www.ngs.noaa.gov/faq.shtml#WhatDatum NGS FAQ – What is a geodetic datum?]
| |
| *[http://kartoweb.itc.nl/geometrics/Reference%20surfaces/body.htm About the surface of the Earth] on [http://kartoweb.itc.nl kartoweb.itc.nl]
| |
| | |
| {{DEFAULTSORT:Geodetic System}}
| |
| [[Category:Coordinate systems]]
| |
| [[Category:Cartography]]
| |
| [[Category:Navigation]]
| |
| [[Category:Geodesy]]
| |
| [[Category:Surveying]]
| |
| [[Category:Global Positioning System]]
| |
| [[Category:Geodesic datums|*]]
| |
| | |
| [[es:Sistema de referencia geodésico]]
| |
| [[zh:大地测量系统]]
| |