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

upuil

MeteoInfo 发表于 2021-2-27 08:54
addcol 改为  addcoldata

可以了，谢谢王老师！

kelly@kitty

请问如果是多个AWX文件，没有t，按文件名列表逐个循环打开文件并插值应该怎么弄啊

kelly@kitty

在08.csv'里存放了逐时的AWX文件名。
tn = readtable('F:\OCF\cloud\08.csv', delimiter=',', format='%s%s%s')
fn =tn['filename']
time=tn['time']
day=tn['day']
显示
java.io.FileNotFoundException: java.io.FileNotFoundException: Invalid file path

MeteoInfo

kelly@kitty 发表于 2021-4-14 17:20
在08.csv'里存放了逐时的AWX文件名。
tn = readtable('F:\OCF\cloud\08.csv', delimiter=',', format='%s% ...

路径分隔符用反斜杠

kelly@kitty

谢谢老师，程序完成了

Finnick

王老师您好，我在利用tostation插值的时候，想对多层进行插值，但是出来的数据都是一样的(不同时刻相同层次的数据是一样得)，我看前面只是在时间维批处理，不知道是否可以进行时间和不同高度层的批处理。

stdata=readtable('station.csv', delimiter=',')
t1=addfile("T.nc","r")
lev=t1['level'][::-1]
tdata=t1['t'][:,::-1,:,:]
time=t1['time'][17：239]
tn = time.dimlen()
ln=lev.dimlen()
tdata2=tdata[17:239,:,'0:60','74:130']
for i in range(0,tn):
    for j in range(0,ln):
        tem=tdata2[i,j,,'0:60','74:130']
        TEMP= tem.tostation(x, y)
        colname =str(j+1)
        print(j)
        stdata.addcoldata(colname,'%d', TEMP)