getconnectedsize.m

识别连通区域

function  ConnectedSize=getConnectedSize

%load a I;

I=imread('f:\2.bmp');

fnum = size(I,2);
snum = size(I,1);

direct_x = zeros(1,8);
direct_y = zeros(1,8);
direct_x = [0,1,-1,-1,-1,1,1,0];
direct_y = [-1,-1,1,0,-1,0,1,1];

count = 1;
for s = 2:snum-1
    for f = 2:fnum-1
       if(I(s,f)~=0)
           c = 0;
           for i = 1:8
                 if(I((s+direct_x(1,i)),(f+direct_y(1,i)))==0) 
                     c = c+1;   
                 end 
           end
           if(c==8)
               count = count+1;
           end
       end
    end
end  %统计孤立点的个数count。

BW2 = bwareaopen(I, 2);
[labeled,numObjects] = bwlabel(BW2,8);

ConnectedSize_Max = max(max(labeled));
Total_num = ConnectedSize_Max+count
ConnectedSize = zeros(1,Total_num);

for c = 1:ConnectedSize_Max
    num = 0;
    for f = 1:fnum
        for s = 1:snum
            if(labeled(s,f)==c)
                num = num+1;
            end
        end
    end
    ConnectedSize(1,c)=num;
end

for c = 1:count
     ConnectedSize(1,ConnectedSize_Max+c)=1;
end