New Xarray accessor for rasters and DEMs  #687

Hi all,

We are planning to build on top of rioxarray for end-users focusing on analysis of EO data by building accessor in our packages GeoUtils and xDEM for raster/vector and DEM analysis, respectively, which are currently built on top of rasterio and geopandas.
See: pydata/xarray#8040 and pydata/xarray#8041

I don't see much overlap for xDEM, except that we could use vertical CRS reprojection through rioxarray if one day we implement it in rasterio, for now not possible (see rasterio/rasterio#2433).

For GeoUtils, part of the package's functionalities are closely linked to georeferencing and aim to simplify the interaction raster/raster or raster/vector, like the reproject_match function of rioxarray. Those did not exist in rasterio, but they have been more of a focus in rioxarray and geocube. In the Raster and Vector classes of GeoUtils, we tried to implement these "convenience" ops for all type of operations and the two objects.

And so I am wondering: where do you think those functionalities should live if we build an xarray accessor? Should we keep them in GeoUtils like is done on top of rasterio (e.g., Raster.reproject(ref), Raster.polygonize(ref)) and implement ds.geo.reproject(ref) and ds.geo.polygonize(ref)) by simply wrapping what is done, or should we try to make those available directly in rioxarray and geocube (but then the functioning of rioxarray would stray quite a bit from that of rasterio)?

More details in our doc for match-reference ops: https://geoutils.readthedocs.io/en/stable/core_match_ref.html

Thanks for all the great work! 😄