Interfaces

In Nipype, interfaces are python modules that allow you to use various external packages (e.g. FSL, SPM or FreeSurfer), even if they themselves are written in another programming language than python. Such an interface knows what sort of options an external program has and how to execute it.

Interfaces vs. Workflows

Interfaces are the building blocks that solve well-defined tasks. We solve more complex tasks by combining interfaces with workflows:

Interfaces Workflows
Wrap *unitary* tasks Wrap *meta*-tasks
  • implemented with nipype interfaces wrapped inside ``Node`` objects
  • subworkflows can also be added to a workflow without any wrapping
  • Keep track of the inputs and outputs, and check their expected types Do not have inputs/outputs, but expose them from the interfaces wrapped inside
    Do not cache results (unless you use [interface caching](advanced_interfaces_caching.ipynb)) Cache results
    Run by a nipype plugin Run by a nipype plugin

    To illustrate why interfaces are so useful, let’s have a look at the brain extraction algorithm BET from FSL. Once in its original framework and once in the Nipype framework.

    BET in the original framework

    Let’s take a look at one of the T1 images we have in our dataset on which we want to run BET.

    from nilearn.plotting import plot_anat
    %matplotlib inline
    import matplotlib.pyplot as plt
    plot_anat('/data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w.nii.gz', title='original',
              display_mode='ortho', dim=-1, draw_cross=False, annotate=False);
    
    ../../_images/basic_interfaces_4_0.png

    In its simplest form, you can run BET by just specifying the input image and tell it what to name the output image:

    bet <input> <output>
    
    %%bash
    
    FILENAME=/data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w
    
    bet ${FILENAME}.nii.gz /output/sub-01_ses-test_T1w_bet.nii.gz
    

    Let’s take a look at the results:

    plot_anat('/output/sub-01_ses-test_T1w_bet.nii.gz', title='original',
              display_mode='ortho', dim=-1, draw_cross=False, annotate=False);
    
    ../../_images/basic_interfaces_8_0.png

    Perfect! Exactly what we want. Hmm… what else could we want from BET? Well, it’s actually a fairly complicated program. As is the case for all FSL binaries, just call it with the help flag -h to see all its options.

    !bet -h
    
    Usage:    bet <input> <output> [options]
    
    Main bet2 options:
      -o          generate brain surface outline overlaid onto original image
      -m          generate binary brain mask
      -s          generate approximate skull image
      -n          don't generate segmented brain image output
      -f <f>      fractional intensity threshold (0->1); default=0.5; smaller values give larger brain outline estimates
      -g <g>      vertical gradient in fractional intensity threshold (-1->1); default=0; positive values give larger brain outline at bottom, smaller at top
      -r <r>      head radius (mm not voxels); initial surface sphere is set to half of this
      -c <x y z>  centre-of-gravity (voxels not mm) of initial mesh surface.
      -t          apply thresholding to segmented brain image and mask
      -e          generates brain surface as mesh in .vtk format
    
    Variations on default bet2 functionality (mutually exclusive options):
      (default)   just run bet2
      -R          robust brain centre estimation (iterates BET several times)
      -S          eye & optic nerve cleanup (can be useful in SIENA)
      -B          bias field & neck cleanup (can be useful in SIENA)
      -Z          improve BET if FOV is very small in Z (by temporarily padding end slices)
      -F          apply to 4D FMRI data (uses -f 0.3 and dilates brain mask slightly)
      -A          run bet2 and then betsurf to get additional skull and scalp surfaces (includes registrations)
      -A2 <T2>    as with -A, when also feeding in non-brain-extracted T2 (includes registrations)
    
    Miscellaneous options:
      -v          verbose (switch on diagnostic messages)
      -h          display this help, then exits
      -d          debug (don't delete temporary intermediate images)
    

    We see that BET can also return a binary brain mask as a result of the skull-strip, which can be useful for masking our GLM analyses (among other things). Let’s run it again including that option and see the result.

    %%bash
    
    FILENAME=/data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w
    
    bet ${FILENAME}.nii.gz /output/sub-01_ses-test_T1w_bet.nii.gz -m
    
    plot_anat('/output/sub-01_ses-test_T1w_bet_mask.nii.gz', title='original',
              display_mode='ortho', dim=-1, draw_cross=False, annotate=False);
    
    /opt/miniconda-latest/envs/neuro/lib/python3.7/site-packages/nilearn/image/resampling.py:512: UserWarning: Casting data from int32 to float32
      warnings.warn("Casting data from %s to %s" % (data.dtype.name, aux))
    
    ../../_images/basic_interfaces_13_1.png

    Now let’s look at the BET interface in Nipype. First, we have to import it.

    BET in the Nipype framework

    So how can we run BET in the Nipype framework?

    First things first, we need to import the BET class from Nipype’s interfaces module:

    from nipype.interfaces.fsl import BET
    

    Now that we have the BET function accessible, we just have to specify the input and output file. And finally, we have to run the command. So exactly like in the original framework.

    skullstrip = BET()
    skullstrip.inputs.in_file = "/data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w.nii.gz"
    skullstrip.inputs.out_file = "/output/T1w_nipype_bet.nii.gz"
    res = skullstrip.run()
    

    If we now look at the results from Nipype, we see that it is exactly the same as before.

    plot_anat('/output/T1w_nipype_bet.nii.gz', title='original',
              display_mode='ortho', dim=-1, draw_cross=False, annotate=False);
    
    ../../_images/basic_interfaces_20_0.png

    This is not surprising, because Nipype used exactly the same bash code that we were using in the original framework example above. To verify this, we can call the cmdline function of the constructed BET instance.

    print(skullstrip.cmdline)
    
    bet /data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w.nii.gz /output/T1w_nipype_bet.nii.gz
    

    Another way to set the inputs on an interface object is to use them as keyword arguments when you construct the interface instance. Let’s write the Nipype code from above in this way, but let’s also add the option to create a brain mask.

    skullstrip = BET(in_file="/data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w.nii.gz",
                     out_file="/output/T1w_nipype_bet.nii.gz",
                     mask=True)
    res = skullstrip.run()
    

    Now if we plot this, we see again that this worked exactly as before. No surprise there.

    plot_anat('/output/T1w_nipype_bet_mask.nii.gz', title='after skullstrip',
              display_mode='ortho', dim=-1, draw_cross=False, annotate=False);
    
    /opt/miniconda-latest/envs/neuro/lib/python3.7/site-packages/nilearn/image/resampling.py:512: UserWarning: Casting data from int32 to float32
      warnings.warn("Casting data from %s to %s" % (data.dtype.name, aux))
    
    ../../_images/basic_interfaces_26_1.png

    Help Function

    But how did we know what the names of the input parameters are? In the original framework, we were able to just run BET, without any additional parameters to get an information page. In the Nipype framework, we can achieve the same thing by using the help() function on an interface class. For the BET example, this is:

    BET.help()
    
    Wraps the executable command ``bet``.
    
    FSL BET wrapper for skull stripping
    
    For complete details, see the `BET Documentation.
    <https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/BET/UserGuide>`_
    
    Examples
    --------
    >>> from nipype.interfaces import fsl
    >>> btr = fsl.BET()
    >>> btr.inputs.in_file = 'structural.nii'
    >>> btr.inputs.frac = 0.7
    >>> btr.inputs.out_file = 'brain_anat.nii'
    >>> btr.cmdline
    'bet structural.nii brain_anat.nii -f 0.70'
    >>> res = btr.run() # doctest: +SKIP
    
    Inputs::
    
            [Mandatory]
            in_file: (a pathlike object or string representing an existing file)
                    input file to skull strip
                    argument: ``%s``, position: 0
    
            [Optional]
            out_file: (a pathlike object or string representing a file)
                    name of output skull stripped image
                    argument: ``%s``, position: 1
            outline: (a boolean)
                    create surface outline image
                    argument: ``-o``
            mask: (a boolean)
                    create binary mask image
                    argument: ``-m``
            skull: (a boolean)
                    create skull image
                    argument: ``-s``
            no_output: (a boolean)
                    Don't generate segmented output
                    argument: ``-n``
            frac: (a float)
                    fractional intensity threshold
                    argument: ``-f %.2f``
            vertical_gradient: (a float)
                    vertical gradient in fractional intensity threshold (-1, 1)
                    argument: ``-g %.2f``
            radius: (an integer)
                    head radius
                    argument: ``-r %d``
            center: (a list of at most 3 items which are an integer)
                    center of gravity in voxels
                    argument: ``-c %s``
            threshold: (a boolean)
                    apply thresholding to segmented brain image and mask
                    argument: ``-t``
            mesh: (a boolean)
                    generate a vtk mesh brain surface
                    argument: ``-e``
            robust: (a boolean)
                    robust brain centre estimation (iterates BET several times)
                    argument: ``-R``
                    mutually_exclusive: functional, reduce_bias, robust, padding,
                      remove_eyes, surfaces, t2_guided
            padding: (a boolean)
                    improve BET if FOV is very small in Z (by temporarily padding end
                    slices)
                    argument: ``-Z``
                    mutually_exclusive: functional, reduce_bias, robust, padding,
                      remove_eyes, surfaces, t2_guided
            remove_eyes: (a boolean)
                    eye & optic nerve cleanup (can be useful in SIENA)
                    argument: ``-S``
                    mutually_exclusive: functional, reduce_bias, robust, padding,
                      remove_eyes, surfaces, t2_guided
            surfaces: (a boolean)
                    run bet2 and then betsurf to get additional skull and scalp surfaces
                    (includes registrations)
                    argument: ``-A``
                    mutually_exclusive: functional, reduce_bias, robust, padding,
                      remove_eyes, surfaces, t2_guided
            t2_guided: (a pathlike object or string representing a file)
                    as with creating surfaces, when also feeding in non-brain-extracted
                    T2 (includes registrations)
                    argument: ``-A2 %s``
                    mutually_exclusive: functional, reduce_bias, robust, padding,
                      remove_eyes, surfaces, t2_guided
            functional: (a boolean)
                    apply to 4D fMRI data
                    argument: ``-F``
                    mutually_exclusive: functional, reduce_bias, robust, padding,
                      remove_eyes, surfaces, t2_guided
            reduce_bias: (a boolean)
                    bias field and neck cleanup
                    argument: ``-B``
                    mutually_exclusive: functional, reduce_bias, robust, padding,
                      remove_eyes, surfaces, t2_guided
            output_type: ('NIFTI' or 'NIFTI_PAIR' or 'NIFTI_GZ' or
                      'NIFTI_PAIR_GZ')
                    FSL output type
            args: (a string)
                    Additional parameters to the command
                    argument: ``%s``
            environ: (a dictionary with keys which are a bytes or None or a value
                      of class 'str' and with values which are a bytes or None or a
                      value of class 'str', nipype default value: {})
                    Environment variables
    
    Outputs::
    
            out_file: (a pathlike object or string representing a file)
                    path/name of skullstripped file (if generated)
            mask_file: (a pathlike object or string representing a file)
                    path/name of binary brain mask (if generated)
            outline_file: (a pathlike object or string representing a file)
                    path/name of outline file (if generated)
            meshfile: (a pathlike object or string representing a file)
                    path/name of vtk mesh file (if generated)
            inskull_mask_file: (a pathlike object or string representing a file)
                    path/name of inskull mask (if generated)
            inskull_mesh_file: (a pathlike object or string representing a file)
                    path/name of inskull mesh outline (if generated)
            outskull_mask_file: (a pathlike object or string representing a file)
                    path/name of outskull mask (if generated)
            outskull_mesh_file: (a pathlike object or string representing a file)
                    path/name of outskull mesh outline (if generated)
            outskin_mask_file: (a pathlike object or string representing a file)
                    path/name of outskin mask (if generated)
            outskin_mesh_file: (a pathlike object or string representing a file)
                    path/name of outskin mesh outline (if generated)
            skull_mask_file: (a pathlike object or string representing a file)
                    path/name of skull mask (if generated)
            skull_file: (a pathlike object or string representing a file)
                    path/name of skull file (if generated)
    
    References:
    -----------
    BibTeX('@article{JenkinsonBeckmannBehrensWoolrichSmith2012,author={M. Jenkinson, C.F. Beckmann, T.E. Behrens, M.W. Woolrich, and S.M. Smith},title={FSL},journal={NeuroImage},volume={62},pages={782-790},year={2012},}', key='JenkinsonBeckmannBehrensWoolrichSmith2012')
    

    As you can see, we get three different pieces of information. First, a general explanation of the class.

    Wraps command **bet**
    
    Use FSL BET command for skull stripping.
    
    For complete details, see the `BET Documentation.
    <http://www.fmrib.ox.ac.uk/fsl/bet2/index.html>`_
    
    Examples
    --------
    >>> from nipype.interfaces import fsl
    >>> from nipype.testing import  example_data
    >>> btr = fsl.BET()
    >>> btr.inputs.in_file = example_data('structural.nii')
    >>> btr.inputs.frac = 0.7
    >>> res = btr.run() # doctest: +SKIP
    

    Second, a list of all possible input parameters.

    Inputs:
    
        [Mandatory]
        in_file: (an existing file name)
            input file to skull strip
            flag: %s, position: 0
    
        [Optional]
        args: (a string)
            Additional parameters to the command
            flag: %s
        center: (a list of at most 3 items which are an integer (int or
             long))
            center of gravity in voxels
            flag: -c %s
        environ: (a dictionary with keys which are a value of type 'str' and
             with values which are a value of type 'str', nipype default value:
             {})
            Environment variables
        frac: (a float)
            fractional intensity threshold
            flag: -f %.2f
        functional: (a boolean)
            apply to 4D fMRI data
            flag: -F
            mutually_exclusive: functional, reduce_bias, robust, padding,
             remove_eyes, surfaces, t2_guided
        ignore_exception: (a boolean, nipype default value: False)
            Print an error message instead of throwing an exception in case the
            interface fails to run
        mask: (a boolean)
            create binary mask image
            flag: -m
        mesh: (a boolean)
            generate a vtk mesh brain surface
            flag: -e
        no_output: (a boolean)
            Don't generate segmented output
            flag: -n
        out_file: (a file name)
            name of output skull stripped image
            flag: %s, position: 1
        outline: (a boolean)
            create surface outline image
            flag: -o
        output_type: ('NIFTI_PAIR' or 'NIFTI_PAIR_GZ' or 'NIFTI_GZ' or
             'NIFTI')
            FSL output type
        padding: (a boolean)
            improve BET if FOV is very small in Z (by temporarily padding end
            slices)
            flag: -Z
            mutually_exclusive: functional, reduce_bias, robust, padding,
             remove_eyes, surfaces, t2_guided
        radius: (an integer (int or long))
            head radius
            flag: -r %d
        reduce_bias: (a boolean)
            bias field and neck cleanup
            flag: -B
            mutually_exclusive: functional, reduce_bias, robust, padding,
             remove_eyes, surfaces, t2_guided
        remove_eyes: (a boolean)
            eye & optic nerve cleanup (can be useful in SIENA)
            flag: -S
            mutually_exclusive: functional, reduce_bias, robust, padding,
             remove_eyes, surfaces, t2_guided
        robust: (a boolean)
            robust brain centre estimation (iterates BET several times)
            flag: -R
            mutually_exclusive: functional, reduce_bias, robust, padding,
             remove_eyes, surfaces, t2_guided
        skull: (a boolean)
            create skull image
            flag: -s
        surfaces: (a boolean)
            run bet2 and then betsurf to get additional skull and scalp surfaces
            (includes registrations)
            flag: -A
            mutually_exclusive: functional, reduce_bias, robust, padding,
             remove_eyes, surfaces, t2_guided
        t2_guided: (a file name)
            as with creating surfaces, when also feeding in non-brain-extracted
            T2 (includes registrations)
            flag: -A2 %s
            mutually_exclusive: functional, reduce_bias, robust, padding,
             remove_eyes, surfaces, t2_guided
        terminal_output: ('stream' or 'allatonce' or 'file' or 'none')
            Control terminal output: `stream` - displays to terminal immediately
            (default), `allatonce` - waits till command is finished to display
            output, `file` - writes output to file, `none` - output is ignored
        threshold: (a boolean)
            apply thresholding to segmented brain image and mask
            flag: -t
        vertical_gradient: (a float)
            vertical gradient in fractional intensity threshold (-1, 1)
            flag: -g %.2f
    

    And third, a list of all possible output parameters.

    Outputs:
    
        inskull_mask_file: (a file name)
            path/name of inskull mask (if generated)
        inskull_mesh_file: (a file name)
            path/name of inskull mesh outline (if generated)
        mask_file: (a file name)
            path/name of binary brain mask (if generated)
        meshfile: (a file name)
            path/name of vtk mesh file (if generated)
        out_file: (a file name)
            path/name of skullstripped file (if generated)
        outline_file: (a file name)
            path/name of outline file (if generated)
        outskin_mask_file: (a file name)
            path/name of outskin mask (if generated)
        outskin_mesh_file: (a file name)
            path/name of outskin mesh outline (if generated)
        outskull_mask_file: (a file name)
            path/name of outskull mask (if generated)
        outskull_mesh_file: (a file name)
            path/name of outskull mesh outline (if generated)
        skull_mask_file: (a file name)
            path/name of skull mask (if generated)
    

    So here we see that Nipype also has output parameters. This is very practical. Because instead of typing the full path name to the mask volume, we can also more directly use the mask_file parameter.

    print(res.outputs.mask_file)
    
    /output/T1w_nipype_bet_mask.nii.gz
    

    Interface errors

    To execute any interface class we use the run method on that object. For FSL, Freesurfer, and other programs, this will just make a system call with the command line we saw above. For MATLAB-based programs like SPM, it will actually generate a .m file and run a MATLAB process to execute it. All of that is handled in the background.

    But what happens if we didn’t specify all necessary inputs? For instance, you need to give BET a file to work on. If you try and run it without setting the input in_file, you’ll get a Python exception before anything actually gets executed:

    skullstrip2 = BET()
    try:
        skullstrip2.run()
    except(ValueError) as err:
        print("ValueError:", err)
    else:
        raise
    
    ValueError: BET requires a value for input 'in_file'. For a list of required inputs, see BET.help()
    

    Nipype also knows some things about what sort of values should get passed to the inputs and will raise (hopefully) informative exceptions when they are violated – before anything gets processed. For example, BET just lets you say “create a mask,” it doesn’t let you name it. You may forget this, and try to give it a name. In this case, Nipype will raise a TraitError telling you what you did wrong:

    try:
        skullstrip.inputs.mask = "mask_file.nii"
    except(Exception) as err:
        if "TraitError" in str(err.__class__):
            print("TraitError:", err)
        else:
            raise
    else:
        raise
    
    TraitError: The 'mask' trait of a BETInputSpec instance must be a boolean, but a value of 'mask_file.nii' <class 'str'> was specified.
    

    Additionally, Nipype knows that, for inputs corresponding to files you are going to process, they should exist in your file system. If you pass a string that doesn’t correspond to an existing file, it will error and let you know:

    try:
        skullstrip.inputs.in_file = "/data/oops_a_typo.nii"
    except(Exception) as err:
        if "TraitError" in str(err.__class__):
            print("TraitError:", err)
        else:
            raise
    else:
        raise
    
    TraitError: The 'in_file' trait of a BETInputSpec instance must be a pathlike object or string representing an existing file, but a value of '/data/oops_a_typo.nii' <class 'str'> was specified.
    

    It turns out that for default output files, you don’t even need to specify a name. Nipype will know what files are going to be created and will generate a name for you:

    skullstrip = BET(in_file="/data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w.nii.gz")
    print(skullstrip.cmdline)
    
    bet /data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w.nii.gz /home/neuro/workshop_weizmann/workshop/nipype/notebooks/sub-01_ses-test_T1w_brain.nii.gz
    

    Note that it is going to write the output file to the local directory.

    What if you just ran this interface and wanted to know what it called the file that was produced? As you might have noticed before, calling the run method returned an object called InterfaceResult that we saved under the variable res. Let’s inspect that object:

    res = skullstrip.run()
    print(res.outputs)
    
    inskull_mask_file = <undefined>
    inskull_mesh_file = <undefined>
    mask_file = <undefined>
    meshfile = <undefined>
    out_file = /home/neuro/workshop_weizmann/workshop/nipype/notebooks/sub-01_ses-test_T1w_brain.nii.gz
    outline_file = <undefined>
    outskin_mask_file = <undefined>
    outskin_mesh_file = <undefined>
    outskull_mask_file = <undefined>
    outskull_mesh_file = <undefined>
    skull_file = <undefined>
    skull_mask_file = <undefined>
    

    We see that four possible files can be generated by BET. Here we ran it in the most simple way possible, so it just generated an out_file, which is the skull-stripped image. Let’s see what happens when we generate a mask. By the way, you can also set inputs at runtime by including them as arguments to the run method:

    res2 = skullstrip.run(mask=True)
    print(res2.outputs)
    
    inskull_mask_file = <undefined>
    inskull_mesh_file = <undefined>
    mask_file = /home/neuro/workshop_weizmann/workshop/nipype/notebooks/sub-01_ses-test_T1w_brain_mask.nii.gz
    meshfile = <undefined>
    out_file = /home/neuro/workshop_weizmann/workshop/nipype/notebooks/sub-01_ses-test_T1w_brain.nii.gz
    outline_file = <undefined>
    outskin_mask_file = <undefined>
    outskin_mesh_file = <undefined>
    outskull_mask_file = <undefined>
    outskull_mesh_file = <undefined>
    skull_file = <undefined>
    skull_mask_file = <undefined>
    

    Nipype knows that if you ask for a mask, BET is going to generate it in a particular way and makes that information available to you.

    Why this is amazing!

    A major motivating objective for Nipype is to streamline the integration of different analysis packages so that you can use the algorithms you feel are best suited to your particular problem.

    Say that you want to use BET, as SPM does not offer a way to create an explicit mask from functional data, but that otherwise, you want your processing to occur in SPM. Although possible to do this in a MATLAB script, it might not be all that clean, particularly if you want your skullstrip to happen in the middle of your workflow (for instance, after realignment). Nipype provides a unified representation of interfaces across analysis packages.

    For more on this, check out the Interfaces and the Workflow tutorial.

    Exercise 1

    Import IsotropicSmooth from nipype.interfaces.fsl and find the FSL command that is being run. What are the mandatory inputs for this interface?

    # write your solution here
    
    from nipype.interfaces.fsl import IsotropicSmooth
    # all this information can be found when we run `help` method. 
    # note that you can either provide `in_file` and `fwhm` or `in_file` and `sigma`
    IsotropicSmooth.help()
    
    Wraps the executable command ``fslmaths``.
    
    Use fslmaths to spatially smooth an image with a gaussian kernel.
    
    Inputs::
    
            [Mandatory]
            fwhm: (a float)
                    fwhm of smoothing kernel [mm]
                    argument: ``-s %.5f``, position: 4
                    mutually_exclusive: sigma
            sigma: (a float)
                    sigma of smoothing kernel [mm]
                    argument: ``-s %.5f``, position: 4
                    mutually_exclusive: fwhm
            in_file: (a pathlike object or string representing an existing file)
                    image to operate on
                    argument: ``%s``, position: 2
    
            [Optional]
            out_file: (a pathlike object or string representing a file)
                    image to write
                    argument: ``%s``, position: -2
            internal_datatype: ('float' or 'char' or 'int' or 'short' or 'double'
                      or 'input')
                    datatype to use for calculations (default is float)
                    argument: ``-dt %s``, position: 1
            output_datatype: ('float' or 'char' or 'int' or 'short' or 'double'
                      or 'input')
                    datatype to use for output (default uses input type)
                    argument: ``-odt %s``, position: -1
            nan2zeros: (a boolean)
                    change NaNs to zeros before doing anything
                    argument: ``-nan``, position: 3
            output_type: ('NIFTI' or 'NIFTI_PAIR' or 'NIFTI_GZ' or
                      'NIFTI_PAIR_GZ')
                    FSL output type
            args: (a string)
                    Additional parameters to the command
                    argument: ``%s``
            environ: (a dictionary with keys which are a bytes or None or a value
                      of class 'str' and with values which are a bytes or None or a
                      value of class 'str', nipype default value: {})
                    Environment variables
    
    Outputs::
    
            out_file: (a pathlike object or string representing an existing file)
                    image written after calculations
    
    References:
    -----------
    BibTeX('@article{JenkinsonBeckmannBehrensWoolrichSmith2012,author={M. Jenkinson, C.F. Beckmann, T.E. Behrens, M.W. Woolrich, and S.M. Smith},title={FSL},journal={NeuroImage},volume={62},pages={782-790},year={2012},}', key='JenkinsonBeckmannBehrensWoolrichSmith2012')
    

    Exercise 2

    Run the IsotropicSmooth for /data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w.nii.gz file with a smoothing kernel 4mm:

    # write your solution here
    
    smoothing = IsotropicSmooth()
    smoothing.inputs.in_file = "/data/ds000114/sub-01/ses-test/anat/sub-01_ses-test_T1w.nii.gz"
    smoothing.inputs.fwhm = 4
    smoothing.inputs.out_file = "/output/T1w_nipype_smooth.nii.gz"
    smoothing.run()
    
    <nipype.interfaces.base.support.InterfaceResult at 0x7f546f0ed4d0>
    

    Exercise 3

    Plot the output of your interface.

    # write your solution here
    
    # we will be using plot_anat from nilearn package
    from nilearn.plotting import plot_anat
    %matplotlib inline
    plot_anat('/output/T1w_nipype_smooth.nii.gz', title='after smoothing',
              display_mode='ortho', dim=-1, draw_cross=False, annotate=False);
    
    ../../_images/basic_interfaces_55_0.png