写一段matlab代码，不要用imfindcircles函数，找到图片中的所有圆形，并计算圆心和半径

### 回答1：

% 设定图片
I = imread('circle.jpg');% 将图片转换成灰度图
Igray = rgb2gray(I);% 对图片进行滤波，以降噪
Iblur = imgaussfilt(Igray,3);% 找到图片中的边缘
Iedge = edge(Iblur);% 将边缘变成极坐标
[r,c] = size(Iedge);
[theta,rho] = cart2pol(1:c,1:r);% 计算极坐标中的梯度
[theta_x,theta_y] = gradient(theta);
[rho_x,rho_y] = gradient(rho);% 将梯度相乘，找到圆心
centers = Iedge.*(theta_x.*rho_y-theta_y.*rho_x);% 找到圆心坐标
[xcenters,ycenters] = find(centers>0);% 计算圆心和半径
for i = 1:length(xcenters)
x = xcenters(i);
y = ycenters(i);
radius = mean(sqrt((x-xcenters).^2+(y-ycenters).^2));
disp(['圆心坐标:(' num2str(x) ',' num2str(y) '), 半径: ' num2str(radius)]);
end

### 回答2：
要找到图片中的所有圆形并计算圆心和半径，可以使用霍夫变换来实现。以下是MATLAB代码示例：

% 读取图片
img = imread('image.jpg');
% 将图片转换为灰度图像
grayImg = rgb2gray(img);
% 对灰度图像进行高斯滤波
filteredImg = imgaussfilt(grayImg, 5);

% 设置霍夫圆变换的参数
sensitivity = 0.85;
radiusRange = [10, 100];

% 霍夫圆变换
[centers, radii, metric] = imfindcircles(filteredImg, radiusRange, 'ObjectPolarity', 'dark', 'Sensitivity', sensitivity);
    
% 展示所有找到的圆形
figure;
imshow(img);
hold on;
viscircles(centers, radii,'EdgeColor','r');

% 显示圆心和半径
for i = 1:length(radii)
    text(centers(i, 1), centers(i, 2), sprintf('(%d, %d)\nRadius: %d', centers(i, 1), centers(i, 2), radii(i)), 'Color', 'red');
end

在这个示例中，首先读取图片并将其转换为灰度图像。然后对灰度图像进行高斯滤波，以平滑图像并减少噪声。接下来使用`imfindcircles`函数进行霍夫圆变换，找到圆形的圆心和半径。设置`Sensitivity`参数可以调整对圆的检测的敏感度，设置`radiusRange`参数限制要检测的圆的半径范围。最后，使用`viscircles`函数将找到的圆形可视化，并使用`text`函数在图像上显示圆心坐标和半径。