pykrige.rk.Krige

class pykrige.rk.Krige(method='ordinary', variogram_model='linear', nlags=6, weight=False, n_closest_points=10, verbose=False)[source]

A scikit-learn wrapper class for Ordinary and Universal Kriging. This works with both Grid/RandomSearchCv for finding the best Krige parameters combination for a problem.

__init__(method='ordinary', variogram_model='linear', nlags=6, weight=False, n_closest_points=10, verbose=False)[source]

Initialize self. See help(type(self)) for accurate signature.

Methods

__init__([method, variogram_model, nlags, …]) Initialize self.
execute(points, *args, **kwargs)
param points:
fit(x, y, *args, **kwargs)
param x:array of Points, (x, y) pairs of shape (N, 2) for 2d kriging
get_params([deep]) Get parameters for this estimator.
predict(x, *args, **kwargs)
param x:array of Points, (x, y) pairs of shape (N, 2) for 2d kriging
score(X, y[, sample_weight]) Return the coefficient of determination R^2 of the prediction.
set_params(**params) Set the parameters of this estimator.
execute(points, *args, **kwargs)[source]
Parameters:
  • points (dict) –
  • Returns
  • -------
  • array (Variance) –
  • array
fit(x, y, *args, **kwargs)[source]
Parameters:
  • x (ndarray) – array of Points, (x, y) pairs of shape (N, 2) for 2d kriging array of Points, (x, y, z) pairs of shape (N, 3) for 3d kriging
  • y (ndarray) – array of targets (N, )
get_params(deep=True)[source]

Get parameters for this estimator.

Parameters:deep (bool, default=True) – If True, will return the parameters for this estimator and contained subobjects that are estimators.
Returns:params – Parameter names mapped to their values.
Return type:mapping of string to any
predict(x, *args, **kwargs)[source]
Parameters:x (ndarray) – array of Points, (x, y) pairs of shape (N, 2) for 2d kriging array of Points, (x, y, z) pairs of shape (N, 3) for 3d kriging
Returns:
Return type:Prediction array
score(X, y, sample_weight=None)[source]

Return the coefficient of determination R^2 of the prediction.

The coefficient R^2 is defined as (1 - u/v), where u is the residual sum of squares ((y_true - y_pred) ** 2).sum() and v is the total sum of squares ((y_true - y_true.mean()) ** 2).sum(). The best possible score is 1.0 and it can be negative (because the model can be arbitrarily worse). A constant model that always predicts the expected value of y, disregarding the input features, would get a R^2 score of 0.0.

Parameters:
  • X (array-like of shape (n_samples, n_features)) – Test samples. For some estimators this may be a precomputed kernel matrix or a list of generic objects instead, shape = (n_samples, n_samples_fitted), where n_samples_fitted is the number of samples used in the fitting for the estimator.
  • y (array-like of shape (n_samples,) or (n_samples, n_outputs)) – True values for X.
  • sample_weight (array-like of shape (n_samples,), default=None) – Sample weights.
Returns:

score – R^2 of self.predict(X) wrt. y.

Return type:

float

Notes

The R2 score used when calling score on a regressor will use multioutput='uniform_average' from version 0.23 to keep consistent with r2_score(). This will influence the score method of all the multioutput regressors (except for MultiOutputRegressor). To specify the default value manually and avoid the warning, please either call r2_score() directly or make a custom scorer with make_scorer() (the built-in scorer 'r2' uses multioutput='uniform_average').

set_params(**params)[source]

Set the parameters of this estimator.

The method works on simple estimators as well as on nested objects (such as pipelines). The latter have parameters of the form <component>__<parameter> so that it’s possible to update each component of a nested object.

Parameters:**params (dict) – Estimator parameters.
Returns:self – Estimator instance.
Return type:object