Freesurfer Abridged: Difference between revisions
Jump to navigation
Jump to search
(One intermediate revision by one other user not shown) | |||
Line 37: | Line 37: | ||
#* Semcat processing | #* Semcat processing | ||
#** In terminal, <code> matlab &</code> | #** In terminal, <code> matlab &</code> | ||
#** In matlab, <code> | #** In matlab, <code>startup_spm12</code> | ||
#** | #** **Check with your local Principle Investigator™ or Graduate Work Horse™ for more details about Semcat bold files** | ||
#* LDT processing | #* LDT processing | ||
#** In terminal, <code>matlab &</code> | #** In terminal, <code>matlab &</code> | ||
#** nii_4D_drop_vols('f.nii', [# vols to drop | #** nii_4D_drop_vols('f.nii', [# vols to drop]) | ||
#** set_TR(2.047, [1 2 3 4 5 6], 'f.nii') | #** set_TR(2.047, [1 2 3 4 5 6], 'f.nii') | ||
#** d=datenum('DD-Mmm-YYYY') (e.g. '05-Oct-1997') | #** d=datenum('DD-Mmm-YYYY') (e.g. '05-Oct-1997') | ||
Line 47: | Line 47: | ||
#** FSTSExtractor('prefix', 'LDT_', 'subject', subj_ID, 'run', i, 'tr', 2.047, 'volumes_dropped', <initial # drop vols>, 'date', d) | #** FSTSExtractor('prefix', 'LDT_', 'subject', subj_ID, 'run', i, 'tr', 2.047, 'volumes_dropped', <initial # drop vols>, 'date', d) | ||
#*** This step must be completed for each f.nii file in each bold folder | #*** This step must be completed for each f.nii file in each bold folder | ||
#*** IT IS VERY IMPORTANT to look at what volumes to drop. It will likely not be the same for each bold file, unless it is a newer LDT participant. You can find more information | #*** IT IS VERY IMPORTANT to look at what volumes to drop. It will likely not be the same for each bold file, unless it is a newer LDT participant. You can find more information [[http://ccnlab.psy.buffalo.edu/wiki/index.php/Freesurfer_BOLD_files]] here. | ||
# Analysis | # Analysis | ||
#* cd to subject folder | #* cd to subject folder |
Latest revision as of 11:55, 11 May 2022
Setting Up
- Note: I am using FS_0184 (semcat) as an example; fill in the parts specific to FS_0184 with the parts specific to your subject
This is an abridged version of the more detailed list provided elsewhere in the wiki.
Structural
- Slicer
- In terminal,
Slicer &
- Load data brainmask.nii
- Save data as brainmask.mgz
- Results: creates brainmask.mgz file
- In terminal,
- Autorecon 1
- In terminal,
autorecon1_SC.sh FS_0184
- Don't bother editing after this step; skip to autorecon 2
- Result: populates mri folder (includes "T1.mgz")
- In terminal,
- Autorecon 2
- In terminal,
autorecon2_SC.sh FS_0184
- Editing
- One terminal,
tksurfer FS_0184 lh inflated
- In another terminal,
tkmedit FS_0184 brainmask.mgz -aux T1.mgz -surfs
- Look for any noticable spikes in the tksurfer brain and edit the brainmask (tkmedit). See wiki page for details; this is the most labor intensive part of anatomical so make sure you're doing it right.
- One terminal,
- Results: creates lh.cortex.label and rh.cortex.label in label folder
- In terminal,
- Autorecon 3
- In terminal,
autorecon3_SC.sh FS_0184
- Editing
- In terminal,
tksurfer -annot aparc FS_0184 lh pial
- If needed, reference
tksurfer -annot aparc fsaverage lh pial
- Repeat with right hemi
- Make sure to EXPORT ANNOTATION! Do not "save overlay!"
- Again, be sure to see the autorecon3 page
- In terminal,
- Result: populates label folder (includes "lh.aparc.annot" and "rh.aparc.annot")
- In terminal,
Functional
- Bold file preparation
- Semcat processing
- In terminal,
matlab &
- In matlab,
startup_spm12
- **Check with your local Principle Investigator™ or Graduate Work Horse™ for more details about Semcat bold files**
- In terminal,
- LDT processing
- In terminal,
matlab &
- nii_4D_drop_vols('f.nii', [# vols to drop])
- set_TR(2.047, [1 2 3 4 5 6], 'f.nii')
- d=datenum('DD-Mmm-YYYY') (e.g. '05-Oct-1997')
- subj_ID='0183'
- FSTSExtractor('prefix', 'LDT_', 'subject', subj_ID, 'run', i, 'tr', 2.047, 'volumes_dropped', <initial # drop vols>, 'date', d)
- This step must be completed for each f.nii file in each bold folder
- IT IS VERY IMPORTANT to look at what volumes to drop. It will likely not be the same for each bold file, unless it is a newer LDT participant. You can find more information [[1]] here.
- In terminal,
- Semcat processing
- Analysis
- cd to subject folder
- mkanalysis.sh
- Contrasts
- cd to subject folder
- Run
mkcontrast.sh
- Pre-processing
- Run
preproc-sess -per-run -fwhm 6 -surface self lhrh -fsd bold -s FS_<subj_ID>
- Run
- Selxavg
- Run
selxavg3-sess -s <subjectname> -analysis <analysis path>
- Analysis path is where your LDT.?h.sm6 folders are.
- Run
- (Optional) Clut and Subparcellation
- Run
clut.sh FS_<subj_ID>
- Run
subpar.py FS_<subj_ID>
- Run
- Blobs