matlab中使用霍夫变换及最小二乘法自动提取不规则圆周上的数据点的代码？

以下是使用霍夫变换和最小二乘法自动提取不规则圆周上的数据点的MATLAB代码：

% 读入图像
img = imread('irregular_circle.png');

% 转为灰度图像
gray = rgb2gray(img);

% 阈值分割
thresh = graythresh(gray);
binary = imbinarize(gray, thresh);

% 边缘检测
edgeImg = edge(binary, 'Canny');

% 霍夫变换
[centers, radii, ~] = imfindcircles(edgeImg, [10 100], ...
    'Sensitivity', 0.95, 'ObjectPolarity', 'bright');

% 显示圆心和半径
figure;
imshow(img);
viscircles(centers, radii, 'EdgeColor', 'b');

% 提取不规则圆周上的数据点
theta = linspace(0, 2 * pi, 100);
x = zeros(length(theta), 1);
y = zeros(length(theta), 1);
for i = 1:length(theta)
    % 计算圆周上的点坐标
    x(i) = centers(1) + radii(1) * cos(theta(i));
    y(i) = centers(2) + radii(1) * sin(theta(i));
end

% 最小二乘法拟合圆周
A = [x y ones(length(theta), 1)];
b = -x.^2 - y.^2;
c = (A' * A) \ (A' * b);
xc = -0.5 * c(1);
yc = -0.5 * c(2);
r = sqrt((c(1)^2 + c(2)^2) / 4 - c(3));

% 显示拟合的圆周和数据点
hold on;
plot(x, y, 'rx');
viscircles([xc yc], r, 'LineStyle', '--');
``