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.
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.
sumber
Jawaban:
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:
dengan hasil ini:
Rotasi saya harus mengikuti tag metadata Origin = SOUTHWEST.
sumber
rotate.bat
dalam 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 ;-)Saya mengunduh 3B42RT.2016031806.7.7day.tif dan .tfw tapi saya tidak melihat ada yang aneh di georeferensi.
Laporan Gdalinfo dari gambar ini adalah
sumber