MATLAB Answers

How can I calculate the thickness for the attached images?

80 views (last 30 days)
Alaa
Alaa on 4 Mar 2015
Commented: Image Analyst on 8 Mar 2015
How can I calculate the thickness for the attached image??
I did segmentation for all the nerves in the image and I labelled each object in the image. What I want is how can I calculate the thickness of each nerve? And how can I highlight a specific part of the nerve if there is an inflation in any nerve?
Your help in this regard is highly appreciated, thank you in advance

  0 Comments

Sign in to comment.

Accepted Answer

Image Analyst
Image Analyst on 4 Mar 2015
Use regionprops to measure area and perimeter. Then assume blobs are long and thin. So area = length*width and perimeter = 2*w+2*l = 2*w+area/w. So
w*P = 2*w^2 + A
or
0 = 2*w^2-w*P + A
Use quadratic formula to solve
w = (P - sqrt(P^2-4*2*A))/(2*2)
w = (P-sqrt(P^2-8*A))/4

  5 Comments

Show 2 older comments
Alaa
Alaa on 8 Mar 2015
Dear Image analysis
I code your answer and it work perfectly, but I don't know how can set a specific color of each region. is there any idea??
In fact, I have the attached image, I completed the segmentation for all the nerves and I extract all the needed information from each nerve. My finall goal is if there is inflation in any nerve how can I detect it and make a final image map where the normal part of the nerve it should be green and the abnormal part let say red color, I used your method and it worked ok, but is there and suggestion how can I set different colors for the both parts either in the grayscale image or binary image?
T
hank you in advance
Image Analyst
Image Analyst on 8 Mar 2015
First you need to get a good binary image like you had before. Maybe you can just threshold it, or maybe use a Frangi filter or maybe use anisotropic diffusion (demo attached).
Then you need to use bwdist() to get the thickness of each place in the nerves. But the thickness will only be along the centerline of the distance transform and you want the entire thickness to be represented by some color, not just the centerline so you need to dilate out the centerline to that it takes up the whole binary image. To do this, use imdilate() which is a local max filter. Dilate it out a lot, like twice as big as the nerve. Then mask the dilated image by the original binary image to crop it to just the nerve boundaries. Now you have a solid color all the way out to the edge and you can simply apply a colormap. So the steps are
  1. process and get a binary image
  2. call bwdist
  3. call imdilate
  4. mask: out = dilatedImage .* binaryImage
  5. colorize: colormap().
If you have trouble, post your code.

Sign in to comment.

More Answers (0)

Categories

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!