|
| 1 | +"""Workflows to resample data.""" |
| 2 | + |
| 3 | +from nipype.interfaces import utility as niu |
| 4 | +from nipype.pipeline import engine as pe |
| 5 | + |
| 6 | + |
| 7 | +def init_resample_raw_wf(bold_file, functional_cache): |
| 8 | + """Resample raw BOLD data to MNI152NLin6Asym:res-2mm space.""" |
| 9 | + from fmriprep.workflows.bold.stc import init_bold_stc_wf |
| 10 | + from niworkflows.engine.workflows import LiterateWorkflow as Workflow |
| 11 | + |
| 12 | + from fmripost_aroma.interfaces.resampler import Resampler |
| 13 | + |
| 14 | + workflow = Workflow(name="resample_raw_wf") |
| 15 | + |
| 16 | + inputnode = pe.Node( |
| 17 | + niu.IdentityInterface(fields=["bold_file", "mask_file"]), |
| 18 | + name="inputnode", |
| 19 | + ) |
| 20 | + inputnode.inputs.bold_file = bold_file |
| 21 | + inputnode.inputs.mask_file = functional_cache["bold_mask"] |
| 22 | + |
| 23 | + outputnode = pe.Node( |
| 24 | + niu.IdentityInterface(fields=["bold_std", "bold_mask_std"]), |
| 25 | + name="outputnode", |
| 26 | + ) |
| 27 | + |
| 28 | + stc_wf = init_bold_stc_wf(name="resample_stc_wf") |
| 29 | + workflow.connect([ |
| 30 | + (inputnode, stc_wf, [ |
| 31 | + ('bold_file', 'inputnode.bold_file'), |
| 32 | + ('mask_file', 'inputnode.mask_file'), |
| 33 | + ]), |
| 34 | + ]) # fmt:skip |
| 35 | + |
| 36 | + resample_bold = pe.Node( |
| 37 | + Resampler(space="MNI152NLin6Asym", resolution="2"), |
| 38 | + name="resample_bold", |
| 39 | + ) |
| 40 | + workflow.connect([ |
| 41 | + (stc_wf, resample_bold, [('outputnode.bold_file', 'bold_file')]), |
| 42 | + (resample_bold, outputnode, [('output_file', 'bold_std')]), |
| 43 | + ]) # fmt:skip |
| 44 | + |
| 45 | + resample_bold_mask = pe.Node( |
| 46 | + Resampler(space="MNI152NLin6Asym", resolution="2"), |
| 47 | + name="resample_bold_mask", |
| 48 | + ) |
| 49 | + workflow.connect([ |
| 50 | + (inputnode, resample_bold_mask, [('mask_file', 'bold_file')]), |
| 51 | + (resample_bold_mask, outputnode, [('output_file', 'bold_mask_std')]), |
| 52 | + ]) # fmt:skip |
| 53 | + |
| 54 | + return workflow |
0 commit comments