图像处理：连通域质心求取

第一个例子：求取单连通域的质心
img=imread('bw.tif');
[m,n]=size(img);

L=bwlabel(img);
sum_x=0;sum_y=0;area=0;
for i=1:m
    for j=1:n
        if L(i,j)==1
            sum_x=sum_x+i;  %所有x*f(x)的和
            sum_y=sum_y+j;  %所有y*f(y)的和
            area=area+1;    %面积求和
        end
    end
end

plot_x=fix(sum_x/area);  
plot_y=fix(sum_y/area);
figure();
imshow(img);
hold on
plot(plot_y,plot_x,'*');
第2个例子：求取单连通域的质心

## 例2：多连通域的质心
I=imread('img.jpg'); 
I_gray=rgb2gray(I);
level=graythresh(I_gray);

[height,width]=size(I_gray);
I_bw=im2bw(I_gray,level);

for i=1:height %%循环中进行反色
for j=1:width
    if I_bw(i,j)==1
        I_bw(i,j)=0;
    else I_bw(i,j)=1;
    end
    
end
end

[L,num]=bwlabel(I_bw,8);  %num是连通域个数
plot_x=zeros(1,num);  %%用于记录质心位置的坐标
plot_y=zeros(1,num);

for k=1:num  %%num个区域依次统计质心位置
    sum_x=0;sum_y=0;area=0;
    for i=1:height
    for j=1:width
       if L(i,j)==k
        sum_x=sum_x+i;
        sum_y=sum_y+j;
        area=area+1;   
       end
    end
    end
    plot_x(k)=fix(sum_x/area);
    plot_y(k)=fix(sum_y/area);
end

figure(1);
imshow(I_bw);
for i=1:num
hold on
plot(plot_y(i) ,plot_x(i), '*');
end