data-raw/create-atlas.R

In ggseg/ggsegFlechsig: Flechsig datasets for the ggseg-plotting tool

library(ggsegExtra)
library(ggseg)
library(ggseg3d)
library(dplyr) # for cleaning the atlas data efficiently
library(tidyr) # for cleaning the atlas data efficiently

# The unique name of the atlas annot, without hemisphere in filename
annot_name <- "flechsig"


for(hemi in c("lh", "rh")){
  tmp <- here::here(sprintf("data-raw/flechsig/%s.fsaverage_Flechsig", hemi))
  freesurfer::mris_convert_annot(
    infile = file.path(freesurfer::fs_subj_dir(), "fsaverage", "surf", "rh.orig"),
    annot = sprintf("%s.label.gii", tmp),
    ext = ".annot",
    outfile = sprintf("%s.annot", tmp)
  )
}

# bash
# sudo cp data-raw/flechsig/lh.fsaverage_flechsig.annot $FREESURFER_HOME/subjects/fsaverage/label/lh.flechsig.annot
# sudo cp data-raw/flechsig/rh.fsaverage_flechsig.annot $FREESURFER_HOME/subjects/fsaverage/label/rh.flechsig.annot



# You might need to convert the annotation file
# convert atlas to fsaverage5
lapply(c("lh", "rh"),
       function(x){
         mri_surf2surf_rereg(subject = "fsaverage",
                             annot = annot_name,
                             hemi = x,
                             output_dir = here::here("data-raw/fsaverage5/"))
       })


# Make  3d ----
flechsig_3d <- make_aparc_2_3datlas(
  annot = annot_name,
  annot_dir = here::here("data-raw/fsaverage5/"),
  output_dir = here::here("data-raw/")
)
ggseg3d(atlas = flechsig_3d)

## fix atlas ----
# you might need to do some alteration of the atlas data,
# like cleaning up the region names so they do not contain
# hemisphere information, and any unknown region should be NA
flechsig_n <- flechsig_3d
flechsig_n <- unnest(flechsig_n, ggseg_3d)
flechsig_n <- mutate(flechsig_n,
                    region = ifelse(grepl("Unknown|\\?", region, ignore.case = TRUE),
                                    NA, region),
                    atlas = "flechsig_3d"
)
flechsig_3d <- as_ggseg3d_atlas(flechsig_n)
ggseg3d(atlas  = flechsig_3d)
ggseg3d(atlas  = flechsig_3d, hemisphere = "left")


# Make palette ----
brain_pals <- make_palette_ggseg(flechsig_3d)
usethis::use_data(brain_pals, internal = TRUE, overwrite = TRUE)
devtools::load_all(".")


# Make 2d polygon ----
flechsig <- make_ggseg3d_2_ggseg(flechsig_3d,
                                 ncores = 15,
                                 steps = 6:7,
                                 tolerance = .1,
                                 output_dir = here::here("data-raw/"))

plot(flechsig)

flechsig %>%
  ggseg(atlas = ., show.legend = TRUE,
        colour = "black",
        mapping = aes(fill=region)) +
  scale_fill_brain("flechsig", package = "ggsegFlechsig", na.value = "black")


usethis::use_data(flechsig, flechsig_3d,
                  internal = FALSE,
                  overwrite = TRUE,
                  compress="xz")


# make hex ----
atlas <- flechsig

p <- ggseg(atlas = atlas,
           hemi = "left",
           view = "lateral",
           show.legend = FALSE,
           colour = "grey30",
           size = .2,
           mapping = aes(fill =  region)) +
  scale_fill_brain2(palette = atlas$palette) +
  theme_void() +
  hexSticker::theme_transparent()

lapply(c("png", "svg"), function(x){
  hexSticker::sticker(p,
                      package = "ggsegFlechsig",
                      filename = sprintf("man/figures/logo.%s", x),
                      s_y = 1.2,
                      s_x = 1,
                      s_width = 1.5,
                      s_height = 1.5,
                      p_family = "mono",
                      p_size = 10,
                      p_color = "grey30",
                      p_y = .6,
                      h_fill = "white",
                      h_color = "grey30"
  )

})

pkgdown::build_favicons(overwrite = TRUE)