QGIS'i koordinatları lat-long yerine long-lat olarak yorumlamak mı?


13

NASA tarafından sağlanan yağış üzerinde TRMM 3B31 dosyalarıyla çalışıyorum. Onları diğer şekil dosyalarıyla birlikte WGS84 projeksiyonunu seçerek QGIS'e yüklediğimde, haritalar dikey olarak görüntülenir, yani QGIS bunun yerine enlem ve boylam olarak yorumlar. İyi koordinatlar heybetli seçerek gdal_translate seçeneği -a_ullr denedim ama sorunu çözmüyor: görüntüler sadece benim gereksinimleri karşılamak için sıkıştırılmış ve gergin ama bu açıkça ne olmasını istiyorum değil.

Bir .csv dosyası yüklediğinizde, enlem ve boylamın ne olması gerektiğini seçebilirsiniz, rasterler için de aynısını yapmanın bir yolu var mı?

193 tarama cihazım olduğu için tekrarlı olarak gerçekleştirilebilecek bir komut önerebilirseniz harika olur.

resim açıklamasını buraya girin

Resim Hindistan'ın bir haritasını ve 3B31 dosyasının ilk bandını, yani yüzey çökelmesini göstermektedir. İki siyah çizgi 40 ve -40 paraleldir. Bunların yatay olmasını ve dolayısıyla yağmur verilerinin Hindistan'la çakışmasını istiyorum.


Bir test görüntüsü çok faydalı olacaktır.
user30184

@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download birçok bağlantı sağlar.
AndreJ

Yanıtlar:


7

Ben den 3b31 HDF dosyasını çalıştı mirador ve karşılaştığınız aynı yönelimi bulundu. Bu yüzden saat yönünün tersine döndürmek için bazı GCP'ler oluşturdum:

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

bu sonuçla:

resim açıklamasını buraya girin

Rotasyonum Origin = SOUTHWEST meta veri etiketine uygun olmalıdır.


Tam da aradığım şey buydu, ancak yazdığınız komutları doğru şekilde kullanamıyorum. Kullanımlarını daha ayrıntılı olarak açıklayabilir misiniz?
Paolo Santini

Bunları rotate.batHDF dosyasıyla aynı klasörde adlandırılmış bir metin dosyasına koyun . Kaynak dosya adını kontrol edin, OSGEO4W Kabuğunu başlatın, cd ... ile klasöre gidin ve toplu işi çalıştırın. İyi eski DOS gezinme her zaman yararlıdır ;-)
AndreJ

3

3B42RT.2016031806.7.7day.tif ve .tfw dosyasını indirdim ancak coğrafi referanslarda tuhaf bir şey görmüyorum.

resim açıklamasını buraya girin

Bu resmin Gdalinfo raporu

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

Görünüşe göre 3b42 tif dosyaları doğru coğrafi referansa sahipken, 3b31 hdf dosyaları sahip değil.
AndreJ

HDF'nin karmaşık olduğu ve GDAL trac.osgeo.org/gdal/wiki/HDF tarafından mükemmel bir şekilde desteklenmediği görülmektedir .
user30184

Mirador'da bulduğunuzlardan bazıları, özelleştirme penceresindeki NetCDF özelliklerini seçerken, coğrafi olarak referanslıdır, diğerleri değildir.
Paolo Santini

@ user30184 uygun coğrafi referans meta veri etiketlerini eklemek veri sağlayıcısına bağlıdır. Sizin için ne ödeme olsun ...
AndreJ
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.