直方图均衡

function Hisequalization()

%读取图像

input_image=imread( './a.bmp' ,   'bmp' );

F=double(input_image);

[m,n]=size(F);

figure;

subplot(221);

imshow(uint8(F));

title( '????' );

subplot(222);

imhist(uint8(F));

His=zeros(1,256); %

for   x=1:m

      for   y=1:n

        M(x,y)=F(x,y);

        a=F(x,y);

        His(a+1)=His(a+1)+1;

      end

end

for   x=1:256

    normHis(x)=His(x)/(m*n); %????????????????

      if (x~=1)

        cumuHis(x)=cumuHis(x-1)+normHis(x); %??????????

      else

        cumuHis(x)=normHis(x);

      end

end

for   i=1:256

    Map(i)=floor(cumuHis(i)*255+0.5); %????

end

 

for   x=1:m

      for   y=1:n

        tmp=M(x,y)+1;

        M(x,y)=Map(tmp);

      end

end

subplot(223);

imshow(uint8(M));

title( '??????' );

subplot(224);

imhist(uint8(M));

%imwrite(uint8(M),'./ttttttt.bmp');

 

 

对于彩色图：

function Hisequalization2()

input_image=imread('./test.png', 'png');

F=double(input_image);

[m,n,o]=size(F);

figure;

subplot(241);

imshow(uint8(F));

title('????');

for z=1:3

    for x=1:m

        for y=1:n

            m(x,y)=F(x,y,z);

        end

    end

    subplot(2,4,1+z);

    imhist(uint8(m));

end

for z=1:3

  

His=zeros(1,256);

for x=1:m

    for y=1:n

        M(x,y,z)=F(x,y,z);

        a=F(x,y,z);

        His(a+1)=His(a+1)+1;

    end

end

sizeof(His);

normHis=zeros(1,256);

for x=1:256

    normHis(x)=His(x)/(m*n);

    if(x~=1)

        cumuHis(x)=cumuHis(x-1)+normHis(x);

    else

        cumuHis(x)=normHis(x);

    end

end

for i=1:256

    Map(i)=floor(cumuHis(i)*255+0.5);

end

 

for x=1:m

    for y=1:n

        tmp=M(x,y,z)+1;

        M(x,y,z)=Map(tmp);

    end

end

for x=1:m

    for y=1:n

       

        m(x,y)=M(x,y,z);

    end

end

subplot(2,4,z+5);

imhest(uint8(m));

end

subplot(225);

imshow(uint8(M));

title('??????');

转载于:https://www.cnblogs.com/wenning/archive/2012/05/07/2487220.html