The Woes of contourf patches, findobj, and Matlab 2014b.

Question

Kellen 2015-1-12

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/169706-the-woes-of-contourf-patches-findobj-and-matlab-2014b

评论： Lauren Johnson 2015-3-27

Hello all, thanks in advance for help.

With the recent update to OS X Yosemite, and of course the required update to MATLAB 2014b, AND the graphics changes to 2014b, I am experiencing some problems with our legacy code.

Background: We have a script called cpatch.m (see below), which [before updating] used

findobj(gca,'type','patches');

on a contourf plot to recolor the patches so that they were a specific color corresponding to a specific contour interval (necessary for the accuracy of our science).

Problem: Contourf in 2014b creates an handle which no longer has any children. It used to be the case that hidden in the children, were the individual patches that one could recolor. Now, these patches do not exist, so one cannot recolor them. In MATLAB:

patches = findobj(gca,'type','patch')

returns

0x0 empty GraphicsPlaceholder array.

and hence is useless.

Ideal Solution: Can someone please tell me how to access the individual patches so that I can change their face color?

Code:

 clev   = 0:.2:1; % These are the contour levels
 colors = bone(length(clev)-1); 
 cpt    = [clev(1:end-1)' clev(2:end)'];
 [C,H] = contourf(rand(10,10),clev);
 patches = findobj(gca,'type','patch');
 for ipatches=1:length(patches)
       c = get(patches(ipatches),'cdata');
       if length(c) <= 1 
           ii = find(cpt(:,1) <= c);
           iwant = max(ii);
           if length(iwant) >= 1
               set(patches(ipatches),'facecolor',colors(iwant,:));
           end
       end
 end

Please help! Thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Kelly Kearney 2015-1-12

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/169706-the-woes-of-contourf-patches-findobj-and-matlab-2014b#answer_164705

在 MATLAB Online 中打开

Can I recommend contourfcmap? I believe it does exactly what your code used to do (assign precise colors to the shaded intervals in a contourf plot), and it works in R2014b.

If you look under the hood of that function, you can find the workaround to extract and recolor. Pre-2014b, contour objects were composed of patches. In 2014b, their underlying polygon components are TriangleStrips, rather than simple patches, and their properties are mostly hidden from the user. But you can still hack into it.

[X,Y,Z] = peaks;
[l, h] = contourf(X,Y,Z,4);
h.FacePrims(1).ColorData = uint8([255 0 0 255]'); % change first interval to red

12 个评论
显示 10更早的评论隐藏 10更早的评论

Kelly Kearney 2015-1-13

在 MATLAB Online 中打开

Interesting. You're right, using datenumbers along the x-axis wreaks havoc with the recolor method in 2014b... I'll have to look into that more. But if you use the calccontour method, neither the dates nor the color-reverting should be an issue.

If you do see color-reverting, the function must be reverting to the default recolor for some reason. Do you have the Mapping Toolbox? If so, perhaps it's a syntax issue? Can you show how you've set up your function call?

I just tested it again with the following example:

z = peaks(100);
t = datenum(2000,1,1) + (1:100);
y = 1:100;
clev = [-5 -3 -2:.5:2 3 5];
cmap = jet(12);
ax(1) = subplot(4,1,1);
[c,h] = contour(t,y,z, clev);
% pcolor(t,y,z); shading flat;
datetick('x');
title('contour');
ax(2) = subplot(4,1,2);
hc = contourfcmap(t,y,z,clev,cmap, 'method', 'calccontour');
datetick('x');
title('countourfcmap, calccontour');
try
    ax(3) = subplot(4,1,3);
    hc = contourfcmap(t,y,z,clev,cmap, 'method', 'recolor');
    datetick('x');
catch
    text(mean(get(ax(3),'xlim')), mean(get(ax(3),'ylim')), ...
        'FAILED','fontsize', 30, 'horiz', 'center');
end
title('contourfcmap, recolor, x=dates');
ax(4) = subplot(4,1,4);
hc = contourfcmap(1:100,y,z,clev,cmap, 'method', 'recolor');
title('contourfcmap, recolor, x=1:nx');
set(ax(1:3), 'xlim', [t(1) t(end)]);
im = frame2im(getframe(gcf));
imwrite(im, 'test1.png', 'png');
export_fig('test2', gcf, '-png', '-nocrop');

Here's the result of the screenshot:

And here's the exported version:

As you can see, I'm replicating the bug with dates and recolor, and also experiencing the color-revert-on-print with the recolor version. But the version using calccontour properly renders both on screen and when printed.

Can you let me know if you see different results running this code snippet? And I'll look into the date bug.

Kellen 2015-1-14

在 MATLAB Online 中打开

Thanks again! Running the code snippets produced identical results.

When I try to run my code, I ran into a bunch of errors, so I ended up having to use no NaN values and set them to a minimal value and allow contourfcmap to use the 'lo' parameter to white them out (which isn't working).

I've attached some made up sample data and a very hacked up version of my code so you can see the results I get. It plots and retains colors which is great, but it seems the plot doesn't line up with the correct values on the x- or y-axes, and the blue box that you will see above the plotted data should be filled with the low value (white), but it is not.

Thanks for all of your help!

Code:

 % Setup Plot
 WD = 7.5; HT=8.5; txtSize=12;
 set(gcf,'DefaultAxesFontSize',txtSize,'DefaultTextFontSize',txtSize)
 set(gcf, 'Units', 'inches', 'PaperPosition', [5,5,WD,HT]);
 set(gcf, 'Position', [5 5 WD HT]);
 lmarg  = .85;  tmarg  = .3;  rmarg  = .5;  bmarg  = 1;  vspace = 0.10;
 height = (HT-4*vspace-tmarg-bmarg)/5;  width  = WD-lmarg-rmarg;
 xshift = lmarg;  yshift = HT - tmarg - height;
 % Load Files
 load ~/Desktop/test_data.mat;
 % Setup axes
 fprintf('Plotting...');
 position = [xshift yshift width height];
 handle = axes('Units','Inches','Position',position);
 % Contour
 HC = contourfcmap(t,y,z,clev,cmap,'lo',[0 0 0],'hi',cmap(end,:),'method','calccontour');
 % Xticks & Text
 datetick('x',28,'keeplimits');
 title('Makeshift Title');
 ylabel('Distance (km)');
 fprintf(' Done!\n');
 % Save Figure
 shg;
 set(gcf,'color','w');
 export_fig ~/Desktop/test1.pdf '-pdf' '-painters'

Result:

Kelly Kearney 2015-1-14

在 MATLAB Online 中打开

I assume the portion that used to be NaNs is that upper triangular bit? When I look at your data, it seems those values are near 0 (-3.1415e-05). But you specify contour levels ranging from -45 to 50, so that value doesn't qualify as low (and instead is being shaded by the light blue you prescribed to the -5-0 interval).

I added the line

z(z == z(end,1)) = -999;

immediately after loading the .mat file, and now when plotted that region will be properly categorized a low.

As for the values lining up, I think you're just referring to the changed axis x-limits. That's just due to differing default behavior between contour and patch; contour sets the axis limits tight to the data, while patch allows a little space around the plotted objects. Adding

axis tight;

before you call datetick should fix that.

With those two small changes, I think the result is what you want:

-

I'll add the NaN problem to my bug list for contourfcmap.

Also, sorry the calccontour method is so slow for this particular example... the polygon-intersection calculation gets pretty resource-intensive with this many contour lines. Might be able to speed that up a lot if I call gpcmex directly rather than using polybool... or perhaps I need to rethink the algorithm there.

Kelly Kearney 2015-3-25

在 MATLAB Online 中打开

If by blocks of data, do you mean the contour lines seem to zig-zag in pixelated manner based on the underlying data resolution?

[x,y,z] = peaks(100);
dn = now;
clev = -8:8;
c = contourcs(x(1,:)+dn,y(:,1),z);
ax(1) = subplot(2,2,1);
contour(x+dn,y,z,clev);
ax(2) = subplot(2,2,2);
contourf(x+dn,y,z,clev);
ax(3) = subplot(2,2,3);
hold on;
arrayfun(@(X) plot(X.X, X.Y), c);
linkaxes(ax, 'x');

If so, then no, I haven't found a way around that behavior, which seems to be a bug in R2014b contour and contourf. Luckily, it does seem to be resolved in 2015a.

With contourfcmap, using the calccontour method will avoid the bug, since it plots lines using the output of contourc, and the zig-zaggy bug doesn't seem to manifest itself in countourc.

Lauren Johnson 2015-3-27