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