# Cell segmentation119

Posted by Steve Eddins,

Blog reader Ramiro Massol asked for advice on segmenting his cell images, so I gave it a try. I'm not a microscopy expert, though, and I invite readers who have better suggestions than mine to add your comments below.

Let's take a look first to see what we have. I'm going to work with a cropped version of the original so that the images aren't too big for the layout of this blog.

Note: you can download the functions imcredit and imoverlay from MATLAB Central.

I = imread('http://blogs.mathworks.com/images/steve/60/nuclei.png');
I_cropped = I(400:900, 465:965);
imshow(I_cropped)
imcredit('Image courtesy of Dr. Ramiro Massol')

Strictly speaking, contrast adjustment isn't usually necessary for segmentation, but it can help the algorithm developer see and understand the image data better. This is a fairly low-contrast image, so I thought it might help. You can adjust the display contrast interactively with imtool, or you can use an automatic method such as adapthisteq. adapthisteq implements a technique called contrast-limited adaptive histogram equalization, or CLAHE. (I always thought "CLAHE" sounded like it must be some Klingon delicacy.)

I_eq = adapthisteq(I_cropped);
imshow(I_eq)

So what happens if we just apply a threshold now?

bw = im2bw(I_eq, graythresh(I_eq));
imshow(bw)

Let's clean that up and then overlay the perimeter on the original image.

bw2 = imfill(bw,'holes');
bw3 = imopen(bw2, ones(5,5));
bw4 = bwareaopen(bw3, 40);
bw4_perim = bwperim(bw4);
overlay1 = imoverlay(I_eq, bw4_perim, [.3 1 .3]);
imshow(overlay1)

Now, I'm not familiar with these cell images, so I don't know exactly what I'm looking at. I assume some of these blobs need more help to be separated properly. One possible approach is called marker-based watershed segmentation. There's a demo of this idea on The MathWorks web site.

With this method, you have to find a way to "mark" at least a partial group of connected pixels inside each object to be segmented. You also have to mark the background.

Let's try to use the bright objects, which I assume are nuclei. The extended maxima operator can be used to identify groups of pixels that are significantly higher than their immediate surrounding.

mask_em = imextendedmax(I_eq, 30);
imshow(mask_em)

Let's clean that up and then overlay it.

mask_em = imclose(mask_em, ones(5,5));
overlay2 = imoverlay(I_eq, bw4_perim | mask_em, [.3 1 .3]);
imshow(overlay2)

Next step: complement the image so that the peaks become valleys. We do this because we are about to apply the watershed transform, which identifies low points, not high points.

I_eq_c = imcomplement(I_eq);

Next: modify the image so that the background pixels and the extended maxima pixels are forced to be the only local minima in the image.

I_mod = imimposemin(I_eq_c, ~bw4 | mask_em);

Now compute the watershed transform.

L = watershed(I_mod);
imshow(label2rgb(L))

I don't know if this is a good segmentation result or not, but I hope some of the methods I've shown will give Dr. Massol some ideas to try.

Get the MATLAB code

Published with MATLAB® 7.2

### Note

kruthika replied on : 1 of 119

hy…im sorry ….amnt here to give a comment…i am working on splittin and merging algorithms..i find there is some error with the splitmerge.m function given in ‘Digital Image Processing using MATLAB’…Could u please help me…i ve just learnt the basics of image segmentation….

thnk u…

Steve replied on : 2 of 119

Kruthika, Rafael Gonzalez knows more about the splitmerge function than I do. You can find contact information for Professor Gonzalez at http://www.imageprocessingplace.com.

ishraq replied on : 3 of 119

hi everbody,,
i have a qustion about coins segmentation ,,
i didn’t find a good way to segment the coins when they are
connected (the edges of coins are connected).

any one can help me please???!?!?!

Steve replied on : 5 of 119

Ishraq – I took a quick peek at your coins image. Possibly something as simple as eroding a thresholded image might work. If not, then you might try marker-controlled watershed segmentation.

ishraq replied on : 6 of 119

Dear Mr Steve ,
Thanks alot for your fast replaying ,
but I still have difficulties in using watershed because
im beginner in using MATLAB,

I like to add 2 new urls for coin images with overlapping between the edges

my problem is about how i can segment these coins in a good way so i can use “bwlabel” and “regionprops”
to recognize the coin types .

i hope from you to help me to find a way to segment these coins because you are expert person .

thanks you alot
ishraq

Steve replied on : 7 of 119

Ishraq – I’m sorry, but I can’t work this problem out for you, or help you learn MATLAB. Go through the Getting Started in MATLAB guide to get some MATLAB knowledge. Then look at the Image Processing Toolbox Users Guide and the product demos on our web site, including the one I already linked to. Look at the material in the documentation and demos on thresholding, morphological operators such as erosion and opening, and watershed segmentation.

Wilson replied on : 8 of 119

Hi there,
I have a query on how to use Matlab for Vehicle Plate segmentation?
I read from the papers that one solution is to gray scale the picture. But how does the system be so clever to segment out the unwanted parts and leave the vehicle plate untouched?
Hope to hear from you guys soon.
Cheers..

subrajeet mohapatra replied on : 9 of 119

Hi,
I have read your book but there is nothing about adaptive histogram equalization.Can u give me some ideas or the algorithm for adaptive histogram equalization.So that i can go to CLAHE.

Steve replied on : 10 of 119

Subrajeet – There is an implementation of CLAHE in the Image Processing Toolbox. See the adapthisteq function.

subrajeet mohapatra replied on : 11 of 119

can i get the complete algorithm

Steve replied on : 12 of 119

Subrajeet – The complete algorithm is summarized at the top of the M-file, plus you can always look at the code. The M-file also contains this reference: Karel Zuiderveld, “Contrast Limited Adaptive Histogram Equalization”, Graphics Gems IV, p. 474-485, code: p. 479-484

subrajeet mohapatra replied on : 13 of 119

Hi,
is there any built in function to measure signal to noise ratio for images ,My friend is telling there is Peak signal to noise ratio (PSNR) im image processing toolbox

Steve replied on : 14 of 119

Subrajeet – There’s no PSNR function in the Image Processing Toolbox, but it’s very easy to implement it directly from the formula.

Fabiano Malhard replied on : 15 of 119

Hi steve,

I am really new at MatLab and I am trying to work with images, reducing its size and later resizing them using one of the commons methods … but i would like to compare the resulting images with the original ones …

Can you help with the SNR function, or there is any other better aproach ???

Thanks

Steve replied on : 16 of 119

Fabiano – There are several different definitions of SNR used for image comparisons, so the details vary depending on the definition you are interested in. But the heart of the computation, typically something like a mean-squared error, is straightforward in MATLAB. If you are comparing A and B, it might look something like this:

mse = mean( (A(:) - B(:)).^2 );

rosani replied on : 17 of 119

Steve,
how to compute a confusion matrix to measure accuracy between two image(rgb,in pixel) from different source.

Steve replied on : 18 of 119

Rosani—I’m confused; what’s a confusion matrix?

bob gold replied on : 19 of 119

Hi my name is Bob Gold a pharmacist in Indiana. I
am interested in developing a method to determine the
number of tablets in a bottle for the purpose of
improving compliance. I am not an expert on cameras .

Do you have any ideas.
Thanks

Rob replied on : 20 of 119

Bob Gold,

I think first of all you’ll need some really wide, flat bottles so the pills are no more than 1 deep in the image. (sorry for the sarcasm, couldn’t resist)

Other than that, you could take a picture of pills in a regular bottle and compute an *estimate* of how many pills are packed in there. Depending on the pill shape, there are many standard packing computations, such as hexagonal-close-packed, body-centered-cubic, etc. Think lattice structures, but all the ones I’ve studied are for spheres. The main failing here is that you won’t get the real number, which it sounds like you need for compliance.

Hope this helps,
Rob

Steve replied on : 21 of 119

Bob—I don’t anything to concrete to suggest beyond what Rob said. It sounds like it might not be feasible except under very particular circumstances.

Rob—Thanks for jumping in.

Steve replied on : 22 of 119

Bob—OK, I do have one thought. I’d be more tempted to use a precision scale instead of a camera and image processing to measure the number of pills. But this kind of quality assurance application is a bit beyond my expertise.

Adnan Butt replied on : 23 of 119

Hello Sir,
I am working in MATLAB image processing and one thing I do have explored that most of the commands extracting quantitaive measures work only on 2D matrices there is not much material on 3D image processing in MATLAB?
Can you guide me in 3d Image processing in MATLAB?
Regards

Steve replied on : 24 of 119

Adnan—There are many functions in MATLAB and the Image Processing Toolbox that work on three-dimensional arrays, including filtering, transforms, morphology, and region measurements. Can you be more specific about what you are looking for?

Ben Taylor replied on : 25 of 119

Improvement: Blobs/particles touching the sides of the images will mess up your size measurements and shape characterization because they are not showing the whole particle. To remove all particles or cells touching the sides I wrote the following simple script that removes any particle touching the border of the image.

%%%%////
function L=borderstrip(L)
display(‘starting border strip’)
border=[L(end,:) L(1,:) L(:,end)’ L(:,1)’]; %these are all the pixels for the borders
background=mode(mode(L)); %this finds the most common pixel for the background which is usually 2

count=0;
for i=1:length(border)
if border(1,i)==2
border(1,i)=0;
else
count=count+1;
blacklist(count)=border(1,i);
end
end
blacklist=unique(blacklist); % this strips out all of the copied black listed pixels to be removed
height=size(L,1);width=size(L,2);
for j=1:length(blacklist)
for h=1:height
for w=1:width
if L(h,w)==blacklist(j)
L(h,w)=background;
end
end
end
end
display(‘ending border strip’)

Steve replied on : 26 of 119

Ben—You can use imclearborder to remove objects touching the border.

Saad Khan replied on : 27 of 119

Hi Steve,
Do you know of any matlab implementation of 3D deblurring (deconvolution)? The matlab image processing toolbox contains four deblurring functions but all of these work with images (2D data).
Thanks!

Steve replied on : 28 of 119

Saad—You can use the deconvolution routines in the Image Processing Toolbox. They are not limited to two dimensions.

VEN replied on : 29 of 119

Respected Sir,
I am Working in Matlab with Some Videos to Identify Moving
Object in that Video. I am Working in Matlab7 in which Video
Processing is not There. Now I am able to Identify the
moving objects Partially. Could you please Suggest me to Segment only that Moving objects In that Video. And also I want to Count the Number of segmented Object. Please help me in this matter.

Thanking You,

maged replied on : 30 of 119

i want functions that make segmantation of nodules in image by matlab
thanks

Shalin replied on : 31 of 119

3D deconvolution:

Hi Steve! Great page… Continuing on Saad’s query, could you help me with pointers on how to represent 3D point spread functions(in my case for optical microscope) and how to employ blind deconvolution.

Thanks
Shalin

Steve replied on : 32 of 119

Shalin—A 3D point spread function is represented as a 3D array containing the impulse response of the blur operator. There is an Image Processing Toolbox function for performing blind deconvolution. Did you have a specific question about it?

Stalin Mohapatra replied on : 33 of 119

Can you plz give me a complete MATLAB code on MEAN Filtering of Noisy Images?
It would be kind enough if you send it to my mail. I’ll be ever grateful to u.

Steve replied on : 34 of 119

Stalin—It sounds like you might just want to use imfilter (or the MATLAB function conv2) with a constant filter, such as ones(5,5)/25.

Adnan replied on : 35 of 119

Hi Steve, I posted the same query a couple of days ago but I cant see my Post now. Any way, There is a camera on road fixed in one position. I want to count the vehicles passed by the camera. But the problem is the segmentation i.e. to ignore any thing except vehicle. Kindly suggest me some proper algorithm for this problem.
Thanking in Anticipation
Regards

Eric replied on : 36 of 119

Hi Steve, I am trying to create circular and elliptical images as quality assurance tests for my image segmentation code. I can do a simple generate xy coordinates from parametric equations command, but I get a real blotchy image. I’ve figured out that it is because I’m not implementing Bresenham style line algorithums. Does matlab have built in Bresenham line drawing algorithums that can be augmented eaisly to draw nice smooth circles and ellipses in pixelated images?

Steve replied on : 37 of 119

Eric—As far as I understand, Bresenham drawing algorithms are for drawing shapes blazingly fast, with an absolute minimum of arithmetic, in low-level code. I’m not sure such a thing is really needed in MATLAB, where we get to do all the math we want. For example, the file toolbox/images/images/private/intline.m draws the same line between two pixels that a Bresenham line-drawing routine but, it uses some simple math to do it. So … what do you mean by “blotchy”?

I would probably “draw” a circle on an image using meshgrid and some sort of distance logic, like this:

[x,y] = meshgrid(linspace(-1,1,200));
bw = hypot(x,y) < = 0.25;

Wander replied on : 38 of 119

Dear Sir,
I want to find out the circularity and cell count of cell images after employing Marker Controlled Watershed transform. (Images are similar to those taken as illustration in this web-page). Can you kindly explain to me how to go ahead in this direction.

Steve replied on : 39 of 119

Wander—Use bwlabel and regionprops.

walid jerbi replied on : 40 of 119

Hi Steve,
I need to merge two regions in a labeled image. I want to know if there is a function that allow mes to do this?
Regards

Steve replied on : 41 of 119

Walid—Can you provide a more specific, detailed description of the operation you want to perform?

sherly replied on : 42 of 119

I am working with counting of cells in an image. I have segmented the cells but i need your help to count the number of cells

deepak replied on : 43 of 119

count the no: of pixels in the each segmented part

Steve replied on : 44 of 119

Sherly—If you have already segmented the cells, then use bwlabel to count them.

Administrator replied on : 46 of 119

Arindam—I am not familiar with cardiac image segmentation.

Asanka replied on : 47 of 119

Is there any way to eliminate closely spaced markers of the wathershed segmentation inorder to avoid over segmentation?

Any suggestions to avoid over segmentation?

Steve replied on : 48 of 119

Asanka—If you are generating your markers using the minima of something like a gradient image, you might try using something like imextendedmin in order to filter out shallow minima.

TMan replied on : 49 of 119

I’m trying to calculate the number distinct objects in a huge (about 6000×6000) binary image B.

[ignore count] = bwlabel(B);

won’t work because of the double precision ‘ignore’-image not fitting my memory.

count = sum(sum(bwmorph(B,’shrink’,Inf)))

is, on the other hand, way too slow.

Anyone having efficient solutions to this problem?

silvia replied on : 50 of 119

hi steve,

I want to measure the area of cell nucleus. Could you explain to me the steps to do that?

Steve replied on : 51 of 119

TMan—You can start with bwlabel and modify it to suit your needs. Look in the file for the line:

numComponents = length(sizes);


That’s the answer you are looking for, and that line occurs before the memory-consuming output label matrix is constructed.

Steve replied on : 52 of 119

Silvia—Once you have the nuclei segmented, you can use bwlabel and regionprops to measure the area. I can’t help you much with the segmentation problem. Image segmentation usually requires some custom algorithm development that depends on the specific characteristics of your data set. Consult texts on microscopy and image processing for some general methods that you might be able to adapt to your data.

JP replied on : 53 of 119

When I copy this image, save it as a PNG, and try to use it, I get an error. Here is my code:

——————–
%test work with imfill
I_cropped = I;
imshow(I_eq)
——————–

All I’m doing is reading in the image, which I copied from the webpage. THen I get this error:

——————-
Error in ==> imfill_test at 4

??? Function ADAPTHISTEQ expected its first input, I, to be two-dimensional.
—————–

I can use adaphisteq if I read in a GIF or TIF with [X,map]=imread… and then use ind2gray, but the example image on this blog will not work for me at all, evne when saved as a TIF, using the same code that reads another TIF. I’m using R2006b. There must be something I’m missing about the definition of the image or something.

JP replied on : 54 of 119

Ok, if anybody else had the same problem as me, just make sure you grab only the top 2-D layer of the PNG after you read it in. On my computer it came in as MxNx3 and made adapthisteq choke. When I said I_cropped=I(:,:,1); it picked up a 2-D matrix and was happy with that input.

Steve replied on : 55 of 119

JP—The second paragraph of this post has a link to the original nuclei.png image for use with this image. The images you see on this web page are automatically captured screen shots from MATLAB, and they are stored in a color format even if the originals are grayscale.

I have modified the first line of code (the call to imread) to show how to read the correct original image directly, using its URL. I hope this will help others avoid similar confusion.

JP replied on : 56 of 119

Thanks, Steve! Pesky color images that aren’t color…

I will keep working to adapt the method to my images. In a nutshell, I have two images of a fuel spray taken with different filters. Each image is zeroed by a dark image, so it is nearly zero outside of the spray itself. The images are stored as matrices and then divided pixel by pixel. Within the spray itself, the resulting data is fine. But outside of the spray, you wind up dividing small noise by small noise and you get pure trash outside of the spray. I’m trying to use the same techniques to outline my spray and make a mask so I can discard everything outside the spray. I’m getting there, but I keep getting tripped up in the difference between matrices, greyscale images and color images. I’ll get it soon. Thanks!

Administrator replied on : 57 of 119

JP—Sounds like an interesting application.

Rudi replied on : 58 of 119

Steve:

Can I use this concept to segment human heads in a crowd scene? But the occlusion problem could be a significant issue for this purpose.

Thanks

Steve replied on : 59 of 119

Rudi—The concept of marker-controlled watershed segmentation is quite general, and might apply to your problem. I would expect the details to be very difficult to work out, however.

Chris replied on : 60 of 119

hi steve,

thanks for this suggestion of the cell segmentation. i have got a question: have you ever implement the pyramid linking algorithm from J. Burt? i am on a similar problem like this cell segmentation. but we’ve got colonies and sometimes the edges are blurred and sometimes they are sharp! if you want i can send you some pictures of these colonies!?!

Thanks
Chris

Steve replied on : 61 of 119

Chris—No, I’m not familiar with that method.

Rupesh Tatiya replied on : 62 of 119

Steve,

I have a similar image segmentation problem at hand where we have the RGB color image of human face and the objective is to segment the eyebrow from the rest of the face. I have tried different algorithms given in various tutorials but they do not work with good enough. So, can you help me in this matter?

Thanks and Regards,

Rupesh Tatiya

Steve replied on : 63 of 119

Rupesh—No, I generally do not have the time available to provide custom algorithm development advice for particular datasets.

Shannon replied on : 64 of 119

i am not sure if this has already been addressed/answered. I used code like yours, and incorporated what you did for segmentation.

I need to create code to count the specific cells once thresholded. like that produced a single number of the cell count, like 18.

i think i need a for loop, but i’m not really strong in matlab, so i was wondering if you knew how to count pixels in order to cell count or what… maybe had sample code or language to steer me in the right direction? THANK YOU!

shannon

Pradeep replied on : 65 of 119

Hi Steve
I am planning to do something on cell segmentation. But I am not an expert. I need your help. I need in areas of calculating more accurate cell density by keeping in mind: Partial Cell, Non-cell area and overlap cell. One more area i want to focus on is auto threshold if you have any ideas regarding this plz let me knw

Thanks

Steve replied on : 66 of 119

Shannon—Try this:

[L, num_cells] = bwlabel(segmented_cell_image);

Steve replied on : 67 of 119

Pradeep—Are your cells reasonably circular? If so, you might consider trying a circular Hough transform. There are many different variations on the literature. There’s at least one available on the MATLAB Central File Exchange.

Pradeep replied on : 68 of 119

thnx steve
i have one more question. why we only use black and white image for cell segmentation. like if we have any color image we first convert it to a grayscale.

Steve replied on : 69 of 119

Pradeep—It is convenient to use a binary (black and white) image to represent the outcome of the fundamental segmentation idea: Choosing which pixels belong to objects of interest, and which pixels belong to the background.

Sam replied on : 70 of 119

Hi ,
How can I identify objects with intermediate intensities in a image(there are bright objects, medium bright objects and the background), when the region around the brighter objects(glow around the bright objects) are of similar intensities to the medium intensity objects and needs to be excluded !

Tilman Sauter replied on : 71 of 119

Dear Sir,

I have developed a 3D image deconvolution code using the approximate method described in “Digital Image Processing, Castleman”. It uses the out of focus microscope images and the out of focus PSF images to determine the ‘noise’, which is then subtracted from the in-focus micrsocope image.

Now I would like to use for example the deconvlucy() function provided by Matlab to do the same thing in order to compare both.

Is it possible to work with out of focus planes in any of matlabs deconvolution functions to achieve 3D image deconvolution?

Tilman

Steve replied on : 72 of 119

Tilman—I don’t have the Castleman reference handy. Can elaborate on what you mean by “work with out of focus planes”?

Tilman Sauter replied on : 73 of 119

Dear Sir,

Thank you for your fast reply. The code I implemented is based on a simple equation:
in-focus-object = in-focus-image – sum(lower-out-of-focus-image convoluted with lower-out-of-focus-PSF + higher-out-of-focus-image convoluted with higher out-of-focus-PSF).
The sum sign defines the number of adjacent planes used, image the obtained image from the microscope and object the deconvoluted image.

The important thing is, that I can use information from the adjacent planes to improve the in-focus-object, whereas with deconvlucy() I can only use the information of the in-focus-image and in-focus-PSF.

Do you think it is possible to use deconvlucy() or another deconvolution function to achieve a comparable result?

With best regards
Tilman

Steve replied on : 74 of 119

Tilman—No, none of the deconvolution methods in the Image Processing Toolbox implement that technique. You would need to modify them yourself.

Steve replied on : 75 of 119

Sam—It’s hard to say without seeing a sample. Sounds like you have a difficult segmentation problem. If the object boundaries are reasonably well-defined, it might be helpful to use the gradient magnitude image.

Steve replied on : 76 of 119

Jon—There is a MATLAB Central search box at the top of the page. You can limit the scope of the search to MATLAB Central blogs. Or you can click on the “Blog archive” link in the side panel. I haven’t written that much about texture analysis, though. You might want to try the texture analysis section in the toolbox documentation.

Jon Hauris replied on : 77 of 119

Steve, I am interested in texture processing. Is there a way to search your excellent archives to see if you have covered this area and where it is located. Thank you, Jon

Giles Kingsley replied on : 78 of 119

Steve,

I have an image of a metal plate from a time-lapse series. There is a grid in the background, which I want to remove. Does this sound like a morphological operation? Throw me a bone here. Edge detection just enhances the grid, and thresholding removes area(my output) from the plate. Any ideas?
-Giles

Steve replied on : 79 of 119

Giles—An opening (IMOPEN) might work. How thick is the grid, and how thick is the object(s) of interest?

IB replied on : 80 of 119

I get the following error when I follow the commands here
“Undefined function or method ‘imoverlay’ for input arguments of type ‘uint8’.”
How could I overcome this error?

Steve replied on : 81 of 119

IB—The third paragraph in this blog post tells you how to get the function imoverlay.

IB replied on : 82 of 119

Steve do you mean the cell segmentation blog? I did not found it in the cell segmentation blog

Steve replied on : 83 of 119

IB—In this blog post, the very one you are commenting on, the third paragraph says “Note: you can download the functions imcredit and imoverlay from MATLAB Central.” The links are provided there.

azhagumani replied on : 84 of 119

i want to expand a region and find out what are all label values it reached.

Steve replied on : 85 of 119

Azhagumani—Form a binary image whose foreground is the region you are interested in, and then dilate it. Then use the dilated image as a logical mask into the original labeled image; that will pull out all the labels adjacent to the original labeled region.

huda replied on : 86 of 119

Hye.. I’m new to this image processing things.. I was told to make an image segmentation of a picture.. The image consists of 3 shapes that overlapping each other.. It was a grayscale image but with 3 different tones of gray colour.. The task was to convert to binary format first and then segment the shape individually.. If during the segmentation, the shape missing some portion, I have to complete the missing portion.. So, can you tell me how to do it in Matlab? And I wonder do I have to set a threshold value?

Steve replied on : 87 of 119

Huda—It’s difficult in general to reconstruct overlapping shapes in a segmentation problem. It’s easier for some specific shapes, though. For example, if your shapes are circles, you could use a circular Hough transform. Hough-like algorithms can often detect a partially occluded shape. The Image Processing Toolbox does not have a circular Hough algorithm, but there’s at least one available on the MATLAB Central File Exchange.

huda replied on : 88 of 119

The shapes are circle, square and rectangle.. Can i used the Hough transform for these shapes? And for the segmentation? How to do it?

Steve replied on : 89 of 119

Huda—The problem sounds kind of contrived. I don’t have any particular suggestions.

Chris Soton replied on : 90 of 119

Hi Steve,

I am MATLAB beginner and I am using a similar image to this example. I have converted my nuclei image to binary and used bwlabel to count the nuclei. However I would like to measure the distance between cell nuclei from a central point on each nuclei (centroid), but I am struggling with the code. Is there any coding help or pointers you could give me with this please?

Thanks,

Chris

Steve replied on : 91 of 119

Chris—You can pass the output of bwlabel to regionprops to get the centroid of each labeled object.

ragu replied on : 92 of 119

hi, steve
i want to develop a general algorithm for image resulting from watershed to display how many pixel occupied in each segment and average intensity of each segment .so that i need variables to sum up the the pixel and average intensity. for limited segment i can develop with specific variables but for generalize i need to get variables which should generate depends on how many segments from watershed output.

Chris Soton replied on : 93 of 119

Hi Steve,

Thanks for the reply, I am sorry but I think I was not clear. I have obtained the centroids of my cells and am struggling to get the distance between the centriods. I would like the code to automatically get the distance between all centriods.

I have tried imdistline but this seems to be a manual code and my attempts to pass the centriod data to this code have failed. Please could you give me any sample code or pointers?

Thanks,

Chris

Steve replied on : 94 of 119

Ragu—You can pass the output from watershed to the function regionprops in order to calculate the desired measurements for each labeled region.

Steve replied on : 95 of 119

Chris—Try something like this:

dist = hypot(centroid1 - centroid2);

Stefan replied on : 96 of 119

Hello Steve!

First of all I’d like to say that I very appreciate your Help here in this Blog. You’re amazing!
I have a similar problem like those discussed here. I have an Image where you can find lots of small cruxes (+), which are all different orientated. I would like to get a contour on every single crux, like you did it on the cells on top of this site. My problem is that i cant isolate the cruxes so i always get connected regions.

The image is here: stud.unileoben.ac.at/~m0335168/RC7.jpg

Thanks a lot!

Steve replied on : 97 of 119

Stefan—I’m sorry, but I don’t have specific suggestions for you.

Chris Soton replied on : 98 of 119

Thanks Steve! Your code combined with another of your examples has helped me greatly.

Ragu replied on : 99 of 119

hi steve
right now i dont have imoverlay function in 6.5.1 version

can you give any function to do the same

or how to develop algorithm

Steve replied on : 100 of 119

kavitha replied on : 101 of 119

Hi Steve.
Its a great page. I got answres for a lot of questions. and in my project i have two images for a same object. I ve segmented those two objects. but i need to calculate the area difference between those objects. how can i do this. pls help me.

Steve replied on : 102 of 119

Kavitha—Use bwlabel and regionprops.

kavitha replied on : 103 of 119

k. steve…. thank u for ur reply. i refered the page which u wrote abt bwlabel.Ur work and discussion are very useful to students like us.. Those examples helped me a lot…

Ragu replied on : 104 of 119

hi steve thanks for the reply

now i have two questions .

1.how many ways are there to make markers.

2.i need to make marker on my gray image before apply to watershed in order to segment the defects(black in color).

where should i make marker either on black parts or gray parts.

Steve replied on : 105 of 119

Ragu—Take a look at the marker-controlled watershed segmentation demo in the Image Processing Toolbox (you can also find it on our web site). You might also be interested in my watershed article in News and Notes from a few years ago.

ragu replied on : 106 of 119

in water marking
fore ground marking is based on the imregionalmin or imregionalmax followed by imclose ,imfill, bwareaopen.

in back ground marking i have little confussion of super imposing fore ground and back ground marker images on gradiant magnitude image rather than original gray image.

Steve replied on : 107 of 119

Ragu—I’m not sure exactly what your question is. Whether to use the original image, the gradient magnitude image, or some other derived image depends on your data set and your chosen method for extracting markers.

ragu replied on : 108 of 119

hi steve thanks for the question,
i have made foreground markers and background markers as per my region of interest.

now should i impose these on original image or gradiant magnitude image.

in example of matlab watermarking it mentioned magnitude gradient.
why?

i belive that imposing markers on original gray image is logically good. if it is not .why?

ragu replied on : 109 of 119

im sorry i forgot to mention that ,my ultimate aim to get segmentation of defects on gray image which is in dark portion ,by applying watershed transform.

for this i need markers to reduce oversegmentation, .

Steve replied on : 110 of 119

Ragu—In my comment #105, I gave a link to my News & Notes article on the watershed transform, which explains the basic concepts of using the transform for segmentation. I encourage you to take a look at it. Here’s an important sentence: “The key behind using the watershed transform for segmentation is this: Change your image into another image whose catchment basins are the objects you want to identify.” Often, the original grayscale image does not have the property that its catchment basins correspond to the objects you are trying to segment. Markers, minima imposition, computing gradient magnitudes, etc., are all techniques for making that happen.

Pratyusha replied on : 111 of 119

HI Steve,
I have a similar problem with segmentation of tumors in an Ultrasound image.
I use a .jpg image file and edge detection techniques, but most of the times it does not work, as the background of an US image is not that clear, i.e., there is a very smooth gradient. I could send/upload some of the files (Images) so you can have a look at them.
Please suggest a suitable edge detection technique for these kind of images.

Thank you!

Steve replied on : 112 of 119

Pratyusha—I don’t have any suggestions for you.

Pratyusha replied on : 113 of 119

Thank you for letting me know. I will try to use a kalman filter for it, or color code the image to highlight the tumor areas.

Francesco replied on : 114 of 119

Goodmorning,
in my engineering degree study i have to make a microstructural characterization of open cell aluminium foams.
After i have done a watershed transform on my 2D slice i have to do a masking of the complement of the binary starting image and the watershed one to reduce the oversegmentation but i don’t know how to make this.

Steve replied on : 115 of 119

Francesco—This blog post has some example code. You might also want to look at the watershed segmentation demo in the Image Processing Toolbox.

Francesco replied on : 116 of 119

thresh=40;
BW=(I>40);
D=bwdist(I1,’euclidean’);
D1=imcomplement(D1);
H=imhmin(D1,13,6);
L=watershed(H);
d=imcomplement(BW)&L;
Now i have reconstructed my open cell foam but if i want to make a geometrical analysis it seems to be that the image has only one region instead of hundreds,because matlab doesn’t recognize that there are many regions.
How i can make to solve this problem??

Steve replied on : 117 of 119

Francesco—It’s your algorithm that doesn’t recognize that there are many regions. It’s hard to say for sure without looking at your image, but I’m not sure I understand why all your steps make sense. You might want to take a look at my News & Notes article about watershed-based segmentation.

deXter Holland replied on : 118 of 119

http://img3.imageshack.us/my.php?image=sample01.jpg

hi all! i need to do cell detection on the image above (identify the presence of those ‘big’ cells). my main problem is the background (noise), i couldn’t remove them properly, so i couldn’t do the usual segmentation easily. for now i’m processing it on grayscale, based on the materials i’m reading.
any suggestion?

Steve replied on : 119 of 119

Dexter—You might try using bwareaopen to remove the small “noise” objects detected.