Fazer o QGIS interpretar coordenadas como latinas longas em vez de latinas?

13

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.

insira a descrição da imagem aqui

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.

Paolo Santini
fonte
Uma imagem de teste seria muito útil.
user30184
@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download fornece muitos links.
Andrej

Respostas:

7

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:

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

com este resultado:

insira a descrição da imagem aqui

Minha rotação deve seguir a marca de metadados Origin = SOUTHWEST.

AndreJ
fonte
Era exatamente isso que eu estava procurando, embora não seja capaz de usar os comandos que você escreveu da maneira correta. Você pode explicar o uso deles em mais detalhes?
Paolo Santini
Coloque-os em um arquivo de texto nomeado rotate.batna 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 ;-) útil
Andrej
3

Eu baixei o 3B42RT.2016031806.7.7day.tif e o .tfw, mas não vejo nada de estranho no georreferenciamento.

insira a descrição da imagem aqui

O relatório Gdalinfo desta imagem é

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
user30184
fonte
Parece que os arquivos 3b42 tif têm georreferenciamento correto, enquanto os arquivos 3b31 hdf não.
Andrej
Parece que o HDF é complicado e não é perfeitamente suportado pelo GDAL trac.osgeo.org/gdal/wiki/HDF .
user30184
Exatamente alguns dos que você encontra no Mirador são georreferenciados, escolhendo as propriedades do NetCDF na janela de personalização, outros não.
Paolo Santini
@ user30184, cabe ao provedor de dados incluir tags de metadados de georreferenciamento adequadas. Você começa o que você pagou ...
Andrej