<div dir="ltr">Hello,<br><br>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].<br><br>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:<br><br><a href="http://www.lfd.uci.edu/~gohlke/pythonlibs/6icuform/GDAL-1.11.2-cp27-none-win_amd64.whl">http://www.lfd.uci.edu/~gohlke/pythonlibs/6icuform/GDAL-1.11.2-cp27-none-win_amd64.whl</a><br><br>The HDF4 file I was using is here:<br><br><a href="ftp://l5eil01.larc.nasa.gov/MISR/MIL2ASLS.002/2011.07.04/MISR_AM1_AS_LAND_P233_O061399_F07_0022.hdf">ftp://l5eil01.larc.nasa.gov/MISR/MIL2ASLS.002/2011.07.04/MISR_AM1_AS_LAND_P233_O061399_F07_0022.hdf</a><br><br>The code I was using was as follows:<br><br>import gdal<br>from gdalconst import *<br>data = gdal.Open('HDF4_EOS:EOS_GRID:' + infile + ':RegParamsLnd:ViewZenAng', GA_ReadOnly).ReadAsArray()<br>print data.shape<br><br>I hope some of you can help me.<br><br>Lucia</div>