MATLAB中灰度图像的圆检测

Question

我的图像如下所示.我的目标是检测第二张图像中显示的圆圈.我用过[centers,radii] = imfindcircles(IM,[100 300]);但却一无所获.

有没有其他方法来检测圆圈？我怎样才能做到这一点？

原始图片:

圆圈:我用油漆画了它.

Answer 1

这是imfindcircles的替代解决方案.基本上对图像进行阈值处理,使用磁盘结构元素对其进行扩展,然后在找到边缘后,使用此处circle_hough的文件交换可用的算法应用Hough变换来检测圆.

这是代码:

clear
clc
close all

A = imread('CircleIm.jpg');

%// Some pre-processing. Treshold image and dilate it.
B = im2bw(A,.85);

se = strel('disk',2);

C = imdilate(B,se);

D = bwareaopen(C,10000);

%// Here imfill is not necessary but you might find it useful in other situations.
E = imfill(D,'holes');

%// Detect edges
F = edge(E);

%// circle_hough from the File Exchange.

%// This code is based on Andrey's answer here:
%https://dsp.stackexchange.com/questions/5930/find-circle-in-noisy-data.

%// Generate range of radii.
 radii = 200:10:250;

h = circle_hough(F, radii,'same');
[~,maxIndex] = max(h(:));
[i,j,k] = ind2sub(size(h), maxIndex);
radius = radii(k);
center.x = j;
center.y = i;

%// Generate circle to overlay
N = 200;

theta=linspace(0,2*pi,N);
rho=ones(1,N)*radius;

%Cartesian coordinates
[X,Y] = pol2cart(theta,rho); 

figure;

subplot(2,2,1)
imshow(B);
title('Thresholded image  (B)','FontSize',16)

subplot(2,2,2)
imshow(E);
title('Filled image (E)','FontSize',16)

subplot(2,2,3)
imshow(F);hold on

plot(center.x-X,center.y-Y,'r-','linewidth',2);

title('Edge image + circle (F)','FontSize',16)

subplot(2,2,4)
imshow(A);hold on
plot(center.x-X,center.y-Y,'r-','linewidth',2);
title('Original image + circle (A)','FontSize',16)

这给出了以下内容:

您可以使用传递给阈值的参数或扩展参数来查看它如何影响结果.

希望有所帮助!