python有没有类似于ncl的linint2_points函数提取某个经纬度数据

搁浅

询问，python程序有没有向ncl的linint2_points函数一样，只要有经纬度信息和风速数据，根据想要提取那个经纬度风速即可提取的函数。或者如果没有的话，python是怎么实现这一功能的呢，求问，谢谢了

werewolf

interpolate.interpn

搁浅

werewolf 发表于 2019-4-23 09:17
interpolate.interpn

好的，多谢，我试一下

搁浅

spd = np.zeros((71,91))
for i in range(71):
    for j in range(91):
        spd[i,j] = struct.unpack('f',f[360 + (i * 91 + j) * 4 : 360 + (i * 91 + j + 1) * 4])[0]
lat = np.arange(-10.,61.,1)
lon = np.arange(60.,151.,1)
lat_,lon_ = np.meshgrid(lon,lat)
data = interpolate.interpn((lon,lat),spd,(116,38),method='linear')
spd是我的二维数据，我想插值提取（38,116）那个点的风速数据，这个函数应该怎么实现呢

搁浅

问题解决
data = interpolate.interpn((lat,lon),spd,(38,116),method='linear')，维度与spd维度相同即可。供其他有用到的人参考。