Code computes histogram equalization not similar to Matlab's histeq function

Question

Mohsin Shah 2018-5-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/399169-code-computes-histogram-equalization-not-similar-to-matlab-s-histeq-function

回答： xiao 2018-5-11

I have a written a simple code for histogram equalization based on this formula: HistEq = (max gray leve/total no of pixels)*cumulative sum(count). I wonder why the result is quite different from Matlab's function histeq? Can someone please figure out where I am wrong?

   clear;
   clc;
  % Input image
Img = imread('lena.bmp');
I = double(Img);
[H, W] = size(I);
total = H * W; 
count = imhist(Img); 
maxL = 255;
csum = 0;  % Initialize to zero
cumuSum = zeros(1, length(count))';
for idx = 1 : length(count)
    csum = csum + count(idx);
    cumuSum(idx) = csum;
end
HistEq = round(maxL/total * cumuSum);
HIm = zeros(H, W);
for idx = 1 : H
    for idy = 1 : W
        HIm(idx, idy) = HistEq(Img(idx, idy));
    end
end