Color My Brain: Part 5

I thought this one was too good to not share.

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

Using SLURM

One of the challenges to learning neuroimaging is the large amount of coding necessary to do it right. Many (like myself) begin with absolutely no coding experience and find themselves lost most of the time. Once I learned how to process a single scan I was pretty proud. Then I was told that I had to process scans in parallel on the supercomputer. I was less proud when I spent two weeks trying to figure that out. Here is a painfully detailed description of the first part of every batch script and some hints for using the supercomputer in general. This is pretty specific to my experiences in the brain imaging and behavior lab (BIABL) at BYU, but may have applicability elsewhere too. Every batch file needs to start with this. Sometimes I’ll copy and paste a large batch script into terminal and miss this part. The batch process will fail and I won’t have any logfiles to figure out the errors. Usually it’s because I missed this....
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