[Gdal-dev] Modis HDF help

[Daniel Victoria]

Hi all,

I bet someone already figure this problem out but I could not find the
solution in the mailing list archives so here goes the question.

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.

Thanks
Daniel