Method for establishing transcranial magnetic stimulation target spot and repetition frequency
Technical Field
The invention relates to a method for establishing transcranial magnetic stimulation targets and repetition frequencies. In particular to a method for establishing a transcranial magnetic stimulation target point according to principal component analysis and establishing a transcranial magnetic stimulation repetition frequency according to the principal frequency of ganglion rhythms at the target point.
Background
Repeated transcranial magnetic stimulation (repetitive transcranial magnetic stimulation, rTMS) is a brain stimulation technique widely used in the study of brain function, brain network, brain circuit, etc.
rTMS target and repetition frequency are important parameters for rTMS brain studies. In rTMS target selection, it is usually determined from functional partitioning of the cerebral cortex, however, a brain functional region is usually large in area and does not give an accurate location of the rTMS target. In terms of rTMS repetition frequency, it is generally believed that high frequency rTMS is used to predispose to cortical excitability and low frequency rTMS is used to inhibit cortical excitability. There is no clear basis for more specific repetition frequency selection. Currently, in existing rTMS studies, it is mostly determined whether the selected target and the repetition frequency are appropriate by observing the effect of rTMS, and it is not possible to compare all targets and all repetition frequencies within the functional region, nor is it possible to determine whether the selected target and the selected repetition frequency are optimal. Thus, for different studies with the same stimulation objective, the possible rTMS targets and repetition frequencies obtained are not consistent, which makes it difficult for subsequent researchers to select appropriate targets and repetition frequencies based on existing studies, and also makes comparison between multiple studies difficult.
Disclosure of Invention
In order to solve the technical problems, the invention aims to provide a method for establishing a transcranial magnetic stimulation target point and a repetition frequency.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method of establishing transcranial magnetic stimulation targets and repetition frequencies comprising the steps of:
s1, recording nerve electrical signals at a plurality of positions of a brain function target area, performing principal component analysis on the nerve electrical signals at all positions, and determining a dominant nerve electrical signal, wherein the position of the nerve electrical signal is an rTMS target;
s2, calculating the dominant frequency of the nerve electric signal at the target point through Fourier transformation;
s3, establishing an rTMS repetition frequency according to the main frequency and the stimulation purpose.
Further, the S1 specifically includes:
s1.1, recording nerve electrical signals of a plurality of leads of a brain function target area by using brain electricity acquisition equipment, and carrying out preliminary pretreatment, wherein the nerve electrical signals after pretreatment are expressed as X= { X 1 ,X 2 ,……,X N Where N is the number of bioelectric signal leads, which will be determined based on the lead spacing of 5mm, the data length of each lead being N, the jth data point of the ith lead being denoted X ij , i=1,2......,N,j=,2......,N;
S1.2, carrying out principal component analysis on nerve electric signals of all leads, and establishing an rTMS target, wherein the method specifically comprises the following steps:
1) Principal component analysis of all leads in the neuroelectric signal, y=a·x, y= { Y 1 ,Y 2 ,……,Y M M is the number of main components, M is less than or equal to N, Y1 is a first main component, A is a main component coefficient matrix, and A= { A 1 ,A 2 ,……,A M }, A 1 N coefficient values are included for the coefficient sequences corresponding to the first principal component, and represent the contribution of the N-lead nerve electric signals to the first principal component;
2) Finding the maximum value in the N first principal component coefficients, wherein the lead corresponding to the maximum value is a nerve electric signal which plays a dominant role in a brain function target area, and determining the position of the lead as an rTMS target point.
Further, the step S2 specifically includes the following steps:
s2.1 records the nerve electrical signal at the rTMS target as x (k), k=1, 2. N is the number of data points, x (k) is transformed into the frequency domain by Fourier transform, and is marked as f (w);
s2.2 calculating the power spectral density, denoted as p (w), p (w) = |f (w) 2 /n;
S2.3 calculating a power spectral density maximum, denoted P, p=max (P (w)); the frequency value corresponding to the maximum value of the power spectral density is marked as W, and P (W) =p is satisfied; w is the dominant frequency of the nerve electrical signal, 1Hz < W <40Hz.
Further, the step S3 specifically includes:
s3.1 when the stimulation objective is to facilitate cortical excitability, 100% W-120% W is used as the rTMS repetition frequency, denoted HF;
s3.2 when the stimulation is aimed at inhibiting cortical excitability, 80% W-100% W may be used as the rTMS repetition frequency, or the rTMS repetition frequency may be determined from the dominant frequency range of the bioelectric signal at the rTMS target, denoted LF.
Further, when the aim of the stimulation is to inhibit cortex excitability, the rTMS repetition frequency is determined according to the main frequency range of the nerve electric signal at the rTMS target, specifically, the main frequency range of the nerve electric signal at the rTMS target is divided into {1Hz-5Hz,5Hz-10Hz,10Hz-15Hz,15Hz-20Hz,20Hz-25Hz,25Hz-30Hz,35Hz-40Hz }, and the corresponding rTMS repetition frequencies are established to be {1Hz,0.9Hz,0.8Hz,0.7Hz,0.6Hz,0.5Hz,0.4Hz }, respectively.
Compared with the prior art, the invention has the beneficial technical effects that:
according to the invention, the rTMS target is determined according to principal component analysis, and the rTMS repetition frequency is established according to the principal frequency of the target nerve electrical signal, so that a reliable basis is provided for the selection of the rTMS target and the repetition frequency, and the rTMS target has good development and application prospects.
Drawings
The invention is further described with reference to the following description of the drawings.
FIG. 1 is a flow chart of a method of the present invention for establishing transcranial magnetic stimulation targets and repetition frequencies.
Detailed Description
As shown in fig. 1, a method for establishing a transcranial magnetic stimulation target and a repetition frequency, and a method for establishing a transcranial magnetic stimulation target and a repetition frequency, comprising the following steps:
s1, recording nerve electrical signals at a plurality of positions of a brain function target area, performing principal component analysis on the nerve electrical signals at all positions, and determining a dominant nerve electrical signal, wherein the position of the nerve electrical signal is an rTMS target;
s1.1, recording nerve electrical signals of a plurality of leads of a brain function target area by using brain electricity acquisition equipment, and carrying out preliminary pretreatment, wherein the nerve electrical signals after pretreatment are expressed as X= { X 1 ,X 2 ,……,X N Where N is the number of bioelectric signal leads, which will be determined based on the lead spacing of 5mm, the data length of each lead being N, the jth data point of the ith lead being denoted X ij , i=1,2......,N,j=,2......,N;
S1.2, carrying out principal component analysis on nerve electric signals of all leads, and establishing an rTMS target, wherein the method specifically comprises the following steps:
1) Principal component analysis of all leads in the neuroelectric signal, y=a·x, y= { Y 1 ,Y 2 ,……,Y M },M Is the number of main components, M is less than or equal to N, Y1 is a first main component, A is a main component coefficient matrix, and A= { A 1 ,A 2 ,……,A M }, A 1 N coefficient values are included for the coefficient sequences corresponding to the first principal component, and represent the contribution of the N-lead nerve electric signals to the first principal component;
2) Finding the maximum value in the N first principal component coefficients, wherein the lead corresponding to the maximum value is a nerve electric signal which plays a dominant role in a brain function target area, and determining the position of the lead as an rTMS target point.
S2, calculating the dominant frequency of the nerve electric signal at the target point through Fourier transformation;
s2.1 records the nerve electrical signal at the rTMS target as x (k), k=1, 2. N is the number of data points, x (k) is transformed into the frequency domain by Fourier transform, and is marked as f (w);
s2.2 calculating the power spectral density, denoted as p (w), p (w) = |f (w) 2 /n;
S2.3 calculating a power spectral density maximum, denoted P, p=max (P (w)); the frequency value corresponding to the maximum value of the power spectral density is marked as W, and P (W) =p is satisfied; w is the dominant frequency of the nerve electrical signal, 1Hz < W <40Hz.
S3, establishing rTMS repetition frequency according to the main frequency combined with the stimulation purpose
S3.1 when the stimulation objective is to facilitate cortical excitability, 100% W-120% W is used as the rTMS repetition frequency, denoted HF;
s3.2 when the stimulation is aimed at inhibiting cortical excitability, 80% W-100% W may be used as the rTMS repetition frequency, or the rTMS repetition frequency may be determined from the dominant frequency range of the bioelectric signal at the rTMS target, denoted LF.
Wherein when the aim of the stimulation is to inhibit cortex excitability, the rTMS repetition frequency is determined according to the main frequency range of the nerve electric signal at the rTMS target, specifically, the main frequency range of the nerve electric signal at the rTMS target is divided into {1Hz-5Hz,5Hz-10Hz,10Hz-15Hz,15Hz-20Hz,20Hz-25Hz,25Hz-30Hz,35Hz-40Hz } and corresponding rTMS repetition frequencies are established as {1Hz,0.9Hz,0.8Hz,0.7Hz,0.6Hz,0.5Hz,0.4Hz } respectively.
Implementation case: in the occipital brain function target area, the ganglion main frequency is about 10Hz, so for the stimulation purpose of facilitating the brain function area, 10Hz rTMS is selected for stimulation, and for the stimulation purpose of inhibiting the brain function area, 0.9Hz or 0.8Hz rTMS is selected for stimulation; in the target region of prefrontal brain function, the ganglion rhythm has a dominant frequency of about 15Hz, so that 15Hz rTMS is selected for stimulation purposes to facilitate the brain function, and 0.8Hz or 0.7Hz rTMS is selected for stimulation purposes to inhibit the brain function.
The above embodiments are only illustrative of the preferred embodiments of the present invention and are not intended to limit the scope of the present invention, and various modifications and improvements made by those skilled in the art to the technical solutions of the present invention should fall within the protection scope defined by the claims of the present invention without departing from the design spirit of the present invention.