- 积分
- 57612
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-21
- 最后登录
- 1970-1-1
|
发表于 2016-3-25 10:31:43
|
显示全部楼层
原版本对于micaps 11类数据的处理有些问题,进行了更新(需要下载最新的MeteoInfo Java 1.3.2R3版)。示例脚本如下:
 - #Plot map
- axesm(position=[0,0,1,1])
- lworld = shaperead('D:/Temp/map/bou2_4p.shp')
- geoshow(lworld, facecolor=[230,230,230])
- #Read micaps 11 data
- for i in range(0, 7, 3):
- fn = 'C:/Temp/16012020.%03i' % i
- print fn
- f = addfile_micaps(fn)
- u = f['U'][0,0,[20,25],[110,120]]
- v = f['V'][0,0,[20,25],[110,120]]
- if i > 0:
- cll() #Clear last add layer
- layer = barbsm(u[::2,::2], v[::2,::2], color='b')
- t = f.gettime(0)
- title(t.strftime('%Y-%m-%d %H:00'))
- outfn = 'C:/Temp/test/barb_' + t.strftime('%Y%m%d%H') + '_%i.png' % i
- savefig(outfn, 800, 400)
|
|