EPSG:1058

Time-dependent Coordinate Frame rotation (geog3D)

Attributes

Data source: EPSG

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

Revision date: 2019-02-21

Remarks: Note the analogy with the Time-dependent Position Vector transformation (code 1055) but beware of the differences! The Position Vector convention is used by IAG. See method codes 1056 and 1057 for similar methods operating between other CRS types.

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.

Transformation of coordinates from one geographic coordinate reference system into another is often carried out as a concatenation of the following operations:

  (geographical to geocentric) + (geocentric to geocentric) + (geocentric to geographic)

The Time-dependent Coordinate Frame rotation (geog3D domain) transformation has 3 steps:

(i) geographic 3D coordinates are converted to geocentric coordinates using EPSG coordinate operation method code 9602;

(ii) the middle step of the concatenated transformation, from geocentric coordinates to geocentric coordinates, uses the Time-dependent Coordinate Frame rotation(geocentric domain) method, EPSG method code 1056;

(iii) finally the geocentric coordinates are converted to geographic 3D using EPSG coordinate operation method code 9602.

Example

Input point: 
Coordinate reference system: WGS 84 (G1674) = ITRF08 (geographic 2D)
  Latitude  =  15 deg 28 min 32.368 sec S
  Longitude =  128 deg 02 min 56.198 sec E
  Ellipsoidal height hs  =  46.244 m
at epoch 2013.90

This transforms to WGS 84 (G1674) = ITRF2008 Cartesian geocentric coordinates:
Xs = -3789470.702 m
Ys =  4841770.411 m
Zs = -1690893.950 m
(at epoch 2013.9)

Transformation parameters ITRF08 to GDA94:
tX = –84.68 mm
dtX = +1.42 mm/yr
tY = –19.42 mm
dtY = +1.34 mm/yr
tZ = +32.01 mm
dtZ = +0.90 mm/yr
rX = -0.4254 msec
drX = +1.5461 msec/yr
rY = +2.2578 msec
drY = +1.1820 msec/yr
rZ = +2.4015 msec
drZ = +1.1551 msec/yr
dS = +0.00971 ppm
ddS = +0.000109 ppm/yr
t0 = 1994.00				

Corrections due to rate of change to each of the 7 transformation parameters for the period (t-t0), taking care to convert the translations to the same units as the source CRS (in this case metres) and the rotations to radians, are:
tX' = –84.68 + [(1.42 * (2013.90-1994.00)] = –56.42 mm = –0.056m
tY' = +0.007 m
tZ' = +0.050 m
rX' = +1.471021E-07 rad
rY' = +1.249830E-07 rad
rZ' = +1.230844E-07 rad
dS' = +0.01188 ppm

from which M = 1.00000001188

Using these parameter values, application of the 7 parameter Coordinate Frame rotation formula results in:
Xt = –3789470.008 m
Yt =  4841770.685 m
Zt = –1690895.103 m
on the GDA94 geocentric coordinate reference system.

   This converts into:
   Latitude  =  15 deg 28 min 32.406 sec S
   Longitude =  128 deg 02 min 56.174 sec E
   Ellipsoidal height =  +46.244m
on the GDA94 geographic 3D coordinate reference system.
MapTiler banner