如何在MATLAB中的bar3图中隐藏零值

时间:2010-01-12 16:14:59

标签: matlab plot data-visualization histogram

我已经使用bar3 plot命令生成了二维直方图(绘图是3D - 几个直方图并排绘制的直方图)。但是,所有零值在x-y平面中显示为平面正方形。有没有办法阻止MATLAB显示值?我已经尝试用NaN替换所有零,但它没有改变关于该情节的任何内容。这是我一直在试验的代码:

x1=normrnd(50,15,100,1); %generate random data to test code
x2=normrnd(40,13,100,1);
x3=normrnd(65,12,100,1);

low=min([x1;x2;x3]);
high=max([x1;x2;x3]);
y=linspace(low,high,(high-low)/4); %establish consistent bins for histogram
z1=hist(x1,y);
z2=hist(x2,y);
z3=hist(x3,y);
z=[z1;z2;z3]';
bar3(z)

正如您所看到的,情节上有很多零值。用NaN替换零后关闭数字并重新绘图似乎没有任何改变:

close
z(z==0)=NaN;
bar3(z)

3 个答案:

答案 0 :(得分:13)

一种解决方案是修改bar3创建的图形对象。首先,您必须获取从bar3返回的句柄:

h = bar3(z);

在您的情况下,h将是一个3元素的句柄向量,每组彩色条一个。然后,下面的代码应该使计数为零的二进制文件不可见:

for i = 1:numel(h)
  index = logical(kron(z(:, i) == 0, ones(6, 1)));
  zData = get(h(i), 'ZData');
  zData(index, :) = nan;
  set(h(i), 'ZData', zData);
end

这是一个例子(带有强制性的自由圈):

enter image description here

工作原理......

如果您的bin计数向量为N-by-1,则bar3将绘制6*N个矩形块(即每个bin的长方体的6个面)。 'ZData'中每组补丁对象的h属性将为(6*N)-by-4,因为每个矩形面有4个角。因此,'ZData'属性的6行的每个簇都是一个bin的6个面的z坐标集。

上面的代码首先创建一个逻辑向量,其中bin计数等于0,然后使用kron函数复制该向量的每个元素6次。这将成为'ZData'属性行的索引,此索引用于将空坐标的补丁的z坐标设置为nan。这将导致补丁无法渲染。


修改

这是一个稍微修改过的代码版本,通过从绘制条形的'ZData' property中获取条形高度使其更加通用,因此它所需要的只是从{{3}返回的句柄}。我还将代码包装在一个函数中(没有错误和输入检查):

function remove_empty_bars(hBars)
  for iSeries = 1:numel(hBars)
    zData = get(hBars(iSeries), 'ZData');  % Get the z data
    index = logical(kron(zData(2:6:end, 2) == 0, ones(6, 1)));  % Find empty bars
    zData(index, :) = nan;                 % Set the z data for empty bars to nan
    set(hBars(iSeries), 'ZData', zData);   % Update the graphics objects
  end
end

答案 1 :(得分:7)

这是一个示例,显示如何隐藏具有零值的条形图。我们从正常的BAR3情节开始:

x = 1:7;
Y = jet(numel(x));
h = bar3(x,Y,'detached');
xlabel x; ylabel y; zlabel z; box on;

before

请注意,变量h包含surface个句柄数组(在这种情况下为3个,每个“组”条形一个。这些组对应于Y的列矩阵,每个都用不同的颜色表示。)

现在隐藏零值的代码:

for i=1:numel(h)
    %# get the ZData matrix of the current group
    Z = get(h(i), 'ZData');

    %# row-indices of Z matrix. Columns correspond to each rectangular bar
    rowsInd = reshape(1:size(Z,1), 6,[]);

    %# find bars with zero height
    barsIdx = all([Z(2:6:end,2:3) Z(3:6:end,2:3)]==0, 2);

    %# replace their values with NaN for those bars
    Z(rowsInd(:,barsIdx),:) = NaN;

    %# update the ZData
    set(h(i), 'ZData',Z)
end

after

说明:

对于每组条形图,创建一个surface图形对象(句柄存储在h(i)中)。它的Z坐标矩阵ZData表示为6*N-by-4矩阵(XDataYDataCData矩阵相同),其中N是数字每组中的矩形条或上例中的7条。

这样,每个矩形用6×4矩阵表示(每个矩阵对应一个X / Y / Z坐标)。例如,一个这样的矩形的坐标看起来像:

>> xx = get(h(3),'XData'); yy = get(h(3),'YData'); zz = get(h(3),'ZData');

>> xx(1:6,:)
ans =
          NaN          2.6          3.4          NaN
          2.6          2.6          3.4          3.4
          2.6          2.6          3.4          3.4
          NaN          2.6          3.4          NaN
          NaN          2.6          3.4          NaN
          NaN          NaN          NaN          NaN

>> yy(1:6,:)
ans =
          NaN          0.6          0.6          NaN
          0.6          0.6          0.6          0.6
          1.4          1.4          1.4          1.4
          NaN          1.4          1.4          NaN
          NaN          0.6          0.6          NaN
          NaN          NaN          NaN          NaN

>> zz(1:6,:)
ans =
          NaN            0            0          NaN
            0            1            1            0
            0            1            1            0
          NaN            0            0          NaN
          NaN            0            0          NaN
          NaN          NaN          NaN          NaN

每个第二列跟踪左侧的点,第三列跟踪右侧的点,当两者连接时将绘制矩形的4个面:

>> surface(xx(1:6,2:3), yy(1:6,2:3), zz(1:6,2:3), cc(1:6,2:3))
>> view(3)

rectangle_surface

第一列和最后一列将通过关闭矩形的边来绘制剩余的两个面。

所有这些矩阵都连接成一个高矩阵,矩形全部使用单个表面对象绘制。这是通过使用NaN值来分隔同一矩形的点内部和差异矩形之间的不同部分来实现的。

所以上面的代码所做的是查找Z高度为零的矩形,并用NaN值替换所有值,这有效地告诉MATLAB不要绘制由这些点形成的表面。

答案 2 :(得分:1)

我的问题不是零值,而是NaN值(在bar3内转换为零值)。 我想继续显示值为零的元素,但不显示值为nan的元素。 我稍微调整了代码,它运行得很好:

for i = 1:numel(h)
  index = logical(kron(isnan(z(:,i)),ones(6,1)));
  zData = get(h(i),'ZData');
  zData(index,:) = nan;
  set(h(i),'ZData',zData);
end

谢谢!