版块推荐
    IE收藏打印复制链接
    «12»Pages: 1/2     Go

    心心水滴论坛SPM5中文教程新编

    -

    选择 0滴只看楼主 wonder  于2008-01-16 21:51     
    总版主 224心级 4197水晶 0钻石30精华; 专业:connectivity/language
    遥想52brain的未来
    -写想法

    在老论坛编译过一个所谓渔版spm5教程,得到了很多滴友的鼓励和支持。可惜种种原因没能完全坚持下来写完。此次响应当前论坛“发好贴,吸引更多滴友参与”的主题管理号召,决心再把这个非常有意义的工作做起来。

    一个人的精力有限,欢迎有能力有时间的滴友来参与一起创作,大家共同汇编成一本。我想编译也不拘泥于什么固定形式,您可以找一章或一节自己感兴趣的来编译,按照自己的理解和喜好来整理叙述。然后别忘了加上您的名字。比如心心水滴论坛...编译。

    今天时间有限,先把老的转过来吧。基本内容还是参照官方的SPM5 manual(http://imaging.mrc-cbu.cam.ac.uk/pdfs/SPM5manual.pdf)来编译。 我想这个帖子大家除了回复自己的编译章节之外就不要另外做其它回复来打断了,想要支持就给辛苦的编译者们送鲜花吧。
    赞好帖- 送花 || 扔蛋-抨滥贴  
    本楼已得:4分。最近评分:
  • 水晶:+1(tearless) 好帖,送鲜花!
  • 水晶:+1(echo2010) 有价值+10
  • 水晶:+1(seashell) 好帖,送鲜花!
  • 搜索更多相关文章:spm 教程 spm5 预处理
    惠皆施来者,来者皆施惠。
    诚信与朋友是人生最大的资产。

    选择 1滴只看该作者 wonder  于2008-01-16 21:54     
    总版主 224心级 4197水晶 0钻石30精华; 专业:connectivity/language
    遥想52brain的未来
    -写想法

    本帖内容由心心水滴论坛(http://52brain.com) wonder 编译自 Functional Imaging Laboratory,Wellcome Department of Imaging Neuroscience,Institute of Neurology, UCL http://www.fil.ion.ucl.ac.uk/spm/ 的英文版本。原理部分为意译,加入了一些自己的理解,力图更直白。但也有一些还不太理解的地方,欢迎大家纠正补充。
    心心水滴论坛版权所有,欢迎大家传阅共享,但是别忘注明来自 http://52brain.com。编译者能力和精力有限,欢迎大家批评指正。


    第一部分:时间维度预处理
    Chapter 1 Slice Timing

    目的:纠正slice(下译为扫描层)之间采集时间的差异。纠正后的新文件名为:a+原文件名。

    原理概述:Sliceorder(扫描层序数) 用于指明扫描层被采集时的顺序,其总数为n,n即一个volume(volume为一个TR时间内的扫描量)里面所包含的扫描层的总数。每个Sliceorder(扫描层系数)同时指明了每个扫描层在图像文件中的存储位置。因此扫描层系数表明了扫描层被采集时的时间顺序。 我们可以用SPM 的 Display 功能来查看这个顺序。移动鼠标的十字准线到 Z 坐标为1的位置时,我们看到的就是一个volume里面第一个slice(扫描层)上的象元(voxel)。

    核磁仪采集数据的时间精度为一个TR,因此我们通常默认一个TR内采集到的数据(一个volume)采自完全相同的时间。但事实并非如此,一个TR内我们需要采集很多扫描层(例如20层)的数据以覆盖较大的脑区(通常是全脑)。这些扫描层是一个接一个按照上面所说的sliceorder(扫描层序数)的顺序来采集的,因此各个扫描层之间的采集时间是有差异的。Slice timing这一步所要做的就是通过一定的算法纠正这一时间差异。

    常规的纠正方法是在保持整段采集信号不变的条件下推前或者推后采集的起始时间。这可以通过简单的移动采集信号的正弦相位来做到。一种常用的数学转换方法:傅立叶转换可以把任何复杂信号转换为由不同频率和相位的简单正弦曲线的线性组合。这样,只要通过加入一个常数到每个频率的相位中,就可以达到移动数据起始和结束时间的目的。

    具体步骤:

    1.1 Data
    预备数个被试或者session的数据。以下所述参数设置将被应用到所有所选数据。

    1.1.1 Sessions
    选择需要纠正的所有图像文件。

    1.2 Number of Slices
    输入扫描层总数。

    1.3 TR
    输入 TR,以秒为单位。

    1.4 TA
    TA必须由使用者手动输入,单位为秒. 一般可以用以下公式计算:TA=TR-(TR/扫描层数). 这里你可以不必计算出结果,直接写上带入了数字的公式就可以了。比如,假如TR为3秒,扫描了20层,则可以直接写为:3-(3/20)。

    1.5 Slice order
    输入扫描顺序。如前述此顺序可通过SPM的Display功能查看。
    ascending (扫描序数从底部到顶部排列,即从1顺序递增到n)
    descending (扫描序数从顶部到底部排列,即从n顺序递减到1)
    interleaved(扫描序数间隔递增或递减,一般顺序为2,4,6,8,10..1,3,5,7,9..n)

    1.6 Reference Slice
    选择参考扫描层(一般可使用默认值),其它扫描层的起始时间都将以此层的起始时间为标准来移动进行校正。
    (心心水滴论坛:http://52brain.com
    赞好帖- 送花 || 扔蛋-抨滥贴  
    惠皆施来者,来者皆施惠。
    诚信与朋友是人生最大的资产。

    选择 2滴只看该作者 wonder  于2008-01-16 21:56     
    总版主 224心级 4197水晶 0钻石30精华; 专业:connectivity/language
    遥想52brain的未来
    -写想法

    第二部分:空间纬维度预处理
              (心心水滴论坛(http://52brain.com)wonder编译制作)
              
              
              Chapter 2  Realign
              被试内的图像数据重排。
              
              2.1 Realign: Estimate
              (重排参数的估计)
              此步骤采用最小方差原理和六参数刚体空间变换来重排从同一个被试上采集到的图像数据[19]。使用者所选取的第一幅图像文件将被作为其它图像重排的参考标准。也就是说,你想要用哪一幅图像作为参考标准,就先选哪一幅图像的文件。参考图像文件不一定非用采集到的第一幅图像,使用最有“代表性”的一幅图像也许更好。
              
              本步骤的目的主要是去除fMRI和PET数据中的运动伪影。图像数据的头文件会被改写以反映数据相对空间位置的变化。此过程的具体参数会在结果窗口中以平移(translation)和旋转(rotation)曲线图显示。每个session的重排参数会被存储到名为 rp*.txt 的文件中。这些参数可以在最后的一般线性模型统计估计中作为混淆因素考虑进去[19]。
              
              2.1.1 Data
              选择一个被试的需要进行此步骤处理的所有sessions。
              
              注:在 coregistration 这一步,首先是对所有的session进行重排,其具体做法是把所选每个session的第一个scan与所选第一个session的第一个scan对齐。然后再把每个session里的其它scan与该session的第一个scan进行对齐。使用此方式进行重排是因为各个session的数据之间可能会有较大差异。
              
              Session
              选择session里所有的scan。
              
              2.1.2 Estimation Options
              这里包括各种注册参数选择项,若对某一个选项不确定,使用软件默认值即可。
              
              Quality
              质量与速度的权衡。选择高质量以最慢的速度给出最精确的结果,低质量以较快的速度给出较不精确的结果。
              此参数的设定实际影响到的是参与参数估计的象元(voxel)的数目。其依据是有些象元(voxel)其实对重排参数的估计贡献不大,可以舍弃。
              
              Separation
              此参数以毫米为单位,表示对参考图像文件进行重采样时采样点之间的间隔。采用点之间间隔越小,结果越精确,运算速度越慢。
              
              Smoothing (FWHM)
              高斯平滑的半高宽值。在估计重排参数之前一般先进行高斯平滑。
              PET数据一般使用 7 mm。
              MRI数据一般使用 5 mm。
              
              Num Passes
              Register to first: 所有图像文件对齐注册到第一幅图像。I
              Register to mean: 使用two-pass处理将所有图像文件对齐注册到所有图像文件的平均图像。
              
              PET 数据一般注册到平均图像。因为PET数据相比fMRI数据噪音更大,文件更少,所以时间的影响更小。
              MRI数据一般注册到第一幅图像。虽然使用two-pass处理可能更精确,但是其对效果的提高与其所损失的运行时间相比得不偿失。
              
              Interpolation
              在估计最佳变换时对数据进行重采样的方法。
              高的degree提供更好的结果,但是也更慢,因为会采样更多的相邻象元(voxel)[52, 53, 54]。
              
              Wrapping
              此参数指示一个volume中数据wrap around in的方向(此处具体理解有待大家补充)。
              
              No wrapping:适用于PET数据或者已进行过空间变换的数据。同时当你不确定自己数据类型时,推荐使用此选项。
              Wrap in Y:适用于没有重排(resilce)过的在Y方向上进行相位编码的MRI数据。
              
              Weighting
              提供一个加权图像,在估计重排参数时对参考图像的每一个象元进行加权。加权系数与标准差成反比。例如当有大量额外的头动(如说话或者特定区域内的严重伪影)时。(此处具体理解有待大家补充)。
              
              (心心水滴论坛(http://52brain.com)编译制作)
              
              2.2 Realign: Reslice
              (据已估计出的参数重排)
              此功能重排上步骤中已进行参数估计和注册的图像文件,使之与参考图像文件达到象元级的匹配精确。重排后的数据被命名为:r+原文件名。
              
              2.2.1 Images
              选择要重排的数据文件。
              
              2.2.2 Reslice Options
              各种重排参数设定,若对某一个选项不确定,使用软件默认值即可。
              
              Resliced images
              All Images (1..n) : 重排所有数据,包括标准参考图像(重排后还是保持原位置不变)。
              Images 2..n : 重排除了标准参考图像之外的所有数据。此选项用于当你以MRI结构像为标准
              重排PET图像数据,而又不想在结果中再生成一个等同的MRI标准结构像时。
              All Images + Mean Image : 重排图像文件之外,另生成一个重排后的平均图像文件。
              Mean Image Only : 只生成重排后的平均图像文件。
              
              Interpolation
              图像文件重采样和重写入的方式。
              .Nearest Neighbour :最快,但不推荐使用。
              Bilinear Interpolation:可用于PET数据,但不是太适用于fMRI数据。
              Fourier Interpolation:此选项仅适用于纯刚体变换,也就是说象元大小必须是相同,并且等方性(正方体)的[17, 14]。
              
              Wrapping
              此参数指示一个volume中数据wrap around in的方向(此处具体理解有待大家补充)。
              
              No wrapping:适用于PET数据或者已进行过空间变换的数据。同时当你不确定自己数据类型时,推荐使用此选项。
              Wrap in Y:适用于没有重排(resilce)过的在Y方向上进行相位编码的MRI数据。
              
              Masking
              因为扫描过程中被试总会或多或少有头动,造成同一个时间系列数据里所采集到的图像的边界不会完全重合。在有些图像还有数据的地方(信号值大于0),其它一些图像已经超出了图像边界(信号值为0)了。在这些信号为0的区域是无法采样数据的,因此SPM只要检测到某一幅图像在某个区域已经超出了边界(即信号为0),就会将其它所有图像的此区域信号值均设为0。此做法相当于取了时间系列数据中所有图像的交集。
              
              2.3 Realign: Estimate %26 Reslice
              将上述参数估计与数据重排合到一起做。全部选项与参数原理均与2.1和2.2中对应项相同。  
              (心心水滴论坛(http://52brain.com)编译制作)
    赞好帖- 送花 || 扔蛋-抨滥贴  
    惠皆施来者,来者皆施惠。
    诚信与朋友是人生最大的资产。

    选择 3滴只看该作者 yubing 于2008-01-17 09:10     
    正式滴友 8心级 29水晶 0钻石2精华; 专业:脑科学
    最近我脑内没啥想法
    -写想法

    谢谢wonder
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 4滴只看该作者 fmri 于2008-01-17 12:25     
    荣誉总版 136心级 292水晶 0钻石10精华; 专业:心理学
    小人物
    -写想法

    :cacakiki21:
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 5滴只看该作者 jia_u 于2008-04-08 20:01     
    待验证滴友 0心级 -12水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    yxh98:
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 6滴只看该作者 yishuihan 于2008-06-21 08:33     
    待验证滴友 0心级 4水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    已经下载,很有用的资料.向你学习.......LZ好好人的.....
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 7滴只看该作者 zshtom 于2008-06-25 16:50     
    正式滴友 8心级 1750水晶 0钻石0精华; 专业:脑科学
    最近我脑内没啥想法
    -写想法

    很好 谢谢!!
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 8滴只看该作者 hsxu 于2008-06-27 17:20     
    新新滴友 2心级 22水晶 0钻石0精华; 专业:生物物理-认知神经心理学
    最近我脑内没啥想法
    -写想法

    非常感谢,看英文真是累啊~还是中文好。
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 9滴只看该作者 sallyZL 于2008-07-15 13:19     
    待验证滴友 1心级 62水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    谢谢楼主

    下了英文的pdf文件,并结合楼主的中文翻译进行了学习,觉得这里的资料对我这样的入门者来说实在是太需要了,再次感谢!希望在不久的将来,我也能在这里参与大家的讨论!
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 10滴只看该作者 multiwavelet 于2008-07-31 15:11     
    新滴友 5心级 5水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    支持一下
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 11滴只看该作者 hbcl5960 于2008-08-10 01:30     
    新滴友 3心级 23水晶 0钻石0精华; 专业:脑科学
    最近我脑内没啥想法
    -写想法

    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 12滴只看该作者 chuan0917 于2008-09-10 08:54     
    新新滴友 2心级 50水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    确实是好东西哦 我英文看得一头雾水 现在雾水少一点了
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 13滴只看该作者 elangbean 于2008-09-11 03:13     
    正式滴友 11心级 238水晶 0钻石1精华; 专业:
    最近我脑内没啥想法
    -写想法

    谢谢wonder,学习ing,希望以后也能来贡献点东东
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 14滴只看该作者 zhuliuhong 于2008-09-12 09:33     
    待验证滴友 1心级 70水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    顶顶顶!!!好东西,就要顶,谢谢楼主。Thanks
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 15滴只看该作者 Jana 于2008-10-07 20:08     
    新滴友 3心级 4水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 16滴只看该作者 simonwy_cas  于2008-10-08 10:10     
    核心滴友 60心级 11水晶 0钻石4精华; 专业:脑科学
    booby
    -写想法

    太强大了
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 17滴只看该作者 simonwy_cas  于2008-10-08 10:11     
    核心滴友 60心级 11水晶 0钻石4精华; 专业:脑科学
    booby
    -写想法

    引用
    原帖由 wonder 于 2008-1-16 21:51 发表
    在老论坛编译过一个所谓渔版spm5教程,得到了很多滴友的鼓励和支持。可惜种种原因没能完全坚持下来写完。此次响应当前论坛“发好贴,吸引更多滴友参与”的主题管理号召,决心再把这个非常有意义的工作做起来。

    一个人的精 ...

    这种方式好啊!目前流行的分工合作写书模式。。。 可惜偶不行,只能大声喊加油咯!
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 18滴只看该作者 yanski  于2008-10-12 14:11     
    正式滴友 34心级 2794水晶 0钻石1精华; 专业:影像医学(fMRI)
    退一步海阔天空
    -写想法

    向楼主学习
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 19滴只看该作者 ivan0909 于2008-10-14 09:12     
    待验证滴友 1心级 36水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    好贴
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 20滴只看该作者 wo657832009 于2008-10-20 18:14     
    待验证滴友 0心级 4水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    好人啊 顶顶顶
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 21滴只看该作者 ivan0909 于2008-10-23 10:59     
    待验证滴友 1心级 36水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    望后续
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 22滴只看该作者 sunshy 于2008-10-23 22:14     
    新滴友 3心级 107水晶 0钻石0精华; 专业:
    最近我脑内没啥想法
    -写想法

    太感谢了,真是入门宝典
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 23滴只看该作者 巧巧宝宝 于2008-10-26 11:28     
    正式滴友 23心级 440水晶 0钻石1精华; 专业:生物医学工程
    最近我脑内没啥想法
    -写想法

    想问:在最后的一般线性模型统计估计中载入 rp*.txt ,与不载入rp*.txt 对结果有什么重大影响吗?谁能说说看?
    引用:
    本步骤的目的主要是去除fMRI和PET数据中的运动伪影。图像数据的头文件会被改写以反映数据相对空间位置的变化。此过程的具体参数会在结果窗口中以平移(translation)和旋转(rotation)曲线图显示。每个session的重排参数会被存储到名为 rp*.txt 的文件中。这些参数可以在最后的一般线性模型统计估计中作为混淆因素考虑进去[19]。
    赞好帖- 送花 || 扔蛋-抨滥贴  

    选择 24滴只看该作者 Carlton 于2008-10-30 22:52     
    正式滴友 20心级 289水晶 0钻石1精华; 专业:
    最近我脑内没啥想法
    -写想法

    回复 24滴 巧巧寶寶 的帖子

    載入 rp*.txt 的目的主要是要移除 "移動因素"造成的信號變化, (remove the variance in the signal cased by the movements)
    Or more precisely, those movment paramters can tell if you have correlated movement with you design.
    或者是說 如果subject有跟實驗刺激有關聯性的移動(例如受到刺激的時候頭就抬一下),你在SPM specify design的時候就可以發現是不是實驗刺激跟movement parameter 有很高的correlation
    這時候 你可能要可慮丟掉這一組數據
    因為你的分析結果產生的Activaton cluster很有可能是移動造成的 而不是真正的BOLD signal
    或者, 你可以把movement parameters 當成confounding variablel 而移除掉他們產生的variance
    但是如果你的實驗刺激跟movement的correlation 太高,你在做GLM parameters estimation的時候 你實驗刺激explaining power會大減

    Anyway,總而言之 如果你發現你的movement parameters ( rp*.txt) 跟design matrix 有很高的correlation, 那就可慮換個subject吧
    赞好帖- 送花 || 扔蛋-抨滥贴