def read_tif(fname): src = gdal.Open(fname, gdal.GA_Update) pdem = src.GetRasterBand(1) gt = src.GetGeoTransform() image = pdem.ReadAsArray() return [gt,image]