clc?
clear?all?
%讀入圖片?
gf=imread（‘5.jpg‘）;
figureimshow（gf）;
title（‘原始圖像‘）;
a=rgb2gray（gf）;
figureimhist（a）?
title（‘原始直方圖‘）;?
gfMin=min（a（:））;?
gfMax=max（a（:））;?
R=gf（::1）;
G=gf（::2）;
B=gf（::3）;
?
%測量圖像尺寸參數(shù)?
[st]=size（a）;?
?
%計算每個灰度級出現(xiàn)的概率?
h=imhist（a）/（s*t）;?
?
%計算灰度Xe，概率p?
for?k=0:255?
????if?sum（h（1:k+1））<（1/2）?
???????continue;?
????elseif?sum（h（1:k+1））==（1/2）?
????????Xe=k;?
????????break;?
????else?
????????Xe=k-1;?
????????break;?
????end?
end?
?????????
p=sum（h（gfMin+1:Xe））;?
%?計算DSIHE直方圖均衡化后的輸出Y（ij）?
for?i=1:s?
????for?j=1:t?
????????x=gf（ij）;?
????????if??x????????????cl=sum（h（gfMin+1:x+1））/p;?
????????????y（ij）=gfMin+（Xe-1-gfMin）*cl;?
????????else?
????????????cu=sum（h（Xe+1:x+1））/（1-p）;?
????????????y（ij）=Xe+（gfMax-Xe）