Czy QGIS interpretuje współrzędne jako długie-długie zamiast lat-długie?

13

Pracuję z plikami TRMM 3B31 dotyczącymi opadów atmosferycznych dostarczonych przez NASA. Kiedy ładuję je do QGIS, wybierając projekcję WGS84 wraz z innymi plikami kształtów, mapy są wyświetlane pionowo, tzn. QGIS interpretuje jako szerokość geograficzną, co zamiast tego jest długością geograficzną i odwrotnie. Próbowałem opcji gdal_translate -a_ullr wybierając narzucanie dobrych współrzędnych, ale to nie rozwiązuje problemu: obrazy są tylko kompresowane i rozciągane, aby spełnić moje wymagania, ale wyraźnie nie tego chciałbym mieć.

Kiedy ładujesz plik .csv, możesz wybrać, jaka ma być szerokość i długość geograficzna, czy istnieje sposób, aby zrobić to samo dla rastrów?

Ponieważ mam raster 193, świetnie byłoby zasugerować polecenie, które można wykonać iteracyjnie.

wprowadź opis zdjęcia tutaj

Zdjęcie pokazuje mapę Indii i pierwsze pasmo pliku 3B31, tj. Opady powierzchniowe. Dwie czarne linie to równoległe linie 40 i -40. Chciałbym, aby były horyzontalne, a zatem dane dotyczące opadów pokrywają się z Indiami.

Paolo Santini
źródło
Bardzo przydatny byłby obraz testowy.
user30184,
@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download zawiera wiele linków.
AndreJ

Odpowiedzi:

7

Próbowałem pliku hdf 3b31 z mirador i znalazłem taką samą orientację, jak napotkałeś. Utworzyłem więc GCP, aby wykonać obrót w lewo:

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

z tym wynikiem:

wprowadź opis zdjęcia tutaj

Moja rotacja powinna być zgodna ze znacznikiem metadanych Origin = SOUTHWEST.

AndreJ
źródło
Właśnie tego szukałem, chociaż nie jestem w stanie prawidłowo używać poleceń, które napisałeś. Czy możesz wyjaśnić ich użycie bardziej szczegółowo?
Paolo Santini,
Umieść je w pliku tekstowym o nazwie rotate.batw tym samym folderze, co plik HDF. Sprawdź nazwę pliku źródłowego, uruchom powłokę OSGEO4W, przejdź do folderu z cd ... i uruchom partię. Dobra stara nawigacja w systemie DOS jest zawsze przydatna ;-)
AndreJ
3

Pobrałem 3B42RT.2016031806.7.7day.tif i .tfw, ale nie widzę nic dziwnego w georeferencji.

wprowadź opis zdjęcia tutaj

Raport Gdalinfo tego obrazu jest

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
użytkownik30184
źródło
Wygląda na to, że pliki tif 3b42 mają poprawną georeferencję, podczas gdy pliki hb 3b31 nie.
AndreJ
Wydaje się raczej, że HDF jest skomplikowany i nie jest w pełni obsługiwany przez GDAL trac.osgeo.org/gdal/wiki/HDF .
user30184,
Dokładnie niektóre z tych, które znajdziesz w Mirador, są georeferencyjne, wybierając właściwości NetCDF w oknie dostosowywania, inne nie.
Paolo Santini,
@ user30184 to do dostawcy danych należy dołączenie odpowiednich tagów metadanych georeferencji. Dostajesz to, za co zapłaciłeś ...
AndreJ