EPSG:1042

Krovak Modified

Attributes

Data source: OGP

Information source: Land Survey Office (ZU), Prague. www.cuzk.cz/zu

Revision date: 2018-08-29

Remarks: Incorporates a polynomial transformation which is defined to be exact and for practical purposes is considered to be a map projection.

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-1042&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>Land Survey Office (ZU), Prague. www.cuzk.cz/zu<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::2010.071&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.100&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::2012.011&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.037&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::8811&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::8833&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::1036&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::8818&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::8819&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: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::8617&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::8618&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::1026&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::1027&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::1028&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::1029&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::1030&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::1031&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::1032&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::1033&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::1034&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::1035&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;epsg:example&gt;</span>For Projected Coordinate Reference System: S-JTSK/05 (Ferro) / Modified Krovak Parameters: Ellipsoid Bessel 1841 a = 6377397.155m 1/f = 299.15281 then e = 0.081696831 e^2 = 0.006674372 Latitude of projection centre = 49°30&#39;00&quot;N = 0.863937979 rad Longitude of Origin = 42°30&#39;00&quot;E of Ferro = 0.741764932 rad Co-latitude of cone axis = 30°17&#39;17.30311&quot; = 0.528627763 rad Latitude of pseudo standard parallel = 78°30&#39;00&quot;N 1.370083463 rad Scale factor on pseudo Standard Parallel (ko) = 0.9999 False Easting = 5000000.00 m False Northing = 5000000.00 m Ordinate 1 of evaluation point Xo = 1089000.00 m Ordinate 2 of evaluation point Yo = 654000.00 m C1 = 2.946529277E-02 C2 = 2.515965696E-02 C3 = 1.193845912E-07 C4 = -4.668270147E-07 C5 = 9.233980362E-12 C6 = 1.523735715E-12 C7 = 1.696780024E-18 C8 = 4.408314235E-18 C9 = -8.331083518E-24 C10 = -3.689471323E-24 Calculated projection constants: A = 6380703.611 B = 1.000597498 gammao = 0.863239103 to = 1.003419164 n = 0.979924705 ro = 1298039.005 Forward calculation for: Latitude = 50°12&#39;32.442&quot;N = 0.876312568 rad Longitude = 34°30&#39;59.179&quot;E of Ferro = 0.602425500 rad U = 0.875596951 V = 0.139422687 S = 1.386275051 D = 0.506554627 theta = 0.496385393 r = 1194731.002 Xp = 1050538.631 Yp = 568990.995 Xr = -38461.369 Yr = -85009.005 dX = -0.077 dY = 0.088 and Southing X = 6050538.71 m Westing Y = 5568990.91 m Reverse calculation for the same Southing and Westing: Xr&#39; = -38461.292 Yr&#39; = -85009.093 dX&#39; = -0.077 dY&#39; = 0.088 Xp&#39; = 1050538.631 Yp&#39; = 568990.995 r&#39; = 1194731.002 theta&#39; = 0.496385393 D&#39; = 0.506554627 T&#39; = 1.386275051 U&#39; = 0.875596951 V&#39; = 0.139422687 lat(iteration 1) = 0.876310603 lat(iteration 2) = 0.876312562 lat(iteration 3) = 0.876312568 Latitude = 0.876312568 rad = 50°12&#39;32.442&quot;N Longitude = 0.294083997 rad = 16°50&#39;59.179&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::1042<span class="nt">&lt;/gml:identifier&gt;</span> <span class="nt">&lt;gml:name&gt;</span>Krovak Modified<span class="nt">&lt;/gml:name&gt;</span> <span class="nt">&lt;gml:remarks&gt;</span>Incorporates a polynomial transformation which is defined to be exact and for practical purposes is considered to be a map projection.<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. From the defining parameters the following constants for the projection may be calculated : A = a(1 - e^2)^0.5 / [1 - e^2 sin^2.(latC)] B = {1 + [e^2 * cos^4(latC) / (1 - e^2)]}^0.5 gammao = asin[sin(latC) / B] to = tan(pi/4 + gammao/2).[(1 + e sin(latC)) / (1 - e sin(latC))]^(e.B/2) / [tan(pi/4 + latC/2)]^B n = sin(latp) ro = kp.A / tan(latp) To derive the projected Southing and Westing coordinates of a point with geographical coordinates (lat, lon) the formulas for the Krovak are: U = 2(atan{to.tan^B(lat/2 + pi/4) / [(1 + e sin(lat)) / (1 - e sin(lat))]^[e.B/2]} - pi/4) V = B(lonO - lon) where lonO and lon must both be referenced to the same prime meridian. T = asin[cos(alphaC).sin(U) + sin(alphaC).cos(U). cos(V)] D = asin[cos(U).sin(V)/cos(T)] theta = n.D r = ro.tan^n(pi/4 + latp/2) / tan^n(T/2 + pi/4) Xp = r.cos(theta) Yp = r.sin(theta) Xr = Xp – Xo Yr = Yp – Yo dX = C1 + C3.Xr – C4.Yr – 2.C6.Xr.Yr + C5.(Xr^2 – Yr^2) + C7.Xr.(Xr^2 – 3.Yr^2) – C8.Yr.(3.Xr^2 – Yr^2) + 4.C9.Xr.Yr.(Xr^2 – Yr^2) + C10.(Xr^4 + Yr^4 – 6.Xr^2.Yr^2) dY = C2 + C3.Yr + C4.Xr + 2.C5.Xr.Yr + C6.(Xr^2 – Yr^2) + C8.Xr.(Xr^2 – 3.Yr^2)+ C7.Yr.(3.Xr^2 – Yr^2) – 4.C10.Xr.Yr.(Xr^2 – Yr^2) + C9.(Xr^4 + Yr^4 – 6.Xr^2.Yr^2) Southing X = FN + Xp – dX Westing Y = FE + Yp – dY The reverse formulas to derive the latitude and longitude of a point from its Southing and Westing values are: Xr&#39; = (Southing – FN) – Xo Yr&#39; = (Westing – FE) – Yo dX&#39; = C1 + C3.Xr&#39; – C4.Yr&#39; – 2.C6.Xr&#39;.Yr&#39; + C5.(Xr&#39;^2 – Yr&#39;^2) + C7.Xr&#39;.(Xr&#39;^2 – 3.Yr&#39;^2) – C8.Yr&#39;.(3.Xr&#39;^2 – Yr&#39;^2) + 4.C9.Xr&#39;.Yr&#39;.(Xr&#39;^2 – Yr&#39;^2) + C10.(Xr&#39;^4 + Yr&#39;^4 – 6.Xr&#39;^2.Yr&#39;^2) dY&#39; = C2 + C3.Yr&#39; + C4.Xr&#39; + 2.C5.Xr&#39;.Yr&#39; + C6.(Xr&#39;^2 – Yr&#39;^2) + C8.Xr&#39;.(Xr&#39;^2 – 3.Yr&#39;^2) + C7.Yr&#39;.(3.Xr&#39;^2 – Yr&#39;^2) - 4.C10.Xr&#39;.Yr&#39;.(Xr&#39;^2 – Yr&#39;^2) + C9.(Xr&#39;^4 + Yr&#39;^4 – 6.Xr&#39;^2.Yr&#39;^2) Xp&#39; = (Southing – FN) + dX&#39; Yp&#39; = (Westing – FE) + dY&#39; r&#39; = [(Yp&#39;)^2 + (Xp&#39;)^2]^(1/2) theta&#39; = atan2[Yp&#39; , Xp&#39;] (see GN7-2 implementation notes in preface for atan2 convention) D&#39; = theta&#39; / sin(latp) T&#39; = 2{atan[((ro / r&#39;)^(1/n)).tan(pi/4 + latp/2)] - pi/4} U&#39; = asin[cos(alphaC).sin(T&#39;) - sin(alphaC).cos(T&#39;).cos(D&#39;)] V&#39; = asin(cos(T&#39;).sin(D&#39;) / cos(U&#39;)) Then latitude lat is found by iteration using U&#39; as the value for lat(j-1) in the first iteration: lat(j) = 2*(atan{tO^(-1/B) tan^(1/B).(U&#39;/2 + pi/4).[(1 + e sin(lat(j-1)) / (1 - e sin(lat(j-1))]^(e/2)} - pi/4) Then lon = lonO - V&#39; / B where lon is referenced to the same prime meridian as lonO.<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::8811&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::8833&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::1036&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::8818&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::8819&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:generalOperationParameter</span> <span class="na">xlink:href=</span><span class="s">&quot;urn:ogc:def:parameter:EPSG::8617&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::8618&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::1026&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::1027&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::1028&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::1029&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::1030&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::1031&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::1032&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::1033&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::1034&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::1035&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-1042">
  <gml:metaDataProperty>
    <epsg:CommonMetaData>
      <epsg:informationSource>Land Survey Office (ZU), Prague. www.cuzk.cz/zu</epsg:informationSource>
      <epsg:revisionDate>2018-08-29</epsg:revisionDate>
      <epsg:changes>
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2010.071" />
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2010.100" />
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2012.011" />
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2017.037" />
        <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::8811" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8833" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1036" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8818" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8819" />
      <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:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8617" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8618" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1026" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1027" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1028" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1029" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1030" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1031" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1032" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1033" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1034" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::1035" />
      <epsg:example>For Projected Coordinate Reference System: S-JTSK/05 (Ferro) / Modified Krovak

Parameters:
Ellipsoid  Bessel 1841   a = 6377397.155m  1/f = 299.15281
    then    e = 0.081696831        e^2 = 0.006674372

Latitude of projection centre = 49°30'00"N =  0.863937979 rad
Longitude of Origin = 42°30'00"E of Ferro = 0.741764932 rad
Co-latitude of cone axis = 30°17'17.30311" = 0.528627763 rad
Latitude of pseudo standard parallel = 78°30'00"N 1.370083463 rad
Scale factor on pseudo Standard Parallel (ko) = 0.9999
False Easting = 5000000.00 m
False Northing = 5000000.00 m
Ordinate 1 of evaluation point Xo = 1089000.00 m
Ordinate 2 of evaluation point Yo = 654000.00 m
C1 = 2.946529277E-02
C2 = 2.515965696E-02
C3 = 1.193845912E-07
C4 = -4.668270147E-07
C5 = 9.233980362E-12
C6 = 1.523735715E-12
C7 = 1.696780024E-18
C8 = 4.408314235E-18
C9 = -8.331083518E-24
C10 = -3.689471323E-24

Calculated projection constants:
A = 6380703.611
B = 1.000597498
gammao = 0.863239103
to = 1.003419164
n = 0.979924705
ro = 1298039.005

Forward calculation for: 
Latitude = 50°12'32.442"N = 0.876312568 rad
Longitude = 34°30'59.179"E of Ferro = 0.602425500 rad

U = 0.875596951
V = 0.139422687
S = 1.386275051
D = 0.506554627
theta = 0.496385393
r = 1194731.002
Xp = 1050538.631
Yp = 568990.995
Xr = -38461.369
Yr = -85009.005
dX = -0.077
dY = 0.088
and	
Southing X = 6050538.71 m
Westing	 Y = 5568990.91 m

Reverse calculation for the same Southing and Westing:

Xr' = -38461.292
Yr' = -85009.093
dX' = -0.077
dY' = 0.088
Xp' = 1050538.631
Yp' = 568990.995
r' = 1194731.002
theta' = 0.496385393
D' = 0.506554627
T' = 1.386275051
U' = 0.875596951
V' = 0.139422687
lat(iteration 1) = 0.876310603 
lat(iteration 2) = 0.876312562
lat(iteration 3) = 0.876312568

Latitude = 0.876312568 rad = 50°12'32.442"N
Longitude = 0.294083997 rad = 16°50'59.179"E</epsg:example>
    </epsg:CoordinateOperationMethodMetaData>
  </gml:metaDataProperty>
  <gml:identifier codeSpace="IOGP">urn:ogc:def:method:EPSG::1042</gml:identifier>
  <gml:name>Krovak Modified</gml:name>
  <gml:remarks>Incorporates a polynomial transformation which is defined to be exact and for practical purposes is considered to be a map projection.</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.

From the defining parameters the following constants for the projection may be calculated :

A = a(1 - e^2)^0.5 / [1 - e^2 sin^2.(latC)]
B = {1 + [e^2 * cos^4(latC) / (1 - e^2)]}^0.5
gammao = asin[sin(latC) / B]
to = tan(pi/4 + gammao/2).[(1 + e sin(latC)) / (1 - e sin(latC))]^(e.B/2) / [tan(pi/4 + latC/2)]^B
n = sin(latp)
ro = kp.A / tan(latp)

To derive the projected Southing and Westing coordinates of a point with geographical coordinates (lat, lon) the formulas for the Krovak are:

U = 2(atan{to.tan^B(lat/2 + pi/4) / [(1 + e sin(lat)) / (1 - e sin(lat))]^[e.B/2]} - pi/4)
V = B(lonO - lon) where lonO and lon must both be referenced to the same prime meridian.
T = asin[cos(alphaC).sin(U) + sin(alphaC).cos(U). cos(V)]
D = asin[cos(U).sin(V)/cos(T)]
theta = n.D
r = ro.tan^n(pi/4 + latp/2) / tan^n(T/2 + pi/4)
Xp = r.cos(theta)
Yp = r.sin(theta)
Xr = Xp – Xo
Yr = Yp – Yo
dX = C1 + C3.Xr – C4.Yr – 2.C6.Xr.Yr + C5.(Xr^2 – Yr^2) + C7.Xr.(Xr^2 – 3.Yr^2) – C8.Yr.(3.Xr^2 – Yr^2) + 4.C9.Xr.Yr.(Xr^2 – Yr^2) + C10.(Xr^4 + Yr^4 – 6.Xr^2.Yr^2)
dY = C2 + C3.Yr + C4.Xr + 2.C5.Xr.Yr + C6.(Xr^2 – Yr^2) + C8.Xr.(Xr^2 – 3.Yr^2)+ C7.Yr.(3.Xr^2 – Yr^2) – 4.C10.Xr.Yr.(Xr^2 – Yr^2) + C9.(Xr^4 + Yr^4 – 6.Xr^2.Yr^2)
Southing   X = FN + Xp – dX
Westing    Y = FE + Yp – dY

The reverse formulas to derive the latitude and longitude of a point from its Southing and Westing values are:

Xr' = (Southing – FN) – Xo
Yr' = (Westing – FE) – Yo
dX' = C1 + C3.Xr' – C4.Yr' – 2.C6.Xr'.Yr' + C5.(Xr'^2 – Yr'^2) + C7.Xr'.(Xr'^2 – 3.Yr'^2) – C8.Yr'.(3.Xr'^2 – Yr'^2) + 4.C9.Xr'.Yr'.(Xr'^2 – Yr'^2) + C10.(Xr'^4 + Yr'^4 – 6.Xr'^2.Yr'^2)
dY' = C2 + C3.Yr' + C4.Xr' + 2.C5.Xr'.Yr' + C6.(Xr'^2 – Yr'^2) + C8.Xr'.(Xr'^2 – 3.Yr'^2) 
+ C7.Yr'.(3.Xr'^2 – Yr'^2) - 4.C10.Xr'.Yr'.(Xr'^2 – Yr'^2) + C9.(Xr'^4 + Yr'^4 – 6.Xr'^2.Yr'^2)
Xp' = (Southing – FN) + dX'
Yp' = (Westing – FE) + dY'
r' = [(Yp')^2 + (Xp')^2]^(1/2)  
theta' = atan2[Yp' , Xp']   (see GN7-2 implementation notes in preface for atan2 convention)
D' = theta' / sin(latp)
T' = 2{atan[((ro / r')^(1/n)).tan(pi/4 + latp/2)] - pi/4}
U' = asin[cos(alphaC).sin(T') - sin(alphaC).cos(T').cos(D')]
V' = asin(cos(T').sin(D') / cos(U'))

Then latitude lat is found by iteration using U' as the value for lat(j-1) in the first iteration:
lat(j) = 2*(atan{tO^(-1/B) tan^(1/B).(U'/2 + pi/4).[(1 + e sin(lat(j-1)) / (1 - e sin(lat(j-1))]^(e/2)} - pi/4) 

Then
lon = lonO - V' / B where lon is referenced to the same prime meridian as lonO.</gml:formula>
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8811" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8833" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1036" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8818" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8819" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8806" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8807" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8617" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8618" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1026" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1027" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1028" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1029" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1030" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1031" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1032" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1033" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1034" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1035" />
</gml:OperationMethod>