专利名称：一种联合fmri和meg的脑区因果连接检测方法
技术领域：
本发明属于图像处理领域，具体涉及一种利用脑功能成像技术的脑区因果连接检测方法，尤其涉及一种联合功能磁共振成像(FMRI)和脑磁图描记术(magneto encephalography, MEG)进行脑区因果连接强度和方向的检测的方法。
背景技术：
随着科学技术的进步，无创性脑功能成像技术日新月异，为我们用多种成像技术进行脑的认知和研究提供了广阔的舞台。单一成像模态，由于成像机理决定了每一种成像技术都具有相应的优点和不足，因此没有一种单一的适合所有实验研究或临床应用的成像方法。目前，功能磁共振成像(Functional Magnetic Resonance Imaging,FMRI)技术以其高空间分辨率、非侵入式等特点在神经疾病诊断治疗方面得到了广泛应用。FMRI技术是通过测量由神经活动引起的脑血流和脑血氧等成分变化而造成的磁共振信号变化来反应人脑内部活动的。由于在代谢水平上测量磁共振信号的变化相较神经元活动具有一定的时间延迟，因此，FMRI技术仅能达到秒级时间分辨率精度。近年来，脑磁图描记术(Magneto enc印halography，MEG)，又称脑磁图仪，逐渐在神经影像领域得到应用。MEG技术可以实时测量脑内神经元之间信号传递时突触后电位变化所产生的磁场信号，因此能够准确测量毫秒级神经元电活动，具有非常强的时间定位能力。但是由于该技术难以准确定位颅内信号源，其空间分辨率较低。因此，如果将FMRI技术与MEG技术的优势结合起来，我们就能达到同时提高空间分辨率和时间分辨率的目的。脑功能一般遵循两个基本组织原则功能集成化和功能特异化。在大空间刻度上， 一个复杂的脑功能可能会由许多功能特异的脑区通过相互作用(集成)来完成；同时，某一特异性脑区也会对许多不同的刺激任务进行表示或加工，通过精细空间刻度上不同的分布式脑活动来对外部不同刺激进行表示。人脑是一个复杂的系统，在不同状态时期各个脑区之间的相互作用会发生变化，反映在脑区间的功能连接强度和方向上。格兰杰因果分析 (Granger causality analysis,GCA)可以用于研究区域间的因果连接关系，为研究揭示人脑信息加工机制提供支持。
发明内容
本发明的目的在于提供一种同时具备高空间分辨率以及高时间分辨率优势的脑区因果连接强度及方向的检测方法。为实现上述目的，本发明提出一种联合FMRI和MEG的脑区因果连接检测方法，包括以下步骤步骤Μ，对扫描脑区获取的FMRI图像进行预处理，并从预处理后的FMRI图像中提取激活区坐标；步骤Sb，对扫描脑区获取的MEG数据进行预处理，并基于提取的FMRI激活区坐标，提取与该激活区坐标对应的脑区的脑区时间序列；步骤&，依据提取的MEG脑区时间序列，计算脑区间的因果连接强度和方向，利用有向网络图显示显著的连接。本发明所述的联合FMRI和MEG的脑区因果连接检测方法，能够有效地利用FMRI 高空间分辨率和MEG高时间分辨率的优势。真实试验数据结果表明，本发明比传统的仅利用FMRI数据更完整准确地检测脑区因果连接，为脑功能数据分析和脑区因果连接检测方法提供了一种新思路，在基于多模态医学成像技术的脑认知加工研究中具有一定的应用价值。
图1为本发明提供的联合FMRI和MEG的脑区因果连接检测的方法流程图。图2是实例中使用本方法检测的脑区因果连接强度图。图3是实例中仅使用单一模态检测的脑区因果连接强度图。
具体实施例方式为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。本发明主要考虑利用功能磁共振成像(FMRI)技术高空间分辨率的优势，定位大脑FMRI图像上的关键神经元位置信息，继而在MEG技术获得上述神经元的高分辨率时间序列，从而在尽可能不丢失大脑信号活动时空二维信息的前提下，利用格兰杰因果分析方法进行脑区因果连接强度及方向的精确检测。图1是本发明所涉及的一种脑区因果连接强度及方向的检测方法的一个具体实施例，其具体涉及联合利用功能磁共振成像(FMRI)和脑磁图描记术(MEG)的脑区因果连接强度及方向的检测方法。该实施例的具体实施步骤如下步骤Μ，对扫描脑区获取的FMRI图像进行预处理，并从预处理后的FMRI图像中提取激活区坐标；1.对FMRI图像进行预处理由于在磁共振成像的扫描过程中，被试的头动、图像内层间扫描时间的差别以及外加磁场的不均勻性等各种因素均会影响采集到的FMRI图像的质量，因此我们需要在保留脑功能图像细节的同时，通过使用脑功能磁共振图像与标准模板进行仿射配准变换方式的预处理，提高脑功能图像的信噪比。所谓仿射配准变换方式是指进行对原始图像缩放，旋转，平移后，根据原始图像和标准模板间的共有特征进行几何配准。所述预处理步骤包括切片扫描时间对齐、图像序列对齐、配准、标准化、平滑、滤波等。在本实施例中，对FMRI图像依次进行图像序列对齐、配准、标准化、空间平滑滤波。所属技术领域的技术人员也可以根据实际需要采用其他预处理方式，或者上述预处理方式的任意组合。2.提取FMRI图像的激活区坐标一般来说，提取激活区坐标的方法有基于一般线性模型(GeneralLinear Model, GLM)的方法，模式分类的方法和基于独立成分分析(Ind印endent Component Analysis,ICA)的方法。其中基于ICA的脑功能激活成分提取模型如下设A e Rnxm为信号混合矩阵；X为脑功能信号的N维观测向量，C为由M(N>M)维统计独立的信号分量构成的信号向量，其中包括与激活相关的脑功能信号分量。ICA分析法可以在混合矩阵A未知的情况下， 寻找线性映射w，从观测的脑功能信号中提取不能被直接观测的原始信号c ( y) :y = wTx
=wAco本发明采用基于数据驱动的ICA方法，首先使用基于最小描述长度(Minimum Description Length, MDL)准则估计独立成分数，继而采用基于互信息最大化的算法计算解混矩阵，分离源信号，获得包含所需激活脑区的独立成分及其Talairach坐标 (Talairach坐标是国际公认的一个标准化的脑图谱坐标系)。步骤Sb，对扫描脑区获取的MEG数据进行预处理，并且基于步骤&提取的FMRI图像的激活区坐标，提取与该激活区坐标对应的脑区的时间维度信息(即脑区时间序列)；1.MEG数据进行预处理由于MEG数据在采集过程中容易受到环境噪声的影响，我们有必要对其进行预处理，以提高数据的信噪比；此外，由于MEG技术的高时间分辨率特性，采集的数据包含了大量的时间维度信息。因此，为保证计算精度与运算时间的平衡，预处理步骤包括数据带通滤波，合适的重复下采样频率等。在本实施例中，对MEG数据依次进行0.5Hz-48Hz的带通滤波和300Hz的下采样频率。所属技术领域的技术人员也可以根据实际需要采用其他预处理方式，或者上述预处理方式的任意组合。2.提取与激活区坐标对应的脑区的时间维度信息如何从获得的MEG数据得到某个局部脑区的时间序列信号，就其本质来说，其实是一个脑磁逆问题，即利用扫描的磁场数据来推断脑内电流源的分布。已有研究证明了导体外的磁场数据无法唯一确定导体内的电流分布，因此脑磁逆问题是不适定的，我们必须在满足条件的解集中通过施加一定限制条件找出合理的解。目前，脑磁逆问题的求解方法主要分为两大类一类基于图像重建技术，其基本思想是将电流分布区域离散成网格，对网络节点处磁源进行重建，重建的磁源参数主要是电流的强度。由于磁场方程中位置参数是非线性的，电流强度参数是线性的，重建时人为地限定了源的位置，因此磁场方程退化为线性方程。另一类为偶极子定位法，该类方法不限定磁源的位置，通常的求解过程是先建立一个目标函数，然后通过调整源的参数使目标函数达到极小，即通过执行一个非线性优化过程来获得磁源参数的解。本发明的实施例采用基于图像重建技术的最小范数解(MinimumNorm Solution, MNS)求解方法。假设磁场方程组为D = LS，其中，D为测得的磁场信号值，S代表节点处电流，L为已知的增益矩阵，取决于网格节点及磁场探测点的位置与头的电导率分布。丽S方法的限制条件是找到满足方程组解中范数最小的S*,即S* = min{ S }，其中，| | ·| 一般取 Frobenius 2 范数，即 | S | = (STS)"2，其中τ表示矩阵的转置。具体求解S*时可以采用基于M00re-Penr0se广义逆L+的方法，即S* = L+D对于广义逆L+的求解，本发明采用基于奇异值分解(Singular ValueDecomposition, SVD)的方法。设L为mXn阶实矩阵(m < η)，其中，m为MEG扫描通道数，η为网格节点数。则对L进行奇异值分解可得
L = U Σ Vt其中U= (U1,u2, -,um)为 mXm 阶正交矩阵，V= (νι;ν2, -,vn)为 nXn 阶正交矩阵，Σ为mXn阶对角矩阵，即Σ= diagOi，δ 2，…，δ ρ)，其中P = m，且S1S δ2彡… 彡δ ρ彡0，若矩阵L的秩rank (L) = r，则上述序列中只有前r个值大于零，S卩δ2彡… ^ δΓ> 5r+1 = ...= δρ = 0那么奇异值分解可表述为
rL = UZVt =YjS1U1V^，其中，Ui(i = 1,2, ...，m)称为奇异值分解的左矢量，Vi (i =
/ 二 1
1,2,…，η)称为奇异值分解的右矢量。由此可求得Moore-Penrose广义逆L+为
r ιL+=Y4-V1U71
i=\由此可根据S* = L+D求出MEG数据所隐含的源空间信息。从而根据ICA分析方法获得的局部激活值最大点体素的坐标，获得该局部脑区的时间序列信息。步骤&，依据提取的MEG脑区时间序列，检测脑区间的因果连接强度和方向，利用有向网络图显示显著的连接。目前，检测脑功能因果连接强度和方向的方法主要有两大类采用动态因果模型(Dynamic Causal Modeling, DCM)和格兰杰因果分析(Granger Causality Analysis, GCA)。两者不同的是，DCM需要预先选定相互作用的区域，并假设这些区域的任意两个之间存在影响，这种预先假设的模型在验证一些有关大脑系统之间的假设时是起作用的，但如果对模型的指定出现偏差，就会导致错误的结论。而GCA则无需事先假设两点之间存在解剖结构的连接性，能更加完整和鲁棒地反应大脑之间的相互连接性。本发明的实施例利用格兰杰因果分析方法检测研究脑区间的连接强度和方向。按照格兰杰因果的定义，若欲判断X是否引起Y，则考察Y的当前值在多大程度上可以由Y的过去值解释，然后考察加入X的滞后值是否能改善解释程度，如果X的滞后值有助于改善对 Y的解释程度，则认为X是Y的格兰杰原因，Y是X的格兰杰结果。格兰杰因果关系检验假设了关于X和Y每一变量的预测的信息全部包含在这些变量的时间序列中。其步骤如下1.将当前的y对其所有的滞后项(yt_1; yt_2，…，yt_q)做有约束回归，即
qyt =Yj^y t-r+ut
/ 二 1其中，y为一个激活区的时间序列，Ut为零均值白噪声，q为模型阶数，通过Akaike 信息准则(Akaike Information Criterion，AIC)计算得到，AIC (q) = 21og(det( Σ ))+2qm2/NΣ = R(O)+ YjA(I)R(I)
/ 二 1其中X(t) = (Xl(t)，X2(t)，…，％(t))T，Xi(t) (i = 1，…，N)为第 i 个激活区的时间序列，R(n) = X(t)XT(t+n)是激活区时间序列矩阵X(t)步长为η的协方差矩阵，Σ为噪声协方差阵，m为激活区数目，N为时间点数目。由此可得有约束残差平方和
权利要求
1.一种联合FMRI和MEG的脑区因果连接检测方法，其特征在于，包括以下步骤 步骤Μ，对扫描脑区获取的FMRI图像进行预处理，并从预处理后的FMRI图像中提取激活区坐标；步骤Sb，对扫描脑区获取的MEG数据进行预处理，并基于提取的FMRI激活区坐标，提取与该激活区坐标对应的脑区的脑区时间序列；步骤&，依据提取的MEG脑区时间序列，计算脑区间的因果连接强度和方向，利用有向网络图显示显著的连接。
2.如权利要求1所述的方法，其特征在于，步骤M中所述对扫描脑区获取的FMRI图像进行预处理，至少包括对扫描脑区获取的FMRI图像进行图像序列对齐、配准、标准化和空间平滑滤波。
3.如权利要求1所述的方法，其特征在于，步骤M中所述从预处理后的FMRI图像中提取激活区坐标，是采用基于数据驱动的ICA方法实现的。
4.如权利要求3所述的方法，其特征在于，所述采用基于数据驱动的ICA方法从预处理后的FMRI图像中提取激活区坐标，包括首先使用基于最小描述长度准则估计独立成分数，继而采用基于互信息最大化的算法计算解混矩阵，分离源信号，获得包含所需激活脑区的独立成分及其Talairach坐标。
5.如权利要求1所述的方法，其特征在于，步骤Sb中所述对扫描脑区获取的MEG数据进行预处理，至少包括对扫描脑区获取的MEG数据进行数据带通滤波和下采样频率。
6.如权利要求1所述的方法，其特征在于，步骤Sb中所述基于提取的FMRI激活区坐标，提取与该激活区坐标对应的脑区的脑区时间序列，包括采用基于图像重建技术的最小范数解方法，对经过预处理的MEG数据进行脑磁逆问题求解，获得对应FMRI图像激活区坐标的脑区时间序列。
7.如权利要求1所述的方法，其特征在于，步骤&中所述依据提取的MEG脑区时间序列计算脑区间的因果连接强度和方向，是利用从MEG数据中提取的脑区时间序列，使用格兰杰因果方法分析时间序列之间的因果连接强度及方向。
全文摘要
本发明涉及一种联合功能磁共振成像(FMRI)和脑磁图描记术(MEG)的脑区因果连接检测方法，该方法包括以下步骤首先，对经过数据预处理的FMRI图像进行激活区坐标的提取；然后，基于提取的FMRI激活区坐标，在经过预处理的MEG数据上提取对应位置区域的脑区时间序列；最后，依据提取的MEG脑区时间序列，计算脑区间的因果连接强度和方向，利用有向网络图显示显著的连接。本发明所述方法是一种有效的联合FMRI和MEG两种成像模态的脑区因果连接检测方法，能够比仅利用FMRI图像的检测方法更完整准确地检测脑区因果连接。
文档编号G06T5/00GK102496159SQ20111039369
公开日2012年6月13日 申请日期2011年12月1日 优先权日2011年12月1日
发明者尤优博, 田捷, 白丽君 申请人:中国科学院自动化研究所