灰度图像处理——连通分量提取（Matlab）

clc;
clear all;
close all;
%提取连通分量测试图像，用来检测食物中的外来物
I=im2double(imread('D:\Gray Files\9-18a.tif'));
%获得图像大小
[M,N]=size(I);
%存放腐蚀后的图像
%设定提取门限，根据经验设定
th=0.808;
%将原图像根据上述门限进行二值处理，大于等于门限为1，小于为0
ind=find(I>=th);
I(ind)=1;
ind=find(I<th);
I(ind)=0;
%获得值为1的x,y坐标
[rows,cols]=find(I==1);
%将x,y坐标点组合成下标矩阵
index=cat(2,rows,cols);
%5*5结构元素
n=5;
B=ones(n,n);
n_l=floor(n/2);
%对边界图进行扩充，四周各加2行、2列0（与结构元素的大小相对应），目的是为了处理边界点
I_pad=padarray(I,[n_l,n_l]);
%将下标矩阵加2
index=index+2;
%保存所有连通对象集
Objs={};
%设定递归上限
set(0,'RecursionLimit',10000);
%递归寻找连通对象
while 1
    %取出下标矩阵中的第一个值
    if isempty(index)
        break;
    end
    %取出下标矩阵中的第一个对象
    p=index(1,:);
    A=FindConnectedPoint(p,I_pad,index,n_l);  
    %将边界线对象保存至对象集
    Objs=cat(1,Objs,A{1,1}-n_l);
    index=A{1,2}; 
    I_pad=A{1,3};        
end
%显示连通对象
J=zeros(M,N);
if ~isempty(Objs)
    for i=1:size(Objs,1)
        Obj=Objs{i,1};
        for j=1:size(Obj,1)
            J(Obj(j,1),Obj(j,2))=1;
        end
    end
end
imshow(J)

FindConnectedPoint函数如下：

function R=FindConnectedPoint(p,Beta_pad,ind,n_l)
    %保存边界坐标对象
    Obj=[];
    %以p为中心点，从边界图中取出3*3大小的块区域
    Block=Beta_pad(p(1,1)-n_l:p(1,1)+n_l,p(1,2)-n_l:p(1,2)+n_l);
    Obj=p;
    %将图像中的p点置零
    Beta_pad(p(1,1),p(1,2))=0;
    %将块区域中的p点置零
    Block(n_l+1,n_l+1)=0;
    A=cell(1,3);
    %将p点从下标矩阵中删除
    for j=1:size(ind,1)
        if ind(j,:)==p
            ind(j,:)=[];
            break;
        end
    end
    %寻找块区域中为1的点的下标
    [rows,cols]=find(Block==1);
    ind_sub=cat(2,rows,cols);
    if ~isempty(ind_sub)
        %确定块中数值为1的点的坐标
        for i=1:size(ind_sub,1) 
            p_next=[];
            if ind_sub(i,1)<=n_l+1
                p_next(1,1)=p(1,1)-abs(n_l+1-ind_sub(i,1));
                if ind_sub(i,2)<=n_l+1                    
                    p_next(1,2)=p(1,2)-abs(n_l+1-ind_sub(i,2));
                else
                    p_next(1,2)=p(1,2)+abs(n_l+1-ind_sub(i,2));                    
                end
            else
                p_next(1,1)=p(1,1)+abs(n_l+1-ind_sub(i,1));
                if ind_sub(i,2)<=n_l+1                    
                    p_next(1,2)=p(1,2)-abs(n_l+1-ind_sub(i,2));
                else
                    p_next(1,2)=p(1,2)+abs(n_l+1-ind_sub(i,2));                    
                end                
            end
            if Beta_pad(p_next(1,1),p_next(1,2))~=0                        
                A=FindConnectedPoint(p_next,Beta_pad,ind,n_l);
                Obj=cat(1,A{1,1},Obj);   
            end                       
            %更新下标矩阵
            ind=A{1,2};
            Beta_pad=A{1,3};
        end
    end
    %返回结果
    R={Obj,ind,Beta_pad};
end