¿Hacer que QGIS interprete las coordenadas como long-lat en lugar de lat-long?

13

Estoy trabajando con los archivos TRMM 3B31 sobre precipitación proporcionados por la NASA. Cuando los cargo en QGIS seleccionando la proyección WGS84 junto con otros archivos de formas, los mapas se muestran verticalmente, es decir, QGIS interpreta como latitud lo que en cambio es longitud y viceversa. He intentado la opción gdal_translate -a_ullr seleccionando imponer buenas coordenadas pero no resuelve el problema: las imágenes solo se comprimen y se estiran para satisfacer mis requisitos, pero esto claramente no es lo que me gustaría tener.

Cuando carga un archivo .csv, puede elegir cuál debe ser la latitud y la longitud, ¿hay alguna manera de hacer lo mismo para los rásteres?

Como tengo 193 ráster, sería genial si puedes sugerir un comando que se pueda ejecutar de forma iterativa.

ingrese la descripción de la imagen aquí

La imagen muestra un mapa de la India y la primera banda del archivo 3B31, es decir, precipitación superficial. Las dos líneas negras son las paralelas 40 y -40. Me gustaría tenerlos horizontales y, por lo tanto, que los datos de lluvia se superpongan con la India.

Paolo Santini
fuente
Una imagen de prueba sería muy útil.
user30184
@ user30184 gis.stackexchange.com/questions/89261/trmm-data-download ofrece muchos enlaces.
AndreJ

Respuestas:

7

Probé el archivo hdf 3b31 de mirador y encontré la misma orientación que encontraste. Así que creé algunos GCP para realizar una rotación en sentido antihorario:

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

con este resultado:

ingrese la descripción de la imagen aquí

Mi rotación debe seguir la etiqueta de metadatos Origin = SOUTHWEST.

AndreJ
fuente
Esto es exactamente lo que estaba buscando, aunque no puedo usar los comandos que escribió de la manera adecuada. ¿Puedes explicar su uso con más detalles?
Paolo Santini
Póngalos en un archivo de texto nombrado rotate.baten la misma carpeta que el archivo HDF. Verifique el nombre del archivo de origen, inicie el Shell OSGEO4W, navegue a la carpeta con cd ... y ejecute el lote. La buena navegación de DOS siempre es útil ;-)
AndreJ
3

Descargué 3B42RT.2016031806.7.7day.tif y .tfw pero no veo nada extraño en la georreferenciación.

ingrese la descripción de la imagen aquí

El informe de Gdalinfo de esta imagen es

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
usuario30184
fuente
Parece que los archivos tif 3b42 tienen georreferenciación correcta, mientras que los archivos hdf 3b31 no.
AndreJ
Parece más bien que HDF es complicado y no es perfectamente compatible con GDAL trac.osgeo.org/gdal/wiki/HDF .
user30184
Exactamente algunos de los que encuentra en Mirador están georreferenciados, eligiendo las propiedades de NetCDF en la ventana de personalización, otros no.
Paolo Santini
@ user30184 depende del proveedor de datos incluir etiquetas de metadatos de georreferenciación adecuadas.
Obtienes