After all processing has been completed, the designated output directory for the qap_functional_temporal.py workflow will contain a set of pdf files that contain the relevant reports for the set of scans undergoing quality assessment. The set of output pdfs includes one pdf file per input scan, e.g.: qap_functional_temporal_sub-01.pdf, which contains the tSNR maps, FD plot and QC metrics for that scan. There will also be a group report pdf in that main output directory, e.g.: qap_functional_temporal_group.pdf, that contains summary metrics for the entire set of scans.

For the individual scan reports:

The tSNR Slice Mosaic

The image showing the tSNR map for a particular scan, e.g.:

tSNR volume, subject sub-01 (session_1_func_1)

is the rendering of the temporal signal-to-noise ratio information from the scan. The tSNR plot is similar to the mean EPI plot generated by the qap_functional_spatial.py workflow, in that both metrics reduce the 4-dimensional scan to a representative 3-dimensional volume that is then split and displayed as a stack of axial (horizontal) slices. For this tSNR plot, the mean of each voxel's timeseries is also computed and is then divided by the timeseries standard deviation. Hence, the tSNR plot shows the voxels in which one would expect to have SNR good enough for statistical analyses. Differences in tSNR are particularly important for comparing the results from region of interest (ROI) analyses, since any observed functional differences might actually be attributable to systematic differences in SNR across the regions being compared. You can learn more about the utility of tSNR plots for fMRI analyses here.

The framewise displacement plot

The image showing the framewise displacement plot for the scan being assessed, e.g.:

FD, subject sub-01 (session_1_func_1)

is the framewise displacement that occurred throughout the scan. This is a temporal motion quality assurance metric and tracks head motions over time, making it easy to determine whether or not the data potentially suffered from significant corruption due to motion. For instance, it is possible to detect if the participant's head was slowly sinking into the cushions in the head coil, or whether the participant was possibly restless or agitated, which would result in several postion changes or movement spikes. The framewise displacement is a frame-by-frame representation of the differences between the BOLD signal intensity value of the n and n+1 timepoints, the n+1 and n+2 timepoints, and so on. The report page for framewise displacement in the functional scan includes both this frame-by-frame plot, as well as a histogram that can be used to visually determine what proportion of timepoints exceeded some pre-set movement threshold (eg: 0.2 mm).

The Temporal Metrics computed on the Functional Scan

The metrics displayed in the Summary Report were computed using the qap_functional_temporal.py workflow and have been displayed as violin plots. Eg:

QC measures (subject sub-01_session_1)

The stars in these plots denote where the score for this particular scan falls in the distribution of all scores for scans that were included as inputs to the the functional-temporal workflow.

The metrics computed are as follows:

  1. DVARS - Spatial standard deviation of the voxelwise temporal derivates
  2. GCOR - Global Correlation
  3. Mean RMSD - Mean Fractional Displacement
  4. Outliers - Mean fraction of outliers per fMRI volume
  5. Quality - Median Distance Index

All metrics are described in more detail in the Taxonomy of QA Measures section of the QAP documentation. Please refer to the QAP website for descriptions of these metrics.

For the group reports:

The violin plots included in the group report, e.g.: QC measures (session_1) are a graphical representation of the columnar values in the qap_functional_temporal.csv file that was created in the main output directory for the workflow. The scores for each metric described above were aggregated to create the distributions that were plotted in both the individual and group reports. Hence, the violin plots in the individual scan reports and the group reports are identical, except that the group reports do not contain any stars denoting individual scans. These group reports are intended to provide the user a means of visually inspecting the overall quality of the temporal data for that group of functional scans.