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

侠女黒丫

你好王老师，运行过程中，出现'NoneType' object has no attribute 'tostation'的错误，是什么原因呢？

侠女黒丫

侠女黒丫 发表于 2017-5-16 10:04
你好王老师，运行过程中，出现'NoneType' object has no attribute 'tostation'的错误，是什么原因呢？

用grd数据ctl插值就不会出现，用nc文件插值就会出现这种情况

MeteoInfo

侠女黒丫 发表于 2017-5-16 10:04
你好王老师，运行过程中，出现'NoneType' object has no attribute 'tostation'的错误，是什么原因呢？

可能是格点数据读取有误，没读出来。

fcb61723062

JMeteoInfo

wuyang0074

王老师，因为是刚刚接触，所以对一些基本的命令语言都不太懂，您能不能把每一行的命令是什么意思都给讲解一下呢，或者给一些学习的资料。谢谢您

MeteoInfo

wuyang0074 发表于 2017-7-7 10:39
王老师，因为是刚刚接触，所以对一些基本的命令语言都不太懂，您能不能把每一行的命令是什么意思都给讲解一 ...

参考这里：http://www.meteothinker.com/docs/meteoinfolab/user_guide.html

论坛也有很多示例脚本程序可以参考：http://bbs.06climate.com/forum.p ... &extra=page%3D1

wuyu

王老师，我运行了如下：
#Open station file as table data
stdata = readtable('F:/MeteoInfo/sample/LonLatStation/Stations.csv', delimiter=',', format='%s%f%f')
x = stdata['Longitude']
y = stdata['Latitude']
#Add data file
f = addfile('F:/MeteoInfo/sample/GrADS/taifenggrd.ctl')
#Get data variable
psv = f['UGRD10m']
#Get time dimension length
tn = psv.dimlen(0)
#Loop
for i in range(0, tn):
    #Get dimension array
    ps = psv[i,'10:50','90:140']
    #Interpolate to stations
    ps_st = ps.tostation(x, y)
    #Add column to table data
    colname = 'UGRD10m' + str(i)
    stdata.addcoldata(colname, '%d', ps_st)

#Save table data to a file
fn = 'f:/test_st.csv'
stdata.savefile(fn)
print 'Finish...'


提示错误为：
run script...
Traceback (most recent call last):
  File "<iostream>", line 6, in <module>
  File "F:\MeteoInfo\pylib\mipylib\dataset\midata.py", line 102, in addfile
    return addfile_grads(fname, False)
  File "F:\MeteoInfo\pylib\mipylib\dataset\midata.py", line 152, in addfile_grads
    meteodata.openGrADSData(fname)
        at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.readDataInfo(GrADSDataInfo.java:841)
        at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.readDataInfo(GrADSDataInfo.java:295)
        at org.meteoinfo.data.meteodata.MeteoDataInfo.openGrADSData(MeteoDataInfo.java:489)
        at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
        at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
        at java.lang.reflect.Method.invoke(Unknown Source)
java.lang.ArrayIndexOutOfBoundsException: java.lang.ArrayIndexOutOfBoundsException: 1

MeteoInfo

wuyu 发表于 2017-8-25 15:48
王老师，我运行了如下：
#Open station file as table data
stdata = readtable('F:/MeteoInfo/sample/Lo ...

看起来是你的GrADS数据文件在读取时出错。如果数据不是很大，你可以把ctl和数据文件附上，我有空了看看。
也可以加入MeteoInfo QQ群交流。

wuyu

taifenggrd.ctl如下：
dset ^taifeng.grd
undef 9.999E+20
ydef 277 linear 0.800000 0.18
xdef 441 linear 81.600000 0.180000
tdef 1 linear 00Z07aug2017 1mo
zdef 1 levels 1000
vars 1

UGRD10m   0  10mu

ENDVARS

stone

马上实验一下