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

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.

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:

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:

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

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

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