我有一个网格化的天气数据集,其尺寸为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) . 那为什么会抛出一个错误?可以解释一下吗?我们将非常感谢您的帮助 .