2014 网站建设/手机百度搜索引擎入口

在之前的文章中，分享了Matlab箭头图的绘制模板：

进一步，如果我们想对每一个箭头赋上颜色，以更加直观地表示其模值的大小，该怎么操作呢？

那么，来看一下模值赋色的箭头图的绘制模板。

先来看一下成品效果：

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

---

模板中最关键的部分内容：

1. 数据准备

此部分主要是读取原始数据。

% 读取数据load data.mat

2. 颜色定义

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

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

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

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

3. 模值赋色的箭头图绘制

使用‘quiver’命令，绘制初始箭头图。

q = quiver(X,Y,U,V);hTitle = title('Visualization of air flow data over North America');hXLabel = xlabel('xaxis');hYLabel = ylabel('yaxis');

4. 细节优化

利用ColortheArrow函数(Suever)，将初始箭头图赋上之前选择的颜色并对线型进行调整：

% 赋色ColortheArrow(q,map)colorbar% 线型调整q.LineWidth = 1;q.ShowArrowHead = 'on';

然后，对坐标轴细节等进行美化：

% 坐标区调整axis tightset(gca, 'Box', 'off', ...                                   % 边框         'LineWidth', 1,...                                  % 线宽         'XGrid', 'off', 'YGrid', 'off', ...                 % 网格         'TickDir', 'out', 'TickLength', [.01 .01], ...      % 刻度         'XMinorTick', 'off', 'YMinorTick', 'off', ...       % 小刻度         'XColor', [.1 .1 .1],  'YColor', [.1 .1 .1])        % 坐标轴颜色% 字体和字号set(gca, 'FontName', 'Arial', 'FontSize', 10)set([hXLabel, hYLabel], 'FontSize', 11, 'FontName', 'Arial')set(hTitle, 'FontSize', 12, 'FontWeight' , 'bold')% 背景颜色set(gcf,'Color',[1 1 1])% 添加上、右框线hold onXL = get(gca,'xlim'); XR = XL(2);YL = get(gca,'ylim'); YT = YL(2);xc = get(gca,'XColor');yc = get(gca,'YColor');plot(XL,YT*ones(size(XL)),'color', xc,'linewidth',1)plot(XR*ones(size(YL)),YL,'color', yc,'linewidth',1)

设置完毕后，以期刊所需分辨率、格式输出图片。

%% 图片输出figW = figureWidth;figH = figureHeight;set(figureHandle,'PaperUnits',figureUnits);set(figureHandle,'PaperPosition',[0 0 figW figH]);fileout = 'test';print(figureHandle,[fileout,'.png'],'-r300','-dpng');

ColortheArrow函数：

function ColortheArrow(q,map)​% author Suever​%// Compute the magnitude of the vectorsmags = sqrt(sum(cat(2, q.UData(:), q.VData(:),reshape(q.WData, numel(q.UData), [])).^2, 2));​%// Get the current colormapcurrentColormap = colormap(map);​%// Now determine the color to make each arrow using a colormap[~, ~, ind] = histcounts(mags, size(currentColormap, 1));​% clims = num2cell(get(gca, 'clim'));% [~, ~, ind] = histcounts(mags, linspace(clims{:}, size(currentColormap, 1)));​%// Now map this to a colormap to get RGBcmap = uint8(ind2rgb(ind(:), currentColormap) * 255);cmap(:,:,4) = 255;cmap = permute(repmat(cmap, [1 3 1]), [2 1 3]);​%// We repeat each color 3 times (using 1:3 below) because each arrow has 3 verticesset(q.Head,'ColorBinding', 'interpolated','ColorData', reshape(cmap(1:3,:,:), [], 4).'); %'​%// We repeat each color 2 times (using 1:2 below) because each tail has 2 verticesset(q.Tail,'ColorBinding', 'interpolated','ColorData', reshape(cmap(1:2,:,:), [], 4).');​end

以上。