How do I overlay contour plots with a different color map for each overlay

I have been able to overlay contours using Matlab 2012b (see the example below) using freezecolors to allow each overlay to have a constant color that is different than the other overlays. But in 2015b freezecolors no longer works. According to this, I shouldn't need freezecolors anymore: http://www.mathworks.com/matlabcentral/fileexchange/7943-freezecolors---unfreezecolors

So I found this, but can't get it to work for what I am trying to do:

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

Any help to get this to work in 2015b would be greatly appreciated!!!

[Carl]

Hi Robert, can you clarify why the second link you posted does not work for you? It seems that you should be able to overlay your contour plots, each on a separate axes with a different colormap.

Answer 1

[KSSV]

Z = peaks(100) ;

% plot valaues 1 
figure
hold on
colors = {'r' 'b' 'g' 'y' 'm' 'k' 'c' 'g' 'r' 'k'} ;
for i = 1:10
    contour(Z,[i i],colors{i},'ShowText','on') ;
end

[Robert]

Thank you KSSV. Your solution applies to a common surface, Z. In the example that I provided, each color represents a separate surface. I.e. Z1 is black, Z2 is red, etc. so that you have multiple surfaces (represented as contour maps) overlaid on each other. But your solution still applies since I can just add for loops for each Z.

[Robert]

Ok. I've selected this as the accepted answer, but now I have another problem. This works if all of the contours can be plotted on the same z-scale, but what if my surfaces have different orders of magnitude? So now I need to have multiple z-scales. Any help?

[KSSV]

I think it shall work even if the orders are different. If it doesn't then by interpolation we can get all surfaces to same orders.

[Robert]

I played around some more with the second link in my original post and got this to work using multiple axes.

nMdls = 2;
Axes = zeros(nMdls);
Axes(1) = axes;
%%Create the first axis and plot the contours
[x,y,z] = peaks;
hold 'on'
for i = 1:2:10
  contour(Axes(1),x,y,z, [i, i],'r', 'ShowText', 'on')
end
%%Create the second axis and plot the contours
Axes(2) = axes;

z2 = z * 1000; %Change the order of magnitude of the orignal z-data hold 'on' for i = 2000:2000:10000 contour(Axes(2),x,y,z2, [i, i],'b', 'ShowText', 'on'); end %% linkaxes(Axes) %% Hide the top axes set(Axes(2), 'visible', 'off'); set(Axes(2), 'XTick', []); set(Axes(2), 'YTick', []);

[Robert]

Sorry, the code should look like this:

nMdls = 2;
Axes = zeros(nMdls);
Axes(1) = axes;
%%Create the first axis and plot the contours
[x,y,z] = peaks;
hold 'on'
for i = 1:2:10
  contour(Axes(1),x,y,z, [i, i],'r', 'ShowText', 'on')
end
%%Create the second axis and plot the contours
Axes(2) = axes;
z2 = z * 1000; %Change the order of magnitude of the orignal z-data
hold 'on'
for i = 2000:2000:10000
  contour(Axes(2),x,y,z2, [i, i],'b', 'ShowText', 'on');
end
linkaxes(Axes)
%%Hide the top axes
set(Axes(2), 'visible', 'off');
set(Axes(2), 'XTick', []);
set(Axes(2), 'YTick', []);

[KSSV]

Using axes is suggested then freezecolors....Good