python二维线性插值

在Python中，二维线性插值通常通过`scipy.interpolate`模块的`griddata`函数或者`interp2d`函数来实现。这两个函数允许你在给定的二维数据集上进行插值，当需要获取数据集中某个不在原点的点的值时。

`griddata`函数接受三个参数：输入坐标（通常是两个数组），对应的数据值，以及一个新的查询点网格。它会返回新的网格对应位置处的插值值。

例如：

from scipy.interpolate import griddata

# 假设我们有数据点 (x, y) 和它们对应的值 z
x = [0, 1, 2, 3]
y = [0, 1, 2, 3]
z = [1, 4, 9, 16]

# 创建一个新的网格，比如均匀间隔的 (0.5, 1.5)
xi, yi = np.mgrid[0:4:10j, 0:4:10j].reshape(2,-1)

# 进行二维线性插值
zi = griddata((x, y), z, (xi, yi), method='linear')

`interp2d`则创建了一个插值表面，并提供了一个类似`numpy`数组的功能，可以直接对网格应用插值：

import numpy as np
from scipy.interpolate import interp2d

# 同样的数据
f = interp2d(x, y, z, kind='linear')  # 'linear' 表示线性插值
zi = f(xi, yi)