写一段matlab代码，实现圆心提取

以下是一段MATLAB代码，用于提取圆心：

% 读取图像
img = imread('image.jpg');

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

% 边缘检测
edge_img = edge(gray_img, 'canny');

% Hough变换
[centers, radii] = imfindcircles(edge_img, [20 50]);

% 绘制圆心
imshow(img);
viscircles(centers, radii, 'EdgeColor', 'b');

% 输出圆心坐标
disp(centers);