直方图统计局部增强

本代码实现了图像局部直方图统计增强，并与全局直方图增强作对比。

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