Difference between revisions of "Cortical Thickness"
(→Freesurfer Segmentation) |
|||
| (4 intermediate revisions by one other user not shown) | |||
| Line 2: | Line 2: | ||
Page describing how we analyze cortical thickness with freesurfer | Page describing how we analyze cortical thickness with freesurfer | ||
| + | |||
| + | ==Setting up Freesurfer and MATLAB== | ||
| + | First make sure that you have the correct lines in your .bashrc file to run freesurfer: | ||
| + | export FREESURFER_HOME=/usr/local/freesurfer | ||
| + | source $FREESURFER_HOME/SetUpFreeSurfer.sh > /dev/null | ||
| + | export SUBJECTS_DIR=/mnt/diskArray/projects/freesurfer | ||
| + | Next make sure that you have the right tool boxes in your matlab search path. This should be done through your startup.m file | ||
| + | addpath(genpath('~/git/yeatmanlab')); | ||
==Create a T1 weighted nifti image for the subject== | ==Create a T1 weighted nifti image for the subject== | ||
Step 1: In a terminal, convert the PAR/REC files to nifti images. You may not need to do this if you have already gone through the [[Anatomy_Pipeline|Anatomy Pipeline]] | Step 1: In a terminal, convert the PAR/REC files to nifti images. You may not need to do this if you have already gone through the [[Anatomy_Pipeline|Anatomy Pipeline]] | ||
| − | cd / | + | cd /mnt/diskArray/projects/MRI/[subid] |
parrec2nii -c -b *.PAR | parrec2nii -c -b *.PAR | ||
Step 2: In MATLAB compute the root mean squared (RMS) image. Once again this might have already been done in the [[Anatomy_Pipeline|Anatomy Pipeline]] so you can re-use that RMS image | Step 2: In MATLAB compute the root mean squared (RMS) image. Once again this might have already been done in the [[Anatomy_Pipeline|Anatomy Pipeline]] so you can re-use that RMS image | ||
| Line 13: | Line 21: | ||
==Freesurfer Segmentation== | ==Freesurfer Segmentation== | ||
Freesurfer is a useful tool for segmenting a T1-weighted image and building a cortical mesh. To segment the subject's T1-weighted image using freesurfer from the command line type: | Freesurfer is a useful tool for segmenting a T1-weighted image and building a cortical mesh. To segment the subject's T1-weighted image using freesurfer from the command line type: | ||
| − | recon-all -i /home/projects/ | + | recon-all -i /home/projects/MRI/[subjid]/[YYYYMMDD]/[subjid]_WIP_MEMP_VBM_SENSE_13_1_MSE.nii.gz -subjid [subid] -all |
| − | We will follow the steps outlined in the [https://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/LongitudinalTutorial Freesurfer Wiki] for analyzing cortical thickness | + | We will follow the steps outlined in the [https://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/LongitudinalTutorial Freesurfer Wiki] for analyzing cortical thickness. The first steps are described [https://surfer.nmr.mgh.harvard.edu/fswiki/LongitudinalProcessing here] |
Latest revision as of 21:15, 3 December 2015
Contents
Page describing how we analyze cortical thickness with freesurfer
Setting up Freesurfer and MATLAB
First make sure that you have the correct lines in your .bashrc file to run freesurfer:
export FREESURFER_HOME=/usr/local/freesurfer source $FREESURFER_HOME/SetUpFreeSurfer.sh > /dev/null export SUBJECTS_DIR=/mnt/diskArray/projects/freesurfer
Next make sure that you have the right tool boxes in your matlab search path. This should be done through your startup.m file
addpath(genpath('~/git/yeatmanlab'));
Create a T1 weighted nifti image for the subject
Step 1: In a terminal, convert the PAR/REC files to nifti images. You may not need to do this if you have already gone through the Anatomy Pipeline
cd /mnt/diskArray/projects/MRI/[subid] parrec2nii -c -b *.PAR
Step 2: In MATLAB compute the root mean squared (RMS) image. Once again this might have already been done in the Anatomy Pipeline so you can re-use that RMS image
T1path = 'Path to t1 weighted image'; T1path = mri_rms(T1path); % Root mean squared image
Freesurfer Segmentation
Freesurfer is a useful tool for segmenting a T1-weighted image and building a cortical mesh. To segment the subject's T1-weighted image using freesurfer from the command line type:
recon-all -i /home/projects/MRI/[subjid]/[YYYYMMDD]/[subjid]_WIP_MEMP_VBM_SENSE_13_1_MSE.nii.gz -subjid [subid] -all
We will follow the steps outlined in the Freesurfer Wiki for analyzing cortical thickness. The first steps are described here
