Python多维插值

callmezoro

以前一直用的python的interp2d插值，这次工作需要一个四维插值，找了一下scipy的interpolate列表，发现一个interpnd的module，打算试一试，但是遇到了报错，我首先自己排查了一下并没有命名为interpnd的函数或者module
之后呢我又help了一下这个函数，并和我常用的interp2d进行了一下对比，结果如图，还是没有找到解决方案我就去看了一下那个interpnd的文件，全是乱码

不知道这个具体问题应该怎么解决

附件是我的数据，程序语句是：
from scipy import interpolate
f=interpolate.interpnd(lon,lat,depth,time,rho,kind='linear')

实际上我用别的三维数据用这个interpnd还是会报错，有没有解决办法啊或者有其他的多维插值的办法

callmezoro

我最后用循环写的二维插值，简化了一部分计算

岫鹿森涌

我也是四维数组插值，写的循环不知道对不对

1. from scipy import interpolate
   
2. d2=np.zeros((684,19,35,180),dtype=float)
   
3. for i in range(0,684):
   
4.     for j in range(0,19):
   
5.         f=interpolate.interp2d(x1,y1,zg11[i,j,:,:],kind='linear')
   
6.     #x2和y2是待插值的经纬度，现在用0.5的分辨率组成的网格
   
7.         lon,lat1 =  np.arange(0,360,2),np.arange(20,90,2)
   
8.     #直接使用f函数把x2和y2上的格点插值出来
   
9.         d2[i,j,:,:] = f(lon,lat1)
   
10.     #输出结果看看
    
11. print(d2.shape)

复制代码