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

MeteoInfo

可以利用脚本程序（MeteoInfoLab环境下）来批量处理格点数据插值为站点数据。示例脚本程序如下，已经注释得很清楚了，这里就不多说了。

1. #Open station file as table data
   
2. stdata = readtable('D:/Temp/China_Prec_2010101420.csv', delimiter=',', format='%s%f%f')
   
3. x = stdata['Longitude']
   
4. y = stdata['Latitude']
   
5. #Add data file
   
6. f = addfile('D:/Temp/GrADS/model.ctl')
   
7. #Get data variable
   
8. psv = f['PS']
   
9. #Get time dimension length
   
10. tn = psv.dimlen(0)
    
11. #Loop
    
12. for i in range(0, tn):
    
13.     #Get dimension array
    
14.     ps = psv[i,'10:60','60:140']
    
15.     #Interpolate to stations
    
16.     ps_st = ps.tostation(x, y)
    
17.     #Add column to table data
    
18.     colname = 'PS' + str(i)
    
19.     stdata.addcoldata(colname, '%d', ps_st)
    
20. 
    
21. #Save table data to a file
    
22. fn = 'D:/Temp/test_st.csv'
    
23. stdata.savefile(fn)
    
24. print 'Finish...'

#Open station file as table data<br /> stdata = readtable('D:/Temp/China_Prec_2010101420.csv', delimiter=',', format='%s%f%f')<br /> x = stdata['Longitude']<br /> y = stdata['Latitude']<br /> #Add data file<br /> f = addfile('D:/Temp/GrADS/model.ctl')<br /> #Get data variable<br /> psv = f['PS']<br /> #Get time dimension length<br /> tn = psv.dimlen(0)<br /> #Loop<br /> for i in range(0, tn):<br />     #Get dimension array<br />     ps = psv[i,'10:60','60:140']<br />     #Interpolate to stations<br />     ps_st = ps.tostation(x, y)<br />     #Add column to table data<br />     colname = 'PS' + str(i)<br />     stdata.addcoldata(colname, '%d', ps_st)<br /> <br /> #Save table data to a file<br /> fn = 'D:/Temp/test_st.csv'<br /> stdata.savefile(fn)<br /> print 'Finish...'

sinohack

再次感谢 老师指教。但是老师，我运行的时候为什么老是报tn=psv.dimlen(0)这一行的错呢？AttributeError: DimVariable instance has no attribute 'dimlen'  另外老师在循环里写的ps=psv[i,(10,60),(60,140)]这是把经纬度选择在（10，60），（60，140）这个范围元素里面的经纬度的意思么？还有就是我试着把tn改成一个数字后再运行，程序就会报ps_st = ps.tostation(x, y) 这一行的错， instance has no attribute 'tostation'。请老师再给讲解一下，谢谢。

MeteoInfo

sinohack 发表于 2015-5-11 21:05
再次感谢 老师指教。但是老师，我运行的时候为什么老是报tn=psv.dimlen(0)这一行的错呢？AttributeError: D ...

1、请下载最新版本，下载地址见帖子1楼
2、(10, 60)和(60,140)是根据经纬度选择数据

sinohack

MeteoInfo 发表于 2015-5-11 22:10
1、请下载最新版本，下载地址见帖子1楼
2、(10, 60)和(60,140)是根据经纬度选择数据

嗯 好的  多谢老师 我这个终于做出来了   太感谢 老师了

yingkr

很有用~~~~~~~~~

liyun5918

学习中！
……

samvelly

{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}

wsjqxjy

,谢谢！

Sunny_H

psv = f['PS']   PS什么意思呢？

MeteoInfo

Sunny_H 发表于 2015-6-19 18:55
psv = f['PS']   PS什么意思呢？

变量名