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