风玫瑰图(Wind rose diagram)是一种特殊的极坐标堆叠图/统计直方图，其能够直观地表示某个地区一段时期内风向、风速的发生频率。

风玫瑰图在建筑规划、环保、风力发电、消防、石油站设计、海洋气候分析等领域都有重要作用，所以在一些顶级期刊中也能够看到它的身影，比如《Nature》：

由于Matlab中未收录风玫瑰图的绘制函数，因此需要大家自行设法解决。

本文使用WindRose工具(Daniel Pereira. Matlab Central, 2023)进行风玫瑰图的绘制，先来看一下成品效果：

特别提示：本期内容『数据+代码』已上传资源群中，加群的朋友请自行下载。有需要的朋友可以关注同名公号 阿昆的科研日常，后台回复关键词【绘图桶】查看加入方式。

---

1. 数据准备

此部分主要是读取原始数据，并初始化绘图参数。

% 读取数据load data.matD = direction;S = speed;

2. 颜色定义

作图不配色就好比做菜不放盐，总让人感觉少些味道。

但颜色搭配比较考验个人审美，需要多加尝试。

这里直接使用TheColor配色工具中的SCI专业配色库：

%% 颜色定义map = TheColor('sci',2068);% map = flipud(map);

获取方式：公众号（阿昆的科研日常）后台回复 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');

也可以尝试其它配色：

以上。