EPSG:1122

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.

This method transforms coordinates between a geographic 3D CRS and a compound CRS consisting of separate geographic 2D and vertical CRSs. The transformation is three-dimensional and reversible.

A transformation using this operation method consists of two parts:
a)	The transformation of coordinates from the geographic 3D to or from the geographic 2D CRS;
b)	The transformation of ellipsoidal height to or from depth.

However, complexities arise if the source 3D and horizontal component of the target compound CRS are based on different geodetic datums. To account for any offset between and non-parallelism of the ellipsoid surfaces, the method is entirely dependent upon a suitable three-dimensional transformation existing between the geographic 3D and geographic 2D CRSs, and this may not always be the case. For this reason the general case is considered to be out of scope here and this method is restricted to circumstances where the source geographic 3D CRS and the target geographic 2D CRS are both based on the same geodetic datum, i.e. the horizontal component of the compound CRS is the two-dimensional subset of the geographic 3D CRS. This geographic 2D CRS is also the CRS used for interpolation of the geoid model file grid.

Forward transformation from geographic 3D to compound
Because of the restriction above, latitude and longitude of the point in the source CRS do not change as a result of the transformation; therefore φt = φs and λt = λs, where the subscripts 's' and ‘t’ indicates the coordinates in the source and target CRSs respectively.

The vertical component of the forward transformation is as for the Geographic3D to Depth (gtx) method, code 1109. The method requires that the ‘Interpolation CRS’ - the horizontal CRS to which the grid containing the hydroid heights (ζ) are referenced - to be specified. That will be the geographic 2D CRS which is both the horizontal subset of the geographic 3D CRS and the horizontal component of the compound CRS. The unit of measure in which the geoid heights are given is part of the file format definition.

The forward case is straightforward:
	φt  =  φs
	λt  =  λs
Then using (φ, λ) for interpolation for the geoid height at the point,
	D  =  ζ – h
where D is depth below tidal surface, ζ is height offset from ellipsoid to tidal surface and h is ellipsoid height.

In a marine context,
	D = (Dobs – hobs) + ζ

Reverse transformation from compound to geographic 3D
For the reverse case the height offset ζ is first interpolated using the geographic 2D coordinates (φ, λ).
Then:
	φt  =  φs
	λt  =  λs
	h  =  H  +  ζ
	h  =  ζ  +  (Dobs –  D)
where h and ζ are heights and D is depth.