#1 Mon 27 December 2021 23:32
- clement.robert
- Juste Inscrit !
- Date d'inscription: 27 Dec 2021
- Messages: 3
Convertir les dataset RGE Alti 5m
Bonjour,
Je développe une application web pour laquel j'ai besoin d'avoir une altitude précise pour une latitude et une longitude donnée.
L'IGN propose un webservice qui fait cela ici https://geoservices.ign.fr/documentatio … rique-rest
Mais ce que je souhaiterais c'est de disposer de genre d'API sur mes serveurs pour ne pas être limité en nombre de requêtes / dépendant du geoservice IGN.
Pour cela j'ai récupéré les RGE Alti 5m en format ASC.
Il y a un projet open-source (opentodata) qui permet de faire ce que je souhaite en "self-hosted".
https://www.opentopodata.org/server/
Là ou ça coince, c'est au niveau de l'ajout du dataset de l'IGN.
Le format ne semble pas correspondre, et je ne sais pas comment convertir les données .asc en un format compatible avec ceux supportés par openToppoData
(.tiff, .hgt, .jp2)
Plus précisément:
"""
Open Topo Data supports all georeferenced raster formats supported by GDAL (e.g, .tiff, .hgt, .jp2).
Datasets can take one of two formats:
A single raster file.
A collection of square raster tiles which follow the SRTM naming convention: the file is named for the lower left corner. So a file named N30W120.tiff would span from 30 to 31 degrees latitude, and -120 to -119 degrees longitude. By default tiles are 1° by 1° and the coordinates are in WGS84, but this can be configured.
If your dataset consists of multiple files that aren't on a nice grid, you can create a .vrt file pointing to the files that Open Topo Data will treat as a single-file dataset. For an example of this process, see the documentation for configuring EMODnet.
"""
Est-ce que l'un d'entre-vous sait si il est possible de convertir cela en un fichier avec un logiciel ou au format / convention de nommage type SRTM ?
Ou à défaut, est-ce qu'il existe une solution similaire au projet opentopodata qui me permettrait d'utiliser le dataset IGN pour me renvoyer via une API HTTP une altitude précise ?
Par avance je vous remercie!
Hors ligne
#2 Tue 28 December 2021 11:11
- ppluvinet
- Participant assidu
- Lieu: VALENCE
- Date d'inscription: 6 Aug 2007
- Messages: 617
Re: Convertir les dataset RGE Alti 5m
Bonjour,
Il doit être possible de convertir les asc en tiff via la librairie GDAL https://gdal.org/programs/gdal_translate.html
GDAL peut être utilisable soit en ligne de commande soit via le logiciel QGIS. Sous QGIS, vous pouvez importer le fichier ASCII et le réexporter dans le format souhaité.
Bonne continuation,
Pascal PLUVINET
Hors ligne
#3 Tue 28 December 2021 15:33
- cquest
- Participant assidu
- Date d'inscription: 6 Jan 2013
- Messages: 874
Re: Convertir les dataset RGE Alti 5m
Une version du RGEALTI 1m repackagée en TIFF COG est disponible sur https://data.cquest.org/ign/rgealti/repack/cog/
Christian Quest - https://amicale.net/@cquest sur Mastodon (terminé twitter/X)
Membre fondateur et porte parole d'OpenStreetMap France
Initiateur de opendatArchives, OpenEventDatabase, Panoramax
Hors ligne
#4 Wed 29 December 2021 09:06
- clement.robert
- Juste Inscrit !
- Date d'inscription: 27 Dec 2021
- Messages: 3
Re: Convertir les dataset RGE Alti 5m
Merci à vous deux pour ces précieuses informations .
De ce que je comprends, comme opentodata a besoin de fichiers tiff qui suivent une convention de nommage bien spécifique (celle de SRTM), il faut que je lui fournisse le fichier .vrt généré avec GDAL pour qu'il fasse correctement le lien entres les tiffs.
Je vais voir si je peux partir de celui que tu as généré: https://data.cquest.org/ign/rgealti/repack/
Je suis un peu occupé cette semaine, mais j'essaierais de faire un retour pour les intéressés quand j'aurais eu le temps de tester.
Hors ligne
#5 Thu 30 December 2021 10:23
- clement.robert
- Juste Inscrit !
- Date d'inscription: 27 Dec 2021
- Messages: 3
Re: Convertir les dataset RGE Alti 5m
J'ai finalement eu le temps de tester et ça fonctionne.
Juste avec la création d'un VRT:
gdalbuildvrt -overwrite -a_srs "EPSG:2154" my.vrt ../ignasc/*.asc
Puis pour optimiser la taille:
gdal_translate -of GTiff -co "TILED=YES" -co COMPRESS=LZW -co BIGTIFF=YES -ot Float32 -a_srs "EPSG:2154" my.vrt my.tiff
Hors ligne