@@ -471,14 +471,11 @@ def _run_interface(self, runtime):
471
471
self .inputs .pre_filter , degree , self .inputs .high_pass_cutoff , TR )
472
472
473
473
if skip_vols :
474
- old_comp , old_basis = components , filter_basis
474
+ old_comp = components
475
475
nrows = skip_vols + components .shape [0 ]
476
476
components = np .zeros ((nrows , components .shape [1 ]),
477
477
dtype = components .dtype )
478
478
components [skip_vols :] = old_comp
479
- filter_basis = np .zeros ((nrows , filter_basis .shape [1 ]),
480
- dtype = filter_basis .dtype )
481
- filter_basis [skip_vols :] = old_basis
482
479
483
480
components_file = os .path .join (os .getcwd (), self .inputs .components_file )
484
481
np .savetxt (components_file , components , fmt = b"%.10f" , delimiter = '\t ' ,
@@ -491,9 +488,12 @@ def _run_interface(self, runtime):
491
488
ncols = filter_basis .shape [1 ] if filter_basis .size > 0 else 0
492
489
header = ['{}{:02d}' .format (ftype , i ) for i in range (ncols )]
493
490
if skip_vols :
494
- ss_cols = np .eye (components .shape [0 ], skip_vols ,
495
- dtype = filter_basis .dtype )
496
- filter_basis = np .hstack ((filter_basis , ss_cols ))
491
+ old_basis = filter_basis
492
+ nrows = filter_basis .shape [0 ] if filter_basis .size > 0 else 0
493
+ filter_basis = np .zeros ((nrows + skip_vols , ncols + skip_vols ),
494
+ dtype = filter_basis .dtype )
495
+ filter_basis [skip_vols :] = old_basis
496
+ filter_basis [:skip_vols , - skip_vols :] = np .eye (skip_vols )
497
497
header .extend (['SteadyState{:02d}' .format (i )
498
498
for i in range (skip_vols )])
499
499
np .savetxt (pre_filter_file , filter_basis , fmt = b'%.10f' ,
0 commit comments