diff --git a/slicetimer.cc b/slicetimer.cc index 98535e195836593ef9928194bbfc46a240b6987e..45937bfbd35791b9ff785e12bf6d79d857eba34a 100644 --- a/slicetimer.cc +++ b/slicetimer.cc @@ -2,7 +2,7 @@ Peter Bannister and Stephen Smith, FMRIB Image Analysis Group - Copyright (C) 2001-2002 University of Oxford */ + Copyright (C) 2001-2003 University of Oxford */ /* CCOPYRIGHT */ @@ -88,10 +88,10 @@ int do_slice_correction() no_volumes = timeseries.tsize(); repeat_time = timeseries.tdim(); if (repeat_time ==0){ - cerr << "Zero TR in file header - fixing ... " ; + // cerr << "Zero TR in file header - fixing ... " ; repeat_time = repeat.value(); } - cerr << "TR = " << repeat_time << endl; + // cerr << "TR = " << repeat_time << endl; slice_spacing = repeat_time / no_slices; if (direction.value() == 1) timeseries. swapdimensions(3,2,1); // Flip z and x @@ -139,8 +139,8 @@ int do_slice_correction() ColumnVector voxeltimeseries = timeseries.voxelts(x_pos,y_pos,slice-1); ColumnVector interpseries = voxeltimeseries; for (int time_step=1; time_step <= no_volumes; time_step++){ - // interpseries(time_step) = interpolate_1d(voxeltimeseries, time_step - offset + recenter); - interpseries(time_step) = kernelinterpolation_1d(voxeltimeseries, time_step - offset + recenter, userkernel, 7); + // interpseries(time_step) = interpolate_1d(voxeltimeseries, time_step - offset - recenter); + interpseries(time_step) = kernelinterpolation_1d(voxeltimeseries, time_step - offset - recenter, userkernel, 7); } timeseries.setvoxelts(interpseries,x_pos,y_pos,slice-1); }