Je travaille avec les fichiers TRMM 3B31 sur les précipitations fournis par la NASA. Lorsque je les charge sur QGIS en sélectionnant la projection WGS84 avec d'autres fichiers de forme, les cartes sont affichées verticalement, c'est-à-dire que QGIS interprète comme latitude ce qui est à la place longitude et vice-versa. J'ai essayé l'option gdal_translate -a_ullr sélectionnant l'imposition de bonnes coordonnées mais cela ne résout pas le problème: les images sont uniquement compressées et étirées pour satisfaire mes exigences mais ce n'est clairement pas ce que j'aimerais avoir.
Lorsque vous chargez un fichier .csv, vous pouvez choisir ce qui doit être la latitude et la longitude, existe-t-il un moyen de faire de même pour les rasters?
Puisque j'ai 193 raster, ce serait bien si vous pouviez suggérer une commande qui peut être exécutée de manière itérative.
L'image montre une carte de l'Inde et la première bande du fichier 3B31, c'est-à-dire les précipitations en surface. Les deux lignes noires sont les parallèles 40 et -40. Je voudrais les avoir horizontaux et donc que les données de pluie chevauchent l'Inde.
la source
Réponses:
J'ai essayé le fichier hdf 3b31 de mirador et j'ai trouvé la même orientation que vous avez rencontrée. J'ai donc créé du GCP pour effectuer une rotation dans le sens inverse des aiguilles d'une montre:
avec ce résultat:
Ma rotation doit suivre la balise de métadonnées Origin = SOUTHWEST.
la source
rotate.bat
dans le même dossier que le fichier HDF. Vérifiez le nom du fichier source, démarrez le shell OSGEO4W, accédez au dossier contenant cd ... et exécutez le lot. Une bonne vieille navigation DOS est toujours utile ;-)J'ai téléchargé 3B42RT.2016031806.7.7day.tif et .tfw mais je ne vois rien d'étrange dans le géoréférencement.
Le rapport Gdalinfo de cette image est
la source