clc
clear?all
%讀入圖片
gf=imread（‘lena.tif‘）;
figureimhist（gf）
title（‘原始直方圖‘）;
gfMin=min（gf（:））;
gfMax=max（gf（:））;

%測量圖像尺寸參數
[st]=size（gf）;

%計算每個灰度級出現的次數
h1=imhist（gf）;

%計算圖像的平均灰度gray_mean
gray_mean=round（sum（gf（:））/（s*t））

????
%計算nlnu
nl=sum（h1（gfMin+1:gray_mean+1））;
nu=sum（h1（gray_mean+2:gfMax+1））;

%?計算BBH直方圖均衡化后的輸出Y（ij）
for?i=1:s
????for?j=1:t
????????x=gf（ij）;
????????if?x<=gray_mean
????????????cl=sum（h1（gfMin+1:x+1））/nl;
????????????y（ij）=gfMin+（gray_mean-gfMin）*cl;
????????else
????????????cu=sum（h1（gray_mean+2:x+1））/nu;
????????????y（ij）=（gray_mean+1）+（gfMax-（gray_mean+1））*cu;
????????end??
????end
end
y=mat2gray（y）;
%顯示BBH均衡化后的圖像及其直方圖
figureimshow（y）
title（‘BBH均衡化后的圖像‘）;
figureimhist（y）
title（‘BBH均衡化后的直方圖‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????863??2013-09-30?09:05??直方圖均衡化\BBHE_1997.m

?????文件????????921??2013-09-24?21:24??直方圖均衡化\DSIHE_1999.m

?????文件???????1245??2013-09-30?09:09??直方圖均衡化\MMBEBHE.m

?????目錄??????????0??2014-04-24?21:34??直方圖均衡化

-----------?---------??----------?-----??----

?????????????????3029????????????????????4