[GRASS-user] TRMM (3B43) handling

[nikos]

Interesting results, read on


On Sun, 2010-11-14 at 16:35 +0100, Markus Neteler wrote:
> On Sun, Nov 14, 2010 at 4:01 PM, nikos <vesnikos at gmail.com> wrote:
> > Hey list,
> >
> >
> > I'm battling to import the 43b3 Product (TRMM perspiration + other
> 
> (http://disc.sci.gsfc.nasa.gov/precipitation/documentation/TRMM_README/TRMM_3B42_readme.shtml)
> 
> > sources)  for more than 5 hours into grass without any success.
> >
> > gdalinfo doesnt list any list any subsets :(
> 
> Does
> gdalinfo --formats
> list the HDF driver?
> 
$ gdalinfo --version
GDAL 1.7.2, released 2010/04/23

$ gdalinfo --formats | grep HDF
  HDF4 (ro): Hierarchical Data Format Release 4
  HDF4Image (rw+): HDF4 Dataset
  HDF5 (ro): Hierarchical Data Format Release 5
  HDF5Image (ro): HDF5 Datase

> For
>  ftp://pps.gsfc.nasa.gov/trmmdata/ByDate/V06/2006/10/01/3B43.061001.6.HDF.gz
> I get this:
> 
> gdalinfo 3B43.061001.6.HDF
> Driver: HDF4/Hierarchical Data Format Release 4
> Files: 3B43.061001.6.HDF
> Size is 512, 512
> Coordinate System is `'
> Corner Coordinates:
> Upper Left  (    0.0,    0.0)
> Lower Left  (    0.0,  512.0)
> Upper Right (  512.0,    0.0)
> Lower Right (  512.0,  512.0)
> Center      (  256.0,  256.0)
> 
> gdalwarp 3B43.061001.6.HDF test.tif
> Input file 3B43.061001.6.HDF has no raster bands.
> 
> This indicates that GDAL 1.7 does not support this HDF subformat of TRMM.
> 
> I see right now this:
> http://trac.osgeo.org/gdal/ticket/3597
> 
> Apparently fixed in GDAL 1.8? Can you update and report here?


GDAL 1.8dev, released 2010/01/19


./gdalinfo  ~/TRMM-3B43/3B43.060801.6.HDF
Driver: HDF4/Hierarchical Data Format Release 4
Files: /home/nikos/greece-hdf/TRMM-3B43/3B43.060801.6.HDF
Size is 512, 512
Coordinate System is `'
Subdatasets:

SUBDATASET_1_NAME=HDF4_SDS:UNKNOWN:"/home/nikos/greece-hdf/TRMM-3B43/3B43.060801.6.HDF":0
  SUBDATASET_1_DESC=[1x1440x400] precipitation (32-bit floating-point)

SUBDATASET_2_NAME=HDF4_SDS:UNKNOWN:"/home/nikos/greece-hdf/TRMM-3B43/3B43.060801.6.HDF":1
  SUBDATASET_2_DESC=[1x1440x400] relativeError (32-bit floating-point)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

----
works

gdal version 1.7.3 doesn't incorporate the fix (based on the release
log) so if anyones reading this should go for 1.8.0+ 

> Markus

Now that I ID the sbs, 

I tried 

./gdal_translate -a_ullr  -180 50 -50 180
HDF4_SDS:UNKNOWN:"/DIR/OF/TRMM_IMAGE/3B43.060801.6.HDF":1 test.tif

which makes an interesting tif, that is not empty, but when I load it
with qgis (use psedocolor) its square(?). The image should have a
1440/400 ratio.

Best regards,

Nikos Ves