¿Lectura de datos GSMaP de JAXA en QGIS?

Question

Estoy intentando procesar los datos de precipitaciones diarias proporcionados por Jaxa

En documentación dice que el formato del archivo dat (digamos: gsmap_nrt.20170514.0.25d.daily.00Z-23Z.dat) es: 4-byte float plain binary,little-endian

¿Cómo puedo leer/procesar estos datos en QGIS o GDAL?

Answer 1

JAXA proporciona un conversor de formato que puede convertir HDF5/HDF4 y sus productos binarios en GeoTiff o KMZ. (Nota: solo para sistemas operativos Windows).

En caso de que el enlace directo a la herramienta no funcione, es el siguiente Sitio web del JAXA Toolkit y el primero (resaltado) es el mismo enlace.

Answer 2

Puede utilizar el controlador vrt de GDAL para leer los datos de origen con un archivo denominado test.vrt así:

<VRTDataset rasterXSize="1440" rasterYSize="480">
  <VRTRasterBand dataType="CFloat32" band="1" subClass="VRTRawRasterBand">
    <SourceFilename relativetoVRT="1">test.dat</SourceFilename>
    <ImageOffset>0</ImageOffset>
    <PixelOffset>4</PixelOffset>
    <LineOffset>5760</LineOffset>
  </VRTRasterBand>
</VRTDataset>

y estos comandos GDAL:

gdal_translate -ot Float32 -of GTiff -a_srs EPSG:4326 -a_ullr 0 60 360 -60 -a_nodata -999.90002441406 test.vrt vrt1.tif
gdalwarp -t_srs WGS84 vrt1.tif vrt2.tif -wo SOURCE_EXTRA=10000 --config CENTER_LONG 0

Esto funciona en todos los sistemas operativos para archivos de acumulación de 24h de 0,25 deg de tamaño de celda. Es necesario adaptar el rasterXSize y rasterYSize para otros archivos de acuerdo con la documentación, así como volver a calcular el LineOffset à rasterXSize * 4

El resultado se puede cargar en QGIS, y tiene el mismo aspecto que el tif creado por la herramienta de conversión de formatos JAXA (probado con datos de 2017-05-01):