我正在检测图像中的圆圈。我返回圆半径和轴的X， Y。我知道如何裁剪1个圆，公式没有问题：X-半径，Y-半径，宽度=2*r，高度=2*r使用im作物。我的问题是当我返回超过1个圆时。我在数组半径数组中得到返回的圆半径。我在centar中得到返回的圆心。当idisp（centar）时，它看起来像这样：

146.4930  144.4943
610.0317  142.1734

当我检查size（centarray）和disp时，我得到：

2 2

所以我知道第一列是X，第二列是Y轴值。所以第一个圆心是146,144。我做了一个只适用于1个圆的循环. "-------" 是我不确定用什么来得到的地方：注意：半径=r

第一圈）

            X = centarray(1)-r;
            Y = centarray(3)-r;
            Width =2*r;
            Width =2*r;

第二圈）

            X = centarray(2);
            Y = centarray(4);
            Width =2*r;
            Width =2*r;

我将如何修改我的代码的 "------" 部分？我也想，如果有3个圆圈，循环将工作，因为我有时从图像中获得多达9个圆圈。

B = imread('p5.tif');
centarray = [];
centarray = [centarray,centers];
radiiarray = [];
radiiarray = [radiiarray,radii];

for j=1:length(radiiarray)

x = centarray((------))-radiiarray(j); %X value to crop
y = centarray((------))-radiiarray(j);    %Y value to crop

width = 2*radiiarray(j);            %WIDTH
height = 2*radiiarray(j);           %HEIGHT

K = imcrop(B, [x y width height]);
end

我的完整代码，不起作用，当我看到值的存储方式时，我意识到为什么…：

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% DETECT + GET X Y WIDTH HEIGHT OF CIRCLES
I = imread('p5.tif');
subplot(2,2,1);imshow(I);title('Original Image'); 
%sharpen edges
B = imsharpen(I);
subplot(2,2,2);imshow(B);title('sharpened edges');


%find circles
Img = im2bw(B(:,:,3));
minRad = 20;
maxRad = 90;
[centers, radii] = imfindcircles(Img, [minRad maxRad], ...
                                    'ObjectPolarity','bright','sensitivity',0.84);
imagesc(Img);
viscircles(centers, radii,'Color','green');
%nuber of circles found


%arrays to store values for radii and centers
centarray = [];
centarray = [centarray,centers];
radiiarray = [];
radiiarray = [radiiarray,radii];
sc = size(centarray);
disp(sc)
disp(centarray)
disp(radiiarray)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%CROP USING VALUE FROM ARRAYS NUMBER OF TIMES THERE ARE CENTERS(number of
%circles)


for j=1:length(radiiarray)


    x = centarray((2*j)-1)-radiiarray(j); %X value to crop
    y = centarray((2*j))-radiiarray(j);    %Y value to crop 
    width = 2*radiiarray(j);            %WIDTH
    height = 2*radiiarray(j);           %HEIGHT
    disp(x)
    disp(y)
    disp(centarray)
    %crop using values
    K = imcrop(B, [x y width height]);
    %togray
    gray = rgb2gray(K);
    subplot(2,2,3);imshow(K);title('cropped before bw');
    Icorrected = imtophat(gray, strel('disk', 15));
    %to black and white
    black = im2bw(Icorrected);
    subplot(2,2,4);imshow(black);title('sharpened edges');
    %read
    results = ocr(black); 
    number = results.Text;
    %display value
    disp(number)

end

任何关于如何创建这种循环的帮助都是值得赞赏的，因为我只是没有更多的想法或找不到答案。

编辑解决方案

嗨，答案是将矩阵视为二维。

for j=1:length(radiiarray)
x=centarray(j,1)
y=centarray(j,2)
width = radiiarray(j)
height = radiiarray(j)
end

随着j增加值现在正确更新。