我有一个包含常规网格数据的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)
插值的输出为array([5.45534467e-10])