how to calculate the diameter of a vessel

hello everyone

I really need help as soon as possible

I want a solution to calculate the diameter of each white line in the image(Vessel Diameter Measurement)

[John D'Errico]

Compute the radius. Multiply by 2.

Really, we are not here to do your homework as soon as possible. IF you could not start it in a timely way, then why should we rush?

Make an effort, as that is the way to learn. Or don't make an effort. In that case, learn to say the phrase "Do ya want fries with that order?"

Answer 1

[Sean de Wolski]

You might find the second part of this webinar helpful:

http://www.mathworks.com/videos/medical-image-processing-with-matlab-81890.html

Answer 2

[Image Analyst]

It's pretty easy - maybe 3 lines of code. First call bwdist to get the distance from the white pixels to the nearest black pixel. Then call bwmorph to skeletonize that so you get the distances along the centerlines. That's the radius so multiply by two to get the diameter. It will be an image looking like a skeleton and to find the diameter you simply look at the value of the pixel at the location you want the diameter of. Basically it's this (untested):

edtImage = 2 * bwdist(~binaryImage);
skeletonImage = bwmorph(binaryImage, 'skel', inf);
diameterImage = edtImage .* double(skeletonImage);

Can't get much simpler, though you may want to add a utility where you let the user call ginput() to report back the diameter at the closest pixel to where they clicked. It's pretty easy, though clever and not obvious (until you know about it), so most novices might embark on some complicated looping process to get cross sections (what a nightmare that would be).

[Image Analyst]

bil - are you still around? Did it work?

[bil bbil]

Sorry for the delay,,I followed the example in medical-image-processing,,I had it calculated the bifurcation points

[Lulu Firdaus]

i got this error.

[Image Analyst]

Try

diameterImage = edtImage .* double(skeletonImage);

[Rey Ezra Langelo]

in addition for anyone who cant see the diameter result, try adding this at the end of your code:

mesh(diameterImage);

so you will get 3D graph, with X and Y as the skeleton's pixel location and Z as the vessel's diameter

[Lulu Firdaus]

it worked, thank you!

[Hasnain Iftikhar]

did you extract thin vessels?

[Image Analyst]

@Hasnain Iftikhar, look here for possible algorithms:

http://www.visionbib.com/bibliography/contentsmedical.html#Medical%20Applications,%20CAT,%20MRI,%20Ultrasound,%20Heart%20Models,%20Brain%20Models

Answer 3

[Nitin]

I would start by first experimenting with regionprops function in matlab

[bil bbil]

you have an idea how to fit in my case??