[发明专利]基于种群多样性控制的脑电源定位方法、系统及设备

说明书

本发明公开了一种基于种群多样性控制的脑电源定位方法、系统及设备，基于患者的个性化MRI，构建患者的个性化真实头模型与源模型。通过基准点法实现EEG‑MRI的粗略配准，通过对电极位置手动微调，实现EEG‑MRI的精细配准，降低空间模型产生的误差。基于配准后的空间信息，建立体积传导模型，计算正向传导矩阵，进而建立源定位数学模型，为降低计算复杂度，对源定位数学模型进行线性与非线性参数分离。面对逆问题的不适定性与病态性，提出一种基于种群多样性控制的量子粒子群算法(DQPSO)，通过控制算法的种群多样性，使算法迭代后期仍然能够保持较高的全局遍历寻优能力，进而提高算法收敛到全局最优的概率，保证高定位精度。

技术领域

本发明涉及生物医学技术领域，具体涉及一种基于种群多样性控制的脑电源定位方法、系统及设备。

背景技术

大脑神经活动时，伴随着电荷、电流的运动，外在表现形式为头皮表面采集到的脑电(EEG)。脑电源定位，就是通过头皮电位的测量数据反向推演脑内神经活动源分布，也即脑电逆问题。常用的方法之一是偶极子定位法，核心思想是假想脑内神经活动分布可以用有限个电流偶极子近似，一般为一到两个。源定位则是通过头皮电位推算出产生颅外电场的等效电流偶极子源的位置与强度信息。脑电逆问题一经提出便受到广泛关注，具备很高的科学价值及临床应用价值，有助于脑科学的研究及脑部疾病的诊断。

目前，临床非侵入性诊断脑部疾病通常需要结合病人的发作状态及多模态医学成像技术进行综合判断，一种是基于脑结构成像(如磁共振成像，MRI；计算机断层扫描，CT)的方法，该方法对形态学异常较为敏感，但对于脑功能异常，敏感性不高，而且脑功能异常往往早于脑结构异常。第二种是基于脑功能成像(如正电子发射体层摄影术，PET；功能磁共振成像，fMRI)的方法，以上两种方法时间分辨率较低，通常在1s左右，而且往往对人体存在一定损伤，如强磁场会对人体带来辐射作用，尤其对于孕妇和儿童。以上方法优点是空间分辨率较高，可达1-10mm左右。第三种是基于神经电生理学原理(EEG、MEG)的方法，该方法相比以上技术，时间分辨率高，通常在1ms左右，同时对人体完全无创无损，但空间分辨率不及以上技术。因此，通过脑电源定位技术，联合个性化的脑结构MRI提供的高空间分辨率和EEG提供的高时间分辨率，能够实现个性化高时空分辨率的病灶定位。

脑电源定位因其不适定性及病态性成为一个极具挑战性的研究课题。不适定性意味着解是不唯一的，无数颅内源分布对应着给定的颅外脑电分布，导致脑电逆问题属于多峰值问题。病态性意味着解空间的敏感性，微小的扰动都会带来解空间的剧烈震荡，而脑电在测量过程中易混入噪声，导致了解空间的复杂性。大多数非线性优化方法在面对多峰值问题时，缺乏跳出局部极值的能力，在面对解空间的震荡性时，缺乏空间遍历能力，容易出现早熟收敛，陷入局部极值，尤其在面对复杂的真实头模型时，定位精度较低。

发明内容

本发明为了解决上述技术问题，提出了如下技术方案：

第一方面，本发明提供了一种基于种群多样性控制的脑电源定位方法，所述方法包括以下步骤：

1)基于患者的个性化MRI，进行脑组织分割，通过三维重建获得患者的个性化头模型与源模型，通过坐标转换配准个性化头模与个性化源模；

2)建立EEG标准电极布置的几何模型，基于基准点法实现EEG标准电极模型与个性化MRI的自动粗略配准，再通过手动微调电极位置实现EEG标准电极模型与个性化MRI的精细配准；

3)构建体积传导空间模型，利用边界元法或有限元法计算正向传导矩阵；

4)建立源定位问题的数学模型并通过参数分离对源定位问题的数学模型进行降维；

5)假设电流偶极子源数目，初始化算法粒子数量和迭代次数，在解空间内均匀初始化粒子位置，利用基于种群多样性的量子粒子群算法：DQPSO进行全局搜索，引入种群多样性指标并设置种群多样性下限，根据种群多样性更新算法的收缩扩张因子，获得源定位结果。