how can imfill an image with a lot of edge?

7 次查看(过去 30 天)
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 个评论
Sean de Wolski
Sean de Wolski 2014-9-8
编辑:Sean de Wolski 2014-9-8
Show us the original image; it might be easier to start with that for the original segmentation.
sara
sara 2014-9-8
hi Sean
my original image is .tif so I attach this in a .zip
thanks

请先登录,再进行评论。

回答(2 个)

Image Analyst
Image Analyst 2014-9-8
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 个评论
sara
sara 2014-9-14
编辑:sara 2014-9-14
hi Image Analyst
I thouth our code is correct for every thing but when I tested it for this image that I attached I understood that the lung result should not have any holes but if you see this image and the result you can see in the bottom and left side of this lung we have a big hole...do you agree with me that if before using ExtractNLargestBlobs we use a circle mask the result can be better?
Image Analyst
Image Analyst 2014-9-14
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.

请先登录,再进行评论。


sara
sara 2014-9-11
hi Image Analyst I use otsu method instead of <45000 in this code and it gave me a good result I attach this.

Community Treasure Hunt

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

Start Hunting!

Translated by