Apply spatial mask across non-spatial dimensions

[colinahill]

odc.geo.xr.mask fails for xarray.DataArray and xarray.Dataset if they have more than the 2 spatial dimensions of latitude and longitude (or equivalent). As other dimensions are common (e.g. time) it would be convenient to allow spatial masking across all other (non-spatial) dimensions.

This also impacts odc.geo.xr.crop, which uses mask.

[robbibt]

Hye @colinahill, that seems problematic - I'm pretty sure I tested this on multi-temporal data originally. Do you have some simple reproducible example I could try this out on?

Edit: we have code below that successfully applies a mask to multitemporal data, so if this no longer works, it's broken recently. What specific error do you get here?
https://knowledge.dea.ga.gov.au/notebooks/How_to_guides/Polygon_drill.html#Mask-data-using-polygon