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