一种基于遗传算法的fMRI约束独立成分分析学习方法

本发明属于脑影像学图像处理，尤其是涉及一种基于遗传算法的fmri约束独立成分分析学习方法。

背景技术：

1、大脑作为人类最重要的器官，它的重要性不仅体现在控制着人类的思维、意识、情感、记忆等各种认知行为，还体现在可以帮助人类实现不同于其他动物的高级认知功能，也是人类迄今为止所知道的最为复杂和精密的系统之一。尽管如此，对大脑神经活动认知机制的探索，一直是国内外科学界力求攻破的一道难题，具有十分重要的研究价值。近年来，随着现代神经影像技术的发展，脑成像技术已成为神经科学中一项最直观的观察手段。其中，功能磁共振成像(functional magnetic resonance imaging,fmri)技术主要利用含氧血红蛋白与脱氧血红蛋白的磁敏感性差异来进行成像，其凭借无侵入、无创伤、无辐射、可重复、可精确定位以及较高时间和空间分辨率等诸多优越特性，广泛应用于认知神经科学、神经心理学、临床医学等各领域的脑科学研究。

2、基于fmri脑功能的多目标约束独立成分分析(constrained independentcomponent analysis,cica)是利用该方法进行脑科学研究的一项关键技术。但是对于多目标cica模型的求解，利用传统线性加权求和方法以及快速不动点算法还存在一些不足，线性加权求和方法需要人为经验设置权重值，快速不动点算法的迭代速度较慢，并且迭代的方向具有不确定性，很容易陷入局部最优的处境。因此，优化和完善多目标cica的求解方法，通过遗传算法自适应地迭代出cica模型的全局最优解集具有重要意义，不仅能够解决人为设置阈值参数的弊端，还能快速自适应地迭代出分离向量矩阵的全局最优解集，提高感兴趣区独立成分提取的精度。然而，当前的遗传算法也可能受到迭代次数和个体遗传过程中的概率性事件的影响，从而导致迭代的解集小幅度变差。

技术实现思路

1、本发明的目的在于提供一种基于遗传算法的fmri约束独立成分分析学习方法。通过cstica方法搭建一个多目标优化框架，然后将每个目标函数中的矩阵通过行或列向量求内积将矩阵迭代转换成多变量单函数进行迭代，最后利用非支配排序遗传算法求解出分离矩阵的全局最优解集，进而还原出相对应的独立源信号。

2、为了达到上述目的，本发明通过以下技术方案实现：

3、一种基于遗传算法的fmri约束独立成分分析学习方法，包含如下步骤：

4、s1、采集若干正常健康被试的静息态fmri脑影像数据；

5、s2、对采集获得的静息态数据进行预处理操作，其中预处理包括时间层校正、头动校正、标准化、平滑、滤波和去线性漂移；

6、s3、根据上述预处理后正常健康被试的静息态fmri脑影像数据，使用cstica方法搭建多目标优化框架，其中时间先验信息和空间先验信息通过主成分分析方法从被试所对应的时间独立成分和空间独立成分中获得；

7、s4、对于每个目标函数中矩阵的迭代，编码上采用两层循环遍历多个行或列向量，再进行向量两两相乘，组成含有多个未知变量的目标函数，从而实现矩阵的迭代优化；

8、s5、利用nsga-ii算法对上述经过处理的多目标函数模型进行求解，利用求解得到的分离矩阵还原相应感兴趣区的独立源信号。

9、在本发明中，预处理包括时间层校正、头动校正、标准化、平滑、滤波和去线性漂移等六个步骤。所有预处理步骤均在dparsf(http://rfmri.org/dparsf)工具上完成。属于现有技术，因此不作详细说明。

10、优选的，s3中的cstica方法包含如下步骤：

11、s31、假设已经获得了m维的混合信号向量x(i)，以及k维的相互独立的源信号s(j)，其中i＝1,2,...,m，j＝1,2,...,k，那么ica模型可表示为：

12、x＝as            (1)

13、其中x＝(x1,x2,…,xm)t表示m维观测信号矩阵，表示m×k维的混合矩阵，s＝(s1,s2,…,sk)t表示k维的源信号矩阵。由于只知道观测信号矩阵x，混合矩阵a和源信号矩阵s是未知数据，所以ica的目标就是通过不断地迭代优化，找到一个m×k维的分离矩阵w，将源信号矩阵s和混合矩阵a分离出来，为了更方便的表示所要求解的分离矩阵w，将上式(1)的矩阵x＝as进行可逆等式转换得到下式(2)：

14、y＝wx             (2)

15、其中y＝(y1,y2,…,ym)t可以看做是源信号s的一个近似估计，是m×k维的分离矩阵；

16、s32、在ica模型的基础上加入先验信息构成一个约束优化问题，其表达形式如式(3)所示：

17、

18、其中j(y)与式(2-3)类似，表示度量输出独立成分的目标函数，不等式g(y)≤0是引入感兴趣成分的先验信息，用来约束独立成分的输出，其中g(y)＝ε(y，r)-ξ，本文采用距离度量函数ε(y，r)＝e(y-r)2，用来度量输出信号y和包含先验信息的参考信号r，ξ表示预先给定的阈值参数；h(y)＝0是等式约束，将ica输出的独立分量单位化，从而保证算法在凸区域内求解，同时防止算法收敛到相同的独立分量；

19、s33、构建多目标cica模型框架，如式(4)：

20、

21、其中，j(wi)代表待估计成分的负熵，y＝zx代表x的白化过程，z表示白化矩阵。ti＝z-1wi表示si对应的时间过程。g(·)表示任意非二次函数，本文采用g(v)＝log(cosh(v))，v是一个零均值单位方差的随机高斯变量，即：

22、

23、其中，rsi表示空间先验信息，ε1(wi)＝e[sirsi]表示用皮尔森相关系数度量si和rsi之间的相似度。rti表示时间先验信息，ε2(wi)＝e[tirti]用来度量ti和rti之间的相似度；为方便运算及代码编写，将等式y＝zx，ti＝z-1wi及e[g(v)]的估计值0.374567代入式(2-7)，即：

24、

25、优选的，s4包含如下步骤：

26、s41、将式(6)中三个目标函数中与未知变量wi相邻进行运算的矩阵拆分成行或列的向量，即在编码过程中对两个运算的矩阵进行两层循环遍历，对未知向量wi进行降维处理；

27、s42、对s41循环遍历出来的向量进行内积运算，最后式(6)会转换成三个一维多变量的目标函数。

28、优选的，s5包含如下步骤：

29、s51、经过步骤4的矩阵处理，利用遗传算法将wi中的i个未知解进行个体编码，每个种群中的个体对应所求解分离向量中的每一个解，以此遵循自然进化准则。设置种群大小以及进化过程中个体选择、交叉和变异的概率，利用非支配排序遗传算法进行种群的迭代；

30、s52、种群迭代完成后，将每一组i个最优解集整合成分离向量wi，利用式和ti＝z-1wi求出空间源信号和时间源信号；

31、综上所述，本发明旨在利用fmri神经成像技术获取脑影像大数据情况下，通过充分利用各个脑区之间的功能连接状态，采用基于多目标时空先验信息的独立成分分析(independent component analysis based on constrained spatiotemporal,cstica)和带精英策略的快速非支配排序遗传算法(non-dominated sorted genetic algorithm-ⅱ)等方法来提取感兴趣区的独立成分。从而为脑皮质功能区病灶的定位(如脑肿瘤治疗、癫痫)、脑部疾病治疗(如：阿尔茨海默病、脑卒中、多发性硬化及帕金森综合症等)以及探索职业脑可塑性重组特征(如海员、飞行员)等提供坚实基础。

32、与现有技术相比，本发明有益效果及显著进步在于：本发明引入非支配排序遗传算法进行多目标约束独立成分模型的求解，有助于自适应且全局最优迭代出研究所需的独立源信号；本发明通过cstica方法和多目标优化框架，避免了ica方法出现成分不确定、检测到的目标成分准确率不理想等问题，提高感兴趣区独立源信号提取的精度；本发明通过遗传算法迭代出的独立源信号可以为脑皮质功能区病灶的定位、脑部疾病治疗以及探索职业脑可塑性重组特征等夯实基础。