data-raw/Olsen-mtl/load_olsen.R

In bbuchsbaum/neuroatlas: useful atlases/parcellations for neuroimaging analysis

setwd("data-raw/Olsen-mtl/")
library(neuroim2)

vol <- read_vol("AVG_MNI_MTL_P33.nii")

full_label <- scan("labels.txt", "")
hemi <- substr(full_label, 1,1)
hemi <- ifelse(hemi == "L", "left", "right")
label <- sapply(full_label, function(x) substr(x, 3, nchar(x)))
names(label) = NULL

#names(labels) <- c("ROINUM", "label", "red", "green", "blue")
paste(as.vector(col2rgb(x)), collapse = " ")
cmap <- rbind(t(col2rgb(cols)), t(col2rgb(cols)))

olsen_mtl <- list(
  name=paste0("Olsen_MTL"),
  atlas=vol,
  cmap=cmap,
  ids=1:length(full_label),
  labels=label,
  orig_labels=full_label,
  hemi=hemi)

class(olsen_mtl) <- c("olsen_mtl", "atlas")
usethis::use_data(olsen_mtl, internal = FALSE, overwrite = TRUE, compress = "xz")