www.gusucode.com > 图像变换及频域滤波 空间域图像增强 图像边缘检测源码程序 > digitalimageprocessprograms/fing.m

    clc;
fid=fopen('d:/img/fing_128.img','r');
f=fread(fid,[128,128],'uchar');
subplot(2,1,1);
imshow(f,[0,255]);
q=zeros(1,256);
for x=1:128
    for y=1:128
        q(f(x,y)+1)=q(f(x,y)+1)+1;
    end
end
s=q./(128*128);
X=0:255;
subplot(2,1,2);
bar(X,s');


figure;
t=zeros(1,256);
t(1)=s(1);
for i=2:256
    t(i)=t(i-1)+s(i);
end
subplot(2,1,1);
bar(X,t');

t0=floor(255*t+0.5);
subplot(2,1,2);
bar(X,t0');


figure;
t1=zeros(1,256);
for i=1:256
    t1(t0(i)+1)=s(i)+t1(t0(i)+1);
end
subplot(2,1,1);
bar(X,t1');


f1=zeros(128,128)
for x=1:128
    for y=1:128
        f1(x,y)=t0(f(x,y)+1);
    end
end
subplot(2,1,2);
imshow(f1,[0,255]);