二维直方图的样条插值 - Python

Question

我有一个 2D 直方图（x,y 坐标和“权重”），我需要进行样条插值并提取数据（所以不仅仅是图形），我如何用 python 做到这一点？ （我在这里附上历史）谢谢！ LightMap X_Z

Answer 1

你可以使用scipy，我从here中拿了下面的例子，稍作修改：

from scipy import interpolate
import matplotlib.pyplot as plt

# define range of x values
x = np.arange(-5.01, 5.01, 0.25)
# define range of z values
z = np.arange(-5.01, 5.01, 0.25)
# create a meshgrid
xx, zz = np.meshgrid(x, z)
# these weights would be given to you in your histogram - here is an example function to generate weights
weights = np.sin(xx**2+zz**2)
# create interpolation object
f = interpolate.interp2d(x, z, weights, kind='cubic')

# generate new ranges of x and z values
xnew = np.arange(-5.01, 5.01, 1e-2)
znew = np.arange(-5.01, 5.01, 1e-2)
# interpolate 
weightsnew = f(xnew, znew)

# plot
plt.imshow(weightsnew)
plt.show()