用matlab计算二值图像的几何矩和中心矩

计算二值图像的几何矩和中心矩可以使用MATLAB中的`regionprops`函数。

示例代码如下：

% 读入二值图像
bw = imread('binary_image.png');

% 计算二值图像的几何矩和中心矩
stats = regionprops('table', bw, 'Centroid', 'Area', 'BoundingBox', 'MajorAxisLength', 'MinorAxisLength', 'Orientation');

% 输出中心矩
% 第一阶中心矩即为质心
x = stats.Centroid(1);
y = stats.Centroid(2);
fprintf('质心坐标为 (%.2f, %.2f)\n', x, y);

% 计算二阶中心矩
mu20 = stats.MajorAxisLength^2/12 + stats.BoundingBox(4)^2/4 - y^2;
mu02 = stats.MinorAxisLength^2/12 + stats.BoundingBox(3)^2/4 - x^2;
mu11 = (stats.BoundingBox(3)*stats.BoundingBox(4))/4 - x*y;
fprintf('二阶中心矩为:\n%.2f %.2f\n%.2f %.2f\n', mu20, mu11, mu11, mu02);

需要注意的是，`regionprops`函数计算的是连通区域的几何矩和中心矩，因此需要保证输入的二值图像中只有一个连通区域。如果有多个连通区域，可以使用`bwlabel`函数将其分割成多个单独的区域，然后分别计算每个区域的几何矩和中心矩。