read gdal raster
adapted doc