Extract non-touching regions from image

1 次查看（过去 30 天）
GIONA CANEVER
GIONA CANEVER 2021-4-6
编辑： DGM 2021-4-6
I've got the following BW image of a flower. I want to extract the non-touching area so for example the isolated region on the right of the image is a separate area from the others.
If I use bwconncomp or bwlabel I obtain 142 regions but there are some regions separated by some white pixels which are considered unite. I can see this by plotting the different area's centroid; as you can see in the second image, the separated area on the right does not have it's own centroid because this area is considered to belong to another.
How can I obtain the separated area and then eventually not consider the regions with dimensions under a certain value?
this is my code
flower = imread('flower.png');
figure
imshow(sestoBW);
flowerGR = im2gray(flower);
flowerBW = imbinarize(flowerGR);
CC = bwconncomp(flowerBW);
stats = regionprops(flowerBW,'Centroid');
centroids = cat(1,stats.Centroid);
hold on
plot(centroids(:,1),centroids(:,2),'b*')
thank you in advance for your help
  3 个评论
显示 1更早的评论
GIONA CANEVER
GIONA CANEVER 2021-4-6
编辑：GIONA CANEVER 2021-4-6
Yes, thank you very much, this is what I was looking for.
Now that I've found all the separated regions, how can I individually extract them? So for example, considering this region on the right, I obtain a new image of the same dimensions of the original one all white except in this area on the right?
something like this
DGM
DGM 2021-4-6
编辑：DGM 2021-4-6
something like this:
inpict=rgb2gray(imread('flower.png'))>0.5;
L=bwlabel(~inpict,8);
% that particular spot is label #90
thisobject=(L==90);
Once you have the label array from the image, you'll have to find the label associated with a particular object. I just did imshow(L) and used the datatip tool.
Once you know what the label is, you can just do a logical test on the label array.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

LO
LO 2021-4-6
编辑：LO 2021-4-6
You can try using imerode on the blobs.
eroded_image = imerode(binary_image,[1 0 1; 1 0 1; 1 0 1]); % this will "shrink" the image on the vertical axis, you can reshape the matrix of ones to set a different rule for erosion: [1 1 1],[1 0 1],[1 1 1] this would be a circular erosion for example
Alternatively you can impose constraints on the distances that centroids should have (all center points having euclidean distance higher than a certain value, are excluded).
Alternatively you could "erode" the image (or threshold the binary image) before using regionprops (if it is not important you conserve the precise sizing of the blobs).
These are some snippets of some code I have been using, perhaps they can help you:
to remove blobs smaller than a certain threshold and if you know the N of blobs you want
CC = bwconncomp(binary_image); % connected shapes in your binary
objects = CC.NumObjects; % number of blobs
while objects > 2; % let's say you expect 2 blobs maximum
numPixels = cellfun(@numel,CC.PixelIdxList);
[smallest,idx] = min(numPixels);
binary_image(CC.PixelIdxList{idx}) = 0;
CC = bwconncomp(binary_image);
objects = CC.NumObjects;
end
to find the coordinates of the points of a given blob, you can use the following code
for i = 1:objects
cc = bwconncomp(binary_image); % find connected components in the image (it should find 2 blobs)
index = cc.PixelIdxList{i}; % indexes of pixels in blob1
[ovlogic,ovidx]=ismember(pixel_idx,index); % find ovarlapping indices in the matrix containing al pixels of all blobs
blob = pixels(ovlogic,:); % pixel coordinates of blob, save this somewhere if you need
end

类别

在 Help Center 和 File Exchange 中查找有关 Image Segmentation and Analysis 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by