rose函数绘制风玫瑰图示例

Aires

搜了下帖子，还没发现有发Matlab绘制风玫瑰图的程序。小A把一小段程序给大家，但其中还有许多地方需要修改~

数据用@lqouc给小A的青岛2012年1月日平均风向数据。如下：
54857.000    36.070   120.330   340.000     7.000
54857.000    36.070   120.330   220.000     3.000
54857.000    36.070   120.330   340.000     7.000
54857.000    36.070   120.330   340.000     6.000
54857.000    36.070   120.330   250.000     1.000
54857.000    36.070   120.330   340.000     3.000
54857.000    36.070   120.330   340.000     4.000
54857.000    36.070   120.330   340.000     2.000
54857.000    36.070   120.330   360.000     3.000
54857.000    36.070   120.330   340.000     4.000
54857.000    36.070   120.330   320.000     3.000
54857.000    36.070   120.330   200.000     5.000
54857.000    36.070   120.330   340.000     2.000
54857.000    36.070   120.330   340.000     2.000
54857.000    36.070   120.330   160.000     1.000
54857.000    36.070   120.330   340.000     2.000
54857.000    36.070   120.330   290.000     2.000
54857.000    36.070   120.330   290.000     2.000
54857.000    36.070   120.330   340.000     6.000
54857.000    36.070   120.330   340.000     8.000
54857.000    36.070   120.330   340.000    10.000
54857.000    36.070   120.330   340.000     7.000
54857.000    36.070   120.330   360.000     2.000
54857.000    36.070   120.330   290.000     7.000
54857.000    36.070   120.330   290.000     4.000
54857.000    36.070   120.330   180.000     3.000
54857.000    36.070   120.330   340.000     2.000
54857.000    36.070   120.330   110.000     4.000
54857.000    36.070   120.330   270.000     1.000
54857.000    36.070   120.330   290.000     3.000

画风玫瑰图可用Matlab里面自带的玫瑰图函数：rose。但是rose函数采用极坐标，即0度为x轴正半轴，逆时针旋转。而风向的定义是y轴正半轴为0度（北风为0度），顺时针旋转。这里用rose对风向进行统计后，用view函数作旋转。
程序如下：
% 用rose函数画风玫瑰图
% 数据为青岛2012年1月日平均风速、风向，第4列为风向
data=importdata('wind.txt');
dir=data(:,4);
dir=dir*pi/180;
x=0:22.5:337.5;
x=x*pi/180;
h=rose(dir,x);
set(h,'linewidth',2);
view(90,-90);
print -dpng -r300 wind_rose.png

图如下：

几点说明：
1.x是为了设置风向的16方位
2.rose中的参数都需要转换为弧度。rose其他用法查看help。
3.view函数通过方位角和仰角确定视角，具体查看help。

rose函数本身默认是这个样子的，想把外围的角度换为16方位的标注，找不到句柄在哪里。知道的希望留个言哈~

2556615693

我是菜鸟啊只能学习了

freekiller

@Aireslz能否将你的实例数据(wind.txt)贴上来。谢谢

Aires

freekiller 发表于 2013-11-11 16:52

                               
登录/注册后可看大图

@Aireslz能否将你的实例数据(wind.txt)贴上来。谢谢

已上传。

freekiller

Aires 发表于 2013-11-11 18:08

                               
登录/注册后可看大图

已上传。

已收到，多谢

张霖逸

默默的为楼主点个赞……

一个不小心

这个好，学习一下

tomming

很好，不过感觉这样画还不如直接用execl画好了。

爆浆炸鸡

正是需要的~~学习一下

bigocean

matlab 一直就有风玫瑰图的函数啊。
figure
load sunspot.dat  % Contains a 2-column vector named sunspot
rose(sunspot(:,2),12)