论文部分内容阅读
作为一种无损伤、不用外源介质和非侵入性的测量技术,核磁共振成像(MRI)为脑功能的研究提供了一种重要的手段。尽管功能核磁共振成像(fMRI)的历史并不长,但是发展却很快,正在逐步成为视觉、语言、工作记忆和其它认知过程脑功能研究的重要工具。本文首先利用非线性动力学的方法检测静息态人脑fMRI数据的时空非线性特性。采用梯度回波快速成像方法在1.5T的核磁共振成像仪上采集九个健康被试的静息态fMRI数据,并通过相关矩阵的特征值谱分析、关联维数分析和时空Lyapunov指数分析来检测静息态人脑fMRI数据的时空非线性特性。通过模拟、调节AR(1)结构纯噪声数据的特征值谱,并将调节后的特征值谱与fMRI数据的特征值谱比较,获得fMRI数据的本征维数的估计。依据所估计的本征维数,用主成分分析从fMRI数据中抽取若干主成分,这样在保证保留准确的相动力学的基础上,尽可能减轻计算量和降低噪声。然后用时间延迟嵌入的方法对所取得的主成分进行相空间重构,并用Grassberger-Procaccia算法估计多变量序列的关联维数。与此同时,对静息态fMRI数据应用耦合映象格子(CML)方法计算时空Lyapunov指数(STLE)及其嵌入维效应和时间演化。非线性统计显著性检验结果表明,在静息态fMRI数据和相应的surrogate数据之间,关联维数和时空Lyapunov指数存在着明显的差异。而且分数关联维数和正的时空Lyapunov指数刻画了静息态人脑fMRI数据的时空非线性动力学特性。这些结果表明,静息态人脑fMRI数据中存在着非线性结构,因此,结果提示静息态fMRI的波动反映了人脑中内在的基本神经活动模式,不能完全归因于噪声。另一方面,时空Lyapunov指数的嵌入维效应及其时间演化也表明,静息态人脑fMRI数据中存在确定性的非线性行为以及动力学特性的短暂稳定性。而正的时空Lyapunov指数还表明静息态人脑中存在时空混沌现象,而时空混沌现象提示体素之间的相关性并不是固定不变的,而随着时间变化,这表明静息态人脑中存在动态功能连接和动态脑功能网络。估计fMRI数据中真实的本征维数(即真实的信号数目)是一个非常重要同时又是非常困难的问题。在本文中,通过构建合适的插值变量,获得了对一阶自回归噪声模型中自回归系数更合理的估计。采用基于一阶自回归噪声模型和三次样条插值相结合的本征维数估计新方法(AR1CSI)、基于一阶自回归噪声模型的本征维数估计方法(AR1)和基于分形的本征维数估计法(FB),分别对仿真数据和静息态fMRI数据进行本征维数估计,结果表明,无论在不同的体素数目、时间长度、信号数目和信噪比还是不同噪声模型的情况下,AR1CSI的本征维数估计方法都比AR1方法及FB方法有更准确的信号数目估计。因此,AR1CSI本征维数估计方法是一种性能更为优异的本征维数估计方法。传统脑功能激活图的检测方法是基于广义线性模型的,尽管该方法取得了许多重要的结果,由于该方法假定神经刺激和fMRI响应之间呈现线性关系,因而存在着明显的缺陷。当今以聚类分析为代表的数据驱动方法受到了日益广泛的重视。本文将基于数据点消息传递(CPMDP)的聚类分析方法引入fMRI数据的分析。为了提高聚类分析结果的可靠性,提出增加一个附加体素的思想。在fMRI数据分析中,用所有体素时间序列与神经刺激标准响应之间的相关系数来确定附加体素。对混合数据集和听觉刺激fMRI数据进行CPMDP聚类分析,并与k-means聚类分析方法相比,结果表明,CPMDP分析不仅不需要预先规定聚类的类数,而且聚类结果也稳定得多,同时由于增加了附加体素,聚类的可靠性也提高了。因此CPMDP聚类方法是一种性能优秀的fMRI数据分析方法,其性能优于k-means聚类分析方法。