Skip to content
Snippets Groups Projects
Commit 06e331c9 authored by Stephen Smith's avatar Stephen Smith
Browse files

err....fixed the timings properly this time...maybe?
parent 4e8eebde
No related branches found
No related tags found
No related merge requests found
......@@ -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);
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment