A kind of H type coil optimization method for brain deep transcranial magnetic stimulation
Technical field
The present invention relates to a kind of H type coil for brain deep transcranial magnetic stimulation.Particularly relate to a kind of H type coil optimization method for brain deep transcranial magnetic stimulation based on true head conductivity model.
Background technology
Transcranial magnetic stimulation technology is a kind of painless, noinvasive, non-invasive technology, and the magnetic signal that the pulse current in coil produces can stimulate to cerebral nerve through skull undampedly.Along with research deeply, the development of technology, the stimulation of shallow-layer can not meet the demand of medical treatment and scientific research, and deep stimulation is arisen at the historic moment.Current coil can only stimulate the peripheral cortex region to 1.5-2cm.Deep stimulates while not increasing excitation, meets the stimulation to cerebral deep structure, reduces the damage to peripheral cortex.
H type coil has complicated structural design, and base portion stimulates major part, the tangential scalp in wire direction; Outstanding loop and link circuit, away from stimulation target, weaken the impact on target area.But current H type coil design is based on Sphere Model, and Sphere Model and true head structure exist larger difference.The accuracy of head model affects the electric field distribution law that transcranial magnetic stimulation coil produces on model.True head mould provides foundation for the design of H type coil, makes coil more level off to truth in the Electric Field Distribution that head produces.
The principle that H type coil utilizes induction field vector superposed, strengthens the induction field in the target tissue of deep; Base portion wire is certain intervals distribution, reduces the stimulus intensity of peripheral cortex with this.Conductor spacing and distribution influence electric field intensity and the deep characteristic of peripheral cortex.The complexity of head and H type coil illustrates that it is far from being enough for only studying the attenuation characteristic of electric field on a certain bar line segment, needs more detailed internal electric field distributed intelligence.Comprise Hippocampus, the isostructural limbic system of corpus amygdaloideum plays an important role in midbrain, information exchange between diencephalon and neopallium structure, be necessary limbic system to separate, set up the true head model comprising limbic system.The present invention adopts the method for finite element to calculate head Electric Field Distribution, and the H type coil conductor spacing of assessment prefrontal, top and sidepiece is on the impact of deep characteristic.
Summary of the invention
Technical problem to be solved by this invention is, a kind of conductor spacing by changing prefrontal, top and sidepiece being provided, using the distribution situation of method simulate electric field at true head conductivity model and the H type coil optimization method for brain deep transcranial magnetic stimulation of deep characteristic of finite element.
The technical solution adopted in the present invention is: a kind of H type coil optimization method for brain deep transcranial magnetic stimulation, use the method simulate electric field of finite element in the distribution situation of true head conductivity model and deep characteristic, optimize the conductor spacing of prefrontal, top and sidepiece, weaken head shallow-layer electric field intensity, increase head deep layer electric field intensity, comprise, following steps:
1) set up by nuclear magnetic resonance or CT data the true head structural model comprising limbic system, in finite element software, set up the true head conductivity model comprising limbic system;
2) the H type coil of tangential described true head structural model outer surface is set up according to the outline of true head structural model; Consider structure and the position characteristics of head edge system, determine prefrontal, top and sidepiece lead model distribution, be based upon the H type coil model at the prefrontal of true head structural model, top and sidepiece different conductor interval;
3) use the Electric Field Distribution of finite element method for simulating H type coil on true head conductivity model, the standard evaluating coil deep characteristic be set, comprise:
(1) in finite element software, pass through stress and strain model, true head conductivity model is divided into limited unit, set up the FEM (finite element) model that can apply to numerical computations; FEM (finite element) model is carried out loading solve, obtain the Electric Field Distribution of H type coil on true head model;
(2) in order to meet the object that deep stimulates, reduce the stimulation to peripheral cortex simultaneously, adopt the ratio of true head conductivity model surface field maximum and true head conductivity model limbic system electric field maximum as the standard evaluated H type coil and design;
4) to evaluate the standard of H type coil design for foundation, the best conductor spacing of H type coil prefrontal, top and sidepiece diverse location is determined.
Step 1) described in the foundation head model that comprises DEEP STRUCTURE comprise:
(1) by nuclear magnetic resonance or the interactively medical image control system of CT data importing, the operation comprising Threshold segmentation, region growth and manually modify, extract the edge of head and limbic system, set up the three-dimensional patch model of head and the three-dimensional patch model of limbic system respectively;
(2) be head physical model and limbic system physical model by three-dimensional for head patch model and the three-dimensional dough sheet model conversation of limbic system in finite element software, described head physical model and limbic system physical model are carried out overlap joint Boolean calculation, set up the true head structural model comprising limbic system, carry out electrical conductivity assignment, set up the true head conductivity model comprising limbic system.
Step 2) described in H type coil model, be adopt pulling mode to set up, the square of the conductor cross-section of H type coil model to be the length of side be 1mm, conductor spacing true head structural model surface 5mm.
Described pulling mode carries out described point line according to true head conductivity model outline in finite element software, form the path on tangential true head conductivity model top layer, set up a cross section perpendicular to path, pull described cross section along path, obtain the H type coil model on tangential true head conductivity model top layer.
Step 3) described in true head conductivity model surface field maximum less with the ratio of true head conductivity model limbic system electric field maximum, represent that H type coil depth characteristic is better.
A kind of H type coil optimization method for brain deep transcranial magnetic stimulation of the present invention, sets up the threedimensional model of head and limbic system by MRI or CT data.The construction features of CONSIDERING EDGE system, changes the conductor spacing of prefrontal, top and sidepiece, and assessment conductor spacing is on the impact of electric field deep characteristic.Due to the complexity of head construction, the wire of different interval produces Electric Field Distribution slightly difference, and therefore employing true head conductivity model surface field maximum and the ratio of limbic system maximum are as the standard evaluating coil deep characteristic.
Accompanying drawing explanation
Fig. 1 is the flow chart of the inventive method;
Fig. 2 is the head model and H type coil model that adopt method of the present invention to set up.
In figure
1:H type coil model 2: frontal lobe top wire
3: prefrontal wire 4: frontal lobe sidepiece wire
Detailed description of the invention
Below in conjunction with embodiment and accompanying drawing, a kind of H coil optimization method for brain deep transcranial magnetic stimulation of the present invention is described in detail.
H type coil is used for stimulating the DEEP STRUCTURE of brain, and cerebral limbic system is made up of limbic lobe and relevant cortex and subcortical structures, relevant with the higher nervous activity such as emotion, spirit, memory.Due to the difference of arrangement of conductors, limbic system produces different effects.In order to reach the darker stimulation degree of depth, reduce the electric field intensity of top layer scalp simultaneously.The H type coil conductor spacing of prefrontal and top and sidepiece is optimized.
As shown in Figure 1, a kind of H type coil optimization method for brain deep transcranial magnetic stimulation of the present invention, use the method simulate electric field of finite element in the distribution situation of true head conductivity model and deep characteristic, to optimize before frontal lobe, above and the conductor spacing of side, weaken head shallow-layer electric field intensity, increase head deep layer electric field intensity, comprise, following steps:
1) set up by nuclear magnetic resonance (MRI) or CT data the true head structural model comprising limbic system, in finite element software, set up the true head conductivity model comprising limbic system;
The head model that described foundation comprises DEEP STRUCTURE comprises:
(1) by nuclear magnetic resonance or the interactively medical image control system (Mimics) of CT data importing, the operation comprising Threshold segmentation, region growth and manually modify, extracts the edge of head and limbic system.After three-dimensional surface modeling, the model obtained is not through any optimization, on the basis not affecting head model structure, smooth operation is adopted to come level and smooth solid surface, simplify operation by dough sheet and reduce tri patch quantity, set up the three-dimensional patch model of head and the three-dimensional patch model of limbic system respectively;
(2) be head physical model and limbic system physical model by three-dimensional for head patch model and the three-dimensional dough sheet model conversation of limbic system in finite element software, described head physical model and limbic system physical model are carried out overlap joint Boolean calculation, ensure the node on two physical model Border faces and unit, set up the true head structural model comprising limbic system, carry out electrical conductivity assignment, set up the true head conductivity model comprising limbic system.
2) the H type coil on tangential described true head structural model surface is set up according to the outline of true head structural model; H type coil is positioned at the anterior cingutate that prefrontal and top conductor line mainly stimulate limbic system; H type coil is positioned at frontal lobe sidepiece wire mainly stimulates the structure such as Hippocampus, corpus amygdaloideum.Consider structure and the position characteristics of head edge system, determine prefrontal, top and sidepiece lead model distribution, the lower edge of H type coil is positioned at the top of head brow ridge, weakens the electric field of H type coil formation to the impact of eyeball.Be based upon the H type coil model at the prefrontal of true head structural model, top and sidepiece different conductor interval.
Described H type coil model, be adopt pulling mode to set up, the square of the conductor cross-section of H type coil model to be the length of side be 1mm, in analog simulation, consider the impact of true coil case and experimenter's hair, H type coil distance true head structural model surface 5mm.Described pulling mode considers that H type coil is positioned at frontal lobe part wire and needs tangential scalp, described point line is carried out according to true head structural model outline in finite element software, form the path on tangential true head conductivity model top layer, set up a cross section perpendicular to path, pull described cross section along path, obtain the H type coil model on tangential true head structural model top layer.
3) use the Electric Field Distribution of finite element method for simulating H type coil on true head conductivity model, the standard evaluating coil deep characteristic be set, comprise:
(1) in finite element software, pass through stress and strain model, true head conductivity model is divided into limited unit, set up the FEM (finite element) model that can apply to numerical computations; FEM (finite element) model is carried out loading solve, obtain the Electric Field Distribution of H type coil on true head model.
(2) in order to meet the object that deep stimulates, reduce the stimulation to peripheral cortex, employing scalp electric field magnitude can as the standard evaluating the modelling of H type coil with the ratio of limbic system electric field magnitude simultaneously.Due to the difference of H type coil conductor spacing, its Electric Field Distribution of bringing out on true head conductivity model slightly difference, therefore utilize the ratio of true head electrical conductivity top layer electric field maximum and true head conductivity model limbic system maximum as deep evaluating characteristics standard, namely evaluate the standard of H type coil design.Described true head conductivity model surface field maximum is less with the ratio of true head conductivity model limbic system electric field maximum, represents that H type coil depth characteristic is better.
4) to evaluate the standard of H type coil design for foundation, the best conductor spacing of H type coil prefrontal, top and sidepiece diverse location is determined.In different H type coil conductor spacings, when the electric field that H type coil produces on true head conductivity model reaches threshold value, if true head conductivity model surface field maximum is close with the ratio of true head conductivity model limbic system electric field maximum, then true head conductivity model surface field smaller is excellent, illustrate under identical depth characteristic, the H type coil model of this conductor spacing is less on the impact of head peripheral cortex.