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

sfz13851581045 · 发表于 2020-10-16 21:07:37

王老师，您好，我有个疑问想请教您一下！请问通过双线性方法插值成站点数据的时候，是插值最靠近该站点四周的格点数据得到的嘛？

MeteoInfo · 发表于 2020-10-16 22:46:05

sfz13851581045 发表于 2020-10-16 21:07
王老师，您好，我有个疑问想请教您一下！请问通过双线性方法插值成站点数据的时候，是插值最靠近该站点四周 ...

是的

sfz13851581045 · 发表于 2020-10-16 23:10:07

感谢王老师的解惑

sfz13851581045 · 发表于 2020-10-30 18:11:48

脚本太棒了，很实用，请问王老师，在保存 table data 的时候可以设置数据的保存格式嘛？

MeteoInfo · 发表于 2020-10-31 10:45:55

sfz13851581045 发表于 2020-10-30 18:11
脚本太棒了，很实用，请问王老师，在保存 table data 的时候可以设置数据的保存格式嘛？

savefile函数有 date_format 和 float_format 参数可以设置数据格式。

upuil · 发表于 2021-2-26 16:43:38

本帖最后由 upuil 于 2021-2-26 16:57 编辑
王老师，请问下面的报错该怎么解决呢?

>>> run script...
Traceback (most recent call last):
File "D:\interplation.py", line 17, in <module>
ps_st = ps.tostation(x,y)
File "D:\MeteoInfo_2.4.5\MeteoInfo\pylib\mipylib\numeric\core\dimarray.py", line 1160, in tostation
r = gdata.data.toStation(x.aslist(), y.aslist())
at org.meteoinfo.data.GridData.toStation(GridData.java:639)
at java.base/jdk.internal.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at java.base/jdk.internal.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:64)
at java.base/jdk.internal.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43)
at java.base/java.lang.reflect.Method.invoke(Method.java:564)
java.lang.ClassCastException: java.lang.ClassCastException: class java.lang.String cannot be cast to class java.lang.Number (java.lang.String and java.lang.Number are in module java.base of loader 'bootstrap')
>>>

复制代码

代码是这样的：

#Open station file as table data
stdata = readtable('D:/stations.txt', delimiter=' ')
x = stdata['lon']
y = stdata['lat']
#Add data file
f = addfile('D:/2019070103.nc')
#Get data variable
psv = f['T2']
#Get time dimension length
tn = psv.dimlen(0)
#Loop
for i in range(0, tn):
#Get dimension array
ps = psv[i,:,:]
#Interpolate to stations
#ps_st = ps.tostation(x, y)
ps_st = ps.tostation(x,y)
#Add column to table data
colname = 'T2' + str(i)
stdata.addcol(colname, '%d', ps_st)
#Save table data to a file
fn = 'D:/test_st.csv'
stdata.savefile(fn)
print 'Finish...'

复制代码

MeteoInfo · 发表于 2021-2-26 17:11:50

upuil 发表于 2021-2-26 16:43
本帖最后由 upuil 于 2021-2-26 16:57 编辑
王老师，请问下面的报错该怎么解决呢?

readtable函数加上format参数，比如format='%2f'指数据中前两列是float类型，具体如何加format要根据数据文件的情况。不加format确实读出来的数组都是字符串类型。

stdata = readtable('D:/station.txt', format='%2f')

upuil · 发表于 2021-2-26 17:20:36

MeteoInfo 发表于 2021-2-26 17:11
readtable函数加上format参数，比如format='%2f'指数据中前两列是float类型，具体如何加format要根据数据 ...

王老师刚那个问题解决了，谢谢您，又出现一个新问题，麻烦您再看一下！

>>> run script...
Traceback (most recent call last):
File "D:\interplation.py", line 22, in <module>
stdata.addcol(colname, '%d', ps_st)
File "D:\MeteoInfo_2.4.5\MeteoInfo\pylib\mipylib\numeric\core\mitable.py", line 297, in addcol
self.data.addColumn(index, colname, dtype)
TypeError: addColumn(): 1st arg can't be coerced to int

复制代码

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

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