求连通域面积matlab

 matlab有自带函数bwreopen()求连通域，并移除小面积连通域，返回一个mask模板不同连通域用不同的label标志，背景用0标志。

function [ mask ] = bwareopenLZ( Lig,minarea)
%UNTITLED3 此处显示有关此函数的摘要
%   此处显示详细说明
row=size(Lig,1);
col=size(Lig,2);

mask=zeros(row,col);     %标记图像
label=1;
queue_head=1;       %队列头
queue_tail=1;       %队列尾
neighbour=[-1 -1;-1 0;-1 1;0 -1;0 1;1 -1;1 0;1 1];  %和当前像素坐标相加得到八个邻域坐标
q2=zeros(row*col,2);%存放坐标位置
for i=2:row-1
    for j=2:col-1       
        if Lig(i,j)==1 && mask(i,j) ==0           
            mask(i,j)=label;
            q2(queue_tail,:)=[i j];
            %q{queue_tail}=[i j];        %用元组模拟队列，当前坐标入列
            queue_tail=queue_tail+1;            
            while queue_head~=queue_tail
                %pix=q{queue_head};
                pix2=q2(queue_head,:);  %pix2是坐标:（x,y）
                for k=1:8               %8邻域搜索
                    %pix1=pix+neighbour(k,:);
                    pix1=pix2+neighbour(k,:);  %邻域坐标
                    if pix1(1)>=2 && pix1(1)<=row-1 && pix1(2) >=2 &&pix1(2)<=col-1
                        if Lig(pix1(1),pix1(2)) == 1 && mask(pix1(1),pix1(2)) ==0  %如果当前像素邻域像素为1并且标记图像的这个邻域像素没有被标记，那么标记
                          mask(pix1(1),pix1(2))=label;
                          %q{queue_tail}=[pix1(1) pix1(2)];
                          q2(queue_tail,:)=[pix1(1) pix1(2)];
                          queue_tail=queue_tail+1;
                        end
                    end
                end                
                queue_head=queue_head+1;  %从头遍历到队列尾部，即遍历8邻域像素的8邻域像素
            end     %此时结束一个连通域的搜索
            %统计该连通域的像素数，小于阈值的连通域进行去除
            area=queue_head-1;
            if(area<minarea)%去除面积小于阈值的连通域
                for i=1:area  %area即该连通域的像素个数 
                    x=q2(i,1);
                    y=q2(i,2);
                    mask(x,y)=0;
                end
            end                
            label=label+1;
            queue_head=1;
            queue_tail=1;            
        end    
    end
end
end