基于密集电极阵列的多目标脑区精确电刺激方法与流程

【技术领域】

本发明涉及一种电刺激方法，具体涉及一种基于密集电极阵列的多目标脑区精确电刺激方法。

背景技术：

经颅直流电刺激(tdcs)是一种利用恒定微小直流电非侵入调控大脑神经元活动的技术，目前对各种脑疾病(脑卒中、癫痫、帕金森病、抑郁症、耳鸣、阿尔茨海默病等)的治疗具有重要的支持作用。其是利用放置在头皮上的电极来注入电流，使神经元静息电位呈现去极化或超极化效应。这种刺激形成的电场很大程度受到个体头部解剖结构和生物物理学的影响，难以在目标脑感兴趣区域(roi)中获得所需的电流密度场。若采用密集电极阵列对整个头皮全方位覆盖，通过一定的优化策略计算出各电极应注入的电流强度，注入不同大小的电流可以将刺激得到的整体的刺激电场分布限制在roi中。而大脑功能的实现通常会涉及多个脑区的交互，同时刺激多个不同的脑区可能能够更好地影响特定的大脑功能，鉴于此dmochowski提出一种基于最小二乘的优化策略，该策略通过头皮各电极记录的eeg信号给与生成这些eeg波形的神经元一定程度的电刺激，但此模式无法指定目标脑区的位置。

技术实现要素：

针对上述方法无法指定目标脑曲位置的问题，本发明提供一种基于密集电极阵列的多目标脑区精确电刺激方法。

本发明是通过以下技术方案实现的，提供一种基于密集电极阵列的多目标脑区精确电刺激方法，包括以下步骤：

s1采集t1加权mri头部影像数据；

s2使用fsl软件对采集的mri头部影响数据进行组织分割，分割成头皮、头骨、脑脊液、脑灰质、脑白质五部分；

s3将s2分割的五部分组织进行人工修正；

s4将n个电极紧贴于头皮表面，取各个电极与头皮接触端向上1mm作为导电胶得到带有电极的头模型；

s5对带有电极的头模型进行有限元三角面片化处理；

s6使用comsol计算出各电极注入单位电流强度时大脑有限元节点处的电场分布矩阵；

s7计算出优化策略。

特别的，所述s4中电极个数n＝64个，所述电极为直径12mm，深度3mm的圆柱形。

特别的，所述s7中优化策略按具体如下方法实施：

s71构建目标函数：

于公式(1)中，i为电极电流强度矩阵，e(r)为节点r处的电场强度，d(r)为节点r处的电流密度的期望方向；

s72简化目标函数：

于公式(2)中，wi为第i个脑区的系数矩阵，所述wi按如下公式进行计算：

于公式(3)中，vm为第m个元素的体积，代表第m个元素内的电流密度期望方向的转置，所述am按如下公式进行计算：

am＝[a1(rm)，a2(rm)，…，an-1(rm)](4)，

于公式(4)中，a为单位电场分布矩阵，n为电极个数，rm表示第m个元素节点位置。

s73施加刺激前，根据各脑区尺寸对其对应的系数矩阵wi进行加权，按如下公式进行加权：

于公式(5)中，m为目标脑感兴趣区域个数，所述αi按如下公式进行计算：

于公式(7)中，其中i＝1，2…m；j＝1，2…m。

特别的，所述s7中，对一部分目标脑区施加高强度刺激、另一部分施加低强度刺激时，对公式(6)引入第i个目标脑感兴趣区域所需要刺激强度的加权系数si，得到如下公式：

特别的，所述s7中，为保证电流输入的安全性，需引入安全限制条件，具体如下：

于公式(8)中，σ为电导率，e为电场分布矩阵，pmax表示安全限制下的功率最大值，若满足上式条件，表示电流安全，表示目标脑感兴趣区域外剩余脑区接受的刺激总功率小于最大允许功率，对进行数值逼近，得到如下公式：

当需要给每个目标脑区施加相同的刺激强度时，按如下公式约束结果精度：

std([w1i，w2i，…，wmi])＜s0(12)。

特别的，所述s7中，当需要给一部分目标脑区施加高强度刺激、另一部分施加低强度刺激时，加入电流分布约束条件，约束条件如下：

本发明提供一种基于密集电极阵列的多目标脑区精确电刺激方法，其可仿真实现多目标脑区精确的tdcs的多种刺激模式，可将电场分布很好的限制在roi中，无需反复挪动电极位置，具有灵活性强，刺激电流安全、精确的优点。

【附图说明】

图1为本发明基于密集电极阵列的多目标脑区精确电刺激方法的流程图；

图2为采用本发明基于密集电极阵列的多目标脑区精确电刺激方法得到的带有电极的头模型；

图3为采用本发明基于密集电极阵列的多目标脑区精确电刺激方法脑部得到强度等级相似刺激的仿真效果图。

【具体实施方式】

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进一步详细说明。

请参阅图1，本发明提供一种基于密集电极阵列的多目标脑区精确电刺激方法，包括以下步骤：

s1采集t1加权mri头部影像数据；

s2使用fsl软件对采集的mri头部影响数据进行组织分割，分割成头皮、头骨、脑脊液、脑灰质、脑白质五部分；

s3使用simpleware软件将s2分割的五部分组织进行人工修正；

s4将64个直径12mm，深度3mm的圆柱形的电极紧贴于头皮表面，取各个电极与头皮接触端向上1mm作为导电胶得到带有电极的头模型，如图2；

s5使用simpleware软件对带有电极的头模型进行有限元三角面片化处理；

s6使用comsol计算出各电极注入单位电流强度时大脑有限元节点处的电场分布矩阵；

s7基于matlab软件使用interior-point搜索算法计算出优化策略。

特别的，所述s7中优化策略按具体如下方法实施：

s71构建目标函数：

于公式(1)中，i为电极电流强度矩阵，e(r)为节点r处的电场强度，d(r)为节点r处的电流密度的期望方向；

s72简化目标函数：

于公式(2)中，wi为第i个脑区的系数矩阵，所述wi按如下公式进行计算：

于公式(3)中，vm为第m个元素的体积，代表第m个元素内的电流密度期望方向的转置，所述am按如下公式进行计算：

am＝[a1(rm)，a2(rm)，…，an-1(rm)](4)，

于公式(4)中，a为单位电场分布矩阵，n为电极个数，rm表示第m个元素节点位置；

s73施加刺激前，根据各脑区尺寸对其对应的系数矩阵wi进行加权，按如下公式进行加权：

于公式(5)中，m为目标脑感兴趣区域个数，所述αi按如下公式进行计算：

于公式(7)中，其中i＝1，2…m；j＝1，2…m。

特别的，所述s7中，对一部分目标脑区施加高强度刺激、另一部分施加低强度刺激时，对公式(6)引入第i个目标脑感兴趣区域所需要刺激强度的加权系数si，得到如下公式：

特别的，所述s7中，为保证电流输入的安全性，需引入安全限制条件，具体如下：

于公式(8)中，σ为电导率，e为电场分布矩阵，pmax表示安全限制下的功率最大值。若满足上式条件，表示电流安全，表示目标脑感兴趣区域外剩余脑区接受的刺激总功率小于最大允许功率，对进行数值逼近，得到如下公式：

当需要给每个目标脑区施加相同的刺激强度时，按如下公式约束结果精度：

std([w1i，w2i，…，wmi])＜s0(12)。

特别的，所述s7中，当需要给一部分目标脑区施加高强度刺激、另一部分施加低强度刺激时，加入电流分布约束条件，约束条件如下：

综上所述，当需要给每个roi施加相同或相似强度的刺激时，基于密集电极阵列的多目标脑区tdsc整体的优化策略①为：

优化策略①可以将各个目标脑区分配到的刺激强度之间的差异控制在一个很小的合理的范围内。

当需要给一部分目标脑区较高强度的刺激，另一部分目标脑区较低强度的刺激时，基于密集电极阵列的多目标脑区tdsc整体的优化策略②为：

以下通过一实施例对本发明进行解释说明，选用基于aal模板的脑区为例，选取三个roi——左侧额上回(sfg)，右侧中央前回(precg)和左侧枕中回(mog)，如表1所示，仿真实现多目标脑区精确tdcs优化方法。

表1

当使用优化策略①时，计算出三个脑区的均衡化系数αi，标准化后得到权值为1：1.1595：1.4411，仿真优化得到结果，如图3所示，经过图3第一列可知三个roi得到了强度等级相似的刺激；在此基础上，给αi分别乘以[0.6,0.3,0.3]，[0.3,0.6,0.3]，[0.3,0.3,0.6]时，使用优化策略②进行优化，sfg脑区，precg脑区以及mog脑区中其中一个脑区将得到较低强度等级的刺激，其余两个脑区将会得到较高强度等级的刺激，如图3第二列至第四列所示。