结合全局-局部图像对比度增强的程序（Global Local Contrast Enhancement）

结合全局-局部图像对比度增强的程序。首先使用全局对比度增强，然后使用局部对比度增强。程序中的参数是经过导师修改过的。文件给出了原论文。对图像有一定的增强效果。

x=imread（‘D:\longlongzhang\tuxiangchuli\jianzhu.jpg‘）;
f1=rgb2gray（x）;
[MN]=size（f1）;
figure（1）
imshow（f1[]）
title（‘原灰度图像‘）
f2=double（f1）/255;
L=3;
L1=（L-1）/2;
Cg=0.07;
C=0.2;
f=zeros（MN）;
gmean=0.0;
for i=1:M
    for j=1:N
        gmean=gmean+f2（ij）;
    end
end
gmean=gmean/（M*N）;
for i=1:M
    for j=1:N
        fx（ij）=（1+Cg）*（f2（ij）-gmean）+0.5;
    end
end
for i=L1+1:M-L1
    for j=L1+1:N-L1
        m=0.0;
        for i1=-L1:L1
            for j1=-L1:L1
                m=m+fx（i+i1j+j1）;
            end
        end
        m=m/（L*L）;
        sigma2=0.0;
        for i1=-L1:L1
            for j1=-L1:L1
                sigma2=sigma2+（fx（i+i1j+j1）-m）^2;
            end
        end
        sigma2=sigma2/（L*L）;
        lsd=sigma2^0.5;
        f（ij）=fx（ij）+C/（lsd+eps）*（fx（ij）-m）;
    end
end
min1=min（f（:））;
max1=max（f（:））;
fx=round（（f-min1）/（max1-min1）*255）;
% f3=mat2gray（f）;
% f4=im2uint8（f3）;
figure（2）
imshow（uint8（fx）[]）
title（‘增强后的灰度图像‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1029  2014-11-11 21:12  Global-Local Contrast Enhancement\main.m

     文件     688076  2014-09-23 16:42  Global-Local Contrast Enhancement\pxc3882365（局部增强）.pdf

     目录          0  2014-11-11 21:28  Global-Local Contrast Enhancement

----------- ---------  ---------- -----  ----

               689105                    3