CN113289249B - 基于密集电极阵列的多目标脑区精确电刺激系统 - Google Patents
基于密集电极阵列的多目标脑区精确电刺激系统 Download PDFInfo
- Publication number
- CN113289249B CN113289249B CN202110607521.3A CN202110607521A CN113289249B CN 113289249 B CN113289249 B CN 113289249B CN 202110607521 A CN202110607521 A CN 202110607521A CN 113289249 B CN113289249 B CN 113289249B
- Authority
- CN
- China
- Prior art keywords
- brain
- formula
- target brain
- stimulation
- intensity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/18—Applying electric currents by contact electrodes
- A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
- A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
- A61N1/36014—External stimulators, e.g. with patch electrodes
- A61N1/36025—External stimulators, e.g. with patch electrodes for treating a mental or cerebral condition
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/18—Applying electric currents by contact electrodes
- A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
- A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
- A61N1/36014—External stimulators, e.g. with patch electrodes
- A61N1/3603—Control systems
- A61N1/36031—Control systems using physiological parameters for adjustment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Geometry (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Hospice & Palliative Care (AREA)
- Psychology (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Social Psychology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Neurology (AREA)
- General Engineering & Computer Science (AREA)
- Developmental Disabilities (AREA)
- Child & Adolescent Psychology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Electrotherapy Devices (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于密集电极阵列的多目标脑区精确电刺激系统,包括以下步骤:采集T1加权MRI头部影像数据;使用FSL软件对采集的MRI头部影响数据进行组织分割,分割成头皮、头骨、脑脊液、脑灰质、脑白质五部分;将分割的五部分组织进行人工修正;将N个电极紧贴于头皮表面,取各个电极与头皮接触端向上1mm作为导电胶得到带有电极的头模型;对带有电极的头模型进行有限元三角面片化处理;计算出各电极注入单位电流强度时大脑有限元节点处的电场分布矩阵;计算出优化策略。该多目标脑区精确电刺激方法无需反复挪动电电极位置,即可实现多目标脑区精确tDCS的多种刺激,灵活性较强。
Description
【技术领域】
本发明涉及一种电刺激方法,具体涉及一种基于密集电极阵列的多目标脑区精确电刺激系统。
【背景技术】
经颅直流电刺激(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按如下公式进行计算:
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中,为保证电流输入的安全性,需引入安全限制条件,具体如下:
于公式(9)中,σ为电导率,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按如下公式进行计算:
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中,为保证电流输入的安全性,需引入安全限制条件,具体如下:
于公式(9)中,σ为电导率,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第二列至第四列所示。
Claims (2)
1.基于密集电极阵列的多目标脑区精确电刺激系统,其特征在于,包括以下步骤:
S1采集T1加权MRI头部影像数据;
S2使用FSL软件对采集的MRI头部影响数据进行组织分割,分割成头皮、头骨、脑脊液、脑灰质、脑白质五部分;
S3将S2分割的五部分组织进行人工修正;
S4将N个电极紧贴于头皮表面,取各个电极与头皮接触端向上1mm作为导电胶得到带有电极的头模型;
S5对带有电极的头模型进行有限元三角面片化处理;
S6使用COMSOL计算出各电极注入单位电流强度时大脑有限元节点处的电场分布矩阵;
S7计算出优化策略,其按具体如下方法实施:
S71构建目标函数:
于公式(1)中,I为电极电流强度矩阵,E(r)为节点r处的电场强度,d(r)为节点r处的电流密度的期望方向;
S72简化目标函数:
于公式(2)中,Wi为第i个脑区的系数矩阵,所述Wi按如下公式进行计算:
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;
对一部分目标脑区施加高强度刺激、另一部分施加低强度刺激时,对公式(6)引入第i个目标脑感兴趣区域所需要刺激强度的加权系数Si,得到如下公式:
为保证电流输入的安全性,需引入安全限制条件,具体如下:
于公式(9)中,σ为电导率,E为电场分布矩阵,pmax表示安全限制下的功率最大值,若满足上式条件,表示电流安全,表示目标脑感兴趣区域外剩余脑区接受的刺激总功率小于最大允许功率,对进行数值逼近,得到如下公式:
当需要给每个目标脑区施加相同的刺激强度时,按如下公式约束结果精度:
std([W1I,W2I,…,WMI])<s0 (12),
当需要给一部分目标脑区施加高强度刺激、另一部分施加低强度刺激时,加入电流分布约束条件,约束条件如下:
2.根据权利要求1所述的基于密集电极阵列的多目标脑区精确电刺激系统,其特征在于,所述S4中电极个数N=64个,所述电极为直径12mm,深度3mm的圆柱形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110607521.3A CN113289249B (zh) | 2021-06-01 | 2021-06-01 | 基于密集电极阵列的多目标脑区精确电刺激系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110607521.3A CN113289249B (zh) | 2021-06-01 | 2021-06-01 | 基于密集电极阵列的多目标脑区精确电刺激系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113289249A CN113289249A (zh) | 2021-08-24 |
CN113289249B true CN113289249B (zh) | 2022-08-09 |
Family
ID=77326629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110607521.3A Active CN113289249B (zh) | 2021-06-01 | 2021-06-01 | 基于密集电极阵列的多目标脑区精确电刺激系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113289249B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113288176B (zh) * | 2021-06-01 | 2023-09-05 | 科悦医疗(苏州)有限公司 | 区分慢波类型的在线慢波增强系统及方法 |
CN114870249B (zh) * | 2022-04-18 | 2023-06-13 | 北京理工大学 | 一种闭环自适应交流电刺激神经网络调控方法及系统 |
CN117594193A (zh) * | 2024-01-17 | 2024-02-23 | 西安电子科技大学 | 基于深度学习的经颅直流电个体化刺激靶点定位方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110507904A (zh) * | 2019-08-22 | 2019-11-29 | 西安八水健康科技有限公司 | 一种计算密集阵列tDCS定向调制最优电极刺激模式的方法 |
CN110559554A (zh) * | 2019-08-22 | 2019-12-13 | 西安八水健康科技有限公司 | 一种融合经颅电刺激和迷走神经刺激的电刺激装置 |
CN110732082A (zh) * | 2019-09-27 | 2020-01-31 | 杭州电子科技大学 | 一种经颅直流电刺激和功能性电刺激的运动功能康复方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9694178B2 (en) * | 2013-10-21 | 2017-07-04 | Neuroelectronics Barcelona S.L. | Method and a system for optimizing the configuration of multisite transcranial current stimulation and a computer-readable medium |
US20160228702A1 (en) * | 2009-04-13 | 2016-08-11 | Research Foundation Of The City University Of New York | Neurocranial Electrostimulation Models, Systems, Devices and Methods |
US20160361534A9 (en) * | 2012-11-16 | 2016-12-15 | Rio Grande Neurosciences, Inc. | Variably configurable, adaptable electrode arrays and effectuating software, methods, and systems |
US10543367B2 (en) * | 2015-06-18 | 2020-01-28 | Taipei Medical University | Transcranial burst electrostimulation apparatus and its applications |
KR102063926B1 (ko) * | 2016-09-23 | 2020-01-08 | 기초과학연구원 | 뇌 자극 방법 |
CN112604163A (zh) * | 2020-12-30 | 2021-04-06 | 杭州电子科技大学 | 一种基于经颅直流电刺激的辅助记忆系统 |
-
2021
- 2021-06-01 CN CN202110607521.3A patent/CN113289249B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110507904A (zh) * | 2019-08-22 | 2019-11-29 | 西安八水健康科技有限公司 | 一种计算密集阵列tDCS定向调制最优电极刺激模式的方法 |
CN110559554A (zh) * | 2019-08-22 | 2019-12-13 | 西安八水健康科技有限公司 | 一种融合经颅电刺激和迷走神经刺激的电刺激装置 |
CN110732082A (zh) * | 2019-09-27 | 2020-01-31 | 杭州电子科技大学 | 一种经颅直流电刺激和功能性电刺激的运动功能康复方法 |
Non-Patent Citations (1)
Title |
---|
Optimized tDCS for Targeting Multiple Brain Regions: An Integrated Implementation;Yu Huang等;《2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC)》;20181029;第3545-3548页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113289249A (zh) | 2021-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113289249B (zh) | 基于密集电极阵列的多目标脑区精确电刺激系统 | |
Seibt et al. | The pursuit of DLPFC: non-neuronavigated methods to target the left dorsolateral pre-frontal cortex with symmetric bicephalic transcranial direct current stimulation (tDCS) | |
Wagner et al. | Noninvasive human brain stimulation | |
Dmochowski et al. | Optimal use of EEG recordings to target active brain areas with transcranial electrical stimulation | |
US8538543B2 (en) | System and method to design structure for delivering electrical energy to tissue | |
Butson et al. | Role of electrode design on the volume of tissue activated during deep brain stimulation | |
Merlet et al. | From oscillatory transcranial current stimulation to scalp EEG changes: a biophysical and physiological modeling study | |
Seo et al. | Relation between the electric field and activation of cortical neurons in transcranial electrical stimulation | |
Nummenmaa et al. | Targeting of white matter tracts with transcranial magnetic stimulation | |
US20160055304A1 (en) | Targeted electrical stimulation | |
Gomez-Tames et al. | Significant group-level hotspots found in deep brain regions during transcranial direct current stimulation (tDCS): A computational analysis of electric fields | |
CN106823137B (zh) | 一种优化神经调控的方法和装置 | |
Bangera et al. | Experimental validation of the influence of white matter anisotropy on the intracranial EEG forward solution | |
Xiao et al. | Theoretical optimization of stimulation strategies for a directionally segmented deep brain stimulation electrode array | |
AU2017281934B2 (en) | Systems and methods for visual analytics of clinical effects | |
Bai et al. | A review of computational models of transcranial electrical stimulation | |
Gomez-Tames et al. | Review on biophysical modelling and simulation studies for transcranial magnetic stimulation | |
Seo et al. | Effect of anatomically realistic full-head model on activation of cortical neurons in subdural cortical stimulation—a computational study | |
Sepulveda et al. | Finite element analysis of current pathways with implanted electrodes | |
CN115779268A (zh) | 一种经颅电刺激的多脑区电极优化方法 | |
Seo et al. | Computational study of subdural cortical stimulation: effects of simulating anisotropic conductivity on activation of cortical neurons | |
Fernandes et al. | Modelling studies of non-invasive electric and magnetic stimulation of the spinal cord | |
Frijns et al. | Integrated use of volume conduction and neural models to simulate the response to cochlear implants | |
Salvador et al. | Personalization of multi-electrode setups in tCS/tES: Methods and advantages | |
Rezaee et al. | A computational pipeline to find lobule-specific electric field distribution during non-invasive cerebellar stimulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |