change raster to readgdal