我有一个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'。我也尝试过Rbf,bisplev但他们也失败了。
kind = 'linear'
kind='cubic'
Rbf
bisplev
所以问题是:是否有一个插值函数可以对大矩阵进行插值?有其他解决方法吗?(或者我是否必须编写一个函数来选择要插入点周围的合适区域并进行调用然后插入interp2d?)
另外:如何处理复数?
由于数据在网格上,因此可以使用RectBivariateSpline。
与复数的工作,你可以插值data.real和data.imag独立(在FITPACK程序IIRC不处理复杂的数据)。
data.real
data.imag