EPSG:7871

COORDINATEOPERATION["NZVD2016 height to Taranaki 1970 height (1)",
    VERSION["LINZ-NZ TARA"],
    SOURCECRS[
        VERTCRS["NZVD2016 height",
            VDATUM["New Zealand Vertical Datum 2016"],
            CS[vertical,1],
                AXIS["gravity-related height (H)",up,
                    LENGTHUNIT["metre",1]],
            ID["EPSG",7839]]],
    TARGETCRS[
        VERTCRS["Taranaki 1970 height",
            VDATUM["Taranaki 1970"],
            CS[vertical,1],
                AXIS["gravity-related height (H)",up,
                    LENGTHUNIT["metre",1]],
            ID["EPSG",5769]]],
    METHOD["Vertical Offset by Grid Interpolation (NZLVD)",
        ID["EPSG",1071]],
    PARAMETERFILE["Vertical offset file","taranaki-1970-to-nzvd2016-conversion.csv"],
    INTERPOLATIONCRS[
        GEOGCRS["NZGD2000",
            DATUM["New Zealand Geodetic Datum 2000",
                ELLIPSOID["GRS 1980",6378137,298.257222101,
                    LENGTHUNIT["metre",1]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433]],
            CS[ellipsoidal,2],
                AXIS["geodetic latitude (Lat)",north,
                    ORDER[1],
                    ANGLEUNIT["degree",0.0174532925199433]],
                AXIS["geodetic longitude (Lon)",east,
                    ORDER[2],
                    ANGLEUNIT["degree",0.0174532925199433]]]],
    OPERATIONACCURACY[0.02],
    USAGE[
        SCOPE["Change of height to a different vertical reference surface."],
        AREA["New Zealand - North Island - Taranaki vertical CRS area."],
        BBOX[-39.92,173.68,-38.41,174.95]],
    ID["EPSG",7871],
    REMARK["Derived at 125 control points. Mean offset 0.286m, standard deviation 0.026m, maximum difference from mean 0.07m."]]

proj4.defs("EPSG:7871","+proj=vgridshift +grids=nz_linz_taraht1970-nzvd2016.tif +multiplier=1");

{
    "$schema": "https://proj.org/schemas/v0.7/projjson.schema.json",
    "type": "Transformation",
    "name": "NZVD2016 height to Taranaki 1970 height (1)",
    "source_crs": {
        "type": "VerticalCRS",
        "name": "NZVD2016 height",
        "datum": {
            "type": "VerticalReferenceFrame",
            "name": "New Zealand Vertical Datum 2016"
        },
        "coordinate_system": {
            "subtype": "vertical",
            "axis": [
                {
                    "name": "Gravity-related height",
                    "abbreviation": "H",
                    "direction": "up",
                    "unit": "metre"
                }
            ]
        },
        "id": {
            "authority": "EPSG",
            "code": 7839
        }
    },
    "target_crs": {
        "type": "VerticalCRS",
        "name": "Taranaki 1970 height",
        "datum": {
            "type": "VerticalReferenceFrame",
            "name": "Taranaki 1970"
        },
        "coordinate_system": {
            "subtype": "vertical",
            "axis": [
                {
                    "name": "Gravity-related height",
                    "abbreviation": "H",
                    "direction": "up",
                    "unit": "metre"
                }
            ]
        },
        "id": {
            "authority": "EPSG",
            "code": 5769
        }
    },
    "interpolation_crs": {
        "type": "GeographicCRS",
        "name": "NZGD2000",
        "datum": {
            "type": "GeodeticReferenceFrame",
            "name": "New Zealand Geodetic Datum 2000",
            "ellipsoid": {
                "name": "GRS 1980",
                "semi_major_axis": 6378137,
                "inverse_flattening": 298.257222101
            }
        },
        "coordinate_system": {
            "subtype": "ellipsoidal",
            "axis": [
                {
                    "name": "Geodetic latitude",
                    "abbreviation": "Lat",
                    "direction": "north",
                    "unit": "degree"
                },
                {
                    "name": "Geodetic longitude",
                    "abbreviation": "Lon",
                    "direction": "east",
                    "unit": "degree"
                }
            ]
        },
        "id": {
            "authority": "EPSG",
            "code": 4167
        }
    },
    "method": {
        "name": "Vertical Offset by Grid Interpolation (NZLVD)",
        "id": {
            "authority": "EPSG",
            "code": 1071
        }
    },
    "parameters": [
        {
            "name": "Vertical offset file",
            "value": "taranaki-1970-to-nzvd2016-conversion.csv",
            "id": {
                "authority": "EPSG",
                "code": 8732
            }
        }
    ],
    "accuracy": "0.02",
    "scope": "Change of height to a different vertical reference surface.",
    "area": "New Zealand - North Island - Taranaki vertical CRS area.",
    "bbox": {
        "south_latitude": -39.92,
        "west_longitude": 173.68,
        "north_latitude": -38.41,
        "east_longitude": 174.95
    },
    "id": {
        "authority": "EPSG",
        "code": 7871
    },
    "remarks": "Derived at 125 control points. Mean offset 0.286m, standard deviation 0.026m, maximum difference from mean 0.07m."
}

<?xml version="1.0" encoding="utf-8"?>
<Map srs="+proj=vgridshift +grids=nz_linz_taraht1970-nzvd2016.tif +multiplier=1">
	<Layer srs="+proj=vgridshift +grids=nz_linz_taraht1970-nzvd2016.tif +multiplier=1">
	</Layer>
</Map>

INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) values ( 7871, 'EPSG', 7871, '+proj=vgridshift +grids=nz_linz_taraht1970-nzvd2016.tif +multiplier=1', 'None');