小编典典

大矩阵的SciPy插值

python

我有一个ndarray(Z),在矩形网格(X,Y)上有一些500000个元素。

现在,我想在x,y的约100个位置插值,这些值不一定在网格上。

我有一些在Matlab中工作的代码:

data = interp2(X,Y,Z, x,y);

但是,当我尝试对scipy.interpolate使用相同的方法时,根据方法的不同会出现各种错误。例如,如果我指定interp2d失败,则会出现MemoryError,如果我指定kind = 'linear',则会出现“ OverflowError:要插入的数据点过多”
kind='cubic'。我也尝试过Rbfbisplev但他们也失败了。

所以问题是:是否有一个插值函数可以对大矩阵进行插值?有其他解决方法吗?(或者我是否必须编写一个函数来选择要插入点周围的合适区域并进行调用然后插入interp2d?)

另外:如何处理复数?


阅读 222

收藏
2021-01-20

共1个答案

小编典典

由于数据在网格上,因此可以使用RectBivariateSpline

与复数的工作,你可以插值data.realdata.imag独立(在FITPACK程序IIRC不处理复杂的数据)。

2021-01-20