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.
źródło
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:
z tym wynikiem:
Moja rotacja powinna być zgodna ze znacznikiem metadanych Origin = SOUTHWEST.
źródło
rotate.bat
w 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 ;-)Pobrałem 3B42RT.2016031806.7.7day.tif i .tfw, ale nie widzę nic dziwnego w georeferencji.
Raport Gdalinfo tego obrazu jest
źródło