-
Notifications
You must be signed in to change notification settings - Fork 195
Open
Description
Hi,
I'm trying to use OrdinaryKriging to extimate value, at frist i used style="grid" and it worked fine. Now I want to speed things up so that it doesnt calculate value for every point, but when using style="points" i get totaly different resoults.
from scipy.interpolate import griddata
gridx = np.arange(-150, 150, 0.5)
gridy = np.arange(-180, 150, 0.5)
z, ss = UK.execute("grid", gridx, gridy)
grid_x, grid_y = np.meshgrid(
np.linspace(-150, 150, z.shape[1]),
np.linspace(-180, 150, z.shape[0])
)
f = griddata(
(grid_x.flatten(), grid_y.flatten()),
z.flatten(),
(x, y),
method='linear'
)
But when i use
z, ss = UK.execute("points", x, y)
then I get:
I can't figure out what I'm doing wrong.
Metadata
Metadata
Assignees
Labels
No labels