Make subset extraction more efficent for large datasets, using bounding box

[astrofrog]

For large cubes, extracting a small rectangular region in image space can be computationally expensive when computing e.g. a profile or histogram because the full mask is computed. We should ideally be able to use the bounding box of selections to first extract a sub-region from the data first, at least for selections where all attributes are pixel components.