1. PETsurfer Kinetic Modeling Tutorial Data Description

The dynamic PET data were completely synthesized based on on average ROI time activity curves (TACs) from a study using [11C]SB207145 which targets the serotonin 5HT4 receptor (see Beliveau, et all, 2016). These ROI TACs where then applied to the anatomical MRIs of subjects not in the PET study to create PET-like dynamic images. The anatomical images are 18 real subjects distributed as part of the FSFAST fMRI tutorial. The anatomical images were gradient unwarped prior to analyzing in recon-all. The data are suitable for analyzing using the MRTM1 and MRTM2 models. The data were synthesized with both tissue fraction effects (TFEs) and partial volume effects (PVEs). The PVE is an isotropic 6mm Gaussian kernel. The data were simulated with slight differences in subjects according to two (arbitrary) groups so that there will be a group effect in the bindind potentials.

2. Getting the Data (not necessary for the Boston FreeSurfer Course)

3. Analyzing an Individual Subject

Change directory into the project folder. There you will see 18 folders called s01, s02, etc. Change directory into sb01. You will see several files:

3.1. Examine Input Data

tac.vol.nii.gz -- this is a multi frame volume holding the time activity curves (TACs) for each voxel
frame.time.dat -- text file indicating the time at which each frame was acquired.
subjectname -- text file with the name of the !FreeSurfer-analyzed anatomy of the subject, uw.fsf01anat in this case.  

Examine the contents of frame.time.dat using any text viewer or editor. Verify that the time file has 38 entries. Also notice that the frames are NOT uniformly distributed in time.

Now examine the dimensions of the TAC with

mri_info tac.vol.nii.gz

Verify that the TAC also has 38 frames (the spatial dimensions are 64x64x30)

3.2. Compute the mean of the TAC and Visualize

Run

mri_concat tac.vol.nii.gz --mean --o tac.vol.mean.nii.gz

Visualize the mean and the TAC volume

freeview  tac.vol.nii.gz tac.vol.mean.nii.gz

Click on the "tac.vol.nii.gz" item in the menu, then show the time series. Click around the brain. Generally, you should see a curve that rises and then slowly returns to 0. Note that freeview displays the curve assuming that the time points are uniformly spaced in time (though they are not, see above).

3.3. Register TAC Mean to the Anatomical

mri_coreg --s uw.fsf01anat --mov tac.vol.mean.nii.gz --reg egistration.lta

This will compute a 6DOF registration to the FreeSurfer anatomical. Check the registration with

tkregisterfv --mov tac.vol.mean.nii.gz --reg registration.lta --surfs

Verify that it is accurate.

3.4. Run Partial Volume Correction

mri_gtmpvc --i tac.vol.nii.gz --o gtm.psf6 --psf 6 - --seg gtmseg.mgz --reg registration.lta \
   --km-ref 8 47 --km-hb 11 12 50 51 --auto-mask 1 .1 \
   --default-seg-merge --no-rescale --update-tt  --mgx .01

This command does a lot so let's break it down: --i tac.vol.nii.gz : this is the multi-frame input that we want to analyze
--o gtm.psf : this is the output folder