XArray

Seasonal Grouping

Extract Time Series (from Location)

time_series = data.isel(x=1000, y=1000).to_pandas().dropna()

Mean Across Multiple Dimensions

data.mean(dim=['lat', 'lon'])