histeqs.asv

基于内容的图像检索算法的研究

function [Image] = HISTEQS( Image )
%HISTEQS 实现直方图均衡化
%程序的功能：实现直方图均衡化的程序。
%程序当中用到的公式：
%     灰度级为[0,L-1]范围的数字图像的直方图是离散函数h(rk)=nk;
%     其中，rk是第k级灰度，nk是图像中灰度级为rk的像素个数
%     灰度rk出现的概率为：p(rk)=nk/n
%     其中，n是图像的像素总数
%程序的参数说明：
%


h=IMHISTS(Image);
p=zeros(1,256);    %一行256列的数组，用于表示相应灰度级出现的概率
s=p;
p=h./numel(Image);
s=cumsum(p); 
s=double(im2uint8(s));

[M,N]=size(Image);
for m=1:M
    for n=1:N
       Image(m,n)=s(double(Image(m,n))+1);
     end
 end
r=double(Image);
s1=r.^1.8;                            %幂变换
Image=im2uint8(mat2gray(s1));