二进制图像圆拟合

时间:2017-12-30 23:13:56

标签: matlab image-processing geometry

我是MATLAB的新手,我正在尝试将圆拟合技术应用于下图。

图像

enter image description here

预期输出

enter image description here

我使用以下代码但没有运气:

[centers, radii] = imfindcircles(W,[30 65],'ObjectPolarity','bright');
figure,imshow(W);% title('Sobel gradient');
viscircles(centers,radii,'Color','b');

请注意,图片

之间的圆圈大小可能会有所不同

任何建议。

2 个答案:

答案 0 :(得分:4)

如此接近!你只需要调整一下灵敏度。范围从0到1,找到的圆圈越多越高。默认情况下为0.85。还得调整半径:

imshow(W) % assuming W is logical
[centers,radii] = imfindcircles(W,[40 60],'ObjectPolarity','bright','Sensitivity',0.925)
viscircles(centers, radii,'EdgeColor','b');

我不认为有可能找到较低的圆圈而不会先得到一些误报 - 最终它不是很圆。

resultimage

答案 1 :(得分:1)

尽管不是一位优秀的数字图像分析师,但我试图用这张图片稍微玩一下。我发现以下方法更加灵活:

img = imread('1GBtO.png');
img = im2bw(img);

rp = regionprops('table',img,'Centroid','Eccentricity','EquivDiameter','MajorAxisLength','MinorAxisLength');
rp.Radii = mean([rp.MajorAxisLength rp.MinorAxisLength],2) / 2;

rp([rp.Eccentricity] > 0.75,:) = [];
rp(rp.EquivDiameter < 10 | rp.EquivDiameter > 30,:) = [];

imshow(img);
hold on;
viscircles(rp.Centroid,rp.Radii);
hold off;

输出远不是理想的输出,但我觉得它比你的第一次试用有所改进:

Output

但我认为您可以通过找到正确的算法来调整半径并通过调整标准来排除不需要的区域来获得所需的结果。