Separating objects in an image

Question

Zoran 2014-8-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/151923-separating-objects-in-an-image

评论： Image Analyst 2017-7-8

Hi,

I'm trying to count objects in an image and the matlab script I wrote works pretty well but some objets are touching (in other image they are even overlapping) and it counts only for one instead of two. I am trying to use the bwdist function associated to the watershed function to separate these objects but it doesn't work pretty well: it cuts my objets at many part where it shouldn't and the counting is much worse.

If somebody could explain me how to proceed in other to separate only the good particles (the 2 touching at the right of my image) I would be grateful.

link for image: download the plastic image

Here is my code (if you run it a lot of figures will appear):

Sorry for writing the comments in french =P

clear all;
k=1; % indice pour la numérotation des images
I=imread('Images particules/DSC_0037.jpg');
figure(k)
k=k+1;
imshow(I);
I_hsv=rgb2hsv(I);
figure(k)
k=k+1;
imshow(I_hsv);
I_h=I_hsv(:,:,1);
I_s=I_hsv(:,:,2);
I_v=I_hsv(:,:,3);
figure(k)
k=k+1;
imshow(I_h)
figure(k)
k=k+1;
imshow(I_s)
figure(k)
k=k+1;
imshow(I_v)
%%Hue
[Gx, Gy] = imgradientxy(I_h);
[Gmag, Gdir] = imgradient(Gx, Gy);
% figure(k);
% k=k+1;
% imshowpair(Gmag, Gdir, 'montage');
I_bw1=Gmag>mean(quantile(Gmag,0.99));
figure(k);
k=k+1;
imshowpair(Gmag,I_bw1,'montage');
%%Saturation
[Gx, Gy] = imgradientxy(I_s);
[Gmag, Gdir] = imgradient(Gx, Gy);
% figure(k);
% k=k+1;
% imshowpair(Gmag, Gdir, 'montage');
I_bw2=Gmag>mean(quantile(Gmag,0.99));
figure(k)
k=k+1;
imshowpair(Gmag,I_bw2,'montage');
%%Variance
[Gx, Gy] = imgradientxy(I_v);
[Gmag, Gdir] = imgradient(Gx, Gy);
% figure(k);
% k=k+1;
% imshowpair(Gmag, Gdir, 'montage');
I_bw3=Gmag>mean(quantile(Gmag,0.99)); % choisir le bon quantile
figure(k)
k=k+1;
imshowpair(Gmag,I_bw3,'montage');
%%Addition images du gradient
I_recomp=I_bw1+I_bw2+I_bw3;
figure(k);
k=k+1;
imshow(I_recomp);
%%Dilatation - fill - erosion
% Element structurant diamond
% Dilatation
SE=strel('octagon',3); % doit être un multiple de 3 !
I_dil=imdilate(I_recomp,SE);
% figure(k)
% k=k+1;
% imshow(I_dil);
% Fill
I_fill=imfill(I_dil,'holes');
% Erosion
I_er=imerode(I_fill,SE);
% figure(k)
% k=k+1;
% imshow(I_er);
%%Elimination du bruit en appliquant un imerode <taille minimale des plastiques en pixels
% Erosion - dilatation
SE=strel('octagon',6); % mesurer la taille maximale d'un plastic en pixel avec imdistline !
I_bruit=imdilate(imerode(I_er,SE),SE);
figure(k)
k=k+1;
imshow(I_bruit);
%%Séparation des particules avec watershed
I_bwdist=-bwdist(~I_bruit);
figure(k);
k=k+1;
imshow(I_bwdist,[]);
I_water=watershed(I_bwdist);
I_bruit(I_water==0)=0;
figure(k);
k=k+1;
imshow(I_bruit);
%%Comptage des particules
cc=bwconncomp(I_bruit);
cc.NumObjects
L=labelmatrix(cc);
RGB_label=label2rgb(L);
figure(k);
k=k+1;
imshow(RGB_label);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Image Analyst 2014-8-22

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/151923-separating-objects-in-an-image#answer_149495

imerode then imdilate is the same as imclose. Anyway, try watershed segmentation.

http://blogs.mathworks.com/steve/2013/11/19/watershed-transform-question-from-tech-support/

http://www.mathworks.com/help/images/examples/marker-controlled-watershed-segmentation.html

http://cmm.ensmp.fr/~beucher/wtshed.html

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Zoran 2014-8-22

在 MATLAB Online 中打开

Thank you for these precious links. Watershed segmentation seems is certainly the answer to my problem.

The problem is I got an over segmentation while applying the function. I tried to eliminate some "holes" that are responsible for that using imextendmin and imimposemin and the result is some what better. But some of my plastic are now cut in new pieces (they should stay compact) and the waterline is cutting my pieces also.

I don't see how I can do to get only one minima per particle. And how can I apply foreground and background marker method to my problem ?

Here's what I changes in my script:

I_bwdist=-bwdist(~I_bruit);
figure(k);
k=k+1;
imshow(I_bwdist,[]);
mask = imextendedmin(I_bwdist,25); % choisir la bonne valeur
figure(k);
k=k+1;
imshowpair(I_bruit,mask,'blend');
I_bwdist2 = imimposemin(I_bwdist,mask);
I_water=watershed(I_bwdist2);
I_bruit(I_water==0)=0;
figure(k);
k=k+1;
imshow(I_bruit);

请先登录，再进行评论。

Answer 2

LAGARI A 2017-7-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/151923-separating-objects-in-an-image#answer_273376

Hai zoran, actually i was searching for the same algorithm. i want to save objects in a image seperately based on its color. object is differentiates based on color. can u help me with this.each object in the image should save as each file.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Image Analyst 2017-7-8

See my File Exchange: http://www.mathworks.com/matlabcentral/fileexchange/?term=authorid%3A31862&sort=downloads_desc, particularly the Image Segmentation Tutorial where I crop out regions of interest, and the color segmentation demos that let you identify regions of certain colors.

请先登录，再进行评论。