x مرکز رو با میانگین x های همه نقاط شکل، و
y مرکز رو با میانگین y های همه نقاط شکل حساب کن.
function [yc, xc] = centermass(image)

[y, x] = find(image == 1);

xc = sum(x) / size(x, 1);
yc = sum(y) / size(y, 1);

end