[qiyuewei]

a=im2double(a);
[M,N]=size(a);
M00=sum(a();
m10=0;
m01=0;

for i=1:M
for j=1:N
m10=m10+(i*a(i,j));
m01=m01+(j*a(i,j));
end

end
x=m10/M00;
y=m01/M00;
c=[x,y];