I’m experimenting with something like xarray for geospatial data.
The idea is to make a common wrapper for various existing spatial data packages with a very minimal interface.
It’s still very much work in progress but does some things xarrays does, like take the mean (or whatever) of a 3d Lat/Lon/Time dataset over the time dimension:
plot(mean(g; dims=Time()))
