CN110547772A - 一种基于脑信号复杂度的个体年龄预测方法 - Google Patents

Abstract

本发明涉及一种脑信号复杂度的个体年龄预测方法，包括：通过近红外成像系统获得被试者脑组织的透射光强数据；将透射光强数据转换为脑生理信息数据，并基于脑生理信息数据计算多尺度熵；计算每个被试者脑区域各检测点对应的多尺度熵值，将测试点按照各脑功能网络进行分类，再将表征同一脑功能网络中所有测试点的多尺度熵值取平均值，得到基于脑功能网络的多尺度熵指数；将每个脑功能网络的多尺度熵指数作为支持向量回归法的输入特征，使用线性核函数得到脑信号复杂度与年龄的相关性方程。上述方法能够通过静息态脑信号复杂度准确预定被试者年龄，测试过程简单，操作方便。

Description

一种基于脑信号复杂度的个体年龄预测方法

技术领域

本发明是关于一种基于脑信号复杂度的个体年龄预测方法，属于脑功能测试领域。

背景技术

人脑信息处理能力的快慢主要与脑功能网络状态密切相关。而脑功能网络状态主要通过脑信号复杂度表征。脑信号复杂度是作用于时间序列的一种指标，一段时间序列通过一个特定的复杂度算法处理，可以得到一个复杂度值，该值表现序列的非规则程度，一般地，不规则变化序列的复杂度较高，像常序列、周期序列这些模式比较简单的序列复杂度值较低。脑信号复杂度通常用标准差和多尺度熵(MSE)的方法来计算。其中，标准差反映的是整体时间序列的变化情况，且用于时间采样频率较低的成像方法；多尺度熵能评估不同时刻下脑信号的变化情况，由于其值是基于多种时间尺度计算得到的，因而多尺度熵应用于时间采样频率较高的成像方法，可以更好地描述脑信号复杂度。

虽然，不同年龄段人群在运动技能、身体协调性以及认知能力等方面存在很大差异。但脑网络状况影响因素众多，并不能直接地将脑网络状况与年龄进行对应。如果可以清楚确定脑网络状况与年龄的相互关系，则不但能进一步解释人脑发育过程，也能为研究发育障碍及儿科首发精神疾病提供理论支持。

发明内容

针对上述问题，本发明的目的是提供一种基于脑信号复杂度的预测个体年龄的方法，通过该方法建立了脑信号复杂度与个体年龄之间的明确关系，使通过脑信号复杂度预测个体年龄成为可能。

为实现上述目的，本发明采取以下技术方案：

本发明提供了一种基于脑信号复杂度的个体年龄预测方法，包括以下步骤：1)通过近红外成像系统获得被试者脑组织的透射光强数据；2)将所述透射光强数据转换为脑生理信息数据，并基于所述脑生理信息数据计算多尺度熵；3)计算每个被试者脑区域各检测点对应的多尺度熵值，将测试点按照各脑功能网络进行分类，再将表征同一所述脑功能网络中所有测试点的所述多尺度熵值取平均值，得到基于所述脑功能网络的多尺度熵指数；4)采用支持向量回归法进行年龄的预测：将每个所述脑功能网络的所述多尺度熵指数作为支持向量回归法的输入特征，使用线性核函数得到所述脑信号复杂度与年龄的相关性方程。

进一步，将被试者按照年龄段分为若干组，计算每组被试者的对应各脑功能网络的所述多尺度熵指数，并对其时间序列进行置换检验，以确认不同组之间被试者的所述多尺度熵指数存在明显差异。

进一步，所述脑生理信息数据包括氧合血红蛋白浓度和脱氧血红蛋白浓度。

进一步，对所述氧合血红蛋白浓度和脱氧血红蛋白浓度去除运动伪迹和生理噪声并进行带通滤波处理，采用所述氧合血红蛋白和脱氧血红蛋白浓度二者中量级大的一者计算多尺度熵。

进一步，所述脑功能网络包括：脑默认网络、额顶控制网络、腹侧注意网络、运动网络、背部注意网络和视觉网络。

进一步，步骤4)中通过留一法交叉验证评估支持向量回归法模型的正确性。

进一步，步骤4)中通过计算实际年龄与预测年龄的皮尔森相关系数评估所述脑信号复杂度与年龄的相关性方程的准确性。

进一步，所述近红外成像系统覆盖被试者脑的额叶、颞叶、顶叶和枕叶区域。

本发明由于采取以上技术方案，其具有以下优点：1、采用静息态脑信号复杂度测试，被试者无须进行知觉和行为的训练就能获取其脑功能数据，测试过程简单，操作方便，不会被试者身体造成伤害，尤其适用于低龄被试者，可用于儿童或青少年的脑信号复杂度检测；2、通过支持向量回归法建立了静息态脑信号复杂度与个体年龄的关系，使通过静息态脑信号复杂度预测个体年龄成为可能；3、通过置换检验、留一法交叉验证实际年龄与预测年龄的皮尔森相关系数,对每组被试者多尺度熵指数的差异和所述脑信号复杂度与年龄的相关性方程的可靠性进行了检验，保证了实验结果的可靠性。

附图说明

图1是本发明一实施例中各脑功能网络不同年龄组的多尺度熵时间序列；

图2是本发明一实施例中由支持向量回归法得出的多尺度熵指数与年龄的相关性方程。

具体实施方式

以下结合附图来对本发明进行详细的描绘。然而应当理解，附图的提供仅为了更好地理解本发明，它们不应该理解成对本发明的限制。在本发明的描述中，需要理解的是，术语仅仅是用于描述的目的，而不能理解为指示或暗示相对重要性。

实施例一

本发明的一个实施例提供了一种基于脑信号复杂度的个体年龄预测方法，包括：

1)通过近红外成像系统获得被试者脑组织的透射光强数据；

2)将透射光强数据转换为脑生理信息数据，并基于脑生理信息数据计算多尺度熵；

3)计算每个被试者脑区域各检测点对应的多尺度熵值，然后根据各检测点的位置计算出对应的各脑功能网络的多尺度熵值，再将表征同一脑功能网络中所有测试点的多尺度熵值取平均值，得到基于脑功能网络的多尺度熵指数；

4)采用支持向量回归法进行年龄的预测：将每个脑功能网络的多尺度熵指数作为支持向量回归法的输入特征，使用线性核函数得到脑信号复杂度与年龄的相关性方程。本实施例通过支持向量回归法建立了静息态脑信号复杂度与个体年龄的关系，使通过静息态脑信号复杂度预测个体年龄成为可能。

其中，步骤1)中获得被试者脑组织的透射光强数据的具体方式如下：首先，使用结构性核磁共振成像(MRI)对被试者的脑进行定位，然后根据定位的结果，在对应位置放置近红外成像系统的发射器和接收器。这里所说的对应位置与测试点对应。近红外成像系统至少应覆盖被试者脑的额叶、颞叶、顶叶和枕叶区域，采样频率设置为50Hz，采样时间约为11分钟。在采样过程中，被试需要闭上双眼，但不能睡着。需说明的是，采样频率和采样时间可以由本领域技术人员根据实际需要确定，本实施例只是给出一个优选值。

将近红外成像系统采集到的投射光强数据根据修正的Beer-Lambert定律转换为氧合血红蛋白(HbO)和脱氧血红蛋白(HbR)浓度，然后去除运动伪迹和生理噪声并进行带通滤波处理，由于HbO浓度信号在近红外成像设备中量级比HbR的浓度信号大，故选取HbO浓度信号用于后续静息态脑信号复杂度分析。

计算基于多尺度熵的脑信号复杂度方法包括如下步骤：

首先，将浓度信号的时间序列{x₁,…,x_i,…,x_N}进行“粗粒化”操作从而得到降采样的时间序列，粗粒化的时间序列y(t)的计算公式为：

N为评估长度，j为粗粒化时间序列的每一个值对应的下标，t为时间尺度大小；本实施例中取t为1-25之间的整数(包括1和25)，t取1时指采用原始时间序列，t取25时指每个窗长中包括25个数据点。

然后，对每个粗粒化的时间序列y(t)通过切比雪夫距离计算其对应的样本熵，样本熵的计算公式为：

m为模式长度，设m＝2，表示使用两连续的数据点进行模式匹配，相似准则r设为0.5，表明到数据点之间的绝对振幅小于时间序列标准的50％时，数据点被认为具有不可区分的振幅。C_m(r)表示满足时间点对i,j占所有时间点对的概率。样本熵通过评估长度为N的时间序列上振幅模型的可预见性来量化信号变化情况。多尺度熵值高表示脑信号复杂度高，多尺度熵值低表示脑信号复杂度低。

根据计算基于多尺度熵的脑信号复杂度方法计算每个被试者脑区域各测试点对应的多尺度熵值，然后根据各测试点位置确定其对应的脑功能网络，此处的脑功能网络包括：脑默认网络、额顶控制网络、腹侧注意网络、运动网络、背部注意网络及视觉网络。这一步就是将测试点按照脑功能网络进行分类，确定每个脑功能网络对应的测试点。再将所有与某一脑功能网络对应的测试点的多尺度熵值取平均得到基于该脑功能网络的多尺度熵指数，以用于后续多尺度熵指数和年龄的相关性的研究。

采用支持向量回归法(SVR)进行年龄的预测，该方法包括将每个脑功能网络的多尺度熵值作为给定样品集进行输入，假定在线性可分训练集中，存在一个超平面f(x_i)＝(w·x)+b实现拟合，那么训练样本可分为两类，即y_i-(w·x)+b＜0的训练样本和y_i-(w·x)+b＞0的训练样本，在回归问题中，理想状态y_i-f(x_i)＝0事实上几乎不可能存在，通常，允许在一个小的正数ε，使|y_i-f(x_i)|≤ε，则认为回归是无损的。选择合适精度的ε，求解对偶化最优问题，计算w，选择拉格朗日乘子a_i>0时，计算b，从而构造线性SVR超平面，该SVR超平面对应的方程即为多尺度熵值与年龄的对应关系方程。在本实施例中ε设为0.001。然后通过留一法交叉验证评估支持向量回归法模型的表现，通过计算实际年龄与预测年龄的皮尔森相关系数评估多尺度熵值与年龄的对应关系的方程的准确性。

为了确定不同的年龄组的脑功能网络即多尺度熵指数是否存在差别，被实施例还包括如下验证步骤：将被试者按照年龄段分为若干组，计算每组被试者的对应各脑功能网络的多尺度熵指数，并对其时间序列进行置换检验，以确认不同组之间被试者的多尺度熵指数存在明显差异。

置换检验过程如下：为了验证不同年龄段的多尺度熵值存在显著差异，设定A组为一个年龄段的多尺度熵值，B组为另一个年龄段的多尺度熵值；零假设为不同的年龄段多尺度熵值不存在差别，在这个检验中，若零假设成立，那么A组和B组的数据分布应该是一样的，即A组和B组服从同一分布。接下来构建检验统计量——A组多尺度熵值的均值与B组多尺度熵值的均值之差。将A组和B组合并到一个集合，随机选择一半的数据作为A组，另一半的数据作为B组。计算并记录第一步中A组和B组的多尺度熵值均值之差；重复随机选择A、B组数据，并求各组均值的步骤999次，重复的次数越多，得到的背景分布越稳定。计算原始A组多尺度熵值的均值与B组多尺度熵值的均值的差值，并求出重复的999次中出现该差值的次数，假设该差值的出现次数为a，显著性值P＝a/999,当显著性值P小于显著性阈值时，零假设不成立，即不同的年龄段多尺度熵值存在差别。通常的显著性阈值选为0.05，但技术人员也可以根据具体试验要求的精确性确定显著性阈值。

实施例二

本实施例通过列举具体实验数据对本发明中技术方案做进一步说明。具体步骤如下：

使用波长为670nm，830nm的近红外成像设备采集107个健康被试者静息态近红外成像光强数据，根据年龄将被试者分为以下三组：6-7岁，8-10岁和11-13岁，然后将红外成像光强数据转换得到HbO浓度数据，然后去除运动伪迹和生理噪声并进行带通滤波处理，使用处理后的HbO浓度数据计算出每个测试点对应的多尺度熵值。

每组被试的样本熵值的时间序列如图1所示，每一组样本熵值均随时间的增加而增加；对于每个时间尺度，年龄大的群组均比年龄小的群组的多尺度熵值大并且不同年龄群组之间没有相互交叠、相互影响的情况。

计算各脑功能网络的多尺度熵指数，然后采用置换检验计算多尺度熵指数与年龄的相关性，得出六个脑功能网络中，除了视觉网络外，其他脑功能网络都呈现年龄与多尺度熵的显著相关。

利用SVR分析方法根据脑功能网络的多尺度熵向量预测被试的年龄并使用留一交叉验证方法进行验证，如图2所示，根据上述多尺度熵值计算实际年龄与预测年龄之间关系得出脑网络相关系数r＝0.56，置换检验p＝0.0000；额顶控制网络相关系数r＝0.41，置换检验显著性值p＝0.0000；腹侧注意网络相关系数r＝0.28，置换检验显著性值p＝0.0014；运动网络相关系数r＝0.32，置换检验显著性值p＝0.0005，以上数据说明脑信号复杂度能较好地预测年龄。

上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。

Claims (8)

1.一种基于脑信号复杂度的个体年龄预测方法，其特征在于，包括以下步骤：

1)通过近红外成像系统获得被试者脑组织的透射光强数据；

2)将所述透射光强数据转换为脑生理信息数据，并基于所述脑生理信息数据计算多尺度熵；

3)计算每个被试者脑区域各检测点对应的多尺度熵值，将测试点按照各脑功能网络进行分类，再将表征同一所述脑功能网络中所有测试点的所述多尺度熵值取平均值，得到基于所述脑功能网络的多尺度熵指数；

4)采用支持向量回归法进行年龄的预测：将每个所述脑功能网络的所述多尺度熵指数作为支持向量回归法的输入特征，使用线性核函数得到所述脑信号复杂度与年龄的相关性方程。

2.如权利要求1所述的基于脑信号复杂度的个体年龄预测方法,其特征在于，将被试者按照年龄段分为若干组，计算每组被试者的对应各脑功能网络的所述多尺度熵指数，并对其时间序列进行置换检验，以确认不同组之间被试者的所述多尺度熵指数存在明显差异。

3.如权利要求1或2所述的基于脑信号复杂度的个体年龄预测方法,其特征在于，所述脑生理信息数据包括氧合血红蛋白浓度和脱氧血红蛋白浓度。

4.如权利要求3所述的基于脑信号复杂度的个体年龄预测方法,其特征在于，对所述氧合血红蛋白浓度和脱氧血红蛋白浓度数据去除运动伪迹和生理噪声并进行带通滤波处理，采用所述氧合血红蛋白和脱氧血红蛋白浓度二者中量级大的一者计算多尺度熵。

5.如权利要求1或2所述的基于脑信号复杂度的个体年龄预测方法,其特征在于，所述脑功能网络包括：脑默认网络、额顶控制网络、腹侧注意网络、运动网络、背部注意网络和视觉网络。

6.如权利要求1或2所述的基于脑信号复杂度的个体年龄预测方法,其特征在于，步骤4)中通过留一法交叉验证评估支持向量回归法模型的正确性。

7.如权利要求1或2所述的基于脑信号复杂度的个体年龄预测方法,其特征在于，步骤4)中通过计算实际年龄与预测年龄的皮尔森相关系数评估所述脑信号复杂度与年龄的相关性方程的准确性。

8.如权利要求1或2所述的基于脑信号复杂度的个体年龄预测方法,其特征在于，所述近红外成像系统覆盖被试者脑的额叶、颞叶、顶叶、枕叶区域。