[Gdal-dev] HDF EOS segmentation fault

[Jan Cermak]

Hi,

When I read in an hdfeos dataset from a MYD06 file I get a segfault of the following type:
*** glibc detected *** python: double free or corruption (fasttop): 0x081b8d90 ***
This occurs either at one of the steps following the data read or at the termination of the program. There does not seem to be a pattern. Here is some example code:

#!/usr/bin/env python
import gdal
DSname='HDF4_EOS:EOS_SWATH:"/home/jancer/data/tmp/MODIS/MYD06_L2.A2006220.1315.005.2006222153137.hdf":mod06:Cloud_Optical_Thickness'
DSobj = gdal.Open(DSname)
data = DSobj.ReadAsArray()
data.shape
halfdata = data/2.

There is a sample file at
ftp://ladsweb.nascom.nasa.gov/allData/5/MYD06_L2/2006/220/MYD06_L2.A2006220.1315.005.2006222153137.hdf


What's my mistake?

Thanks
Jan