hydromt_wflow.workflows.spatial_interpolation#

hydromt_wflow.workflows.spatial_interpolation(forcing: ~xarray.core.dataarray.DataArray, ds_like: ~xarray.core.dataset.Dataset, interp_type: str, nnearest: int | None = 4, p: float | None = 2, remove_missing: bool | None = False, cov: str | None = '1.0 Exp(10000.)', src_drift: ~numpy.ndarray | None = None, trg_drift: ~numpy.ndarray | None = None, mask_name: str | None = None, logger: ~logging.Logger | None = <Logger hydromt_wflow.workflows.forcing (WARNING)>) → DataArray[source]#

Interpolate spatial forcing data from station observations to a regular grid.

This workflow uses the wradlib.ipol.interpolate function in wradlib. It wraps the following interpolation types into a single function: - “nearest”: nearest-neighbour interpolation, also works with a single station. - “idw”: inverse-distance weighing using 1 / distance ** p. - “linear”: linear interpolation using scipy.interpolate.LinearNDInterpolator, may result in missing values when station coverage is limited. - “ordinarykriging”: interpolate using Ordinary Kriging, see wradlib documentation for a full explanation: wradlib.ipol.OrdinaryKriging <https://docs.wradlib.org/en/latest/generated/wradlib.ipol.OrdinaryKriging.html> - “externaldriftkriging”: Kriging interpolation including an external drift, see wradlib documentation for a full explanation: wradlib.ipol.ExternalDriftKriging <https://docs.wradlib.org/en/latest/generated/wradlib.ipol.ExternalDriftKriging.html>

Parameters:

  • forcing (xr.DataArray) – GeoDataArray with the forcing data with time and index of the point data.

  • ds_like (xr.Dataset) – Target dataset defining the grid for interpolation.

  • interp_type (str) – Interpolation method. Options: “nearest”, “idw”, “linear”, “ordinarykriging”, “externaldriftkriging”.

  • nnearest (int, optional) – Maximum number of neighbors for interpolation. Default is 4.

  • p (float, optional) – Power parameter for IDW interpolation. Default is 2.

  • remove_missing (bool, optional) – Whether to mask NaN values in the input data. Default is False.

  • cov (str, optional) – Covariance model for Kriging. Default is ‘1.0 Exp(10000.)’.

  • src_drift (np.ndarray, optional) – External drift values at source points (stations).

  • trg_drift (np.ndarray, optional) – External drift values at target points (grid).

  • mask_name (str, optional) – Name of the mask variable in ds_like, by default None.

Returns:

Interpolated forcing data on the targeted grid.

Return type:

xr.DataArray

See also

wradlib.ipol.interpolate <https://docs.wradlib.org/en/latest/ipol.html#wradlib.ipol.interpolate>