MeteoInfoLab脚本示例：风玫瑰图

MeteoInfo

利用极坐标可以绘制风玫瑰图，不过二者的区别还是挺大的。极坐标的0度在东边，逆时针递增；风玫瑰图的0度在北边，顺时针递增。为了方便绘制风玫瑰图增加了windrose函数，输入数据是风向和风速以及其它一些设置。

示例脚本用了Origin网站的一个示例数据（http://www.originlab.com/doc/Tutorials/Windrose-Graph），数据文件： windrose.txt (379 Bytes, 下载次数: 314)

2017-3-27 23:57 上传

点击文件名下载附件

1. fn = r'D:\Temp\ascii\windrose.txt'
   
2. ncol = numasciicol(fn)
   
3. nrow = numasciirow(fn)
   
4. a = asciiread(fn,shape=(nrow,ncol))
   
5. ws=a[:,0]
   
6. wd=a[:,1]
   
7. 
   
8. n = 16
   
9. wsbins = arange(0., 21.1, 4)
   
10. cols = makecolors(['k','y','r','b','g'], alpha=0.7)
    
11. rtickloc = [0.05,0.1,0.15,0.18]
    
12. ax, bars = windrose(wd, ws, n, wsbins, rmax=0.18, colors=cols, rtickloc=rtickloc)
    
13. colorbar(bars, shrink=0.6, vmintick=True, vmaxtick=True, xshift=10, \
    
14.     label='m/s', labelloc='bottom')
    
15. title('Windrose example')

fn = r'D:\Temp\ascii\windrose.txt'<br /> ncol = numasciicol(fn)<br /> nrow = numasciirow(fn)<br /> a = asciiread(fn,shape=(nrow,ncol))<br /> ws=a[:,0]<br /> wd=a[:,1]<br /> <br /> n = 16<br /> wsbins = arange(0., 21.1, 4)<br /> cols = makecolors(['k','y','r','b','g'], alpha=0.7)<br /> rtickloc = [0.05,0.1,0.15,0.18]<br /> ax, bars = windrose(wd, ws, n, wsbins, rmax=0.18, colors=cols, rtickloc=rtickloc)<br /> colorbar(bars, shrink=0.6, vmintick=True, vmaxtick=True, xshift=10, \<br />     label='m/s', labelloc='bottom')<br /> title('Windrose example')

︶ㄣ安定■丶

点完赞再睡～～

qdsonglinyy

每天都有新内容，赞啊，我是铁粉

自在飞花

谢谢楼主

lijian

谢谢楼主{:5_231:}

Amon8859

谢谢楼主

小1龙2虾3

要是能一起计算静风频率就更好了

静言_GRMC

可以把两个方位中间的空格去掉吗
两条数据紧挨着

MeteoInfo

静言_GRMC 发表于 2017-6-11 23:38
可以把两个方位中间的空格去掉吗
两条数据紧挨着

可以用width参数，缺省是0.5，width=1就没有空了。

wanghp

{:eb502:}{:eb502:}