MATLAB 彩色风羽图

Liujiaxin

在前人绘制风廓线雷达风羽图的基础上，通过for循环以及增加颜色索引
N_color=32;
mcp=colormap(parula(N_color));%生成颜色
P=sqrt(u^2+v^2);
[~,~,P_color] = histcounts(P,linspace(P_min,P_max,N_color));%颜色索引
mcp(P_color,:);
%check for scale and wind barb property specification
wbproperties = '''color'',mcp(P_color,:),''linewidth'',1.1';
实现彩色风羽图的绘制。

以下为一随机生成的风速的示例程序
clc;close all;clear all;
[x,y] = meshgrid(-20:4:20,-10:4:10);
px=rand(size(x,1),size(x,2)).*10;
py=rand(size(x,1),size(x,2)).*10;
scale=1.2;
[m,n]=size(x);
scale_x = scale*(max(max(x))-min(min(x)))/n;
scale_y = scale*(max(max(y))-min(min(y)))/m;
P=sqrt(px.^2+py.^2);
P_min=floor(min(min(P))/5)*5;
P_max=ceil(max(max(P))/5)*5;
figure;
for i=1:size(x,2)
    for j=1:size(x,1)

        windbarb3(x(j,i),y(j,i),px(j,i),py(j,i),scale_x,scale_y,P_min,P_max,scale);
        %原程序是四个风羽一起画，这里是一个一个的画，先不要增加水平及垂直方向的尺度scale_x scale_y，P_min,P_max是颜色索引的最小最大值

    end
end
colormap;
colorbar;
caxis([[P_min,P_max]]);%颜色条范围