With xarray-regrid it is possible to regrid between two rectilinear grids. The following methods are supported:
- Linear
- Nearest-neighbor
- Conservative
- Cubic
- "Most common value", as well as other zonal statistics (e.g., variance or median).
All regridding methods can operate lazily on Dask arrays.
Note that "Most common value" is designed to regrid categorical data to a coarse resolution. For regridding categorical data to a finer resolution, please use "nearest-neighbor" regridder.
Regridding is a common operation in earth science and other fields. While xarray does have some interpolation methods available, these are not always straightforward to use. Additionally, methods such as conservative regridding, or taking the most common value, are not available in xarray.
For a minimal install:
pip install xarray-regrid
To improve performance in certain cases:
pip install xarray-regrid[accel]
which includes optional extras such as:
dask
: parallelization over chunked datasparse
: for performing conservative regridding using sparse weight matricesopt-einsum
: optimized einsum routines used in conservative regridding
Benchmarking varies across different hardware specifications, but the inclusion of these extras can often provide significant speedups.
The xarray-regrid routines are accessed using the "regrid" accessor on an xarray Dataset:
import xarray_regrid
ds = xr.open_dataset("input_data.nc")
ds_grid = xr.open_dataset("target_grid.nc")
ds.regrid.linear(ds_grid)
For examples, see the benchmark notebooks and the demo notebooks.
The benchmark notebooks contain comparisons to more standard methods (CDO, xESMF).
To be able to run the notebooks, a conda environment is required (due to ESMF and CDO).
You can install this environment using the environment.yml
file in this repository.
Micromamba is a lightweight version of the much faster "mamba" conda alternative.
micromamba create -n environment_name -f environment.yml
This package was developed under Netherlands eScience Center grant NLESC.OEC.2022.017.