Question

我正在检测图像中的圆圈。我返回圆的半径和轴的X，Y。我知道如何裁剪1圈没有问题的公式： X-radius, Y-radius, width=2*r,height=2*r使用imcrop。我的问题是当我返回超过1圈。我在数组radiiarray中返回圆圈半径。我在centarray返回了圈子中心。当我disp(centarray)时，它看起来像这样：

146.4930  144.4943
610.0317  142.1734

当我检查size(centarray)和disp时，我得到了：

2 2

所以我理解第一列是X，第二列是Y轴值。所以第一圈中心将是146,144。我制作了一个仅适用于1圈的循环。 “-------”是我不确定要使用的地方：注意：radius = r

第1圈）

            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现在正确地增加值更新。

Answer 1

答案是将矩阵视为二维。

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

因为j现在正确地增加值更新。

感谢@beaker的评论！这就是为什么我想出来的

循环遍历数组矩阵matlab

1 个答案: