< Neuroimaging Data Processing

Overview of FSL

FSL is a comprehensive library of analysis tools for FMRI, MRI and DTI brain imaging data. It could be used either by FSL main GUI or the commond-lines. The advantage of the GUI lies in its simplicity but the drawback is the less flexibility comparing to the command-lines. FSLView as the display tool could be completely separate use from processing and analysis. Currently FSL only accepts the input files in NIFTI format, and the DICOM files need to be converted to NIFTI after acquisition and before processing. The convertion from DICOM to NIFIT could be fulfilled by several methods, such as dcm2nii from mricron or dcmstack. In the GUI there are nine popular modules are defined, and depending on the input file types, they fall into functional-oriented, structural-oriented or diffusion-oriented modules. The structural-oriented functions include "BET brain extraction", "FAST segmentation", "FLIRT: linear registration"; the functional-oriented functions invole "SUSAN noise reduction","FEAT FMRI analysis" and "MELODIC ICA"; the DTI-oriented functions incorporate "FDT diffusion". Besides these functional modules, there is a simulator called "POSSIUM MRI simulator", which is a software tool to produce realistic simulated MRI and FMRI images or time series. POSSUM (Physics-Oriented Simulated Scanner for Understanding MRI) includes tools for the pulse sequence generation, signal generation, noise addition and image reconstruction.

The GUI of FSL with nine modules

How to launch FSL

Suppose all the environmental vairables are settled well for FSL, there are two ways to lauch FSL, either from the user GUI or command-line.

User GUI

Just type fsl in the console and get access to the main GUI, from where you could reach to the modules by clicking on the options.

 fsl

Depending on your operating system, you might need to do some additional steps prior to type in this starting command into your console/terminal

operating systempreparing steps
Linux/Unixto tell your terminal, that the following commands are referred to the fsl-toolbox, your need to type in
 FSL

in capital letters. Afterwards you will get a message, that your terminal will accept fsl-commands.

WindowsFSL is not running under Windows. Therefore you will need additional software, to emulate another operating system. FSL will run perfectly under a vitual machine like VMware Player[1] or VirtualBox[2]
Mac OS Xno prior steps needed

If you want to jump to specific modules directly, then just type in the name starting with a capital letter, e.g. Melodic for Melodic ICA

 Melodic
Command-line

If you want to launch a function from the command line, just type the function name in lowercase, along with the necessary arguments to the function.

 melodic -i inputfile

FSL function overviews

Structural image processing
FunctionGUI/commandAimOnline document
BETBothBrain extractionhttp://fsl.fmrib.ox.ac.uk/fsl/fslwiki/BET
FASTBothTissue segmentationhttp://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FAST
FIRSTCommandSubcortical segmentationhttp://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FIRST
FLIRTBothLinear registrationhttp://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FLIRT
FNIRTCommandNonlinear registrationhttp://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FNIRT
FUGUECommandEPI distortion correctionhttp://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FUGUE
SIENACommandAtrophy analysishttp://fsl.fmrib.ox.ac.uk/fsl/fslwiki/SIENA
FSL-VBMCommandGrey matter densityhttp://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FSLVBM
Functional MRI data analysis

In FSL, the functional images could be analyzed in two integration workflows, rather than step-by-step as in SPM and AFNI. Analysis for an experiment can be set up in less than 1 minute and finished in 5-20 minutes for the first-level session, producing a web page analysis report, including colour activation images and time-course plots of data vs model.

FEAT is based on general linear modelling (GLM). It allows to describe the experimental design; then a model is created that should fit the data, telling you where the brain has activated in response to the stimuli. In FEAT, the GLM method used on first-level (time-series) data is known as FILM. FILM uses a robust and accurate nonparametric estimation of time series autocorrelation to prewhiten each voxel's time series; this gives improved estimation efficiency compared with methods that do not pre-whiten.

GUI for FEAT with six tabs

MELODIC is model-free analysis which uses Independent Component Analysis (ICA) to decompose a single or multiple 4D data sets into different spatial and temporal components. MELODIC can pick out different activation and artefactual components without any explicit time series model being specified.

GUI for MELODIC with six tabs
FSLUTILS

FSLUTILS is a set of useful command-line utilities which allow the conversion, processing etc. of Analyze and Nifti format data sets. Many of them work on both 3D and 4D data. For each of these programs, type just the program name to get the usage help. Those commands could be broadly divided into four categories depending on their aims.

1. Mathematical manipulation of images

FunctionAimExamplary command
fslccRun cross-correlations between every volume in one 4D data set with every volume in anotherfslcc input1.nii input2.nii (measure similarities in ICA outputs)
fslfftOutputs the Fast-Fourier Transform (or inverse) for a complex input volumefslfft input.vol output.vol
fslmathsSimple but powerful program to allow mathematical manipulation of imagesfslmaths func_data.nii -bptf 25.0 -1 func_data_tempfilt.nii (bandpass temporal filtering)
fslmeantsOutput the average timeseries of a set of voxelsfslmeants -i input.nii -o mean.txt -m mask (preliminary step for seed-based correlelation)
fslroiExtract region of interest (ROI) from an imagefslroi input.nii output_roi.nii 0 1 (extract the first volume)

2. File combination and splitting

FunctionAimExamplary command
fslcomplexAllow 3D or 4D complex image files to be split or constructed from corresponding real componentsfslcomplex -complexsplit source dest 1 3
fslinterleaveInterleave two inputs to form a combined image, and this combination is in Z axis onlyfslinterleave input1 input2 output
fslmergeConcatenate image files into a single output. This concatenation can be in time, X, Y or Z axisfslmerge -t outfile infile1 infile2 (concatenate files in time series)
fslsliceSplit a 3D file into lots of 2D files along Z-axis.fslslice volfile slicefile
fslsplitSplit a 4D file into lots of 3D files in either time, X,Y,or Z axisfslsplit -t input outfileprefix (split volumes along time series)

3. Header-related utilities

FunctionAimExamplary command
fslcpgeomCopy certain parts of the header information from one image to another with identical file type.fslcpgeom source dest
fslcreatehdCreates a new image header along with a zero intensity data image.fslcreatehd
fslhdReport every field of an Analyze or Nifti headerfslhd input.nii
fslinfoReport a basic subset of an Analyze or Nifti header.fslinfo input.nii
fslvalReport a particular parameter from an image headerfslval input.nii

4. Orientation-related Utilities

FunctionAimExamplary command
fslreorient2stdReorient an image to match the orientation of the standard template images (MNI152), but this only rotated to 90, 180 or 270 degree and can't substitute coregistrationfslreorient2std input output
fslorientAdvanced tool that reports or sets the orientation information in a filefslorient -forceradiological input
fslswapdimAdvanced tool that re-orders the data storage to permit changes between axial, sagittal and coronal slicingfslswapdim
FSLView

The FSLView could be launched from terminal by typingː

 fslview --help

1.Cursor view

FSL cursorview

Box-1ː show and control the cursor position in voxels

Box-2ː show and control the cursor position in mm

Box-3ː shows and controls the current volume in a 4D dataset

Box-4ː shows the value in the voxel at the current cursor position

2.FSL image overlay

FSL image overlay

Box-1ː The layer marked with chevrons to the right of it is the "main layer". This layer is the first one loaded during any given session and cannot be removed from the list as much of the viewer's display capabilities are determined from this layer's attributes.

Box-2ː Visibilty" checkbox controls whether or not a given layer is visible; "Lock" checkbox determines if a given layer can be edited or not; "Transparency" slider determines how the selected layer blends with the layers below it

Resources and references

  1. http://www.vmware.com/products/player/
  2. https://www.virtualbox.org/
This article is issued from Wikibooks. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.