EPSG:9810

Polar Stereographic (variant A)

Attributes

Data source: OGP

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

Revision date: 2018-08-29

Remarks: Latitude of natural origin must be either 90 degrees or -90 degrees (or equivalent in alternative angle unit).

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-9810&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>2018-08-29<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::2003.220&quot;</span> <span class="nt">/&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::2010.054&quot;</span> <span class="nt">/&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.024&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::8801&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::8802&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::8805&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::8806&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::8807&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;epsg:example&gt;</span>For Projected Coordinate Reference System: WGS 84 / UPS North Parameters: Ellipsoid: WGS 84 a = 6378137.0 metre 1/f = 298.2572236 then e = 0.081819191 Latitude of natural origin (latO): 90°00&#39;00.000&quot;N =1.570796327 rad Longitude of origin (longO): 0°00&#39;00.000&quot;E=0.0 rad Scale factor at natural origin (ko): 0.994 False easting (FE) 2000000.00 metre False northing (FN) 2000000.00 metre Forward calculation for: Latitude (lat) =73°N =1.274090354 rad Longitude (lon) =44°E =0.767944871 rad t = 0.150412808 rho = 1900814.564 whence E = 3320416.75 m N = 632668.43 m Reverse calculation for the same Easting and Northing (3320416.75 E, 632668.43 N) first gives: rho&#39; = 1900814.566 t&#39; = 0.150412808 chi = 1.2722090 Then Latitude (lat) = 73°00&#39;00.000&quot;N Longitude (lon) = 44°00&#39;00.000&quot;E<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::9810<span class="nt">&lt;/gml:identifier&gt;</span> <span class="nt">&lt;gml:name&gt;</span>Polar Stereographic (variant A)<span class="nt">&lt;/gml:name&gt;</span> <span class="nt">&lt;gml:remarks&gt;</span>Latitude of natural origin must be either 90 degrees or -90 degrees (or equivalent in alternative angle unit).<span class="nt">&lt;/gml:remarks&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. For the forward conversion from latitude and longitude, for the south pole case E = FE + rho * sin(lon – lonO) N = FN + rho * cos(lon – lonO) where t = tan(pi/4 + lat/2) / {[(1 + e sin(lat)) / (1 – e sin(lat))]^(e/2)} rho = 2*a*ko*t / {[(1+e)^(1+e) (1–e)^(1–e)]^0.5} For the north pole case, rho and E are found as for the south pole case but t = tan(pi/4 – lat/2) * {[(1 + e sin(lat)) / (1 – e sin(lat))]^(e/2)} N = FN – rho * cos(lon – lonO) For the reverse conversion from easting and northing to latitude and longitude, lat = chi + (e^2/2 + 5e^4/24 + e^6/12 + 13e^8/360) sin(2 chi) + (7e^4/48 + 29e^6/240 + 811e^8/11520) sin(4 chi) + (7e^6/120 + 81e^8/1120) sin(6 chi) + (4279e^8/161280) sin(8 chi) where rho&#39; = [(E-FE)^2 + (N – FN)^2]^0.5 t&#39; =rho&#39; {[(1+e)^(1+e) * (1– e)^(1-e)]^0.5} / (2 a ko) and for the south pole case chi = 2 atan(t&#39; ) – pi/2 but for the north pole case chi = pi/2 - 2 atan(t&#39;) Then for for both north and south cases if E = FE, lon = lonO else for the south pole case lon = lonO + atan2[(E – FE),(N – FN)] and for the north pole case lon = lonO + atan2[(E – FE),(FN – N)] (see GN7-2 implementation notes in preface for atan2 convention)<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::8801&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::8802&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::8805&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::8806&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::8807&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-9810">
  <gml:metaDataProperty>
    <epsg:CommonMetaData>
      <epsg:informationSource>EPSG guidance note #7-2, http://www.epsg.org</epsg:informationSource>
      <epsg:revisionDate>2018-08-29</epsg:revisionDate>
      <epsg:changes>
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2003.220" />
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2010.054" />
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2017.024" />
      </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::8801" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8802" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8805" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8806" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8807" />
      <epsg:example>For Projected Coordinate Reference System: WGS 84 / UPS North

Parameters:
Ellipsoid: WGS 84
a = 6378137.0 metre
1/f = 298.2572236
then e = 0.081819191

Latitude of natural origin (latO): 90°00'00.000"N =1.570796327 rad
Longitude of origin (longO): 0°00'00.000"E=0.0 rad
Scale factor at natural origin (ko): 0.994
False easting (FE) 2000000.00 metre
False northing (FN) 2000000.00 metre

Forward calculation for: 
Latitude (lat) =73°N =1.274090354 rad
Longitude (lon) =44°E =0.767944871 rad

t  = 0.150412808
rho = 1900814.564
whence
E = 3320416.75 m
N =  632668.43 m

Reverse calculation for the same Easting and Northing (3320416.75 E, 632668.43 N) first gives:
rho' = 1900814.566
t'  = 0.150412808
chi  = 1.2722090

Then
Latitude (lat) = 73°00'00.000"N
Longitude (lon) = 44°00'00.000"E</epsg:example>
    </epsg:CoordinateOperationMethodMetaData>
  </gml:metaDataProperty>
  <gml:identifier codeSpace="IOGP">urn:ogc:def:method:EPSG::9810</gml:identifier>
  <gml:name>Polar Stereographic (variant A)</gml:name>
  <gml:remarks>Latitude of natural origin must be either 90 degrees or -90 degrees (or equivalent in alternative angle unit).</gml:remarks>
  <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.

For the forward conversion from latitude and longitude, for the south pole case

E = FE + rho * sin(lon – lonO)
N = FN + rho * cos(lon – lonO)
where
t = tan(pi/4 + lat/2) / {[(1 + e sin(lat)) / (1 – e sin(lat))]^(e/2)}
rho = 2*a*ko*t / {[(1+e)^(1+e) (1–e)^(1–e)]^0.5}

For the north pole case, 
rho and E are found as for the south pole case but
t  = tan(pi/4 – lat/2) *  {[(1 + e sin(lat)) / (1 – e sin(lat))]^(e/2)}
N = FN – rho * cos(lon – lonO)

For the reverse conversion from easting and northing to latitude and longitude,
lat = chi + (e^2/2 + 5e^4/24 + e^6/12 + 13e^8/360) sin(2 chi) 
+ (7e^4/48 + 29e^6/240 + 811e^8/11520) sin(4 chi)
+ (7e^6/120 +  81e^8/1120) sin(6 chi)  + (4279e^8/161280) sin(8 chi)

where rho'  = [(E-FE)^2  + (N – FN)^2]^0.5
t'   =rho' {[(1+e)^(1+e) * (1– e)^(1-e)]^0.5} / (2 a ko)
and for the south pole case
chi  = 2 atan(t' ) – pi/2 
but for the north pole case
chi  =  pi/2 - 2 atan(t')

Then for for both north and south cases if E = FE, lon = lonO
else for the south pole case
lon = lonO + atan2[(E – FE),(N – FN)]
and for the north pole case
lon = lonO + atan2[(E – FE),(FN – N)] (see GN7-2 implementation notes in preface for atan2 convention)</gml:formula>
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8801" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8802" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8805" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8806" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8807" />
</gml:OperationMethod>