clc;
clear?all;
image=imread（‘5.bmp‘）;
im1=rgb2gray（image）;
figure（10）imshow（im1）;

im1=medfilt2（im1[3?3]）;?%二維中值濾波處理圖像噪聲%
figure（12）imshow（im1）;

BW?=?edge（im1‘sobel‘）;?%邊緣檢測
figure（13）imshow（BW）;

[imximy]=size（BW）;
msk=[0?0?0?0?0;
?????0?1?1?1?0;
?????0?1?1?1?0;
?????0?1?1?1?0;
?????0?0?0?0?0;];
B=conv2（double（BW）double（msk））;?%二維卷積操作??Smoothing??image?to?reduce?the?number?of?connected?components
figure（1）imshow（B）;

grd=B;
bg1=imclose（grdstrel（‘rectangle‘[519]））;%取矩形框的閉運算
figure（2）imshow（bg1）;title（‘圖像閉運算[519]‘）;%輸出閉運算的圖像
bg3=imopen（bg1strel（‘rectangle‘[519]））;%取矩形框的開運算
figure（3）imshow（bg3）;title（‘圖像開運算[519]‘）;%輸出開運算的圖像
bg2=imopen（bg3strel（‘rectangle‘[191]））;%取矩形框的開運算
figure（4）imshow（bg2）;title（‘圖像開運算[191]‘）;%輸出開運算的圖像

[xy]=size（bg2）;???????????????
II=double（bg2）;?


%將I變為對應區塊的函數?區域照度優化！
x11=x/2;
x12=uint16（x11）;
x13=double（x12）;
y11=y/4;
y12=uint16（y11）;
y13=double（y12）;
y23=3*y13;
for?i=x13:x
????for?j=y13:y23
????????I（ij）=II（ij）;
????end
end

figure（14）imshow（II）;
figure（15）imshow（I）;
%取邊緣????二值化?????????????????
z0=max（max（I））;??%取I中最大值??????????????
z1=min（min（I））;??%取I中最小值????????????????
T=（z0+z1）/2;??????????????????????
TT=0;
S0=0;?n0=0;
S1=0;?n1=0;
allow=0.5;???????????????????????
d=abs（T-TT）;
count=0;????????????????????????

[aabb]=size（I）
while（d>=allow）?????????????????
????count=count+1;
????for?i=1:aa
????????for?j=1:bb
????????????if?（I（ij）>=T）
????????????????S0=S0+I（ij）;
????????????????n0=n0+1;
????????????end
????????????if?（I（ij）????????????????S1=S1+I（ij）;
????????????????n1=n1+1;
????????????end
????????end
????end?
????T0=S0/n0;
????T1=S1/n1;
????TT=（T0+T1）/2;
????d=abs（T-TT）;
????T=TT;
end

Seg=zeros（xy）;

for?i=1:aa
????for?j=1:bb
????????if（I（ij）>=1.9*T）%照度影響明顯?降低必然明顯影響準確度
????????????Seg（ij）=1;???????????????
????????end
????end
end
figure（16）imshow（Seg）;
SI=1-Seg;????????
figure（17）imshow（SI）;
se1=strel（‘square‘3）;?%選擇3腐蝕?為使得邊界精細?????????????
SI1=imerode（SIse1）;??
figure（17）imshow（SI1）;????
BW=SI-SI1;??????????????????????????
figure（5）;
imshow（BW）;title（‘BW‘）?????


%?增加顆粒
hsize=[33];
H?=?fspecial（‘average‘hsize）;?%預定義濾波器
smoothed?=?imfilter（BWH‘replicate‘）;
figure（6）;
imshow（smoothed）;?title（‘3*3?Smoothed?Image‘）;
%填充封閉區域
BW1?=?bwfill（BW‘holes‘）;
figure（7）;
imshow（BW1）;
%取邊緣
BW2?=?bwmorph（BW1‘remove‘）;
figure（8）;
imshow（BW2）;
%close?pot?去除較小的白點
tt2=bwarea（BW1）;?%計算BW的面積
tt3=tt2/2
if（tt3>=10000）
???????J2=bwareaopen（BW110000）;
?????elseif（tt3>=9000）?
?????????J2?=?bwareaopen（BW19000）;
?????????elseif（tt3>=8000）?
?????????J2?=?bwareaopen（BW18000）;
?????elseif（tt3>=7000）?
?????????J2?=?bwareaopen（BW17000）;
?????elseif（tt3>=6000）?
?????????J2?=?bwareaopen（BW16000）;
?????elseif（tt3>=5000）?
?????????J2?=?bwareaopen（BW15000）;
?????elseif（tt3>=4000）?
?????????J2?=?bwareaopen（BW14000）;
?????elseif（tt3>=3000）
?????????J2?=?bwarea

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

?????文件???????4545??2011-12-24?15:16??r1.m

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

?????????????????4545????????????????????1