Generate 3D histogram while capping the height of bars

Question

Ken 2024-4-1，16:59

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2101361-generate-3d-histogram-while-capping-the-height-of-bars

评论： Adam Danz 2024-4-10，16:56

I'm generating 3D histograms. It often happens that, given the nature of our data, one bar towers above all the rest. Since the histogram automatically scales to the height of the tallest bar, this makes it hard to see the trends that we would really like to see in the shorter bars. So I would like to generate a histogram (color-coded by height, that's really nice) with the height of individual bars capped. In other words, if I set the cap at 1000 and the computed height of a bar is 4000, it would display with a height of 1000.

It doesn't look like hist3 has a zlimit option. However, this answer contains the following fascinating line of code, a feature I don't see documented in the Matlab documentation, that calculates the matrix of bin values rather than plotting the actual histogram. I have tested, it actually does.

cnt = hist3(X,{bin_x,bin_y});

It would be very easy to apply a floor function to the matrix thus calculated.

If this is documented, could someone please point me to the place? And regardless, can I call hist3 in some other way to plot the contents of this matrix once generated and modified?

2 个评论
显示无隐藏无

Dyuman Joshi 2024-4-1，17:09

编辑：Dyuman Joshi 2024-4-1，17:17

That feature is documented in the documentation of hist3 - check the second last syntax.

You could use max() to compare and cap the number accordingly.

Yes, you can call hist3() after manipulating the data, if that is what you meant.

Ken 2024-4-1，21:29

编辑：Ken 2024-4-1，21:30

It looks like I cannot call hist3 on the modified output, because it is not prepared to accept a matrix of bin values, but only a two column matrix of X, Y pairs. I can call bar3, but unfortunately it seems to handle color very differently than hist3 does. 'FaceColor', 'interp' is not accepted, alas!

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Steven Lord 2024-4-1，22:31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2101361-generate-3d-histogram-while-capping-the-height-of-bars#answer_1434511

在 MATLAB Online 中打开

Rather than trying to make hist3 or histogram2 truncate the bins, consider using zlim to control the limits of the Z axis.

x = randn(1, 1e5);

y = randn(1, 1e5);

histogram2(x, y, FaceColor = 'flat'); % Note the center bars have height > 400

figure

histogram2(x, y, FaceColor = 'flat')

zlim([0 400]) % Limit the axes to showing only up to 400

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Adam Danz 2024-4-1，23:06

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2101361-generate-3d-histogram-while-capping-the-height-of-bars#answer_1434526

编辑：Adam Danz 2024-4-1，23:14

在 MATLAB Online 中打开

How to cap the height of a 2D histogram without affecting the max color for each bar

This is tricky because hist3 does not return the handle to the surface object it creates. Follow these steps to clip the height of the 3D bars while maintaining the bar color indicating the original heights (if I understood the question correctly).

Note that this produces a misleading plot where bar height and bar color no longer agree which may cause confusion or misinterpretation of the results. Assuming there is good reason to do this, it should be clearly indicated in text that bar height is capped.

% Produce demo data
rng default
x = randn(1000,1);
y = randn(1000,1);
% Plot the 3D historam
figure()
hist3([x,y],'NBins',[5,5],'CDataMode','auto','FaceColor','interp')
colorbar()

Note that the tallest bar has a height of 256. We will cap the height at z=200.

% Define the max height

zClip = 200;

% Get the surface handle - this assumes there is only 1 surface in the axes

ax = gca();

s = findobj(ax,'type','surface');

% Prevent CData from updating when the ZData changes. This will also

% prevent CLim from updating.

s.CDataMode = 'manual';

% Replace ZData larger than zClip

s.ZData = min(s.ZData, zClip);

subtitle(sprintf('Bar height is capped at %.1f',zClip))

Now the 3D histogram is capped at z = 200 while the color data remains at its original values showing a max height of 256.

2 个评论
显示无隐藏无

Ken 2024-4-9，11:28

Wow, snagging the surface object out of the axes? Pretty snazzy!

But actually I would like the colors to rescale. I am dealing with a situation in which one tall bar can completely dominate the picture, thereby making almost all other bars almost the same color.