绘制线并切断循环区域

Question

运行以下代码后，我得到了以下图片。

file='grayscale.png';
I=imread(file);
bw = im2bw(I);
bw = bwareaopen(bw,870);
imwrite(bw,'noiseReduced.png')
subplot(2,3,1),imshow(bw);
[~, threshold] = edge(bw, 'sobel');
fudgeFactor = .5;
im = edge(bw,'sobel', threshold * fudgeFactor);
subplot(2,3,2), imshow(im), title('binary gradient mask');

se = strel('disk',5);
closedim = imclose(im,se);
subplot(2,3,3), imshow(closedim), title('Connected Cirlces');
cc = bwconncomp(closedim);

S = regionprops(cc,'Centroid'); //returns the centers S(2) for innercircle
numPixels = cellfun(@numel,cc.PixelIdxList);
[biggest,idx] = min(numPixels);
im(cc.PixelIdxList{idx}) = 0;
subplot(2,3,4), imshow(im), title('Inner Cirlces Only');
c = S(2);

我的目标是在圆形物体周围绘制一个红色圆圈（参见图像）并从原始图像“I”中剪切圆形区域（区域）并将裁剪区域保存为图像或执行其他任务。我该怎么办？

Answer 1

白色像素的凸包将为您提供相当好的圆近似值。你可以找到中心作为船体区域的质心，半径作为从中心到船体顶点的平均距离。

Answer 2

或者，您可以使用包含所有点的最小r来优化/拟合圆圈：

bw = imread('http://i.stack.imgur.com/il0Va.png');
[yy xx]=find(bw);

现在，让p为参数化圆的三个向量：p(1), p(2)是中心的x-y坐标，其半径为p(3)。然后我们想要最小化r（即p(3)）：

obj = @(p) p(3);

受制于圈内的所有点

con = @(p) deal((xx-p(1)).^2+(yy-p(2)).^2-p(3).^2, []);

使用fmincon进行优化：

[p, fval] = fmincon(obj, [mean(xx), mean(yy), size(bw,1)/4], [],[],[],[],[],[],con);

产量

p =
471.6397  484.4164  373.2125

绘制结果

imshow(bw,'border','tight');
colormap gray;hold on;
t=linspace(-pi,pi,1000);
plot(p(3)*cos(t)+p(1),p(3)*sin(t)+p(2),'r', 'LineWidth',1);

你可以生成一个与bw大小相同的二进制掩码，其中true在圈内，false在外面

msk = bsxfun(@plus, ((1:size(bw,2))-p(1)).^2, ((1:size(bw,1)).'-p(2)).^2 ) <= p(3).^2;

面具看起来像：