Back to list of all tutorials | Back to course page | Previous | Next

TRACULA outputs

Remember...

For each new terminal that you open, you must do:

export SUBJECTS_DIR=$TUTORIAL_DATA/diffusion_tutorial/fs
cd $TUTORIAL_DATA/diffusion_tutorial

Outputs from quality assessment

The diffusion MRI contrast is designed to detect the microscopic motion of water molecules. Unfortunately this makes diffusion MRI particularly sensitive to head motion. Head motion does not result only in misalignment between consecutive DWIs, which could be corrected by registering the DWIs to each other. It can also alter the intensities of the DWIs, and this cannot be corrected by registration. The quality assessment step in TRACULA computes four measures of head motion that summarize the above effects. To see an example of these motion measures for one scan, run:

cat trc/subject1/dmri/dwi_motion.txt

This file contains four values: the average volume-by-volume translation, the average volume-by-volume rotation, the percent of slices with excessive intensity drop-out, and the average drop-out score for slices with excessive intensity drop-out. Excessive intensity drop-out can be identified by slices that are dark, or even completely black, adjacent to normal looking slices.

AvgTranslation AvgRotation PercentBadSlices AvgDropoutScore
0.304789 0.00400526 0 1

Our subject was pretty well-behaved. These motion measures can be used to ensure that groups of subjects are matched with respect to head motion, or to introduce head motion as a nuissance regressor in group analyses. For more information, see Yendiki et al. 2014.

If you want a more detailed break down of motion measures during the scan, run:

cat trc/subject1/dmri/dwi_motion_byvol.txt

This file contains motion measures for each DWI volume of this subject:

TranslationX TranslationY TranslationZ RotationX RotationY RotationZ PercentBadSlices AvgDropoutScore
0 0 0 0 0 0 0 1
0.00424115 -0.00629262 -0.00592024 -0.00034372 -0.000342477 4.26509e-05 0 1
-0.028909 -0.142125 0.00304372 -0.00240405 0.000228201 -0.000238103 0 1
0.0320273 -0.0183681 -0.00179705 0.000278721 -0.000196739 0.00143512 0 1
-0.0587013 0.156494 0.081299 -0.000440066 0.000288482 -0.00193024 0 1
0.0599447 0.155258 -0.024998 0.000163482 -5.79141e-06 0.00301999 0 1
-0.0945436 -0.440727 0.0647819 -0.00143778 -0.000311254 -0.00360836 0 1
0.0549738 -0.0785657 0.0209674 0.000369791 7.12276e-06 0.00315994 0 1
-0.0651216 0.259283 0.0244383 0.000830467 0.000287222 -0.00311584 0 1
...

Outputs from tensor fit

Even though TRACULA relies on the ball-and-stick model instead of the tensor model to reconstruct tracts, the outputs of the tensor fit, especially anisotropy and diffusivity maps, are still useful for group analysis. In addition to the tract-based analysis that can be performed on the tract-based averages of anisotropy and diffusivity, the same maps could also be used for ROI-based or even voxel-based group analyses.

Outputs in native diffusion space

The maps produced by the diffusion tensor fit can be found in the dmri directory which was created when trac-all -prep was run. Let's list these files for subject1:

ls trc/subject1/dmri/dtifit*

The contents of these maps are explained below:

1. dtifit_FA.nii.gz - Fractional anisotropy
2. dtifit_MD.nii.gz - Mean diffusivity
3. dtifit_MO.nii.gz - Mode of the anisotropy
4. dtifit_S0.nii.gz - Non-diffusion weighted (b=0) image
5. dtifit_L1.nii.gz - Primary eigenvalue of the tensor
6. dtifit_L2.nii.gz - Secondary eigenvalue of the tensor
7. dtifit_L3.nii.gz - Tertiary eigenvalue of the tensor
8. dtifit_V1.nii.gz - Primary eigenvector of the tensor
9. dtifit_V2.nii.gz - Secondary eigenvector of the tensor
10. dtifit_V3.nii.gz - Tertiary eigenvector of the tensor

One can use freeview to view the FA map for subject1, overlaid with the subject's anatomical segmentation mapped from T1 to diffusion space. This is a good way to check the quality of the intra-subject registration. You can do this with the command below:

freeview trc/subject1/dmri/dtifit_FA.nii.gz \
         trc/subject1/dlabel/diff/aparc+aseg.bbr.nii.gz:colormap=lut:opacity=.2 &

Note: The above assumed that bbregister was used for the intra-subject registration. If FLIRT was used, the transformed FA map would be called aparc+aseg.flt.nii.gz instead.

Hover your mouse on different voxels throughout the image and notice the FA value change in the bottom toolbar. Look for areas of low FA and areas of high FA. Scroll through the slices using the Page Up and Page Down buttons (or up and down arrows). Feel free to change the contrast by holding down Shift and the left mouse button while dragging the mouse across the viewer window. When you are done, hit the X in the top right to close freeview.

To see the correspondence between the FA and anatomical segmentation in diffusion space and check the accuracy of the intra-subject registration, toggle between the two images by using the 'Alt + C' hot key.

Alternatively, you can adjust the opacity:

You can make one of the views bigger by choosing the 1x1 button ( ). Use the orientation buttons to switch the plane ( , , , ). Put your cursor at the edge of a gyrus or ventricle as a landmark before toggling back and forth between the two volumes. Scroll through the slices and try it again on a different section of the brain. Notice how accurate the registration is - the subject's anatomy should line up well between the FA map and the segmentation. When you are done, be sure to close freeview again.

Tractography outputs

Visualizing the probability distribution of selected white-matter pathways

To visualize the posterior distribution of the inferior longitudinal fasiculus (ILF), do the following:

freeview -v trc/subject1/dmri/dtifit_FA.nii.gz:opacity=.6 \
            trc/subject1/dpath/rh.ilf_avg16_syn_bbr/path.pd.nii.gz:colormap=heat:isosurface=1:isosurface_color=magenta:name=rh.ilf \
            trc/subject1/dpath/lh.ilf_avg16_syn_bbr/path.pd.nii.gz:colormap=heat:isosurface=1:isosurface_color=yellow:name=lh.ilf

Note: The above assumed that bbregister was used for the intra-subject registration and affine registration to MNI space was used for the inter-subject alignment. If anything else was used, bbr and mni would have to be adjusted accordingly.

Click on the 3D button ( )to view isosurfaces of the left and right ILF. Use the 1x1 button ( ) to make it fill the window. You can rotate the image by grabbing it with the mouse and dragging. You can also move the slice planes by grabbing those and dragging them out of the way. Hold down the right mouse button and push the mouse forward for a smooth zoom inwards.

Here we chose to threshold the isosurface display at a very low value (1, which corresponds to the lowest probability voxels). Click on either rh.ilf or lh.ilf, and you can adjust the isosurface threshold to view higher or lower probability parts of the ILF.

Try to make the tracts look similar to the screenshot below.

You can also adjust the heatmap thresholds for the posterior distribution of the ILF. You'll want to go back to the non-3D views ( , , ) to see the effect.

Close freeview when you are done.

Visualizing the end points of selected white-matter pathways on the cortical surface

freeview -f $SUBJECTS_DIR/subject1/surf/lh.inflated:curvature_method=binary:label=trc/subject1/dpath/lh.uf_avg16_syn_bbr/endpt1.surf.label:label_color=green:label=trc/subject1/dpath/lh.uf_avg16_syn_bbr/endpt2.surf.label:label_color=yellow

Visualizing the probability distributions of all white-matter pathways simultaneously

TRACULA also outputs a merged 4D volume, which is a concatenation of the probability distributions of all the pathways that were reconstructed. To visualize all pathways at once, open this merged 4D volume in freeview using the -tv option:

freeview -v  trc/subject1/dmri/dtifit_FA.nii.gz:opacity=.6 \
         -tv trc/subject1/dpath/merged_avg16_syn_bbr.mgz &

Move around in freeview to find the location in this screenshot, or look for an area that interests you.

The lower-left corner of the freeview display gives the color lookup table for all the tracts. The tracts in the merged volume are displayed at 20% of their maximum threshold, by default. Hence the posterior distribution of some of the tracts may not be displayed in its entirety. The threshold for each tract can be adjusted in the 3D and non-3D views by selecting a specific tract from the lookup table and adjusting the threshold button at the bottom of the freeview diplay panel (between 1 and 100).

When you are done, close freeview.

Now that you've seen what the output looks like, you can take a look at the resulting stats for each tract by clicking on the Next button below.

Summary

By the end of this page, you should know how to:

• Check head motion measures

• View the maps output by the tensor fit (e.g., FA) in freeview

• Overlay these with the anatomical segmentation to check the intra-subject registration
• Visualize and manipulate tractography outputs in freeview

Back to list of all tutorials | Back to course page | Previous | Next