clc;
clear;
X=imread（‘D:\人臉識別論文\程序\程序\第四章（自商圖像）\自商圖像（可以用）\002.bmp‘）;
figure
subplot（221）;
imshow（X）

[mn]=size（X）;
allmean=mean（mean（X））;
sum1=0;sum2=0;
for?i=1:m
????for?j=1:n
????????if?X（ij）????????????sum1=sum1+X（ij）;
????????????sum2=sum2+1;
????????end
????end
end
threshold=sum2/sum1;
A=zeros（mn）;
for?i=1:m
????for?j=1:n
????????if?X（ij）>threshold
????????????A（ij）=1;
????????else
????????????A（ij）=0;
????????end
????end
end
subplot（222）;
imshow（A）
%對圖像進行二值化

%利用模板確定圖像中比較暗的部分
M=damb（mn）;
Afan=zeros（mn）;
for?i=1:m
????for?j=1:n
????????if?A（ij）==0?&&?M（ij）==1
????????????Afan（ij）=1;
????????else
????????????Afan（ij）=0;
????????end
????end
end
%斷開連接處理
se1=strel（‘square‘8）;
Afan=imerode（Afanse1）;
Afan=imdilate（Afanse1）;?
subplot（223）;
imshow（Afan）

[Afannum]=bwlabeln（Afan4）;
area1=0.01*m*n
area2=0.07*m*n
B=zeros（mn）;
for?i=1:num
????l=length（find（Afan==i））;
????if?l>=area1?&&?l<=area2
????????B（find（Afan==i））=1;
????end
end
%確定出需要調整的暗區

X（find（B==1））=50;
subplot（224）;
imshow（X）
imwrite（X‘D:\abcdefj.bmp‘‘bmp‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-11-12?10:33??SQI\
?????文件????????1214??2014-10-05?17:01??SQI\alonedark.m
?????文件????????1416??2014-10-05?17:20??SQI\alonelight.m
?????文件?????????320??2014-10-05?11:01??SQI\damb.m
?????文件?????????950??2014-10-23?21:45??SQI\iterative.m
?????文件????????2099??2014-10-19?15:47??SQI\SQI.m