MeteoInfoLab脚本：TXT站点数据绘制等值线填色图

静言_GRMC

辞觉 发表于 2017-12-21 10:22
请教一下按照您的脚本运行的，但是我的图为什么站点位置不对呢？有的都跑到省外去了

确定经纬度没问题的话，那就可能地图不太对

afan1996

感谢分享(*^▽^*)~

yesyesnotion

请问“gtemp,gx,gy = griddata((lon, lat), pm, xi=(x, y), method='idw', radius=0.8)”中gx和gy是有用的变量吗？后边也没有出现这两个量

yesyesnotion

请问下这个错误怎么解决？纠结了一天了，没有找到答案

falconet

很好！！学习了！

lifetime

levs = [0.1,10,25,50,100,200]
cols = [(165,243,141),(61,185,63),(99,184,249),(0,0,254),(243,5,238),(129,0,64)]
layer = contourfm(x, y, gtemp,levs,colors=cols)
slayer = scatterm(lon, lat, rain, levs, colors=cols, size=8)
运行的时候，提示：
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "C:\Program Files (x86)\MeteoInfo\pylib\mipylib\plotlib\miplot.py", line 2121, in scatterm
    r = gca.scatter(*args, **kwargs)
  File "C:\Program Files (x86)\MeteoInfo\pylib\mipylib\plotlib\mapaxes.py", line 776, in scatter
    ls = plotutil.setlegendscheme_point(ls, **kwargs)
  File "C:\Program Files (x86)\MeteoInfo\pylib\mipylib\plotlib\plotutil.py", line 536, in setlegendscheme_point
    kwargs['color'] = colors
IndexError: index out of range: 6

请问是什么原因？

MeteoInfo

lifetime 发表于 2018-10-26 19:38
levs = [0.1,10,25,50,100,200]
cols = [(165,243,141),(61,185,63),(99,184,249),(0,0,254),(243,5,238), ...

cols的长度应该比levs多1，也就是cols要有7个元素。

ymy_00000

fname='H:/ymy/Precipitation6h/2015/6/SURF_WEA_PRE_6HOUR-2015060100.txt'
ncol = numasciicol(fname)
nrow = numasciirow(fname)

#Read station data
a = asciiread(fname,shape=(nrow,ncol))
lon = a[:,6]
lat = a[:,5]
rain= a[:,8]
#griddata function - interpolate
x = arange(75, 135, 0.28169015)
y = arange(15, 55, 0.28169015)
prg = griddata((lon, lat), rain, xi=(x, y), method='idw', radius=0.5)[0]
prg.savedata('H:/ymy/6小时降水/IDW/2015/6/rain_test.txt')

运行的时候提示：
>> run script...
Traceback (most recent call last):
  File "<iostream>", line 3, in <module>
  File "E:\ymy\MetoeInfo\MeteoInfo_1.6.2\MeteoInfo\pylib\mipylib\dataset\midata.py", line 440, in numasciicol
    ncol = ArrayUtil.numASCIICol(filename, delimiter, headerlines)
        at java.io.FileInputStream.open0(Native Method)

        at java.io.FileInputStream.open(Unknown Source)

        at java.io.FileInputStream.<init>(Unknown Source)

        at java.io.FileInputStream.<init>(Unknown Source)

        at org.meteoinfo.data.ArrayUtil.numASCIICol(ArrayUtil.java:171)

        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.io.FileNotFoundException: java.io.FileNotFoundException: H:\ymy\Precipitation6h\2015\6\SURF_WEA_PRE_6HOUR-2015060100.txt (???????????)

我想请问下这个是什么原因导致的！

MeteoInfo

ymy_00000 发表于 2018-12-17 14:01
fname='H:/ymy/Precipitation6h/2015/6/SURF_WEA_PRE_6HOUR-2015060100.txt'
ncol = numasciicol(fname) ...

ava.io.FileNotFoundException  这个应该不难理解吧

ymy_00000

MeteoInfo 发表于 2018-12-17 14:31
ava.io.FileNotFoundException  这个应该不难理解吧

可我的文件名和文件路径是正确的，所以我不太清楚这个错误是什么原因导致的