removing some connected component
38 views (last 30 days)
In Image Processing, ''imclearborder'' is used to remove connected components that touch the borders of the image.
But i want to remove some specific, for example in my condition it will be like this '' if (49,90) pixel is in any connected component then remove it like imclearborder does.
And also sometimes i may not want to remove all components but only those which touch right border of the image or up border etc.
How can i make my special adjustments? Is it concerned with imclearborder(IM,4) OR imclearborder(IM,8)?
Thanks in advance.
David Legland on 13 Jul 2012
Edited: David Legland on 13 Jul 2012
if your input is a label image, you can use the following syntax:
img(img == ind) = 0;
This will detect which pixels are labeled with label IND, and set them to 0, which is the usual label for background. You can convert binary image to label image by using the bwlabel function.
If you want to remove the region around the point (40, 50):
ind = img(50, 40); % beware of index ordering
img(img == ind) = 0;
If you want to remove labels touching right border:
inds = unique(img(:, end));
inds(inds==0) = 0;
img(ismember(img, inds)) = 0;
Not that the bwconncomp function may be useful too. The output structure is different, and may be more efficient in some cases.
More Answers (2)
Image Analyst on 13 Jul 2012
See my image segmentation tutorial: http://www.mathworks.com/matlabcentral/fileexchange/?term=authorid%3A31862 It shows you how to filter your blobs based on things like their area or intensity (or whatever) and then give you back the filtered labeled image without those blobs using the ismember() function. They key snippet of code is below:
% Now I'll demonstrate how to select certain blobs based using the ismember function.
% Let's say that we wanted to find only those blobs
% with an intensity between 150 and 220 and an area less than 2000 pixels.
% This would give us the three brightest dimes (the smaller coin type).
allBlobIntensities = [blobMeasurements.MeanIntensity];
allBlobAreas = [blobMeasurements.Area];
% Get a list of the blobs that meet our criteria and we need to keep.
allowableIntensityIndexes = (allBlobIntensities > 150) & (allBlobIntensities < 220);
allowableAreaIndexes = allBlobAreas < 2000; % Take the small objects.
keeperIndexes = find(allowableIntensityIndexes & allowableAreaIndexes);
% Extract only those blobs that meet our criteria, and
% eliminate those blobs that don't meet our criteria.
% Note how we use ismember() to do this.
keeperBlobsImage = ismember(labeledImage, keeperIndexes);
% Re-label with only the keeper blobs kept.
labeledDimeImage = bwlabel(keeperBlobsImage, 8); % Label each blob so we can make measurements of it
% Now we're done. We have a labeled image of blobs that meet our specified criteria.
subplot(3, 3, 7);
title('"Keeper" blobs (3 brightest dimes in a re-labeled image)');