MeteoInfoLab脚本示例：站点数据风向杆图

MeteoInfo

绘制风向杆图的函数是barbs和barbsm（绘制在地图上），需要站点的经/纬度，风场U/V分量或者风速/风向。函数里的isuv参数默认是True，表明风场用的是U/V分量，如果isuv的值为False则表示风场用的是风向/风速数据。这里示例从MICAPS 1类数据中读出站点经/纬度和风速/风向数据并绘制站点风向杆图。

1. #Set data folders
   
2. basedir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
   
3. datadir = os.path.join(basedir, 'sample/MICAPS')
   
4. mapdir = os.path.join(basedir, 'map')
   
5. #Read shape files
   
6. lworld = shaperead(os.path.join(mapdir, 'country1.shp'))
   
7. lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
   
8. #Read station data
   
9. f = addfile_micaps(os.path.join(datadir, '10101414.000'))
   
10. lon = f['Longitude'][:]
    
11. lat = f['Latitude'][:]
    
12. ws = f['WindSpeed'][:]
    
13. wd = f['WindDirection'][:]
    
14. #Plot
    
15. axesm(bgcolor=(204,255,255))
    
16. geoshow(lworld, edgecolor='k', facecolor=(255,251,195))
    
17. geoshow(lchina, edgecolor='k')
    
18. layer = barbsm(lon, lat, wd, ws, color='b', isuv=False)
    
19. yticks([20,30,40,50])
    
20. title('Station barb plot example')
    
21. xlim(72, 136)
    
22. ylim(16, 55)

#Set data folders<br /> basedir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'<br /> datadir = os.path.join(basedir, 'sample/MICAPS')<br /> mapdir = os.path.join(basedir, 'map')<br /> #Read shape files<br /> lworld = shaperead(os.path.join(mapdir, 'country1.shp'))<br /> lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))<br /> #Read station data<br /> f = addfile_micaps(os.path.join(datadir, '10101414.000'))<br /> lon = f['Longitude'][:]<br /> lat = f['Latitude'][:]<br /> ws = f['WindSpeed'][:]<br /> wd = f['WindDirection'][:]<br /> #Plot<br /> axesm(bgcolor=(204,255,255))<br /> geoshow(lworld, edgecolor='k', facecolor=(255,251,195))<br /> geoshow(lchina, edgecolor='k')<br /> layer = barbsm(lon, lat, wd, ws, color='b', isuv=False)<br /> yticks([20,30,40,50])<br /> title('Station barb plot example')<br /> xlim(72, 136)<br /> ylim(16, 55)

风格而才

风格而才 发表于 2022-2-22 11:10
麻烦问一下王老师，EC再分析资料能画带风向杆的风场图么，有没有参考的例子呀，找了好久没找到。micaps资料 ...

不好意思  找到了！希望对其他人有帮助。speed = sqrt(u*u+v*v)
layer = barbsm(u, v, speed,color='k', size = 10)

Nancy_fq

王老师给力！难道是传说的沙发……

qxjy009

谢谢王老师，我现在还有两个问题请教，
一个是您这边的isuv,os等颜色跟其他字符不同，我的没有，不过也能出这个风向杆图，但是风向杆颜色跟您一样设置 color='b'，出来是还是黑色的，改成其他的也没变化
第二个就是我加载了一个地形图，它的色标我自己设置的lgs文件，然后脚本中是这样写的
ldx = shaperead('D:/javasoft/sample/high-hz.grd')
geoshow(ldx, facecolor='D:/javasoft/sample/Legend/out1.lgs')
但是出图没有按我设的色标，这个加载颜色的命令是不是有问题

MeteoInfo

qxjy009 发表于 2016-12-15 10:56
谢谢王老师，我现在还有两个问题请教，
一个是您这边的isuv,os等颜色跟其他字符不同，我的没有，不过也能 ...

建议使用最新版本。

用makelegend函数载入色标文件：
ls = makelegend('D:/javasoft/sample/Legend/out1.lgs')

在geoshow函数中用symbolspec参数：
geoshow(ldx, symbolspec=ls)

风格而才

麻烦问一下王老师，EC再分析资料能画带风向杆的风场图么，有没有参考的例子呀，找了好久没找到。micaps资料我知道能画来，但是EC资料只有uv风分量，我只会画带箭头的矢量图，带风向杆的该怎么画呢。