Hello,
I have noticed some problem of getting correct slice timing of multiband fMRI data generated by the custom sequence built by CMRR of University of Minnesota. For reference of their slice timing, here is a
link. We use Siemens Trio scanner with 32 channel coil.
The command I used for to3d was like this:
to3d -epan -prefix empty_run.01 -assume_dicom_mosaic -oblique_origin -time:zt 58 274 2000 FROM_IMAGE -anatparent inplane.02+orig. -zpad 1 ...
Running "3dinfo -slice_timing" for the generated empty_run.01+orig gave me output like this:
0.000000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|-0.666000|0.000000
Running it for any specific sub-brick gave me this:
0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000|0.000000
But if I run "dicom_hdr -slice_times" for any dicom file other than the first volume, it gives me correct timing that is consistent with what the document (linked above) says. Notice that the document mentioned a bug that the timing stored for the first volume was always wrong (simply speaking, all time stamps stored for odd slices are larger than the actual TR).
So I am not sure if this problem is because to3d reads the slice timing information from the first dicom image and gets confused? If so, can you make it read from the last image?
The forum does not allow me to upload dicom file. If you need sample dicom file please let me know how I can upload.
Thanks!
Mingbo