Calculation method for carrying out target precise stimulation aiming at time interference stimulation

Abstract

The invention provides a calculation method for performing target precise stimulation for time interference stimulation. The calculation method comprises the following steps: (1) segmenting a head magnetic resonance image according to different tissues; (2) gridding the image based on the segmentation result of the magnetic resonance image, and calculating a current guide field matrix corresponding to each electrode by using the electrodes; (3) calculating to obtain the electric field vector intensity of all gridding coordinates when the electrode pair is used for electrical stimulation; (4) calculating the electric field vector intensity in the target area and the envelope electric field intensity of all grid coordinates; (5) calculating the average value and the maximum value of the envelope electric field intensity of all the grid coordinates of all the electrode pair combinations in the target area, and screening all the electrode pair combinations to obtain a first screened electrode pair combination; and (6) calculating the average value and the maximum value of the envelope electric field intensity of the gray matter area when the electrode pair is subjected to electrical stimulation, and screening the first screening electrode pair combination to obtain an optimal electrode pair combination.

Description

Technical Field

[0001] This invention belongs to the field of neurostimulation technology, and in particular relates to a calculation method for precise targeted stimulation of time-interference stimuli. Background Technology

[0002] Temporal interference stimulation (TIS) is a novel neuromodulation technique based on transcranial electrical stimulation (TCS). This technique uses two high-frequency alternating current stimuli of similar frequencies, superimposed and coupled to form a low-frequency envelope electric field, which stimulates neurons. Based on the strong penetrability of high-frequency signals to human tissues and the high controllability of low-frequency signals to neurons, TIS can overcome the limitation of existing TCS which only stimulates the surface of the cerebral cortex, enabling stimulation of deep brain nuclei and tissues.

[0003] Due to the complexity of the human brain's structure, and the fact that the envelope signal formed by high-frequency alternating current transmitted through the scalp to the deep brain can only be calculated through physical field simulation, the simulated physical fields obtained from different scalp stimulation points vary significantly. Furthermore, to achieve optimal stimulation of neurons in the target region, the envelope electric field strength in the target region needs to be sufficiently large, while minimizing the electric field strength in non-target regions. Currently, a significant amount of time is required to search for the optimal electrode stimulation points to obtain the best stimulation effect. Summary of the Invention

[0004] The purpose of this invention is to address the above problems by proposing a calculation method for precise target stimulation based on temporal interference stimuli. Under a defined target stimulation area, the optimal scalp electrode stimulation point can be quickly found, maximizing the intensity of the envelope signal in the target area and minimizing the intensity in other areas, thereby achieving the effect of precise target stimulation.

[0005] To achieve the above objectives, this invention proposes a calculation method for precise target stimulation of temporal interference stimuli, which includes the following steps:

[0006] (1) Segment the head magnetic resonance images according to different tissues;

[0007] (2) The segmented image is meshed, and the current guiding field matrix corresponding to each electrode is calculated according to the electrode position of the 10-10 standard leads;

[0008] (3) Four different electrode combinations are combined to form 6 electrode pair combinations. All electrode combinations in the 10-10 standard leads are obtained by arranging and combining them. When all electrode pairs are electrically stimulated, the electric field vector intensity of all gridded coordinates is calculated.

[0009] (4) The combination of multiple grid coordinates is called the target area. Using the method in step (3), the electric field vector intensity in the target area is calculated when the electrode pairs in all 6 electrode pair combinations are electrically stimulated. Based on this electric field vector intensity, the envelope electric field intensity of all grid coordinates in the target area of ​​all 6 electrode pair combinations is calculated.

[0010] (5) Based on the envelope electric field intensity of all grid coordinates, calculate the average and maximum values ​​of the envelope electric field intensity of all 6 electrode pairs in the target area at all grid coordinates. Based on the above two electric field intensities, screen all 6 electrode pairs in the total electrode combination to obtain the first screened electrode pair combination.

[0011] (6) Calculate the average and maximum values ​​of the envelope electric field strength of the gray matter region when the electrode pairs in the first screening electrode pair combination are electrically stimulated, and based on the average and maximum values ​​of the envelope electric field strength of the gray matter region and the target region, screen the first screening electrode pair combination to obtain the optimal electrode pair combination.

[0012] Furthermore, in step (1), the head image is segmented into brain white matter, brain gray matter, cerebrospinal fluid, skull, and scalp tissue.

[0013] Furthermore, the specific method for step (3) is as follows:

[0014] Four different electrodes are used as electrode combinations. All electrode combinations in the 10-10 standard leads are obtained by arranging and combining them. The four electrodes are paired up to form two electrode pairs, resulting in a total of six electrode pair combinations. For any electrode pair, the electric field vector intensity of all gridded coordinates when using the electrode pair for electrical stimulation is obtained using the current guiding field matrix calculated in step (2).

[0015] Furthermore, in step (4), the formula (1) for the envelope electric field strength is as follows:

[0016] (1)

[0017] in, and Let be the electric field vector intensity of two electrode pairs in an electrode pair combination under the s coordinate, where s represents the coordinate value; Representing two electric fields and The spatial angle, That is, the value of the envelope electric field strength formed by this coordinate system.

[0018] Furthermore, in step (5), the specific method for obtaining the first screening electrode pair combination is as follows:

[0019] (5.1) For all electrode pair combinations, retain the electrode pair combinations with the maximum envelope electric field strength in the target region that is greater than the 80th percentile.

[0020] (5.2) The first screening electrode pair combination is obtained by taking the average value of the envelope electric field strength of the electrode pair combination retained in step (5.1) in the target area according to the preset conditions.

[0021] Furthermore, the preset condition is: retain electrode pair combinations with an average envelope electric field strength greater than 0.4 V / m or a strength greater than the 90th percentile.

[0022] Furthermore, the specific method for step (6) is as follows:

[0023] (6.1) Based on the maximum value of the envelope electric field strength of the first screening electrode pair combination in the gray matter region, retain electrode pair combinations with a strength less than 10% of the quantiles;

[0024] (6.2) Based on the average value of the envelope electric field intensity of the electrode pair combination retained in step (6.1) in the gray matter region, the retained intensity is less than that of the electrode pair combination at the 10th quantile.

[0025] (6.3) Based on the electrode pair combinations retained in step (6.2), the average value of the envelope electric field intensity of each electrode pair in the gray matter region is called E. GM The average value of the envelope electric field strength in the target region is called E. ROI Calculate E ROI / E GM , keep E ROI / E GM Electrode pair combinations greater than the 95th quantile;

[0026] (6.4) Based on the electrode pair combinations retained in step (6.3), the electrode pair combination with the maximum average value of the envelope electric field intensity in the target region is the optimal combination.

[0027] Furthermore, the present invention proposes an apparatus comprising a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the steps of any of the aforementioned methods for calculating precise target stimulation for temporal interference stimuli.

[0028] Furthermore, the present invention proposes a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of any of the aforementioned methods for calculating precise target stimulation for time-interference stimuli.

[0029] Beneficial effects: Compared with the prior art, the technical solution of the present invention has the following beneficial technical effects:

[0030] The method proposed in this invention comprehensively considers the electric field intensity of the target region and the electric field intensity of the gray matter region, ensuring that the electric field intensity is maximized in the target region and minimized in other regions, thereby achieving rapid search for the location of the optimal electrode stimulation point. Attached Figure Description

[0031] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the embodiments are briefly described below.

[0032] Figure 1 This is a flowchart of the method of the present invention. Detailed Implementation

[0033] To more clearly illustrate the technical solution of the present invention, the present invention will be further described below in conjunction with the accompanying drawings and embodiments. The following embodiments will help those skilled in the art to further understand the present invention, but do not limit the present invention in any way. It should be noted that those skilled in the art can make several modifications and improvements without departing from the concept of the present invention. These all fall within the protection scope of the present invention.

[0034] like Figure 1 As shown, this invention proposes a method for calculating precise target stimulation for time-interference stimuli, which includes the following steps:

[0035] (1) Segment the head magnetic resonance images according to different tissues;

[0036] (2) The segmented image is meshed, and the current guiding field matrix corresponding to each electrode is calculated according to the electrode position of the 10-10 standard leads;

[0037] (3) Four different electrode combinations are combined to form 6 electrode pair combinations. All electrode combinations in the 10-10 standard leads are obtained by arranging and combining them. When all electrode pairs are electrically stimulated, the electric field vector intensity of all gridded coordinates is calculated.

[0038] (4) The combination of multiple grid coordinates is called the target area. Using the method in step (3), the electric field vector intensity in the target area is calculated when the electrode pairs in all 6 electrode pair combinations are electrically stimulated. Based on this electric field vector intensity, the envelope electric field intensity of all grid coordinates in the target area of ​​all 6 electrode pair combinations is calculated.

[0039] (5) Based on the envelope electric field intensity of all grid coordinates, calculate the average and maximum values ​​of the envelope electric field intensity of all 6 electrode pairs in the target area at all grid coordinates. Based on the above two electric field intensities, screen all 6 electrode pairs in the total electrode combination to obtain the first screened electrode pair combination.

[0040] (6) Calculate the average and maximum values ​​of the envelope electric field strength of the gray matter region when the electrode pairs in the first screening electrode pair combination are electrically stimulated, and based on the average and maximum values ​​of the envelope electric field strength of the gray matter region and the target region, screen the first screening electrode pair combination to obtain the optimal electrode pair combination.

[0041] Furthermore, the specific method for step (3) is as follows:

[0042] Four different electrodes are used as electrode combinations. All electrode combinations in the 10-10 standard leads are obtained by arranging and combining them. The four electrodes are paired up to form two electrode pairs, resulting in a total of six electrode pair combinations. For any electrode pair, the electric field vector intensity of all gridded coordinates when using the electrode pair for electrical stimulation is obtained using the current guiding field matrix calculated in step (2).

[0043] Furthermore, in step (4), the formula (1) for the envelope electric field strength is as follows:

[0044] (1)

[0045] in, and Let be the electric field vector intensity of two electrode pairs in an electrode pair combination under the s coordinate, where s represents the coordinate value; Representing two electric fields and The spatial angle, That is, the value of the envelope electric field strength formed by this coordinate system.

[0046] Furthermore, in step (5), the specific method for obtaining the first screening electrode pair combination is as follows:

[0047] (5.1) For all electrode pair combinations, retain the electrode pair combinations with the maximum envelope electric field strength in the target region that is greater than the 80th percentile.

[0048] (5.2) The first screening electrode pair combination is obtained by taking the average value of the envelope electric field strength of the electrode pair combination retained in step (5.1) in the target area according to the preset conditions.

[0049] Furthermore, the preset condition is: retain electrode pair combinations with an average envelope electric field strength greater than 0.4 V / m or a strength greater than the 90th percentile.

[0050] Furthermore, the specific method for step (6) is as follows:

[0051] (6.1) Based on the maximum value of the envelope electric field strength of the first screening electrode pair combination in the gray matter region, retain electrode pair combinations with a strength less than 10% of the quantiles;

[0052] (6.2) Based on the average value of the envelope electric field intensity of the electrode pair combination retained in step (6.1) in the gray matter region, the retained intensity is less than that of the electrode pair combination at the 10th quantile.

[0053] (6.3) Based on the electrode pair combinations retained in step (6.2), the average value of the envelope electric field intensity of each electrode pair in the gray matter region is called E. GM The average value of the envelope electric field strength in the target region is called E. ROI Calculate E ROI / E GM , keep E ROI / E GM Electrode pair combinations greater than the 95th quantile;

[0054] (6.4) Based on the electrode pair combinations retained in step (6.3), the electrode pair combination with the maximum average value of the envelope electric field intensity in the target region is the optimal combination.

[0055] Furthermore, the present invention proposes an apparatus comprising a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the steps of any of the aforementioned methods for calculating precise target stimulation for temporal interference stimuli.

[0056] Furthermore, the present invention proposes a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of any of the aforementioned methods for calculating precise target stimulation for time-interference stimuli.

[0057] The specific embodiments of the present invention have been described above. It should be understood that the present invention is not limited to the specific embodiments described above, and those skilled in the art can make various modifications or variations within the scope of the claims, which do not affect the essence of the present invention.

Claims

1. A calculation method for precise target stimulation of temporally interfering stimuli, characterized in that, The method includes the following steps: (1) Segment the head magnetic resonance images according to different tissues; (2) The segmented image is meshed, and the current guiding field matrix corresponding to each electrode is calculated according to the electrode position of the 10-10 standard leads; (3) Four different electrode combinations are combined to form 6 electrode pair combinations. All electrode pair combinations in the 10-10 standard leads are obtained by arranging and combining them. When all electrode pairs are electrically stimulated, the electric field vector intensity of all gridded coordinates is calculated. (4) The combination of multiple grid coordinates is called the target area. Using the method in step (3), the electric field vector intensity in the target area is calculated when the electrode pairs in all 6 electrode pair combinations are electrically stimulated. Based on this electric field vector intensity, the envelope electric field intensity of all grid coordinates in the target area of ​​all 6 electrode pair combinations is calculated. (5) Based on the envelope electric field intensity of all grid coordinates, calculate the average and maximum values ​​of the envelope electric field intensity of all 6 electrode pair combinations in the target area. Based on the above two electric field intensities, screen all 6 electrode pair combinations to obtain the first screened electrode pair combination. (6) Calculate the average and maximum values ​​of the envelope electric field strength of the gray matter region when the electrode pairs in the first screening electrode pair combination are electrically stimulated, and based on the average and maximum values ​​of the envelope electric field strength of the gray matter region and the target region, screen the first screening electrode pair combination to obtain the optimal electrode pair combination.

2. The method for calculating precise target stimulation for time-interference stimuli according to claim 1, characterized in that, The specific method for step (3) is as follows: Four different electrodes are used as electrode combinations. All electrode combinations in the 10-10 standard leads are obtained by arranging and combining them. The four electrodes are paired up to form two electrode pairs, resulting in a total of six electrode pair combinations. For any electrode pair, the electric field vector intensity of all gridded coordinates when using the electrode pair for electrical stimulation is obtained using the current guiding field matrix calculated in step (2).

3. The method for calculating precise target stimulation for time-interference stimuli according to claim 1, characterized in that, In step (4), the formula (1) for the envelope electric field strength is as follows: (1) in, and Let be the electric field vector intensity of two electrode pairs in an electrode pair combination under the s coordinate, where s represents the coordinate value; Representing two electric fields and The spatial angle, That is, the value of the envelope electric field strength formed by this coordinate system.

4. The method for calculating precise target stimulation for time-interference stimuli according to claim 1, characterized in that, The specific method for obtaining the first screening electrode pair combination in step (5) is as follows: (5.1) For all electrode pair combinations, retain the electrode pair combinations with the maximum envelope electric field strength in the target region that is greater than the 80th percentile. (5.2) The first screening electrode pair combination is obtained by taking the average value of the envelope electric field strength of the electrode pair combination retained in step (5.1) in the target area according to the preset conditions.

5. The method for calculating precise target stimulation for time-interference stimuli according to claim 4, characterized in that, The preset condition is: retain electrode pair combinations with an average envelope electric field strength greater than 0.4 V / m or a strength greater than the 90th percentile.

6. The method for calculating precise target stimulation for time-interference stimuli according to claim 1, characterized in that, The specific method for step (6) is as follows: (6.1) Based on the maximum value of the envelope electric field strength of the first screening electrode pair combination in the gray matter region, retain electrode pair combinations with a strength less than 10% of the quantiles; (6.2) Based on the average value of the envelope electric field intensity of the electrode pair combination retained in step (6.1) in the gray matter region, the retained intensity is less than that of the electrode pair combination at the 10th quantile. (6.3) Based on the electrode pair combinations retained in step (6.2), the average value of the envelope electric field intensity of each electrode pair in the gray matter region is called E. GM The average value of the envelope electric field strength in the target region is called E. ROI Calculate E ROI / E GM , keep E ROI / E GM Electrode pair combinations greater than the 95th quantile; (6.4) Based on the electrode pair combinations retained in step (6.3), the electrode pair combination with the maximum average value of the envelope electric field intensity in the target region is the optimal combination.

7. An apparatus, characterized in that, The device includes a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the steps of the calculation method for targeted precise stimulation of any of the time-interference stimuli as described in claims 1-6.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the steps of the calculation method for precise target stimulation for any of the time-interference stimuli described in claims 1-6.

Images