1、统计各灰度级的像素数
nk(1),nk(2),nk(3),...,nk(256)
MATLAB中输入如下语句：
%求nk
nk=zeros(1,256);
for i=0:255
num=find(I==(i+1));%I为原图像的矩阵
nk(i+1)=length(num);
end
2、求概率pr
MATLAB输入：
%求pr（rk）=nk/MN
pr=zeros(1,256);
for i=1:256
pr(i)=nk(i)/(width*height);
end
3、对概率从前至后依次累加
s(1) s(2) ...s(256)
MATLAB格式：
s=zeros(1,256);
for i=1:L256
for j=1:i
s(i)=s(i)+pr(j);
end
end
4、将分布函数值量化到[0,255]中对应的灰度级
MATLAB语言：
sk=zeros(1,256);
for i=1:256
sk(i)=round((256-1)*sk(i));
end
5、求灰度级出现的概率
MATLAB语言：
for i=0:255
pr(i+1)=sum(s(find(sk==i)));
end
pr %显示pr值
6、将原图像各灰度级的点的灰度改为新的灰度
MATLAB语言：
J=I;
for i=1:256
J(find(J==(i-1)))=sk(i);
end
