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

Luo丹王道

王老师，我改成了
#Open station file as table data
stdata = readtable('E:/Temp/duan.csv', delimiter=',', format='%s%f%f')
x = stdata['Longitude']
y = stdata['Latitude']
#Add data file
f = addfile('E:/CN025/temp.nc')
#Get data variable
psv = f['temp']
#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 = 'temp' + str(i)
    stdata.addcol(colname, '%d', ps_st)
#Save table data to a file
fn = 'E:/Temp/duantem.csv'
stdata.savefile(fn)
print 'Finish...'
查了下您发的别的关于打开netCDF数据的帖子，还是不懂怎么读，请您帮帮忙~~~

MeteoInfo

Luo丹王道 发表于 2017-1-13 08:49
王老师，我改成了
#Open station file as table data
stdata = readtable('E:/Temp/duan.csv', delimiter ...

有错误信息吗？

Luo丹王道

run script...
Longitude
Latitude
SLF4J: Failed to load class "org.slf4j.impl.StaticLoggerBinder".
SLF4J: Defaulting to no-operation (NOP) logger implementation
SLF4J: See http://www.slf4j.org/codes.html#StaticLoggerBinder for further details.
temp

Luo丹王道

MeteoInfo 发表于 2017-1-13 08:58
有错误信息吗？

错误信息是
run script...
Longitude
Latitude
SLF4J: Failed to load class "org.slf4j.impl.StaticLoggerBinder".
SLF4J: Defaulting to no-operation (NOP) logger implementation
SLF4J: See http://www.slf4j.org/codes.html#StaticLoggerBinder for further details.
temp
Finish...
TypeError: __len__() should return an int

MeteoInfo

Luo丹王道 发表于 2017-1-13 09:21
错误信息是
run script...
Longitude

你先在MeteoInfoLab桌面环境中运行，调试好了再用Dos命令运行。

Luo丹王道

MeteoInfo 发表于 2017-1-13 10:19
你先在MeteoInfoLab桌面环境中运行，调试好了再用Dos命令运行。

meteoinfolab的桌面环境指的是MeteoInfo软件吗？
我用MeteoInfo导入netCDF数据，然后点“显示数据”，出来的数据都是错误的值，但是点绘图显示的又是正确的值范围，所以接下来的to station data也做不了~~~
我觉得我上面的DOS命令应该就是读入netCDF数据那几行有错误，可是就是不知道怎么写，所以卡在那了，请您帮帮我，谢谢！

MeteoInfo

Luo丹王道 发表于 2017-1-13 15:28
meteoinfolab的桌面环境指的是MeteoInfo软件吗？
我用MeteoInfo导入netCDF数据，然后点“显示数据”，出 ...

在MeteoInfo目录中运行MeteoInfoLab.exe文件，参考此贴：MeteoInfoLab脚本汇总贴
http://bbs.06climate.com/forum.p ... 151&fromuid=106
(出处: 气象家园)

vortexroc

请问如何把插值后得到的数组（比如示例中的ps_st变量）保留1位小数？

MeteoInfo

vortexroc 发表于 2017-2-22 16:59
请问如何把插值后得到的数组（比如示例中的ps_st变量）保留1位小数？

目前还没有小数位数的选项，以后会考虑加上。

dicemark

请问这个问题怎么解决，不是很懂，软件是在官网下的最新版，csv文件是上次你发的那个。请教一下
run script...
Traceback (most recent call last):
  File "<iostream>", line 6, in <module>
  File "C:\Users\ZQ\Desktop\MeteoInfo\pylib\mipylib\dataset\dimvariable.py", line 231, in dimlen
    return self.dims[idx].getLength()
AttributeError: DimVariable instance has no attribute 'dims'