R/skelConnectivity.R

In AlessioVeneziano/indianaBones: Analysis of bone 3D morphology

#' Skeleton connectivity
#'
#' Calculates the connectivity of the topological skeleton as the number of branches connected to each node.
#'
#' @param skel topological skeleton, output of \code{readAmiraSkeleton}
#'
#' @return A vector of the number of branches connected to each node. The vector has length equal to the number of nodes in the topological skeleton.
#'
#' @author Alessio Veneziano, Fabio Alfieri
#'
#' @references
#' Veneziano A, Cazenave M, Alfieri F, Panetta D, Marchi D. 2021. Novel strategies for the characterization of cancellous bone morphology: Virtual isolation and analysis. American Journal of Physical Anthropology.
#'
#' @examples
#' #Calculate skeleton connectivity
#' data(exampleSkeleton)
#' conn<-skelConnectivity(exampleSkeleton)
#'
#' @export

skelConnectivity<-function(skel){
  skVC<-skel[[2]]
  skEC<-skel[[3]]
  nodeC<-as.vector(table(c(skEC)+1))

  return(nodeC)
}