基于希尔伯特黄变换的脑磁图动态功能连接构建方法技术

本发明专利技术涉及基于希尔伯特黄变换法的脑磁图动态脑功能连接构建方法。该方法包括：1)脑磁图数据的预处理和源重构，脑区代表信号的选取；2)脑区信号的希尔伯特黄变换和脑区相干矩阵构建、瞬时脑网络构建；3)网络分析，包括瞬时脑网络分析、平均脑网络分析和聚类脑网络分析。采用希尔伯特黄变换元素共轭相乘相干变换构建的瞬时网络连接，完整保留大脑活动固有的非线性信息，且充分利用脑磁图高时间分辨率的特点，在时间上动态显示大脑活动过程，为了解神经疾病发生产生的影响、阐述认知过程、脑部手术进行术前/术后评估、脑机接口探索等领域提供帮助。提供帮助。提供帮助。

【技术实现步骤摘要】
基于希尔伯特黄变换的脑磁图动态功能连接构建方法


[0001]本专利技术涉及脑磁图功能网络构建
，尤其涉及一种基于静息态脑磁图数据的动态脑功能连接构建方法。

技术介绍

[0002]人类大脑由大约10^11个神经元细胞组成，神经元之间由约10^15个神经突触相互连接进而形成一个极为复杂的网络。从传统的细胞生物学角度去了解人脑功能、神经疾病原理、大脑发育与老化机制、认知的产生过程是远远不够的。迅速发展的医学影像分析技术证实了脑网络的存在，越来越多的神经解剖和神经生理研究聚焦于脑网络。然而，很多研究仅专注于网络连接的空间特征，忽略了网络的时间维度信息。
[0003]大脑活动是大脑响应外界刺激和适应外界活动的动态变化过程，该变化通常发生在毫秒级尺度。脑磁图是一种非侵入式脑功能信号采集设备,其获取的颅外磁信号可反映神经活动发生的位置和时间过程。脑磁图记录信号时间的精度可达毫秒，与大脑处理信息的时间尺度较为一致，因此，脑磁图为探究动态功能脑网络提供了潜在的有效工具。如能在脑磁图时间分辨率下构建并分析动态功能脑网络，将可能推动神经外科手术评估、神经系统疾病发病机制等研究的发展。
[0004]脑磁图动态功能连接方法中,以滑窗法为典型代表，该方法描述连续时间窗内功能连接的波动。时间窗内的功能连接构建，通常以信号间存在线性相关为前提，比如计算信号间幅度的相关性、考虑信息流方向的格兰杰因果关系、信号间的相干性等；或者考虑了大脑活动固有的非线性信息,比如相位耦合的相位锁相值、相位滞后指数等。上述方法中，在限定窗宽的时间尺度下,一定程度上丢失了大脑活动随时间变化的细节信息。同时，大量研究表明，大脑震荡产生的节律可在不同频带上形成脑网络，由于窗宽的限制，动态分析的频率范围将受限于窗周期。
[0005]以小波变换为代表的时频分析法很好地解决了基于滑窗法进行动态功能连接构建的不足。针对滑窗法中由于窗宽导致的动态分析频率范围受限的问题，小波变换借助丰富的小波基将信号分解至多频段，在研究过程中，我们不仅可以在单频带中评估功能连接，还可扩展至多频带甚至跨频带。在时间分辨率方面，小波相干支持信号在逐样本级别上的动态功能连接构建。
[0006]然而，尽管小波变换优于其他常见的时频变换方法(如短时傅里叶变换等)，但其在信号分解方面依赖于小波基函数的选取。由于基函数长度的限制，采用小波相干法构建功能连接时，无法敏锐地捕获信号的快速波动信息，在一定程度上也将导致信号能量谱的泄露。此外，加之小波变换的非自适应性无法体现信号本身的性质，因此，对于非线性非稳定且时空分辨率高的脑磁图信号来说，小波变换并不是最优的分析方法。综上所述，寻找一种能充分体现脑磁图信号本身特点、且不丢失信号信息和分辨率的方法，将在脑磁图动态功能连接构建上具有重要意义。

技术实现思路

[0007]针对现有方法特别是小波变换在脑磁图动态功能连接构建中存在的非自适应性、信息泄露、动态分辨率不足的问题，本专利技术提出一种基于希尔伯特黄变换的脑磁图动态脑功能连接构建方法。通过自适应信号分解，以瞬时相干来建立动态功能连接，可消除小波变换中基函数的选择对信号的限制，充分发挥脑磁图高时间分辨率的优势。此方法构建的动态功能连接能在时间上保留大脑活动过程的细节，且原理简单计算方便。为实现上述目的，本专利技术采用如下技术方案：
[0008]一种基于希尔伯特黄变换法的脑磁图动态功能连接构建方法，包括以下步骤：
[0009]步骤1，获取被试的脑磁图数据及相应被试的MRIT1图像。
[0010]步骤2，对原始脑磁图信号进行预处理，采用被试MRIT1图像构建头模型和源模型。
[0011]步骤2.1，对原始脑磁图数据进行检查，去除由于传感器故障，导致出现>800fT的高振幅或无振幅的通道，并用邻近通道的平均值/加权平均值来代替；采用0.5HZ高通滤波去除基线漂移，50Hz陷波器去除工频干扰，采用独立成分分析法去除信号中含有的肌电，眼电，心电这些干扰信号；对数据以等时长为间隔进行分割。
[0012]步骤2.2，对MRIT1图像进行分割提取脑表面，并建立头的体积传导模型；对MRIT1图像进行切片，构成基于皮质表面信息描述偶极子位置和方向的源模型。
[0013]步骤3，基于Beamformer法，对经过预处理的脑磁图信号、头模型、源模型进行源重构，创建通道信号与源的映射关系。
[0014]步骤4，根据解剖标记模板将大脑分为116个脑区，每个脑区有数量不等的经过源映射的通道信号。
[0015]步骤5，对各个脑区提取功率最大的信号作为脑区代表信号。
[0016]步骤6，对脑区代表信号进行希尔伯特黄变换。
[0017]步骤6.1，对脑区代表信号进行经验模态分解(EMD)。根据公式X(t)＝MF1(t)+IMF2(t)+
…
+MF
N
(t)+(t)，每个代表信号X(t)将自适应分解得到数量不一的本征模态函数IMF和一个残差函数r(t)，本专利技术中使用EMD需要满足的条件与其他人描述一致，通过计算各个本征模态函数与脑区代表信号的相关性，一致选取相关性较强的前几个本征模态函数进行下一步计算。
[0018]步骤6.2，将选取得到的本征模态函数进行希尔伯特变换。以g(t)代表当前进行希尔伯特变换的IMF(t)，希尔伯特变换公式为g(t)的希尔伯特分析变换值z(t)为(α(t)代表瞬时幅度，代表瞬时相位)。
[0019]步骤7，脑磁图动态功能连接构建。
[0020]步骤7.1，对经过希尔伯特黄变换的本征模态函数信号矩阵进行元素共轭相乘构建交叉谱其中和表示信号X1(t)和信号X2(t)经过希尔伯特黄变换后形成的本征模态函数信号矩阵。
[0021]步骤7.2，对交叉谱进行平滑和归一化其中E表示平滑，E
′
为E的导数，得到相干矩阵。可选地，计算交叉谱各时间点的相位值得到信号间的相位锁相矩阵。
[0022]步骤7.3，对任意两两感兴趣区求其相干矩阵/相位锁相矩阵，取出同一时间点的相干值/相位锁相值组成N
×
N的瞬时相干网络/瞬时相位锁相网络，将各个时间点的瞬时相干网络/瞬时相位锁相网络按时间顺序排列，即得到大脑动态功能连接网络。
[0023]与现有技术小波变换相比，本专利技术的效果是：
[0024]我们模拟了一对在时间和频率上具有动态一致性的时间序列，以研究与小波变换相比，我们的方法在脑磁图信号中的有效性。每个时间序列包含两个不同的频率部分，我们将这对时间序列模拟信号命名为1和2。实验中引入了显著性检验，以确保我们的连接性是有效的。在图7(a)中，两个模拟信号具有相似的功率谱分布；在图7(b)中，图谱分为两部分，左列是小波变换的结果，模拟信号1和2的小波谱为前两行所示，第三行表示它们的交叉谱，第四行是95％置信区间的结果。希尔伯特谱在右栏显示，其结果与小波结果相似。比较图7(b)的左右列，其结果显示基于小波变换和希尔伯特黄变换的频谱在时间<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于希尔伯特黄变换法的脑磁图动态功能连接构建方法，其特征在于，包括以下步骤：步骤1，获取被试的脑磁图数据及相应被试的MRIT1图像；步骤2，对原始脑磁图信号进行预处理，采用被试MRIT1图像构建头模型和源模型；步骤2.1，对原始脑磁图数据进行检查，去除由于传感器故障，导致出现&gt;800fT的高振幅或无振幅的通道，并用邻近通道的平均值/加权平均值来代替；采用0.5HZ高通滤波去除基线漂移，50Hz陷波器去除工频干扰，采用独立成分分析法去除信号中含有的肌电，眼电，心电这些干扰信号；对数据以等时长为间隔进行分割；步骤2.2，对MRIT1图像进行分割提取脑表面，并建立头的体积传导模型；对MRIT1图像进行切片，构成基于皮质表面信息描述偶极子位置和方向的源模型；步骤3，基于Beamformer法，对经过预处理的脑磁图信号、头模型、源模型进行源重构，创建通道信号与源的映射关系；步骤4，根据解剖标记模板将大脑分为116个脑区，每个脑区有数量不等的经过源映射的通道信号；步骤5，对各个脑区提取功率最大的信号作为脑区代表信号；步骤6，对脑区代表信号进行希尔伯特黄变换；步骤6.1，对脑区代表信号进行经验模态分解(EMD)；根据公式X(t)＝IMF1(t)+IMF2(t)+

【专利技术属性】
技术研发人员：杨春兰，刘露，李志梅，任洁钏，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：