how can imfill an image with a lot of edge?

30 views (last 30 days)
sara
sara on 8 Sep 2014
Commented: Image Analyst on 20 May 2017
hi
I want to imfill an image of lung to have an image that shows me just left and right lung.. I obtain boundary of lung wit this cod that I attach
if true
function [perimImage] = filling( I )
s = uint8(255 * mat2gray(I));
s = im2bw(s);
perimImage = xor(s,imerode(s, true(3)));
figure, imshow(perimImage)
end end
when I use this code my result is like the image that I put.
I think it is a good result but when I use imfill, the result is this:
do you have any suggestion for me to have an image like this after using imfill?
be black or white is not important for me I just want to segment lung to the result that I say.. thanks
  2 Comments
sara
sara on 8 Sep 2014
hi Sean
my original image is .tif so I attach this in a .zip
thanks

Sign in to comment.

Answers (2)

Image Analyst
Image Analyst on 8 Sep 2014
Draw lines around the border of your binary image:
binaryImage(:,1) = true;
binaryImage(:,end) = true;
binaryImage(1,:) = true;
binaryImage(end,:) = true;
Now do a fill to fill the region between the border and the body.
binaryImage = imfill(binaryImage, 'holes');
Now do an imclearborder to get rid of all that and leave just the lung edges
binaryImage = imclearborder(binaryImage);
Now do a fill again to get the lungs filled.
lungsMask = imfill(binaryImage, 'holes');
And you're basically done. If there is any small blobs remaining, use my ExtractNLargestBlobs function (attached) to extract just the 2 largest blobs, which will be the lungs.
  14 Comments
Image Analyst
Image Analyst on 14 Sep 2014
It's not a "hole" but more of a "bay" because that region is actually a white region that is an incursion from the bright outer body into the lung. If you want to smooth the boundary you need to use active contours. This will "cut off" that bay. I'm attaching an example.

Sign in to comment.


sara
sara on 11 Sep 2014
hi Image Analyst I use otsu method instead of <45000 in this code and it gave me a good result I attach this.
  2 Comments
Image Analyst
Image Analyst on 20 May 2017
Who are you asking? I gave my code above in this post: click here

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!