Estou trabalhando com os arquivos TRMM 3B31 sobre precipitação fornecidos pela NASA. Quando os carrego no QGIS, selecionando a projeção WGS84 junto com outros arquivos de forma, os mapas são exibidos verticalmente, ou seja, o QGIS interpreta como latitude, o que é longitude e vice-versa. Tentei a opção gdal_translate -a_ullr, selecionando a imposição de boas coordenadas, mas isso não resolve o problema: as imagens são compactadas e esticadas apenas para satisfazer meus requisitos, mas claramente não é o que eu gostaria de ter.
Ao carregar um arquivo .csv, você pode escolher qual deve ser a latitude e a longitude, existe uma maneira de fazer o mesmo com rasters?
Como tenho 193 raster, seria ótimo se você pudesse sugerir um comando que possa ser executado iterativamente.
A imagem mostra um mapa da Índia e a primeira banda do arquivo 3B31, ou seja, precipitação superficial. As duas linhas pretas são as paralelas 40 e -40. Gostaria de tê-los na horizontal e, portanto, os dados de chuva se sobrepõem à Índia.
fonte
Respostas:
Eu tentei o arquivo hb 3b31 do mirador e encontrei a mesma orientação que você encontrou. Então, criei um GCP para executar uma rotação no sentido anti-horário:
com este resultado:
Minha rotação deve seguir a marca de metadados Origin = SOUTHWEST.
fonte
rotate.bat
na mesma pasta que o arquivo HDF. Verifique o nome do arquivo de origem, inicie o OSGEO4W Shell, navegue até a pasta com cd ... e execute o lote. O bom e velho navegação DOS é sempre ;-) útilEu baixei o 3B42RT.2016031806.7.7day.tif e o .tfw, mas não vejo nada de estranho no georreferenciamento.
O relatório Gdalinfo desta imagem é
fonte