fMRI图像预处理

fMRI图像简介

脑成像主要分为两种：结构性脑成像与功能性脑成像。并且存在一些不同的类别来执行每种类别的成像。
结构性脑成像：它涉及脑结构的研究，也涉及疾病和损伤的诊断。例如，如果你在一次事故中然后得了中风或类似的疾病，你可能会去做一个结构成像看看其影响。
执行方式包括:计算机轴向断层扫描(CAT)，磁共振成像(MRI)，正电子发射断层扫描(PET)。
功能性脑成像：它可用于研究认知和情感过程，是在一段时间内多次测量测量大脑某块区域
执行方式包括：PET ,fMRI ,EEG ,MEG
fMRI与MRI的区别：fMRI需要采集一段时间内的MRI图像然后看这段时间的脑成像变化。

fMRI图像从每个体素中提取到活动强度信息，每个体素都有一个空间的位置，在连续的时间观察这个体素，我们就可以看到这个体素的活动强度随时间的变化，从而提取出来关于时间的活动强度信息，从而关联我们的任务或者被试的活动。fMRI从基础上说是每个体素的时序变化，这些信号代表血氧依赖水平(BOLD 对比度)，测量血红蛋白与脱氧血红蛋白
的比例，测量该组织的新陈代谢速度，间接反映神经元活动。
血流响应函数(hemodynamic response function,HRF)，代表神经活动触发的fmri响应。

fMRI数据的处理过程

⁃ 数据获重建过程
⁃ 预处理过程：时间校准，头动校准，配准与标准化，空间平滑处理
⁃ 数据分析过程：定位大脑活动、连接强度分析、预测

fMRI的数据结构

时间分辨率(temporal resolution)：时间分辨率由扫描一张图片所需的时间决定，事实上扫描一张图片的时间是由一个扫描时的参数TR决定的。他决定了我们区分观察不同时间点上大脑变化的能力。
空间分辨率(spacial resolution)：空间分辨率决定了我们区分不同区域微小变化的能力（图片是否清晰）。
在磁场等条件一样时，fMRI扫描出的图像的时间分辨率越高/低，对应的空间分辨率就越低/高，很多时候我们需要根据不同的情况在两者之间找到一个平衡点。

结构图像(Structural (T1) images)有着极其高的空间分辨率。我们可以通过它区分不同的组织。在这个图里面你可以分辨灰质和白质，也可以找出解剖上的边界等等各种结构上的微小细节。然而结构图像的时间分辨率却极其的低，事实上它就是一个静态的图像。
功能图像(Functional (T2*) images)的时间分辨率相对很高。通过它我们可以研究不同时间段里大脑发生了什么变化，例如在这里有两个情况下扫描得到的“大脑功能图像序列”，情况A是用指头敲击，情况B是静息状态。这样我们就可以分析两种状态相比之下大脑的信号有什么区别。但其空间分辨率较低，相对结构图像来说很模糊。

在一个实验中，通常我们都会选取很多被试（subject）去做相同的任务（task）。每一个被试会参加很多轮扫描（sessions），在一轮里又会被扫描很多次（run）。我们通常会让被试重复完成几次实验任务。每一次扫描都会得到很多数据。每一次有很多体素，每个体素都有一个自己的强度值。这也是最基础的测量值。所以整个实验的数据就是这样一层层下来的，一个实验有很多被试，每个被试会参加很多轮扫描。每次扫描又会扫描很多次。每次扫描又会包含很多体素。
通过改变我们激发原子核的间隔时间TR和激发后到我们开始收集数据的时间TE，就可以控制哪些组织的特征得到强调。

MRI的目标是构建一幅图像或者是一个与空间定位对应的数字矩阵，这幅图像描述了样本的原子核的一些特性的空间分布，这可能是原子核的密度或者它们所属的组织的弛豫时间。
在fMRI中我们还有另一个对比，叫做T2*，它是T2和主磁场相互作用产生的结果。因为T2很重要，所以在这里稍微再讲一下T2。
因为脱氧血红蛋白会抑制MR信号，因此脱氧血红蛋白浓度提高就会导致fMRI信号的减弱，我们把这种信号的减弱称为T2衰减，它是由于脱氧血红蛋白导致主磁场的紊乱而加快的信号衰减。实际上BOLD fMRI测量的就是在富含氧合血红蛋白与富含脱氧血红蛋白的血液之间其T2信号的差异。而扫描器可以设置成消除或者强化这些磁场的紊乱，后者就是BOLD fMRI的基础。
T2和T2的区别:简单来说，因为不可避免的整体磁场的不均匀和氢原子核与周围其他原子核自旋的交互作用，使得氢原子核的移相更快，从而导致了信号的衰减加速，这样的信号衰减就是T2衰减。换句话说，T2衰减就是没有对氢原子核的移相进行任何矫正的一种信号衰减。而T2衰减是比T2衰减更慢的信号衰减，T2是一种要对磁场的不均匀进行矫正才能获得的一种信号。
根据主要对T1、T2或T2敏感度不同，我们可以构建出不同的图像。因为T1和T2会因组织不同而有所差异，它们可以表示出CSF与灰质和白质之间的边界。而T2对血流和氧合作用比较敏感，它可用于对大脑功能的成像。

fMRI信号、噪音与血流

fMRI往往用于研究大脑的具体功能，扫出来的是功能图像，也叫做T2*权重图像。虽然它的空间分辨率比较低，但是时间分辨率很高，可以在很短的时间内扫出一叠功能图像。这样子就可以研究实验操作究竟是如何影响大脑的MRI信号的。
一个fMRI实验相当于包含了一个MRI图像序列。我们可以从中看出血氧信号随着时间的变化。例如这里有一个MRI的图像序列，我们来着重看看大脑中的一个体素。每个体素都有着自己的空间位置和一个对应的值。如果我们把每个时间点这个体素对应的值单独拿出来，我们就得到了一个时间序列。也就是说每个体素都对应了一个自己的时间序列。

我们其实就是想看看这个时间序列是如何随着我们的实验操作而变化的。大多数fMRI实验中我们用到的这个值都是BOLD（Blood Oxygenation Level Dependent）比例指数，BOLD指数指的是有氧血红蛋白（Oxygenated Hemoglobin）的含量和脱氧血红蛋白（Deoxygenated Hemoglobin）含量的比值。

fMRI的噪音来源

BOLD fMRI包含多种来源的噪音，与设备本身和被试本身相关
系统内*电子的热运动
磁场和其梯度的不稳定性
头动及其对磁场的交互影响
生理影响：心跳、呼吸，co2浓度

噪声如何在数据内出现

高频的spikes
图像的伪影和畸变
低频的漂移和周期性的波动

如何减少噪音和伪影的影响

在获取过程中：
我们需要高质量的控制以确保扫描器本身工作很正常
合理安排实验的序列
我们还可以用一些特殊的序列，比如自旋回波(spin-echo)，同时多层采样，z-shimming(减少某个区域的伪影)
尽可能减少头动
在数据分析中：需要检查一下我们得到的数据

异常值/伪影的鉴别和校正

通过预处理来调整头动和漂移
一些数据处理的手段(回归，层级建模)
还可以建模并去除低频与周期性波动的内容

在检视你的数据过程中，你需要着重关注以下几点：

1. 覆盖范围(FOV)与先前的计划是否一致；
2. RF噪音与畸变的图片
3. 瞬态梯度伪影
4. 重影
5. 磁敏感性伪影
6. 任务相关的头动

fMRI的数据预处理

刚采集的原始图像数据会经过一系列的预处理步骤。这些步骤主要是分辨并去除伪影（Artifact），或者是检验一些模型所需的假设是否成立。具体来说，有三个目的：
尽量减少因为数据采集（Data acquisition）和生理学特性（Physiological）导致的误差
我们想检验一下模型的统计假设，并且尽可能的做一些变换让数据符合这些假设
将不同个体的数据的脑区位置标准化（Standardize the location of brain region）以便于进行组间分析（Group analysis）。只有这样组间分析才具有较好的效度（Validity）和灵敏度（Sensitivity）

预处理步骤包括：可视化（Visualization）、去伪影（Artifact removal）、时间配准（Slice time correction）、头动校正（Motion correction）、生理噪音校正（Correction for physiological effect）、结构功能配准（Co-registration）、标准化（Normalization）和时空间滤波（Spatial and temporal filtering）。
一、Visualization and Artifact Removal
可视化的去伪影。在任何情况下这都一定是预处理的第一步。我们会探索性的看看刚采集的原始图像数据（Raw imagine data）有哪些问题。找到了问题才知道要用什么模型去去伪影。比如有些时候数据会出现一些异常陡峭的峰波（Spike）或者是缓慢的偏倚（Drift）。有很多种方法，例如主成分分析法（PCA，Principle components analysis），可以用来探测这些异常的峰波。

二、Slice Time Correction
在扫描一次完整大脑（Brain volume）的周期（TR，Repetition time）内，我们会扫好几片脑片（Slice），由于我们一个时间点只能扫描一张脑片，如果是按照顺序一张一张扫下去的话，每张脑片之间的扫面时间点都会有区别。比如最顶端的脑片的扫描时间会相对低端的脑片延迟2s。例如在这张图里。脑片1、2、3分别在不同的时间点上被扫描。如果我们不做时间回归直接处理的话，由于我们采样的时间点不相同，三个脑片的信号看起来很不一样，尽管实际上他们来源于同一个信号。

因此我们需要用插值（Interpolation）等方法来获得三个脑片相同时间点的数据。所谓插值就是通过与未知点相邻的已知点的信号值来预测位置点的值。我们通常用的插值函数类有线性函数和三角函数。当然，你也可以用傅里叶变换。求出傅里叶级数后平移相位就好。（其实个人认为这个就是三角函数插值）
三、Motion Correction
头动（Head motion）也是一个很严重的干扰。我们在做后续分析的时候需要假设每个体素时时刻刻对应的都是大脑的同一个位置。可是扫描的时候被试的头肯定会有微小的移动。这样的话一个体素可能在下一秒就对应了脑中另外一个截然不同的位置。因此我们会用刚性变换（Rigid body transformation）把所有的图像中的脑都固定在同一个靶位置（Target imagine）（我们假设头动过程中脑一直是个刚体）。
一个刚体变换包含了六个*度，即关于X,Y,Z轴三个方向的平移（Translation）以及围绕这三个轴的旋转（Rotation）。
下面这张图说明了以上的几种基本变换，当然有时候我们还可能用到缩放变换（Scaling）和切变换（Shearing）其实还有很多的线性变换可以用的。例如：
相似变换（平移，旋转和等比缩放，七个*度）
仿射变换（关于三个轴的平移，旋转，缩放和切变换，12个*度）

除此之外还有非线性变换（Warping method），即方程对图像的处理是非线性（Nonlinear）的。
不管怎样我们变换的最终目的是使处理过的图像与靶图像尽量的重合。我们会有一个代价函数来描述处理过的图像与靶图像的相似度和重合度。例如我们可以使用所有体素数据值的差的平方和或者互信息（Mutual information）。
注：信息论里一种有用的信息度量，它可以看成是一个随机变量中包含的关于另一个随机变量的信息量，或者说是一个随机变量由于已知另一个随机变量而减少的不肯定性（百度百科-互信息）

一般来说靶图像选取的是第一张图像或者所以图像关于时间的平均图像（Mean imagine）。例如这张图，我们有一个靶图像，还有一个等着我们变换的图像（Input imagine）。我们希望经过一些列的平移很旋转使蓝色的大脑和绿色的重合在一起。这就是头动校正。

注：软件计算出的旋转角度和平移距离的曲线，横坐标是输入图像（Input imagine）的序号（按时间排列）
四、Co-registration
接下来的事情就是结构图像与功能图像配准（Co-registration）了。我们希望能把功能图像上的点定位在有着较高分辨率的结构图像上。

虽然我们分析的是较低分辨率的功能图像，但当我们将其定位在结构图象上时，我们可以发现一些更细节的问题。配准也为我们之后要讲到的标准脑（Standard coordinate system）提供了很大的帮助。
一般来说我们会在实验的每轮扫描的开始扫一张结构图像。虽然说结构图像和功能图像是在同一论扫描里扫得的。但往往它们往往是没法重合覆盖在一起的。一是因为他们测量不是一个东西（编者注：一个是组织一个是血流信号），二是因为测量过程中会有头动。
结构图像和功能图像有这么几个区别：

信号的数量级不一样，没法直接对比
形状不一样。至少要用仿射变换才能配在一起。在这里我们一般使用Mutual Information 来作为代价函数

注：软件SPM算出的配准变换方式。
五、Normalization

图1为fMRI数据预处理流程
在这个模块中我们会讨论将图像与脑图谱模板（atlas template）配准及平滑（smoothing）的过程。通常来说，每个人的大脑都是不同的，两个个体的大脑尺寸差异可以达到30%，大脑形状也可能存在显著差别，所以标准化（normalization）所做的正是通过拉伸、压缩以及卷绕使得扫描得到的大脑与标准大脑模板一致。这对于组分析（group analysis）而言很重要，因为我们想要对比不同的大脑，当我们关注某一个体素（voxel）时，我们希望能够在所有被试个体的扫描图像中看到这个体素，这就需要所有的大脑“匹配”至同一个标准脑空间（standard brain space）。所以高分辨率的T1结构像现在又被配准至标准空间（normalized space），如下图2所示。

图2将高分辨率T1结构像配准至标准空间
注：从图 1的流程图中可以看出，课程中所讲的标准化是先将高分辨率的T1结构像与标准模板（此处的标准模板是高分辨率的结构像标准模板）进行配准，再直接将这些配准参数应用于功能像上（因为功能像之前已经与结构像进行过一次配准），这也是目前比较普遍的一种做法，因为T1结构像分辨率高，大脑内部结构可精准对齐。但还有另外一种做法，可以直接将功能像与标准模板（此处的标准模板是低分辨率的功能像标准模板）进行配准，这种做法在SPM软件中比较常见，不需要借助T1结构像，但在对齐效果上较差。
标准化的优点是：
可以以一致的方式报告和解释空间位置。在标准空间内讨论某一坐标时，我们总能知道它在哪里
实验结果可以被推广到更大的群体，因为我们可以对不同个体进行比较
实验结果可以跨研究进行比较，前提是每个研究在标准化过程中采用同一大脑模板
不同个体的实验结果可以进行平均
缺点是：
会降低空间分辨率
引入潜在误差

目前主要有两个脑图谱模板被广泛使用。
Talairach图谱，由Talairach建立，随后Talairach和Tournoux在1988年对其进行了更新。Talairach图谱是基于一位60岁女性尸体的单个大脑半球建立的，在该图谱中，原点被设定在前联合（Anterior Commissure），水平向被设定为前联合与后联合（Posterior Commissure）的连线方向。
MNI脑图谱（Montreal Neurological Institute Brain Atlas）也越来越常用。该图谱集合了许多正常对照组的MRI扫描数据，在现行的ICBM-152模板中就集合了152个右利手个体的脑区图。相较于基于单个个体建立的脑图谱（如Talairach图谱），MNI图谱更具代表性。
进行标准化的方法有很多，这里作简要介绍：
标定点法（Landmark-based methods）：根据大脑的解剖结构特征，找到不同个体的大脑标定点，并将其与标准模板的标定点对齐（如Talairach图谱的标记点就包括前联合、后联合、颅正中矢状平面以及在每一方向上的大脑外边界）。
基于体积配准法（volume-based registration）：利用线性和非线性变换实现（与结构像和功能像配准时所采用的方法相似）。
计算解剖学法（computational anatomy methods），采用微分同胚变换。
表层法（surface-based methods），主适用于皮层表面（you sort of blow up the brain and do normalization that way）。

六、Spatial and temporal filtering空间滤波smoothing
在统计分析前对采集到的数据进行空间平滑（smoothing）是很普遍的，这可以提高信噪比（signal to noise ratio），消除伪影（artifacts），使数据更符合分布假设（distributional assumptions）。例如，我们经常假设数据呈正态分布，通过对许多体素取平均值的平滑过程，可以增加其呈正态分布的可能性。
空间滤波的优点是：
通过模糊残留的解剖结构差异，克服标准化过程的缺陷。标准化过程并不完美，有时会造成图像上出现一些解剖结构差异（失匹配现象），这就需要通过模糊的方式来消除差异。
提高信噪比。例如，有一片区域被**，我们对这片区域进行均值滤波，可以保留其信号强度但降低噪声强度（因为噪声多为正态分布，有正有负，取均值后噪声强度降低，均值滤波实质是一种低通滤波，保留低频的信号成分，滤掉高频噪声成分），从而提高信噪比。
增强统计分析的效力。
满足高斯随机场理论的分析要求。高斯随机场理论在多重比较中广泛应用。
缺点是：
降低图像分辨率，丢失信息。

对图像进行空间滤波时，我们经常使用高斯核（Gaussian kernel），该平滑核的尺寸（平滑范围）由半高宽（full width at half maximum, FWHM）决定。如下图 3所示，半高宽是指：分布函数在其峰值一半位置上的峰宽。对于高斯分布而言，其半高宽与标准差存在关系：

图3半高宽
根据匹配滤波器（matched filter）的设计原则，与信号尺度相匹配的滤波器可以为我们提供最大信噪比，所以如果我们知道**区的确切尺度，就应该以该尺度进行平滑，这样可以保留信号幅度并降低噪声。但通常平滑核的尺寸是在观察**区之前就根据先验知识确定的，所以我们不大可能使得该尺寸与**区尺度完全匹配。另外，我们通常采用同样尺寸的平滑核对整个大脑图像进行平滑，但在不同脑区，**区的尺度不尽相同。为了避免这样的问题，可以采用自适应平滑（adaptive smoothing）的方法，涉及非平稳空间的Gaussian-Markov随机场，实现空间和时间的平滑。
注：通常来说，平滑核的选取在不同的实验要求下有不同的标准，增加平滑核尺寸可导致检测到更大范围的**区域，但会减弱对较小**区域的检测能力，如下图 4所示。一般选取两倍于体素尺寸的平滑核，即若体素大小为3×3×3，一般选取FWHM=6mm。

图 4 平滑核尺寸对**检测造成影响的例子

最后，再来回顾一下fMRI数据预处理的流程（如下图 5所示），扫描得到的结构像和功能像都需要进行预处理。结构像，通常是T1加权，需要与功能像进行配准，并配准至标准空间。功能像，通常是T2*加权，需要进行去噪，扫描层时间校正，以及头动校正之后，将标准化过程中结构像的配准参数应用于功能像上，使功能像也配准至标准空间（如MNI脑图谱）。最后，功能像进行平滑，就可以开始数据分析了。

图 5 fMRI数据预处理流

转载申明：

链接：https://zhuanlan.zhihu.com/p/22385985
来源：知乎
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。