Reseau Geodesique de la Polynesie Francaise
EPSG
Gouvernement de la Polynésie Française, Service de l'Urbanisme, Section topographie. Also www.shom.fr, technical report 002/199.
2020-03-14
Geodesy.
SHOM report gives scale difference as 0.999 999 9907 (wrt unity). Transformation is to original Transit definition of WGS 84. It is consistent with later WGS 84 realisations G730, G873 and G1150 to no better than 1m.
Position Vector transformation (geog2D domain)
(0.072, -0.507, -0.245, 0.0183, -0.0003, 0.007, -0.0093)
Definition: OGC Well Known Text 2 (2019)
COORDINATEOPERATION["RGPF to WGS 84 (1)", VERSION["URB-Pyf"], SOURCECRS[ GEOGCRS["RGPF", DATUM["Reseau Geodesique de la Polynesie Francaise", 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",4687]]], TARGETCRS[ GEOGCRS["WGS 84", ENSEMBLE["World Geodetic System 1984 ensemble", MEMBER["World Geodetic System 1984 (Transit)"], MEMBER["World Geodetic System 1984 (G730)"], MEMBER["World Geodetic System 1984 (G873)"], MEMBER["World Geodetic System 1984 (G1150)"], MEMBER["World Geodetic System 1984 (G1674)"], MEMBER["World Geodetic System 1984 (G1762)"], MEMBER["World Geodetic System 1984 (G2139)"], ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]], ENSEMBLEACCURACY[2.0]], 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",4326]]], METHOD["Position Vector transformation (geog2D domain)", ID["EPSG",9606]], PARAMETER["X-axis translation",0.072, LENGTHUNIT["metre",1], ID["EPSG",8605]], PARAMETER["Y-axis translation",-0.507, LENGTHUNIT["metre",1], ID["EPSG",8606]], PARAMETER["Z-axis translation",-0.245, LENGTHUNIT["metre",1], ID["EPSG",8607]], PARAMETER["X-axis rotation",0.0183, ANGLEUNIT["arc-second",4.84813681109536E-06], ID["EPSG",8608]], PARAMETER["Y-axis rotation",-0.0003, ANGLEUNIT["arc-second",4.84813681109536E-06], ID["EPSG",8609]], PARAMETER["Z-axis rotation",0.007, ANGLEUNIT["arc-second",4.84813681109536E-06], ID["EPSG",8610]], PARAMETER["Scale difference",-0.0093, SCALEUNIT["parts per million",1E-06], ID["EPSG",8611]], OPERATIONACCURACY[0.5], USAGE[ SCOPE["Geodesy."], AREA["French Polynesia - onshore and offshore. Includes Society archipelago, Tuamotu archipelago, Marquesas Islands, Gambier Islands and Austral Islands."], BBOX[-31.24,-158.13,-4.52,-131.97]], ID["EPSG",8828], REMARK["SHOM report gives scale difference as 0.999 999 9907 (wrt unity). Transformation is to original Transit definition of WGS 84. It is consistent with later WGS 84 realisations G730, G873 and G1150 to no better than 1m."]]
Definition: PROJ.4
+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=GRS80 +step +proj=helmert +x=0.072 +y=-0.507 +z=-0.245 +rx=0.0183 +ry=-0.0003 +rz=0.007 +s=-0.0093 +convention=position_vector +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1
Definition: JavaScript (Proj4js)
proj4.defs("EPSG:8828","+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=GRS80 +step +proj=helmert +x=0.072 +y=-0.507 +z=-0.245 +rx=0.0183 +ry=-0.0003 +rz=0.007 +s=-0.0093 +convention=position_vector +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1");
Definition: JSON
{ "$schema": "https://proj.org/schemas/v0.7/projjson.schema.json", "type": "Transformation", "name": "RGPF to WGS 84 (1)", "source_crs": { "type": "GeographicCRS", "name": "RGPF", "datum": { "type": "GeodeticReferenceFrame", "name": "Reseau Geodesique de la Polynesie Francaise", "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": 4687 } }, "target_crs": { "type": "GeographicCRS", "name": "WGS 84", "datum_ensemble": { "name": "World Geodetic System 1984 ensemble", "members": [ { "name": "World Geodetic System 1984 (Transit)", "id": { "authority": "EPSG", "code": 1166 } }, { "name": "World Geodetic System 1984 (G730)", "id": { "authority": "EPSG", "code": 1152 } }, { "name": "World Geodetic System 1984 (G873)", "id": { "authority": "EPSG", "code": 1153 } }, { "name": "World Geodetic System 1984 (G1150)", "id": { "authority": "EPSG", "code": 1154 } }, { "name": "World Geodetic System 1984 (G1674)", "id": { "authority": "EPSG", "code": 1155 } }, { "name": "World Geodetic System 1984 (G1762)", "id": { "authority": "EPSG", "code": 1156 } }, { "name": "World Geodetic System 1984 (G2139)", "id": { "authority": "EPSG", "code": 1309 } } ], "ellipsoid": { "name": "WGS 84", "semi_major_axis": 6378137, "inverse_flattening": 298.257223563 }, "accuracy": "2.0", "id": { "authority": "EPSG", "code": 6326 } }, "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": 4326 } }, "method": { "name": "Position Vector transformation (geog2D domain)", "id": { "authority": "EPSG", "code": 9606 } }, "parameters": [ { "name": "X-axis translation", "value": 0.072, "unit": "metre", "id": { "authority": "EPSG", "code": 8605 } }, { "name": "Y-axis translation", "value": -0.507, "unit": "metre", "id": { "authority": "EPSG", "code": 8606 } }, { "name": "Z-axis translation", "value": -0.245, "unit": "metre", "id": { "authority": "EPSG", "code": 8607 } }, { "name": "X-axis rotation", "value": 0.0183, "unit": { "type": "AngularUnit", "name": "arc-second", "conversion_factor": 4.84813681109536e-06 }, "id": { "authority": "EPSG", "code": 8608 } }, { "name": "Y-axis rotation", "value": -0.0003, "unit": { "type": "AngularUnit", "name": "arc-second", "conversion_factor": 4.84813681109536e-06 }, "id": { "authority": "EPSG", "code": 8609 } }, { "name": "Z-axis rotation", "value": 0.007, "unit": { "type": "AngularUnit", "name": "arc-second", "conversion_factor": 4.84813681109536e-06 }, "id": { "authority": "EPSG", "code": 8610 } }, { "name": "Scale difference", "value": -0.0093, "unit": { "type": "ScaleUnit", "name": "parts per million", "conversion_factor": 1e-06 }, "id": { "authority": "EPSG", "code": 8611 } } ], "accuracy": "0.5", "scope": "Geodesy.", "area": "French Polynesia - onshore and offshore. Includes Society archipelago, Tuamotu archipelago, Marquesas Islands, Gambier Islands and Austral Islands.", "bbox": { "south_latitude": -31.24, "west_longitude": -158.13, "north_latitude": -4.52, "east_longitude": -131.97 }, "id": { "authority": "EPSG", "code": 8828 }, "remarks": "SHOM report gives scale difference as 0.999 999 9907 (wrt unity). Transformation is to original Transit definition of WGS 84. It is consistent with later WGS 84 realisations G730, G873 and G1150 to no better than 1m." }
Definition: MapServer - MAPfile
PROJECTION "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=GRS80" "step" "proj=helmert" "x=0.072" "y=-0.507" "z=-0.245" "rx=0.0183" "ry=-0.0003" "rz=0.007" "s=-0.0093" "convention=position_vector" "step" "inv" "proj=cart" "ellps=WGS84" "step" "proj=pop" "v_3" "step" "proj=unitconvert" "xy_in=rad" "xy_out=deg" "step" "proj=axisswap" "order=2,1" END
Definition: Mapnik
<?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=GRS80 +step +proj=helmert +x=0.072 +y=-0.507 +z=-0.245 +rx=0.0183 +ry=-0.0003 +rz=0.007 +s=-0.0093 +convention=position_vector +step +inv +proj=cart +ellps=WGS84 +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=GRS80 +step +proj=helmert +x=0.072 +y=-0.507 +z=-0.245 +rx=0.0183 +ry=-0.0003 +rz=0.007 +s=-0.0093 +convention=position_vector +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1"> </Layer> </Map>
Definition: SQL (PostGIS)
INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) values ( 8828, 'EPSG', 8828, '+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=GRS80 +step +proj=helmert +x=0.072 +y=-0.507 +z=-0.245 +rx=0.0183 +ry=-0.0003 +rz=0.007 +s=-0.0093 +convention=position_vector +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1', 'None');