Lassen Sie QGIS Koordinaten als Long-Lat statt Lat-Long interpretieren?

13

Ich arbeite mit den TRMM 3B31-Dateien über die von der NASA bereitgestellten Niederschläge. Wenn ich sie in QGIS lade und die WGS84-Projektion zusammen mit anderen Formdateien auswähle, werden die Karten vertikal angezeigt, dh QGIS interpretiert als Breite, was stattdessen Länge ist, und umgekehrt. Ich habe versucht, die Option gdal_translate -a_ullr auszuwählen, um gute Koordinaten zu erhalten, aber dies löst das Problem nicht: Die Bilder werden nur komprimiert und gestreckt, um meine Anforderungen zu erfüllen, aber dies ist eindeutig nicht das, was ich gerne hätte.

Wenn Sie eine CSV-Datei laden, können Sie auswählen, welcher Breitengrad und welcher Längengrad verwendet werden soll. Gibt es eine Möglichkeit, dasselbe für Raster zu tun?

Da ich ein Raster von 193 habe, wäre es großartig, wenn Sie einen Befehl vorschlagen könnten, der iterativ ausgeführt werden kann.

Bildbeschreibung hier eingeben

Das Bild zeigt eine Karte von Indien und das erste Band der 3B31-Datei, dh Oberflächenniederschlag. Die beiden schwarzen Linien sind die 40 und -40 parallel. Ich hätte sie gerne horizontal und daher überschneiden sich die Regendaten mit Indien.

Paolo Santini
quelle
Ein Testbild wäre sehr nützlich.
user30184
@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download enthält viele Links.
AndreJ

Antworten:

7

Ich habe die 3b31-hdf-Datei von mirador ausprobiert und die gleiche Ausrichtung wie bei Ihnen festgestellt. Also habe ich ein GCP erstellt, um eine Drehung gegen den Uhrzeigersinn durchzuführen:

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

mit diesem Ergebnis:

Bildbeschreibung hier eingeben

Meine Drehung sollte dem Metadaten-Tag Origin = SOUTHWEST folgen.

AndreJ
quelle
Dies ist genau das, wonach ich gesucht habe, obwohl ich die von Ihnen geschriebenen Befehle nicht richtig verwenden kann. Können Sie ihre Verwendung näher erläutern?
Paolo Santini
Fügen Sie sie in eine Textdatei ein, die sich rotate.batim selben Ordner wie die HDF-Datei befindet. Überprüfen Sie den Namen der Quelldatei, starten Sie die OSGEO4W-Shell, navigieren Sie zum Ordner mit der CD ... und führen Sie den Stapel aus. Gute alte DOS-Navigation ist immer nützlich ;-)
AndreJ
3

Ich habe 3B42RT.2016031806.7.7day.tif und .tfw heruntergeladen, aber in der Georeferenzierung sehe ich nichts Seltsames.

Bildbeschreibung hier eingeben

Gdalinfo Bericht von diesem Bild ist

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
quelle
Es scheint, dass die 3b42-TIF-Dateien über die richtige Georeferenzierung verfügen, die 3b31-HDF-Dateien dagegen nicht.
AndreJ
Es scheint eher so, als ob HDF kompliziert ist und von GDAL trac.osgeo.org/gdal/wiki/HDF nicht perfekt unterstützt wird .
user30184
Genau einige von denen, die Sie in Mirador finden, sind georeferenziert. Wählen Sie die NetCDF-Eigenschaften im Anpassungsfenster aus, andere nicht.
Paolo Santini
@ user30184 Es ist Sache des Datenanbieters, geeignete Georeferenzierungs-Metadaten-Tags einzuschließen. Sie bekommen, wofür Sie bezahlt haben ...
AndreJ