Import TRMM: Difference between revisions
(new, based on comments from Markus Metz) |
(added rotation trick from Markus Metz) |
||
Line 1: | Line 1: | ||
TRMM data are usually delivered in a specific HDF format which is (partially) supported from GDAL 1.8+ (so, {{cmd|r.in.gdal}} would in theory work). | TRMM data are usually delivered in a specific HDF format which is (partially) supported from GDAL 1.8+ (so, {{cmd|r.in.gdal}} would in theory work, see also below for rotation problem). | ||
However, it is easier to get binary grids from here | However, it is easier to get instead binary grids from here: | ||
ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/3B43_V6/ | ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/3B43_V6/ | ||
Line 15: | Line 15: | ||
Optionally convert then to mm/hr to some other time unit. | Optionally convert then to mm/hr to some other time unit. | ||
=== Getting rotated maps in TRMM-HDF? === | |||
The TRMM 3B43 HDF files are rotated by 90° (clockwise), so north is right and west is up. The included datasets "precipitation" and "error" can be dumped from the HDF files to a binary or ASCII file with the "hdp" tool which should be on any system that has the HDF4 library installed. | |||
Rotate "precipitation" with | |||
hdp dumpsds -r 4 -d -o <outfile> -b <TRMM_3B43.hdf> | |||
and "error" with | |||
hdp dumpsds -r 5 -d -o <outfile> -b <TRMM_3B43.hdf> | |||
These are rotated and need to be rotated back. The (properly oriented) binary grids are much easier to handle than the original HDF files. | |||
[[Category: FAQ]] | [[Category: FAQ]] |
Revision as of 11:54, 17 November 2010
TRMM data are usually delivered in a specific HDF format which is (partially) supported from GDAL 1.8+ (so, r.in.gdal would in theory work, see also below for rotation problem).
However, it is easier to get instead binary grids from here:
ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/3B43_V6/
The format is explained here. Hence, to import gridded TRMM into a Lat-long location with is done via
# generally r.in.bin -f -b north=50 south=-50 east=180 west=-180 rows=400 cols=1440 ....
e.g.
r.in.bin -f -b input=3B43.100101.6.precipitation.bin output=3B43.100101.6.precipitation \ north=50 south=-50 east=180 west=-180 rows=400 cols=1440
Optionally convert then to mm/hr to some other time unit.
Getting rotated maps in TRMM-HDF?
The TRMM 3B43 HDF files are rotated by 90° (clockwise), so north is right and west is up. The included datasets "precipitation" and "error" can be dumped from the HDF files to a binary or ASCII file with the "hdp" tool which should be on any system that has the HDF4 library installed.
Rotate "precipitation" with
hdp dumpsds -r 4 -d -o <outfile> -b <TRMM_3B43.hdf>
and "error" with
hdp dumpsds -r 5 -d -o <outfile> -b <TRMM_3B43.hdf>
These are rotated and need to be rotated back. The (properly oriented) binary grids are much easier to handle than the original HDF files.