基于MATLAB的风速风向风玫瑰图绘制方法与应用研究

基于matlab根据风速风向画风玫瑰图。

最近在分析气象数据的时候，总得处理风速风向这种CP组合。今儿给大伙儿整点实用的——用MATLAB画风玫瑰图。这玩意儿本质上就是个极坐标版的直方图，但配上彩色扇区确实比表格数据直观多了。

基于matlab根据风速风向画风玫瑰图。

先搞点模拟数据热热身。假设我们有个月的观测记录，每天24组数据：

direction = rand(720,1)*360; % 随机生成720个0-360度的风向 speed = 3 + randn(720,1)*2; % 风速均值3m/s的正态分布 speed(speed<0) = 0; % 风速不能为负

接下来是核心操作——数据分箱。风向要分成16个扇区（每个22.5度），风速分成五个等级：

theta_edges = linspace(0, 360, 17); % 16个风向区间边界 speed_edges = [0 2 5 8 12 20]; % 风速分级标准 % 双重分箱操作 [~, dir_bins] = histc(direction, theta_edges); [~, spd_bins] = histc(speed, speed_edges);

这里有个骚操作，用accumarray统计联合分布：

count_matrix = accumarray([dir_bins, spd_bins], 1, [16,5]); count_matrix(:,end) = []; % 剔除超出风速范围的数据 total = sum(count_matrix,2); frequency = count_matrix ./ total; % 各风速等级占比

重点来了，极坐标可视化部分。注意MATLAB原生的polarplot有点难用，咱们换个姿势玩：

figure('Color','w') ax = polaraxes(); theta = deg2rad(linspace(0, 360, 17)); % 转换成弧度 % 堆叠式柱状图 for k = 1:size(frequency,2) h = polarhistogram('BinEdges', theta, 'BinCounts', frequency(:,k)*100,... 'FaceColor', jet(5)(k,:), 'EdgeColor','none'); hold on end % 调教图形细节 ax.ThetaDir = 'clockwise'; % 风向顺时针增加 ax.ThetaZeroLocation = 'top'; % 0度指向正北 rlim([0 40]) % 半径范围 ax.RTickLabel = {'10%','20%','30%','40%'}; legend({'0-2m/s','2-5m/s','5-8m/s','8-12m/s','>12m/s'},'Location','southoutside')

运行完你会看到个五彩斑斓的玫瑰，不过可能遇到两个坑：一是颜色映射不匹配，这时候检查jet(5)是不是生成了五个颜色；二是图例位置可能跑偏，把Location参数改成'eastoutside'试试。

进阶玩法可以加个风速统计表：

uistack(ax,'bottom') % 把极坐标轴置底 annotation('textbox',[0.7 0.6 0.2 0.3],'String',... sprintf('最大风速: %.1fm/s\n平均风速: %.1fm/s',max(speed),mean(speed)),... 'EdgeColor','none','FontSize',10)

实测发现用histogram2也能快速出图，但自定义程度低：

histogram2(direction, speed, theta_edges, speed_edges,... 'DisplayStyle','tile','ShowEmptyBins','on') set(gca,'ThetaDir','clockwise','ThetaZeroLocation','top')

最后友情提示：处理真实数据时记得清洗异常值。曾经有个项目因为传感器故障导致出现3600m/s的风速，画出来的玫瑰图直接变成仙人掌，血泪教训啊...（完）