EPSG:9836

Geocentric/topocentric conversions

Attributes

Data source: OGP

Information source: EPSG guidance note #7-2, http://www.epsg.org

Revision date: 2017-06-13

MapTiler banner

Export

Definition: OGP XML

<div class="syntax"><pre><span class="cp">&lt;?xml version=&quot;1.0&quot; encoding=&quot;UTF-8&quot;?&gt;</span> <span class="nt">&lt;gml:OperationMethod</span> <span class="na">xmlns:epsg=</span><span class="s">&quot;urn:x-ogp:spec:schema-xsd:EPSG:1.0:dataset&quot;</span> <span class="na">xmlns:gml=</span><span class="s">&quot;http://www.opengis.net/gml/3.2&quot;</span> <span class="na">xmlns:xlink=</span><span class="s">&quot;http://www.w3.org/1999/xlink&quot;</span> <span class="na">gml:id=</span><span class="s">&quot;iogp-method-9836&quot;</span><span class="nt">&gt;</span> <span class="nt">&lt;gml:metaDataProperty&gt;</span> <span class="nt">&lt;epsg:CommonMetaData&gt;</span> <span class="nt">&lt;epsg:informationSource&gt;</span>EPSG guidance note #7-2, http://www.epsg.org<span class="nt">&lt;/epsg:informationSource&gt;</span> <span class="nt">&lt;epsg:revisionDate&gt;</span>2017-06-13<span class="nt">&lt;/epsg:revisionDate&gt;</span> <span class="nt">&lt;epsg:changes&gt;</span> <span class="nt">&lt;epsg:changeID</span> <span class="na">xlink:href=</span><span class="s">&quot;urn:ogc:def:change-request:EPSG::2017.018&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;/epsg:changes&gt;</span> <span class="nt">&lt;epsg:show&gt;</span>true<span class="nt">&lt;/epsg:show&gt;</span> <span class="nt">&lt;epsg:isDeprecated&gt;</span>false<span class="nt">&lt;/epsg:isDeprecated&gt;</span> <span class="nt">&lt;/epsg:CommonMetaData&gt;</span> <span class="nt">&lt;/gml:metaDataProperty&gt;</span> <span class="nt">&lt;gml:metaDataProperty&gt;</span> <span class="nt">&lt;epsg:CoordinateOperationMethodMetaData&gt;</span> <span class="nt">&lt;epsg:isOperationReversible&gt;</span>true<span class="nt">&lt;/epsg:isOperationReversible&gt;</span> <span class="nt">&lt;epsg:signReversal</span> <span class="na">changeSign=</span><span class="s">&quot;false&quot;</span> <span class="na">xlink:href=</span><span class="s">&quot;urn:ogc:def:parameter:EPSG::8837&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;epsg:signReversal</span> <span class="na">changeSign=</span><span class="s">&quot;false&quot;</span> <span class="na">xlink:href=</span><span class="s">&quot;urn:ogc:def:parameter:EPSG::8838&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;epsg:signReversal</span> <span class="na">changeSign=</span><span class="s">&quot;false&quot;</span> <span class="na">xlink:href=</span><span class="s">&quot;urn:ogc:def:parameter:EPSG::8839&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;epsg:example&gt;</span>For Geocentric CRS = WGS 84 (EPSG CRS code 4978) and Topocentric origin Xo = 3652 755.3058 m Topocentric origin Yo = 319 574.6799 m Topocentric origin Zo = 5201 547.3536 m Ellipsoid parameters: a = 6378137m.0 1/f = 298.2572236 First calculate additional ellipsoid parameters: e^2 = 0.006694380 eta = 0.006739497 b = 6356752.314 Next, derive Po, Lo from Xo,Yo,Zo by the formulas given in method 9602: p = 3666708.2376 q = 0.9583523313 Po = 0.9599310885 rad Lo = 0.0872664625 rad Forward calculation for point with geocentric coordinates: X= 3771 793.968 Y= 140 253.342 Z= 5124 304.349 gives topocentric coordinates U= -189 013.869 V= -128 642.040 W= -4 220.171<span class="nt">&lt;/epsg:example&gt;</span> <span class="nt">&lt;/epsg:CoordinateOperationMethodMetaData&gt;</span> <span class="nt">&lt;/gml:metaDataProperty&gt;</span> <span class="nt">&lt;gml:identifier</span> <span class="na">codeSpace=</span><span class="s">&quot;IOGP&quot;</span><span class="nt">&gt;</span>urn:ogc:def:method:EPSG::9836<span class="nt">&lt;/gml:identifier&gt;</span> <span class="nt">&lt;gml:name&gt;</span>Geocentric/topocentric conversions<span class="nt">&lt;/gml:name&gt;</span> <span class="nt">&lt;gml:formula&gt;</span>Note: These formulas have been transcribed from EPSG Guidance Note #7-2. Users are encouraged to use that document rather than the text which follows as reference because limitations in the transcription will be avoided. First it is necessary to derive ellipsoidal values Po, Lo of the topocentric origin from their geocentric values Xo, Yo, Zo through the reverse formulas given in method 9602. Then topocentric coordinates [U, V, W] are computed as follows: U = – (X-Xo) sin Lo + (Y-Yo) cos Lo V = – (X–Xo) sin Po cos Lo – (Y–Yo) sin Po sin Lo + (Z–Zo) cos Po W = (X–Xo) cos Po cos Lo + (Y–Yo) cos Po sin Lo + (Z–Zo) sin Po The reverse formulas to calculate geocentric coordinates from topocentric coordinates are: X = Xo – U sin Lo – V sin Po cos Lo + W cos Po cos Lo Y = Yo + U cos Lo – V sin Po sin Lo + W cos Po sin Lo Z = Zo + V cos Po + W sin Po<span class="nt">&lt;/gml:formula&gt;</span> <span class="nt">&lt;gml:generalOperationParameter</span> <span class="na">xlink:href=</span><span class="s">&quot;urn:ogc:def:parameter:EPSG::8837&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;gml:generalOperationParameter</span> <span class="na">xlink:href=</span><span class="s">&quot;urn:ogc:def:parameter:EPSG::8838&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;gml:generalOperationParameter</span> <span class="na">xlink:href=</span><span class="s">&quot;urn:ogc:def:parameter:EPSG::8839&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;/gml:OperationMethod&gt;</span> </pre></div>
<?xml version="1.0" encoding="UTF-8"?>
 <gml:OperationMethod xmlns:epsg="urn:x-ogp:spec:schema-xsd:EPSG:1.0:dataset" xmlns:gml="http://www.opengis.net/gml/3.2" xmlns:xlink="http://www.w3.org/1999/xlink" gml:id="iogp-method-9836">
  <gml:metaDataProperty>
    <epsg:CommonMetaData>
      <epsg:informationSource>EPSG guidance note #7-2, http://www.epsg.org</epsg:informationSource>
      <epsg:revisionDate>2017-06-13</epsg:revisionDate>
      <epsg:changes>
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2017.018" />
      </epsg:changes>
      <epsg:show>true</epsg:show>
      <epsg:isDeprecated>false</epsg:isDeprecated>
    </epsg:CommonMetaData>
  </gml:metaDataProperty>
  <gml:metaDataProperty>
    <epsg:CoordinateOperationMethodMetaData>
      <epsg:isOperationReversible>true</epsg:isOperationReversible>
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8837" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8838" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8839" />
      <epsg:example>For Geocentric CRS = 	WGS 84 (EPSG CRS code 4978)
and		
Topocentric origin Xo 	=	3652 755.3058 m
Topocentric origin Yo 	=	  319 574.6799 m
Topocentric origin Zo	=	5201 547.3536 m

Ellipsoid parameters:	a = 6378137m.0  1/f = 298.2572236

First calculate additional ellipsoid parameters:
e^2 = 0.006694380
eta = 0.006739497
b = 6356752.314

Next, derive Po, Lo from Xo,Yo,Zo by the formulas given in method 9602:
p	=	3666708.2376	
q	=	0.9583523313	
Po	=	0.9599310885	rad
Lo	=	0.0872664625	rad

Forward calculation for point with geocentric coordinates:
X= 3771 793.968
Y=  140 253.342
Z= 5124 304.349

gives topocentric coordinates
U= -189 013.869
V= -128 642.040
W=    -4 220.171</epsg:example>
    </epsg:CoordinateOperationMethodMetaData>
  </gml:metaDataProperty>
  <gml:identifier codeSpace="IOGP">urn:ogc:def:method:EPSG::9836</gml:identifier>
  <gml:name>Geocentric/topocentric conversions</gml:name>
  <gml:formula>Note: These formulas have been transcribed from EPSG Guidance Note #7-2. Users are encouraged to use that document rather than the text which follows as reference because limitations in the transcription will be avoided.

First it is necessary to derive ellipsoidal values Po, Lo of the topocentric origin from their geocentric values Xo, Yo, Zo through the reverse formulas given in method 9602.

Then topocentric coordinates [U, V, W] are computed as follows:

U = – (X-Xo) sin Lo + (Y-Yo) cos Lo
V = – (X–Xo) sin Po cos Lo – (Y–Yo) sin Po sin Lo + (Z–Zo) cos Po
W = (X–Xo) cos Po cos Lo + (Y–Yo) cos Po sin Lo + (Z–Zo) sin Po 


The reverse formulas to calculate geocentric coordinates from topocentric coordinates are:

X = Xo – U sin Lo  – V sin Po cos Lo + W cos Po cos Lo
Y = Yo + U cos Lo – V sin Po sin Lo  + W cos Po sin Lo 
Z = Zo + V cos Po + W sin Po</gml:formula>
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8837" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8838" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8839" />
</gml:OperationMethod>