Far interpretare le coordinate QGIS come long-lat anziché lat-long?


13

Sto lavorando con i file TRMM 3B31 sulle precipitazioni fornite dalla NASA. Quando li carico su QGIS selezionando la proiezione WGS84 insieme ad altri file di forma, le mappe vengono visualizzate verticalmente, ovvero QGIS interpreta come latitudine ciò che invece è longitudine e viceversa. Ho provato l'opzione gdal_translate -a_ullr selezionando imponendo buone coordinate ma non risolve il problema: le immagini sono solo compresse e allungate per soddisfare le mie esigenze, ma chiaramente non è quello che vorrei avere.

Quando carichi un file .csv puoi scegliere quale deve essere la latitudine e quale la longitudine, c'è un modo per fare lo stesso con i raster?

Dal momento che ho 193 raster sarebbe fantastico se puoi suggerire un comando che può essere eseguito in modo iterativo.

inserisci qui la descrizione dell'immagine

L'immagine mostra una mappa dell'India e la prima banda del file 3B31, ovvero le precipitazioni di superficie. Le due linee nere sono il parallelo 40 e -40. Mi piacerebbe averli orizzontali e quindi i dati sulla pioggia si sovrappongono all'India.


Un'immagine di prova sarebbe molto utile.
user30184

@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download fornisce molti collegamenti.
AndreJ,

Risposte:


7

Ho cercato il file 3b31 HDF da mirador ed ho trovato lo stesso orientamento che hai incontrato. Quindi ho creato un GCP per eseguire una rotazione in senso antiorario:

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

con questo risultato:

inserisci qui la descrizione dell'immagine

La mia rotazione dovrebbe seguire il tag metadati Origin = SOUTHWEST.


Questo è esattamente quello che stavo cercando, anche se non sono in grado di utilizzare i comandi che hai scritto nel modo corretto. Puoi spiegarne l'uso in modo più dettagliato?
Paolo Santini,

Inseriscili in un file di testo denominato rotate.batnella stessa cartella del file HDF. Controllare il nome del file di origine, avviare OSGEO4W Shell, accedere alla cartella con cd ... ed eseguire il batch. La buona vecchia navigazione DOS è sempre utile ;-)
AndreJ

3

Ho scaricato 3B42RT.2016031806.7.7day.tif e .tfw ma non vedo nulla di strano nella georeferenziazione.

inserisci qui la descrizione dell'immagine

Il rapporto Gdalinfo di questa immagine è

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

Sembra che i file tif 3b42 abbiano una corretta georeferenziazione, mentre i file hdf 3b31 no.
AndreJ,

Sembra piuttosto che HDF sia complicato e non perfettamente supportato da GDAL trac.osgeo.org/gdal/wiki/HDF .
user30184

Alcuni di quelli che trovi in ​​Mirador sono georeferenziati, scegliendo le proprietà NetCDF nella finestra di personalizzazione, altri no.
Paolo Santini,

@ user30184 spetta al fornitore di dati includere tag di metadati di georeferenziazione adeguati. Ottieni ciò per cui hai pagato ...
AndreJ,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.