一种基于功能磁共振成像的脑区效应连接分析系统

摘要

本发明公开了一种基于功能磁共振成像的脑区效应连接分析系统，包括实验范式设计模块、核磁成像处理模块、大脑图像分区模块、效应连接检测模块和分析结果显示模块。实验范式设计模块对实验对象参与的实验内容和实验流程进行设计；核磁成像处理模块通过参数配置单元对磁共振设备采集到的图像进行预处理并保存；大脑图像分区模块加载和储存大脑分区模板，对核磁成像进行区域划分；效应连接检测模块进行定向线性因果关系分析，计算效应连接值和连接指标；分析结果显示模块将效应连接分析结果和相关的脑功能指标反馈给医护人员和被试。本发明在功能磁共振成像质量评估、脑功能诊断与调节、认知功能研究、精神疾病治疗等领域具有重要的应用价值。

说明书

技术领域

本发明涉及一种医学影像分析系统，具体是一种基于功能磁共振成像的脑区效应连接分析系统，属于生物医学信息处理技术领域。

背景技术

人脑是自然界中最复杂的系统之一。一个成年人的大脑中，大约有10¹¹个神经元细胞，这些数量巨大的神经元细胞通过大约10¹⁵个突触交错连接，形成了一个高度复杂的脑神经系统。虽然不同神经元、不同脑区对应不同的功能特异性，但即使是完成一项非常简单的任务，也需要大脑多个不同的区域相互作用、相互联系，构成一个网络来协调工作。

近年来，一些科学家充分认识到了构建人脑网络动态模型的重要性，开展了一系列针对大脑功能连接的研究，力图在广泛的时空尺度上将大脑连接成动态的功能网络，借助先进的技术手段进一步挖掘该网络的连接规律。脑电图（Electroencephalogram, EEG）是目前最敏感的脑功能监测方法，空间分辨率只能达到厘米级别，并且极易受各种因素干扰，多数情况下不能作为大脑功能网络研究的唯一依据。

脑影像学技术能够从微观和宏观水平上分析局部和整体大脑功能网络特性的变化。基于血液氧合依赖水平（Blood Oxygen Level Dependent, BOLD）的功能磁共振成像（Functional Magnetic Resonance Imaging, fMRI）技术利用磁共振造影来测量神经元活动所引发的氧合血红蛋白含量变化，不仅时间分辨率更高，其空间分辨率也可达到毫米水平。作为一种非损伤性脑成像技术，磁共振成像在脑功能研究中发挥了不可替代的作用。借助磁共振成像技术可以将对大脑的研究扩展至记忆力、注意力等，识别研究对象感觉到的声音、图像或者阅读的词语、执行的动作，甚至能够实现对个人内心世界的揭示。

脑功能网络从整体的角度克服了以往孤立研究单个脑区的缺陷，但是目前的脑网络研究大多集中在对无向网络功能连接的研究上。无向连接分析起来比较简单，但是忽略了大脑神经系统中非常重要的方向和权重信息。因此，构建有向的功能网络来刻画不同脑区神经活动之间的因果关系，更深入细致的理解大脑的结构组织模式和功能活动规律成为大脑网络研究的当务之急。

在大脑网络研究中，具有方向性的功能连接称为效应连接。效应连接体现了空间上分离的神经单元（脑区）的神经活动在时间上的关联性或统计依赖关系，通过量化神经单元之间的直接或间接影响，反映神经活动的动态过程以及实验因素对神经活动的调节作用，更接近于真实的脑功能机制。大脑功能连接模式的研究目前只是停留在体验式实验的设计与分析阶段，对脑功能网络连接机制、尤其是对脑区效应连接机制的模型化和机理性研究相对较少。同时，在当前的认知研究和临床应用中，也非常缺乏与先进的功能磁共振成像技术相配套、满足实际应用需要的综合性脑区效应连接分析系统。

发明内容

针对现有技术存在的不足之处和实际应用的需要，本发明要解决的问题是：

提供一种基于功能磁共振成像脑区效应连接分析系统，实现对大脑功能区域之间的效应连接进行判定和记录。

为了达到上述目的，本发明采取以下技术方案：

一种基于功能磁共振成像的脑区效应连接分析系统，包括实验范式设计模块、核磁成像处理模块、大脑图像分区模块、效应连接检测模块和分析结果显示模块：

所述实验范式设计模块用于对实验对象（被试者）参与的实验内容和流程进行设计；

所述核磁成像处理模块包括图像读取单元、格式转换单元、图像预处理单元和参数配置单元，通过参数配置界面对磁共振设备采集到的大脑三维结构像和时间层功能像进行读取和格式转换，执行预处理操作；

所述大脑图像分区模块用于加载和储存各种标准化的大脑分区模板，将操作人员选定的模板与磁成像处理模块预处理后的核磁图像进行匹配，划分图像为若干大脑区域； 

所述效应连接检测模块包括时间序列提取单元、因果连接分析单元和连接指标计算单元，用于提取大脑功能磁共振图像中不同标准分区对应的时间序列，对一组实验对象的时间序列进行单样本T检验（One sample T-test），将得到的组分析结果进行定向线性因果关系分析，计算不同脑区时间序列之间的效应连接值和效应连接指标；

所述分析结果显示模块将效应连接分析结果和相关的脑功能指标反馈给医护人员、被试和操作人员。

进一步，所述核磁成像处理模块将读取到的磁共振图像由DICOM格式转换为NIFTI格式。

进一步，所述核磁成像处理模块执行的图像预处理操作包括功能像头动校正、结构像配准、结构像灰质分割、结构像和功能像标准化、功能像平滑等预处理操作。

进一步，所述效应连接检测模块采用多变量自回归（Multivariate Autoregressive, MVAR）模型测量两组时间序列之间的定向线性因果关系。根据脑区时间序列建立的p阶MVAR模型表达式可表示为：

，

其中，u[n]、v[n]及w[n]分别为时间序列x[n]、y[n]及其组合向量（矩阵z[n]）的自回归估计残差。在上述MVAR模型中，残差的方差或方差矩阵可表示为：

，，。

根据MVAR模型得到的定向线性因果关系可以表示为：

。

其中，F_x→y表示以x[n]为参考区域指向y[n]的线性定向影响；F_y→x表示从y[n]到x[n]的线性定向影响。

如果x[n]的过去值能够预测y[n]的当前值，F_x→y的值就会明显大于0 (T₁ << M₁)；反之，F_1→2值近似等于0 (T₂ ≈ M₂)，F_y→x的情况类似。F_x→y与F_y→x差值的符号代表因效连接方向，差值大小表示因效连接的权值；由于T₁ ≤ M₁、T₂ ≤ M₂，F_x→y和F_y→x总是非负，需要将它们的值进行归一化；然后采用错误发现率校正方法（FDR）确定其显著性水平的阈值，并将结果映射到每个脑区，形成大脑因效连接网络。

进一步，所述分析结果显示模块显示的脑功能指标主要包括脑区因果连接图、加权因果连接图、出\入度、加权出\入度、聚类系数、因果密度等。

采用上述技术方案后，本发明的有益技术效果是：

（1）通过量化磁共振成像时间序列之间的作用程度，根据节点间的有向连接强度判定和描述不同脑区在统计意义上的因果关系；

（2）通过分析体现生理和心理信息传递方向的有向连接机制，探寻大脑功能特征参数变化的起因，并完善大脑各个区域相互作用所形成的拓扑结构；

（3）通过分析脑区之间协同工作模式和调控与被调控关系，发现大脑特定区域在体现抑郁症发病机理和症状表现的相应有向连接中所起到的作用，以更好的解释患者大脑内的信息加工和精神表达机制；

（4）在功能磁共振成像质量评估、脑功能诊断与调节、认知功能研究、精神疾病治疗等多个领域具有重要的应用价值。

附图说明

图1是本发明一种基于功能磁共振成像的脑区效应连接分析系统的结构示意图；

图2是本发明中核磁成像处理模块的结构示意图；

图3是本发明中效应连接检测模块的结构示意图；

图4是本发明一种基于功能磁共振成像的脑区效应连接分析系统的操作流程图。

具体实施方式

下面结合附图和具体实施例对本发明作进一步的说明。

本发明一种基于功能磁共振成像的脑区效应连接分析系统的结构如图1所示，包括实验范式设计模块、核磁成像处理模块、大脑图像分区模块、效应连接检测模块和分析结果显示模块。本实施例中，系统平台采用Windows操作系统下的图形化集成开发环境Visual C# 2010实现，数据存储与处理载体为64位计算机，人机交互界面载体为液晶显示屏或触摸屏。

操作人员通过实验范式设计模块对实验对象（被试）参与的实验内容和实验流程进行设计。本实施例中，系统采取随机、对照和重复的组块实验方案，激活状态和控制状态两个组块交替进行。

如图2所示，核磁成像处理模块包括图像读取单元、格式转换单元、图像预处理单元和参数配置单元。操作人员首先通过参数配置单元进行参数设置后，在图像读取单元中对磁共振设备采集到的大脑三维结构像和时间层功能像进行读取，再在格式转换单元中将读取到的结构像和功能像由DICOM格式转换为NIFTI格式，然后在图像预处理单元中进行头动校正、配准、分割结构像、标准化、平滑等预处理操作，最后核磁成像处理传递至效应连接检测模块。本实施例中，核磁成像处理模块执行的预处理操作步骤依次为：功能像头动校正、结构像配准、结构像灰质分割、结构像和功能像标准化、功能像平滑。 

大脑图像分区模块用于加载和储存各种标准化的大脑分区模板，如AAL分区模板、Brodmann分区模板、CH2分区模板等。操作人员将选定的模板与磁成像处理模块得到的图像进行匹配，划分图像为若干大脑区域，还可以依据先验知识选取部分脑区或一些功能相关脑区组成的社团组合，最终将大脑图像分区结果传递至效应连接检测模块。本实施例中，系统默认选择被划分成90个脑区的AAL大脑分区模板进行全脑效应连接分析。

如图3所示，效应连接检测模块包括时间序列提取单元、因果连接分析单元和连接指标计算单元。操作人员对效应连接检测模块的检测选项进行选择，通过时间序列提取单元提取大脑不同标准分区对应的时间序列，再通过因果连接分析单元对一组实验对象的时间序列进行统计分析，将得到的组分析结果进行定向线性因果关系分析，计算不同脑区时间序列之间的效应连接值，然后通过连接指标计算单元计算出脑区效应连接的各项指标；再将计算结果传递至分析结果显示模块。本实施例中，系统提取预处理后磁共振图像中各脑区的体素平均值作为大脑血氧依赖水平(BOLD)响应的时间序列。

效应连接检测模块采用MVAR模型测量两组时间序列之间的定向线性因果关系。根据脑区的时间序列建立的p阶MVAR模型表达式可表示为：

，

其中，u[n]、v[n]及w[n]分别为时间序列x[n]、y[n]及其组合向量（矩阵z[n]）的自回归估计残差。

上述MVAR模型中，残差的方差或方差矩阵可表示为：

，，。

根据MVAR模型得到的定向线性因果关系可以表示为：

。

其中，F_x→y表示以x[n]为参考区域指向y[n]的线性定向影响；F_y→x表示从y[n]到x[n]的线性定向影响。

如果x[n]的过去值能够预测y[n]的当前值，F_x→y的值就会明显大于0 (T₁ << M₁)；反之，F_1→2值近似等于0 (T₂ ≈ M₂)，F_y→x的情况类似。F_x→y与F_y→x差值的符号代表因效连接方向，差值大小表示因效连接的权值；由于T₁ ≤ M₁、T₂ ≤ M₂，F_x→y和F_y→x总是非负，需要将它们的值进行归一化；然后采用错误发现率校正方法（FDR）确定其显著性水平的阈值，并将结果映射到每个脑区，形成大脑因效连接网络。

当x[n]过去值的相关信息丢失时，以前能被x[n]过去值解释的z[n]的方差，现在只能用y[n]的过去值来解释，这就导致了F_y→x测量值的波动。本实施例中，为了避免测量偏差的干扰，采用F_x→y和F_y→x的差值来考察时间序列之间的因果关系。当F_x→y—F_y→x>0时，表示x[n]对y[n]存在正向影响；F_x→y—F_y→x<0表示y[n]对x[n]的影响为正向。

操作人员可以对分析结果显示模块的显示选项进行选择，分析结果显示模块将效应连接分析结果和相关的脑功能指标反馈给医护人员和被试和操作人员。分析结果显示模块显示的大脑功能因效连接指标主要包括脑区因果连接图、加权因果连接图、出\入度、加权出\入度、聚类系数、因果密度等。本实施例中，系统显示的因效连接指标为加权因果图、加权出\入度和因果密度。

如图4所示，本发明一种基于功能磁共振成像的脑区效应连接分析方法包括以下步骤：

（1）操作人员通过实验范式设计模块对实验对象（被试）参与的实验内容和实验流程进行设计；

（2）操作人员通过核磁成像处理模块中的参数配置界面设置参数，对磁共振设备采集到的图像进行读取和格式转换，再进行头动校正、配准、分割结构像、标准化、平滑等预处理；

（3）操作人员通过大脑图像分区模块选定一种标准化的大脑分区模板，并与预处理后的核磁图像进行匹配，划分图像为若干大脑区域；

（4）操作人员通过效应连接检测模块提取大脑功能磁共振图像中不同标准分区对应的时间序列；

（5）操作人员通过效应连接检测模块对该组实验对象的脑区时间序列进行统计分析，计算不同脑区时间序列之间的效应连接值和效应连接指标；  

（6）操作人员通过分析结果显示模块将效应连接分析结果和相关的脑功能指标反馈给医护人员和被试，如果医护人员要求更换其它脑区标准模板进行效应连接检测分析，转至步骤（3）；否则执行步骤（7）；

（7）操作人员收到医护人员关于实验重新开始的通知，转至步骤（1）；否则执行步骤（8）；

（9）结束本次分析。

本发明还可以按照如下方式进行，以达到更有效或更有针对性的效应连接分析效果：

（1）操作人员可以在大脑图像分区模块存储或加载的标准化大脑分区模板中，依据先验知识选取部分脑区或一些功能相关脑区组成的社团组合，并与核磁成像处理模块预处理后的图像进行匹配。在大脑分区结果传递至效应连接检测模块后，提取所选定区域对应的时间序列进行效应连接分析。

（2）操作人员除了采用大脑图像分区模块自带的标准化大脑分区模板，还可以在大脑结构像中自由选取若干个感兴趣区域（Region of Interest, ROI），并与核磁成像处理模块预处理后的图像进行匹配，在效应连接检测模块中提取各个感兴趣区域对应的时间序列进行效应连接分析。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。