风玫瑰图(Wind rose diagram)是一种特殊的极坐标堆叠图/统计直方图，其能够直观地表示某个地区一段时期内风向、风速的发生频率。
风玫瑰图在建筑规划、环保、风力发电、消防、石油站设计、海洋气候分析等领域都有重要作用，所以在一些顶级期刊中也能够看到它的身影，比如《Nature》：
由于Matlab中未收录风玫瑰图的绘制函数，因此需要大家自行设法解决。
本文使用WindRose工具(Daniel Pereira. Matlab Central, 2023)进行风玫瑰图的绘制，先来看一下成品效果：
特别提示：本期内容『数据+代码』已上传资源群中，加群的朋友请自行下载。有需要的朋友可以在公众号后台回复关键词【绘图桶】查看加群方式。
1. 数据准备
此部分主要是读取原始数据，并初始化绘图参数。
% 读取数据
load data.mat
D = direction;
S = speed;
2. 颜色定义
作图不配色就好比做菜不放盐，总让人感觉少些味道。
但颜色搭配比较考验个人审美，需要多加尝试。
这里直接使用TheColor配色工具中的SCI专业配色库：
%% 颜色定义
map = TheColor('sci',2068);
% map = flipud(map);
(点击上图查看TheColor具体功能)
获取方式：公众号（阿昆的科研日常）后台回复 TC
3. 风玫瑰图绘制
利用‘WindRose’工具，绘制风玫瑰图。
Options = {'anglenorth',0,...
'angleeast',90,...
'labels',{'N (0°)','NE (45°)','E (90°)','SE (135°)','S (180°)','SW (225°)','W (270°)','NW (315°)'},...
'freqlabelangle','auto',...
'nspeeds',7,...
'ndirections',18,...
'lablegend','Wind speeds in m/s',...
'legendtype',2,...
'legendvariable','W_s',...
'min_radius',0.1,...
'titlestring','Wind Rose Plot',...
'height',12*50,...
'width',15*50,...
'cmap',map};
figureHandle = WindRose(D,S,Options);
其中，‘Options’为自定义绘图参数集，包括标签、标题、图例、配色、组数等等，可以根据自己的需要进行设置。
绘制完成后，按所需格式、分辨率导出图像：
%% 图片输出
print('test.png','-r300','-dpng');
也可以尝试其它配色：
以上。
转自：“阿昆的科研日常”微信公众号
如有侵权，请联系本站删除！
