You could do this with your original EPI datasets:
3dcalc -a epi_r1+orig -expr t -prefix tt_epi_time
This replaces every voxel with the time it was acquired because 't' is a reserved variable for that purpose in 3dcalc. Then you could just apply the same warping to cardinal with 3dWarp on that time dataset. There is the question of what to do with interpolation there. It may be most appropriate to use Nearest Neighbor interpolation rather than computing an interpolated value.
Edited 1 time(s). Last edit at 10/17/2018 07:19PM by Daniel Glen.