条纹二值化，细化，轮廓提取源码

[filename pathname] = uigetfile（‘*.bmp;*.pgm;*.tif;*.jpg;*.*‘ ‘Load input image‘）;
 filename = [pathname filename]; 
 im_original=imread（filename）;
 [nxnyl]=size（im_original）; 
 if l>1
     image_I=rgb2gray（im_original）;%将原图象转换为灰度图象
 else
     image_I= im_original;
 end
 figure（1）;imshow（image_I）;title（‘原图像‘）
 



 image_I=double（image_I）;
 image_I12=image_I;
  for i=1:3
 image_I1=imfilter（image_Ifspecial（‘gaussian‘102））;
  end
 figure（3）;imshow（uint8（image_I1））;title（‘lvbo图像‘） 
 

im_graylap=image_I1;
T1=0.5*（double（min（im_graylap（:）））+double（max（im_graylap（:））））;
done=0