叶片轮廓提取代码

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;