equalisasi.m

来自「to flip image in matlab」· M 代码 · 共 43 行

function Equalisasi
global IMG;
image=double(IMG);
s=size(image);
result=zeros(s(1),s(2),s(3));
volume=s(1)*s(2);
for i=1:256
    histR(i)=0;
    histG(i)=0;
    histB(i)=0;
end
for i=1:s(1)
    for j=1:s(2)
        histR(image(i,j,1)+1)=histR(image(i,j,1)+1)+1;
        histG(image(i,j,2)+1)=histG(image(i,j,2)+1)+1;
        histB(image(i,j,3)+1)=histB(image(i,j,3)+1)+1;
    end
end
for i=1:256
    sumR=0;
    sumG=0;
    sumB=0;
for j=1:i
        sumR=sumR+histR(j)/volume;
        sumG=sumG+histG(j)/volume;
        sumB=sumB+histB(j)/volume;
end
    histR_eq(i)=round(255*sumR);
    histG_eq(i)=round(255*sumG);
    histB_eq(i)=round(255*sumB);
end
for i=1:s(1)
    for j=1:s(2)
        result(i,j,1)=histR_eq(image(i,j,1)+1);
        result(i,j,2)=histG_eq(image(i,j,2)+1);
        result(i,j,3)=histB_eq(image(i,j,3)+1);
    end
end
res=uint8(result);
IMG=res;
imshow(IMG)