EPSG Dataset :: v12.003

Vertical Offset by Grid Interpolation (asc)

Coordinate Operation Method Details [VALID]

Name:

Vertical Offset by Grid Interpolation (asc)

Code:

1085

Operation is Reversible:

Yes

Formula:

The relationship between two gravity-related coordinate reference systems may be available through gridded data files of offsets (sometimes called height differences).

The vertical offset at a point is first interpolated within the grid of values using bi-linear interpolation. For interpolation within the offset grid, horizontal positions in the Interpolation CRS are used.

The interpolated value is then applied as a correction parameter A. The interpolated offset is added to the height in the source CRS:
H2 = H1 + A1>2
where
H2 = value in the target vertical coordinate reference system.
H1 = value in the source vertical coordinate reference system.
A1>2 is the offset to be applied for the transformation from CRS 1 to CRS 2. Its value for the forward calculation is the value of the origin of the source CRS 1 in the target CRS 2.

For the reverse transformation from CRS 2 to CRS 1, the same formula is used but with the sign of the offset A1>2 reversed:
H1 = H2 + (–A1>2)

Care! In the EPSG Dataset, offsets are additive (for the forward transformation). See EPSG Guidance Note 7-2, Implementation Notes - Offsets.

Example:

For coordinate transformation: NZVD2016 height to Dunedin 1958 height (2), code 9314:

For a point at 44.42°S, 168.92°E (using CRS code 4167, NZGD2000, for interpolation) with NZVD2016 height of 50.000 metres, to find its Dunedin 1958 height:

First obtain the offsets at each of the surrounding grid nodes:
NW corner 44.4000°S, 168.9000°E, offset = 0.3028m
NE corner 44.4000°S, 168.9333°E, offset = 0.3029m
SE corner 44.4333°S, 168.9333°E, offset = 0.3055m
SW corner 44.4333°S, 168.9000°E, offset = 0.3051m

Then using bi-linear interpolation for 44.42°S, 168.92°E, offset = 0.304m

Then Dunedin 1958 height = 50.000 + 0.304 = 50.304m.

For the reverse calculation to find the NZVD2016 height from a point with Dunedin 1958 height of 50.304m:

First obtain the offsets at each of the surrounding grid nodes again using coordinates referenced to NZGD2000 for interpolation:
NW corner 44.4000°S, 168.9000°E, offset = 0.3028m
NE corner 44.4000°S, 168.9333°E, offset = 0.3029m
SE corner 44.4333°S, 168.9333°E, offset = 0.3055m
SW corner 44.4333°S, 168.9000°E, offset = 0.3051m

Then using bi-linear interpolation for 44.42°S, 168.92°E, offset = 0.304m

Then NZVD2016 height = 50.304 - 0.304 = 50.000m.

Method
Parameters:

Parameter Name	Parameter Code	Sign reversal	Parameter Description
Vertical offset file	8732	Yes	The name of the [path and] file containing vertical offsets, the differences in gravity-related height.
EPSG code for Interpolation CRS	1048	No	The EPSG code for the CRS that should be used to interpolate gridded data.

Meta Data

Information Source:

IOGP

Data Source:

EPSG

Revision Date:

July 4, 2020

Change ID:

[2020.057]

Wildcard	Example
% or *	Any wildcard characters
? or _	One wildcard character
[ - ]	[ - ] Range e.g. N[a-c] will return all geodetic objects that contains a string starting with the letter 'N' or 'n' followed by any of the letters contained within the range specified within the square bracket, e.g. letters a, b and c.
\|	Or e.g. E[a\|e]M will return all geodetic objects that contains a string starting with the letter E(e), ends with the letter M(m) and also contain either the letter a or letter e.