Hello, The gdalinfo.exe correctly shows the subdataset of a HDF4 file with dimension [180x8x32x9]. But the same subdataset read from gdal in python is only [180x8x32].
I was using the gdal binary from Christoph Gohlke in windows, and according to him, the GDAL binaries are built with HDF-4.2.10: http://www.lfd.uci.edu/~gohlke/pythonlibs/6icuform/GDAL-1.11.2-cp27-none-win_amd64.whl The HDF4 file I was using is here: ftp://l5eil01.larc.nasa.gov/MISR/MIL2ASLS.002/2011.07.04/MISR_AM1_AS_LAND_P233_O061399_F07_0022.hdf The code I was using was as follows: import gdal from gdalconst import * data = gdal.Open('HDF4_EOS:EOS_GRID:' + infile + ':RegParamsLnd:ViewZenAng', GA_ReadOnly).ReadAsArray() print data.shape I hope some of you can help me. Lucia
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
