24小时热门版块排行榜

返回列表

bucheron

金虫 (小有名气)

应助: 21 (小学生)
金币: 652
散金: 16
红花: 5
帖子: 178
在线: 189.9小时
虫号: 2183449
注册: 2012-12-13
性别: GG
专业: 燃烧学

[交流] 怎样在同一个图中画两个等高线图？

目的：存在两组等高线，两组等高线有交集。将两个等高线画在同一个图中，每组等高线使用不同的colormap，两组等高线的交集应该是两组颜色的重叠。

实现方法一

使用contour，不用填充颜色的contourf，以便第二组等高线不会覆盖住第一组的交集部分。使用例如colormap([cool(64);gray(64)])命令，将colormap变成两种色彩集的合并。

例子1

clear all % clear all the data in the workspace
close all % close all the current figures
clc % clear the command window
%-----------------------------
% generate the data
[X,Y,Z1] = peaks(64); % X,Y,Z data of the first group
Z2=X.^2+Y.^2; % Z data of the second group
%%
% first method
% use two axes in the same figure
% set the figure's properties
h=figure
scrsz = get(0,'ScreenSize');
set(h,'Position',[scrsz(4).*(1/8) scrsz(4).*(1/20) scrsz(3)*3/5 scrsz(4).*(4/5)])
%--------------------------------------------------------------------------
% the first axes
hAxes(1)=axes('Unit','pixels','position',[100 100 400 400]);
%--------
% combine two colormap in the same one, here we use 'jet' and 'gray'
colormap([jet(64);gray(64)])
% the jet color will be used two index the first contour and the gray will
% be used for the second one.
%--------
[c1,h1] = contour(X(1,

,Y(:,1),Z1); % plot the first contour
set(h1,'ShowText','on')
cmin1=min(Z1(

);
cmax1=max(Z1(

);
% use caxis to set the limit of the colorbar, the first one will occupy the
% low half one
caxis(hAxes(1),[cmin1,2*cmax1-cmin1])
% then define the axis
set(hAxes(1),'YColor','k','Box','on');
set(hAxes(1),'FontName','Helvetica','FontSize',20,'LineWidth',1)
xlabel(hAxes(1),'$ X $ [-]','Color','k','Interpreter','LaTex','FontSize',20);
ylabel(hAxes(1),'$ Y $ [-]','Color','k','Interpreter','LaTex','FontSize',20);
set(hAxes(1),'xlim',[min(X(1,

), max(X(1,

)],'YAxisLocation','left','Color','w');
set(hAxes(1),'ylim',[min(Y(:,1)), max(Y(:,1))]);
%--------------------------------------------------------------------------
% the second axes
hAxes(2)=axes('Unit','pixels','position',get(hAxes(1),'position'));
[c2,h2] = contour(X(1,

,Y(:,1),Z2); % plot the second contour
cmin2=min(Z2(

);
cmax2=max(Z2(

);
% use caxis to set the limit of the colorbar, the second one will occupy the
% higher half one
caxis(hAxes(2),[2*cmin2-cmax2,cmax2])
set(hAxes(2),'visible', 'off') % the the visable of axes 2 off

实现方法二
原则上来说，在一个figure中只能使用一种colormap，所以如果不像第一种方法定义colorbar，两组等高线的色调是一样的。解决这个问题可以采用freezecolor方法，即第一个contour采用一中colormap，然后freeze，对第二个contour采用另外一个colormap，也freeze。
具体可以参考
http://ycc.math.fju.edu.tw/research/biology/b2012fMRI/fMRICODE/simtb_v18/simtb_v18/display/freezeColors/demo/html/freezeColors_pub.html

回复此楼