histograms for ND data

Versión (2.44 KB) por Diederick
Create histograms for ND data using ND bins

549 descargas

Actualizado 21 Jun 2011

Ver licencia

histnd Histogram count of ND data with ND bins

histmat = histnd(x, y, ..., xedges, yedges, ...)
Extract ND histogram data containing the number of events
of [x, y, ...] tuples that fall in each bin of the ND-grid defined by
xedges, yedges, .... The edges are passed to histc internally and should
therefore conform to histc's input restrictions: the edge-vectors should
be monotonically non-decreasing.

[histmat, nOOF, OOFidx] = histnd(x, y, ..., xedges, yedges, ...)
If any values are outside of the range of the edges, they are not
counted. The number of those cases and their linear index in the input
data is however returned in the second and third output arguments.


events = 1000000;
x1 = sqrt(0.05)*randn(events,1)-0.5; x2 = sqrt(0.05)*randn(events,1)+0.5;
y1 = sqrt(0.05)*randn(events,1)+0.5; y2 = sqrt(0.05)*randn(events,1)-0.5;
x= [x1;x2]; y = [y1;y2];

For linearly spaced edges:
xedges = linspace(-1,1,64); yedges = linspace(-1,1,64);
histmat = histnd(x, y, xedges, yedges);
figure; pcolor(xedges,yedges,histmat'); colorbar ; axis square tight ;

For nonlinearly spaced edges:
xedges_ = logspace(0,log10(3),64)-2; yedges_ = linspace(-1,1,64);
histmat_ = histnd(x, y, xedges_, yedges_);
figure; pcolor(xedges_,yedges_,histmat_'); colorbar ; axis square tight ;

3D data
x = 3.*randn(640000,1);
y = 1.*randn(640000,1);
z = 1.*randn(640000,1);
histmat = histnd(x,y,z,linspace(min(x),max(x),20),linspace(min(y),max(y),20),linspace(min(z),max(z),20));
% make 3D hist, color of points indicates count
[xp,yp,zp] = meshgrid(linspace(min(x),max(x),20),linspace(min(y),max(y),20),linspace(min(z),max(z),20));
% cut away histogram positions where count is 0
qzero = histmat==0;
histmat(qzero) = [];
xp(qzero) = [];
yp(qzero) = [];
zp(qzero) = [];
% draw points
figure;%('Renderer','OpenGL') % might need the openGL renderer to handle so many points
ax = scatter3(xp(:),yp(:),zp(:),'.');
% color them according to count
cdata = histmat(:)./max(histmat);
xlabel('X'), ylabel('Y'), zlabel('Z')

Citar como

Diederick (2023). histograms for ND data (https://www.mathworks.com/matlabcentral/fileexchange/31889-histograms-for-nd-data), MATLAB Central File Exchange. Recuperado .

Compatibilidad con la versión de MATLAB
Se creó con R2011a
Compatible con cualquier versión
Compatibilidad con las plataformas
Windows macOS Linux
Más información sobre Histograms en Help Center y MATLAB Answers.

Inspirado por: 2D Histogram Calculation

Community Treasure Hunt

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

Start Hunting!
Versión Publicado Notas de la versión