12
$\begingroup$

I have an image:

enter image description here

I tried to extract the outline of the bubbles, but I met with some problems. Firstly I tried the method used in Mathematica documentation for Analyze Segmented Cells In An Image, but it didn't work well, I got an irregular figure:

enter image description here

Then I tried the function EdgeDetect. Things got a little bit better, but it's still not what I want. Could you please help me with my problem?

$\endgroup$
5
  • $\begingroup$ One question is.., what do you want out of it? The number of bubbles? $\endgroup$ – Öskå Jul 18 '14 at 17:46
  • $\begingroup$ I want the outline of each bubble. $\endgroup$ – nix Jul 18 '14 at 17:49
  • 1
    $\begingroup$ I've got something like this not perfect, but better :) $\endgroup$ – Öskå Jul 18 '14 at 18:30
  • $\begingroup$ @Öskå I think it's much better than my picture! Could you please share your method with us? $\endgroup$ – nix Jul 19 '14 at 2:26
  • $\begingroup$ Rahul's answer is even better ;o) Quite the same method though :) $\endgroup$ – Öskå Jul 19 '14 at 9:06
17
$\begingroup$
image = Import["http://i.stack.imgur.com/UNXEV.jpg"];
n = TotalVariationFilter[image, 0.05]
b = LocalAdaptiveBinarize[n, 50]
s = DeleteSmallComponents[ColorNegate[DeleteSmallComponents[b, 20^2]], 20^2]
(d = DistanceTransform[s]) // ImageAdjust
m = MaxDetect[d, 10]
(w = WatershedComponents[ColorNegate[d], m]) // Colorize
c = ComponentMeasurements[w, {"Centroid", "EquivalentDiskRadius"}, #2 < 50 &];
Show[image, Graphics[{Thick, Red, Circle @@@ c[[All, 2]]}], ImageSize -> 1280]

enter image description here

(Click for bigger.)

Run the code to see the intermediate steps.

Thanks to nikie for the second half of this answer.

LocalAdaptiveBinarize is new in version 10, but all it does is "creates a binary image by replacing values above the mean of the range-$r$ neighborhood with $1$ and others with $0$." So you can replace it with something like, say,

b = Binarize[ImageSubtract[n, MeanFilter[n, 50]], 0]

which gives similar results, though not identical for some reason. Possibly there's a more elegant way to do it.

$\endgroup$
3
  • $\begingroup$ Do you think that the same method could work with Polygons? I asked this question hoping to get the same kind of answer :) $\endgroup$ – Öskå Jul 19 '14 at 9:07
  • $\begingroup$ Thank you for your answer! I think // ImageAdjust // Colorize are of no use during the processing, since d and w aren't changed by these 2 commands. So you just put them here to show the effect of ImageAdjust and Colorize? What's more, for users using Mathematica below 10: Note that function LocalAdaptiveBinarize appears only in Mathematica 10, which means that maybe you cannot use LocalAdaptiveBinarize. For concrete solution, Rahul has showed above. $\endgroup$ – nix Jul 20 '14 at 7:42
  • $\begingroup$ I put the // ImageAdjust and // Colorize there to visualize what the intermediate steps are doing, because you can't tell anything from looking at d and w directly. I've updated my answer with a replacement for LocalAdaptiveBinarize. $\endgroup$ – user484 Jul 20 '14 at 8:56

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.