- 
                Notifications
    
You must be signed in to change notification settings  - Fork 197
 
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

