matlab分割大圆和小圆

结果：

 过程：

 

 

 

 代码：

close all; clear all; clc; 
 
I = imread('Fig0943.tif'); 
se = strel('disk', 26); 
se1 = strel('disk', 50); 
se2 = strel('disk', 30);
se3 = strel('disk', 35);
K = imclose(I, se); 
figure; 
subplot(121), imshow(I); title('原图');
subplot(122), imshow(K); title('闭运算去小圆');
 
 
N = imdilate(K, se2);%膨胀去除边缘黑色
L=imerode(N,se3);
figure; 
subplot(121), imshow(N); title('膨胀去除边缘黑色');
subplot(122), imshow(L);title('腐蚀');
 
 
J = imopen(L, se1); 
J= imbinarize(J);
M = edge(J, 'log'); 
figure; 
subplot(121), imshow(J);title('开运算并且二值化'); 
subplot(122),imshow(M); title('提取边界');
 
M = im2uint8(M);%将图像格式转换为unit8
 
K=imadd(I,M);
figure; 
imshow(K); title('相加后结果图');