[GIS] Reading JAXA GSMaP data in QGIS

gdaljaxaqgis

I'm trying to process the daily rainfall data provided by Jaxa

The documentation says that the format of the dat file (say: gsmap_nrt.20170514.0.25d.daily.00Z-23Z.dat) is: 4-byte float plain binary,little-endian

How can I read/process this data in QGIS or GDAL?

Best Answer

JAXA provides a format converter which can convert HDF5/HDF4 and their binary products into either GeoTiff or KMZ. (NB: only for Windows OS).

Just in case the above direct link to the tool does not work, this is JAXA Toolkit website and the first one (highlighted) is the same link.

Related Solutions

[GIS] Reading palsar k&c 50 data in GDAL or other open source application

Not sure if the below is suitable. It's free (but not opensource as far as I can tell):

http://earth.eo.esa.int/polsarpro/about.html

It supports

ALOS PALSAR: dual pol and quad pol (level 1.1 and 1.5 - JAXA CEOS and ERSDAC-Vexcel)

It includes a tool to:

Create BMP files from binary data (modulus, phase, real part, imaginary part) with different colourmaps (grey, jet, hsv);

From bitmap you can then convert to anything.

You could also send an mail to the GDAL list - http://lists.osgeo.org/mailman/listinfo/gdal-dev

[GIS] Reading data from raster band with Java and GDAL

There are a few ways to read the data. Forgive me, but I am not familiar with the java bindings and syntax. One method is to use GDALRasterBand->ReadBlock():

http://gdal.org/classGDALRasterBand.html#a09e1d83971ddff0b43deffd54ef25eef

That reads in the internal block defined by the driver, which is one scanline for the aaigrid driver. See the example in the documentation. To do this manually, similar to your example, use GDALRasterBand->RasterIO(), which appears to be wrapped in java with ReadRaster(). Loop over the raster GDALDataset->GetRasterYSize() times and read in one line:

double[]data = new double[nXSize];
for(int i = 0; i < nYSize; i++)
{
    band.ReadRaster(0, i, nXsize, 1, data);
    //do something with your data
}

You can change the size of your input window for whatever you would like to get from your data. For example, as above, use a 3x3 window on the ReadRaster() call to calculate slope and aspect. I just found a link to the javadoc here:

http://gdal.org/java/

Best Answer

Related Solutions

[GIS] Reading palsar k&c 50 data in GDAL or other open source application

[GIS] Reading data from raster band with Java and GDAL

Related Question