The Brain Imaging Data Structure (BIDS) format

• 
  project/
  
•   └── subject
  
•     └── session (optional)
  
•       └── acquisition

The Brain Imaging Data Structure (BIDS) format

• 
• Dicom header information gets stored in .json files
• BIDS specification PDF

Heudiconv

• Heudiconv uses heuristics to find the dicoms corresponding to anatomical, functional, diffusion-weighted etc. sequences
• These heuristics need to be adapted to each study
• E.g. "fmri_1" could be an fMRI session which has "fMRI" in the description and consists of 13650 dicoms
• "fmri_2" could have "fMRI" in the description and consist of 16044 dicoms
• Use what you know about your sequences to build a heuristic!

Example heuristic for heudiconv

import os

def create_key(template, outtype=('nii.gz',), annotation_classes=None):

    if template is None or not template:
        raise ValueError('Template must be a valid format string')
    return template, outtype, annotation_classes

def infotodict(seqinfo):

    print(seqinfo)

    t1w = create_key('sub-{subject}/anat/sub-{subject}_T1w')
    dwi = create_key('sub-{subject}/dwi/sub-{subject}_run-{item:01d}_dwi')
    task1 = create_key('sub-{subject}/func/sub-{subject}_task-task1_rec-{rec}_run-{item:01d}_bold')
    task2 = create_key('sub-{subject}/func/sub-{subject}_task-task2_rec-{rec}_run-{item:01d}_bold')

    info = {t1w: [], dwi: [], task1: [], task2: []}

    for s in seqinfo:
        if (s.dim3 == 176) and ('t1' in s.series_description or 'T1' in s.series_description) and not s.is_derived: # and (s.dim4 == 1)
            info[t1w] = [s.series_id]  # assign if a single series meets criteria
        if (s.dim3 == 3900) and (s.dim4 == 1) and ('dti' in s.series_description or 'DTI' in s.series_description):
            info[dwi].append(s.series_id)  # append if multiple series meet criteria
        if (s.dim3 == 13650) and ('fmri_1' in s.series_description):
            if s.is_motion_corrected:  # exclude non motion corrected series
                info[task1].append({'item': s.series_id, 'rec': 'corrected'})
            else:
                info[task1].append({'item': s.series_id, 'rec': 'uncorrected'})
	      if (s.dim3 == 16044) and ('fmri_2' in s.series_description):
	          if s.is_motion_corrected:  # exclude non motion corrected series
                info[task2].append({'item': s.series_id, 'rec': 'corrected'})
            else:
                info[task2].append({'item': s.series_id, 'rec': 'uncorrected'})
    return info

Containers vs. Virtual Machines

• 
• Source: Infoworld

Docker containers

• Example: fmriprep bundles e.g. Freesurfer, FSL, AFNI, ANTs, among other tools
• You can make your own container when you have finished an analysis to archive everything needed to reproduce, e.g. package versions

Containers

• 
• Source: Gorgolewski et al., 2017, PLoS Computational Biology

Singularity containers

• Docker usually needs certain admin rights
• Singularity has been developed as a high performance cluster implementation for container based processing
• Singularity is installed on Akalla

Singularity workflow

1. Obtain the relevant container in the version that you want
2. Either directly from Singularity Hub (no conversion necessary) or from Docker Hub (docker2singularity conversion necessary)
3. Conversion with docker2singularity requires a local installation of Docker
4. Run the container

Container maintenance

• Containers can be archived in a central location on the server so that everyone has access to all of them
• Containers can be quite large (e.g. >10 GB per container for fmriprep), so deleting old/unused containers is important

Quality control of imaging data with mriqc

singularity run mriqc_container input_folder output_folder \
participant --no-sub -w work_folder --verbose-reports

Quality control of imaging data with mriqc

• 
• Source: Esteban et al., 2017, PLoS ONE

Quality control of imaging data with mriqc

• Overview of image quality metrics (IQMs) such as Framewise Displacement (FD), DVARS, SNR etc. in your raw data
• Both for anatomical as well as functional MRI data

Quality control of imaging data with mriqc

Preprocessing your data with fmriprep

singularity run fmriprep_container input_folder \
output_folder participant -w work_folder

Preprocessing your data with fmriprep - T1

Preprocessing your data with fmriprep - fMRI

Preprocessing your data with fmriprep

• Example output