spatial_means function

Computes spatial empirical means using a dataframe or a stars object