Построение результата интерполяции - PullRequest
0 голосов
/ 04 января 2019

У меня есть файл h5, содержащий данные регулярной сетки.Я использовал код, по которому я могу легко получить интерполированное значение для трех заданных значений.Я использовал функцию RegularGridInterpolator для интерполяции.Теперь я хочу составить график, чтобы проверить правильность интерполяции или нет.Но я не понимаю, как я могу это сделать.Может кто-нибудь помочь мне сделать это, пожалуйста?Вот мой код:

import numpy as np   
import h5py
from scipy.interpolate import RegularGridInterpolator
f = h5py.File('file.h5', 'r')  
list(f.keys())
dset = f[u'data']
dset.shape
dset.value.shape
dset[0:63,0:63,0:63]
x = np.linspace(-10, 320, 64) 
y = np.linspace(-10, 320, 64)
z = np.linspace(-10, 320, 64)
my_interpolating_function = RegularGridInterpolator((x, y, z), dset.value)
pts = np.array([7.36970468e-09, -4.54271563e-09, 1.51802701e-09])
my_interpolating_function(pts)

Результатом интерполяции является массив ([5.45534467e-10])

...