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 <[email protected]> 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 _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
