matlab最小二乘圆法求圆度

以下是使用MATLAB中的最小二乘圆法求圆度的步骤：

1.读取二值图像并显示：

img = imread('image.jpg');
imshow(img);

2.使用边缘检测算法检测图像的边缘坐标：

BW = edge(img,'canny');
imshow(BW);
[x] = find(BW);

3.使用最小二乘圆法拟合圆：

x_m = mean(x);
y_m = mean(y);
u = x - x_m;
v = y - y_m;
Suu = sum(u.^2);
Suv = sum(u.*v);
Svv = sum(v.^2);
Suuv = sum(u.*u.*v) ;
Suvv = sum(u.*v.*v);
A = [Suu Suv; Suv Svv];
B = [-Suuv/2 - Suvv/2]';
a = A\B;
xc = a(1);
yc = a(2);
R = sqrt((u-xc).^2 + (v-yc).^2);

4.计算目标像素面积及其圆度：

area = sum(BW(:));
perimeter = sum(R(:));
circularity = 4*pi*area/perimeter^2;

其中，`circularity`即为所求的圆度。