豆粑粑 matlab scatterbar3 bar3 画3D的bar图
论坛里别人的问题。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
clear all
clc
clf
=peaks(20);
X=x(:);
Y=y(:);
Z=z(:);
scatterbar3(X,Y,Z,0.16)
colormap(jet)
axis tight
view(-30,40)
grid on
box on
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=;
% x=;
% z=;
% scatterbar3(x,y,z,1)
% colorbar
%
% % EXAMPLE 2:
% =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
=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)
end
end
end
zlim=;
if zlim(1)>0,zlim(1)=0;end
if zlim(2)<0,zlim(2)=0;end
axis()
caxis()
function drawbar(x,y,z,z_max,width)
edge_col='w';
% fac_alp=0.6;
fac_alp = abs(z/z_max)*0.65+0.3;
h(1)=patch([-width -width width width]+x,[-width width width -width]+y,,'b','edgecolor',edge_col,'facealpha',fac_alp);
h(2)=patch(width.*[-1 -1 1 1]+x,width.*[-1 -1 -1 -1]+y,z.*,'b','edgecolor',edge_col,'facealpha',fac_alp);
h(3)=patch(width.*[-1 -1 -1 -1]+x,width.*[-1 -1 1 1]+y,z.*,'b','edgecolor',edge_col,'facealpha',fac_alp);
h(4)=patch([-width -width width width]+x,[-width width width -width]+y,,'b','edgecolor',edge_col,'facealpha',fac_alp);
h(5)=patch(width.*[-1 -1 1 1]+x,width.*+y,z.*,'b','edgecolor',edge_col,'facealpha',fac_alp);
h(6)=patch(width.*+x,width.*[-1 -1 1 1]+y,z.*,'b','edgecolor',edge_col,'facealpha',fac_alp);
set(h,'facecolor','flat','FaceVertexCData',z,'facealpha',fac_alp)
页:
[1]