Faire QGIS interpréter les coordonnées comme long-lat au lieu de lat-long?

13

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.

entrez la description de l'image ici

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.

Paolo Santini
la source
Une image de test serait très utile.
user30184
@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download donne beaucoup de liens.
AndreJ

Réponses:

7

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:

gdal_translate -of VRT HDF4_SDS:UNKNOWN:"3b31.hdf":0 0.vrt
gdal_translate -of VRT -gcp 0 0 -180 -40 -gcp 160 0 -180 40 -gcp 0 720 180 -40 -gcp 160 720 180 40 0.vrt 1.vrt
gdalwarp -r bilinear -t_srs EPSG:4326 1.vrt 0.tif

avec ce résultat:

entrez la description de l'image ici

Ma rotation doit suivre la balise de métadonnées Origin = SOUTHWEST.

AndreJ
la source
C'est exactement ce que je cherchais, bien que je ne puisse pas utiliser correctement les commandes que vous avez écrites. Pouvez-vous expliquer leur utilisation plus en détail?
Paolo Santini
Placez-les dans un fichier texte nommé rotate.batdans 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 ;-)
AndreJ
3

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.

entrez la description de l'image ici

Le rapport Gdalinfo de cette image est

gdalinfo 3B42RT.2016031806.7.7day.tif
Driver: GTiff/GeoTIFF
Files: 3B42RT.2016031806.7.7day.tif
       3B42RT.2016031806.7.7day.tfw
Size is 1440, 480
Coordinate System is `'
Origin = (-180.000000000000000,60.000000000000000)
Pixel Size = (0.250000000000000,-0.250000000000000)
Metadata:
  TIFFTAG_DATETIME=2016:03:18 10:25:21
  TIFFTAG_DOCUMENTNAME=./image/3B42RT.2016031806.7.7day.tif
  TIFFTAG_IMAGEDESCRIPTION=IDL TIFF file
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=IDL 7.1.1, ITT Visual Information Solutions
  TIFFTAG_XRESOLUTION=100
  TIFFTAG_YRESOLUTION=100
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-180.0000000,  60.0000000)
Lower Left  (-180.0000000, -60.0000000)
Upper Right ( 180.0000000,  60.0000000)
Lower Right ( 180.0000000, -60.0000000)
Center      (   0.0000000,   0.0000000)
Band 1 Block=1440x2 Type=UInt16, ColorInterp=Gray
user30184
la source
Il semble que les fichiers tif 3b42 aient un géoréférencement correct, contrairement aux fichiers hdf 3b31.
AndreJ
Il semble plutôt que HDF est compliqué et n'est pas parfaitement pris en charge par GDAL trac.osgeo.org/gdal/wiki/HDF .
user30184
Exactement certains de ceux que vous trouvez dans Mirador sont géoréférencés, en choisissant les propriétés NetCDF dans la fenêtre de personnalisation, d'autres non.
Paolo Santini
@ user30184, c'est au fournisseur de données d'inclure les balises de métadonnées de géoréférencement appropriées. Vous en avez pour votre
argent