点列表上的 2D 插值 Python

Question

我有一个插值问题。它不应该太复杂，但我找不到任何有效的解决方案。

我正在研究 2D 机翼截面，我知道机翼的每个坐标点 (x_inf,y_inf) 上的压力（标量），并且我想将此压力内插到机翼的另一个离散化点（(x_profil ,y_profile) 坐标）。这是我拥有的配置文件的图片（（x_inf，y_inf），绿色）和我想要插值的网格（（x_profil，y_profil），红色）。 Points for interpolation。 数据 x_inf、y_inf、p_inf 是相同大小的 numpy 数组（它们是从特定文件中提取的）。 x_profil, y_profil 也是相同大小的 numpy 数组（但与 _inf 数据不同）。

我首先尝试了 interp2d 函数，但结果是一个大小为插值点大小平方的数组。

pressure=interp2d(x_inf,y_inf,p_inf)
p_profil=pressure(x_profil,y_profil)

我还尝试使用 interp1d 仅在 x 轴上进行插值，但这也不起作用。将插值类型设置为“最近”或“零”有效，但插值中有“孔”，如下图所示：presure interpolation，其中绿色点是输入数据，红色点是插值.

pressure=interp1d(x_inf,p_inf,kind='zero',bounds_error=False)
p_profil=pressure(x_profil)

我正在使用 Python 2.7.10、scipy 0.16.1 和 numpy 1.9.2，在带有 Enthought Canopy 的 Windows 7 上运行。

有人知道我该如何解决我的问题吗？

提前非常感谢！

Answer 1

听起来好像是 interp2d 函数认为您正在将其传递给常规网格。来自 numpy 文档：

> 定义数据点坐标的数组。 如果点位于规则网格上，x 可以指定列坐标，y 可以指定行坐标，例如：

x, y : array_like x = [0,1,2]; y = [0,3]； z = [[1,2,3], [4,5,6]] 否则，x 和 y 必须指定每个点的完整坐标，例如：

x = [0,1,2,0,1,2]; y = [0,0,0,3,3,3]； z = [1,2,3,4,5,6] 如果 x 和 y 是多维的，则在使用前将它们展平。

z : array_like 在数据点处插值的函数值。如果 z 是多维数组，则在使用前将其展平。如果 x 和 y 指定列和行坐标，则展平 z 数组的长度为 len(x)*len(y)，如果 x 和 y 指定，则为 len(z) == len(x) == len(y)每个点的坐标。

听起来你需要强制 x 和 y 明确表示完整坐标