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!