Suite

Gdal_translate : erreur de projection ?

Gdal_translate : erreur de projection ?


J'utilise les routines de la bibliothèque GDAL pour créer une hiérarchie d'images en mosaïque à utiliser avec Google Earth. Les images ne sont que des fichiers PNG créés à l'aide de Matplotlib. Pour mosaïquer les images, nous utilisonsgdal_translate,gdalwarpetgdal2tiles.py. Lors de l'utilisation de l'installation Macports de GDAL (installée viainstaller le port sudo gdal, tout fonctionne parfaitement et la hiérarchie en mosaïque est créée et chargée dans Google Earth. Mais pour des raisons que je ne comprends pas, l'installation d'Anaconda nous renvoie une erreur.

Voici les résultats degdalinfosur un fichier PNG type que nous produisons :

% gdalinfo plot.png"> cross-posting cette question, mais je viens de découvrir ce site d'échange de piles et je ne sais pas encore comment migrer ma question).


Il y a environ une semaine, l'emballage Anaconda de GDAL inclut désormais tous les fichiers de projection. La commande d'installation pour l'empaquetage Anaconda de GDAL est

% conda installer gdal

Les fichiers de projection sont stockés sous$ANACONDA/gdal/share$ donc réglageGDAL_DATA' vers ce chemin est tout ce qui est nécessaire pour éviter les problèmes que j'avais ci-dessus.

% export GDAL_DATA=$ANACONDA/share/gdal

Plus besoin de compter sur MacPorts, Frameworks, HomeBrew, gdal.org et ainsi de suite pour ces fichiers supplémentaires.


(Note : ceci est une réponse partielle !)

J'ai compris que les fichiers .vrt créés par l'installation de MacPorts et l'installation d'Anaconda étaient différents. La chaîne "WKT" pour la version Anaconda contenait une chaîne supplémentaire "TOWGS84[0,0,0,0,0,0,0]", c'est-à-dire

Plot.vrt (Anaconda)

GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], TOWGS84[0,0,0,0,0,0, 0], AUTORITE["EPSG","6326"]], PRIMEM["Greenwich",0, AUTORITE["EPSG","8901"]], UNITÉ["degré",0.0174532925199433, AUTORITE["EPSG"," 9108"]], AUTORITÉ["EPSG","4326"]]

MacPorts (plot.vrt)

GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTORITE["EPSG","7030"]], AUTORITE["EPSG","6326"]], PRIMEM[ "Greenwich",0, AUTORITE["EPSG","8901"]], UNITÉ["degré",0.0174532925199433, AUTORITE["EPSG","9122"]], AUTORITE["EPSG","4326"]]

J'ai donc créé un fichier "pro.prj" qui contient la chaîne WKT (la version Anaconda, avec la chaîne supplémentaire TOWGS84[0,0,0,0,0,0,0,0]"), et l'ai transmis à 'gdal_translate'

% gdal_translate -of VRT -a_srs proj.prj -gcp 0 0 -120 0 -gcp 2400 0 -60 0 -gcp 2400 2400 -60 -60 -90 plot.png">

Problème résolu! Avec le recul, le problème était très simple.

L'installation de MacPorts fonctionnerait correctement (bien que gdal2tiles.py ait émis des avertissements) sans la variable d'environnement GDAL_DATA définie. J'en étais donc venu à croire que mon problème était suffisamment "simple" pour ne pas nécessiter de fichiers de données supplémentaires.

Cependant, la version Anaconda ne fonctionnerait pas sans GDAL_DATA correctement défini. Et pour une raison quelconque, les fichiers de données installés avec Macports ne fonctionnent pas pour Anaconda (les noms de fichiers sont tous les mêmes, mais il existe de subtiles différences de contenu).

Anaconda fonctionnera avec ses propres fichiers de données :

export GDAL_DATA=/Bibliothèque/Frameworks/GDAL.framework/Versions/Current/Resources/gdal

Notez que cela fonctionne également pour MacPorts.

Et pour mémoire, même si MacPorts ne semblait pas avoir besoin d'un paramètre GDAL_DATA pour mon problème simple, le paramètre

exporter GDAL_DATA=/opt/local/share/gdal

supprime tous les avertissements ennuyeux que gdal2tiles.py émet avec Macports. Encore une fois, ces fichiers de données ne fonctionnent PAS avec Anaconda, pour des raisons que je ne pouvais même pas commencer à comprendre !


Voir la vidéo: 041614 Change map projection