EPSG:1501

COORDINATEOPERATION["NAD83 to NAD83(HARN) (28)",
    VERSION["NGS-Usa OR WA"],
    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(HARN)",
            DATUM["NAD83 (High Accuracy Reference Network)",
                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",4152]]],
    METHOD["NADCON",
        ID["EPSG",9613]],
    PARAMETERFILE["Latitude difference file","wohpgn.las"],
    PARAMETERFILE["Longitude difference file","wohpgn.los"],
    OPERATIONACCURACY[0.05],
    USAGE[
        SCOPE["Geodesy."],
        AREA["United States (USA) - Oregon and Washington."],
        BBOX[41.98,-124.79,49.05,-116.47]],
    ID["EPSG",1501],
    REMARK["Accuracy 67% confidence level. Uses NADCON method which expects longitudes positive west; EPSG GeogCRSs NAD83 (code 4269) and NAD83(HARN) (code 4152) have longitudes positive east. May be taken as approximate transformation NAD83-WGS 84 - see code 1737."]]

proj4.defs("EPSG:1501","+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=hgridshift +grids=us_noaa_WO.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(HARN) (28)",
    "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(HARN)",
        "datum": {
            "type": "GeodeticReferenceFrame",
            "name": "NAD83 (High Accuracy Reference Network)",
            "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": 4152
        },
        "remarks": "In CONUS, AK, HI and PRVI replaces NAD83 for applications with an accuracy of better than 1m. Replaced by NAD83(FBN) in CONUS, American Samoa and Guam / NMI, by NAD83(NSRS2007) in Alaska, by NAD83(PA11) in Hawaii and by NAD83(HARN Corrected) in PRVI."
    },
    "method": {
        "name": "NADCON",
        "id": {
            "authority": "EPSG",
            "code": 9613
        }
    },
    "parameters": [
        {
            "name": "Latitude difference file",
            "value": "wohpgn.las",
            "id": {
                "authority": "EPSG",
                "code": 8657
            }
        },
        {
            "name": "Longitude difference file",
            "value": "wohpgn.los",
            "id": {
                "authority": "EPSG",
                "code": 8658
            }
        }
    ],
    "accuracy": "0.05",
    "scope": "Geodesy.",
    "area": "United States (USA) - Oregon and Washington.",
    "bbox": {
        "south_latitude": 41.98,
        "west_longitude": -124.79,
        "north_latitude": 49.05,
        "east_longitude": -116.47
    },
    "id": {
        "authority": "EPSG",
        "code": 1501
    },
    "remarks": "Accuracy 67% confidence level. Uses NADCON method which expects longitudes positive west; EPSG GeogCRSs NAD83 (code 4269) and NAD83(HARN) (code 4152) have longitudes positive east. May be taken as approximate transformation NAD83-WGS 84 - see code 1737."
}

<?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=us_noaa_WO.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=us_noaa_WO.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 ( 1501, 'EPSG', 1501, '+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=hgridshift +grids=us_noaa_WO.tif +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1', 'None');