Add support for reading a subset of raster cells to from_gdal

[kordejong]

This is convenient for performing scalability experiments. Something like this:

array = from_gdal(name, center, shape, partition_shape);
array = from_gdal(name, south_west, north_east, partition_shape);