Membuat QGIS menginterpretasikan koordinat sebagai panjang-panjang dan bukan panjang-panjang?

13

Saya bekerja dengan file TRMM 3B31 tentang presipitasi yang disediakan oleh NASA. Ketika saya memuatnya pada QGIS memilih proyeksi WGS84 bersama-sama dengan file bentuk lain peta ditampilkan secara vertikal, yaitu QGIS mengartikan sebagai garis lintang apa yang bukan garis bujur dan sebaliknya. Saya telah mencoba opsi gdal_translate -a_ullr memilih penerapan koordinat yang baik tetapi tidak menyelesaikan masalah: gambar hanya dikompresi dan diregangkan untuk memenuhi persyaratan saya, tetapi ini jelas bukan yang ingin saya miliki.

Ketika Anda memuat file .csv Anda dapat memilih apa yang harus menjadi garis lintang dan garis bujur apa, apakah ada cara untuk melakukan hal yang sama untuk raster?

Karena saya memiliki 193 raster, akan lebih bagus jika Anda dapat menyarankan perintah yang dapat dilakukan iteratif.

masukkan deskripsi gambar di sini

Gambar menunjukkan peta India dan pita pertama dari file 3B31, yaitu presipitasi permukaan. Dua garis hitam adalah paralel 40 dan -40. Saya ingin mereka horisontal dan karenanya data hujan tumpang tindih dengan India.

Paolo Santini
sumber
Gambar uji akan sangat berguna.
user30184
@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download memberi banyak tautan.
AndreJ

Jawaban:

7

Saya mencoba file hdb 3b31 dari mirador dan menemukan orientasi yang sama seperti yang Anda temui. Jadi saya membuat beberapa GCP untuk melakukan rotasi berlawanan arah jarum jam:

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

dengan hasil ini:

masukkan deskripsi gambar di sini

Rotasi saya harus mengikuti tag metadata Origin = SOUTHWEST.

AndreJ
sumber
Ini persis apa yang saya cari, meskipun saya tidak dapat menggunakan perintah yang Anda tulis dengan cara yang tepat. Bisakah Anda menjelaskan penggunaannya secara lebih rinci?
Paolo Santini
Masukkan mereka ke dalam file teks bernama rotate.batdalam folder yang sama dengan file HDF. Periksa nama file sumber, mulai OSGEO4W Shell, navigasikan ke folder dengan cd ... dan jalankan batch. Navigasi DOS lama yang baik selalu berguna ;-)
AndreJ
3

Saya mengunduh 3B42RT.2016031806.7.7day.tif dan .tfw tapi saya tidak melihat ada yang aneh di georeferensi.

masukkan deskripsi gambar di sini

Laporan Gdalinfo dari gambar ini adalah

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
pengguna30184
sumber
Tampaknya file tif 3b42 memiliki georeferensi yang benar, sedangkan file hdf 3b31 belum.
AndreJ
Tampaknya HDF rumit dan tidak sepenuhnya didukung oleh GDAL trac.osgeo.org/gdal/wiki/HDF .
user30184
Tepatnya, beberapa yang Anda temukan di Mirador mengalami georeferensiasi, memilih properti NetCDF di jendela penyesuaian, yang lain tidak.
Paolo Santini
@ user30184 terserah penyedia data untuk menyertakan tag metadata georeferensi yang tepat. Anda mendapatkan apa yang Anda bayar ...
AndreJ