put func-to-anat linear xfm invert node in one place to reduce repeated calculation.

Author: HechengJin0

CPAC/nuisance/nuisance.py

@@ -628,6 +628,7 @@ def create_regressor_workflow(nuisance_selectors,
         'lat_ventricles_mask_file_path',
         'functional_brain_mask_file_path',
         'func_to_anat_linear_xfm_file_path',
+        'anat_to_func_linear_xfm_file_path',
         'mni_to_anat_linear_xfm_file_path',
         'anat_to_mni_initial_xfm_file_path',
         'anat_to_mni_rigid_xfm_file_path',
@@ -667,6 +668,7 @@ def create_regressor_workflow(nuisance_selectors,
 
         "Transformations": {
             "func_to_anat_linear_xfm": (inputspec, "func_to_anat_linear_xfm_file_path"),
+            "anat_to_func_linear_xfm": (inputspec, "anat_to_func_linear_xfm_file_path"),
             "mni_to_anat_linear_xfm": (inputspec, "mni_to_anat_linear_xfm_file_path"),
             "anat_to_mni_initial_xfm": (inputspec, "anat_to_mni_initial_xfm_file_path"),
             "anat_to_mni_rigid_xfm": (inputspec, "anat_to_mni_rigid_xfm_file_path"),

CPAC/nuisance/utils/__init__.py

@@ -324,22 +324,18 @@ def generate_summarize_tissue_mask(nuisance_wf,
             mask_to_epi.inputs.interp = 'nearestneighbour'
 
             if regressor_selector['extraction_resolution'] == "Functional":
-                # invert func2anat matrix to get anat2func_linear_xfm
-                anat2func_linear_xfm = pe.Node(interface=fsl.ConvertXFM(), name='{}_anat2func_linear_xfm'.format(prev_mask_key))
-                anat2func_linear_xfm.inputs.invert_xfm = True
-                nuisance_wf.connect(*(
-                    pipeline_resource_pool['Transformations']['func_to_anat_linear_xfm'] + 
-                    (anat2func_linear_xfm, 'in_file')
-                ))
                 # apply anat2func matrix
                 mask_to_epi.inputs.apply_xfm = True
                 mask_to_epi.inputs.output_type = 'NIFTI_GZ'
                 nuisance_wf.connect(*(
                     pipeline_resource_pool['Functional_mean'] +
                     (mask_to_epi, 'reference')
                 ))
-                nuisance_wf.connect(anat2func_linear_xfm, 'out_file', 
-                                    mask_to_epi, 'in_matrix_file')
+                nuisance_wf.connect(*(
+                    pipeline_resource_pool['Transformations']['anat_to_func_linear_xfm'] + 
+                    (mask_to_epi, 'in_matrix_file')
+                ))
+
             else:
                 resolution = regressor_selector['extraction_resolution']
                 mask_to_epi.inputs.apply_isoxfm = resolution

CPAC/pipeline/cpac_pipeline.py

@@ -2069,6 +2069,18 @@ def build_workflow(subject_id, sub_dict, c, pipeline_name=None, num_ants_cores=1
                     'inputspec.func_to_anat_linear_xfm_file_path'
                 )
 
+                # invert func2anat matrix to get anat2func_linear_xfm
+                anat2func_linear_xfm = pe.Node(interface=fsl.ConvertXFM(), name='anat_to_func_linear_xfm_{0}_{1}'.format(regressors_selector_i, num_strat))
+                anat2func_linear_xfm.inputs.invert_xfm = True
+                workflow.connect(
+                    node, out_file,
+                    anat2func_linear_xfm, 'in_file')
+                workflow.connect(
+                    anat2func_linear_xfm, 'out_file', 
+                    regressor_workflow,
+                    'inputspec.anat_to_func_linear_xfm_file_path'
+                )
+
                 node, out_file = new_strat.get_leaf_properties()
                 workflow.connect(
                     node, out_file,