豆粑粑，matab 我画你个饼 ，3D等值面

meatball1982

之前，陌生女人的提问，
解答过程中觉得思路还不错（当时是在车上，临时的一个小程序。
今天，突然想到可以用patch的方法画各个3D的等值面。
效果，还不错。

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. %% outline
   
6. % plot BING
   
7. % contact meatball1982@163.com
   
8. % if u really, really got problem with this m file.
   
9. 
   
10. 
    
11. %% main
    
12. % generate data
    
13. [x,y,z]=peaks(100);
    
14. % get the contour line
    
15. [C,h]=contour(x,y,z,18);
    
16. 
    
17. % get the structure, close the contour
    
18. S=contourdata(C);
    
19. clf
    
20. 
    
21. 
    
22. %% plot BING
    
23. for i = 1: length(S)
    
24.     s_lev(i,1)=S(i).level;
    
25. end
    
26. 
    
27. tm_s = unique(s_lev);
    
28. n_col = length(tm_s);
    
29. [loca,locb]=ismember(s_lev,tm_s); % get the lin color position in col_mm
    
30. col_mm=jet(n_col);
    
31. 
    
32. %% outplut
    
33. % the normal one
    
34. subplot(1,2,1)
    
35. surf(x,y,z,'edgecolor','none')
    
36. colormap(col_mm)
    
37. view(-20,20)
    
38. axis tight
    
39. colorbar
    
40. % grace BING plot
    
41. subplot(1,2,2)
    
42. hold on
    
43. for i = 1:length(S)
    
44.     x_lin = S(i).xdata;
    
45.     y_lin = S(i).ydata;
    
46.     z_va  = S(i).level;
    
47.     z_lin = z_va.*ones(size(x_lin));
    
48.     alp_va = abs(z_va./max(abs(z(:))))/2+0.5;
    
49. %     h=patch(x_lin,y_lin,z_lin,col_mm(locb(i),:),...
    
50. %                                          'edgecolor',col_mm(locb(i),:)); % no transpancy
    
51.     h=patch(x_lin,y_lin,z_lin,col_mm(locb(i),:),...
    
52.                                           'edgecolor',col_mm(locb(i),:),...
    
53.                                           'facealpha',alp_va); % transpancy
    
54. end
    
55. axis tight
    
56. view(-20,20)
    
57. grid on

复制代码