MeteoInfoLab脚本示例：站点数据绘制等值线

气象万千2011

MeteoInfo 发表于 2015-7-11 18:12
MeteoInfoLab是用Jython对MeteoInfo库的功能封装了很多类似NCL和MatLab函数，编写脚本更加简洁、方便。可 ...

嗯，谢谢王老师，这个软件我还需要多多学习

ldepn

运行的时候出现这样一个错误
run script...
Traceback (most recent call last):
  File "<iostream>", line 16, in <module>
AttributeError: 'NoneType' object has no attribute 'griddata'

MeteoInfo

ldepn 发表于 2015-7-17 12:36
运行的时候出现这样一个错误
run script...
Traceback (most recent call last):

估计是你的站点数据没获取正确（是空变量None Type），当然就无法格点化了（griddata）。最好贴出你的脚本程序和数据文件以便查找问题。

ldepn

MeteoInfo 发表于 2015-7-17 12:57
估计是你的站点数据没获取正确（是空变量None Type），当然就无法格点化了（griddata）。最好贴出你的脚 ...

#Set data folders
basedir = 'E:/Baidudownload/MeteoInfo/MeteoInfo(Java)'
datadir = os.path.join(basedir, 'sample/MICAPS')
mapdir = os.path.join(basedir, 'map')
#Read shape files
bou2_layer = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
bou1_layer = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))
china_layer = shaperead(os.path.join(mapdir, 'china.shp'))
city_layer = shaperead(os.path.join(mapdir, 'res1_4m.shp'))
#Read station
dataf = addfile_micaps(os.path.join(datadir, '10101414.000'))
pr = f.stationdata('Precipitation6h')
#griddata function - interpolate
x = arange(75, 135, 0.5)
y = arange(18, 55, 0.5)
prg = pr.griddata((x, y), method='idw', radius=3)
#Plot
axesm()
geoshow(bou2_layer, edgecolor='lightgray')
geoshow(bou1_layer, facecolor=(0,0,255))
geoshow(city_layer, facecolor='r', size=4, labelfield='NAME', fontname=u'楷体', fontsize=16, yoffset=15)
geoshow(china_layer, visible=False)
levs = [0.1, 1, 2, 5, 10, 20, 25, 50, 100]
cols = [(255,255,255),(170,240,255),(120,230,240),(200,220,50),(240,220,20),(255,120,10),(255,90,10), \
(240,40,0),(180,10,0),(120,10,0)]
layer = contourfm(prg, levs, colors=cols)
masklayer(china_layer, [layer])
colorbar(layer)
xlim(72, 136)
ylim(16, 55)
text(95, 52, u'全国降水量实况图', fontname=u'黑体', fontsize=16)
text(95, 50, u'(2010-10-14 08:00 至 2010-10-14 14:00)', fontname=u'黑体', fontsize=14)
#Add south China Sea
sc_layer = bou1_layer.clone()
axesm(position=[0.14,0.18,0.15,0.2], axison=False)
geoshow(sc_layer, facecolor=(0,0,255))
xlim(106, 123)
ylim(2, 23)

MeteoInfo

ldepn 发表于 2015-7-17 14:02
#Set data folders
basedir = 'E:/Baidudownload/MeteoInfo/MeteoInfo(Java)'
datadir = os.path.join( ...

#Read station
dataf = addfile_micaps(os.path.join(datadir, '10101414.000'))
pr = f.stationdata('Precipitation6h')

应该是：

#Read station data
f = addfile_micaps(os.path.join(datadir, '10101414.000'))
pr = f.stationdata('Precipitation6h')

ldepn

MeteoInfo 发表于 2015-7-17 14:09
#Read station
dataf = addfile_micaps(os.path.join(datadir, '10101414.000'))
pr = f.stationdata( ...

噢噢，我粗心了。

风荷1

请问王老师，脚本中把图片输出到文件的语句是什么呢?

MeteoInfo

风荷1 发表于 2015-7-31 09:45
请问王老师，脚本中把图片输出到文件的语句是什么呢?

参考此帖：发布MeteoInfo Java 1.2.2
http://bbs.06climate.com/forum.p ... 995&fromuid=106
(出处: 气象家园)

风荷1

谢谢答复。。。。。

iij

老师：您好，我用cressman插值，XNum、YNum、RadList设置有没有什么规定还是随意设置的呢？我有99个站点数据，XNum、YNum设置的15，RadList设置的10,8,6,4，得到的结果总感觉怪怪的，不是很正确。
还想请问一下老师，C#调用出图，color自定义怎么定义呢？谢谢老师。D:\5.PNG