myhisteq.m

该程序实现对直方图均衡化

function pimproved=myhidteq(pmat)
%直方图均衡化，输入为原始图像的灰度值，输出为均衡化的图像灰度值

pmax=max(max(pmat));   %原始图像最大的灰度值
[row col]=size(pmat);    %原始图像的大小
s=int(log2(pmax));   %原始图像灰度级个数L-1
%nf用来统计各灰度级的像素数，pin_grade用来记录图像各点所在的灰度级
nf=zeros(s+1);   
pin_grade=zeros(row,col);  
for m=1:s+1
    for i=1:row
        for j=1:col
            if(pmat(i,j)>=2^(m-1))&(pmat(i,j)<2^m)
                nf(m)=nf(m)+1;   %统计各灰度级的像素数
                pin_grade(i,j)=m-1;   %记录图像各点所在的灰度级
            end
        end
    end
end
%统计原始图像直方图个灰度级的频数
pf=zeros(s+1);
for m=1:s+1
    pf(m)=nf(m)/sum(nf);
end
%计算累计分布函数
cf=cumsum(pf);
%灰度级映射时用到的常数gmax和gmin
gmax=2^(s+1)-1;
gmin=0;
%灰度级映射公式
g=int((gmax-gmin)*cf+gmin+0.5);
%统计映射后各灰度级的像素数目
ng=zeros(s+1);
for m=1:s+1
    for i=1:s+1
        if(g(i)==m)
            ng(m)=ng(m)+nf(i);
        end
    end
end
%输出图像直方图的频数
pg=ng./sum(ng);
%计算输出图像各点的灰度值
pimproved=pmat+(2^g(pin_grade)-2^pin_grade);