使用scipy.interp2d插值获全相同错误值，如何修复代码？

核心原因 #

你的代码出现Z值全相同的问题，本质是输入数据格式不符合interp2d的要求：interp2d仅支持两种输入模式——要么传入1D的x/y轴数组+对应二维网格的Z值，要么传入严格规则网格的二维X/Y数组（X每行重复x轴值，Y每列重复y轴值）。如果输入的X/Y是散点或非规则二维数组，interp2d会错误解析坐标，导致插值失效。

以下分两种场景给出修复代码：

场景1：原始数据是规则网格 #

如果你的3*3数据本身是规则网格点（比如x轴为[0,1,2]，y轴为[0,1,2]，对应网格坐标），只需将X/Y提取为1D轴数组，再传入interp2d：

import numpy as np
from scipy.interpolate import interp2d
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

# 定义规则网格的轴与对应Z值
x = np.array([0, 1, 2])
y = np.array([0, 1, 2])
X, Y = np.meshgrid(x, y)  # 生成二维网格坐标
Z = np.array([[1, 3, 2], [4, 6, 5], [7, 9, 8]])  # 自定义3*3 Z数据

# 正确创建插值函数：传入1D x/y轴 + 二维Z网格
f = interp2d(x, y, Z, kind='linear')

# 生成5*5的新插值点
x_new = np.linspace(x.min(), x.max(), 5)
y_new = np.linspace(y.min(), y.max(), 5)
Z_new = f(x_new, y_new)

# 绘制3D曲面
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
X_new, Y_new = np.meshgrid(x_new, y_new)
ax.plot_surface(X_new, Y_new, Z_new, cmap='viridis')
plt.show()

场景2：原始数据是散点（非规则网格） #

如果你的3*3数据是随机散点，interp2d不适用，改用CloughTocher2DInterpolator（专门处理散点的二维插值工具）：

import numpy as np
from scipy.interpolate import CloughTocher2DInterpolator
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

# 散点形式的3*3测试数据
X_scatter = np.array([0, 1, 2, 0.5, 1.5, 2.5, 1, 2, 3])
Y_scatter = np.array([0, 0, 0, 1, 1, 1, 2, 2, 2])
Z_scatter = np.array([1, 2, 3, 4, 5, 6, 7, 8, 9])

# 创建散点插值器
points = np.column_stack((X_scatter, Y_scatter))
f = CloughTocher2DInterpolator(points, Z_scatter)

# 生成5*5的新网格
x_new = np.linspace(X_scatter.min(), X_scatter.max(), 5)
y_new = np.linspace(Y_scatter.min(), Y_scatter.max(), 5)
X_new, Y_new = np.meshgrid(x_new, y_new)
Z_new = f(X_new, Y_new)

# 绘制3D曲面+原始散点对比
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(X_new, Y_new, Z_new, cmap='viridis')
ax.scatter(X_scatter, Y_scatter, Z_scatter, color='red', s=50)
plt.show()

额外推荐：使用官方替代工具 #

scipy已标记interp2d为过时，推荐用RegularGridInterpolator处理规则网格插值，用法更清晰：

from scipy.interpolate import RegularGridInterpolator

# 规则网格轴与Z值
x = np.array([0, 1, 2])
y = np.array([0, 1, 2])
Z = np.array([[1, 3, 2], [4, 6, 5], [7, 9, 8]])

# 创建插值器
f = RegularGridInterpolator((x, y), Z)

# 生成新网格坐标并计算插值
x_new = np.linspace(x.min(), x.max(), 5)
y_new = np.linspace(y.min(), y.max(), 5)
X_new, Y_new = np.meshgrid(x_new, y_new)
points_new = np.column_stack((X_new.ravel(), Y_new.ravel()))
Z_new = f(points_new).reshape(X_new.shape)

# 绘图逻辑同前
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(X_new, Y_new, Z_new, cmap='viridis')
plt.show()

内容的提问来源于stack exchange，提问作者M_Sea