Geographic Coordinate Conversion

In geodesy, conversion among different geographic coordinate systems is made necessary by the different geographic coordinate systems in use across the world and over time.

Coordinate conversion is composed of a number of different types of conversion: format change of geographic coordinates, conversion of coordinate systems, or transformation to different geodetic datums. Geographic coordinate conversion has applications in cartography, surveying, navigation and geographic information systems.

In geodesy, geographic coordinate conversion is defined as translation among different coordinate formats or map projections all referenced to the same geodetic datum. A geographic coordinate transformation is a translation among different geodetic datums. Both geographic coordinate conversion and transformation will be considered in this article.

This article assumes readers are already familiar with the content in the articles geographic coordinate system and geodetic datum.

Change of units and format

Informally, specifying a geographic location usually means giving the location's latitude and longitude. The numerical values for latitude and longitude can occur in a number of different units or formats:

There are 60 minutes in a degree and 60 seconds in a minute. Therefore, to convert from a degrees minutes seconds format to a decimal degrees format, one may use the formula

    Geographic Coordinate Conversion .

To convert back from decimal degree format to degrees minutes seconds format,

    Geographic Coordinate Conversion 

where Geographic Coordinate Conversion  and Geographic Coordinate Conversion  are just temporary variables to handle both positive and negative values properly.

Coordinate system conversion

A coordinate system conversion is a conversion from one coordinate system to another, with both coordinate systems based on the same geodetic datum. Common conversion tasks include conversion between geodetic and earth-centered, earth-fixed (ECEF) coordinates and conversion from one type of map projection to another.

From geodetic to ECEF coordinates

Geographic Coordinate Conversion 
The length PQ, called the prime vertical radius, is Geographic Coordinate Conversion . The length IQ is equal to Geographic Coordinate Conversion . Geographic Coordinate Conversion .

Geodetic coordinates (latitude Geographic Coordinate Conversion , longitude Geographic Coordinate Conversion , height Geographic Coordinate Conversion ) can be converted into ECEF coordinates using the following equation:

    Geographic Coordinate Conversion 

where

    Geographic Coordinate Conversion 

and Geographic Coordinate Conversion  and Geographic Coordinate Conversion  are the equatorial radius (semi-major axis) and the polar radius (semi-minor axis), respectively. Geographic Coordinate Conversion  is the square of the first numerical eccentricity of the ellipsoid. Geographic Coordinate Conversion  is the flattening of the ellipsoid. The prime vertical radius of curvature Geographic Coordinate Conversion  is the distance from the surface to the Z-axis along the ellipsoid normal.

Properties

The following condition holds for the longitude in the same way as in the geocentric coordinates system:

    Geographic Coordinate Conversion 

And the following holds for the latitude:

    Geographic Coordinate Conversion 

where Geographic Coordinate Conversion , as the parameter Geographic Coordinate Conversion  is eliminated by subtracting

    Geographic Coordinate Conversion 

and

    Geographic Coordinate Conversion 

Furthermore:

    Geographic Coordinate Conversion 

Orthogonality

The orthogonality of the coordinates is confirmed via differentiation:

    Geographic Coordinate Conversion 

where

    Geographic Coordinate Conversion 

(see also "Meridian arc on the ellipsoid").

From ECEF to geodetic coordinates

The conversion of ECEF coordinates to longitude is:

    Geographic Coordinate Conversion .

where atan2 is the quadrant-resolving arc-tangent function. The geocentric longitude and geodetic longitude have the same value; this is true for Earth and other similar shaped planets because they have a large amount of rotational symmetry around their spin axis (see triaxial ellipsoidal longitude for a generalization).

The conversion for the latitude and height involves a circular relationship involving N, which is a function of latitude:

    Geographic Coordinate Conversion ,
    Geographic Coordinate Conversion .

It can be solved iteratively, for example, starting with a first guess h≈0 then updating N. More elaborate methods are shown below. The procedure is, however, sensitive to small accuracy due to Geographic Coordinate Conversion  and Geographic Coordinate Conversion  being maybe 106 apart.

Newton–Raphson method

The following Bowring's irrational geodetic-latitude equation, derived simply from the above properties, is efficient to be solved by Newton–Raphson iteration method:

    Geographic Coordinate Conversion 

where Geographic Coordinate Conversion  and Geographic Coordinate Conversion  as before. The height is calculated as:

    Geographic Coordinate Conversion 

The iteration can be transformed into the following calculation:

    Geographic Coordinate Conversion 

where Geographic Coordinate Conversion 

The constant Geographic Coordinate Conversion  is a good starter value for the iteration when Geographic Coordinate Conversion . Bowring showed that the single iteration produces a sufficiently accurate solution. He used extra trigonometric functions in his original formulation.

Ferrari's solution

The quartic equation of Geographic Coordinate Conversion , derived from the above, can be solved by Ferrari's solution to yield:

    Geographic Coordinate Conversion 
The application of Ferrari's solution

A number of techniques and algorithms are available but the most accurate, according to Zhu, is the following procedure established by Heikkinen, as cited by Zhu. It is assumed that geodetic parameters Geographic Coordinate Conversion  are known

    Geographic Coordinate Conversion 

Note: arctan2[Y, X] is the four-quadrant inverse tangent function.

Power series

For small e2 the power series

    Geographic Coordinate Conversion 

starts with

    Geographic Coordinate Conversion 

Geodetic to/from ENU coordinates

To convert from geodetic coordinates to local tangent plane (ENU) coordinates is a two-stage process:

  1. Convert geodetic coordinates to ECEF coordinates
  2. Convert ECEF coordinates to local ENU 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 Geographic Coordinate Conversion  and an aircraft at Geographic Coordinate Conversion , then the vector pointing from the radar to the aircraft in the ENU frame is

    Geographic Coordinate Conversion 

Note: Geographic Coordinate Conversion  is the geodetic latitude; the geocentric latitude is inappropriate for representing vertical direction for the local tangent plane and must be converted if necessary.

From ENU to ECEF

This is just the inversion of the ECEF to ENU transformation so

    Geographic Coordinate Conversion 

Conversion across map projections

Conversion of coordinates and map positions among different map projections reference to the same datum may be accomplished either through direct translation formulas from one projection to another, or by first converting from a projection Geographic Coordinate Conversion  to an intermediate coordinate system, such as ECEF, then converting from ECEF to projection Geographic Coordinate Conversion . The formulas involved can be complex and in some cases, such as in the ECEF to geodetic conversion above, the conversion has no closed-form solution and approximate methods must be used. References such as the DMA Technical Manual 8358.1 and the USGS paper Map Projections: A Working Manual contain formulas for conversion of map projections. It is common to use computer programs to perform coordinate conversion tasks, such as with the DoD and NGA supported GEOTRANS program.

Datum transformations

Geographic Coordinate Conversion 
The different possible paths for transforming geographic coordinates from datum Geographic Coordinate Conversion  to datum Geographic Coordinate Conversion 

Transformations among datums can be accomplished in a number of ways. There are transformations that directly convert geodetic coordinates from one datum to another. There are more indirect transforms that convert from geodetic coordinates to ECEF coordinates, transform the ECEF coordinates from one datum to another, then transform ECEF coordinates of the new datum back to geodetic coordinates. There are also grid-based transformations that directly transform from one (datum, map projection) pair to another (datum, map projection) pair.

Helmert transformation

Use of the Helmert transform in the transformation from geodetic coordinates of datum Geographic Coordinate Conversion  to geodetic coordinates of datum Geographic Coordinate Conversion  occurs in the context of a three-step process:

  1. Convert from geodetic coordinates to ECEF coordinates for datum Geographic Coordinate Conversion 
  2. Apply the Helmert transform, with the appropriate Geographic Coordinate Conversion  transform parameters, to transform from datum Geographic Coordinate Conversion  ECEF coordinates to datum Geographic Coordinate Conversion  ECEF coordinates
  3. Convert from ECEF coordinates to geodetic coordinates for datum Geographic Coordinate Conversion 

In terms of ECEF XYZ vectors, the Helmert transform has the form (position vector transformation convention and very small rotation angles simplification)

    Geographic Coordinate Conversion 

The Helmert transform is a seven-parameter transform with three translation (shift) parameters Geographic Coordinate Conversion , three rotation parameters Geographic Coordinate Conversion  and one scaling (dilation) parameter Geographic Coordinate Conversion . The Helmert transform is an approximate method that is accurate when the transform parameters are small relative to the magnitudes of the ECEF vectors. Under these conditions, the transform is considered reversible.

A fourteen-parameter Helmert transform, with linear time dependence for each parameter,: 131-133  can be used to capture the time evolution of geographic coordinates dues to geomorphic processes, such as continental drift and earthquakes. This has been incorporated into software, such as the Horizontal Time Dependent Positioning (HTDP) tool from the U.S. NGS.

Molodensky-Badekas transformation

To eliminate the coupling between the rotations and translations of the Helmert transform, three additional parameters can be introduced to give a new XYZ center of rotation closer to coordinates being transformed. This ten-parameter model is called the Molodensky-Badekas transformation and should not be confused with the more basic Molodensky transform.: 133-134 

Like the Helmert transform, using the Molodensky-Badekas transform is a three-step process:

  1. Convert from geodetic coordinates to ECEF coordinates for datum Geographic Coordinate Conversion 
  2. Apply the Molodensky-Badekas transform, with the appropriate Geographic Coordinate Conversion  transform parameters, to transform from datum Geographic Coordinate Conversion  ECEF coordinates to datum Geographic Coordinate Conversion  ECEF coordinates
  3. Convert from ECEF coordinates to geodetic coordinates for datum Geographic Coordinate Conversion 

The transform has the form

    Geographic Coordinate Conversion 

where Geographic Coordinate Conversion  is the origin for the rotation and scaling transforms and Geographic Coordinate Conversion  is the scaling factor.

The Molodensky-Badekas transform is used to transform local geodetic datums to a global geodetic datum, such as WGS 84. Unlike the Helmert transform, the Molodensky-Badekas transform is not reversible due to the rotational origin being associated with the original datum.: 134 

Molodensky transformation

The Molodensky transformation converts directly between geodetic coordinate systems of different datums without the intermediate step of converting to geocentric coordinates (ECEF). It requires the three shifts between the datum centers and the differences between the reference ellipsoid semi-major axes and flattening parameters.

The Molodensky transform is used by the National Geospatial-Intelligence Agency (NGA) in their standard TR8350.2 and the NGA supported GEOTRANS program. The Molodensky method was popular before the advent of modern computers and the method is part of many geodetic programs.

Grid-based method

Geographic Coordinate Conversion 
Magnitude of shift in position between NAD27 and NAD83 datum as a function of location.

Grid-based transformations directly convert map coordinates from one (map-projection, geodetic datum) pair to map coordinates of another (map-projection, geodetic datum) pair. An example is the NADCON method for transforming from the North American Datum (NAD) 1927 to the NAD 1983 datum. The High Accuracy Reference Network (HARN), a high accuracy version of the NADCON transforms, have an accuracy of approximately 5 centimeters. The National Transformation version 2 (NTv2) is a Canadian version of NADCON for transforming between NAD 1927 and NAD 1983. HARNs are also known as NAD 83/91 and High Precision Grid Networks (HPGN). Subsequently, Australia and New Zealand adopted the NTv2 format to create grid-based methods for transforming among their own local datums.

Like the multiple regression equation transform, grid-based methods use a low-order interpolation method for converting map coordinates, but in two dimensions instead of three. The NOAA provides a software tool (as part of the NGS Geodetic Toolkit) for performing NADCON transformations.

Multiple regression equations

Datum transformations through the use of empirical multiple regression methods were created to achieve higher accuracy results over small geographic regions than the standard Molodensky transformations. MRE transforms are used to transform local datums over continent-sized or smaller regions to global datums, such as WGS 84. The standard NIMA TM 8350.2, Appendix D, lists MRE transforms from several local datums to WGS 84, with accuracies of about 2 meters.

The MREs are a direct transformation of geodetic coordinates with no intermediate ECEF step. Geodetic coordinates Geographic Coordinate Conversion  in the new datum Geographic Coordinate Conversion  are modeled as polynomials of up to the ninth degree in the geodetic coordinates Geographic Coordinate Conversion  of the original datum Geographic Coordinate Conversion . For instance, the change in Geographic Coordinate Conversion  could be parameterized as (with only up to quadratic terms shown): 9 

    Geographic Coordinate Conversion 

where

    Geographic Coordinate Conversion  parameters fitted by multiple regression
    Geographic Coordinate Conversion 
    Geographic Coordinate Conversion  scale factor
    Geographic Coordinate Conversion  origin of the datum, Geographic Coordinate Conversion 

with similar equations for Geographic Coordinate Conversion  and Geographic Coordinate Conversion . Given a sufficient number of Geographic Coordinate Conversion  coordinate pairs for landmarks in both datums for good statistics, multiple regression methods are used to fit the parameters of these polynomials. The polynomials, along with the fitted coefficients, form the multiple regression equations.

See also

References

Tags:

Geographic Coordinate Conversion Change of units and formatGeographic Coordinate Conversion Coordinate system conversionGeographic Coordinate Conversion Datum transformationsGeographic Coordinate ConversionCartographyGeodesyGeodetic datumGeographic coordinate systemGeographic information systemsNavigationSurveying

🔥 Trending searches on Wiki English:

Reba McEntireSteve JobsWashington, D.C.Murder of Grace MillaneJames VI and IThe Gentlemen (2024 TV series)J. Robert OppenheimerGame of ThronesBoeingBob DylanKorean WarThe Acolyte (TV series)ArgylleWordleLondonAtomic bombings of Hiroshima and NagasakiBrazil national football teamCrocus City Hall attackShohei OhtaniGiancarlo Esposito2026 FIFA World Cup qualificationWikiLeaksNarendra ModiThe Rookie (TV series)Rachin RavindraLensey NamiokaBob MarleyMurder MubarakWales national football teamMillennials2024 Indian Premier LeagueJoe KeeryList of dates for EasterCordarrelle PattersonYellowstone (American TV series)Jake Paul vs. Mike TysonDua LipaKatrina LawRonna McDanielIppei MizuharaThe Bachelor (American TV series)Cody RhodesMao ZedongKen MilesLisa LillienMarilyn MonroeHelldivers 2Candace OwensSocial mediaCatherine, Princess of WalesAndre Jin CoquillardCrocus City Hall2026 FIFA World Cup qualification (UEFA)Angelina JolieI-35W Mississippi River bridgeBad Boys Ride or DiePoor Things (film)Sagrada FamíliaDan + ShayArtificial intelligenceThe Bachelorette (American TV series)Justin BieberMexicoAnatomy of a FallRussian invasion of UkraineSri LankaThe Regime (miniseries)Derek DraperAlbert EinsteinMiranda CosgroveSpain national football teamBaltimore Harbor TunnelList of Marvel Cinematic Universe filmsKate GarrawayXHamsterGilbert's syndrome🡆 More