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