User Tools

Site Tools


processing_stream_-_quality_assessment

This is an old revision of the document!


Processing stream - Quality assessment

Motivation

There are several things that can go wrong during the acquisition phase (scanning) of a subject. Among those are several that impact the usability of a subject's dataset. While any given project is still in the stage of data collection (subjects are still being scanned), there is always the chance to decide that a particular subject might introduce too much noise into the eventually performed group statistic and should be discarded (and in this case replaced by another subject).

Relevant for that decision could be one of the following issues:

  • the subject had to exit the scanner before the experiment was completed → usually such a dataset needs to be discarded
  • the subject couldn't restrain from moving their head during the experiment → depending on how difficult it is to find a replacement subject, it is advised to discard such a dataset
  • the scanner produced disproportionally strong noise in the data → if possible, such a dataset should also be discarded

Of course there are still many other possible reasons to discard any given subject (e.g. a score on a questionnaire/behavioral measure indicates that the subject does not fall into the distribution of the examined population of subjects), but especially the second and third issue mentioned above can be detected even before entering a subject's dataset into any given group analysis.

Requirements

To run the fMRI quality checking function, the images need to be in one of the functional imaging data formats currently supported by the xff class (Analyze/NIftI, FMR/STC, VTC).

Usage

For additional options, please consult the fmriquality reference manual page.

The most basic (and pre-configured) way of running fmriquality is by simply passing in the filename(s) or object of the run to check:

% using Analyze files
qas = fmriquality(findfiles(sessionfolder, '*.img'));
 
% alternatively, for FMR
fmr = xff('*.fmr', 'Select FMR for which you want to check the data quality...');
qas = fmriquality(fmr);

The returned variable is of type struct and contains at least the following fields:

   .Dims         1x4 array, size
   .Filename     the first filename given
   .Masks        automatically detected masks (foreground, background, etc)
   .Raw          mean, stdev, and null-voxel image
   .TempFiltered re-created summary values/maps after applying temporal filtering
   .Quality      summary images trying to capture overall quality measures (SNR, CNR, etc.)
   .TC           diverse time courses

The output of fmriquality can be passed on to fmriqasheet, which in turn creates a new figure and displays part of the information in the structure, which can be used to decide on whether or not a subject would likely introduce too much noise/bias at the group level.

processing_stream_-_quality_assessment.1274933303.txt.gz · Last modified: 2010/05/27 04:08 by jochen