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 see through skull undampedly and stimulate to cerebral nerve.Along with research deeply, the development of technology, the stimulation of shallow-layer can not meet the demand of medical treatment and scientific research, deep stimulation is arisen at the historic moment.Current coil can only stimulate the peripheral cortex region to 1.5-2cm.Deep stimulates in 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 is to stimulate major part, the tangential scalp of wire direction; Outstanding loop and link circuit, away from stimulation target, weaken the impact on target area.But current H type coil designs based on Sphere Model, there is larger difference in Sphere Model and true head structure.The accuracy of head model affects the electric field distribution law that transcranial magnetic stimulation coil produces on model.True head mould is for the design of H type coil provides foundation, and the Electric Field Distribution that coil is produced at head more levels off to truth.
H type coil utilizes the vector superposed principle of induction field, strengthens the induction field in the target tissue of deep; Base portion wire is certain intervals and distributes, and reduces the stimulus intensity of peripheral cortex with this.Conductor spacing and distribution influence electric field intensity and the deep characteristic of peripheral cortex.It is far from being enough that the attenuation characteristic of electric field on a certain line segment only studied in the complexity explanation of head and H type coil, needs more detailed internal electric field distributed intelligence.Comprise Hippocampus, the information exchange of the isostructural limbic system of corpus amygdaloideum between midbrain, diencephalon and neopallium structure plays an important role, be necessary limbic system to separate, set up comprise limbic system true head model.The present invention adopts the method for finite element to calculate head Electric Field Distribution, the impact of the H type coil conductor spacing of assessment prefrontal, top and sidepiece on deep characteristic.
Summary of the invention
Technical problem to be solved by this invention is, provide a kind of by changing the conductor spacing of prefrontal, top and sidepiece, use the method simulate electric field of finite element at the H type coil optimization method for brain deep transcranial magnetic stimulation of distribution situation and the deep characteristic of true conductivity model.
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 distribution situation and the deep characteristic of true head conductivity model, 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 that comprises limbic system, in finite element software, set up the true head conductivity model that comprises limbic system;
2) set up the H type coil of tangential described true head structural model outer surface 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 prefrontal, top and the sidepiece different conductor interval of true head structural model;
3) use the Electric Field Distribution of finite element method for simulating H type coil on true head conductivity model, the standard of evaluating coil deep characteristic be set, comprising:
(1) in finite element software, divide by grid, true head conductivity model is divided into limited unit, foundation can apply to the FEM (finite element) model of numerical computations; In FEM (finite element) model, load and solve, obtain the Electric Field Distribution of H type coil on true head model;
(2) object stimulating in order to meet deep, reduce the stimulation to peripheral cortex simultaneously, adopt the standard of evaluating the design of H type coil that is compared to of true head conductivity model surface field maximum and true head conductivity model limbic system electric field maximum;
4) taking the standard of evaluating H type coil design as foundation, determine the best conductor spacing of H type coil prefrontal, top and sidepiece diverse location.
Step 1) the described foundation head model that comprises DEEP STRUCTURE comprises:
(1) by nuclear magnetic resonance or the interactively medical image control system of CT data importing, comprise Threshold segmentation, region growing and the operation of manually modifying, extract the edge of head and limbic system, set up respectively the three-dimensional patch model of head and the three-dimensional patch model of limbic system;
(2) in finite element software, be head physical model and limbic system physical model by three-dimensional head patch model and the three-dimensional dough sheet model conversation of limbic system, described head physical model and limbic system physical model are overlapped to Boolean calculation, the true head structural model that foundation comprises limbic system, carry out electrical conductivity assignment, set up the true head conductivity model that comprises limbic system.
Step 2) described H type coil model, be to adopt pulling mode to set up, the conductor cross-section of H type coil model is that the length of side is the square of 1mm, conductor spacing true head structural model surface 5mm.
Described pulling mode is to carry 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 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.
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, the conductor spacing of change prefrontal, top and sidepiece, the impact of assessment conductor spacing on electric field deep characteristic.Due to the complexity of head construction, it is slightly variant that the wire of different interval produces Electric Field Distribution, therefore adopts the standard of evaluating coil deep characteristic that is compared to of true head conductivity model surface field maximum and limbic system maximum.
Brief description of the drawings
Fig. 1 is the flow chart of the inventive method;
Fig. 2 is head model and the H type coil model that adopts 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 structure under limbic lobe and relevant cortex and cortex, relevant with higher nervous activities such as emotion, spirit, memories.Due to the difference of arrangement of conductors, in limbic system, produce different effects.In order to reach the darker stimulation degree of depth, reduce the electric field intensity of top layer scalp simultaneously.H type coil conductor spacing to 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 distribution situation and the deep characteristic of true conductivity model, 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 that comprises limbic system, in finite element software, set up the true head conductivity model that comprises limbic system;
The head model that described foundation comprises DEEP STRUCTURE comprises:
(1) by nuclear magnetic resonance or the interactively medical image control system of CT data importing (Mimics), comprise Threshold segmentation, region growing and the operation of manually modifying, extract the edge of head and limbic system.After three-dimensional surface modeling, the model obtaining does not pass through any optimization, do not affecting on the basis of head model structure, adopt smooth operation to come level and smooth solid surface, simplify operation by dough sheet and reduce tri patch quantity, set up respectively the three-dimensional patch model of head and the three-dimensional patch model of limbic system;
(2) in finite element software, be head physical model and limbic system physical model by three-dimensional head patch model and the three-dimensional dough sheet model conversation of limbic system, described head physical model and limbic system physical model are overlapped to Boolean calculation, ensure node and unit on two physical model Border faces, the true head structural model that foundation comprises limbic system, carry out electrical conductivity assignment, set up the true head conductivity model that comprises limbic system.
2) set up the H type coil on tangential described true head structural model surface according to the outline of true head structural model; H type coil is positioned at prefrontal and top conductor line mainly stimulates the anterior cingutate of 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 impact on eyeball of electric field that H type coil forms.Be based upon the H type coil model at prefrontal, top and the sidepiece different conductor interval of true head structural model.
Described H type coil model, be to adopt pulling mode to set up, the conductor cross-section of H type coil model is that the length of side is the square of 1mm, in analog simulation, consider the impact of true coil case and experimenter's hair, H type coil is apart from the surperficial 5mm of true head structural model.Described pulling mode is to consider that H type coil is positioned at frontal lobe part wire and needs tangential scalp, in finite element software, carry out described point line according to true head structural model outline, 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 of evaluating coil deep characteristic be set, comprising:
(1) in finite element software, divide by grid, true head conductivity model is divided into limited unit, foundation can apply to the FEM (finite element) model of numerical computations; In FEM (finite element) model, load and solve, obtain the Electric Field Distribution of H type coil on true head model.
(2) object stimulating in order to meet deep reduces the stimulation to peripheral cortex simultaneously, adopts scalp electric field magnitude and the ratio of limbic system electric field magnitude to can be used as the standard of evaluating the modelling of H type coil.Due to the difference of H type coil conductor spacing, its Electric Field Distribution of bringing out on true head conductivity model is slightly variant, therefore utilize the deep evaluating characteristics standard that is compared to of true head electrical conductivity top layer electric field maximum and true head conductivity model limbic system maximum, 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) taking the standard of evaluating H type coil design as foundation, determine the best conductor spacing of H type coil prefrontal, top and sidepiece diverse location.In different H type coil conductor spacings, the electric field producing on true head conductivity model at H type coil reaches threshold value, if true head conductivity model surface field maximum approaches with the ratio of true head conductivity model limbic system electric field maximum, 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.