clear
close?all
I=imread（‘D:\畢設\測試圖片3\test0066.jpg‘）;?
subplot（331）imshow（I）?title（‘原始圖像‘）;
I1=rgb2gray（I）;
subplot（332）imshow（I1）title（‘灰度圖像‘）;?
Q1=medfilt2（I1[77]）;?
subplot（333）imshow（Q1）;title（‘去噪圖像‘）;??
Z1=histeq（Q1）;?subplot（334）;
imhist（Q1）;title（‘直方圖‘）;?
subplot（335）;
imhist（Z1）;
title（‘均衡化‘）;
temp=60;?
[mn]=size（Q1）;?
Z2=zeros（mn）;?
for?i=1:m????
????for?j=1:n??????????
????????if?Q1（ij）>30&Q1（ij）<140;??????????????
????????????Z2（ij）=100;????
????????else??Z2（ij）=0;????????????
????????end
????end
end
subplot（336）;?
imshow（Z2）;
title（‘閾值點‘）;
Z1=double（Z1）;?
f=fft2（Z1）;?
TO1=fftshift（f）;?
[mn]=size（f）;?
d0=10;?
r1=0.5;?
rh=2;?
c=4;?
n1=floor（m/2）;?
n2=floor（n/2）;?
for?i=1:m;?????
????for?j=1:n?????????
????????d=sqrt（（i-n1）^2+（j-n2）^2）;?????????
????????h=（rh-r1）*（1-exp（-c*（d.^2/d0.^2）））+r1;????????
????????TO1（ij）=h*TO1（ij）;?????
????end
end
TO1=ifftshift（TO1）;
TO1=uint8（real（ifft2（TO1）））;?
subplot（337）;
imshow（TO1）;
title（‘同態濾波圖像‘）;
E1=im2bw（TO10.70）;?
subplot（338）;
imshow（E1）;
title（‘二值化處理圖像‘）;
BY1=edge（TO1‘canny‘）;
subplot（339）imshow（BY1）title（‘邊緣檢測后圖像‘）;
se=[1;1;1];??
FS1=imerode（BY1se）;?
figuresubplot（331）imshow（FS1）title（‘腐蝕后邊緣圖像‘）;?
se=strel（‘rectangle‘[2525]）;
TC1=imclose（FS1se）;
subplot（332）imshow（TC1）title（‘填充后圖像‘）;
XT1=bwareaopen（TC12000）;
subplot（333）imshow（XT1）t