求形心的这么写好像不可以啊,版主帮忙看一下啊
im=imread('b.bmp'); [m,n]=size(im); sumxc=0; sumyc=0; sum=0;
for i=1:m, for j=1:n
if pic(i,j)>0
sumxc=sumxc+pic(i,j)*i; sumyc=sumyc+pic(i,j)*j; sum=sum+pic(i,j);
end
end; end
xc=sumxc/sum; yc=sumyc/sum;
figure,imshow(im); hold on; plot(xc,yc,'r+');