不使用hist（）Matlab的灰度标准化直方图

Question

: a grayscale img in [0..255] ：[0..255]中的灰度img
: img histogram normalized - an array 1X256 divided by total number of pixels ：img histogram标准化 - 数组1X256除以像素总数

这是我的解决方案：

function [h] = histImage(img)
    h=zeros(1,256)
    for i=1:size(h,2)
       h(i) = length(find(img==i));
    end
    h = h./sum(h);

有没有更好的方法呢？

Answer 1

“更好”总是在旁观者的眼中。 无论如何，这是使用accumarray进行上述操作的accumarray ：

%# each pixel contributes 1/nPixels to the counts of the histogram
%# to use graylevels as indices, we have to add 1 to make it 1...256

nPix = numel(img);
h = accumarray(img(:)+1,ones(nPix,1)/nPix,[256 1],@sum,0);