@@ -266,7 +266,7 @@ def run_dmriprep_pe(dwi_file, dwi_file_AP, dwi_file_PA, bvec_file, bval_file,
266
266
eddy .inputs .repol = True
267
267
eddy .inputs .niter = 1 # TODO: change back to 5 when running for real
268
268
269
- def drop_eddy_outliers (outlier_report , threshold ):
269
+ def drop_outliers_fn (outlier_report , threshold ):
270
270
"""Get list of scans that exceed threshold for number of outliers
271
271
272
272
Parameters
@@ -310,18 +310,18 @@ def num_outliers(scan, outliers):
310
310
311
311
return drop_scans
312
312
313
- drop_outliers = pe .Node (niu .Function (
313
+ drop_outliers_node = pe .Node (niu .Function (
314
314
input_names = ["outlier_report" , "threshold" ],
315
315
output_names = ["drop_scans" ],
316
- function = drop_eddy_outliers ),
317
- name = "drop_outliers "
316
+ function = drop_outliers_fn ),
317
+ name = "drop_outliers_node "
318
318
)
319
319
320
- drop_outliers .inputs .threshold = 1
320
+ drop_outliers_node .inputs .threshold = 1
321
321
wf .connect (prep , "fsl_eddy.out_outlier_report" ,
322
- drop_outliers , "outlier_report" )
322
+ drop_outliers_node , "outlier_report" )
323
323
324
- def save_outlier_list (drop_scans ):
324
+ def save_outlier_list_fn (drop_scans ):
325
325
"""Save list of outlier scans to file
326
326
327
327
Parameters
@@ -340,14 +340,15 @@ def save_outlier_list(drop_scans):
340
340
np .savetxt (outpath , drop_scans , fmt = "%d" )
341
341
return outpath
342
342
343
- save_drop_scans = pe .Node (niu .Function (
343
+ save_outlier_list_node = pe .Node (niu .Function (
344
344
input_names = ["drop_scans" ],
345
345
output_names = ["outpath" ],
346
- function = save_outlier_list ),
347
- name = "save_drop_scans "
346
+ function = save_outlier_list_fn ),
347
+ name = "save_outlier_list_node "
348
348
)
349
349
350
- wf .connect (drop_outliers , "drop_scans" , save_drop_scans , "drop_scans" )
350
+ wf .connect (drop_outliers_node , "drop_scans" ,
351
+ save_outlier_list_node , "drop_scans" )
351
352
352
353
merge = pe .Node (fsl .Merge (dimension = 't' ), name = "mergeAPPA" )
353
354
merge .inputs .in_files = [dwi_file_AP , dwi_file_PA ]
@@ -468,7 +469,7 @@ def get_orig(subjects_dir, sub):
468
469
wf .connect (prep , "fsl_eddy.out_shell_alignment_parameters" ,
469
470
datasink , "dmriprep.qc.@eddyparamsshellalign" )
470
471
471
- wf .connect (save_drop_scans , "outpath" , datasink , "dmriprep.qc.@droppedscans" )
472
+ wf .connect (save_outlier_list_node , "outpath" , datasink , "dmriprep.qc.@droppedscans" )
472
473
473
474
wf .connect (get_tensor , "out_file" , datasink , "dmriprep.dti.@tensor" )
474
475
wf .connect (get_tensor , "fa_file" , datasink , "dmriprep.dti.@fa" )
0 commit comments