scatterbar3.m

function scatterbar3(X,Y,Z,width)
%SCATTERBAR3   3-D scatter bar graph.
%   SCATTERBAR3(X,Y,Z,WIDTH) draws 3-D bars of height Z at locations X and Y with width WIDTH.
%
%   X, Y and Z must be of equal size.  If they are vectors, than bars are placed
%   in the same fashion as the SCATTER3 or PLOT3 functions.
%
%   If they are matrices, then bars are placed in the same fashion as the SURF
%   and MESH functions.
%
%   The colors of each bar read from the figure's colormap according to the bar's height.
%
%   NOTE:  For best results, you should use the 'zbuffer' renderer.  To set the current
%   figure renderer to 'zbuffer' use the following command:
%
%       set(gcf,'renderer','zbuffer')
%
%    % EXAMPLE 1:
%    y=[1 2 3 1 2 3 1 2 3];
%    x=[1 1 1 2 2 2 3 3 3];
%    z=[1 2 3 6 5 4 7 8 9];
%    scatterbar3(x,y,z,1)
%    colorbar
%
%    % EXAMPLE 2:
%    [X,Y]=meshgrid(-1:0.25:1);
%    Z=2-(X.^2+Y.^2);
%    scatterbar3(X,Y,Z,0.2)
%    colormap(hsv)
%
%    % EXAMPLE 3:
%    t=0:0.1:(2*pi);
%    x=cos(t);
%    y=sin(t);
%    z=sin(t);
%    scatterbar3(x,y,z,0.07)

% By Mickey Stahl - 2/25/02
% Engineering Development Group
% Aspiring Developer


% mod by : mm
% mod on : 2018.04.06
% contact me if you get problem about this function
% meatball1982@163.com

[r,c]=size(Z);

z_max=max(abs(Z(:)));
for j=1:r,
    for k=1:c,
        if ~isnan(Z(j,k))
%             drawbar(X(j,k),Y(j,k),Z(j,k),z_max,width/2)
              drawbar(X(j,k),Y(j,k),Z(j,k),z_max,width/2)
        end
    end
end



zlim=[min(Z(:)) max(Z(:))];
if zlim(1)>0,zlim(1)=0;end
if zlim(2)<0,zlim(2)=0;end
axis([min(X(:))-width max(X(:))+width min(Y(:))-width max(Y(:))+width zlim])
caxis([min(Z(:)) max(Z(:))])

function drawbar(x,y,z,z_max,width)
edge_col='w';
% fac_alp=0.6;
fac_alp = abs(z/z_max)*0.2+0.2;
h(1)=patch([-width -width width width]+x,[-width width width -width]+y,[0 0 0 0],'b','edgecolor',edge_col,'facealpha',fac_alp);
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);
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);
h(4)=patch([-width -width width width]+x,[-width width width -width]+y,[z z z z],'b','edgecolor',edge_col,'facealpha',fac_alp);
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);
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);
set(h,'facecolor','flat','FaceVertexCData',z,'facealpha',fac_alp)