X=imread（‘I.bmp‘）;
%X=rgb2gray（X）;
figureimshow（X）;
title（‘原始圖像‘）;

Y=histeq（X）;

figureimshow（Y）;
title（‘全局均衡化圖像‘）;
E=6.0;k0=0.3;k1=0.2;k2=0.3;
[mn]=size（X）;
X1=X;
Mg=sum（X1（:））/（m*n）;
Vg=sum（sum（（X1-Mg）.^2））/（m*n）;
for?i=2:m-1
????for?j=2:n-1
????????X2=X1（i-1:i+1j-1:j+1）;
????????Ms=sum（X2（:））/9;
????????Vs=sum（sum（（X2-Ms）.^2））/9;
????????if（Ms<=k0*Mg&&k1*Vg<=Vs&&Vs????????????X1（ij）=E*X1（ij）;
????????end
????end
end
imwrite（X1‘x.bmp‘）;
figureimshow（X1）;
title（‘直方圖統計局部增強圖像‘）;

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

?????文件????????547??2018-05-10?10:35??histeqarea.m

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

??????????????????547????????????????????1