clc;clear;
%導入圖片
srcDir=uigetdir（‘‘）;
allnames=struct2cell（dir（‘*.bmp‘））;?
[klen]=size（allnames）;
%批量處理
for?i=1:len
name=allnames{1i};
a=imread（name）;?
a1=a（::1）;
b=im2bw（a10.65）;
st=strel（‘disk‘8）;
c=imopen（bst）;
d=bwareaopen（c30000）;
d=~d;
e=medfilt2（d）;
f=bwareaopen（e50000）;
%擴充平移
if?i~=16&i~=17&i~=18&i~=48&i~=138
m=zeros（800200）;
f=[f?m];
else
???f0=zeros（200800）;
???f=[f;f0];
end
%去除葉柄
j=500;?k=sum（f（j:）==1）;
while?k>=80
????j=j+1;
????k=sum（f（j:）==1）;
???if??j==800
???????j=680;k=79;
??