Region Growing

Segmentation by growing a region from seed point using intensity mean measure

Description

Simple but effective example of "Region Growing" from a single seed point.

The region is iteratively grown by comparing all unallocated neighbouring pixels to the region. The difference between a pixel's intensity value and the region's mean, is used as a measure of similarity. The pixel with the smallest difference measured this way is allocated to the region.

This process stops when the intensity difference between region mean and new pixel becomes larger than a certain treshold

• Acknowledgements: This submission has inspired the following:; Region Growing (2D/3D grayscale)
• MATLAB release: MATLAB 7.5 (R2007b)

[Xin Liu]

It's very helpful. Thanks!

[jichao zhao]

That is excellent one. I am wondering whether you can extend your work to colorful graph, not just gray scale.

[balamurugan samy]

useful

[Roberto Mantia]

It does what it should do. Very good !

[Zonglei Zhen]

region growing is a very useful method for sgementation. This program is a good start of it.

[Paolo L]

very good, simple to read, simple to understand, it works fine for a start, useful for student exercises

[ryan hartadi]

it works very good, so simple. need to de developed in other area

[rasool shaik]

good

[Ida Haggstrom]

 

[ucd puri]

Hi All, This is an excellent piece of code. I am just an intermediate matlab programmer and have a couple of questions on this code.
(1) It works perect on normal images. But i have a dicom int16 medical image, when i give it as an input, i get an output image with only one pixel (the seed pixel that i selected) . I changes the code a little bit for example instad of imread, i used dicomread, and in place of imshow(I,[]) i used imshow(I,'displayRange',[]).
(2) the file explaination reads '"This process stops when the intensity difference between region mean and
 new pixel become larger than a certain treshold (t)" but this t is 0.2 which you have defined in the file, what if i want to use an grey intensity level of say 600 as threshold. As my image is int16 dicom, it has 2^16 , so it grey sclae ranges from -32768 to + 32768.
thanks a lot for help in advance.
tanuj

[ucd puri]

Hi, I figured out what was the problem, my data set was in int16 format so, it was eaisy to make amendments, but i am not sure if my region is segmented correctly. Regards tanuj

[Jag]

thanks, seems like a handy piece of code.
On R2008b I had to change line 70 to

dist = abs(neg_list(1:neg_pos,3)-double(reg_mean));

to get it to run with the sample image.

[Torres Chhetri]

its was a great help, would you like to extend your for taking random seed pixels for a colorful image. if u do that, it would be so great.

[siti 85]

cannot download.. why?

[jack lumber]

 

[taccini]

Hello All,
first i want to thank Dirk about his work and all the matlab community because there are very helpful for me.
well i need to use the region growing for 3D matrix.
this code could be adapted or any one have some idea?
thanks.

[anoop R]

Very useful piece of code.

[Natarajan]

code works except as Ron did . got to change the line 70.

However the i didn't get a segmented/thresholded output for 16 bit grey scale! Thanks

Would be great if someone could give an idea of whats the issue!

[Jess]

No need to change line 70 if use im2double with imread as given in the example (line 21 - comment section)

[Esther Kneepkens]

Thank you a lot!

[Hank Will]

This looks pretty cool! But I wonder processing larger images with this routine (e.g. satellite data). Any suggestions? In my case I don't get any error, but processing a 10 Bit panchromatic imagery subset of just 15 Mb leads to no valuable results.
Would highly appreciate any thoughts about workflows for matlab for this purpose!
Thanks!

[Aravind]

can somebody explain how to segment a carotid artery ?

[myriam robert]

I have a MRI coup from a 3D data matrix "Y" so it's Y( : , : , 14) .
the algorithm doesn't work on this coup. In the result I have only one white pixel ( the seedpoint).
I change the line 70 but it doesn't change anything.
can you help me please?

[myriam robert]

I change reg_maxdist to 10 and it's ok.
Does anybody have find how use the region growing for 3D matrix ?
thank you

[paramveer sran]

i have downloaded this piece o code but when i tried to run on the same image i.e medtest.png i got the following error. could you plz tell me where is the problem and how to solve this
>> regiongrowing
??? Maximum recursion limit of 500 reached. Use set(0,'RecursionLimit',N)
to change the limit. Be aware that exceeding your available stack space can
crash MATLAB and/or your computer.

Error in ==> isstruct at 99
thanx

[michiko]

very helpful!Thank you very much!!!

[Xin He]

very helpful thanks!

[saif]

thank you

[kalarmago]

What are x and y, I believe image(x,y) coordinate, but when I tested other x,y values (eg, x=308,y=244, I(x,y)=0.73) then the segmented result image is not desired.

[fsdm fsdm]

 

[Sivan Ran]

Hi,
How can I make it work for 3D images?

Thank you

[Deepak]

Thq Dirk-Jan Kroon for such short code.I am working dicom image.I convert it into JPG.i have a problem over run time.

??? Undefined function or method 'regiongrowing' for input
arguments of type 'double'.

I even defined im2double() fnc.but i got an error...