[gdal-dev] Unable to read data correctly in gdal python

[Lucia Leon]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20150423/a7ff952f/attachment.html>