search

GeoStat-Framework / PyKrige

Kriging Toolkit for Python
https://pykrige.readthedocs.io
BSD 3-Clause "New" or "Revised" License
740 stars 186 forks source link

3D kriging takes very long time to execute #218

Closed askar-INTERA closed 2 years ago

askar-INTERA commented 2 years ago

Hi I am doing 3d ordinary kriging using 50 points that yield the values that I want to interpolate over 448 points only. This takes forever to be executed and I am not sure why, can you please help.

HK = pd.read_csv ('HK_3D.csv') PP = pd.read_csv ('PP_HK_3D.csv') parameters={'sill': 15.0, 'range': 30, 'nugget': 0.0} ok = OrdinaryKriging3D(PP['x'], PP['y'], PP['z'], PP['k'], variogram_model = "exponential", variogram_parameters=parameters, enable_plotting=True, verbose=True) k3d1, ss3d = ok.execute("grid", HK['x'], HK['y'], HK['z'])

MuellerSeb commented 2 years ago

That is strange. Have tried to do the same with GSTools?

I would also need some more information.