Jose Gomez-Dans wrote:
Running gdalinfo -stats on the new filename proceeds reporting the
expected information, band by band... slowly but surely, until it
reaches band 33 (it's band 33 everytime, even with different datasets):
Band 33 Block=128x128 Type=Int16, ColorInterp=Undefined
ERROR 4: `HDF4_EOS:EOS_GRID:"../BlahBlah.hdf":blah:blah' does not exist
in the file system,
and is not recognised as a supported dataset
name.
ERROR 1: IReadBlock failed at X offset 0, Y offset 0
I can read individual bands OK (although I haven't tried looping through
them), including band 33. I have stacks of memory, and data size should
not be an issue.
Jose,
I suspect the HDF library is using a table of file accessors that is only
32 large,so when GDAL tries to open a 33rd file things fail. Your workaround
of translating to geotiff and working with that should be fine. If you really
wanted to fix this you could likely modify a constant in the hdf library for
the table size and recompile it but that could be somewhat involved.
Best regards,
--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up | Frank Warmerdam, [email protected]
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush | Geospatial Programmer for Rent
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
