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

Sunny_H

MeteoInfo 发表于 2015-6-19 19:53
变量名

好强大啊 谢谢！

Luo丹王道

按照您的脚本修改的脚本：
#Open station file as table data
stdata = readtable('E:/Temp/station.csv')
x = stdata['Longitude']
y = stdata['Latitude']
#Add data file
f = addfile('E:/CMIP5_Historical_pr_1901_2005_1x1.ctl')
#Get data variable
psv = f['pr']
#Get time dimension length
tn = psv.dimlen(0)
#Loop
for i in range(0, tn):
    #Get dimension array
    ps = psv[i,(10,60),(60,140)]
    #Interpolate to stations
    ps_st = ps.tostation(x, y)
    #Add column to table data
    colname = 'pr' + str(i)
    stdata.addcolumn(colname, '%d', ps_st)
#Save table data to a file
fn = 'E:/Temp/test_st.csv'
stdata.savefile(fn)
print 'Finish...'

错误是
run script...
File Format Error!
Longitude
Traceback (most recent call last):
  File "<iostream>", line 3, in <module>
  File "D:\work software\MeteoInfo_Java_1.2.7R1_Files\MeteoInfo\pylib\mipylib\midata.py", line 54, in __getitem__
    coldata = self.data.getColumnData(key)
        at org.meteoinfo.table.ColumnData.<init>(ColumnData.java:32)
        at org.meteoinfo.table.DataTable.getColumnData(DataTable.java:504)
        at org.meteoinfo.table.DataTable.getColumnData(DataTable.java:485)
        at org.meteoinfo.data.TableData.getColumnData(TableData.java:312)
        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.NullPointerException: java.lang.NullPointerException
王老师，这要怎么改呢？

MeteoInfo

Luo丹王道 发表于 2015-7-10 16:26
按照您的脚本修改的脚本：
#Open station file as table data
stdata = readtable('E:/Temp/station.csv' ...

你按照类似这样的格式准备站点数据文件：

Stid,Longitude,Latitude
57562,111.3600006,29.57999992
57584,113.0800018,29.37999916
52633,98.41000366,38.81000137
57662,111.6800003,29.04999924
59096,114.4800034,24.37000084
59072,112.3799973,24.78000069
57297,114.0500031,32.13000107
58111,116.3799973,33.86000061
57774,112.1600037,27.45000076
50877,129.5800018,46.29999924

Luo丹王道

MeteoInfo 发表于 2015-7-10 16:58
你按照类似这样的格式准备站点数据文件：

Stid,Longitude,Latitude

王老师，我下载的站点数据文件是您以前帖子的China_Prec_2010101420.csv
可是还是这个错误，呜呜~~~

MeteoInfo

Luo丹王道 发表于 2015-7-10 17:57
王老师，我下载的站点数据文件是您以前帖子的China_Prec_2010101420.csv
可是还是这个错误，呜呜~~~

我最近对MeteoInfoLab的更新可能影响了这个脚本的运行，我抽空测试一下。

Luo丹王道

MeteoInfo 发表于 2015-7-10 18:07
我最近对MeteoInfoLab的更新可能影响了这个脚本的运行，我抽空测试一下。

好的，麻烦王老师了~请您尽快测试，我下周赶着交数据，呜呜~

MeteoInfo

Luo丹王道 发表于 2015-7-10 22:02
好的，麻烦王老师了~请您尽快测试，我下周赶着交数据，呜呜~

已经更新了软件，可以在这里下载：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435

脚本程序会在1楼更新。

Luo丹王道

MeteoInfo 发表于 2015-7-10 23:50
已经更新了软件，可以在这里下载：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435

...

王老师，终于搞定了~~~非常非常非常感谢您啊！！！好佩服您，您太牛啦！！！
O(∩_∩)O谢谢~

lovechang1314

如有ctl中有多个时次，怎么加循环呢？

MeteoInfo

lovechang1314 发表于 2015-7-19 19:47
如有ctl中有多个时次，怎么加循环呢？

脚本里有对时间的循环