EPSG:1071

Vertical Offset by Grid Interpolation (NZLVD)

Attributes

Data source: EPSG

Information source: IOGP

Revision date: 2019-03-22

Formula

The relationship between NZVD16 height and thirteen local gravity-related coordinate reference systems (LVDs) for New Zealand is 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 grids, horizontal positions in the NZGD2000 geographic 2D CRS are used.

The interpolated value is then applied as an offset:
•	If the height in a LVD is required when a NZVD2016 height is given, add the interpolated offset to the NZVD2016 height (LVD height = NZVD2016 height + interpolated offset).
•	If the NZVD2016 height is required when a LVD height is given, subtract the interpolated offset from the LVD height (NZVD2016 height = LVD height - interpolated offset).

Care! In the EPSG Dataset, offsets are additive (for the forward transformation). See EPSG Guidance Note 7-2, Implementation Notes - Offsets. The LINZ convention assumes that the offset is subtracted, and grid file names reflect that convention. The offset files therefore have names which conflict with the source and target CRSs of the EPSG transformation.

Example

For a point in the Dunedin 1958 height area at 44.42°S, 168.92°E (NZGD2000) 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:
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.
MapTiler banner