脑功能磁共振成像数据处理算法及应用研究

发布时间：2020-06-09 18:39

【摘要】： 高空间分辨率的功能磁共振成像(Functional magnetic resonance imaging, fMRI)技术是人们认识脑的一种重要的无损检测手段,并且在相关领域受到了极大关注。fMRI脑功能信息提取技术是将fMRI与大脑认知、神经科学和临床应用相结合的关键。目前该技术仍然不完善,需要进一步的发展。论文以脑功能活动区域定位、聚类的fMRI数据处理算法和脑功能不对称性等问题为核心,提出了利用fMRI时间序列的频率信息、相位谱(Phase spectrum,PS)信息、卷积功率谱(Convolution power spectrum,CPS)信息和仿射聚类(Affinity propagation clustering,APC)等方法来探测脑功能活动,发展了相应的抗干扰能力强的时频fMRI信号分析法,并进一步提出了利用定量分析功率谱变化来研究运动皮层功能不对称性的思想。通过仿真与实际数据的分析证明了这些方法和技术创新的有效性和应用价值。分列如下: 时域fMRI数据处理方法常忽略时间序列的频率信息,本论文提出了一种利用时间序列在希尔伯特(Hilbert)空间的频率信息去探测功能活动区的新方法。它的主要思想是利用任务状态与非任务状态的频率熵信息(Frequency entropy information,FEI)的差作为判断脑激活的指标。通过仿真确定了建议方法的有效性。比较统计参数图(Statistical parametric mapping,SPM)、支持向量机(Support vector machine,SVM)和FEI等算法的接受者工作特性曲线(Receiver operating characteristic,ROC),频率信息方法显示出优良性能,实际的fMRI数据处理也证实该方法能有效探测脑功能活动。提出了利用fMRI时间序列的相位谱信息来探测脑功能活动的新方法。基本思想是将信号特征频率上的相位用作脑激活指标,并且通过仿真和实际fMRI数据验证和确认了这种方法的有效性。常规fMRI数据分析方法以构造体素的血氧水平依赖性(Blood oxygen level-dependent,BOLD)响应时间函数为目标,但是功率谱分析理论完全聚焦在理解作用系统的动力学能量变化。因此,针对fMRI数据处理,论文提出了一种新的基于先验镜像信号的卷积能谱分析模型。首先,抑制fMRI信号中的噪声,计算测量的fMRI数据与预先设计的实验模式之间的卷积;然后将卷积信号的功率谱密度(Power spectrum density, PSD)分析作为定量分析BOLD信号变化的指标。并且对比常规的功率谱、SPM和SVM方法,仿真和实际fMRI块实验数据分析结果揭示CPS方法能有效地探测脑功能活动。对于具有复杂特性的fMRI时间序列,论文的结果证明抗噪声能力强的CPS方法是揭示脑功能信息的有用的分析工具。聚类分析是一种很好的的数据驱动方法。然而对于fMRI时间序列数据,庞大的计算量使它很难实用。基于此,论文提出了将有效、快速的新聚类算法—仿射聚类APC用于分析脑功能活动的fMRI大数据集。该算法不用随机选择初始的类代表点,它将所有的数据点都作为潜在的类代表点,通过最小化能量函数与信息传递架构(Message passing architecture),得到最优化的类代表点与它们对应的类。四个仿真数据集和三个真实fMRI数据集(包括块设计与事件相关(Event-related)实验)的分析结果揭示APC能有效地探测脑功能活动并能区分不同的响应模式。性能参数指标平均均方误差(Average squared error)揭示在块设计与事件相关实验中APC算法明显优于k中心聚类(k-centers clustering)。研究结果表明利用fMRI数据,APC算法能有效地探测和聚类脑功能活动区。本论文提出了利用定量分析BOLD信号的功率变化来研究脑功能不对称性的思想。六个右利手被实验者参与了fMRI实验,论文分析了双手运动与单手运动的实验数据。功率谱方法表明在任务与休息状态之间,大脑右侧运动皮层相对左侧运动皮层,BOLD信号有更大的功率差,这说明其右侧运动区域有更多相关神经细胞被诱发。此外,通过对比BOLD信号幅度变化,对用功率谱方法测量功能不对称性进行了合理性分析。对于脑功能不对称性分析,功率谱方法确信是有效的定量分析方法。
【图文】：

时间序列,头像,子区域

产生一个 128×128 个体素的仿真 fMRI 脑图像，如图2-1 所示。图 2-1 由五个不同尺寸与形状的活动子区域（471 个体素），11026个体素的时不变黑边和信号随机变化的 4887 个体素（灰质，白质与脑室）构成。图 2-1 仿真的 fMRI 头像块设计与事件相关（Event-related）实验：两个仿真数据集 BS1 和 BS2（图 2-2（c）和图 2-3（c））分别构建五个已知位置的活动子区域的数据。设计数据集 BS1（图 2-2（c））用来说明对仿真块数据的血液动力学响应的形状变化。三种不同的活动模式分别叠加到区域“A”和“E”（“A_E”），区域“B”和“C”（“B_C”）与区域“D”。空间分离但时间模式相同的区域“A” 和“E” 被合并为相同的激活源区（“A_E”），类似的在我们的研究中，“B” 和“C” 也被合并为相同的激活源区（“B_C”）。子区域体素的时间序列由信号与高斯噪声混合构成，SNR=0.5。被设计的方块信号如图 2-2（c）第一排所示。数据集 BS2（图 2-3（c））设计来说明事件相关实验。活动模式分别叠加到 5 个激活区域，并且 SNR 与 BS1 一致。为了检验FEI的效果，仿真数据也用SVM方法来进行对比分析。我们用线性核SVM

设计实验,功能,范围,方法

2-2 在 SNR = 0.5 下仿真块设计实验的探测结果：（a）SPM 计算的功能SVM（中）和 FEI 方法（右）（p-value <0.05，FWE 校验）;（b）三种线（p-value 范围从 0.00001 到 0.05，，FWE 校验）;（c）BS1 的三种实
【学位授予单位】：电子科技大学
【学位级别】：博士
【学位授予年份】：2010
【分类号】：R445.2

【引证文献】