freesurfer_lut.R

In neuroconductor/freesurfer: Wrapper Functions for 'Freesurfer'

read_lut = function(){ 
  lut_file = file.path(fs_dir(), "FreeSurferColorLUT.txt")
  lut = readLines(lut_file)
  lut = grep("^(#.*)?$", lut, 
             value = TRUE, invert = TRUE) 
  lut = lut[ !(lut %in% "") ]
  
  str_to_matrix <- function(strings) {
    strings = strsplit(strings, split = " +")
    strings = do.call("rbind", strings)
    return(strings)
  }  
  
  mat = str_to_matrix(lut)
  colnames(mat) = c("index", "label", 
                    "R", "G", "B", "A")
  mat = data.frame(mat, stringsAsFactors = FALSE)
  mat$index = as.numeric(mat$index)
  mat$R = as.numeric(mat$R)
  mat$G = as.numeric(mat$G)
  mat$B = as.numeric(mat$B)
  mat$A = as.numeric(mat$A)
  
  return(mat) 
}

fs_lut = read_lut()
fs_lut = dplyr::arrange(fs_lut, index)
stopifnot(!any(duplicated(fs_lut$index)))
save(fs_lut, file = "data/fs_lut.rda", compression_level = 9)