[GRASS-user] TRMM (3B43) handling

Markus Neteler neteler at osgeo.org
Sun Nov 14 14:41:47 EST 2010


On Sun, Nov 14, 2010 at 7:42 PM, nikos <vesnikos at gmail.com> wrote:
> 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)
...
>> 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


...excellent!

> 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+

(we could ask backport to 1.7)

>> 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.

Wait - AFAIK you need to use "gdalwarp" on it, not gdal_translate.

cheers
Markus


More information about the grass-user mailing list