Multivariable neural decoding method based on OPM-MEG and RSA

Abstract

The invention discloses a multivariable neural decoding method based on OPM-MEG and RSA, and belongs to the technical field of biomedical engineering. The method comprises the following steps: acquiring resting state structural MRI data and an OPM-MEG scanning image for joint registration; playing auditory stimulation, completing OPM-MEG data acquisition, and preprocessing the data to obtain segmented data; averaging different test times of the same auditory stimulation on a segment data level, constructing a neural response mode vector, and calculating the Pearson correlation of all possible pairwise stimulation pairs at time points to obtain a representation similarity matrix and a neural response mode similarity average result; determining time points with differences one by one; a sequence of source activation times for each subject is estimated, and activation regions are determined using T-tests. According to the method, the decoding representation process of the brain can be disclosed by calculating the distance of the neural reaction mode in the representation space, so that the response mode of the brain under a specific task or condition can be comprehensively understood.

Description

Technical Field

[0001] The invention belongs to the technical field of biomedical engineering, and in particular relates to a multivariate neural decoding method based on OPM-MEG and RSA. Background Art

[0002] In the field of neuroscience, the brain converts external stimuli into intentions through decoding. By decoding neural activity, we can understand how the brain processes information and generates thoughts and behaviors. This is of great significance to many fields, including neuroscience, psychology, medicine, engineering, and computer science.

[0003] In recent years, a variety of measurement modes have been developed to decode brain neural activity. In non-invasive brain function measurement, optically pumped magnetometer-based magnetoencephalography (OPM-MEG) can directly measure the magnetic activity of neurons and capture the dynamic changes of the brain at the millisecond level. It has become an important brain function imaging technology and has broad application prospects.

[0004] In addition to measurement patterns, methods for decoding neural activity have also made significant progress, mainly including two categories: univariate analysis and multivariate analysis. Among them, multivariate analysis can simultaneously consider the activity patterns of multiple brain regions, integrate the overall information of the brain, and more comprehensively understand the brain's activity patterns under specific tasks or conditions.

[0005] Representation Similarity Analysis (RSA), as a method of multivariate analysis, provides a powerful explanatory framework through the core concept of "representation space", in which the neural response to each stimulus is represented as a high-dimensional pattern vector.

[0006] The traditional application scenario of RSA is fMRI data in the visual field. Compared with the visual field, the current research on brain decoding in the auditory field using RSA has two problems: on the one hand, there is a lack of description of how the human brain processes sound. Sound is the basis of language and communication and contains social and biological related information. Exploring how the brain processes sound and what characteristics can better understand the process of sound understanding and recognition is crucial to understanding the functions of the language and auditory systems. Although a variety of multivariate analysis methods have been developed, only a few studies have applied these methods to time series neuroimaging data to date; on the other hand, some existing studies on RSA decoding in the auditory field focus more on locating brain areas with different response characteristics to stimuli, but this brings many problems, one of which is the fuzzy interpretation of adjacent time points. If a faster change in brain activity occurs between two time points, fMRI may not be able to distinguish such changes, but blur them into an overall activity pattern, which may lead to underestimation of the differences between adjacent time points in similarity analysis, omission of certain dynamic process information, and limit the interpretation of pattern similarity. Thus, although RSA is a promising technique for the analysis of fMRI data, it does not necessarily allow for specific conclusions to be drawn about how information is encoded in the underlying activation patterns, whether the effects are multidimensional or unidimensional, or any other substantive differences. Summary of the invention

[0007] To solve the above technical problems, the present invention provides a multivariate neural decoding method based on OPM-MEG and RSA. OPM-MEG can provide high-temporal and spatial resolution data on brain activity and provide rich samples of neuronal pattern information, while RSA can capture the subtle interactions between different neurons in the brain, analyze the patterns in these data, and reveal the representation of brain activity at different stages. Therefore, the present invention can combine high-resolution neurophysiological data and high-dimensional analysis methods to understand the process of neural representation more finely and comprehensively.

[0008] To achieve the above object, the present invention adopts the following technical solutions:

[0009] A multivariate neural decoding method based on OPM-MEG and RSA, comprising the following steps:

[0010] Step S1, acquiring resting-state structural MRI data and performing segmentation processing;

[0011] Step S2, scanning the subject wearing the OPM-MEG sensor helmet using a structured light scanner to obtain a scanned image and perform co-registration with the segmented MRI data;

[0012] Step S3, playing auditory stimulation and completing OPM-MEG data acquisition;

[0013] Step S4, preprocessing the collected OPM-MEG data to obtain segmented data;

[0014] Step S5, averaging different trials of the same auditory stimulus at the segmented data level and constructing a neural response pattern vector, and calculating the Pearson correlation of all possible pairwise stimulus pairs at each time point;

[0015] Step S6, obtaining a representation similarity matrix in the time dimension according to the result of step S5, averaging the non-diagonal elements of the representation similarity matrix, performing the operation of step S6 on all subjects, and obtaining an average result of the similarity of the neural response patterns;

[0016] Step S7, performing a non-parametric Wilcoxon rank sum test at each time point to determine the time points with differences;

[0017] Step S8, estimating the source activation time series of each subject through MNE-Python;

[0018] Step S9, projecting the result of step S8 onto a common brain template and averaging the results among all subjects, and using a T-test to determine the activated area.

[0019] In a second aspect, the present invention provides an electronic device, comprising: one or more processors; a memory for storing one or more programs; wherein, when the one or more programs are executed by the one or more processors, the one or more processors implement the aforementioned multivariate neural decoding method based on OPM-MEG and RSA.

[0020] In a third aspect, the present invention provides a computer-readable storage medium having executable instructions stored thereon, which, when executed by a processor, enables the processor to implement the aforementioned multivariate neural decoding method based on OPM-MEG and RSA.

[0021] The beneficial effects of the present invention are:

[0022] (1) Direct neurophysiological data with high spatial and temporal resolution: On the one hand, the OPM-MEG used in the present invention can simultaneously measure neuronal activities at multiple locations in the brain with high spatial and temporal resolution, providing rich samples of neuronal pattern information and more finely capturing the information and processes represented by a large number of neuronal groups. Even small angular differences between adjacent dipoles will produce separable field patterns in MEG statistics. At the same time, the propagation of the magnetic field is less affected by the skull and tissue, so MEG has higher positioning accuracy. On the other hand, OPM-MEG directly measures the magnetic activity of neurons and can provide direct observation of brain activity.

[0023] (2) Comprehensive consideration of the activity patterns of a large number of neuronal groups: The RSA adopted in the present invention can simultaneously utilize the spatial and temporal richness of multi-channel data and comprehensively consider multiple brain regions, rather than being limited to the activities of a single brain region. In the representation space, the neural response to each stimulus is represented as a high-dimensional neural response pattern vector, which can reveal the brain's decoding representation process by calculating the distance of neural response patterns in the representation space, thereby providing a more comprehensive understanding of the brain's response patterns under specific tasks or conditions.

[0024] (3) Dynamic evolution of neural decoding: The present invention combines OPM-MEG with RSA, which can not only locate brain areas with different response characteristics to stimuli, but also understand the process-level information contained in these brain areas and reveal the continuous changes in neural response patterns at different stages of the decoding process. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] Figure 1 This is a flow chart of a multivariate neural decoding method based on OPM-MEG and RSA of the present invention;

[0026] Figure 2 A schematic diagram of a similarity matrix representing an embodiment of the present invention;

[0027] Figure 3 This is a schematic diagram of the similarity results of neural response patterns according to an embodiment of the present invention;

[0028] Figure 4 It is a schematic diagram of source localization results of pattern similarity changes according to an embodiment of the present invention. DETAILED DESCRIPTION

[0029] The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

[0030] The present invention provides a multivariate neural decoding method based on OPM-MEG and RSA, such as Figure 1 As shown, the specific steps are as follows:

[0031] Step S1: Acquire resting structural MRI data and perform segmentation processing. The specific implementation steps are as follows:

[0032] Step S1-1, setting MRI scanning parameters through a magnetic resonance imaging console;

[0033] Step S1-2, the subject lies on the scanning bed and fixes the head;

[0034] Step S1-3, monitoring the subject's state in real time through the built-in monitoring system of the MRI device, and collecting MRI scan image data in a static state;

[0035] Step S1-4: segment the acquired MRI scan image data to obtain an anatomical label map containing segmentation, surface reconstruction, and morphological measurement information of various brain structures.

[0036] Step S2: Scan the subject wearing the OPM-MEG sensor helmet based on a structured light scanner to obtain a scanned image and perform co-registration with the segmented MRI data. The specific implementation steps are as follows:

[0037] Step S2-1, after the subject wears the rigid helmet with the OPM-MEG sensor, a structured light scanner is used to scan the subject's head to obtain a point cloud image of the relative position between the subject's head and the OPM-MEG sensor;

[0038] Step S2-2, using a region growing method to separate the helmet point cloud and the face point cloud obtained in the point cloud image;

[0039] Step S2-3, registering the rigid helmet model and the helmet point cloud based on a random sampling consistency initial registration algorithm;

[0040] Step S2-4, locating the nose tip region in the scanned head point cloud obtained in S2-1 and the nose tip region in the MRI structural image obtained in S1-4 and registering them using a random sampling consistency initial registration algorithm;

[0041] Step S2-5: Based on the iterative closest point algorithm, the registration results obtained in steps S2-3 and S2-4 are precisely registered to minimize the sum of squared Euclidean distances between the nose tip area corresponding to the structured light scan and the MRI structure image. The calculation formula is as follows;

[0042] ,

[0043] in, and is a pair of corresponding points between the MRI structural image and the structured light scanning point cloud, with a total of For the corresponding points, is a rotation matrix used to describe the rotation transformation of MRI structural images and structured light scanning point clouds in three-dimensional space, and T is a translation matrix used to describe the translation transformation of MRI structural images and structured light scanning point clouds in three-dimensional space;

[0044] Step S3: Play auditory stimulation and complete OPM-MEG data collection. The specific implementation steps are as follows:

[0045] Step S3-1, the subject is provided with a helmet having an OPM-MEG sensor placed thereon and adjusted to a suitable position;

[0046] Step S3-2: Play auditory stimulation and collect data.

[0047] Step S4: pre-process the collected OPM-MEG data to obtain segmented data. The specific implementation steps are as follows:

[0048] Step S4-1, perform 1-30 Hz offline bandpass filtering on the collected OPM-MEG data, plot the power spectrum density, and check and delete bad tracks based on the power spectrum density;

[0049] Step S4-2, segmenting the data according to the recorded stimulation synchronization trigger signal, extracting the data of 100 ms before the baseline and 1000 ms after the baseline;

[0050] Step S4-3, performing artifact removal on the data based on independent component analysis;

[0051] Step S4-4, perform visual inspection, judge and remove abnormal segmented data;

[0052] Step S4-5: The segmented OPM-MEG signal obtained after processing in the above steps is as follows:

[0053] ,

[0054] in, is the basic EEG signal, t is the time, is the kth component of the basic EEG signal, Q is the number of components, A Based on magnetic brain signal The instantaneous amplitude of the kth component, Based on magnetic brain signal The instantaneous phase of the kth component, e is a constant, and j is an imaginary unit.

[0055] Step S5: average different trials of the same auditory stimulus at the segmented data level and construct a neural response pattern vector, and calculate the Pearson correlation of all possible pairwise stimulus pairs at each time point. The specific implementation steps are as follows:

[0056] Step S5-1, averaging different trials of the same sound stimulation to obtain a neural response pattern vector corresponding to the sound stimulation, where the dimension of the neural response pattern vector is the number of OPM-MEG sensors;

[0057] Step S5-2, calculate the Pearson correlation R between the two sound stimulus pairs, where an R value of 1 represents that the neural response patterns under the two conditions are completely correlated, and an R value of 0 represents that the neural response patterns under the two conditions are completely uncorrelated:

[0058] ,

[0059] in, and are the brain nerve signal strengths measured by all sensors at a certain time t before and after the baseline, that is, the observed values ​​of the neural response pattern vector, and are the average values ​​of brain neural signal strength measured by all sensors at all times before and after the baseline, respectively. The numerator is the covariance, which measures whether the pattern distribution before and after the baseline tends to be consistent, and the denominator is the product of the standard deviation for normalization;

[0060] Step S5-3: Repeat steps S5-1 and S5-2 for all possible stimulus pairs, and obtain the correlation of all stimulus pairs at the current moment.

[0061] Step S6: According to the result of step S5, a representation similarity matrix in the time dimension is obtained, and the non-diagonal elements of the representation similarity matrix are averaged within the subject and between the subjects to obtain the average result of the similarity of the neural response pattern. The specific implementation steps are as follows:

[0062] Step S6-1, using the correlation of all stimulus pairs at the current moment in step S5-3 to form a representation similarity matrix at the current moment, wherein the representation similarity matrix is ​​a square matrix symmetric about the diagonal, the dimension of which is the number of sound stimuli, and each row or column represents the similarity between one stimulus and all stimuli;

[0063] Step S6-2, repeating step S5-3 for each subject at each time point, to obtain a series of representation similarity matrices for each subject along the time axis;

[0064] Step S6-3, average the non-diagonal elements in the lower left corner of each similarity matrix, calculate the average Pearson correlation R value at each time point, repeat the non-diagonal element averaging operation for all subjects, average these averages across subjects, and obtain a time series of neural response pattern similarity.

[0065] Step S7: Perform non-parametric Wilcoxon rank sum test at each time point to determine the time points with differences. The specific implementation steps are as follows:

[0066] Step S7-1, calculating the average neural response pattern similarity of each subject 100 ms before the baseline;

[0067] Step S7-2, combining the average neural response pattern similarity of each subject at a certain moment after the baseline with the average neural response pattern similarity results 100 ms before the baseline into a whole;

[0068] Step S7-3, sorting the merged whole according to the value from small to large;

[0069] Step S7-4, assigning a rank r to each sorted neural response pattern similarity;

[0070] Step S7-5: Use the U statistic as the statistic for the Wilcoxon rank sum test, and the calculation formula is:

[0071] ,

[0072] ,

[0073] ,

[0074] in, and are the data sample sizes at a certain moment before and after the baseline, and Respectively represent the rank of each sample after merging and sorting at a certain moment before and after the baseline;

[0075] Step S7-6: Calculate the rank sum of the pre-baseline and post-baseline data respectively to obtain the Wilcoxon rank sum statistic U, and obtain the corresponding significance level according to the Wilcoxon rank sum test distribution table. value;

[0076] Step S7-7, set the difference level to 0.05, and search the distribution table to find the corresponding rank sum critical value;

[0077] Step S7-8, if the rank sum obtained in S7-6 is greater than the rank sum critical value in S7-7, it is considered that the similarity of the neural response pattern at this moment is different from the baseline, otherwise it is considered that there is no difference;

[0078] Step S7-9, repeat step S7-8 at each time point to obtain a series of time points with differences.

[0079] Step S8: Estimate the source activation time series of each subject through MNE-Python. The specific implementation steps are as follows:

[0080] Step S8-1, obtaining the position of the subject's OPM-MEG sensor and the corresponding brain model data from step S2-1;

[0081] Step S8-2, obtaining the subject source space information from step S1-4;

[0082] Step S8-3, calculating a forward solution by combining the position of the OPM-MEG sensor, the subject's brain model data, and the subject's source space information to obtain the contribution of the current source in the brain to the scalp magnetic field;

[0083] Step S8-4: Based on the dynamic statistical parameter mapping algorithm, estimate the source space activity of the significant time interval obtained in step S7-9, and standardize the result into a Z value.

[0084] Step S9: Project the result of step S8 onto the common brain template and average it among the subjects, and use T test to determine the activated area. The specific implementation steps are as follows:

[0085] Step S9-1, projecting the source activity Z values ​​of all subjects in step S8-4 onto the common brain template and averaging them;

[0086] Step S9-2: For the Z value obtained in step S9-1, use the T test to calculate the statistical value of each vertex and time point and further obtain the corresponding test significance level. value;

[0087] Step S9-3: Use FDR correction to avoid multiple comparison problems and adjust the test significance level Correction is performed and the significance level of the correction test is The value replaces the original data value;

[0088] Step S9-4: The values ​​were projected onto the common brain template and the source activation distribution of the time interval obtained in step S7-9 was viewed based on the MNE visualization tool.

[0089] Figure 2 The schematic result of representing the similarity matrix is ​​shown, where N represents the number of sound stimuli, so the dimension of the matrix is ​​N. 2 Each value in the matrix represents the correlation between the sound stimulus pairs in the corresponding row and column. This is a diagonally symmetric square matrix. The closer the value is to 1, the more correlated the response patterns of the two stimuli are. The closer the value is to 0, the less correlated the response patterns of the two stimuli are. Each row or column represents the correlation between the stimulus and all stimuli.

[0090] Figure 3 The figure shows the changes in the similarity of neural response patterns during auditory processing. The curve represents the average similarity value at each time point, reflecting the changes in the similarity of neural response patterns over time. The light shade represents the standard deviation, and the horizontal line above the curve marks the time points that are different from the baseline.

[0091] Figure 4 The sources of the above-mentioned pattern similarity changes were displayed from the lateral, medial and ventral perspectives. It can be seen that the bilateral temporal lobes and parahippocampal gyrus of the brain contribute to the changes in the similarity of neural response patterns.

[0092] In a second aspect, the present invention provides an electronic device, comprising: one or more processors; a memory for storing one or more programs; wherein, when the one or more programs are executed by the one or more processors, the one or more processors implement the aforementioned multivariate neural decoding method based on OPM-MEG and RSA.

[0093] In a third aspect, the present invention provides a computer-readable storage medium having executable instructions stored thereon, which, when executed by a processor, enables the processor to implement the aforementioned multivariate neural decoding method based on OPM-MEG and RSA.

[0094] The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A multivariate neural decoding method based on OPM-MEG and RSA, characterized in that: The following steps are involved: Step S1, acquiring resting-state structural MRI data and performing segmentation processing; Step S2, scanning the subject wearing the OPM-MEG sensor helmet using a structured light scanner to obtain a scanned image and perform co-registration with the segmented MRI data; Step S3, playing auditory stimulation and completing OPM-MEG data acquisition; Step S4, preprocessing the collected OPM-MEG data to obtain segmented data; Step S5, averaging different trials of the same auditory stimulus at the segmented data level and constructing a neural response pattern vector, and calculating the Pearson correlation of all possible pairwise stimulus pairs at each time point; Step S6, obtaining a representation similarity matrix in the time dimension according to the result of step S5, averaging the non-diagonal elements of the representation similarity matrix, performing the operation of step S6 on all subjects, and obtaining an average result of the similarity of the neural response patterns; Step S7, performing a non-parametric Wilcoxon rank sum test at each time point to determine the time points with differences; Step S8, estimating the source activation time series of each subject through MNE-Python; Step S9, projecting the result of step S8 onto a common brain template and averaging the results among all subjects, and using a T-test to determine the activated area.

2. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 1, characterized in that: The step S1 comprises: Step S1-1, setting MRI scanning parameters through a magnetic resonance imaging console; Step S1-2, the subject lies on the scanning bed and fixes the head; Step S1-3, monitoring the subject's state in real time through the built-in monitoring system of the MRI device, and collecting MRI scan image data in a static state; Step S1-4: segment the acquired MRI scan image data to obtain an anatomical label map containing segmentation, surface reconstruction, and morphological measurement information of various brain structures.

3. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 1, characterized in that: The step S2 comprises: Step S2-1, after the subject wears the rigid helmet with the OPM-MEG sensor, a structured light scanner is used to scan the subject's head to obtain a point cloud image of the relative position between the subject's head and the OPM-MEG sensor; Step S2-2, using a region growing method to separate the helmet point cloud and the face point cloud obtained in the point cloud image; Step S2-3, registering the rigid helmet model and the helmet point cloud based on a random sampling consistency initial registration algorithm; Step S2-4, locating the nose tip area in the point cloud image obtained in step S2-1 and the nose tip area in the MRI scan image obtained in step S1-4, and aligning them using a random sampling consistency initial alignment algorithm; Step S2-5, based on the iterative closest point algorithm, the registration results obtained in step S2-3 and step S2-4 are precisely registered to minimize the sum of squares of the Euclidean distances between the point cloud image and the nose tip area corresponding to the MRI scan image; ， in, and is a pair of corresponding points between the MRI scan image and the point cloud image, with a total of For the corresponding points, represents the rotation matrix and T represents the translation matrix.

4. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 1, characterized in that: The S4 includes: Step S4-1, perform 1-30 Hz offline bandpass filtering on the collected OPM-MEG data, plot the power spectrum density, and check and delete bad tracks based on the power spectrum density; Step S4-2, segmenting the data according to the recorded auditory stimulation synchronization trigger signal, extracting the data of 100 ms before the baseline and 1000 ms after the baseline; Step S4-3, performing artifact removal based on independent component analysis on the data extracted in step S4-2; Step S4-4: Perform visual inspection, judge and remove abnormal segmented data, and obtain segmented OPM-MEG signal data.

5. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 4, characterized in that: The segmented OPM-MEG signal is: ， in, is the basic EEG signal, t is the time, is the kth component of the basic EEG signal, Q is the number of components, A Based on magnetic brain signal The instantaneous amplitude of the kth component, Based on magnetic brain signal The instantaneous phase of the kth component, e is a constant, and j is an imaginary unit.

6. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 1, characterized in that: The S5 includes: Step S5-1, averaging different trials of the same sound stimulation to obtain a neural response pattern vector corresponding to the sound stimulation, where the dimension of the neural response pattern vector is the number of OPM-MEG sensors; Step S5-2, calculate the Pearson correlation R between the two sound stimulus pairs, where an R value of 1 represents that the neural response patterns under the two conditions are completely correlated, and an R value of 0 represents that the neural response patterns under the two conditions are completely uncorrelated: ， in, and are the brain nerve signal strengths measured by all sensors at time t before and after the baseline, i.e., the observed values ​​of the neural response pattern vector, and They are the average values ​​of brain neural signal strength measured by all OPM-MEG sensors at all times before and after the baseline; Step S5-3: Repeat steps S5-1 and S5-2 for all possible pairs of sound stimuli to obtain the correlation of all stimulus pairs at the current moment.

7. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 1, characterized in that: The S6 includes: Step S6-1, using the correlation of all stimulus pairs at the current moment in step S5-3 to form a representation similarity matrix at the current moment, wherein the representation similarity matrix is ​​a square matrix symmetrical about the diagonal, the dimension of which is the number of sound stimuli, and each row or column represents the similarity between one stimulus and all stimuli; Step S6-2, repeating step S5-3 for each subject over time, to obtain a series of representation similarity matrices for each subject along the time axis; Step S6-3, average the non-diagonal elements in the lower left corner of each similarity matrix, calculate the average Pearson correlation R value at each time point, repeat the non-diagonal element averaging operation for all subjects, average the obtained average values ​​among the subjects, and obtain the time series of neural response pattern similarity.

8. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 1, characterized in that: The S7 includes: Step S7-1, calculating the average neural response pattern similarity of each subject 100 ms before the baseline; Step S7-2, combining the average neural response pattern similarity of each subject at a certain moment after the baseline with the average neural response pattern similarity 100 ms before the baseline; Step S7-3, sorting the merged whole according to the value from small to large; Step S7-4, assigning a rank r to each sorted neural response pattern similarity; Step S7-5: Use the U statistic as the statistic for the Wilcoxon rank sum test, and the calculation formula is: ， ， ， in, and are the data sample sizes at a certain moment before and after the baseline, and Respectively represent the rank of each sample after merging and sorting at a certain moment before and after the baseline; Step S7-6: Calculate the rank sum of the pre-baseline and post-baseline data respectively to obtain the Wilcoxon rank sum statistic U, and obtain the corresponding significance level according to the Wilcoxon rank sum test distribution table. value; Step S7-7, set the difference level to 0.05, and search the distribution table to find the corresponding rank and critical value; Step S7-8, if the rank sum obtained in step S7-6 is greater than the rank sum critical value in step S7-7, it is considered that the similarity of the neural response pattern at this moment is different from the baseline, otherwise it is considered that there is no difference; Step S7-9, repeat step S7-8 at each time point to obtain a series of time points with significant differences.

9. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 1, characterized in that: The S8 includes: Step S8-1, obtaining the position of the subject's OPM-MEG sensor and the corresponding brain model data from step S2-1; Step S8-2, obtaining the subject source space information from step S1-4; Step S8-3, calculating a forward solution by combining the position of the OPM-MEG sensor, the subject's brain model data, and the subject's source space information to obtain the contribution of the current source in the brain to the scalp magnetic field; Step S8-4: Based on the dynamic statistical parameter mapping algorithm, estimate the source space activity of the significant difference time point obtained in step S7-9, and standardize the result into a Z value.

10. The multivariate neural decoding method based on OPM-MEG and RSA according to claim 1, characterized in that: The S9 includes: Step S9-1, projecting the Z values ​​corresponding to the source activities of all subjects in S8-4 onto the common brain template and averaging them; Step S9-2: Use T test to calculate the statistical value of each vertex and time point and obtain the corresponding test significance level value; Step S9-3: Test significance level The corrected test significance level was obtained after the value was corrected value, replacing the significance level in the original data value; Step S9-4: Set the corrected test significance level of step S9-3 to The values ​​were projected onto the common brain template and the spatial distribution of source activations at the time points obtained in steps S7-9 was viewed using the MNE-based visualization tool.

11. An electronic device, characterized in that: include: one or more processors; A memory for storing one or more programs; Wherein, when one or more programs are executed by the one or more processors, the one or more processors implement the multivariate neural decoding method based on OPM-MEG and RSA as described in any one of claims 1-10.

12. A computer-readable storage medium, characterized in that: Executable instructions are stored thereon, which, when executed by a processor, enable the processor to implement the multivariate neural decoding method based on OPM-MEG and RSA as described in any one of claims 1-10.

Images