数组的质心,可以看为将数组值作为像素点灰度值时的重心。
在文章Matlab质心算法中提到了质心的算法。
第一种算法源自文章Matlab质心算法:
% B为待求质心的数组
x = 0;
y = 0;
sum1 = 0;
for i = 1 : size(B, 1)
for j = 1 : size(B, 2)
y = y + B(i, j) * i;
x = x + B(i, j) * j;
sum1 = sum1 + B(i, j);
end
end
x = x / sum1;
y = y / sum1;
文章Matlab质心算法中也给出了另一种算法:求每一行/列的质心,然后再求数组的质心。
但在本文中有一个更简单的算法,代码更简洁:
% B为待求质心的数组
x1 = sum(sum(B .* [1 : size(B, 2)])) / sum(B(:));
y1 = sum(sum(B' .* [1 : size(B, 1)])) / sum(B(:));
欢迎大家批评指正。