Reference slice in slice timing in SPM

1 min read

I saw at least two web pages saying that the reference slice can be chosen during slice timing correction. For example, I disagree.

Let’s take the first volume as example. Assume TR=2 and number of slices are 31. During data analysis (GLM or ROI analysis), we assume that the first volume is at time 0. In reality, the first acquired slice in this volume is at time 0, the 2nd is at time 2/31, the 3rd 4/31, and so on. The 31st slice is acquired at time 60/31. So except for the first acquired slice, the timing of all other slices should be corrected (or shifted) to time 0. That means, only the 1st acquired slice can be the “reference slice”. Reference slice is not a choice in data preprocessing step.

In SPM’s convention, slice #1 is the most bottom slice which may or may not be the first acquired slice depending on the slice order. Assume the slice order is descending (meaning data are acquired from top to bottom). Then the first acquired slice is slice #31; and the last acquired is slice #1. So in this case, the reference slice should be 31, and the slice order is 31, 30, 29, …, 1.

Choosing a wrong reference slice can affect the timing of your signal. Below I plotted the fMRI signal time series in motor cortex in a finger tapping experiment (tapping for 15s, rest 20s). The top plot is for reference slice = 1, the middle plot is for reference slice = 31. For the bottom plot, I did not do slice timing correction at all. I simply plot the time series in this region (which is at slice #23, or 0.5 seconds away from #31). The vertical line indicates the onset of tapping instruction. Click to zoom the figure.

You can see from the above plot that the rising timing of the signal can be off by 2s if difference reference slice is chosen. Does this matter? It depending on if timing is important to your experiment. For a block design, it’s probably fine; but for event design, or when the question you address is about timing, then this is important.

(The only preprocessing done on the images are reslicing to the first volume. No normalizing, smoothing etc ).

Subscribe to my blog
Don't want to miss new papers in your field? Check out Stork we developed:

fNIRS Journal Club

We have an online fNIRS Journal Club in the end of each month. For the last journal club on March 21st, the presenter was...
Xu Cui
22 sec read


想了解近红外fnirs技术的同学,可以考虑参加我们周六中午的知识竞猜活动。因为心仪科技的慷慨赞助,我们今年的红包比往年大几倍,从5元涨到20元,最后还有大题,奖励100元的大红包。欢迎大家踊跃参加,也欢迎邀请朋友来。 时间:北京时间本周六(2020/2/15)中午12点,在“NIRS知识红包竞猜”群。 不了解心仪的同学,可以趁此机会了解一下: 上海心仪电子科技有限公司是一家专业从事心理学、认知神经科学和脑科学产品技术服务商,目前公司代理了德国ANT便携式脑电、美国ISS公司的Imagent近红外、韩国OBELAB便携式近红外、MAG&More公司的PowerMag系列经颅磁,为脑功能成像和调控提供全套解决方案。
Xu Cui
2 sec read

After fNIRS experiment …

Xu Cui
13 sec read

7 Replies to “Reference slice in slice timing in SPM”

  1. Agree. I also found that using a different reference slice changed my GLM results (event-related analysis).

  2. well, I actually disagree on this issue. First of all, the ref slice selection is a subject preference. If you select the 1st or last slice as ref slice, you also need to do corresponding event timing adjustment when you specify the design matrix in SPM. This latter procedure ensures that the ref slice maps with the event timing well.
    Also, give the slice-timing correction tends to get more inaccurate with longer time delays, when TR = 2s, for example, people tended to choose the middle slice so that the maximum correction is 1 sec not 2 sec if you choose the 1st or last slice as ref slice.

  3. @Jian
    Yes, if you do event timing adjustment then you can choose whatever ref slice. But it seems many people are not aware of this.

    about 2nd point: I thought about this and think there is no “better” ref point (slice). If you choose the middle slice as ref, then you need to adjust the 1st and last slice by 1s; if you choose the 1st slice as ref, you need to adjust the middle slice by 1s, but you don’t need to adjust the last one by 2s, because the last slice is very close to the 1st slice in the next TR so you only need to adjust by TR/n second (n is number of slices). The periodicity of the scanning, I think, makes no slice a better ref slice in terms of accuracy of temporal interpolation.

  4. @xu
    I am glad we can communicate this way! :), again, congrats for the new daddy!
    as to the 2nd point, I still think choose the middle slice as ref slice would save people a lot of hassle since every slices with that specific TR is aligned well, and you only need to remember one lag. Instead, if you choose the 1st slice as ref slice, and decide to align the last slice with the first slice (ref slice) of next TR, then you need to worry about that for this specific slice, you have a different ref point, and I don’t think SPM or any conventional Statistic package allow you to deal with slices that are “acquired” at different time and automatically adjust that.

  5. for the event timing adjustment in the design matrix, do you mean the hrf with time derivative? But, I was suggested not to use this option…

  6. @Ting Xiang
    Some people adjust the onsets of events in their designs, but the correct SPM way for doing this is adjusting microtime resolution and onset. Sorry for posting at such an outdated discussion, but someone may occasionly find it on google just like me 🙂

Leave a Reply

Your email address will not be published. Required fields are marked *