EPSG:9653

Complex polynomial of degree 4

Attributes

Data source: OGP

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

Revision date: 2017-06-13

Remarks: Coordinate pairs treated as complex numbers. This exploits the correlation between the polynomial coefficients and leads to a smaller number of coefficients than the general polynomial of degree 4.

MapTiler banner

Export

Definition: OGP XML

<?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-9653">
  <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>false</epsg:isOperationReversible>
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8619" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8620" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8621" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8622" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8694" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8695" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8624" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8625" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8626" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8627" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8628" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8629" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8630" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8631" />
      <epsg:example>For transformation RD / Netherlands New to ED50 / UTM zone 31N,

Eo1 = 155000
No1 = 463000
Eo2 = 663395.607
No2 = 5781194.380
A1 = -51.681
A2 = 3290.525
A3 = 20.172
A4 = 1.133
A5 = 2.075
A6 = 0.251
A7 = 0.075
A8 = -0.012

For source coordinate system E1=200000  N1=500000, then
E2 =707155.557  N2 = 5819663.128.</epsg:example>
    </epsg:CoordinateOperationMethodMetaData>
  </gml:metaDataProperty>
  <gml:identifier codeSpace="IOGP">urn:ogc:def:method:EPSG::9653</gml:identifier>
  <gml:name>Complex polynomial of degree 4</gml:name>
  <gml:remarks>Coordinate pairs treated as complex numbers.  This exploits the correlation between the polynomial coefficients and leads to a smaller number of coefficients than the general polynomial of degree 4.</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.

The relationship between two projected coordinate reference systems may be approximated more elegantly by a single polynomial regression formula written in terms of complex numbers. The advantage is that the dependence between the ‘A’ and ‘B’ coefficients (for U and V) is taken into account in the formula, resulting in fewer coefficients for the same order polynomial. A fourth-order polynomial in complex numbers is used in The Netherlands for transforming coordinates referenced to the Amersfoort / RD system to and from ED50 / UTM.

mT.(dX + i. dY) = (A1 + i. A2).(U + i.V) + (A3 + i. A4).(U + i.V)^2 + (A5 + i. A6).(U + i.V)^3 + (A7 + i.A8).(U + i.V)^4

where U = mS.(XS - XS0)
           V = mS.(YS - YS0)
and mS, mT are the scaling factors for the coordinate differences in the source and target coordinate reference systems.

The polynomial to degree 4 can alternatively be expressed in matrix form.

Then
XT  = XS - XS0 + XT0 + dX
YT  = YS - YS0 + YT0 + dY

where
XT , YT      are coordinates in the target coordinate reference system,
XS , YS      are coordinates in the source coordinate reference system,
XS0 , YS0   are coordinates of the evaluation point in the source coordinate reference system,
XT0 , YT0   are coordinates of the evaluation point in the target coordinate reference system.

Note that the zero order coefficients of the general polynomial, A0 and B0, have apparently disappeared.  In reality they are absorbed by the different coordinates of the source and of the target evaluation point, which in this case, are numerically very different because of the use of two different projected coordinate reference systems for source and target.

The transformation parameter values (the coefficients) are not reversible.  For the reverse transformation a different set of parameter values are required, used within the same formulas as the forward direction.</gml:formula>
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8619" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8620" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8621" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8622" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8694" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8695" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8624" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8625" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8626" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8627" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8628" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8629" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8630" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8631" />
</gml:OperationMethod>
<?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-9653">
  <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>false</epsg:isOperationReversible>
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8619" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8620" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8621" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8622" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8694" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8695" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8624" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8625" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8626" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8627" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8628" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8629" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8630" />
      <epsg:signReversal changeSign="false" xlink:href="urn:ogc:def:parameter:EPSG::8631" />
      <epsg:example>For transformation RD / Netherlands New to ED50 / UTM zone 31N,

Eo1 = 155000
No1 = 463000
Eo2 = 663395.607
No2 = 5781194.380
A1 = -51.681
A2 = 3290.525
A3 = 20.172
A4 = 1.133
A5 = 2.075
A6 = 0.251
A7 = 0.075
A8 = -0.012

For source coordinate system E1=200000  N1=500000, then
E2 =707155.557  N2 = 5819663.128.</epsg:example>
    </epsg:CoordinateOperationMethodMetaData>
  </gml:metaDataProperty>
  <gml:identifier codeSpace="IOGP">urn:ogc:def:method:EPSG::9653</gml:identifier>
  <gml:name>Complex polynomial of degree 4</gml:name>
  <gml:remarks>Coordinate pairs treated as complex numbers.  This exploits the correlation between the polynomial coefficients and leads to a smaller number of coefficients than the general polynomial of degree 4.</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.

The relationship between two projected coordinate reference systems may be approximated more elegantly by a single polynomial regression formula written in terms of complex numbers. The advantage is that the dependence between the ‘A’ and ‘B’ coefficients (for U and V) is taken into account in the formula, resulting in fewer coefficients for the same order polynomial. A fourth-order polynomial in complex numbers is used in The Netherlands for transforming coordinates referenced to the Amersfoort / RD system to and from ED50 / UTM.

mT.(dX + i. dY) = (A1 + i. A2).(U + i.V) + (A3 + i. A4).(U + i.V)^2 + (A5 + i. A6).(U + i.V)^3 + (A7 + i.A8).(U + i.V)^4

where U = mS.(XS - XS0)
           V = mS.(YS - YS0)
and mS, mT are the scaling factors for the coordinate differences in the source and target coordinate reference systems.

The polynomial to degree 4 can alternatively be expressed in matrix form.

Then
XT  = XS - XS0 + XT0 + dX
YT  = YS - YS0 + YT0 + dY

where
XT , YT      are coordinates in the target coordinate reference system,
XS , YS      are coordinates in the source coordinate reference system,
XS0 , YS0   are coordinates of the evaluation point in the source coordinate reference system,
XT0 , YT0   are coordinates of the evaluation point in the target coordinate reference system.

Note that the zero order coefficients of the general polynomial, A0 and B0, have apparently disappeared.  In reality they are absorbed by the different coordinates of the source and of the target evaluation point, which in this case, are numerically very different because of the use of two different projected coordinate reference systems for source and target.

The transformation parameter values (the coefficients) are not reversible.  For the reverse transformation a different set of parameter values are required, used within the same formulas as the forward direction.</gml:formula>
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8619" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8620" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8621" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8622" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8694" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8695" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8624" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8625" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8626" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8627" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8628" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8629" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8630" />
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::8631" />
</gml:OperationMethod>