GDAL Dataset and Band abstraction for simple tiled raster calculations (AKA "map algebra").
Can handle rasters with different extents as long as they overlap.
Cellsizes must be the same.