How to bin 2d data?

36 次查看(过去 30 天)
Yeping Sun
Yeping Sun 2016-8-6
编辑: Thorsten 2016-8-9
Dear all,
I have a 2-dimension coordinate set (x,y) which contains 5000 points and I want to construct a 2d map to reflect the point density distribution. firstly I want to divide both of the ranges of the two coordinate sets (-3 to2 for x; and -2.5 to 3 for y) into like 37 bins,so there will be totally 37x37 bins. Then I want to assign each point into these bins.To do that I need to construct a indicator function δkij (k=1:5000; i=1:37; j=1:37) which denotes whether point k falls in bin (i,j): if point k falls into bin (i,j), then δkij=1, otherwiseδkij=0. Then the histogram at bin (i,j) is the sum of δkij (k=1:5000). At last I wish to plot a contour map or hotmap to show the hisrogram. How does matlab realize these computation?
All the best!
Yeping Sun
  1 个评论
Azzi Abdelmalek
Azzi Abdelmalek 2016-8-6
Your question is not clear, post a short example, give a brief explanation on what you want, then post the expected result

请先登录,再进行评论。

采纳的回答

Image Analyst
Image Analyst 2016-8-6
编辑:Image Analyst 2016-8-6
Try (the badly named) hist3() to construct a 2-D histogram:
% Create sample data.
numPoints = 5000;
x = -3 + 5 * rand(numPoints, 1);
y = -2.5 + 5.5 * rand(numPoints, 1);
% Now construct the histogram
xy = [x, y];
subplot(1,2,1);
hist3(xy, [37, 37]); % Plots bar chart.
subplot(1,2,2);
counts = hist3(xy, [37, 37]);
imshow(counts, []); % Show as image
colormap(hot(256));
colorbar;
axis on;
  8 个评论
Yeping Sun
Yeping Sun 2016-8-8
Hi! I actually want the axes on, but I want them to be in the range of -3 to 3. Could you more specific when you said "you can send in the values to imshow() with the xaxis and y axis options"? I've tried
imshow(heatMap, [-3 3]);
it doesn't work.
By adding your code for the kernel thing, now the image becomes (the attached). It colorbar is still mad.
Image Analyst
Image Analyst 2016-8-9
SOrry, I meant XData and YData, like this:
imshow(heatMap, 'XData', [-3 3], 'YData', [-3,3]);
but you should really use histogram2() and get the BinEdges property and use that instead of [-3, 3] because [-3,3] may not be correct if the histogram was automatic. However if you passed in the edges, then you know what the edges are and you can use that.

请先登录,再进行评论。

更多回答(2 个)

Steven Lord
Steven Lord 2016-8-6
Use the histogram2 function. The DisplayStyle 'tile' looks like it's exactly what you're looking for.

Thorsten
Thorsten 2016-8-9
编辑:Thorsten 2016-8-9
Have a look at this:
D = dlmread('/Users/hansen/Downloads/PC1-PC2.txt');
Nbins = 37;
H = hist3(D, {linspace(-3,3,Nbins) linspace(-3,3,Nbins)});
imshow(H, [])
colormap hot
colorbar
The values of H are in the range [0, 196], and that is the range of the colorbar.
For a contour plot, use:
contour(H)

类别

Help CenterFile Exchange 中查找有关 Histograms 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by