【MATLAB】一个图像处理的小作业

马上就要读研究生了，希望自己养成写博客的好习惯，能及时整理学到的东西，同时也共享给大家。

今天做了一个图像处理的小作业，题目如下：

显微应用中一个预处理步骤是从两组或更多组重叠的类似颗粒（见下图）中分离出单个独立的一种圆颗粒。假设所有颗粒的大小相同，提出一种产生3幅图像的形态学算法，这3幅图像分别仅由如下物体组成：
（a）仅与图像边界融合在一起的颗粒
（b）仅彼此重叠的颗粒
（c）没有重叠的颗粒

分析

       题目要求分开三类颗粒：与边界重叠的、彼此重叠的（两个颗粒）以及没有重叠的单个颗粒。显然可以采取面积法。

解答

clc;clear;

%读入图像
f=imread('FigP0936.tif');
[a,b]=size(f);

%二值化
f=im2bw(f);

%开运算
se=strel('disk',3');%圆盘型结构元素
f=imopen(f,se);%开运算

subplot(141);imshow(f);
title('二值化图像','fontname','Microsoft YaHei UI Light','FontSize',16);

% L为标记矩阵，n为找到连接分量的总数
[L,n]=bwlabel(f);


%仅与边界重叠的颗粒
f1=f;
for k=1:a
    if(L(k,1)~=0)
        I= L==L(k,1);
        f1(I)=0;
        I=0;
    end
end
for k=1:a
    if(L(k,b)~=0)
        I= L==L(k,b);
        f1(I)=0;
        I=0;
    end
end
for k=1:b
    if(L(1,k)~=0)
        [I]=find(L==L(1,k));
        f1(I)=0;
        I=0;
    end
end
for k=1:b
    if(L(a,k)~=0)
        [I]=find(L==L(a,k));
        f1(I)=0;
        I=0;
    end
end
f1=f-f1;
subplot(142);imshow(f1);

title('仅与边界重叠的','fontname','Microsoft YaHei UI Light','FontSize',16);

%仅与彼此重叠的颗粒，面积在300~450间
f2=f;
for k=1:n
    [r,c]=find(L==k);
    if( size(r,1)>300 && size(r,1)<450 )
        for i=1:size(r,1)
            f2(r(i),c(i))=0;
        end
    end 
end
f2=f-f2;

subplot(143);imshow(f2);
title('仅与彼此重叠的','fontname','Microsoft YaHei UI Light','FontSize',16);

%没有重叠的颗粒
f3=f;
for k=1:n
    [r,c]=find(L==k);
    if( size(r,1)<300 )
        for i=1:size(r,1)
            f3(r(i),c(i))=0;
        end
    end 
end
f3=f-f3;
f3=f3-f1;%排除与边界重叠的颗粒

subplot(144);imshow(f3);
title('没有重叠的','fontname','Microsoft YaHei UI Light','FontSize',16);

最后得到的结果如下