How can I fit a circle to a region of pixels?

Question

Kelsey Hilton 2016-10-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/306198-how-can-i-fit-a-circle-to-a-region-of-pixels

回答： Gianluca Iori 2018-7-24

2D_trabecularBoneSlice.png

I have a series of 2D CT slices of trabecular bone. I have binarized them and detected the boundaries of the bone, and now I would like to do a bit more processing. I am trying to fit the largest possible circle into each slice in the black pixel regions. I've included an image so you can see what I mean.

I've tried something like this:

stats = regionprops('table',bw,'Centroid','EquivDiam');

where bw is my binary image without the boundaries drawn in.

Any help would be great! I'm struggling to find other examples which match mine.

Thanks in advance!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Massimo Zanetti 2016-10-7

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/306198-how-can-i-fit-a-circle-to-a-region-of-pixels#answer_237919

You can manage by exploiting the "shape measurements" and other properties that function regionprops allows you to use:

https://it.mathworks.com/help/images/ref/regionprops.html

3 个评论
显示 1更早的评论隐藏 1更早的评论

Kelsey Hilton 2016-10-7

test.png

I've attached an example of what I get if I use that option. Maybe it'll help.

Massimo Zanetti 2016-10-7

在 MATLAB Online 中打开

You are considering "on" pixels the white ones. Try negate the input of bwconncomp. For example instead of using:

C = bwconncomp(BW)

use

C = bwconncomp(~BW)

This way you will consider the black areas.

请先登录，再进行评论。

Answer 2

Bert 2016-10-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/306198-how-can-i-fit-a-circle-to-a-region-of-pixels#answer_237925

在 MATLAB Online 中打开

I'm not familiar with this toolbox. But have you tried simply scanning all point? you could do something like this:

% create indexes
ix = repmat([1:size(bw,2)],size(bw,1),1);
iy = repmat([1:size(bw,1)]',1,size(bw,1));
% all in vector
bw = bw(:);
ix = ix(:);
iy = iy(:);
% init loop
center = [];
radius = 0;
% indexes of all black points, possibally the center
I = 1:length(bw);
I = I(bw==0);
for i = I
    % all distances from thi point to white ones
    distances = sqrt( (ix(i)-ix(bw==1)).^2 + (iy(i)-iy(bw==1)).^2 );
    % smallest distance
    distance = min(distances);
    if distance>radius % this point has the biggest minimal distance to a white one yet
        radius = distance;
        center = [iy(i) ix(i)];
    end
end