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

upuil · 发表于 2021-2-26 17:22:45

upuil 发表于 2021-2-26 17:20
王老师刚那个问题解决了，谢谢您，又出现一个新问题，麻烦您再看一下！

这里stdata.addcol(colname, '%d', ps_st) ，%d该怎么选呢，我试了%s，%f还是一样报错!

MeteoInfo · 发表于 2021-2-27 08:54:43

upuil 发表于 2021-2-26 17:22
这里stdata.addcol(colname, '%d', ps_st) ，%d该怎么选呢，我试了%s，%f还是一样报错!

addcol 改为 addcoldata

upuil · 发表于 2021-2-27 15:39:19

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

可以了，谢谢王老师！

kelly@kitty · 发表于 2021-4-14 11:41:23

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

kelly@kitty · 发表于 2021-4-14 17:20:24

在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 · 发表于 2021-4-14 21:02:27

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

路径分隔符用反斜杠

kelly@kitty · 发表于 2021-4-15 16:07:23

谢谢老师，程序完成了

Finnick · 发表于 2022-2-6 18:31:41

本帖最后由 Finnick 于 2022-2-6 19:45 编辑

王老师您好，我在利用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)

		自动登录	找回密码
密码			立即注册