豆粑粑 matlab scatterbar3 bar3 画3D的bar图

meatball1982

论坛里别人的问题。
http://www.ilovematlab.cn/thread-540258-1-1.html

how to use bar3 to plot the 3D data with special x,y,z.

Got a function on file exchange.
scatterbar3
Good work, but not great. I modify it by the facealpha option, the big abs(z), get a clear face. the smaller z will be more transparent.

Use peaks function as a example

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. [x,y,z]=peaks(20);
   
6. X=x(:);
   
7. Y=y(:);
   
8. Z=z(:);
   
9. scatterbar3(X,Y,Z,0.16)
   
10. colormap(jet)
    
11. axis tight
    
12. view(-30,40)
    
13. grid on
    
14. box on
    

复制代码

1. function scatterbar3(X,Y,Z,width)
   
2. %SCATTERBAR3   3-D scatter bar graph.
   
3. %   SCATTERBAR3(X,Y,Z,WIDTH) draws 3-D bars of height Z at locations X and Y with width WIDTH.
   
4. %
   
5. %   X, Y and Z must be of equal size.  If they are vectors, than bars are placed
   
6. %   in the same fashion as the SCATTER3 or PLOT3 functions.
   
7. %
   
8. %   If they are matrices, then bars are placed in the same fashion as the SURF
   
9. %   and MESH functions.
   
10. %
    
11. %   The colors of each bar read from the figure's colormap according to the bar's height.
    
12. %
    
13. %   NOTE:  For best results, you should use the 'zbuffer' renderer.  To set the current
    
14. %   figure renderer to 'zbuffer' use the following command:
    
15. %
    
16. %       set(gcf,'renderer','zbuffer')
    
17. %
    
18. %    % EXAMPLE 1:
    
19. %    y=[1 2 3 1 2 3 1 2 3];
    
20. %    x=[1 1 1 2 2 2 3 3 3];
    
21. %    z=[1 2 3 6 5 4 7 8 9];
    
22. %    scatterbar3(x,y,z,1)
    
23. %    colorbar
    
24. %
    
25. %    % EXAMPLE 2:
    
26. %    [X,Y]=meshgrid(-1:0.25:1);
    
27. %    Z=2-(X.^2+Y.^2);
    
28. %    scatterbar3(X,Y,Z,0.2)
    
29. %    colormap(hsv)
    
30. %
    
31. %    % EXAMPLE 3:
    
32. %    t=0:0.1:(2*pi);
    
33. %    x=cos(t);
    
34. %    y=sin(t);
    
35. %    z=sin(t);
    
36. %    scatterbar3(x,y,z,0.07)
    
37. 
    
38. % By Mickey Stahl - 2/25/02
    
39. % Engineering Development Group
    
40. % Aspiring Developer
    
41. 
    
42. 
    
43. % mod by : mm
    
44. % mod on : 2018.04.06
    
45. % contact me if you get problem about this function
    
46. % meatball1982@163.com
    
47. 
    
48. [r,c]=size(Z);
    
49. 
    
50. z_max=max(abs(Z(:)));
    
51. for j=1:r,
    
52.     for k=1:c,
    
53.         if ~isnan(Z(j,k))
    
54.             drawbar(X(j,k),Y(j,k),Z(j,k),z_max,width/2)
    
55.         end
    
56.     end
    
57. end
    
58. 
    
59. 
    
60. 
    
61. zlim=[min(Z(:)) max(Z(:))];
    
62. if zlim(1)>0,zlim(1)=0;end
    
63. if zlim(2)<0,zlim(2)=0;end
    
64. axis([min(X(:))-width max(X(:))+width min(Y(:))-width max(Y(:))+width zlim])
    
65. caxis([min(Z(:)) max(Z(:))])
    
66. 
    
67. function drawbar(x,y,z,z_max,width)
    
68. edge_col='w';
    
69. % fac_alp=0.6;
    
70. fac_alp = abs(z/z_max)*0.65+0.3;
    
71. h(1)=patch([-width -width width width]+x,[-width width width -width]+y,[0 0 0 0],'b','edgecolor',edge_col,'facealpha',fac_alp);
    
72. h(2)=patch(width.*[-1 -1 1 1]+x,width.*[-1 -1 -1 -1]+y,z.*[0 1 1 0],'b','edgecolor',edge_col,'facealpha',fac_alp);
    
73. h(3)=patch(width.*[-1 -1 -1 -1]+x,width.*[-1 -1 1 1]+y,z.*[0 1 1 0],'b','edgecolor',edge_col,'facealpha',fac_alp);
    
74. h(4)=patch([-width -width width width]+x,[-width width width -width]+y,[z z z z],'b','edgecolor',edge_col,'facealpha',fac_alp);
    
75. h(5)=patch(width.*[-1 -1 1 1]+x,width.*[1 1 1 1]+y,z.*[0 1 1 0],'b','edgecolor',edge_col,'facealpha',fac_alp);
    
76. h(6)=patch(width.*[1 1 1 1]+x,width.*[-1 -1 1 1]+y,z.*[0 1 1 0],'b','edgecolor',edge_col,'facealpha',fac_alp);
    
77. set(h,'facecolor','flat','FaceVertexCData',z,'facealpha',fac_alp)
    

复制代码