Restricting the maximum boundary of region growing algorithm

2 views (last 30 days)
Hi,
I am trying to implement a multilayer region growing algorithm for breast tumor segmentation. This is a reference of a previous question that I asked.
TM25.jpg is the original image. I want to restrict the boundary of the region growing algorithm. The boundary is defined as the rectangular region in e2.jpg. I want to manually select the seed point inside the rectangular ROI. However, after implementing the multilayer region growing algorithm , I get the result shown in e1.jpg which is not correct. Any suggestions would be appreciated.
I=im2double(imread('TM25.jpg'));
I=rgb2gray(I);
imshow(I)
k = waitforbuttonpress;
point1 = get(gca,'CurrentPoint'); % button down detected
finalRect = rbbox; % return figure units
point2 = get(gca,'CurrentPoint'); % button up detected
point1 = point1(1,1:2); % extract x and y
point2 = point2(1,1:2);
p1 = min(point1,point2); % calculate locations
offset = abs(point1-point2); % and dimensions
x = [p1(1) p1(1)+offset(1) p1(1)+offset(1) p1(1) p1(1)];
y = [p1(2) p1(2) p1(2)+offset(2) p1(2)+offset(2) p1(2)];
hold on
axis manual
plot(x,y)
[x,y]=getpts;x=round(x);y=round(y);
a=imgaussfilt(I,2);
b=adapthisteq(a);
m=regiongrowing_MLT(b,x,y,12);
m=imfill(m,'holes');
figure,imshow(m)

Accepted Answer

Image Analyst
Image Analyst on 27 Mar 2021
You can call regionprops() on the regions and ask it for the perimeter. Take appropriate action if the perimeter is more than you like
props = regionprops(mask, 'Perimeter');
allPerimeters = [props.Perimeter];
if any(allPerimeters) > maxAllowablePerimeter
% Do something...
end
  7 Comments

Sign in to comment.

More Answers (0)

Products


Release

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!