线性插值的代码python_python中带有数据的线性插值

我有一个网格化的天气数据集，其尺寸为33 X 77 X 77.第一个维度是时间，其余分别是Lat和Lon . 我需要每次将数据插入(线性或最近邻居)到不同的点(lat＆lon)并将其写入csv文件 . 我已经使用了scipy的interp2d函数，它成功了一步 . 因为我有很多地方，所以我不想随着时间的推移而循环 .

下面显示的是我写的一段代码，任何人都可以提出一个更好的方法来完成任务吗？

import sys ; import numpy as np ; import scipy as sp ; from scipy.interpolate import interp2d ;import datetime ; import time ; import pygrib as pg ;

grb_f=pg.open('20150331/gfs.20150331.grb2') lat=tmp[0].data(lat1=4,lat2=42,lon1=64,lon2=102)[1] ; lat=lat[:,0];

lon=tmp[0].data(lat1=4,lat2=42,lon1=64,lon2=102)[2] ; lon=lon[0,:] ;

temp=np.empty((0,lon.shape[0]))

for i in range(0,tmp.shape[0]):

dat=tmp[i].data(lat1=4,lat2=42,lon1=64,lon2=102)

temp=np.concatenate([temp,dat[0]-273.15],axis=0)

temp1=temp.reshape(tmp.shape[0],lat.shape[0],lon.shape[0])

x=77 ; y=28 #(many points)

f=interp2d(lon,lat, temp1[0,:,:],kind='linear',copy=False,bounds_error=True ) ; Z=f(x,y)

编辑::

我没有制作3D矩阵，而是垂直添加数据并制作大小为2541 X 77且纬度为2541 X 1的lat和lon的数据矩阵.interp2d函数给出无效长度误差 .

f=interp2d(lon,lat, temp1[0,:,:],kind='linear',copy=False,bounds_error=True )

“非矩形网格的输入z的长度无效”)

ValueError：非矩形网格的输入z的长度无效

我的x，y，z矩阵的长度是相同的(2541,2541,2541) . 那为什么会抛出一个错误？可以解释一下吗？我们将非常感谢您的帮助 .