hydrobox.geostat.ordinary_kriging

hydrobox.geostat.ordinary_kriging(variogram, grid_resolution, exact=True, cond_err='nugget', pseudo_inv=True, pseudo_inv_type='pinv', return_type='plot', **kwargs)

Use a scikit-gstat Variogram to estimate spatial properties of a sample. Uses this variogram to interpolate using kriging. The Kriging is done with the Ordinary class. Refer to the docs of Ordinary to learn about the parameters.

Note

Right now, there are only very limited possibilities to specify an interpolation grid. Only setting the resolution of the result like: 50x50 is possible by passing the integer 50. More fine-grained control will be added with a future release.

Parameters

  • variogram (skgstat.Variogram) – Variogram used for kriging

  • grid_resolution (int) – Resoultion of the interpolation grid. The resolution will be used in all input data dimensions, which can lead to non-quadratic grid cells.

  • exact (bool) – If True (default), the input data will be matched exactly. Refer to Ordinary for more info.

  • cond_err (str, float, list) – Measurement error, or variogram nugget. Refer to Ordinary for more info.

  • pseudo_inv (bool) – If True, the Kriging is more robust, but also slower. Refer to Ordinary for more info.

  • pseudo_inv_type (str) – Type of matrix inversion used if pseudo_inv is True. Refer to Ordinary for more info.

  • return_type (str) – Specify how the result should be retuned. Can be the kriging class itself ('object'), the interpolated grid ('grid') or a plot of the grid ('plot').

Return type

Union[List[ndarray], Figure]

Returns

  • results (numpy.ndarray, numpy.ndarray) – Interpolation grid and kriging error grid

  • fig (plotly.graph_objects.Figure, matplotlib.Figure) – Figure of the result plot.