EPSG:5021

COORDINATEOPERATION["Porto Santo 1995 to PTRA08 (1)",
    VERSION["IGP-Prt MadArch"],
    SOURCECRS[
        GEOGCRS["Porto Santo 1995",
            DATUM["Porto Santo 1995",
                ELLIPSOID["International 1924",6378388,297,
                    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",4663]]],
    TARGETCRS[
        GEOGCRS["PTRA08",
            DATUM["Autonomous Regions of Portugal 2008",
                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",5013]]],
    METHOD["Geocentric translations (geog2D domain)",
        ID["EPSG",9603]],
    PARAMETER["X-axis translation",-503.229,
        LENGTHUNIT["metre",1],
        ID["EPSG",8605]],
    PARAMETER["Y-axis translation",-247.375,
        LENGTHUNIT["metre",1],
        ID["EPSG",8606]],
    PARAMETER["Z-axis translation",312.582,
        LENGTHUNIT["metre",1],
        ID["EPSG",8607]],
    OPERATIONACCURACY[2.0],
    USAGE[
        SCOPE["Transformation of coordinates at 2m level of accuracy."],
        AREA["Portugal - Madeira, Porto Santo and Desertas islands - onshore."],
        BBOX[32.35,-17.31,33.15,-16.23]],
    ID["EPSG",5021],
    REMARK["Derived at 34 points in May 2009. Residuals at 25 test points within 0.5m horizontal and 2m vertical. Info source also provides a Position Vector tfm of similar accuracy. See codes 5022 and 5023 for preferred tfms for islands of Porto Santo and Maderia."]]

proj4.defs("EPSG:5021","+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=push +v_3 +step +proj=cart +ellps=intl +step +proj=helmert +x=-503.229 +y=-247.375 +z=312.582 +step +inv +proj=cart +ellps=GRS80 +step +proj=pop +v_3 +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": "Porto Santo 1995 to PTRA08 (1)",
    "source_crs": {
        "type": "GeographicCRS",
        "name": "Porto Santo 1995",
        "datum": {
            "type": "GeodeticReferenceFrame",
            "name": "Porto Santo 1995",
            "ellipsoid": {
                "name": "International 1924",
                "semi_major_axis": 6378388,
                "inverse_flattening": 297
            }
        },
        "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": 4663
        }
    },
    "target_crs": {
        "type": "GeographicCRS",
        "name": "PTRA08",
        "datum": {
            "type": "GeodeticReferenceFrame",
            "name": "Autonomous Regions of Portugal 2008",
            "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": 5013
        }
    },
    "method": {
        "name": "Geocentric translations (geog2D domain)",
        "id": {
            "authority": "EPSG",
            "code": 9603
        }
    },
    "parameters": [
        {
            "name": "X-axis translation",
            "value": -503.229,
            "unit": "metre",
            "id": {
                "authority": "EPSG",
                "code": 8605
            }
        },
        {
            "name": "Y-axis translation",
            "value": -247.375,
            "unit": "metre",
            "id": {
                "authority": "EPSG",
                "code": 8606
            }
        },
        {
            "name": "Z-axis translation",
            "value": 312.582,
            "unit": "metre",
            "id": {
                "authority": "EPSG",
                "code": 8607
            }
        }
    ],
    "accuracy": "2.0",
    "scope": "Transformation of coordinates at 2m level of accuracy.",
    "area": "Portugal - Madeira, Porto Santo and Desertas islands - onshore.",
    "bbox": {
        "south_latitude": 32.35,
        "west_longitude": -17.31,
        "north_latitude": 33.15,
        "east_longitude": -16.23
    },
    "id": {
        "authority": "EPSG",
        "code": 5021
    },
    "remarks": "Derived at 34 points in May 2009. Residuals at 25 test points within 0.5m horizontal and 2m vertical. Info source also provides a Position Vector tfm of similar accuracy. See codes 5022 and 5023 for preferred tfms for islands of Porto Santo and Maderia."
}

<?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=push +v_3 +step +proj=cart +ellps=intl +step +proj=helmert +x=-503.229 +y=-247.375 +z=312.582 +step +inv +proj=cart +ellps=GRS80 +step +proj=pop +v_3 +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=push +v_3 +step +proj=cart +ellps=intl +step +proj=helmert +x=-503.229 +y=-247.375 +z=312.582 +step +inv +proj=cart +ellps=GRS80 +step +proj=pop +v_3 +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 ( 5021, 'EPSG', 5021, '+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=push +v_3 +step +proj=cart +ellps=intl +step +proj=helmert +x=-503.229 +y=-247.375 +z=312.582 +step +inv +proj=cart +ellps=GRS80 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1', 'None');