[Francophone] Statut registre IGNF dans PostGIS

[Even Rouault]

Salut Didier,

Merci pour ces infos très précieuses.

Pour l'instant, dans le cadre de mes travaux pour PROJ 6, j'avais fait il y a 
quelques jours une moulinette d'ingestion du fichier "IGNF" qui contient les 
chaines PROJ.4 pour les mettre dans un formalisme compatible de la 
structuration ISO-19111 (telle qu'employée par la base de données EPSG), qui 
sera utilisé pour la nouvelle base de données de CRS et de transformations de 
coordonnées

Le résultat est là
https://github.com/rouault/proj.4/blob/iso19111/data/sql/ignf.sql

Mais il est clair qu'il aurait été mieux de repartir
de https://registre.ign.fr/ign/IGNF/IGNF/ si j'avais su. Pour une prochaine 
itération...

Plus de détails et contexte général dans
https://github.com/OSGeo/proj.4/pull/1149
http://even.rouault.free.fr/proj_cpp_api/rfc-2.html#database

Exemple:

- Requête de la définition d'IGNF:LAMB1 en chaine WKT2

$ projinfo IGNF:LAMB1 -o WKT2_2018

WKT2_2018 string: 
PROJCRS["Lambert I",
    BASEGEOGCRS["Nouvelle Triangulation Francaise Paris grades",
        DATUM["Nouvelle Triangulation Francaise Paris grades",
            ELLIPSOID["Clarke 1880 (IGN)",6378249.2,293.466021293627,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Paris",2.5969213,
            ANGLEUNIT["grad",0.0157079632679489]]],
    CONVERSION["Conversion for Lambert I",
        METHOD["Lambert Conic Conformal (1SP)",
            ID["EPSG",9801]],
        PARAMETER["Latitude of natural origin",49.5,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",0.99987734,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",600000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",200000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["easting (X)",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["northing (Y)",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["unknown"],
        AREA["World"],
        BBOX[-90,-180,90,180]],
    ID["IGNF","LAMB1"]]

- Requête des transformations possibles de IGNF:LAMB1 vers EPSG:4326 sous 
forme de pipeline de transformation compatible PROJ5

$ projinfo -s IGNF:LAMB1 -t EPSG:4326 -o PROJ
-------------------------------------
Operation n°1:

PROJ string: 
+proj=pipeline +step +inv +proj=lcc +lat_1=49.5 +lat_0=49.5 +lon_0=0 
+k_0=0.99987734 +x_0=600000 +y_0=200000 +ellps=clrk80ign +pm=paris +step 
+proj=hgridshift +grids=ntf_r93.gsb +step +proj=unitconvert +xy_in=rad 
+xy_out=deg +step +proj=axisswap +order=2,1

-------------------------------------
Operation n°2:

PROJ string: 
+proj=pipeline +step +inv +proj=lcc +lat_1=49.5 +lat_0=49.5 +lon_0=0 
+k_0=0.99987734 +x_0=600000 +y_0=200000 +ellps=clrk80ign +pm=paris +step 
+proj=cart +ellps=clrk80ign +step +proj=helmert +x=-168 +y=-60 +z=320 +step 
+inv +proj=cart +ellps=WGS84 +step +proj=unitconvert +xy_in=rad +xy_out=deg 
+step +proj=axisswap +order=2,1

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com