Discussion :

[lhapaipai]

Bonjour,

J'aimerais pour un projet d'une web app offline utiliser la BD Ortho disponible via les geoservices de l'ign : https://geoservices.ign.fr/bdortho
Si l'on télécharge l'orthophoto d'un département, on obtient environ 500 rasters d'à peu près 100Mo chacun. et un fichier mosaique.vrt qui regroupe tous ces fichiers. Les rasters sont projetés en Lambert 93.

J'aurais donc besoin de générer des tuiles XYZ jusqu'au zoom 19 à l'échelle d'un département dans la projection Web Mercator.

Pour le moment j'ai essayé naïvement de faire

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
gdal2tiles.py -s EPSG:2154 --xyz --zoom=7-19 -w leaflet mosaique.vrt output

mais bien sûr je fais exploser mon ordinateur sans même générer la moindre tuile

Auriez-vous des pistes sur la manière de m'y prendre sachant que je n'ai à ma disposition qu'un ordinateur de bureau (processeur i7, 32Go de ram).

Merci et bonne soirée

edit: ou bien savez-vous si l'ign propose cette orthophoto en EPSG:3857 en téléchargement direct dans les géoservices ? elle est bien déjà disponible derrière le protocole WMTS.

edit2: peut-être que je pourrais déjà séparer les 2 étapes reprojection / construction des tuiles.
en reprojetant un fichier virt composé de 3 images dans une image. puis en appelant gdal2tiles sur cette image résultante mais sans reprojection j'arrive à générer mes tuiles (alors que si j'appelais gdal2tiles sur le fichier virt directement mon ordinateur ramait et ne produisait rien)
existe-t-il une possibilité d'appeler gdalwarp sur mon fichier mosaique.vrt et de générer en sortie non pas une seule image mais 500 images reprojetées de 25000px de côté environ ?

[lhapaipai]

Du coup j'ai résolu mon problème,
si cela peut intéresser quelqu'un,


d'abord j'ai calculé les coordonnées des bbox associées à toutes les tuiles de zoom 13 (dans la projection web mercator) nécessaires à la couverture de mon département.
j'ai ensuite reprojeté ma grosse mosaique.vrt dans chacune des bbox.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
gdalwarp ./mosaique_POSIX.vrt -s_srs EPSG:2154 -t_srs EPSG:3857 -te <ma-bbox> output/<z>-<x>-<y>.jp2

On se retrouve avec toujours autant d'images (voire plus) mais cette fois-ci dans la bonne projection. il est alors possible d'appeler gdal2tiles.py sur chaque petite image, le processus n'est donc plus bloquant.

j'ai donc appelé gdal2tiles.py sur chacune de ces images en allant du zoom 13 à 19 (pour les zooms d'avant je referai une mosaique.vrt avec mes petites tuiles de zoom 13.
pourquoi zoom 13 : ça fait des images de grosso modo 100Mo donc pratique à manipuler.
pourquoi ce choix de bbox : car lors de l'appel de gdal2tiles les tuiles viendront parfaitement se découper dans la bbox -> ça évite d'avoir des zones blanches.