[Gdal-dev] Modis HDF help

[Andrey Kiselev]

On Mon, Sep 04, 2006 at 11:08:20PM -0300, Daniel Victoria wrote:
> I'm trying to import Modis BRDF Albedo (Mod43B) using gdal. The Modis
> web page tells it's a 22 band image (2 quality, 10 albedo for diferent
> bands and white/black body) and sure enough, issuing gdalinfo on the
> file will give me 2 subdatasets:
> 
> SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"MOD43B3.A2003193.h21v10.004.2003226054122.hdf":MOD_Grid_BRDF:Albedo
>  SUBDATASET_1_DESC=[1200x1200x10x2] Albedo (16-bit integer)
> SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:"MOD43B3.A2003193.h21v10.004.2003226054122.hdf":MOD_Grid_BRDF:Albedo_Quality
>  SUBDATASET_2_DESC=[1200x1200x2] Albedo_Quality (32-bit unsigned integer)
> 
> Where subdataset 1 has 10x2 bands of 1200 x 1200. But then, when I
> issue gdalinfo on the first subdataset:
> 
> gdalinfo 
> HDF4_EOS:EOS_GRID:"MOD43B3.A2003193.h21v10.004.2003226054122.hdf":MOD_Grid_BRDF:Albedo
> 
> I get 1200 bands of 1200x1 instead of 20 bands of 1200x1200? So, what
> am I doing wrong? Modis reprojection tool will list the 22 bands but
> then I don't know how to use MRT in bash mode to import 60 images in
> one swift stroke.


Daniel,

There is nothing wrong in your actions, but the HDF4 driver does not
proprly work with 4D datasets. Could you make available for download one
of those datasets?

Regards,
Andrey

-- 
Andrey V. Kiselev
ICQ# 26871517