CN104740781A - 基于神经纤维走向的矢量经颅磁刺激方法 - Google Patents
基于神经纤维走向的矢量经颅磁刺激方法 Download PDFInfo
- Publication number
- CN104740781A CN104740781A CN201510172215.6A CN201510172215A CN104740781A CN 104740781 A CN104740781 A CN 104740781A CN 201510172215 A CN201510172215 A CN 201510172215A CN 104740781 A CN104740781 A CN 104740781A
- Authority
- CN
- China
- Prior art keywords
- coil
- vector
- target
- stimulation
- nerve fiber
- 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.)
- Pending
Links
- 239000013598 vector Substances 0.000 title claims abstract description 138
- 210000004126 nerve fiber Anatomy 0.000 title claims abstract description 123
- 238000011491 transcranial magnetic stimulation Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000000638 stimulation Effects 0.000 claims abstract description 90
- 239000000835 fiber Substances 0.000 claims abstract description 57
- 210000003792 cranial nerve Anatomy 0.000 claims abstract description 50
- 230000005684 electric field Effects 0.000 claims abstract description 21
- 210000004556 brain Anatomy 0.000 claims description 45
- 210000004761 scalp Anatomy 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 10
- 210000003710 cerebral cortex Anatomy 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 238000002598 diffusion tensor imaging Methods 0.000 claims description 6
- 230000004936 stimulating effect Effects 0.000 claims description 6
- 239000003990 capacitor Substances 0.000 claims description 4
- 238000009792 diffusion process Methods 0.000 claims description 4
- 230000005674 electromagnetic induction Effects 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 210000005036 nerve Anatomy 0.000 abstract description 9
- 230000007383 nerve stimulation Effects 0.000 abstract description 3
- 230000002490 cerebral effect Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 210000003050 axon Anatomy 0.000 description 2
- 230000007830 nerve conduction Effects 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 208000021384 Obsessive-Compulsive disease Diseases 0.000 description 1
- 208000018737 Parkinson disease Diseases 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 230000001054 cortical effect Effects 0.000 description 1
- 210000001787 dendrite Anatomy 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 206010015037 epilepsy Diseases 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 230000028161 membrane depolarization Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 208000028173 post-traumatic stress disease Diseases 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Landscapes
- Magnetic Treatment Devices (AREA)
Abstract
一种基于神经纤维走向的矢量经颅磁刺激方法,包括:确定目标脑神经纤维的刺激方向;将线圈产生的刺激最强的点与目标刺激点进行位置叠加,同时将线圈产生的磁场最强的方向与目标刺激方向相重合;分别计算基于目标神经纤维走向的脑神经纤维各空间体素点感应电场矢量,以及感应电流矢量。本发明能够提高神经刺激的空间分辨率,使用者可以根据实验要求,定位目标刺激神经,通过线圈与神经位置匹配,达到可预期的刺激效果。本发明还可以估算基于线圈与神经纤维位置关系的感应电场矢量及感应电流矢量,使用者根据线圈与大脑刺激区域的位置关系,计算出基于目标神经纤维走向的感应电场矢量和感应电流矢量。
Description
技术领域
本发明涉及一种经颅磁刺激方法。特别是涉及一种基于神经纤维走向的矢量经颅磁刺激方法。
背景技术
经颅磁刺激(TMS)是近30年发展起来的一种新的刺激技术,由于磁刺激具有无创,操作简便,相比于电刺激技术具有电安全性好等优势,在神经、精神、信息等领域,尤其是抑郁症,强迫症,帕金森症,癫痫以及创伤后应激障碍等的治疗方面,具有很好的应用前景,成为近年来研究的热点。
经颅磁刺激装置是通过大功率电容向电感线圈放电,短时间的在线圈上产生变化的电流。根据电磁感应定律,变化的电流会在线圈周围产生变化的磁场。
将这种变化的磁场作用于大脑,透过头骨、脑脊液、大脑皮层等组织,到达到神经纤维处。变化的磁场在神经纤维上会感生出感应电场。当感应电场达到一定数值时,会改变神经纤维内外的电位差异,称之为去极化作用。这种作用在神经纤维上进行传递,并通过神经树突在神经之间进行传递,从而达到改变神经传导的刺激效果。
多年来,经颅磁刺激仪器在功能和刺激效果上都有了很大发展。现有的技术是基于大脑的核磁共振图像进行刺激。通过扫描人脑的核磁共振图像,重建大脑皮质的三维图像,使用定位导航技术确定要刺激的脑区进行刺激。但其存在如下不足:
1.对于刺激区域的分辨率太低,神经刺激主要是针对神经轴突的刺激,通过改变轴突内外电位来达到神经传导的目的,而定位与皮质层面的刺激器其刺激位置分辨率太低。
2.刺激效果难以精确评估。线圈放电产生的磁场会在大脑刺激区域的神经纤维上产生感应电场,感应电场强度直接关系到刺激效果。感应电场强度的计算与线圈放电电流,线圈与大脑刺激区域位置关系有关,现有技术很难实现基于神经纤维走向的精确磁刺激。
发明内容
本发明所要解决的技术问题是,提供一种能够提高神经刺激的空间分辨率的基于神经纤维走向的矢量经颅磁刺激方法。
本发明所采用的技术方案是:一种基于神经纤维走向的矢量经颅磁刺激方法,包括如下步骤:
1)确定目标脑神经纤维的刺激方向;
2)将线圈产生的刺激最强的点与目标刺激点进行位置叠加,同时将线圈产生的磁场最强的方向与目标刺激方向相重合;
3)分别计算基于目标神经纤维走向的脑神经纤维各空间体素点感应电场矢量,以及感应电流矢量。
步骤1)所述的脑神经纤维的刺激方向,首先获取人脑扩散张量成像数据,根据纤维示踪技术,得到目标神经纤维;在目标神经纤维上提取一定数量的空间体素点,并对所述的空间体素点进行线性拟合,得到初次拟合矢量方向;逐次增加相同数量的像素点再次进行拟合,每拟合一次都要计算拟合后的矢量方向与初始拟合矢量方向之间的夹角,当夹角角度大于设定值时,拟合终止,取终止前一次的矢量方向作为最佳拟合方向,确定为目标神经纤维的刺激方向。
所述的线性拟合,是首先在目标神经纤维上确定目标刺激点,并以目标刺激点为起点,沿目标神经纤维,以大脑皮层后联合至前联合方向或者前联合至后联合方向,取若干空间体素点,并对由空间体素点组成的纤维采用空间求回归直线法进行线性拟合,得到拟合方向。
所述的计算拟合后的矢量方向与初次拟合矢量方向之间的夹角,是指首先在笛卡尔坐标系下的图像中,分别提取两个矢量上任意一点坐标及目标刺激点坐标,然后依据坐标点分别计算出两个矢量的单位矢量,最后计算单位矢量夹角作为两个矢量的夹角。
步骤2)包括:
(1)线圈采用八字线圈,使用三维定位导航技术对线圈形状和人脑轮廓进行定位,使用笛卡尔坐标***一线圈轮廓及人脑头皮、目标脑神经纤维上各点坐标;
(2)在笛卡尔坐标系下的目标脑神经纤维上各个空间体素坐标数据集中,选定目标刺激点坐标a1(x1,y1,z1),计算得到目标刺激点投影到人脑头皮上点的坐标a'1(x'1,y'1,z'1),并将线圈产生磁场最强的点与a'1进行位置叠加,确保目标刺激点的刺激磁场最强,所述最强的点是八字线圈中心位置点;
(3)做基于八字线圈平面的三个正交方向x、y、z,根据线圈形状及电磁感应原理,线圈y方向产生的磁场最强,将线圈y方向与目标脑神经纤维的刺激方向相重合,确保线圈沿目标刺激方向进行刺激。
步骤(3)所述的将线圈y方向与目标脑神经纤维的刺激方向相重合,包括:
基于八字线圈平面的三个正交方向x、y、z存在三个正交面xy、yz、xz,设A为目标神经纤维的刺激方向,A到正交面xy的投影为到正交面yz的投影为到正交面xz的投影为将线圈y方向与A之间的夹角用γ表示,线圈y方向与之间的夹角用α表示,线圈y方向与的夹角用β表示,线圈y方向与的夹角用θ表示;
线圈绕x轴方向的转动,调整夹角β的角度;线圈绕y轴方向的转动,调整夹角θ的角度;线圈绕z轴方向的转动,调整夹角α的角度,通过线圈的调整,将其中任意两个夹角的角度调整为0°,然后再调节第三个夹角的角度,来调节线圈y方向与A之间的夹角γ,此时第三个夹角的角度与夹角γ的角度一一对应,当夹角γ角度为0°时,线圈y方向与A方向相同。
步骤3)包括:
(1)在确保步骤2)中的第(2)步和第(3)步实现的条件下,提取目标神经纤维上各空间体素坐标a1(x1,y1,z1),a2(x2,y2,z2)…an(xn,yn,zn),设线圈周长为m,则提取线圈单位长度矢量的坐标d1(xd1,yd1,zd1),d2(xd2,yd2,zd2)…dm(xdm,ydm,zdm),计算目标神经纤维任一个空间体素点到线圈单位长度矢量的距离r为
(2)将放电电路相关参数电阻R、线圈电感L、充电电容容值C、充电电压U代入RLC电路的放电电流公式:计算得到电流变化率最大值
其中,
(3)将和距离r代入感应电场值公式计算得出基于目标刺激方向的目标脑神经纤维上各个空间体素的感应电场矢量;
其中,为线圈电流变化率最大值,为线圈上的一个单位长度,μ0为真空磁导率,r为目标神经纤维任一个空间体素点到线圈单位长度矢量上的距离,L为线圈周长;
(4)根据扩散张量成像原理可知,目标神经纤维上各空间体素点都存在三个方向正交的扩散张量分量x、y、z,且其中一个方向的分量幅值最大,设为x方向,所述x方向与目标神经纤维走向一致,基于目标脑神经纤维走向的各空间体素点感应电场矢量是在x方向分量上的投影
(5)将投影代入感应电流公式计算得出基于目标脑神经纤维走向的各空间体素点感应电流矢量;
其中,σ为目标神经纤维的电导率。
所述的目标脑神经纤维,是单一脑神经纤维或者多束脑神经纤维。
所述的线圈,是单一线圈或者多个线圈的神经纤维的矢量磁刺激,当是多个线圈刺激时,多个线圈同时针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激。
所述的刺激,是指同步、顺序或者随机的脑神经纤维矢量磁刺激;
同步是指多个线圈同时针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激;
顺序是指多个线圈按一定顺序针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束并按一定的线圈顺序进行基于神经纤维走向的矢量经颅磁刺激;
随机是指多个线圈随机针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束随机进行基于神经纤维走向的矢量经颅磁刺激。
本发明的基于神经纤维走向的矢量经颅磁刺激方法,具有的优点和积极效果:
1.提高神经刺激的空间分辨率。本发明采用了基于神经纤维走向的刺激方式对于刺激效果更加准确且具有方向性,使用者可以根据实验要求,定位目标刺激神经,通过线圈与神经位置匹配,达到可预期的刺激效果。并且这样的刺激效果,可以实现对于神经层面的脑区关联性的研究。
2.设计还可以估算基于线圈与神经纤维位置关系的感应电场矢量及感应电流矢量,使用者根据线圈与大脑刺激区域的位置关系,计算出基于目标神经纤维走向的感应电场矢量和感应电流矢量。
附图说明
图1是本发明基于神经纤维走向的矢量经颅磁刺激方法流程示意图;
图2是本发明的神经纤维拟合示意图;
图3是本发明的神经纤维目标刺激点与线圈中心点匹配示意图;
图4是本发明的神经纤维目标刺激方向与线圈产生的磁场最强的方向夹角示意图;
图5是本发明的线圈放电电路示意图;
图6是基于目标脑神经纤维走向的感应电场与感应电流生成示意图。
具体实施方式
下面结合实施例和附图对本发明的基于神经纤维走向的矢量经颅磁刺激方法做出详细说明。
如图1所示,本发明的基于神经纤维走向的矢量经颅磁刺激方法,包括如下步骤:
1)确定目标脑神经纤维的刺激方向;
所述的脑神经纤维的刺激方向,首先获取人脑扩散张量成像(diffusion tensorimaging,DTI)数据,根据纤维示踪技术,得到目标神经纤维;在目标神经纤维上提取一定数量的空间体素点,并对所述的空间体素点进行线性拟合,得到初次拟合矢量方向;逐次增加相同数量的像素点再次进行拟合,每拟合一次都要计算拟合后的矢量方向与初始拟合矢量方向之间的夹角,当夹角角度大于设定值时,拟合终止,取终止前一次的矢量方向作为最佳拟合方向,确定为目标神经纤维的刺激方向。
所述的目标神经纤维,是指扫描受试者大脑的磁共振-扩散张量图像,在目标脑神经纤维束附近选取一个目标点作为种子点,利用纤维示踪技术(包括:基于向量场法、“解扩散方程”法和“随机走”法)得到目标脑神经纤维的三维图像。
所述的线性拟合,是首先在目标神经纤维上确定目标刺激点P(附图2所示),并以目标刺激点为起点,沿目标神经纤维,以大脑皮层后联合(Posterior Commissure)至前联合(AnteriorCommissure)方向或者前联合至后联合方向,取若干空间体素点,并对由空间体素点组成的纤维采用空间求回归直线法进行线性拟合,得到初始拟合矢量。拟合具体方法如下:
设提取的空间体素点的坐标数据集a1(x1,y1,z1),a2(x2,y2,z2)…an(xn,yn,zn),设拟合矢量方程为z=Ax+By+C,设定每个空间体素点的z轴坐标到拟合矢量的z轴坐标之间距离为d,则有
d1=z1-C-Ax1-By1……dn=zn-C-Axn-Byn。当取最小值时,为最佳拟合矢量。即:
根据式(1),可以得到参数A、B和C,最终确定拟合矢量方程z=Ax+By+C。
所述的计算拟合后的矢量方向与初次拟合矢量方向之间的夹角,是指首先在笛卡尔坐标系下的图像中,分别提取两个矢量上任意一点坐标及目标刺激点坐标,然后依据坐标点分别计算出两个矢量的单位矢量,最后计算单位矢量夹角作为两个矢量的夹角;(附图2所示)。
两条拟合矢量之间夹角的计算方法如下:
在两条拟合矢量和上,各提取任意一点坐标(xA1 yA1 zA1)和(xA2 yA2 zA2)和目标刺激点P坐标(xP yP zP),求两条拟合矢量的单位矢量:
为的单位矢量,向量坐标(xT yT zT)为:
为的单位矢量,向量坐标(xR yR zR)为:
求单位向量与的夹角为:
所述的最佳拟合方向,例如:(附图2所示)目标刺激点为P,夹角阈值为5°,提取20个像素点进行拟合,结果作为矢量第二次再提取20个像素点(总共40个像素点)进行拟合,结果作为矢量且与夹角θ(1)<5°,第三次再提取20个像素点(总共60个像素点)进行拟合,结果作为矢量且与夹角θ(2)<5°第四次提取20个像素点(总共80个像素点)进行拟合,结果作为矢量且与夹角θ(3)>5°,则选取作为最终的拟合矢量,其方向为最终拟合方向。
2)将线圈产生的刺激最强的点与目标刺激点进行位置叠加,同时将线圈产生的磁场最强的方向与目标刺激方向相重合;包括:
(1)线圈采用八字线圈,使用三维定位导航技术对线圈形状和人脑轮廓进行定位,使用笛卡尔坐标***一线圈轮廓及人脑头皮、目标脑神经纤维上各点坐标;
(2)在笛卡尔坐标系下的目标脑神经纤维上各个空间体素坐标数据集a1(x1,y1,z1),a2(x2,y2,z2)…an(xn,yn,zn)中,选定目标刺激点坐标a1(x1,y1,z1),计算得到目标刺激点投影到人脑头皮上点的坐标a'1(x'1,y'1,z'1),并将线圈产生磁场最强的点与a'1进行位置叠加,确保目标刺激点的刺激磁场最强,所述最强的点是八字线圈中心位置点(附图3所示);
所述的目标刺激点投影到人脑头皮上点的坐标a'1(x'1,y'1,z'1)的选取方法为:
根据头皮定位结果,可以得到头皮上各点的坐标a'1(x'1,y'1,z'1),a'2(x'2,y'2,z'2)…a'n(x'n,y'n,z'n)。计算目标刺激点a1(x1,y1,z1)到头皮各点距离,即:
距离最短的点作为目标刺激点在头皮上的投影点a'1(x'1,y'1,z'1);
(3)做基于八字线圈平面的三个正交方向x、y、z(附图4所示),根据线圈形状及电磁感应原理,线圈y方向产生的磁场最强,将线圈y方向与目标脑神经纤维的刺激方向相重合,确保线圈沿目标刺激方向进行刺激。
所述的将线圈y方向与目标脑神经纤维的刺激方向相重合,包括:
基于八字线圈平面的三个正交方向x、y、z(附图4所示)存在三个正交面xy、yz、xz,设A为目标神经纤维的刺激方向,A到正交面xy的投影为到正交面yz的投影为到正交面xz的投影为将线圈y方向与A之间的夹角用γ表示,线圈y方向与之间的夹角用α表示,线圈y方向与的夹角用β表示,线圈y方向与的夹角用θ表示;
线圈绕x轴方向的转动,调整夹角β的角度;线圈绕y轴方向的转动,调整夹角θ的角度;线圈绕z轴方向的转动,调整夹角α的角度,通过线圈的调整,将其中任意两个夹角的角度调整为0°,然后再调节第三个夹角的角度,来调节线圈y方向与A之间的夹角γ,此时第三个夹角的角度与夹角γ的角度一一对应,当夹角γ角度为0°时,线圈y方向与A方向相同。
3)分别计算基于目标神经纤维走向的脑神经纤维各空间体素点感应电场矢量,以及感应电流矢量。包括:
(1)在确保步骤2)中的第(2)步和第(3)步实现的条件下,提取目标神经纤维上各空间体素坐标a1(x1,y1,z1),a2(x2,y2,z2)…an(xn,yn,zn),设线圈周长为m,则提取线圈单位长度矢量的坐标d1(xd1,yd1,zd1),d2(xd2,yd2,zd2)…dm(xdm,ydm,zdm),计算目标神经纤维任一个空间体素点到线圈单位长度矢量的距离,即:
(2)(附图5所示)将放电电路相关参数电阻R、线圈电感L、充电电容容值C、充电电压U代入RLC电路的放电电流公式,即:
计算得到电流变化率最大值其中,
(3)将和距离r代入感应电场值公式,即:
计算得出基于目标刺激方向的目标脑神经纤维上各个空间体素的感应电场矢量;
其中,为线圈电流变化率最大值,为线圈上的一个单位长度,μ0为真空磁导率,r为目标神经纤维任一个空间体素点到线圈单位长度矢量上的距离,L为线圈周长;
(4)根据扩散张量成像原理可知,目标神经纤维上各空间体素点都存在三个方向正交的扩散张量分量x、y、z(附图6所示),且其中一个方向的分量幅值最大,设为x方向,所述x方向与目标神经纤维走向一致,基于目标脑神经纤维走向的各空间体素点感应电场矢量是在x方向分量上的投影(图6所示);
(5)将投影代入感应电流公式计算得出基于目标脑神经纤维走向的各空间体素点感应电流矢量;
其中,σ为目标神经纤维的电导率;
在本发明的一种基于神经纤维走向的矢量经颅磁刺激方法中:
所述的目标脑神经纤维,可以是单一脑神经纤维或者多束脑神经纤维;即在提取目标脑神经纤维的空间体素提取过程中,提取一束或者多束(相互存在交叉或者不存在交叉)的脑神经纤维束,进行基于神经纤维走向的矢量磁刺激。
所述的线圈,是单一线圈或者多个线圈的神经纤维的矢量磁刺激,当是多个线圈刺激时,多个线圈同时针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激。
所述的刺激,是指同步、顺序或者随机的脑神经纤维矢量磁刺激;其中:
同步是指多个线圈同时针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激;
顺序是指多个线圈按一定顺序针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束并按一定的线圈顺序进行基于神经纤维走向的矢量经颅磁刺激;
随机是指多个线圈随机针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束随机进行基于神经纤维走向的矢量经颅磁刺激。
Claims (10)
1.一种基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,包括如下步骤:
1)确定目标脑神经纤维的刺激方向;
2)将线圈产生的刺激最强的点与目标刺激点进行位置叠加,同时将线圈产生的磁场最强的方向与目标刺激方向相重合;
3)分别计算基于目标神经纤维走向的脑神经纤维各空间体素点感应电场矢量,以及感应电流矢量。
2.根据权利要求1所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,步骤1)所述的脑神经纤维的刺激方向,首先获取人脑扩散张量成像数据,根据纤维示踪技术,得到目标神经纤维;在目标神经纤维上提取一定数量的空间体素点,并对所述的空间体素点进行线性拟合,得到初次拟合矢量方向;逐次增加相同数量的像素点再次进行拟合,每拟合一次都要计算拟合后的矢量方向与初始拟合矢量方向之间的夹角,当夹角角度大于设定值时,拟合终止,取终止前一次的矢量方向作为最佳拟合方向,确定为目标神经纤维的刺激方向。
3.根据权利要求2所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,所述的线性拟合,是首先在目标神经纤维上确定目标刺激点,并以目标刺激点为起点,沿目标神经纤维,以大脑皮层后联合至前联合方向或者前联合至后联合方向,取若干空间体素点,并对由空间体素点组成的纤维采用空间求回归直线法进行线性拟合,得到拟合方向。
4.根据权利要求2所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,所述的计算拟合后的矢量方向与初次拟合矢量方向之间的夹角,是指首先在笛卡尔坐标系下的图像中,分别提取两个矢量上任意一点坐标及目标刺激点坐标,然后依据坐标点分别计算出两个矢量的单位矢量,最后计算单位矢量夹角作为两个矢量的夹角。
5.根据权利要求1所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,步骤2)包括:
(1)线圈采用八字线圈,使用三维定位导航技术对线圈形状和人脑轮廓进行定位,使用笛卡尔坐标***一线圈轮廓及人脑头皮、目标脑神经纤维上各点坐标;
(2)在笛卡尔坐标系下的目标脑神经纤维上各个空间体素坐标数据集中,选定目标刺激点坐标a1(x1,y1,z1),计算得到目标刺激点投影到人脑头皮上点的坐标a'1(x'1,y'1,z'1),并将线圈产生磁场最强的点与a'1进行位置叠加,确保目标刺激点的刺激磁场最强,所述最强的点是八字线圈中心位置点;
(3)做基于八字线圈平面的三个正交方向x、y、z,根据线圈形状及电磁感应原理,线圈y方向产生的磁场最强,将线圈y方向与目标脑神经纤维的刺激方向相重合,确保线圈沿目标刺激方向进行刺激。
6.根据权利要求5所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,步骤(3)所述的将线圈y方向与目标脑神经纤维的刺激方向相重合,包括:
基于八字线圈平面的三个正交方向x、y、z存在三个正交面xy、yz、xz,设A为目标神经纤维的刺激方向,A到正交面xy的投影为到正交面yz的投影为到正交面xz的投影为将线圈y方向与A之间的夹角用γ表示,线圈y方向与之间的夹角用α表示,线圈y方向与的夹角用β表示,线圈y方向与的夹角用θ表示;
线圈绕x轴方向的转动,调整夹角β的角度;线圈绕y轴方向的转动,调整夹角θ的角度;线圈绕z轴方向的转动,调整夹角α的角度,通过线圈的调整,将其中任意两个夹角的角度调整为0°,然后再调节第三个夹角的角度,来调节线圈y方向与A之间的夹角γ,此时第三个夹角的角度与夹角γ的角度一一对应,当夹角γ角度为0°时,线圈y方向与A方向相同。
7.根据权利要求1所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,步骤3)包括:
(1)在确保步骤2)中的第(2)步和第(3)步实现的条件下,提取目标神经纤维上各空间体素坐标a1(x1,y1,z1),a2(x2,y2,z2)…an(xn,yn,zn),设线圈周长为m,则提取线圈单位长度矢量的坐标d1(xd1,yd1,zd1),d2(xd2,yd2,zd2)…dm(xdm,ydm,zdm),计算目标神经纤维任一个空间体素点到线圈单位长度矢量的距离r为 i∈1~n,dj∈d1~dm;
(2)将放电电路相关参数电阻R、线圈电感L、充电电容容值C、充电电压U代入RLC电路的放电电流公式:计算得到电流变化率最大值
其中,
(3)将和距离r代入感应电场值公式计算得出基于目标刺激方向的目标脑神经纤维上各个空间体素的感应电场矢量;
其中,为线圈电流变化率最大值,为线圈上的一个单位长度,μ0为真空磁导率,r为目标神经纤维任一个空间体素点到线圈单位长度矢量上的距离,L为线圈周长;
(4)根据扩散张量成像原理可知,目标神经纤维上各空间体素点都存在三个方向正交的扩散张量分量x、y、z,且其中一个方向的分量幅值最大,设为x方向,所述x方向与目标神经纤维走向一致,基于目标脑神经纤维走向的各空间体素点感应电场矢量是在x方向分量上的投影
(5)将投影代入感应电流公式计算得出基于目标脑神经纤维走向的各空间体素点感应电流矢量;
其中,σ为目标神经纤维的电导率。
8.根据权利要求1所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,所述的目标脑神经纤维,是单一脑神经纤维或者多束脑神经纤维。
9.根据权利要求1所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,所述的线圈,是单一线圈或者多个线圈的神经纤维的矢量磁刺激,当是多个线圈刺激时,多个线圈同时针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激。
10.根据权利要求1所述的基于神经纤维走向的矢量经颅磁刺激方法,其特征在于,所述的刺激,是指同步、顺序或者随机的脑神经纤维矢量磁刺激;
同步是指多个线圈同时针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激;
顺序是指多个线圈按一定顺序针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束并按一定的线圈顺序进行基于神经纤维走向的矢量经颅磁刺激;
随机是指多个线圈随机针对一束脑神经纤维束进行基于神经纤维走向的矢量经颅磁刺激,或者是每个线圈对应一束脑神经纤维束随机进行基于神经纤维走向的矢量经颅磁刺激。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510172215.6A CN104740781A (zh) | 2015-04-10 | 2015-04-10 | 基于神经纤维走向的矢量经颅磁刺激方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510172215.6A CN104740781A (zh) | 2015-04-10 | 2015-04-10 | 基于神经纤维走向的矢量经颅磁刺激方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104740781A true CN104740781A (zh) | 2015-07-01 |
Family
ID=53581385
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510172215.6A Pending CN104740781A (zh) | 2015-04-10 | 2015-04-10 | 基于神经纤维走向的矢量经颅磁刺激方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104740781A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105617532A (zh) * | 2016-02-16 | 2016-06-01 | 深圳英智科技有限公司 | 经颅磁刺激定位方法及*** |
CN106296808A (zh) * | 2016-08-15 | 2017-01-04 | 浙江工业大学 | 一种交互式脑纤维选择与可视化方法 |
CN108804852A (zh) * | 2018-06-28 | 2018-11-13 | 东汉新能源汽车技术有限公司 | 纤维定向的预测方法及装置 |
CN109011157A (zh) * | 2018-07-03 | 2018-12-18 | 中国科学院电工研究所 | 经颅磁刺激线圈定位装置 |
CN109999348A (zh) * | 2019-03-20 | 2019-07-12 | 上海市精神卫生中心(上海市心理咨询培训中心) | 一种基于弥散张量成像实现经颅磁刺激深部脑区的方法 |
CN113440731A (zh) * | 2021-08-11 | 2021-09-28 | 成都理工大学 | 用于颅脑电磁刺激的不对称异型线圈及颅脑电磁刺激*** |
CN115282488A (zh) * | 2022-06-21 | 2022-11-04 | 中国科学院自动化研究所 | 一种磁刺激仪 |
-
2015
- 2015-04-10 CN CN201510172215.6A patent/CN104740781A/zh active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105617532A (zh) * | 2016-02-16 | 2016-06-01 | 深圳英智科技有限公司 | 经颅磁刺激定位方法及*** |
CN105617532B (zh) * | 2016-02-16 | 2018-05-15 | 深圳英智科技有限公司 | 经颅磁刺激定位方法及*** |
CN106296808A (zh) * | 2016-08-15 | 2017-01-04 | 浙江工业大学 | 一种交互式脑纤维选择与可视化方法 |
CN106296808B (zh) * | 2016-08-15 | 2019-04-09 | 浙江工业大学 | 一种交互式脑纤维选择与可视化方法 |
CN108804852A (zh) * | 2018-06-28 | 2018-11-13 | 东汉新能源汽车技术有限公司 | 纤维定向的预测方法及装置 |
CN109011157A (zh) * | 2018-07-03 | 2018-12-18 | 中国科学院电工研究所 | 经颅磁刺激线圈定位装置 |
CN109999348A (zh) * | 2019-03-20 | 2019-07-12 | 上海市精神卫生中心(上海市心理咨询培训中心) | 一种基于弥散张量成像实现经颅磁刺激深部脑区的方法 |
CN113440731A (zh) * | 2021-08-11 | 2021-09-28 | 成都理工大学 | 用于颅脑电磁刺激的不对称异型线圈及颅脑电磁刺激*** |
CN113440731B (zh) * | 2021-08-11 | 2022-02-15 | 成都理工大学 | 用于颅脑电磁刺激的不对称异型线圈及颅脑电磁刺激*** |
CN115282488A (zh) * | 2022-06-21 | 2022-11-04 | 中国科学院自动化研究所 | 一种磁刺激仪 |
CN115282488B (zh) * | 2022-06-21 | 2023-09-12 | 中国科学院自动化研究所 | 一种磁刺激仪 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104740781A (zh) | 基于神经纤维走向的矢量经颅磁刺激方法 | |
Laakso et al. | Effects of coil orientation on the electric field induced by TMS over the hand motor area | |
Miranda et al. | Modeling the current distribution during transcranial direct current stimulation | |
Nummenmaa et al. | Comparison of spherical and realistically shaped boundary element head models for transcranial magnetic stimulation navigation | |
EP3106204B1 (en) | Intracerebral current simulation method and device thereof, and transcranial magnetic stimulation system including intracerebral current simulation device | |
Suh et al. | Influence of anisotropic conductivity in the skull and white matter on transcranial direct current stimulation via an anatomically realistic finite element head model | |
Sack et al. | Optimizing functional accuracy of TMS in cognitive studies: a comparison of methods | |
Laakso et al. | Computational analysis shows why transcranial alternating current stimulation induces retinal phosphenes | |
CN103917271B (zh) | 确定经颅磁刺激对脑部的效果的方法及设备 | |
Krieg et al. | Computational and experimental analysis of TMS-induced electric field vectors critical to neuronal activation | |
US9883812B2 (en) | Enhanced multi-core beamformer algorithm for sensor array signal processing by combining data from magnetoencephalography | |
De Geeter et al. | Effective electric fields along realistic DTI-based neural trajectories for modelling the stimulation mechanisms of TMS | |
Lee et al. | Electric field model of transcranial electric stimulation in nonhuman primates: correspondence to individual motor threshold | |
Gomez et al. | Uncertainty quantification in transcranial magnetic stimulation via high-dimensional model representation | |
KR102577945B1 (ko) | 국소적 경두개 자기 자극 장치 및 방법 | |
Soldati et al. | Computational errors of the induced electric field in voxelized and tetrahedral anatomical head models exposed to spatially uniform and localized magnetic fields | |
CN105142721B (zh) | 利用结构有效性指标来规划导航经颅磁刺激 | |
Laakso et al. | Computational dosimetry of induced electric fields during realistic movements in the vicinity of a 3 T MRI scanner | |
CN112704486B (zh) | 基于电磁仿真计算的tms线圈位姿图谱生成方法 | |
Gomez–Tames et al. | TMS activation site estimation using multiscale realistic head models | |
Salvador et al. | How electrode montage affects transcranial direct current stimulation of the human motor cortex | |
Salvador et al. | High permeability cores to optimize the stimulation of deeply located brain regions using transcranial magnetic stimulation | |
CN116099128A (zh) | 一种基于功能性核磁共振实现经颅磁精准刺激*** | |
Lueckel et al. | Functional connectivity-and E-field-optimized TMS targeting: method development and concurrent TMS-fMRI validation | |
JP2022519520A (ja) | 10-20システム基盤の位置情報提供方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150701 |
|
RJ01 | Rejection of invention patent application after publication |