Skel2Graph 3D
MATLAB code to derive the network graph of a 3D voxel skeleton
This function converts a 3D binary voxel skeleton into a network graph described by nodes and edges. The input is a 3D binary image containing a one-dimensional voxel skeleton, generated e.g. using the "Skeleton3D" thinning function available on MFEX. The output is the adjacency matrix of the graph, and the nodes and links of the network as MATLAB structure.
Usage:
[A,node,link] = Skel2Graph(skel,THR)
where "skel" is the input 3D binary image, and "THR" is a threshold for the minimum length of branches, to filter out skeletonization artifacts.
A is the adjacency matrix with the length of the links as matrix entries, and node/link are the structures describing node and link properties.
A node has the following properties:
- idx List of voxel indices of this node
- links List of links connected to this node
- conn List of destinations of links of this node
- comX,comY,comZ Center of mass of all voxels of this node
- ep 1 if node is endpoint (degree 1), 0 otherwise
A link has the following properties:
- n1 Node where link starts
- n2 Node where link ends
- point List of voxel indices of this link
A second function, "Graph2Skel3D.m", converts the network graph back into a cleaned-up voxel skeleton image.
An example of how to use these functions is given in the script "Test_Skel2Graph3D.m", including a test image. In this example, it is also demonstrated how to iteratively combine both conversion functions in order to obtain a completely cleaned skeleton graph.
Any comments, corrections or suggestions are highly welcome. If you include this in your own work, please cite our publicaton [1].
Philip Kollmannsberger 09/2013, 01/2016
philipk@gmx.net
[1] Kollmannsberger, Kerschnitzki et al.,
"The small world of osteocytes: connectomics of the lacuno-canalicular network in bone."
New Journal of Physics 19:073019, 2017.
Citation pour cette source
Kollmannsberger, Kerschnitzki et al., "The small world of osteocytes: connectomics of the lacuno-canalicular network in bone." New Journal of Physics 19:073019, 2017.
Compatibilité avec les versions de MATLAB
Plateformes compatibles
Windows macOS LinuxCatégories
- MATLAB > Graphics > 2-D and 3-D Plots >
Tags
Remerciements
Inspiré par : Skeleton3D
A inspiré : Calculate the network graph of a 3D voxel skeleton.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Découvrir Live Editor
Créez des scripts avec du code, des résultats et du texte formaté dans un même document exécutable.
Les versions qui utilisent la branche GitHub par défaut ne peuvent pas être téléchargées
Version | Publié le | Notes de version | |
---|---|---|---|
1.22.0.1 | updated description |
|
|
1.22.0.0 | Linked to Skeleton3D |
|
|
1.21.0.0 | Linked to github |
||
1.2.0.0 | Fixed typo
|
||
1.1.0.0 | Function returned an error if the skeleton volume was not padded with zeros. This is now corrected. |
||
1.0.0.0 |