每天一点matlab——提取图像中的小圆圈

%图像变为bmp文件
clear all;clc;
I=imread('I.bmp'); %读取图片
imshow(I);title('源图像');
J=im2bw(I,graythresh(I));
figure,imshow(J);title('初始二值图像');
J=~J;
figure,imshow(J);title('取反目标白色');
%除去边界白色非目标
J=imclearborder(J,4);%imclearborder:Suppress light structures connected to image border,抑制和图像边界相连的亮对象
figure;imshow(J),title('目标图像');
%腐蚀膨胀操作，去除小的白色点
se=strel('disk',1);
J=imerode(J,se);
J=imdilate(J,se);
figure,imshow(J);title('最终图像');