Which 3dTshift method does fmriprep use?
It seems that fmriprep uses the default method in 3dTshift to perform slice time correction based on the call to it from line 121 in stc.py copied below. According to the current docs for 3dTshift, that would be Fourier. Is this correct or does fmriprep somehow use a different method from 3dTshift? slice_timing_correction = pe.Node( TShift( outputtype='NIFTI_GZ', tr=f'{metadata["RepetitionTime"]}s', slice_timing=metadata['SliceTiming'], slice_encoding_direction=metadata.get('SliceEncodingDirection', 'k'), tzero=tzero,