Intersecting shapefiles and raster data at scale - good design patterns?

Thanks for sharing Kevin!

Our general recommendation is that packages like xagg shouldn’t do I/O at all. Consume and produce standard data containers and leave the I/O to xarray / pandas etc.