Method and device for modeling anisotropic conductivity head model by combining cortex excitability
Technical Field
The invention relates to the field of transcranial magnetic stimulation head model modeling, in particular to an anisotropic conductivity head model modeling method and device combined with cortex excitability.
Background
Transcranial magnetic stimulation is a non-invasive neurostimulation technology and is widely applied to researches on brain functions, brain networks, brain circuits and the like. The technology utilizes a variable magnetic field generated by pulse current to act on brain tissue so as to generate an induced electric field, and when the induced electric field in a brain tissue region reaches or exceeds a certain threshold value, the neural activity of the region can be activated or inhibited.
TMS simulation research is an important means for predicting TMS effect and optimizing TMS coil and stimulation parameters. The accurate establishment of the TMS head model is the basis for carrying out TMS intracranial field distribution simulation operation. The current TMS head model is usually constructed by adopting individualized MRI T1 and DTI, mainly focuses on structural information such as sulcus, fasciculus and the like, and reflects the difference of different individuals on brain structures. For the conductivity assignment of the head model, the conductivity assignment is uniform among different individuals and does not reflect the difference of cortical excitability of different individuals.
Therefore, we propose an anisotropic conductivity head model modeling method and apparatus combining cortical excitability to solve the above problems.
Disclosure of Invention
The invention aims to solve the defect that the difference of the excitation of different individual cortex cannot be reflected in the prior art, and provides a method and a device for modeling an anisotropic conductivity head model in combination with the excitation of the cortex.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method of modeling an anisotropic conductivity head model in conjunction with cortical excitability, the method comprising:
establishing an anisotropic conductivity head model based on MRI and DTI images;
acquiring an individual cortex excitability comprehensive index through an array target TMS experiment;
and constructing an anisotropic conductivity head model containing cortex excitability according to the anisotropic conductivity head model and the individual cortex excitability comprehensive index.
Further, the establishing of the anisotropic conductivity head model based on MRI and DTI images specifically includes:
MRI is carried out by adopting a T1 preparation to obtain an MRI T1 image, and the MRI T1 image is subjected to image correction, tissue segmentation and cortical reconstruction to obtain a real structural head model;
decoupling, surface net separating and body net separating operations are carried out on the head model to obtain an isotropic head model;
extracting a diffusion tensor according to the DTI data, and converting the diffusion tensor into a conductivity tensor;
and combining the conductivity tensor with the isotropic head model to obtain an anisotropic conductivity head model.
Further, the obtaining of the individual cortical excitability comprehensive index through the array target TMS experiment specifically comprises:
selecting a stimulation target area as a primary motor cortex, moving a TMS coil to each target point, implementing single TMS on each target point, and recording a motion threshold and a motion evoked potential of each target point during single TMS;
calculating a target area density mean value PAVE according to the motion threshold value, calculating discrete values of the motion threshold value and the target area density mean value PAVE, determining a target area density center according to the discrete values, and obtaining the individual cortical excitability comprehensive index according to the motion threshold value and the motion evoked potential at the target area density center.
Furthermore, the method for obtaining the individual cortical excitability comprehensive index through the array target TMS experiment is characterized by comprising the following steps:
the motion-evoked potentials include a motion-evoked potential amplitude and a motion-evoked potential latency.
Further, the constructing an anisotropic conductivity head model containing cortical excitability according to the anisotropic conductivity head model and the individual cortical excitation comprehensive index specifically includes:
setting default conductivity of each layer of brain tissue in the anisotropic head model, and calculating the mean value of the individual cortical excitation comprehensive index as a total cortical excitation comprehensive index;
and acquiring an individual conductivity coefficient according to the individual cortex excitation comprehensive index and the overall cortex excitation comprehensive index, and acquiring an anisotropic conductivity head model of cortex excitability according to the default conductivity and the individual conductivity coefficient.
Anisotropic conductivity head model modeling apparatus incorporating cortical excitability, comprising:
the first establishing module is used for establishing an anisotropic conductivity head model based on MRI and DTI images;
the acquisition module is used for acquiring an individual cortical excitability comprehensive index through an array target TMS experiment;
and the second establishing module is used for establishing an anisotropic conductivity head model containing cortex excitability according to the anisotropic conductivity head model and the individual cortex excitation comprehensive index. .
Further, the first establishing module specifically includes:
the isotropic head model establishing module is used for carrying out MRI by adopting a T1 preparation to obtain an MRI T1 image, and carrying out image correction, tissue segmentation and cortical reconstruction on the MRI T1 image to obtain a real structural head model;
decoupling, surface net separating and body net separating operations are carried out on the head model to obtain an isotropic head model;
extracting a diffusion tensor according to the DTI data, and converting the diffusion tensor into a conductivity tensor;
and the isotropic conductivity head model establishing module is used for combining the conductivity tensor and the isotropic conductivity head model to obtain an anisotropic conductivity head model.
Further, the acquiring module specifically includes:
the first acquisition module selects the stimulation target area as a primary motor cortex, moves the TMS coil to each target point, implements single TMS on each target point, and records the motion threshold and the motion evoked potential of each target point during single TMS;
the second acquisition module is used for calculating a target area density mean value PAVE according to the motion threshold value, calculating discrete values of the motion threshold value and the target area density mean value PAVE, determining a target area density center according to the discrete values, and acquiring the individual reflection cortical excitability comprehensive index according to the motion threshold value and the motion evoked potential at the target area density center.
Further, the obtaining module is characterized in that:
the motion-evoked potentials include a motion-evoked potential amplitude and a motion-evoked potential latency.
Further, the second establishing module specifically includes:
the third acquisition module is used for setting default conductivity of each layer of brain tissue in the anisotropic head model and calculating the mean value of the individual cortical excitation comprehensive index to serve as a total cortical excitation comprehensive index;
and the anisotropic conductivity head model establishing module is used for acquiring an individualized conductivity coefficient according to the individual cortical excitation comprehensive index and the overall cortical excitation comprehensive index and acquiring the anisotropic conductivity head model of cortical excitability according to the default conductivity and the individualized conductivity coefficient.
The method takes a motion threshold, a motion evoked potential amplitude and a motion evoked potential latency in an array target TMS experiment as cortical excitation comprehensive indexes; and realizing individual conductivity assignment according to the cortical excitation comprehensive index, and constructing an anisotropic conductivity head model containing individual cortical excitation. The transcranial magnetic stimulation head model modeling method combines functional information on the basis of structural information, reflects individual cortical excitation, is closer to a real brain, can further improve the authenticity and reliability of TMS modeling simulation, has important research significance in simulating TMS effect, optimizing TMS coil and stimulation parameters, and has good development and application prospects.
Drawings
FIG. 1 shows a flow chart of an anisotropic conductivity head model modeling method in combination with cortical excitability in example 1 of the present invention;
FIG. 2 is a flowchart showing a specific example of the method for modeling an anisotropic conductivity head model in combination with cortical excitability in example 1 of the present invention;
fig. 3 is a block diagram showing an anisotropic conductivity head model modeling apparatus incorporating cortical excitability in embodiment 2 of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments.
Examples1
The embodiment provides a better anisotropic conductivity head model modeling method combined with cortex excitability, which is used in an anisotropic conductivity head model modeling device combined with cortex excitability. As shown in fig. 1 and 2, the method comprises the following steps:
s101: establishing an anisotropic conductivity head model based on MRI and DTI images;
further, MRI is carried out by adopting a T1 preparation to obtain an MRI T1 image, and the MRI T1 image is subjected to image correction, tissue segmentation and cortical reconstruction to obtain a real structural head model;
decoupling, surface net separating and body net separating operations are carried out on the head model to obtain an isotropic head model;
extracting a diffusion tensor according to the DTI data, and converting the diffusion tensor into a conductivity tensor;
the conductivity tensor is combined with the isotropic head model to obtain an anisotropic conductivity head model.
S102: acquiring an individual cortex excitability comprehensive index through an array target TMS experiment;
further, selecting a stimulation target area as a primary motor cortex, moving a TMS coil to each target point, performing single TMS on each target point, and recording a motion threshold and a motion evoked potential of each target point during single TMS, wherein the motion evoked potential comprises a motion evoked potential amplitude and a motion evoked potential latency;
specifically, the method comprises the following steps: selecting the stimulation target area as primary motor cortex, wherein the target area is rectangular and longnWide ismThe number of target points isNThen, thenN=n*mThe target point coordinates are expressed as (i, j),i=1…n,j=1…m(ii) a The TMS coil is controlled by the mechanical arm to move to each target point, single TMS is carried out on each target point, the motion threshold value of each target point during single TMS is recorded, and the motion threshold value is counted asP ij Recording the motion evoked potential amplitude and the motion evoked potential latency of each target spot;
calculating a target area density mean value PAVE according to the motion threshold value, calculating discrete values of the motion threshold value and the target area density mean value PAVE, determining a target area density center according to the discrete values, and obtaining an individual reflection cortical excitability comprehensive index according to the motion threshold value and the motion evoked potential at the target area density center;
specifically, the method comprises the following steps: calculating the density mean value of the target area according to the motion threshold value of each target point
PAVE,
,
Calculating the motion threshold of each target point
P ij And
PAVEand calculating the absolute value of the difference value, and recording the absolute value as a discrete value
E ij ,
,
Discrete valueE ij The target point corresponding to the minimum value of the target area is determined as a target area density center, a motion threshold value, a motion-evoked potential amplitude and a motion-evoked potential latency at the target area density center are used for evaluating cortical excitability, and the motion threshold value is counted asPC,The motion-evoked potential amplitude is recorded asAC,The motor evoked potential latency is recordedTCDefining a comprehensive index reflecting individual cortical excitabilityg,g=AC/(PC*TC)
S103: constructing an anisotropic conductivity head model containing cortex excitability according to the anisotropic conductivity head model and the individual cortex excitability comprehensive index;
further, setting default conductivity of each layer of brain tissue in the anisotropic head model, and calculating the mean value of the individual cortical excitation comprehensive index as the overall cortical excitation comprehensive index;
specifically, default electrical conductivities of layers of brain tissue in the anisotropic head model are sets 0 Calculating the cortical excitation comprehensive index of all individualsgIs a mean value ofGAs a global cortical stimulation composite indicator, with a default conductivitys 0 And (7) corresponding.
And acquiring an individual conductivity coefficient according to the individual cortex excitation comprehensive index and the overall cortex excitation comprehensive index, and acquiring an anisotropic conductivity head model of cortex excitability according to the default conductivity and the individual conductivity coefficient.
Specifically, the cortical excitation comprehensive index of the individualgCombined index of total cortical excitationGAs an individualized conductivity coefficientk,k=g/GCalculating the product of the individualized conductivity coefficient and the default conductivity as the conductivity estimate containing the individualized cortical excitabilitys i =k*s 0 And finally constructing an anisotropic head model containing the individual TMS cortex excitability.
Example 2
The embodiment provides an anisotropic conductivity head model modeling device which is better combined with cortical excitability. As shown in fig. 3, the apparatus includes:
s201: the first establishing module is used for establishing an anisotropic conductivity head model based on MRI and DTI images;
further, an isotropic head model building module adopts a T1 preparation to carry out MRI to obtain an MRI T1 image, and carries out image correction, tissue segmentation and cortical reconstruction on the MRI T1 image to obtain a real structural head model;
decoupling, surface net separating and body net separating operations are carried out on the head model to obtain an isotropic head model;
extracting a diffusion tensor according to the DTI data, and converting the diffusion tensor into a conductivity tensor;
the isotropic conductivity head model building module is used for combining the conductivity tensor and the isotropic conductivity head model to obtain an anisotropic conductivity head model;
s202: the acquisition module is used for acquiring an individual cortical excitability comprehensive index through an array target TMS experiment;
further, the first acquisition module selects a stimulation target area as a primary motor cortex, moves the TMS coil to each target point, performs single TMS on each target point, and records a motor threshold and a motor evoked potential of each target point during single TMS, wherein the motor evoked potential comprises a motor evoked potential amplitude and a motor evoked potential latency;
the second acquisition module is used for calculating a target area density mean value PAVE according to the motion threshold value, calculating a discrete value of the motion threshold value and the target area density mean value PAVE, determining a target area density center according to the discrete value, and acquiring an individual reflection cortical excitability comprehensive index according to the motion threshold value and the motion evoked potential at the target area density center.
S203: the second establishing module is used for establishing an anisotropic conductivity head model containing cortex excitability according to the anisotropic conductivity head model and the individual cortex excitation comprehensive index;
further, a third acquisition module sets default conductivity of each layer of brain tissue in the anisotropic head model, and calculates a mean value of the individual cortical excitation comprehensive index as a total cortical excitation comprehensive index;
and the anisotropic conductivity head model establishing module is used for acquiring an individualized conductivity coefficient according to the individual cortical excitation comprehensive index and the overall cortical excitation comprehensive index and acquiring an anisotropic conductivity head model of cortical excitability according to the default conductivity and the individualized conductivity coefficient.
Embodiments of the present invention further provide a non-transitory computer storage medium storing computer-executable instructions, where the computer-executable instructions may implement any one of the embodiments. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a Flash Memory (Flash Memory), a Hard Disk (Hard Disk Drive, HDD), a Solid-State Drive (SSD), or the like; the storage medium may also comprise a combination of memories of the kind described above.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a read-only memory (ROM), a Random Access Memory (RAM), or the like.
The above description is only a few preferred embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can combine, substitute, or change the technical solutions and their inventive concepts within the technical scope of the present invention.