EndPoints in 3D skeleton

9 views (last 30 days)
Hege
Hege on 6 Apr 2021
Commented: Hege on 12 Apr 2021
I have a 3D Skeletonized Micro-CT Fiber model which obtains completely seperated 7 volumes(fibers). I would like to find the endpoints of each and every fiber and I used below commands to find the end points.
B1=bwskel(testvol,'Minbranchlength',20);
BP = bwmorph3(B1,'endpoints');
Endpoints_Cor=regionprops3(BP,'VoxelList');
stats4=regionprops3(B1,'Volume','VoxelList');
However, my model obtains 7 volumes which are completely seperated. But it gives only 13 endpoints(But I am expecting 14 end points.7*2=14). I noticed one of the endpoint of my third volume does not obtain in the end point list. only one end point is in the list. Could you please explain how I should obtain all the endpoints in each and every volume and why I am missing one endpoint of the third volume? And also, is there anyway to categorise which endpoints belongs to which volume?
Your answers are highly appreciated!

Accepted Answer

Image Analyst
Image Analyst on 7 Apr 2021
bwmorph() has an 'endpoints' option. The image should be a skeleton. bwmorph() also has a skeleton option but I think bwskel() has a better algorithm.
  5 Comments
Image Analyst
Image Analyst on 8 Apr 2021
You can do a brute force loop to find the distance of every point to every other point.
Hege
Hege on 12 Apr 2021
Thank you so much Image Analyst! I really appreciate your answer.
What I did was, found the minimum distance for each and every voxel and rearranged the regionporops3() data as I want using a one ending point of each volume.
Thank you so much!

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!