Вернуть z значений 3D scipy.spatial. Delaunay

Question

У меня есть следующий код, который создает поверхность:

u = np.array(mat) 
v = np.array(stk) 

x = u
y = v
z = np.array(vol)

tri = Delaunay(np.array([u,v]).T)

Я хотел бы сравнить поверхность с моим полным набором данных и найдите СКО. Однако для этого мне нужно найти значения z точек на поверхности, которые не являются вершинами симплексов. Как мне go сделать это?

Answer 1

Вы можете использовать scipy.interpolate.LinearNDInterpolator.

Пример использования:

interpolator = LinearNDInterpolator(tri, z)
interpolator(200, 200)

или с массивом точек:

interpolator(np.array([x,y]).T)