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]);