python 三維散點插值 griddata

#三維點插值
#在三維空間中，利用實際點的值推算出網格點的值
import numpy as np 

point_grid =np.array([[0.0,0.0,0.0],[0.4,0.4,0.4],[0.8,0.8,0.8],[1.0,1.0,1.0]])#網格點坐標

def func(x, y, z):
    return x*(1-x)*np.cos(4*np.pi*x) * (np.sin(4*np.pi*y**2)**2)*z

points = np.random.rand(10, 3)#實際點坐標
values = func(points[:,0], points[:,1],points[:,2])#實際點的值

from scipy.interpolate import griddata
grid_z0 = griddata(points, values, point_grid, method='nearest')#插值計算，計算出網格點的值