Isotropic Semivariogram in pykrige

[gustavowillam]

I have been using pykrige and would like to know how to calculate an isotropic semivariogram. In the pykrige.ok.OrdinaryKriging, OrdinaryKriging method, the parameter (anisotropy_angle), when using 0, calculates the isotropic semivariogram or an anisotropic semivariogram with an angle 0 degrees ?

[MuellerSeb]

To calculate a Semi-variogram, you could use GSTools, which is compatible with pykrige.

Here you can find a tutorial on how to estimate and fit an isotropic semi-variogram.

The variogram can be used with pykrige afterwards in the following way:

cov_model = gstools.Gaussian(dim=2)
cov_model.fit_variogram(bin_center, gamma)  # see the tutorial
OK = OrdinaryKriging(data, **cov_model.pykrige_kwargs)

The interface to pykrige will be updated with PR #125

[MuellerSeb]

Some clarification:
The anisotropy_angle_* are only useful/meaningful when anisotropy_scaling_* is given.
The anisotropy_scaling_* define the ratios between the length-scale in main direction and the length scales in transversal directions. If there is no anisotropy, the rotation angles are meaningless.

[MuellerSeb]

Closing due to inactivity. Feel free to re-open.