C++ code to process hdf files with gdal