How to plot two contour filled plots of two different datasets using different colormap and caxis range in the same axis?

Question

ANKAN SARKAR 2023-11-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2049527-how-to-plot-two-contour-filled-plots-of-two-different-datasets-using-different-colormap-and-caxis-ra

评论： ANKAN SARKAR 2023-11-28

Dear all,

I want to plot two contour filled plots of two different variables with different colormaps and caxis ranges in the same axis. I tried with 'hold on' to plot two contour filled plots in the same axis, but it didn't work. Maybe it doesn't work with two contour filled plots?? I have gone through a previous conversation here (https://www.matlab.com/matlabcentral/answers/827960-setting-color-scale-of-two-datasets-in-the-same-axis-in-matlab), and modfied the code as per my case. But, I had no luck. I want the kind of figure like the figure shown in Ding et al. (2020).

I am attaching this figure as a sample for your convenience.

Can anyone please help me to solve the issue? That will be really helpful for me. My modifed MATLAB code is given below:

figure(1)
% one axes
contourf(X,Y,A,100,'linecolor','none')
cb1 = colorbar;
caxis([-2 2])
xlim([20 26])
ylim([100 1000])
h1 = gca;
h2 = axes;
% the other axes
contourf(X,Y,B,100,'linecolor','none')
% set(h2,'Color','none')
caxis([-0.25 0.25])
xlim([20 26])
ylim([100 1000])
cb2 = colorbar;
box on
set(h1,'yscale','log','YDir','reverse','Position', [0.45,0.48,0.26,0.42],'YTick','','YTickLabel','','XTick','','XTickLabel','',...
    'FontName', 'Arial','FontSize',18,'FontWeight','bold','LineWidth', 1.5)
set(h2,'yscale','log','YDir','reverse','Position', [0.45,0.48,0.26,0.42],'YTick','','YTickLabel','','XTick','','XTickLabel','',...
    'FontName', 'Arial','FontSize',18,'FontWeight','bold','LineWidth', 1.5)
% store the position and limits because adjusting the colorbar will mess these up
axpos = get(h1,'position');
xl = get(h1,'xlim');
yl = get(h1,'ylim');
set(cb1,'Position',[0.718,0.483,0.007,0.199],'YTick',-2:0.5:2,'FontName', 'Arial','FontSize',10.5,'FontWeight','bold','LineWidth', 1.5)
set(cb2,'Position',[0.718,0.698,0.007,0.199],'YTick',-0.25:0.1:0.25,'FontName', 'Arial','FontSize',10.5,'FontWeight','bold','LineWidth', 1.5)
% reassert axes geometry so they match
set(h1,'position',axpos,'xlim',xl,'ylim',yl);
set(h2,'position',axpos,'color','none','xlim',xl,'ylim',yl);
% set one colormap for each axes
cmap1 =jet(40);  % Specify the desired number of colors
colormap(h1,cmap1)
cmap2 =parula(60);  % Specify the desired number of colors
colormap(h2,cmap2)

4 个评论
显示 2更早的评论隐藏 2更早的评论

ANKAN SARKAR 2023-11-20

Hello @Dyuman Joshi,

I meant that when I added hold on in the code for overlaying two contourf plot on the same axis, it didn't happen. Earlier, I used 'hold on' successfully when I wanted to plot a contour filled plot for one variable and contours for another variable. But, 'hold on' is not working when I am trying to plot two contour filled plot on the same axis.

And, I am using WRF output (netcdf file). Since dataset is large and very diifficult to share, I am giving a brief description of my dataset and what I want. I want to plot WRF variable as latitude (x-axis) vs pressure (y-axis) similar to the attached figure. X is latitude , Y is Pressure, and A or B is the WRF variable. All three have same size (34*415). And If I run each contourf plot individually, it is coming. But, I want both of the contour filled plot to be plotted on the same axis. But, with the given code, maybe one contourf plot is getting overlapped with each other. Maybe, that's why it's not coming. So, is there anything wrong with my code? Can you please tell me?That will be very helpful.

Thannks and regards,

Ankan

DGM 2023-11-20

I'm a bit distracted at the moment, but as I'm sure you've found out by now, you can't do this in a single parent axes (or at least not reasonably). As in the linked thread, you need to use two overlaid axes and you'd then also need to manipulate the transparency of the upper contour() object accordingly.

Manipulating the contour transparency is a bit of an obstacle, depending on the version.

https://www.mathworks.com/matlabcentral/answers/1864708-add-transparency-to-a-contourf#answer_1113543

Contour objects had a scalar FaceAlpha property added in R2022b, so that doesn't help you either way. The first two options may be useful.

There may be other simpler options available via strategic selection of the LevelList property of the contour object. Having relevant example data would help.

ANKAN SARKAR 2023-11-20

Hello @DGM,

I understood that two seperate axes should be created for the two ontourf and then these axes should be overlaid. That's why I used the following lines in my MATLAB code:

h1 = gca;

h2 = axes;

I also tried linkaxes and linkprop to link those axes as sugggested in the following two links:

https://in.mathworks.com/matlabcentral/answers/101346-how-do-i-use-multiple-colormaps-in-a-single-figure-in-r2014a-and-earlier

https://in.mathworks.com/matlabcentral/answers/194554-how-can-i-use-and-display-two-different-colormaps-on-the-same-figure

But I had no luck. Any help will be greatly appreciated.

Thanks and regards,

Ankan

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

DGM 2023-11-21

编辑：DGM 2023-11-21

在 MATLAB Online 中打开

b.png
a.png

I just threw together two different examples using the ways I mentioned. Neither is ideal. Manipulating contour() objects is just problematic. Some part of my brain is telling me that there's something I'm missing, but that's probably just the sleep deprivation talking.

This example demonstrates how truncating the LevelList can be used to open up the periphery of a contour. This only works if the "background" of the upper contour at its periphery corresponds to its minimal values (like a 2D gaussian or something)

% make some fake Z data from some images

A = imread('a.png');

B = imread('b.png');

A = im2double(A)*50 - 25;

B = im2double(B)*100;

% the background axes

contourf(A,'edgecolor','none');

cb1 = colorbar;

h1 = gca;

h2 = axes;

% the foreground axes

[~,hc2] = contourf(B,'edgecolor','none');

set(h2,'color','none') % make the axes background itself invisible

cb2 = colorbar;

% make the base of the upper contour transparent by truncating LevelList

% this approach really only works if you're truncating the lower levels

% is there a workaround? i don't know.

hc2.LevelList(1) = [];

% store the position and limits because adjusting

% the colorbar will mess these up

axpos = get(h1,'position');

xl = get(h1,'xlim');

yl = get(h1,'ylim');

% adjust the colorbars so they don't overlap

barl = 0.48;

cb1.Position(2) = cb1.Position(2)+cb1.Position(4)*(1-barl);

cb1.Position(4) = cb1.Position(4)*barl;

cb2.Position(4) = cb2.Position(4)*barl;

% reassert axes geometry so they match

set(h1,'position',axpos,'xlim',xl,'ylim',yl);

set(h2,'position',axpos,'color','none','xlim',xl,'ylim',yl);

% set one colormap for each axes

colormap(h1,turbo)

colormap(h2,summer)

This uses manipulation of undocumented properties to control transparency. This is ugly, and will probably become problematic. Changes to the figure will cause the transparency to be reverted. Other than setting up an update function, is there a way to avoid this problem? I don't know.

figure

% make some fake Z data from some images

A = imread('a.png');

B = imread('b.png');

A = im2double(A)*50 - 25;

B = im2double(B)*100;

% the background axes

contourf(A,'edgecolor','none');

cb1 = colorbar;

h1 = gca;

h2 = axes;

% the foreground axes

[~,hc2] = contourf(B,'edgecolor','none');

set(h2,'color','none') % make the axes background itself invisible

cb2 = colorbar;

% store the position and limits because adjusting

% the colorbar will mess these up

axpos = get(h1,'position');

xl = get(h1,'xlim');

yl = get(h1,'ylim');

% adjust the colorbars so they don't overlap

barl = 0.48;

cb1.Position(2) = cb1.Position(2)+cb1.Position(4)*(1-barl);

cb1.Position(4) = cb1.Position(4)*barl;

cb2.Position(4) = cb2.Position(4)*barl;

% reassert axes geometry so they match

set(h1,'position',axpos,'xlim',xl,'ylim',yl);

set(h2,'position',axpos,'color','none','xlim',xl,'ylim',yl);

% set one colormap for each axes

colormap(h1,turbo)

colormap(h2,summer)

% manipulate transparency via undocumented descendant objects

% this is copied from 1659975, though you may not need the flexibility

% need to force the object to be initialized

% otherwise it may take a second or two before it exists

% all these changes are volatile and will need to be reasserted if the figure is changed

drawnow

filltriangles = hc2.FacePrims;

% set the fill objects to transparent

for k = 1:length(filltriangles)

% Have to set this. The default is 'truecolor' which ignores alpha.

filltriangles(k).ColorType = 'truecoloralpha';

% The 4th element is the 'alpha' value. First 3 are RGB. Note, the

% values expected are in range 0-255.

%filltriangles(k).ColorData(4) = 128;

end

% if you want to make the contour background regions fully transparent

filltriangles(1).ColorData(4) = 0;

With any luck, this junk answer will invoke Cunningham's law.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

ANKAN SARKAR 2023-11-28

Hello @DGM,

Thank you for your answer. I finally solved my problem by replacing the zero and umwanted values with NaN. So, the MATLAB code I shared before is working.

请先登录，再进行评论。