Hi- Sorry for stupidly ignoring gdal on the gdal list. Brain cramp. In any case, yes I would recommend this. from osgeo import gdal g = gdal.Open("test.tif") dataarray = g.ReadAsArray() ...