본문 바로가기
카테고리 없음

Gps Coordinate Conversion Formula

ocnepgypfacolen 2020. 8. 31. 13:19



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.

Convert Coordinates - Calculate a position in a variety of formats. For lat/lon, the letters 'N, S, E, and W can come before or after the coordinates. So long as the coordinate can be understood, punctuation, spaces, and °' are optional. Acceptable formats include. Jan 18, 2019  Add GPS coordinates to address. If you are looking to add GPS coordinates to an address, please go to this page. How to write a GPS coordinate in Excel. GPS coordinates are generally written using the symbols ° (degrees), ' (minutes), and ' (seconds). But if you keep your coordinate as a string, no calculation is possible.

In geodesy, geographic coordinate conversion is defined as translation among different coordinate formats or map projections all referenced to the same geodetic datum.[1] 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.

Gps

  • 2Coordinate system conversion
    • 2.2From ECEF to geodetic coordinates
      • 2.2.2Ferrari's solution
    • 2.3Geodetic to/from ENU coordinates
  • 3Datum transformations

Change of units and format[edit]

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:[2]

  • sexagesimal degree: degrees, minutes, and seconds : 40° 26′ 46″ N 79° 58′ 56″ W
  • degrees and decimal minutes: 40° 26.767′ N 79° 58.933′ W
  • decimal degrees: 40.446° N 79.982° W

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

decimaldegrees=degrees+minutes/60+seconds/3600{displaystyle {rm {{decimal degrees}={rm {{degrees}+{rm {{minutes}/60+{rm {{seconds}/3600}}}}}}}}}.

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

degrees=⌊decimaldegrees⌋minutes=⌊60×(decimaldegrees−degrees)⌋seconds=3600×(decimaldegrees−degrees)−60×minutes{displaystyle {begin{aligned}{rm {degrees}}&=lfloor {rm {{decimal degrees}rfloor }}{rm {minutes}}&=lfloor 60times ({rm {{decimal degrees}-{rm {{degrees})rfloor }}}}{rm {seconds}}&=3600times ({rm {{decimal degrees}-{rm {{degrees})-60times {rm {minutes}}}}}}end{aligned}}}

where the notation ⌊x⌋{displaystyle lfloor xrfloor } means take the integer part of x{displaystyle x} and is called a floor function.

Coordinate system conversion[edit]

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 ECEF coordinates and conversion from one type of map projection to another.

From geodetic to ECEF coordinates[edit]

The length PQ is called prime vertical radius (N(ϕ){displaystyle ,N(phi )}). The length IQ is equal to e2N(ϕ){displaystyle ,e^{2}N(phi )}. R = (X,Y,Z){displaystyle ,(X,Y,Z)}.

Geodetic coordinates (latitude ϕ{displaystyle phi }, longitude λ{displaystyle lambda }, height h{displaystyle h}) can be converted into ECEF coordinates using the following equation:[3]

X=(N(ϕ)+h)cos⁡ϕcos⁡λY=(N(ϕ)+h)cos⁡ϕsin⁡λZ=(b2a2N(ϕ)+h)sin⁡ϕ{displaystyle {begin{aligned}X&=left(N(phi )+hright)cos {phi }cos {lambda }Y&=left(N(phi )+hright)cos {phi }sin {lambda }Z&=left({frac {b^{2}}{a^{2}}}N(phi )+hright)sin {phi }end{aligned}}}

where

N(ϕ)=a2a2cos2⁡ϕ+b2sin2⁡ϕ=a1−e2sin2⁡ϕ,{displaystyle N(phi )={frac {a^{2}}{sqrt {a^{2}cos ^{2}phi +b^{2}sin ^{2}phi }}}={frac {a}{sqrt {1-e^{2}sin ^{2}phi }}},}

and a{displaystyle a} and b{displaystyle b} are the equatorial radius (semi-major axis) and the polar radius (semi-minor axis), respectively. e2=1−b2a2{displaystyle e^{2}=1-{frac {b^{2}}{a^{2}}}} is the square of the first numerical eccentricity of the ellipsoid. The prime vertical radius of curvatureN(ϕ){displaystyle ,N(phi )} is the distance from the surface to the Z-axis along the ellipsoid normal (see 'Radius of curvature on the Earth'). The parameter h{displaystyle h} is eliminated by subtracting

X2+Y2cos⁡ϕ=N+h{displaystyle {frac {sqrt {X^{2}+Y^{2}}}{cos phi }}=N+h}

and

Zsin⁡ϕ=b2a2N+h,{displaystyle {frac {Z}{sin phi }}={frac {b^{2}}{a^{2}}}N+h,}

Gps coordinates to map

so the following equation holds:

pcos⁡ϕ−Zsin⁡ϕ−e2N(ϕ)=0,{displaystyle {frac {p}{cos phi }}-{frac {Z}{sin phi }}-e^{2}N(phi )=0,}

where p=X2+Y2{displaystyle p={sqrt {X^{2}+Y^{2}}}}.

The orthogonality of the coordinates is confirmed via differentiation:

(dXdYdZ)=(−sin⁡λ−sin⁡ϕcos⁡λcos⁡ϕcos⁡λcos⁡λ−sin⁡ϕsin⁡λcos⁡ϕsin⁡λ0cos⁡ϕsin⁡ϕ)(dEdNdU),(dEdNdU)=((N(ϕ)+h)cos⁡ϕ000M(ϕ)+h0001)(dλdϕdh),{displaystyle {begin{aligned}{begin{pmatrix}dXdYdZend{pmatrix}}&={begin{pmatrix}-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{pmatrix}}{begin{pmatrix}dEdNdUend{pmatrix}},{begin{pmatrix}dEdNdUend{pmatrix}}&={begin{pmatrix}left(N(phi )+hright)cos phi &0&00&M(phi )+h&00&0&1end{pmatrix}}{begin{pmatrix}dlambda dphi dhend{pmatrix}},end{aligned}}}

where

M(ϕ)=a(1−e2)(1−e2sin2⁡ϕ)3/2{displaystyle M(phi )={frac {a(1-e^{2})}{left(1-e^{2}sin ^{2}phi right)^{3/2}}}}

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

From ECEF to geodetic coordinates[edit]

The conversion of ECEF coordinates to geodetic coordinates (such WGS84) involves more trigonometry but is sensitive to small accuracy due to Rn and h being maybe 10^6 apart,[4][5] but the longitude is the same as that of the geocentric, λ{displaystyle ,lambda }.

There are several methods that solve the equation; two are shown.

Newton–Raphson method[edit]

The following Bowring's irrational geodetic-latitude equation[6] is efficient to be solved by Newton–Raphson iteration method:[7][8]

κ−1−e2aκp2+(1−e2)Z2κ2=0,{displaystyle kappa -1-{frac {e^{2}akappa }{sqrt {p^{2}+(1-e^{2})Z^{2}kappa ^{2}}}}=0,}

where κ=pZtan⁡ϕ{displaystyle kappa ={frac {p}{Z}}tan phi } and p=x2+y2.{displaystyle p={sqrt {x^{2}+y^{2}}}.} The height is calculated as:

h=e−2(κ−1−κ0−1)p2+Z2κ2,{displaystyle h=e^{-2}(kappa ^{-1}-{kappa _{0}}^{-1}){sqrt {p^{2}+Z^{2}kappa ^{2}}},}
κ0=(1−e2)−1.{displaystyle kappa _{0}=left(1-e^{2}right)^{-1}.}

The iteration can be transformed into the following calculation:

κi+1=ci+(1−e2)Z2κi3ci−p2=1+p2+(1−e2)Z2κi3ci−p2,{displaystyle 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}}},}

where ci=(p2+(1−e2)Z2κi2)3/2ae2.{displaystyle c_{i}={frac {left(p^{2}+left(1-e^{2}right)Z^{2}kappa _{i}^{2}right)^{3/2}}{ae^{2}}}.}

The constant κ0{displaystyle ,kappa _{0}} is a good starter value for the iteration when h≈0{displaystyle happrox 0}. Bowring showed that the single iteration produces a sufficiently accurate solution. He used extra trigonometric functions in his original formulation.

Ferrari's solution[edit]

The quartic equation of κ{displaystyle kappa }, derived from the above, can be solved by Ferrari's solution[9][10] to yield:

ζ=(1−e2)z2/a2,ρ=(p2/a2+ζ−e4)/6,s=e4ζp2/(4ρ3a2),t=1+s+s(s+2)3,u=ρ(1+t+1/t),v=u2+e4ζ,w=e2(u+v−ζ)/(2v),κ=1+e2(u+v+w2+w)/(u+v).{displaystyle {begin{aligned}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}/(4rho ^{3}a^{2}),[6pt]t&={sqrt[{3}]{1+s+{sqrt {s(s+2)}}}},[6pt]u&=rho (1+t+1/t),[6pt]v&={sqrt {u^{2}+e^{4}zeta }},[6pt]w&=e^{2}(u+v-zeta )/(2v),[6pt]kappa &=1+e^{2}({sqrt {u+v+w^{2}}}+w)/(u+v).end{aligned}}}
The application of Ferrari's solution[edit]

A number of techniques and algorithms are available but the most accurate, according to Zhu,[11] is the following procedure established by Heikkinen,[12] as cited by Zhu. It is assumed that geodetic parameters {a,b,e}{displaystyle {a,b,e}} are known

r=X2+Y2e′2=(a2−b2)/b2F=54b2Z2G=r2+(1−e2)Z2−e2(a2−b2)c=e4Fr2G3s=1+c+c2+2c3P=F3(s+1s+1)2G2Q=1+2e4Pr0=−(Pe2r)1+Q+12a2(1+1/Q)−P(1−e2)Z2Q(1+Q)−12Pr2U=(r−e2r0)2+Z2V=(r−e2r0)2+(1−e2)Z2z0=b2ZaVh=U(1−b2aV)ϕ=arctan⁡[Z+e′2z0r]λ=arctan⁡2[Y,X]{displaystyle {begin{matrix}r&=&{sqrt {X^{2}+Y^{2}}}e'^{2}&=&(a^{2}-b^{2})/b^{2}F&=&54b^{2}Z^{2}G&=&r^{2}+(1-e^{2})Z^{2}-e^{2}(a^{2}-b^{2})c&=&{frac {e^{4}Fr^{2}}{G^{3}}}s&=&{sqrt[{3}]{1+c+{sqrt {c^{2}+2c}}}}P&=&{frac {F}{3left(s+{frac {1}{s}}+1right)^{2}G^{2}}}Q&=&{sqrt {1+2e^{4}P}}r_{0}&=&{frac {-(Pe^{2}r)}{1+Q}}+{sqrt {{frac {1}{2}}a^{2}left(1+1/Qright)-{frac {P(1-e^{2})Z^{2}}{Q(1+Q)}}-{frac {1}{2}}Pr^{2}}}U&=&{sqrt {(r-e^{2}r_{0})^{2}+Z^{2}}}V&=&{sqrt {(r-e^{2}r_{0})^{2}+(1-e^{2})Z^{2}}}z_{0}&=&{frac {b^{2}Z}{aV}}h&=&Uleft(1-{frac {b^{2}}{aV}}right)phi &=&arctan left[{frac {Z+e'^{2}z_{0}}{r}}right]lambda &=&arctan 2[Y,X]end{matrix}}}

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

Power Series[edit]

For small e2 the power series

κ=∑i≥0αie2i{displaystyle kappa =sum _{igeq 0}alpha _{i}e^{2i}}

starts with

α0=1;{displaystyle alpha _{0}=1;}
α1=aZ2+p2;{displaystyle alpha _{1}={frac {a}{sqrt {Z^{2}+p^{2}}}};}
α2=aZ2Z2+p2+2a2p22(Z2+p2)2.{displaystyle alpha _{2}={frac {aZ^{2}{sqrt {Z^{2}+p^{2}}}+2a^{2}p^{2}}{2(Z^{2}+p^{2})^{2}}}.}

Geodetic to/from ENU coordinates[edit]

To convert from geodetic coordinates to local 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[edit]

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 {Xr,Yr,Zr}{displaystyle {X_{r},Y_{r},Z_{r}}} and an aircraft at {Xp,Yp,Zp}{displaystyle {X_{p},Y_{p},Z_{p}}} then the vector pointing from the radar to the aircraft in the ENU frame is

[xyz]=[−sin⁡λrcos⁡λr0−sin⁡ϕrcos⁡λr−sin⁡ϕrsin⁡λrcos⁡ϕrcos⁡ϕrcos⁡λrcos⁡ϕrsin⁡λrsin⁡ϕr][Xp−XrYp−YrZp−Zr]{displaystyle {begin{bmatrix}xyzend{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}}}

Note:ϕ{displaystyle phi } is the geodetic latitude. A prior version of this page showed use of the geocentric latitude (ϕ′{displaystyle phi ^{prime }}). 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:

tan⁡ϕ′=ZrXr2+Yr2=N(ϕ)(1−f)2+hN(ϕ)+htan⁡ϕ{displaystyle 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 }

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.'

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).

tan⁡λ=YrXr{displaystyle tan lambda ={frac {Y_{r}}{X_{r}}}}

Note: Unambiguous determination of ϕ{displaystyle phi } and λ{displaystyle lambda } requires knowledge of which quadrant the coordinates lie in.

From ENU to ECEF[edit]

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

[XYZ]=[−sin⁡λ−sin⁡ϕcos⁡λcos⁡ϕcos⁡λcos⁡λ−sin⁡ϕsin⁡λcos⁡ϕsin⁡λ0cos⁡ϕsin⁡ϕ][xyz]+[XrYrZr]{displaystyle {begin{bmatrix}XYZend{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}xyzend{bmatrix}}+{begin{bmatrix}X_{r}Y_{r}Z_{r}end{bmatrix}}}

Conversion across map projections[edit]

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 A{displaystyle A} to an intermediate coordinate system, such as ECEF, then converting from ECEF to projection B{displaystyle B}. 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[13] and the USGS paper Map Projections: A Working Manual[14] 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.[15]

Datum transformations[edit]

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 the 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.

Grid-based method[edit]

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.[16] 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).[17] 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.[18][19]

Molodensky transformation[edit]

The Molodensky transformation converts directly between geodetic coordinate systems of different datums without the intermediate step of converting to geocentric coordinates (ECEF).[20] 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.[21] The Molodensky method was popular before the advent of modern computers and the method is part of many geodetic programs.

Multiple regression equations[edit]

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.[22] The standard NIMA TM 8350.2, Appendix D,[23] lists MRE transforms from several local datums to WGS 84, with accuracies of about 2 meters.[24]

The MREs are a direct transformation of geodetic coordinates with no intermediate ECEF step. Geodetic coordinates ϕB,λB,hB{displaystyle phi _{B},lambda _{B},h_{B}} in the new datum B{displaystyle B} are modeled as polynomials of up to the ninth degree in the geodetic coordinates ϕA,λA,hA{displaystyle phi _{A},lambda _{A},h_{A}} of the original datum A{displaystyle A}. For instance, the change in ϕB{displaystyle phi _{B}} could be parameterized as (with only up to quadratic terms shown)[22]:9

Δϕ=a0+a1U+a2V+a3U2+a4UV+a5V2+⋯{displaystyle Delta phi =a_{0}+a_{1}U+a_{2}V+a_{3}U^{2}+a_{4}UV+a_{5}V^{2}+cdots }

where

ai=parametersfittedbymultipleregressionU=K(ϕA−ϕm)V=K(λA−λm)K=scalefactorϕm,λm=originofthedatumA{displaystyle {begin{aligned}a_{i}&={rm {parameters fitted by multiple regression}}U&=K(phi _{A}-phi _{m})V&=K(lambda _{A}-lambda _{m})K&={rm {scale factor}}phi _{m},lambda _{m}&={rm {{origin of the datum }A}}end{aligned}}}

with similar equations for Δλ{displaystyle Delta lambda } and Δh{displaystyle Delta h}. Given a sufficient number of (A,B){displaystyle (A,B)} 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.

Helmert transformation[edit]

Use of the Helmert transform in the transformation from geodetic coordinates of datum A{displaystyle A} to geodetic coordinates of datum B{displaystyle B} occurs in the context of a three-step process:[25]

  1. Convert from geodetic coordinates to ECEF coordinates for datum A{displaystyle A}
  2. Apply the Helmert transform, with the appropriate A→B{displaystyle Ato B} transform parameters, to transform from datum A{displaystyle A} ECEF coordinates to datum B{displaystyle B} ECEF coordinates
  3. Convert from ECEF coordinates to geodetic coordinates for datum B{displaystyle B}

In terms of ECEF XYZ vectors, the Helmert transform has the form[25]

[XBYBZB]=[cxcycz]+(1+s×10−6)⋅[1−rzryrz1−rx−ryrx1]⋅[XAYAZA].{displaystyle {begin{bmatrix}X_{B}Y_{B}Z_{B}end{bmatrix}}={begin{bmatrix}c_{x}c_{y}c_{z}end{bmatrix}}+(1+stimes 10^{-6})cdot {begin{bmatrix}1&-r_{z}&r_{y}r_{z}&1&-r_{x}-r_{y}&r_{x}&1end{bmatrix}}cdot {begin{bmatrix}X_{A}Y_{A}Z_{A}end{bmatrix}}.}

The Helmert transform is a seven-parameter transform with three translation (shift) parameters cx,cy,cz{displaystyle c_{x},c_{y},c_{z}}, three rotation parameters rx,ry,rz{displaystyle r_{x},r_{y},r_{z}} and one scaling (dilation) parameter s{displaystyle s}. 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.[26]

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

Molodensky-Badekas transformation[edit]

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.[26]: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 A{displaystyle A}
  2. Apply the Molodensky-Badekas transform, with the appropriate A→B{displaystyle Ato B} transform parameters, to transform from datum A{displaystyle A} ECEF coordinates to datum B{displaystyle B} ECEF coordinates
  3. Convert from ECEF coordinates to geodetic coordinates for datum B{displaystyle B}

The transform has the form[30]

Convert Gps Coordinates To Address

[XBYBZB]=[XAYAZA]+[ΔXAΔYAΔZA]+[1−rzryrz1−rx−ryrx1]⋅[XA−XA0YA−YA0ZA−ZA0]+ΔS[XA−XA0YA−YA0ZA−ZA0].{displaystyle {begin{bmatrix}X_{B}Y_{B}Z_{B}end{bmatrix}}={begin{bmatrix}X_{A}Y_{A}Z_{A}end{bmatrix}}+{begin{bmatrix}Delta X_{A}Delta Y_{A}Delta Z_{A}end{bmatrix}}+{begin{bmatrix}1&-r_{z}&r_{y}r_{z}&1&-r_{x}-r_{y}&r_{x}&1end{bmatrix}}cdot {begin{bmatrix}X_{A}-X_{A}^{0}Y_{A}-Y_{A}^{0}Z_{A}-Z_{A}^{0}end{bmatrix}}+Delta S{begin{bmatrix}X_{A}-X_{A}^{0}Y_{A}-Y_{A}^{0}Z_{A}-Z_{A}^{0}end{bmatrix}}.}

where (XA0,YA0,ZA0){displaystyle (X_{A}^{0},Y_{A}^{0},Z_{A}^{0})} is the origin for the rotation and scaling transforms and ΔS{displaystyle Delta S} 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.[26]:134

See also[edit]

References[edit]

  1. ^Roger Foster; Dan Mullaney. 'Basic Geodesy Article 018: Conversions and Transformations'(PDF). National Geospatial Intelligence Agency. Retrieved 4 March 2014.
  2. ^'Coordinate transformer'. Ordnance Survey Great Britain. Retrieved 4 March 2014.
  3. ^B. Hofmann-Wellenhof; H. Lichtenegger; J. Collins. GPS - theory and practice. Section 10.2.1. p. 282. ISBN3-211-82839-7.
  4. ^R. Burtch, A Comparison of Methods Used in Rectangular to Geodetic Coordinate Transformations.[permanent dead link]
  5. ^Featherstone, W. E.; Claessens, S. J. (2008). 'Closed-Form Transformation between Geodetic and Ellipsoidal Coordinates'. Stud. Geophys. Geod. 52 (1): 1–18. doi:10.1007/s11200-008-0002-6.
  6. ^Bowring, B. R. (1976). 'Transformation from Spatial to Geographical Coordinates'. Surv. Rev. 23 (181): 323–327. doi:10.1179/003962676791280626.
  7. ^Fukushima, T. (1999). 'Fast Transform from Geocentric to Geodetic Coordinates'. J. Geod. 73 (11): 603–610. doi:10.1007/s001900050271. (Appendix B)
  8. ^Sudano, J. J. (1997). 'An exact conversion from an earth-centered coordinate system to latitude, longitude and altitude'. doi:10.1109/NAECON.1997.622711.Cite journal requires |journal= (help)
  9. ^Vermeille, H., H. (2002). 'Direct Transformation from Geocentric to Geodetic Coordinates'. J. Geod. 76 (8): 451–454. doi:10.1007/s00190-002-0273-6.
  10. ^Gonzalez-Vega, Laureano; PoloBlanco, Irene (2009). 'A symbolic analysis of Vermeille and Borkowski polynomials for transforming 3D Cartesian to geodetic coordinates'. J. Geod. 83 (11): 1071–1081. doi:10.1007/s00190-009-0325-2.
  11. ^Zhu, J. (1994). 'Conversion of Earth-centered Earth-fixed coordinates to geodetic coordinates'. IEEE Transactions on Aerospace and Electronic Systems. 30: 957–961. doi:10.1109/7.303772.
  12. ^Heikkinen, M. (1982). 'Geschlossene formeln zur berechnung räumlicher geodätischer koordinaten aus rechtwinkligen koordinaten'. Z. Vermess. (in German). 107: 207–211.
  13. ^'TM8358.2: The Universal Grids: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS)'(PDF). National Geospatial-Intelligence Agency. Retrieved 4 March 2014.
  14. ^Snyder, John P. (1987). Map Projections: A Working Manual. USGS Professional Paper: 1395.
  15. ^'MSP GEOTRANS 3.3 (Geographic Translator)'. NGA: Coordinate Systems Analysis Branch. Retrieved 4 March 2014.
  16. ^'ArcGIS Help 10.1: Grid-based methods'. ESRI. Retrieved 5 March 2014.
  17. ^'NADCON/HARN Datum ShiftMethod'. bluemarblegeo.com. Retrieved 5 March 2014.
  18. ^'NADCON - Version 4.2'. NOAA. Retrieved 5 March 2014.
  19. ^Mulcare, Donald M. 'NGS Toolkit, Part 8: The National Geodetic Survey NADCON Tool'. Professional Surveyor Magazine. Archived from the original on 6 March 2014. Retrieved 5 March 2014.
  20. ^'ArcGIS Help 10.1: Equation-based methods'. ESRI. Retrieved 5 March 2014.
  21. ^'Datum Transformations'. National Geospatial-Intelligence Agency. Retrieved 5 March 2014.
  22. ^ abUser's Handbook on Datum Transformations Involving WGS 84(PDF) (Report). Special Publication No. 60 (3rd ed.). Monaco: International Hydrographic Bureau. August 2008. Retrieved 2017-01-10.
  23. ^'DEPARTMENT OF DEFENSE WORLD GEODETIC SYSTEM 1984 Its Definition and Relationships with Local Geodetic Systems'(PDF). National Imagery and Mapping Agency (NIMA). Retrieved 5 March 2014.
  24. ^Taylor, Chuck. 'High-Accuracy Datum Transformations'. Retrieved 5 March 2014.
  25. ^ ab'Equations Used for Datum Transformations'. Land Information New Zealand (LINZ). Retrieved 5 March 2014.
  26. ^ abcd'Geomatics Guidance Note Number 7, part 2 Coordinate Conversions and Transformations including Formulas'(PDF). International Association of Oil and Gas Producers (OGP). Archived from the original(PDF) on 6 March 2014. Retrieved 5 March 2014.
  27. ^Bolstad, Paul. GIS Fundamentals, 4th Edition(PDF). Atlas books. p. 93. ISBN978-0-9717647-3-6. Archived from the original(PDF) on 2016-02-02.
  28. ^'Addendum to NIMA TR 8350.2: Implementation of the World Geodetic System 1984 (WGS 84) Reference Frame G1150'(PDF). National Geospatial-Intelligence Agency. Retrieved 6 March 2014.
  29. ^'HTDP - Horizontal Time-Dependent Positioning'. U.S. National Geodetic Survey (NGS). Retrieved 5 March 2014.
  30. ^'Molodensky-Badekas (7+3) Transformations'. National Geospatial Intelligence Agency (NGA). Retrieved 5 March 2014.

Excel Gps Coordinate Conversion Formula

Retrieved from 'https://en.wikipedia.org/w/index.php?title=Geographic_coordinate_conversion&oldid=916881374'





티스토리툴바