f=imread('12.jpg');
imshow(f);
f=im2double(f);
h=fspecial('gaussian',25,15)
g=imfilter(f,h,'replicate');
figure,imshow(g);%平滑處理
g=im2bw(g,1.5*graythresh(g));
figure,imshow(g);


s=bwmorph(g,'skel',Inf);
figure,imshow(s);
s1=bwmorph(s,'spur',8);
figure,imshow(s1);