Center of mass and total mass of a matrix
tot_mass = sum(A(:));
[ii,jj] = ndgrid(1:size(A,1),1:size(A,2));
R = sum(ii(:).*A(:))/tot_mass;
C = sum(jj(:).*A(:))/tot_mass;
out = [tot_mass,R,C]
OR
x = 1 : size(A, 2); % Columns.
y = 1 : size(A, 1); % Rows.
[X, Y] = meshgrid(x, y)
meanA = mean(A(:))
centerOfMassX = mean(A(:) .* X(:)) / meanA
centerOfMassY = mean(A(:) .* Y(:)) / meanA
% Method 2: using regionprops
props = regionprops(true(size(A)), A, 'WeightedCentroid')