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 and coordinate systems must be the same.