[Francophone] Statut registre IGNF dans PostGIS

[Didier Richard]

Bonjour,

(désolé pour mon temps de réponse)

Quelques petites informations si ça peut aider étant à l'origine des registres IGN (xml) et des fichiers comme IGNF, la grille ntf_r93.gsb, fichiers conçus avec l'aide du Service de la Géodésie de l'IGN et de la cellule normalisation de l'époque.

Even m'avait interpellé justement sur la mise à jour du registre il y a quelques mois ... j'ai pris quelques informations depuis, mais pas d'action réelle.

Le registre des systèmes de référence de coordonnées, ainsi que des méthodes de passage d'un système à un autre est diffusé à l'adresse : http://librairies.ign.fr/geoportail/resources/IGNF.xml 
C'est un fichier qui s'appuie sur la norme ISO 19139 ; il décrit les systèmes, les paramètres et les méthodes de transformation entre système et qui sont de la responsabilité de l'IGN. Il fait aussi le lien entre nomenclature IGN et quant elle existe la nomenclature EPSG (enfin, OGP maintenant). Il décrit aussi les transformations (interpolation géocentrique, conversion altimétrique, lambert conique conforme sécant à deux parallèles, etc ...). Du lourd !D
Il existe aussi un registre en ligne : https://registre.ign.fr/ign/IGNF/IGNF/ qui permet aux logiciels qui consomment le registre statique IGNF.xml de naviguer dans les systèmes, paramètres, etc ... registre mis en ligne par l'équipe Géoportail.

Ce registre XML est régulièrement mis à jour (dernière version : mars 2018).

Sur la même page du site de la géodésie (Cf. https://geodesie.ign.fr/index.php?page=documentation#titre3), on trouve aussi un lien vers un script SQL d'import dans PostGIS, fichier que j'avais conçu à l'époque (2008 !-), donc plus à jour ...
Ce qui est important dans ce fichier SQL, c'est en fait la dernière colonne de la table spatial_ref_sys qui contient la chaîne PROJ4 qui fait donc référence au registre IGNF que l'on trouve dans le projet PROJ (sans le 4 maintenant).
C'était la seule façon de faire en sorte que sous PostGIS on puisse effectuer des transformations correctes entre NTF et RGF93 (passage de Lambert II étendue à Lambert 93 par exemple) avec les grilles. La transformation via les codes EPSG étant incorrecte car ne prenant pas en compte la grille d'interpolation bilinéaire nécessaire pour corriger le changement de géoïde.

Idéalement, il faudrait usiner un script pour fabriquer automagiquement les registres IGNF (PROJ, WKT et PostGIS). J'ai un tel script XSLT qui s'occupe uniquement de PROJ et qui date de ... 2008 (mais la structure du IGNF.xml n'a pas changé, donc ça devrait le faire encore) !

Le véritable enjeu de la transformation du registre IGNF (xml ISO 19139) en quelque chose de consommable par PROJ et ces consommateurs, c'est celui des grilles car il en existe un paquet ... Elles sont toutes accessibles via https://geodesie.ign.fr/index.php?page=documentation#titre6 qui donne aussi bien la grille NTF RGF93 (Cf. https://geodesie.ign.fr/contenu/fichiers/gr3df97a.txt ) que les grilles de conversion altimétriques :
* http://geodesie.ign.fr/contenu/fichiers/documentation/grilles/outremer/ggm00.txt
* http://geodesie.ign.fr/contenu/fichiers/documentation/grilles/outremer/ggg00_ls.txt
* ggg00_lsv2.mnt
* http://geodesie.ign.fr/contenu/fichiers/documentation/grilles/outremer/ggg00_sb.txt 
* ...
Elles sont toutes là : https://geodesie.ign.fr/index.php?page=grilles

Le hic, c'est :

1) IGNF.xml en référence certaines complètement (URL), certaines par le nom (e.g. ggg00_lsv2.mnt)
2) IGNF.xml en référence certaines par un format (.txt), certaines par un autre format (.mnt)
    (Cf. https://geodesie.ign.fr/contenu/fichiers/documentation/grilles/notices/Grilles-MNT-TXT_Formats.pdf)
3) IGNF donne des noms qui ne sont pas très "parlants" !

Il faut donc transformer ces grilles dans un format consommable par PROJ (c'est ce que l'IGN a fait à l'époque pour la grille NTF <-> RGF93) ... et fabriquer ainsi les bonnes chaînes WKT et PROJ (et mettre à jour le registre IGNF dans PROJ avec les quelques nouveaux systèmes sur les îles entre autres).

Voilà, si on peut effectivement avancer collectivement sur la prise en compte de tout cela, ça serait hyper-cool !

@+
--
RICHARD Didier - Chef du service ValiLab
http://fr.linkedin.com/pub/didier-richard/98/2a3/a8/ - https://www.osgeo.org/member/didier/
6/8 avenue Blaise Pascal - BP Champs-sur-Marne - 77455 MARNE-LA-VALLÉE CEDEX 2
Tél : +33 (0) 1 43 98 83 23


________________________________________
De : Francophone [francophone-bounces at lists.osgeo.org] de la part de Even Rouault [even.rouault at spatialys.com]
Envoyé : mardi 25 septembre 2018 10:33
À : Régis Haubourg
Cc : Francophone local chapter discussions
Objet : Re: [Francophone] Statut registre IGNF dans PostGIS

On mardi 25 septembre 2018 10:07:04 CEST Régis Haubourg wrote:
> Merci Even pour les précisions, la cible pour PostGIS3 me parait assez
> claire.
>
> Plus je creuse en amont et moins je comprends à quel point l'utilisation de
> grille est décrit par les standards EPSG/OGC.

Quelques liens "utiles":

- OGC Topic 2 - Spatial referencing by coordinates / ISO 19111:2017
http://portal.opengeospatial.org/files/?artifact_id=39049
C'est la modélisation théorique des CRS et des notions de transformation de coordonnées

- WKT 2: http://docs.opengeospatial.org/is/12-063r5/12-063r5.html
Une implémentation basée sur cette modélisation abstraite

- EPSG Guidance Note 7-2
https://www.iogp.org/bookstore/product/coordinate-conversions-and-transformation-including-formulas/
qui décrit les différents types de transformation possibles

>
> Un exemple pour OSGB 1936
>  -  http://spatialreference.org/ref/epsg/4277/ ne mentionne aucune grille.

Comme son nom l'indique :-), spatialreference.org n'est ni un site de référence
ni maintenu depuis une bonne dizaine d'années.
Ceci dit il est normal qu'une définition de CRS ne contienne pas de grille.
Un CRS existe en tant que lui-même. Il n'est pas défini par rapport à un autre.
La notion de noeu TOWGS84 ou EXTENSION["PROJ4_GRIDS", "..."] dans le WKT1 ou
dans les chaines PROJ.4 est un "hack" pour faciliter la vie en passant
par un système WGS84 pivot, ce qui est loin d'être idéal car induit des pertes
de précision.
Cf diapos 26 et 27 de
https://download.osgeo.org/gdal/presentations/GDAL%202.3%20(FOSS4G-E%202018)%20-%2020%20year%20already%20and%20heading%20to%20the%20cloud.pdf
pour quelques limitations de cette approche dite "early binding" vs
l'approche "late binding" où on a d'une part les définitions de CRS et d'autre
part un catalogue de transformations dispoibles pour des couples (src_srs, dst_srs).

>  -  EPSG.io, maintenu par Klokantech,  mentionne  qu'il y a une grille de
> type NTV2 qui permet des reprojection à 1m.Le nom de la grille est fourni
> dans la page, mais n'est disponible dans aucune des chaînes de définition
> EPSG / OGC / Proj, ESRI, etc..

Tu parles de OSTN02_NTv2.gsb ?
C'est discuté dans https://github.com/OSGeo/proj-datumgrid/issues/21
la grille est dispo
https://www.ordnancesurvey.co.uk/business-and-government/help-and-support/navigation-technology/os-net/ostn02-ntv2-format.html

PROJ a maintenu un dépot dédié avec plusieurs paquets de grilles: un paquet
mondial de taille modeste avec les grilles les plus courantes et basse résolution,
et des paquets continentaux/régionaux avec des grilles plus précises:
https://github.com/OSGeo/proj-datumgrid/

Mais oui de manière générale, il peut y avoir des écarts dans le nommage des grilles ou le format effectivement
utilisé. D'où mon idée dans mes travaux à venir d'avoir une sorte de table
d'alias pour faire l'adaptation entre le nom/format du registre officiel et
celui effectivement utilisé.

>
> Est-ce que cela signifie qu'aucun organisme ne standardise cette notion de
> grille?
> Est-ce que cela implique que c'est à l'utilisateur, comme dans Geoserver
> [0], d'ajouter lui même les transformations qu'il souhaite?
> Est-ce que cela nous autorise arbitrairement à modifier les définitions
> EPSG Lambert pour ajouter la grille ign ?

C'est ce qui est fait dans le fichier IGNF de PROJ avec les chains +proj= qui
contiennent un +nadgrids.

> Cela serait bien le plus simple
> pour attendre un postGIS3 plus propre au niveau de la gestion des registres
> et des transformations.
>
> Régis (en croisade pour une utilisation des grilles simple et à la portée
> de tous)

Bonne chance !

--
Spatialys - Geospatial professional services
http://www.spatialys.com
_______________________________________________
Francophone mailing list
Francophone at lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/francophone