第一个例子:求取单连通域的质心
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