EPSG:1070

Point motion by grid (Canada NTv2_Vel)

Attributes

Data source: OGP

Information source: Natural Resources Canada.

Revision date: 2018-08-29

Remarks: Interpolation within the grid is in the horizontal component of the source CRS.

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-1070&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>Natural Resources Canada.<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::2016.060&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::2018.034&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>false<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::1050&quot;</span> <span class="nt">/&gt;</span> <span class="nt">&lt;epsg:example&gt;</span>For geodetic station Centennial Monument NCC100 with ellipsoidal coordinates referenced to NAD83(CSRS)v6 at coordinate epoch 2010.00: lat = 45°25&#39;45.714920&quot;N lon = 75°42&#39;05.960075&quot;W h = 39.524 m t = 2010.00 (years) whose coordinates are required at coordinate epoch 2002.00. First interpolate within the grid file for velocities at lat = 45°25&#39;45.714920&quot;N, lon = 75°42&#39;05.960075&quot;W: VN = -0.00156 m/yr VE = 0.00177 m/yr V_up = 0.00202 m/yr NAD83(CSRS)v6 uses the GRS 1980 ellipsoid for which a = 6378137.0m and e^2 = 0.006694380. Then: V_lat = -2.44979E-10 rad/yr V_lon = 3.94759E-10 rad/yr and lat = 0.792892001 + (2002.00-2010.00) * -2.44979E-10 = 0.792892003 radians = 45°25&#39;45.715324&quot;N lon = -1.321243139 + (2002.00-2010.00) * 3.94759E-10 = -1.321243142 radians = 75°42&#39;05.960726&quot;W h = 39.524 + (2002.00-2010.00) * 0.00202 = 39.508m t = 2002.00<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::1070<span class="nt">&lt;/gml:identifier&gt;</span> <span class="nt">&lt;gml:name&gt;</span>Point motion by grid (Canada NTv2_Vel)<span class="nt">&lt;/gml:name&gt;</span> <span class="nt">&lt;gml:remarks&gt;</span>Interpolation within the grid is in the horizontal component of the source CRS.<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. The NTv2_Vel grid contains point velocities resolved into north, east, and up components. First the three velocities at a point are interpolated within the NTv2_Vel grid using bi-linear interpolation within the source and target CRS. Then ellipsoidal coordinates of a point at coordinate epoch t1 may be calculated at any other coordinate epoch t2 from: [Lat(t2)] = [Lat(t1)] + (t2-t1) * [V_lat] [Lon(t2)] = [Lon(t1)] + (t2-t1) * [V_lon] [h(t2)] = [h(t1)] + (t2-t1) * [V_up] where V describes the velocities of the ellipsoidal coordinates. These are given as a rate in linear units (metres per year) resolved into north, east and up components VN, VE and Vh. The north and east components are converted into latitude and longitude components by: V_lat = VN / (rho+h) V_lon = VE / [(nu+h) cos(lat)] where rho is radius of curvature of the CRS&#39;s ellipsoid in the plane of the meridian at latitude lat rho = a(1 – e^2)/(1 – e^2sin^2(lat))^3/2 nu is radius of curvature of the ellipsoid perpendicular to the meridian at latitude lat nu = a /(1 – e^2sin^2(lat))^1/2<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::1050&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-1070">
  <gml:metaDataProperty>
    <epsg:CommonMetaData>
      <epsg:informationSource>Natural Resources Canada.</epsg:informationSource>
      <epsg:revisionDate>2018-08-29</epsg:revisionDate>
      <epsg:changes>
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2016.060" />
        <epsg:changeID xlink:href="urn:ogc:def:change-request:EPSG::2018.034" />
      </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::1050" />
      <epsg:example>For geodetic station Centennial Monument NCC100 with ellipsoidal coordinates referenced to NAD83(CSRS)v6 at coordinate epoch 2010.00: 
lat = 45°25'45.714920"N
lon = 75°42'05.960075"W
h = 39.524 m
t = 2010.00 (years)				
whose coordinates are required at coordinate epoch 2002.00. 

First interpolate within the grid file for velocities at lat = 45°25'45.714920"N, lon = 75°42'05.960075"W:
VN  = -0.00156 m/yr
VE  =  0.00177 m/yr
V_up = 0.00202 m/yr

NAD83(CSRS)v6 uses the GRS 1980 ellipsoid for which a = 6378137.0m and e^2 = 0.006694380.

Then:
	V_lat = -2.44979E-10 rad/yr
	V_lon =  3.94759E-10 rad/yr
and
	lat =  0.792892001 + (2002.00-2010.00) * -2.44979E-10 =  0.792892003 radians = 45°25'45.715324"N
	lon = -1.321243139 + (2002.00-2010.00) *  3.94759E-10 = -1.321243142 radians = 75°42'05.960726"W
	h   =       39.524 + (2002.00-2010.00) * 0.00202 = 39.508m
	t = 2002.00</epsg:example>
    </epsg:CoordinateOperationMethodMetaData>
  </gml:metaDataProperty>
  <gml:identifier codeSpace="IOGP">urn:ogc:def:method:EPSG::1070</gml:identifier>
  <gml:name>Point motion by grid (Canada NTv2_Vel)</gml:name>
  <gml:remarks>Interpolation within the grid is in the horizontal component of the source CRS.</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 NTv2_Vel grid contains point velocities resolved into north, east, and up components.

First the three velocities at a point are interpolated within the NTv2_Vel grid using bi-linear interpolation within the source and target CRS.

Then ellipsoidal coordinates of a point at coordinate epoch t1 may be calculated at any other coordinate epoch t2 from:

[Lat(t2)] = [Lat(t1)] + (t2-t1) * [V_lat]
[Lon(t2)] = [Lon(t1)] + (t2-t1) * [V_lon]
[h(t2)] = [h(t1)] + (t2-t1) * [V_up]

where V describes the velocities of the ellipsoidal coordinates. These are given as a rate in linear units (metres per year) resolved into north, east and up components VN, VE and Vh. The north and east components are converted into latitude and longitude components by:
	V_lat = VN / (rho+h)
	V_lon = VE / [(nu+h) cos(lat)] 
where 
	rho is radius of curvature of the CRS's ellipsoid in the plane of the meridian at latitude lat 
	rho = a(1 – e^2)/(1 – e^2sin^2(lat))^3/2
	nu is radius of curvature of the ellipsoid perpendicular to the meridian at latitude lat 
	nu = a /(1 – e^2sin^2(lat))^1/2</gml:formula>
  <gml:generalOperationParameter xlink:href="urn:ogc:def:parameter:EPSG::1050" />
</gml:OperationMethod>