Dask and Xarray have revolutionized the analysis of large scale geospatial data. They make it easy for researchers to intuitively query large multi-dimensional datasets, and then run those queries at large scale.

import xarray as xr
ds = xr.open_dataset(
hourly_wind = (
    .sel(lon=slice(180, None))
    .coarsen(lat=2, lon=2, boundary='trim')

More Examples

For more in-depth examples consider the following: