Dear all: I tried it with every kind of gdal binaries available: from Christoph Gohlke, from Tamas Szekeres, and from Osgeo4w.
But all the gdal binaries report wrong dimension of the HDF dataset. What should I do? Is the problem inherent to gdal source code? Any of yours idea is appreciated. Thanks. Lucia ---------- Forwarded message ---------- From: Lucia Leon <leoneslu...@gmail.com> Date: Thu, Apr 23, 2015 at 8:52 PM Subject: Unable to read data correctly in gdal python To: gdal-...@lists.osgeo.org 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
_______________________________________________ osgeo4w-dev mailing list osgeo4w-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/osgeo4w-dev