@@ -490,25 +490,25 @@ def _run_interface(self, runtime):
490
490
data = np .concatenate ([vol .get_data ().reshape (vol .get_shape ()[:3 ] + (- 1 ,)) for vol in vollist ], axis = 3 )
491
491
if isdefined (self .inputs .regress_poly ):
492
492
timepoints = img .get_shape ()[- 1 ]
493
- X = np .ones ((timepoints ,1 ))
493
+ X = np .ones ((timepoints , 1 ))
494
494
for i in range (self .inputs .regress_poly ):
495
- X = np .hstack ((X ,legendre (i + 1 )(np .linspace (- 1 , 1 , timepoints ))[:, None ]))
495
+ X = np .hstack ((X , legendre (i + 1 )(np .linspace (- 1 , 1 , timepoints ))[:, None ]))
496
496
betas = np .dot (np .linalg .pinv (X ), np .rollaxis (data , 3 , 2 ))
497
- datahat = np .rollaxis (np .dot (X [:,1 :],
497
+ datahat = np .rollaxis (np .dot (X [:, 1 :],
498
498
np .rollaxis (betas [1 :, :, :, :], 0 , 3 )),
499
499
0 , 4 )
500
500
data = data - datahat
501
501
img = nb .Nifti1Image (data , img .get_affine (), img .get_header ())
502
- nb .save (img , self ._gen_output_file_name ('detrended' ))
502
+ nb .save (img , self ._gen_output_file_name ('detrended' ))
503
503
meanimg = np .mean (data , axis = 3 )
504
504
stddevimg = np .std (data , axis = 3 )
505
- tsnr = meanimg / stddevimg
505
+ tsnr = meanimg / stddevimg
506
506
img = nb .Nifti1Image (tsnr , img .get_affine (), img .get_header ())
507
- nb .save (img , self ._gen_output_file_name ())
507
+ nb .save (img , self ._gen_output_file_name ())
508
508
img = nb .Nifti1Image (meanimg , img .get_affine (), img .get_header ())
509
- nb .save (img , self ._gen_output_file_name ('mean' ))
509
+ nb .save (img , self ._gen_output_file_name ('mean' ))
510
510
img = nb .Nifti1Image (stddevimg , img .get_affine (), img .get_header ())
511
- nb .save (img , self ._gen_output_file_name ('stddev' ))
511
+ nb .save (img , self ._gen_output_file_name ('stddev' ))
512
512
return runtime
513
513
514
514
def _list_outputs (self ):
0 commit comments