R/eigSeg.R

In stnava/itkImageR: Examples using Rcpp to interface R and ITK

eigSeg <- function(mask = NA, imgList = NA) {
  if (typeof(mask) != "S4") {
    print(args(eigSeg))
    return(1)
  }
  maskvox <- (mask > 0)
  maskseg <- antsImageClone(mask)
  maskseg[maskvox] <- 0
  if (length(imgList) > 0) {
    if (typeof(imgList) == "list") 
      mydata <- imageListToMatrix(imgList, mask)
    if (typeof(imgList) == "character") 
      mydata <- imagesToMatrix(imgList, mask)
    segids <- apply(abs(mydata), 2, which.max)
    segmax <- apply(abs(mydata), 2, max)
    maskseg[maskvox] <- (segids * (segmax > 1e-09))
    print(max(segmax))
    print(max(segids))
    return(maskseg)
  } else print("No images in list")
}