[FR]: Consider adding `RasterIndex` & `CRSIndex` by default to returned Xarray objects

### Feature Summary

I think it makes sense for Xee to automatically add a RasterIndex (https://rasterix.readthedocs.io/en/latest/raster_index/creating.html) where it makes sense (nearly always?) to the returned Xarray object.

You could also add a `CRSIndex` from xproj: https://xproj.readthedocs.io/en/latest/

### Use Cases

Easier handling of raster data. By using RasterIndex and CRSIndex, many classes of bugs are avoided by default