EPSG:9887

COORDINATEOPERATION["NAD83 to NAD83(CSRS)v2 (2)",
    VERSION["SK PMC-Can SK"],
    SOURCECRS[
        GEOGCRS["NAD83",
            DATUM["North American Datum 1983",
                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]],
            ID["EPSG",4269]]],
    TARGETCRS[
        GEOGCRS["NAD83(CSRS)v2",
            DATUM["North American Datum of 1983 (CSRS) version 2",
                ELLIPSOID["GRS 1980",6378137,298.257222101,
                    LENGTHUNIT["metre",1]],
                ANCHOREPOCH[1997]],
            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]],
            ID["EPSG",8237]]],
    METHOD["NTv2",
        ID["EPSG",9615]],
    PARAMETERFILE["Latitude and longitude difference file","SK83-98.gsb"],
    OPERATIONACCURACY[1.5],
    USAGE[
        SCOPE["Transformation of coordinates at 1m to 2m level of accuracy."],
        AREA["Canada - Saskatchewan."],
        BBOX[49,-110,60.01,-101.34]],
    ID["EPSG",9887],
    REMARK["Can be taken as an approximate transformation NAD83 to WGS 84 - see code 1697."]]

proj4.defs("EPSG:9887","+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=hgridshift +grids=ca_nrc_SK83-98.tif +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1");

{
    "$schema": "https://proj.org/schemas/v0.7/projjson.schema.json",
    "type": "Transformation",
    "name": "NAD83 to NAD83(CSRS)v2 (2)",
    "source_crs": {
        "type": "GeographicCRS",
        "name": "NAD83",
        "datum": {
            "type": "GeodeticReferenceFrame",
            "name": "North American Datum 1983",
            "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": 4269
        },
        "remarks": "Longitude is POSITIVE EAST. The adjustment included connections to Greenland and Mexico but the system was not adopted there. For applications with an accuracy of better than 1m replaced by NAD83(HARN) in the US and PRVI and by NAD83(CSRS) in Canada."
    },
    "target_crs": {
        "type": "GeographicCRS",
        "name": "NAD83(CSRS)v2",
        "datum": {
            "type": "GeodeticReferenceFrame",
            "name": "North American Datum of 1983 (CSRS) version 2",
            "anchor_epoch": 1997,
            "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": 8237
        },
        "remarks": "Adopted by the Canadian federal government from 1998-01-01 and by the provincial governments of British Columbia, New Brunswick, Prince Edward Island and Quebec. Replaces NAD83(CSRS96). Replaced by NAD83(CSRS)v3 (code 8240). Longitudes are POSITIVE EAST."
    },
    "method": {
        "name": "NTv2",
        "id": {
            "authority": "EPSG",
            "code": 9615
        }
    },
    "parameters": [
        {
            "name": "Latitude and longitude difference file",
            "value": "SK83-98.gsb",
            "id": {
                "authority": "EPSG",
                "code": 8656
            }
        }
    ],
    "accuracy": "1.5",
    "scope": "Transformation of coordinates at 1m to 2m level of accuracy.",
    "area": "Canada - Saskatchewan.",
    "bbox": {
        "south_latitude": 49,
        "west_longitude": -110,
        "north_latitude": 60.01,
        "east_longitude": -101.34
    },
    "id": {
        "authority": "EPSG",
        "code": 9887
    },
    "remarks": "Can be taken as an approximate transformation NAD83 to WGS 84 - see code 1697."
}

<?xml version="1.0" encoding="utf-8"?>
<Map srs="+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=hgridshift +grids=ca_nrc_SK83-98.tif +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1">
	<Layer srs="+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=hgridshift +grids=ca_nrc_SK83-98.tif +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1">
	</Layer>
</Map>

INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) values ( 9887, 'EPSG', 9887, '+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=hgridshift +grids=ca_nrc_SK83-98.tif +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1', 'None');