Color My Brain: Part 2

Another picture of my brain poked and prodded until ready to color

Comments

Post a Comment

Popular posts from this blog

Making a Brain Template

Image
My last post was on how to obtain a representative sample to create a brain template. This post shows how to create the actual template. This step is fairly straightforward because most of the heavy lifting is already performed by antsMultivariateTemplateContstruction.sh . Here’s the entirety of the script so that you have an idea of how it all looks together. #!/bin/bash #SBATCH --time=75:00:00 #SBATCH --ntasks=1 #SBATCH --nodes=1 #SBATCH --mem-per-cpu=32768M #SBATCH -o /fslhome/username/logfiles/dataset/output_temp.txt #SBATCH -e /fslhome/username/logfiles/dataset/error_temp.txt #SBATCH -J "ucsdtemp" #SBATCH --mail-user=LongbranchPennywhistle@gmail.com #SBATCH --mail-type=BEGIN #SBATCH --mail-type=END #SBATCH --mail-type=FAIL export PBS_NODEFILE=`/fslapps/fslutils/generate_pbs_nodefile` export PBS_JOBID=$SLURM_JOB_ID export PBS_O_WORKDIR="$SLURM_SUBMIT_DIR" export PBS_QUEUE=batch export OMP_NUM_THREADS...
Read more

Image Manipulation and Statistics

The following represents the current extent to which the iClass package I’ve been developing has been tested. The goal is to impliment it in ANTsR eventually when I’m not in the midst of medical school classes. For now this can be used if ANTsR is already installed and the following code is run to download the package from github. ## install.packages("devtools") devtools::install_github("Tokazama/iClass") library(ANTsR) library(h5) home <- '/Volumes/SANDISK/datasets/ucsd/' ucsd <- read.csv(paste(home, 'spreadsheets/ucsdWOna.csv', sep = ""))[, -1] iGroup The first thing I want to do is create an object that represents my image information in a convenient way. I can do this using the iGroup class. The following demonstrates how I do so with whole brain morphometry images. wblist <- c() boolwb <- rep(FALSE, nrow(ucsd)) for (i in 1:nrow(ucsd)) { tmppath <- pas...
Read more

Visualizing Cortical Thickness In a Population

Image
Today my goal is to develop some methods for comparing an individual's neuroanatomy to a population's and make some graphs that help convey my findings. To demonstrate this I will be using the Pediatric Template of Brain Perfusion (PTBP), which has already been processed and can be found here . I'll start off by loading some libraries I'll need and the CSV file. The CSV file contain demographics data and anatomically specific information about the cortical thickness, fractional anistrophy, cerebral blood flow, and blood oxygenation level. For now I'm just going to focus on cortical thickness (found by the prefix “Thick*”). Picture of amygdala from recent paper on the effects of cannabis suppressPackageStartupMessages(library(dplyr)) suppressPackageStartupMessages(library(ggplot2)) home &lt;- '/Users/zach8769/' ptbp &lt;- read.csv(paste(home, 'Desktop/ptbp/data/ptbp_summary_demographics.csv', sep = ""), header = TRUE) thick_pt...
Read more