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.
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.
fuente
Respuestas:
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:
con este resultado:
Mi rotación debe seguir la etiqueta de metadatos Origin = SOUTHWEST.
fuente
rotate.bat
en 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 ;-)Descargué 3B42RT.2016031806.7.7day.tif y .tfw pero no veo nada extraño en la georreferenciación.
El informe de Gdalinfo de esta imagen es
fuente