% %图像形态学处理
SE=strel('arbitrary',[0 1 0;1 1 1;0 1 0]);
imagen=imerode(imagen,SE);
imagen=imdilate(imagen,SE);
%imagen =~im2bw(imagen,threshold); %将图像二值化
%**************************************************
% Remove all object containing fewer than 30 pixels
imagen = bwareaopen(imagen,30); %排除小于30像素的连通区域
subplot(3,2,2);
imshow(imagen);
title('二值化结果');