www.gusucode.com > 关于海航matlab和lingo的训练题 > center.m

    %3.求质心
function [xmean, ymean] = center(bw)
bw=bwfill(bw,'holes');
area = bwarea(bw);


 [m n] =size(bw);
bw=double(bw);
xmean =0; ymean = 0;
for i=1:m,
for j=1:n,
xmean = xmean + j*bw(i,j);
ymean = ymean + i*bw(i,j);
end;
end;
if(area==0)
    xmean=0;
    ymean=0;
else 
    xmean = xmean/area;
    ymean = ymean/area;
    xmean = round(xmean);
  ymean = round(ymean);
    
end