EPSG:1087

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.

In France the national mapping agency (IGN) has promolgated a transformation between the old geographic 2D coordinate reference system NTF and the modern 3-dimensional system RGF93 which uses geocentric translations interpolated from a grid file. This is a two-step operation which first requires interpolation of a gridded dataset in the geographic 2D domain to obtain the geocentric translations, then the transformation is made using the Geocentric translations method. It is described in IGN document NTG-88. The method has also been used in New Caledonia. In summary:
•	The grid file nodes are referenced to the modern (RGF93 in France) geographic 2D CRS (the Interpolation CRS).
•	Within the grid file the sense of the parameter values is from old to modern (for France, NTF to RGF93), that is the offsets are additive when the old CRS (in France, NTF) is the source CRS and the modern CRS (in France, RGF93) is the target CRS.

For the forward transformation, old CRS to modern CRS, because the interpolation CRS is not the source CRS, in principle iteration is required. The source CRS coordinates are used for the initial bi-linear interpolation of the grid to obtain preliminary geocentric translation offsets. These are applied to calculate provisional target CRS coordinates:
  X't = Xs + tX
  Y't = Ys + tY
  Z't = Zs + tZ
This process is iterated until the changes in target coordinates are insignificant. 

An approximation to an accuracy of 1cm may be made if a suitable 'standard transformation' is available. This avoids iteration. The steps for the forward transformation then are:
•	Convert source CRS coordinates to interpolation CRS coordinates by using the standard transformation.	
•	Using these, interpolate within the grid file to obtain the three geocentric translations (tX, tY, tZ) applicable at the point.
•	Using the Geocentric translation (geog 2D domain) method,  transform coordinates of the point in the source CRS to the target CRS:
  Xt = Xs + tX
  Yt = Ys + tY
  Zt = Zs + tZ

For the reverse transformation from modern CRS (in France, RGF93) to old CRS (in France, NTF) where the modern CRS is the Interpolation CRS, the steps are:
•	Using the modern (RGF93) geographic coordinates, interpolate within the grid file to obtain the three geocentric translations (tX, tY, tZ) applicable at the point.
•	Transform modern CRS (RGF93) coordinates to old CRS (NTF) coordinates, taking into account the sense of the geocentric translation parameter values: 
  X(old) = X(modern)+(– tX) which for France may be expressed as X(NTF) = X(RGF93)+(– tX)
  Y(old) = Y(modern)+(– tY) which for France may be expressed as Y(NTF) = Y(RGF93)+(– tY)
  Z(old) = Z(modern)+(– tZ) which for France may be expressed as Z(NTF) = Z(RGF93)+(– tZ)