Skip to content

Commit

Permalink
ENH: Remove ICA-AROMA from workflow and docs
Browse files Browse the repository at this point in the history
  • Loading branch information
effigies committed Mar 13, 2023
1 parent d14b811 commit f8f64bb
Show file tree
Hide file tree
Showing 26 changed files with 39 additions and 794 deletions.
2 changes: 0 additions & 2 deletions .circleci/config.yml
Original file line number Diff line number Diff line change
Expand Up @@ -551,7 +551,6 @@ jobs:
mkdir -p /tmp/${DATASET}/derivatives_partial
sudo cp -a /tmp/${DATASET}/work /tmp/${DATASET}/work_partial
sudo cp -a /tmp/${DATASET}/work /tmp/${DATASET}/work_bids
sudo rm -rf /tmp/${DATASET}/work_partial/fmriprep_*_wf/single_subject_01_wf/func_preproc_task_mixedgamblestask_run_02_wf/ica_aroma_wf
# Nipype sometimes fails to pick up when the base directory changes
# This is a cheap workflow, so let's not fuss
sudo rm -rf /tmp/${DATASET}/work /tmp/${DATASET}/work_bids/fmriprep_*_wf/single_subject_01_wf/anat_preproc_wf/surface_recon_wf/gifti_surface_wf || true
Expand Down Expand Up @@ -608,7 +607,6 @@ jobs:
--output-layout legacy \
--sloppy --write-graph --use-syn-sdc --mem-mb 14336 \
--output-spaces MNI152NLin2009cAsym fsaverage5 fsnative MNI152NLin6Asym anat \
--aroma-melodic-dimensionality 2 --use-aroma \
--nthreads 4 --cifti-output --project-goodvoxels -vv
- store_artifacts:
path: /tmp/ds005/derivatives_partial
Expand Down
4 changes: 0 additions & 4 deletions .circleci/ds005_legacy_partial_fasttrack_outputs.txt
Original file line number Diff line number Diff line change
Expand Up @@ -17,11 +17,9 @@ fmriprep/sub-01/fmap/sub-01_run-02_fmapid-auto00000_desc-magnitude_fieldmap.nii.
fmriprep/sub-01/fmap/sub-01_run-02_fmapid-auto00000_desc-preproc_fieldmap.json
fmriprep/sub-01/fmap/sub-01_run-02_fmapid-auto00000_desc-preproc_fieldmap.nii.gz
fmriprep/sub-01/func
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_AROMAnoiseICs.csv
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_boldref.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_desc-confounds_timeseries.json
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_desc-confounds_timeseries.tsv
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_desc-MELODIC_mixing.tsv
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_from-scanner_to-boldref_mode-image_xfm.txt
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_from-scanner_to-T1w_mode-image_xfm.txt
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_from-T1w_to-scanner_mode-image_xfm.txt
Expand Down Expand Up @@ -49,8 +47,6 @@ fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_d
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-brain_mask.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-preproc_bold.json
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-preproc_bold.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-smoothAROMAnonaggr_bold.json
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-smoothAROMAnonaggr_bold.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-T1w_boldref.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-T1w_desc-aparcaseg_dseg.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-T1w_desc-aseg_dseg.nii.gz
Expand Down
4 changes: 0 additions & 4 deletions .circleci/ds005_legacy_partial_outputs.txt
Original file line number Diff line number Diff line change
Expand Up @@ -71,11 +71,9 @@ fmriprep/sub-01/fmap/sub-01_run-02_fmapid-auto00000_desc-magnitude_fieldmap.nii.
fmriprep/sub-01/fmap/sub-01_run-02_fmapid-auto00000_desc-preproc_fieldmap.json
fmriprep/sub-01/fmap/sub-01_run-02_fmapid-auto00000_desc-preproc_fieldmap.nii.gz
fmriprep/sub-01/func
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_AROMAnoiseICs.csv
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_boldref.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_desc-confounds_timeseries.json
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_desc-confounds_timeseries.tsv
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_desc-MELODIC_mixing.tsv
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_from-scanner_to-boldref_mode-image_xfm.txt
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_from-scanner_to-T1w_mode-image_xfm.txt
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_from-T1w_to-scanner_mode-image_xfm.txt
Expand Down Expand Up @@ -103,8 +101,6 @@ fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_d
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-brain_mask.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-preproc_bold.json
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-preproc_bold.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-smoothAROMAnonaggr_bold.json
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-MNI152NLin6Asym_desc-smoothAROMAnonaggr_bold.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-T1w_boldref.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-T1w_desc-aparcaseg_dseg.nii.gz
fmriprep/sub-01/func/sub-01_task-mixedgamblestask_run-02_space-T1w_desc-aseg_dseg.nii.gz
Expand Down
3 changes: 1 addition & 2 deletions .vscode/launch.json
Original file line number Diff line number Diff line change
Expand Up @@ -34,7 +34,6 @@
"/tmp/ds005/derivatives",
"participant",
"--sloppy", "--write-graph",
"--use-aroma",
"--skull-strip-template", "OASIS30ANTs:res-1",
"--output-space", "T1w", "template", "fsaverage5", "fsnative",
"--template-resampling-grid", "native",
Expand Down Expand Up @@ -98,4 +97,4 @@
"justMyCode": false
}
]
}
}
8 changes: 0 additions & 8 deletions Dockerfile
Original file line number Diff line number Diff line change
Expand Up @@ -223,14 +223,6 @@ WORKDIR $ANTSPATH
RUN curl -sSL "https://dl.dropbox.com/s/gwf51ykkk5bifyj/ants-Linux-centos6_x86_64-v2.3.4.tar.gz" \
| tar -xzC $ANTSPATH --strip-components 1

# Installing and setting up ICA_AROMA
WORKDIR /opt/ICA-AROMA
RUN curl -sSL "https://github.com/oesteban/ICA-AROMA/archive/v0.4.5.tar.gz" \
| tar -xzC /opt/ICA-AROMA --strip-components 1 && \
chmod +x /opt/ICA-AROMA/ICA_AROMA.py
ENV PATH="/opt/ICA-AROMA:$PATH" \
AROMA_VERSION="0.4.5"

WORKDIR /opt
RUN curl -sSLO https://www.humanconnectome.org/storage/app/media/workbench/workbench-linux64-v1.5.0.zip && \
unzip workbench-linux64-v1.5.0.zip && \
Expand Down
1 change: 0 additions & 1 deletion REFERENCES.md
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,6 @@
| MCFLIRT | https://doi.org/10.1006/nimg.2002.1132 | https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/MCFLIRT |
| SUSAN | https://doi.org/10.1023/A:1007963824710 | https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/SUSAN |
| MELODIC | | https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/MELODIC |
| ICA-AROMA | http://www.sciencedirect.com/science/article/pii/S1053811915001822 | https://github.com/rhr-pruim/ICA-AROMA/ |
| PRELUDE & FUGUE | https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FUGUE/Guide | https://nipype.readthedocs.io/en/latest/interfaces/generated/workflows.dmri/fsl.utils.html#cleanup-edge-pipeline |
| **AFNI** | | https://doi.org/10.1006/cbmr.1996.0014; https://doi.org/10.1016/j.neuroimage.2011.08.056
| 3dvolreg | | https://afni.nimh.nih.gov/pub/dist/doc/program_help/3dvolreg.html |
Expand Down
2 changes: 1 addition & 1 deletion docs/_static/sbatch.slurm
Original file line number Diff line number Diff line change
Expand Up @@ -40,7 +40,7 @@ subject=$( sed -n -E "$((${SLURM_ARRAY_TASK_ID} + 1))s/sub-(\S*)\>.*/\1/gp" ${BI
find ${LOCAL_FREESURFER_DIR}/sub-$subject/ -name "*IsRunning*" -type f -delete

# Compose the command line
cmd="${SINGULARITY_CMD} /data /data/${DERIVS_DIR} participant --participant-label $subject -w /work/ -vv --omp-nthreads 8 --nthreads 12 --mem_mb 30000 --output-spaces MNI152NLin2009cAsym:res-2 anat fsnative fsaverage5 --use-aroma --fs-subjects-dir /fsdir"
cmd="${SINGULARITY_CMD} /data /data/${DERIVS_DIR} participant --participant-label $subject -w /work/ -vv --omp-nthreads 8 --nthreads 12 --mem_mb 30000 --output-spaces MNI152NLin2009cAsym:res-2 anat fsnative fsaverage5 --fs-subjects-dir /fsdir"

# Setup done, run the command
echo Running task ${SLURM_ARRAY_TASK_ID}
Expand Down
2 changes: 1 addition & 1 deletion docs/faq.rst
Original file line number Diff line number Diff line change
Expand Up @@ -202,7 +202,7 @@ In addition to the ``--output-spaces`` that you specify, *fMRIPrep* will
internally require the ``MNI152NLin2009cAsym`` template.
If the ``--skull-strip-template`` option is not set, then ``OASIS30ANTs``
will be used.
Finally, both the ``--cifti-output`` and ``--use-aroma`` arguments require ``MNI152NLin6Asym``.
Finally, the ``--cifti-output`` argument requires ``MNI152NLin6Asym``.
To do so, follow the next steps.

1. By default, a mirror of *TemplateFlow* to store the resources will be
Expand Down
1 change: 0 additions & 1 deletion docs/installation.rst
Original file line number Diff line number Diff line change
Expand Up @@ -80,7 +80,6 @@ the ``fmriprep`` package:
- AFNI_ (version 22.3.06)
- `C3D <https://sourceforge.net/projects/c3d/>`_ (version 1.3.0)
- FreeSurfer_ (version 7.3.2)
- `ICA-AROMA <https://github.com/maartenmennes/ICA-AROMA/archive/v0.4.5.tar.gz>`_ (version 0.4.5)
- `bids-validator <https://github.com/bids-standard/bids-validator>`_ (version 1.8.0)
- `connectome-workbench <https://www.humanconnectome.org/software/connectome-workbench>`_ (version 1.5.0)

Expand Down
53 changes: 7 additions & 46 deletions docs/outputs.rst
Original file line number Diff line number Diff line change
Expand Up @@ -30,10 +30,9 @@ upcoming `BEP 011`_ and `BEP 012`_).
.. important::
In order to remain agnostic to any possible subsequent analysis,
*fMRIPrep* does not perform any denoising (e.g., spatial smoothing) itself.
There are two exceptions to this principle (described in their corresponding
sections below):
There are exceptions to this principle (described in its corresponding
section below):

- ICA-AROMA's *non-aggressive* denoised outputs, and
- CompCor regressors, which are calculated after temporal high-pass filtering.

Layout
Expand Down Expand Up @@ -228,19 +227,11 @@ These :abbr:`TSV (tab-separated values)` tables look like the example below,
where each row of the file corresponds to one time point found in the
corresponding :abbr:`BOLD (blood-oxygen level dependent)` time series::

csf white_matter global_signal std_dvars dvars framewise_displacement t_comp_cor_00 t_comp_cor_01 t_comp_cor_02 t_comp_cor_03 t_comp_cor_04 t_comp_cor_05 a_comp_cor_00 a_comp_cor_01 a_comp_cor_02 a_comp_cor_03 a_comp_cor_04 a_comp_cor_05 non_steady_state_outlier00 trans_x trans_y trans_z rot_x rot_y rot_z aroma_motion_02 aroma_motion_04
682.75275 0.0 491.64752000000004 n/a n/a n/a 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 -0.00017029 -0.0 0.0 0.0 0.0
669.14166 0.0 489.4421 1.168398 17.575331 0.07211929999999998 -0.4506846719 0.1191909139 -0.0945884724 0.1542023065 -0.2302324641 0.0838194238 -0.032426848599999995 0.4284323184 -0.5809158299 0.1382414008 -0.1203486637 0.3783661265 0.0 0.0 0.0207752 0.0463124 -0.000270924 -0.0 0.0 -2.402958171 -0.7574011893
665.3969 0.0 488.03 1.085204 16.323903999999995 0.0348966 0.010819676200000001 0.0651895837 -0.09556632150000001 -0.033148835 -0.4768871111 0.20641088559999998 0.2818768463 0.4303863764 0.41323714850000004 -0.2115232212 -0.0037154909000000004 0.10636180070000001 0.0 0.0 0.0 0.0457372 0.0 -0.0 0.0 -1.341359143 0.1636017242
662.82715 0.0 487.37302 1.01591 15.281561 0.0333937 0.3328022893 -0.2220965269 -0.0912891436 0.2326688125 0.279138129 -0.111878887 0.16901660629999998 0.0550480212 0.1798747037 -0.25383302620000003 0.1646403629 0.3953613889 0.0 0.010164 -0.0103568 0.0424513 0.0 -0.0 0.00019174 -0.1554834655 0.6451987913

Finally, if ICA-AROMA is used, the MELODIC mixing matrix and the components classified as noise
are saved::

sub-<subject_label>/
func/
sub-<subject_label>_[specifiers]_AROMAnoiseICs.csv
sub-<subject_label>_[specifiers]_desc-MELODIC_mixing.tsv
csf white_matter global_signal std_dvars dvars framewise_displacement t_comp_cor_00 t_comp_cor_01 t_comp_cor_02 t_comp_cor_03 t_comp_cor_04 t_comp_cor_05 a_comp_cor_00 a_comp_cor_01 a_comp_cor_02 a_comp_cor_03 a_comp_cor_04 a_comp_cor_05 non_steady_state_outlier00 trans_x trans_y trans_z rot_x rot_y rot_z
682.75275 0.0 491.64752000000004 n/a n/a n/a 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 -0.00017029 -0.0 0.0
669.14166 0.0 489.4421 1.168398 17.575331 0.07211929999999998 -0.4506846719 0.1191909139 -0.0945884724 0.1542023065 -0.2302324641 0.0838194238 -0.032426848599999995 0.4284323184 -0.5809158299 0.1382414008 -0.1203486637 0.3783661265 0.0 0.0 0.0207752 0.0463124 -0.000270924 -0.0 0.0
665.3969 0.0 488.03 1.085204 16.323903999999995 0.0348966 0.010819676200000001 0.0651895837 -0.09556632150000001 -0.033148835 -0.4768871111 0.20641088559999998 0.2818768463 0.4303863764 0.41323714850000004 -0.2115232212 -0.0037154909000000004 0.10636180070000001 0.0 0.0 0.0 0.0457372 0.0 -0.0 0.0
662.82715 0.0 487.37302 1.01591 15.281561 0.0333937 0.3328022893 -0.2220965269 -0.0912891436 0.2326688125 0.279138129 -0.111878887 0.16901660629999998 0.0550480212 0.1798747037 -0.25383302620000003 0.1646403629 0.3953613889 0.0 0.010164 -0.0103568 0.0424513 0.0 -0.0 0.00019174

**Multi-echo derivatives**.
For multi-echo datasets, the output ``_bold`` series are "optimally combined" by
Expand Down Expand Up @@ -529,31 +520,6 @@ the brain's outer edge (*crown*) mask.
The procedure essentially follows the initial proposal of the approach by Patriat et al.
[Patriat2017]_ and is described in our ISMRM abstract [Provins2022]_.

**AROMA confounds**.
:abbr:`AROMA (Automatic Removal Of Motion Artifacts)` is an :abbr:`ICA (independent components analysis)`
based procedure to identify confounding time series related to head-motion [Prium2015]_.
ICA-AROMA can be enabled with the flag ``--use-aroma``.

- ``aroma_motion_XX`` - the motion-related components identified by ICA-AROMA.

.. danger::
If you are already using AROMA-cleaned data (``~desc-smoothAROMAnonaggr_bold.nii.gz``),
do not include ICA-AROMA confounds during your design specification or denoising procedure.

Additionally, as per [Hallquist2013]_ and [Lindquist2019]_, when using AROMA-cleaned data
most of the confound regressors should be recalculated (this feature is a work-in-progress,
follow up on `#1905 <https://github.com/nipreps/fmriprep/issues/1905>`__).
Surprisingly, `our simulations
<https://github.com/nipreps/fmriprep-notebooks/blob/9933a628dfb759dc73e61701c144d67898b92de0/05%20-%20Discussion%20AROMA%20confounds%20-%20issue-817%20%5BJ.%20Kent%5D.ipynb>`__
(with thanks to JD. Kent) suggest that using the confounds as currently calculated by
*fMRIPrep* --before denoising-- would be just fine.

.. caution::
*Nonsteady-states* (or *dummy scans*) in the beginning of every run
are dropped **before** ICA-AROMA is performed.
Therefore, any subsequent analysis of ICA-AROMA outputs must drop the same
number of *nonsteady-states*.

Confounds and "carpet"-plot on the visual reports
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
The visual reports provide several sections per task and run to aid designing
Expand Down Expand Up @@ -698,11 +664,6 @@ the following invocation::
.. [Power2016] Power JD, A simple but useful way to assess fMRI scan qualities. NeuroImage. 2016.
doi:`10.1016/j.neuroimage.2016.08.009 <http://doi.org/10.1016/j.neuroimage.2016.08.009>`_
.. [Prium2015] Pruim RHR, Mennes M, van Rooij D, Llera A, Buitelaar JK, Beckmann CF.
ICA-AROMA: A robust ICA-based strategy for removing motion artifacts from fMRI data.
Neuroimage. 2015 May 15;112:267–77.
doi:`10.1016/j.neuroimage.2015.02.064 <https://doi.org/10.1016/j.neuroimage.2015.02.064>`_.
.. [Provins2022] Provins C et al., Quality control and nuisance regression of fMRI, looking out
where signal should not be found. Proc. Intl. Soc. Mag. Reson. Med. 31, London (UK). 2022
doi:`10.31219/osf.io/hz52v <https://doi.org/10.31219/osf.io/hz52v>`_.
Expand Down
6 changes: 3 additions & 3 deletions docs/spaces.rst
Original file line number Diff line number Diff line change
Expand Up @@ -138,11 +138,11 @@ Modifiers are not allowed when providing nonstandard spaces.

Preprocessing blocks depending on standard templates
""""""""""""""""""""""""""""""""""""""""""""""""""""
Some modules of the pipeline (e.g., the ICA-AROMA denoising, the generation of
HCP compatible *grayordinates* files, or the *fieldmap-less* distortion correction)
Some modules of the pipeline (e.g., the generation of HCP compatible
*grayordinates* files, or the *fieldmap-less* distortion correction)
operate in specific template spaces.
When selecting those modules to be included (using any of the following flags:
``--use-aroma``, ``--cifti-outputs``, ``--use-syn-sdc``) will modify the list of
``--cifti-outputs``, ``--use-syn-sdc``) will modify the list of
*internal* spaces to include the space identifiers they require, should the
identifier not be found within the ``--output-spaces`` list already.
In other words, running *fMRIPrep* with ``--output-spaces MNI152NLin6Asym:res-2
Expand Down
2 changes: 1 addition & 1 deletion docs/usage.rst
Original file line number Diff line number Diff line change
Expand Up @@ -180,7 +180,7 @@ Because processing this emerging type of datasets (*densely sampled neuroimaging

.. attention::
When the intention is to combine the *anatomical fast-track* with some advanced options that involve
standard spaces (e.g., ``--use-aroma`` or ``--cifti-output``), please make sure you include the
standard spaces (e.g., ``--cifti-output``), please make sure you include the
``MNI152NLin6Asym`` space to the ``--output-spaces`` list in the first invocation of *fMRIPrep*
(or *sMRIPrep*) from which the results are to be reused.
Otherwise, a warning message indicating that *fMRIPrep*'s expectations were not met will be issued,
Expand Down
52 changes: 2 additions & 50 deletions docs/workflows.rst
Original file line number Diff line number Diff line change
Expand Up @@ -567,56 +567,8 @@ segmentation, the `discover_wf` sub-workflow calculates potential
confounds per volume.

Calculated confounds include the mean global signal, mean tissue class signal,
tCompCor, aCompCor, Frame-wise Displacement, 6 motion parameters, DVARS, spike regressors,
and, if the ``--use-aroma`` flag is enabled, the noise components identified by ICA-AROMA
(those to be removed by the "aggressive" denoising strategy).
Particular details about ICA-AROMA are given below.

ICA-AROMA
~~~~~~~~~
:py:func:`~fmriprep.workflows.bold.confounds.init_ica_aroma_wf`

ICA-AROMA denoising is performed in ``MNI152NLin6Asym`` space, which is automatically
added to the list of ``--output-spaces`` if it was not already requested by the user.
The number of ICA-AROMA components depends on a dimensionality estimate made by
FSL MELODIC.
For datasets with a very short TR and a large number of timepoints, this may result
in an unusually high number of components.
By default, dimensionality is limited to a maximum of 200 components.
To override this upper limit one may specify the number of components to be extracted
with ``--aroma-melodic-dimensionality``.
Further details on the implementation are given within the workflow generation
function (:py:func:`~fmriprep.workflows.bold.confounds.init_ica_aroma_wf`).

*Note*: *non*-aggressive AROMA denoising is a fundamentally different procedure
from its "aggressive" counterpart and cannot be performed only by using a set of noise
regressors (a separate GLM with both noise and signal regressors needs to be used).
Therefore instead of regressors, *fMRIPrep* produces *non*-aggressive denoised 4D NIFTI
files in the MNI space:

``*space-MNI152NLin6Asym_desc-smoothAROMAnonaggr_bold.nii.gz``

Additionally, the MELODIC mix and noise component indices will
be generated, so non-aggressive denoising can be manually performed in the T1w space with ``fsl_regfilt``, *e.g.*::

fsl_regfilt -i sub-<subject_label>_task-<task_id>_space-T1w_desc-preproc_bold.nii.gz \
-f $(cat sub-<subject_label>_task-<task_id>_AROMAnoiseICs.csv) \
-d sub-<subject_label>_task-<task_id>_desc-MELODIC_mixing.tsv \
-o sub-<subject_label>_task-<task_id>_space-T1w_desc-AROMAnonaggr_bold.nii.gz

*Note*: The non-steady state volumes are removed for the determination of components in melodic.
Therefore ``*MELODIC_mixing.tsv`` may have zero padded rows to account for the volumes not used in
melodic's estimation of components.

A visualization of the AROMA component classification is also included in the HTML reports.

.. figure:: _static/aroma.svg

Maps created with maximum intensity projection (glass brain) with a black
brain outline.
Right hand side of each map: time series (top in seconds),
frequency spectrum (bottom in Hertz).
Components classified as signal in green; noise in red.
tCompCor, aCompCor, Frame-wise Displacement, 6 motion parameters, DVARS, and
spike regressors.

.. _bold_t2s:

Expand Down
Loading

0 comments on commit f8f64bb

Please sign in to comment.