格点插值为站点数据批量处理

lixing123

学习中学习中学习中学习中

Jackluo

谢谢楼主分享~

刘思然

大家有没有打开grib文件插值到站点或者是nc文件插值到站点的呢？有的贴出来参考下呀，试了好久都不对

MeteoInfo

刘思然 发表于 2017-12-6 21:45
大家有没有打开grib文件插值到站点或者是nc文件插值到站点的呢？有的贴出来参考下呀，试了好久都不对

参考此贴一楼的脚本

Luo丹王道

王老师，您好！我用温度的NC数据插值，运行都是正常的，可是就是得出的结果都是负的多位数或者正的多位数，这是什么原因呢？

Luo丹王道

我的运行脚本如下：
#Open station file as table data
stdata = readtable('E:/Temp/fujianstations.csv', delimiter=',', format='%s%f%f')
x = stdata['Longitude']
y = stdata['Latitude']
#Add data file
f = addfile('E:/gradsfj/max24temp.ctl')
#Get data variable
psv = f['tem']
#Get time dimension length
tn = psv.dimlen(0)
#Loop
for i in range(0, tn):
    #Get dimension array
    ps = psv[i,(20,30),(110,125)]
    #Interpolate to stations
    ps_st = ps.tostation(x, y)
    #Add column to table data
    colname = 'tem' + str(i)
    stdata.addcol(colname, '%d', ps_st)
#Save table data to a file
fn = 'E:/Temp/t4.csv'
stdata.savefile(fn)
print 'Finish...'

MeteoInfo

Luo丹王道 发表于 2017-12-24 06:38
我的运行脚本如下：
#Open station file as table data
stdata = readtable('E:/Temp/fujianstations.csv ...

请仔细看看本贴一楼的代码。

Luo丹王道

MeteoInfo 发表于 2017-12-24 20:25
请仔细看看本贴一楼的代码。

可以了，NC文件可以直接打开~

qiaody

感谢楼主！非常有用

零陵路

感谢楼主！

　　大连妇科医院哪家较好 http://www.0411bh.com

　　大连男科医院哪家好 http://www.lnbohaink.com