[go: up one dir, main page]

CN112168167A - Method and electronic device for detecting single auditory evoked signal in magnetoencephalography - Google Patents

Method and electronic device for detecting single auditory evoked signal in magnetoencephalography Download PDF

Info

Publication number
CN112168167A
CN112168167A CN202011026442.5A CN202011026442A CN112168167A CN 112168167 A CN112168167 A CN 112168167A CN 202011026442 A CN202011026442 A CN 202011026442A CN 112168167 A CN112168167 A CN 112168167A
Authority
CN
China
Prior art keywords
signal
sensor
sample
signals
brain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011026442.5A
Other languages
Chinese (zh)
Other versions
CN112168167B (en
Inventor
郭弘
吴腾
彭翔
张建玮
冯雨龙
肖伟
孙晨曦
吴玉龙
张相志
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking University
Original Assignee
Peking University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Peking University filed Critical Peking University
Priority to CN202011026442.5A priority Critical patent/CN112168167B/en
Publication of CN112168167A publication Critical patent/CN112168167A/en
Application granted granted Critical
Publication of CN112168167B publication Critical patent/CN112168167B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Physiology (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明提供提供了一种脑磁图单次听觉诱发信号检测方法及电子装置,包括:将待检测脑磁信号转换为一个二维矩阵数据,依次以各脑磁传感器为选定传感器,计算选定传感器与设定邻域内各传感器的相关系数,并根据选定传感器二维矩阵数据、设定邻域内各脑磁传感器的二维矩阵数据及相关系数,计算各选定传感器的增加信号;根据增加信号及各脑磁传感器的位置,得到脑磁信号视图;提取脑磁信号视图的信号空间分布特征,并对所述信号空间分布特征进行分类,判断待检测信号中包含的单次听觉诱发信号或噪声信号。本发明能够最大程度上保留原始信号的信息,而不会影响信号中包含的其余的信息,检测更加准确,计算也更加简便。

Figure 202011026442

The present invention provides a method and an electronic device for detecting a single auditory evoked signal of a magnetoencephalogram, comprising: converting the magnetoencephalography signal to be detected into a two-dimensional matrix data, sequentially taking each magnetoencephalography sensor as a selected sensor, and calculating the selected sensor. Determine the correlation coefficient between the sensor and each sensor in the set neighborhood, and calculate the increase signal of each selected sensor according to the two-dimensional matrix data of the selected sensor, the two-dimensional matrix data and the correlation coefficient of each brain magnetism sensor in the set neighborhood; Add the signal and the position of each EEG sensor to obtain a view of the EEG signal; extract the signal spatial distribution features of the EEG view, classify the signal spatial distribution features, and determine the single auditory evoked signal contained in the signal to be detected or noise signal. The present invention can retain the information of the original signal to the greatest extent without affecting the rest of the information contained in the signal, so that the detection is more accurate and the calculation is simpler.

Figure 202011026442

Description

Magnetoencephalogram single auditory evoked signal detection method and electronic device
Technical Field
The invention relates to the field of magnetoencephalogram feature signal identification in the field of biological feature identification, in particular to a magnetoencephalogram single auditory evoked signal automatic detection method and an electronic device.
Background
The magnetoencephalogram can detect weak magnetic signals emitted by the human brain by using an extremely sensitive magnetic detector, such as an atomic optical pump magnetometer, a superconducting quantum interferometer and the like, extract corresponding information from the signals and analyze the brain state of a tested person. The current application of magnetoencephalography is mainly two-fold: one aspect is the partitioning and study of different functional areas of the brain network. On the other hand, the compound has important clinical application in the diagnosis and treatment of some brain functional diseases, such as psychosis, brain injury and autism. In particular, in epileptic conditions, the magnetoencephalogram can approximate the location of an epileptic focus, thereby allowing noninvasive and non-invasive diagnosis of the patient.
In the study of human brain function using magnetoencephalogram, the main study object is the signals emitted by the brain when the subject is subjected to different stimuli, and these signals are called event-induced signals. For example, when the testee is stimulated by different visual, auditory, sensory and olfactory senses, different functional areas of the brain will send out corresponding event inducing signals. However, since brain magnetic signals are very weak, and various noises exist. If the time-induced signal is observed, the testee needs to be stimulated for multiple times, the event-induced signal is measured for multiple times, then the event-induced signal obtained through multiple measurement is stored for signal averaging, so that the noise signal which randomly appears in the event-induced signal is suppressed, and the regular signal pattern is strengthened. On the one hand, however, the method of averaging signals yields results that depend on the common trends and "components" present in the observed waveforms. For each stimulation, the waveform, the appearance time, the appearance position, the frequency and other characteristics of each corresponding evoked signal are changed due to the differences of the attention, the tension state and the excitation state of the testee. Thus, part of the information of the event-induced signal in a single measurement is lost in averaging the signal. In addition, in a single evoked signal, there are a number of other brain region-emanating signals as well as ocular and electrocardiographic artifact signals, which also occur randomly, playing a role as noise in the detection of event-evoked signals. Averaging also results in loss of information for these signals. On the other hand, if there is significant noise in a single measurement, the averaging operation also introduces noise into the averaged observed signal. The averaged signal from the averaging operation can only provide a pattern of event-induced signals that is not representative of any single event-induced signal.
Currently, methods for detecting single-event evoked signals are mainly proposed based on electroencephalography, such as separating feature quantities from single-event evoked signals by using independent component analysis, signal space projection, and the like. But the characteristics of such isolation are not obvious due to the particularity and diversity of the brain signal sources. With the improvement of the detection technology, the waveform of the signal naturally becomes the entry point of the next single event induced signal detection technology, and people construct statistics by using the frequency characteristics, the amplitude characteristics and the phase characteristics of the single event induced signal appearing in the electroencephalogram and detect the single event induced signal by using a statistical method. However, there is a difference between the magnetoencephalogram and the event-induced signal features detected in the electroencephalogram, and the useful information about the signal is not the same. Especially, the magnetoencephalogram generally has more sensors due to the positioning requirement of a signal source, and the spatial distribution of detection signals can be obtained. While the electroencephalogram measures the potential difference between the two sensors. Therefore, based on the above reasons, we propose a new algorithm for detecting single event evoked signals in magnetoencephalogram, especially single auditory evoked signals.
Disclosure of Invention
Aiming at the defects of the traditional single event evoked signal detection method in the magnetoencephalogram, the invention provides a magnetoencephalogram single auditory evoked signal detection method and an electronic device.
The technical scheme of the invention is as follows:
a method for detecting a single auditory evoked signal of a magnetoencephalogram comprises the following steps:
1) converting the brain magnetic signals to be detected into two-dimensional matrix data of M x N, wherein M is the number of the brain magnetic sensors, and N is a set time length;
2) sequentially taking each brain magnetic sensor as a selected sensor, calculating the correlation coefficient of the selected sensor and each sensor in a set neighborhood, and calculating an increasing signal of each selected sensor according to the two-dimensional matrix data of the selected sensor, the two-dimensional matrix data of each brain magnetic sensor in the set neighborhood and the correlation coefficient;
3) obtaining a brain magnetic signal view according to the added signals and the positions of the brain magnetic sensors;
4) and extracting the signal space distribution characteristics of the brain magnetic signal view, classifying the signal space distribution characteristics, and judging a single auditory evoked signal or a noise signal contained in the signal to be detected.
Further, the increase signal for each selected sensor is calculated by the following strategy:
1) selecting the sensor X in the corresponding set neighborhood with the correlation coefficient larger than the threshold cjSignal GjAnd correspondingly setting the sensor X in the neighborhoodjSignal G ofjBy a correlation coefficient TjFor weighting, superimposed on the selected sensor XiSignal G ofiThe obtained signal
Figure BDA0002702249320000021
As selected sensor XiNew detection signals, wherein I is more than or equal to 1 and less than or equal to I, I is the number of the selected sensors, J is more than or equal to 1 and less than or equal to J, and J is the number of the sensors in the set neighborhood with the correlation coefficient larger than the threshold value c;
2) if none of the correlation coefficients is greater than the threshold c, the increase signal G for the selected sensor is selectedi'=Gi
Further, the brain magnetic signal view is a 2D view; obtaining the brain magnetic signal view by the following steps:
1) calculating the average power of each brain magnetic sensor strengthening signal;
2) and coding the average power by using colors, and drawing the average power to a spatial position corresponding to the brain magnetic sensor to obtain a brain magnetic signal view.
Further, extracting signal space distribution characteristics of the brain magnetic signal view through a trained two-dimensional deep convolution neural network, and classifying the signal space distribution characteristics; wherein a two-dimensional deep convolutional neural network is trained by:
a) respectively converting a plurality of single auditory evoked sample signals and sample noise signals into sample data with labels;
b) sequentially taking each magnetoencephalography sensor as a selected sensor, calculating sample correlation coefficients of the selected sensor and each magnetoencephalography sensor in a set sample neighborhood, and calculating a sample increase signal of each selected sensor according to the sample data of the selected sensor, the sample data of each magnetoencephalography sensor in the set sample neighborhood and the sample correlation coefficients;
c) adding signals and the positions of the brain magnetic sensors according to the samples to obtain a plurality of single auditory evoked sample brain magnetic signal views or sample noise views;
d) and training the two-dimensional deep convolutional neural network through an auditory evoked sample brain magnetic signal view and a sample noise view to obtain the trained two-dimensional deep convolutional neural network.
Further, the two-dimensional deep convolutional neural network comprises a google lenet network, an AlexNet network or a two-dimensional deep convolutional neural network with the node number and the weight meeting the requirements.
Further, preprocessing the single auditory evoked sample signal before converting the single auditory evoked sample signal into sample data with a label; the pretreatment comprises the following steps: removing ambient noise signals from the single auditory evoked sample signal, removing eye movement and cardiac artifact noise signals, and removing other signals emanating from the brain.
Further, the ambient noise signal includes: power frequency noise; the method for removing power frequency noise comprises the following steps: filtering removal at 50Hz or 60Hz is used.
Further, the method for removing the eye movement and the cardiac artifact noise signal comprises the following steps: a signal space projection algorithm or an independent component analysis algorithm is used.
Further, other signals emitted by the brain include: an alpha signal and a beta signal; the method for removing the alpha signal and the beta signal comprises the following steps: filtering removal at 10Hz and 20Hz was used.
A storage medium having a computer program stored therein, wherein the computer program is arranged to perform the above-mentioned method when executed.
An electronic device comprising a memory having a computer program stored therein and a processor arranged to run the computer to perform the method as described above.
The invention has the beneficial effects that:
a new algorithm for single auditory evoked signal detection in magnetoencephalography based on artificial neural networks is proposed;
compared with the traditional signal averaging method, the single auditory evoked signal detection does not need to store data, and can be realized synchronously with data detection, namely in real time;
a single auditory evoked signal detection is able to preserve the information of the original signal to the maximum extent possible, without affecting the rest of the information contained in the signal, compared to traditional methods of signal averaging;
compared with the traditional algorithm for inducing signals by single events in electroencephalogram, the algorithm utilizes the signal space distribution information contained in the electroencephalogram signals, so that the detection is more accurate, and the calculation is simpler and more convenient.
Drawings
FIG. 1 is a flow chart of an algorithm for automatic detection of a single-measure auditory evoked signal in accordance with an embodiment of the present invention.
FIG. 2 is a flow chart of a single-measurement auditory evoked signal enhancement algorithm in an embodiment of the present invention.
Fig. 3 is a 2D view rendered with a raw acquired single auditory evoked signal and an enhanced single auditory evoked signal in an embodiment of the present invention.
Fig. 4 is a confusion matrix for a single auditory evoked signal test on training set data and test set data using google lenet in an embodiment of the present invention.
Detailed Description
In order to better explain the technical scheme of the invention, the invention is further described in detail with reference to the accompanying drawings and specific embodiments. In order to solve the existing technical problems, the invention provides a single auditory evoked signal detection method in a magnetoencephalogram based on an artificial neural network, which specifically comprises the following steps:
the invention relates to a method for automatically detecting single auditory evoked signals in a magnetoencephalogram, which comprises the following steps as shown in figure 1:
1) the known single auditory evoked signals are signal pre-processed in order to remove noise therefrom. Wherein each single auditory evoked signal is a data set in the form of a two-dimensional matrix with a size of M N, M being the number of magnetoencephalography sensors; and N is a set time length.
2) Calculating a correlation coefficient for the signal preprocessed in step 1). The neighborhood radius size is defined as r. For each sensor, a neighborhood is defined, and for ease of description, the sensor is referred to as the selected sensor, and the sensors within its neighborhood are referred to as neighborhood sensors. And respectively calculating a correlation coefficient between each neighborhood sensor detection signal and the selected sensor detection signal, namely a neighborhood correlation coefficient.
3) As shown in fig. 2, the signal preprocessed in step 1) is subjected to linear superposition enhancement. When the correlation coefficient between the neighborhood signal and the selected signal calculated in the step 2) is larger than the threshold c, weighted linear superposition is carried out on the signals of the neighborhood signal and the selected signal, then the superposed signal is divided by the number of the superposed neighborhood signals to serve as a new selected signal, and the weight is the correlation coefficient between the neighborhood signal and the selected signal. But if the number of the superimposed domain signals is 0, no operation is performed on the selected signal.
4) And drawing a 2D view by using the enhanced signal obtained in the step 3). The average power of the enhancement signal is calculated, color coded, and combined with the sensor position information, a 2D view can be rendered.
5) The spatial distribution enhancement operation described above is performed on a known single auditory evoked signal. A 2D view of a single auditory evoked signal is obtained as training set 1. Wherein each single auditory evoked signal is a data set in the form of a two-dimensional matrix with a size of M N, M being the number of magnetoencephalography sensors; and N is a set time length.
6) The spatial distribution enhancement operation described above is performed on a known noise signal. A 2D view of the noise signal is obtained as training set 2. Each noise signal is a data set in a two-dimensional matrix form with the size of M x N, and M is the number of the brain magnetic sensors; and N is a set time length. The measurement environment of the noise signal is the same as the single auditory evoked signal measurement environment.
7) Training a single auditory evoked signal detection model by using the training set 1 and the training set 2 obtained in the steps 5) and 6) to form a training set total; the single auditory evoked signal detection model comprises a feature learning module and a classification output module; the feature learning module comprises a two-dimensional deep convolutional neural network and is used for learning the spatial distribution features of signals presented in a training set, namely a 2D view, extracting the spatial distribution features of the signals and constructing a corresponding deep convolutional neural network model; and then, the classification output module divides the training set data into auditory evoked signals and noise signals according to the spatial distribution characteristics of the single auditory evoked signals extracted from the characteristic learning module.
8) The above-described spatial distribution enhancement operation is performed on the unknown signal to obtain a 2D view thereof. Each unknown signal is a data set in a two-dimensional matrix form with the size of M x N, and M is the number of the brain magnetic sensors; and N is a set time length.
9) And respectively inputting the 2D views of the unknown signals into the trained single auditory evoked signal detection model to obtain corresponding auditory evoked signal detection results.
10) And determining which data are auditory evoked signals and which data are noise signals in the unknown magnetoencephalogram signals according to the auditory evoked signals detection result obtained in the step 9).
Further, the noise includes some environmental noise signals, eye movement and cardiac artifact noise signals, and other signals emitted by the brain.
Further, the environmental noise signal, such as power frequency noise, can be removed by filtering at 50Hz or 60 Hz.
Further, the eye movement and the cardiac artifact signals can be removed by a signal space projection algorithm or an independent component analysis method.
Further, other signals emitted from the brain, such as the alpha signal and the beta signal, can be removed by filtering at 10Hz and 20 Hz.
Further, the neighborhood radius r should be dependent on the device that measures the magnetoencephalogram.
Further, the correlation coefficient threshold c should be determined according to the apparatus for measuring the magnetoencephalogram.
Further, the spatial resolution of the 2D view is related to the number of sensors that the apparatus for measuring the magnetoencephalogram has.
The following is an application example of the present invention:
first there is a known single measurement of the auditory evoked signal. The example signals used in this embodiment are brain magnetic signals measured by SQUID-MEG equipment of CTF, canada, which has 274 sensors, so 274 sets of signals are measured. The measurement of a single auditory evoked signal lasts approximately 0.3 seconds, with a sampling rate of 600Hz, and thus a data length of 0.3 x 600. The single auditory evoked signal fragment is a data matrix of 274 x 180. A single measurement of the auditory evoked signal is employed. There are 197 such single auditory evoked signals in total.
There is also a known noise signal, which should be the same duration as the auditory evoked signal, and a data matrix with data size of 274 x 180. And measured in the same environment as the single auditory evoked signal. There are 200 random noise signals.
The known auditory evoked signals and the noise signals are subjected to signal preprocessing. Including 60Hz, 120Hz, 180Hz, 10Hz, 20Hz filtering operations. Filtering at 60Hz and multiples thereof is intended to remove power frequency noise (the segment of data is measured in europe), filtering at 10Hz is intended to remove the alpha signal, and filtering at 20Hz is intended to remove the beta signal. Meanwhile, signal space projection operation is carried out on the eye movement artifact and the heartbeat artifact signals so as to eliminate the influence of the eye movement artifact and the heartbeat artifact on the magnetoencephalogram signals.
The known single auditory evoked signal and the noise signal are subjected to a signal enhancement processing operation. Each of the 274 sensors is first treated as a selected sensor in turn. The neighborhood of the selected sensor is defined by centering on the location of the selected sensor and by 0.022 x 1.7 meters as radius, and the sensors within the neighborhood of the selected sensor may be referred to as neighborhood sensors. The correlation coefficient between the signal detected by the selected sensor and the signal detected by each of the neighboring sensors is calculated separately. If the correlation coefficient is greater than or equal to 0.8, the signals detected by the adjacent sensors are superposed on the signals detected by the selected sensors by the weight of the correlation coefficient. If the correlation coefficient of the two is less than 0.8, no action is taken. And (4) taking the number of the final superposed signals as the number of the field signals, and taking the obtained signals as the signals detected by the final selected sensor. But if the number of the superimposed domain signals is 0, no operation is performed on the selected signal. Fig. 3 is a 2D view rendered with a single acoustic evoked signal from a raw acquisition and a single acoustic evoked signal after enhancement.
The radius of the above-mentioned field is selected to be 0.022 × 1.7 meters, and according to the distribution of the sensors, it can be known that there are 7-10 field sensors in the radius of the neighborhood, and it can be known through statistics that this parameter is relatively in accordance with this algorithm in this embodiment. In addition, the correlation coefficient threshold is also statistically selected to be 0.8 in the present embodiment.
The 2D view is rendered using the single auditory evoked signal and the noise signal subjected to the signal enhancement operation described above. Since the spatial distribution of the sensors is known. The enhancement signal can be color coded according to its spatial distribution, mapped to the corresponding locations, resulting in a 2D view describing the spatial distribution of the signal with a resolution of 224 x 224 pixels. Finally, 197 2D views of the single auditory evoked signal and 200 2D views of the noise signal can be obtained. There are 397 2D views as our training set.
And training the pre-training network GoogLeNet by using the obtained training set. The structure of the pre-training network GoogLeNet is adjusted correspondingly, so that the recognition of single auditory evoked signals can be completed. The main adjustment is as follows: the number of filters of the last fully connected layer is modified to be equal to the number of classifications (i.e., single auditory evoked signal versus noise signal). Meanwhile, the learning rate factor of the full connection layer is increased, so that the data can be converged more quickly; the number of output classes of the designated network in the classification layer is modified to 2. The classification layer is replaced with a new classification layer without class labels. When training the network, it will automatically set the output class of the layer; and modifying the last dropout layer, wherein the dropout layer is mainly used for randomly setting the parameter transmitted from the previous layer to be 0 so as to prevent the network from being overfitted. Here the probability of setting the layer setting element to 0 is adjusted from 0.5 to 0.6. The network is then trained using a training set.
The trained network can be used for detecting single auditory evoked signals of original training data and also can be used for detecting single auditory evoked signals existing in unknown brain magnetic signals. In an embodiment, the trained network is used to detect the original training data, i.e. the single auditory evoked signals in 397 2D views, with a detection success rate of 100%. A single auditory evoked signal in the unknown data is also detected. The unknown data contains 199 single auditory evoked signals and 200 noise signals. The final detection success rate was also 100%. As shown in fig. 4, is a confusion matrix representing the detection results.
The above embodiments are only intended to illustrate the technical solution of the present invention and not to limit the same, and a person skilled in the art can modify the technical solution of the present invention or substitute the same without departing from the spirit and scope of the present invention, and the scope of the present invention should be determined by the claims.

Claims (10)

1.一种脑磁图单次听觉诱发信号检测方法,其步骤包括:1. A magnetoencephalography single-shot auditory evoked signal detection method, the steps comprising: 1)将待检测脑磁信号转换为一M*N的二维矩阵数据,其中M为脑磁传感器数目,N为设定时间长度;1) Convert the magnetic brain signal to be detected into a two-dimensional matrix data of M*N, wherein M is the number of magnetic brain sensors, and N is the set time length; 2)依次以各脑磁传感器为选定传感器,计算选定传感器与设定邻域内各传感器的相关系数,并根据选定传感器二维矩阵数据、设定邻域内各脑磁传感器的二维矩阵数据及相关系数,计算各选定传感器的加强信号;2) Take each brain magnetism sensor as the selected sensor in turn, calculate the correlation coefficient between the selected sensor and each sensor in the set neighborhood, and set the two-dimensional matrix of each brain magnetism sensor in the neighborhood according to the two-dimensional matrix data of the selected sensor. data and correlation coefficients to calculate the enhanced signal of each selected sensor; 3)根据加强信号及各脑磁传感器的位置,得到脑磁信号视图;3) According to the enhanced signal and the position of each magnetoencephalography sensor, the view of the magnetoencephalography signal is obtained; 4)提取脑磁信号视图的信号空间分布特征,并对所述信号空间分布特征进行分类,判断待检测信号中包含的单次听觉诱发信号或噪声信号。4) Extracting the signal spatial distribution feature of the magnetic brain signal view, classifying the signal spatial distribution feature, and judging the single auditory evoked signal or noise signal contained in the signal to be detected. 2.如权利要求1所述的方法,其特征在于,通过以下策略计算各选定传感器的增加信号:2. The method of claim 1, wherein the increased signal of each selected sensor is calculated by the following strategy: 1)选取相关系数大于阈值c的相应设定邻域内传感器Xj信号Gj,并将该相应设定邻域内传感器Xj的信号Gj以相关系数Tj为权重,叠加到选定传感器Xi的信号Gi上,得到的信号
Figure FDA0002702249310000011
为选定传感器Xi的新信号,其中,1≤i≤I,I为选定传感器的数量,1≤j≤J,J为相关系数大于阈值c的设定邻域内传感器的数量;
1) Select the signal G j of the sensor X j in the corresponding set neighborhood with the correlation coefficient greater than the threshold c, and superimpose the signal G j of the sensor X j in the corresponding set neighborhood with the correlation coefficient T j as the weight, and superimpose it on the selected sensor X i on the signal Gi , the resulting signal
Figure FDA0002702249310000011
is the new signal of the selected sensor X i , where 1≤i≤I, I is the number of selected sensors, 1≤j≤J, J is the number of sensors in the set neighborhood whose correlation coefficient is greater than the threshold c;
2)若相关系数都不大于阈值c,则选定传感器的增加信号Gi'=Gi2) If the correlation coefficients are not greater than the threshold value c, the increase signal G i '=G i of the selected sensor.
3.如权利要求1所述的方法,其特征在于,所述脑磁信号视图为2D视图;通过以下步骤得到所述脑磁信号视图:3. The method of claim 1, wherein the magnetic brain signal view is a 2D view; and the magnetic brain signal view is obtained by the following steps: 1)计算每个脑磁传感器加强信号的平均功率;1) Calculate the average power of the enhanced signal of each brain sensor; 2)将该平均功率用颜色编码,绘制到脑磁传感器对应的空间位置,得到脑磁信号视图。2) The average power is color-coded and drawn to the corresponding spatial position of the magnetic brain sensor to obtain a view of the magnetic brain signal. 4.如权利要求1所述的方法,其特征在于,通过一个训练好的二维深度卷积神经网络,提取脑磁信号视图的信号空间分布特征,并对所述信号空间分布特征进行分类;其中通过以下步骤训练二维深度卷积神经网络:4. The method according to claim 1, wherein, by a trained two-dimensional deep convolutional neural network, the signal space distribution feature of the encephalomagnetic signal view is extracted, and the signal space distribution feature is classified; where a two-dimensional deep convolutional neural network is trained by the following steps: a)将若干的单次听觉诱发样本信号与样本噪声信号,分别转换为带有标签的样本数据;a) Convert several single auditory evoked sample signals and sample noise signals into labeled sample data respectively; b)依次以各脑磁传感器为选定传感器,计算选定传感器与设定样本邻域内各脑磁传感器的样本相关系数,并根据选定传感器样本数据、设定样本邻域内各脑磁传感器的样本数据及样本相关系数,计算各选定传感器的样本增加信号;b) Take each brain magnetism sensor as the selected sensor in turn, calculate the sample correlation coefficient between the selected sensor and each brain magnetism sensor in the set sample neighborhood, and set the sample data of each brain magnetism sensor in the sample neighborhood according to the selected sensor sample data. Sample data and sample correlation coefficient, calculate the sample increase signal of each selected sensor; c)根据各样本增加信号及各脑磁传感器位置,得到若干的单次听觉诱发样本脑磁信号视图或样本噪声视图;c) according to the increase of the signal of each sample and the position of each EEG sensor, to obtain several single auditory evoked sample EEG signal views or sample noise views; d)通过听觉诱发样本脑磁信号视图与样本噪声视图,对二维深度卷积神经网络进行训练,得到训练好的二维深度卷积神经网络。d) Train the two-dimensional deep convolutional neural network through the auditory evoked sample encephalomagnetism signal view and the sample noise view, and obtain a trained two-dimensional deep convolutional neural network. 5.如权利要求4所述的方法,其特征在于,所述二维深度卷积神经网络包括GoogLeNet网络、AlexNet网络或节点数与权重符合要求的二维深度卷积神经网络。5. The method according to claim 4, wherein the two-dimensional deep convolutional neural network comprises a GoogLeNet network, an AlexNet network or a two-dimensional deep convolutional neural network with the required number of nodes and weights. 6.如权利要求4所述的方法,其特征在于,将单次听觉诱发样本信号转换为带有标签的样本数据之前,对单次听觉诱发样本信号进行预处理;所述预处理包括:去除单次听觉诱发样本信号中的环境噪声信号、去除眼动与心动伪迹噪声信号和去除脑部发出的其他信号。6. The method according to claim 4, wherein before converting the single auditory evoked sample signal into labeled sample data, preprocessing is performed on the single auditory evoked sample signal; the preprocessing comprises: removing the Ambient noise from single auditory evoked samples, removal of eye movement and cardiac artifact noise, and removal of other signals from the brain. 7.如权利要求4所述的方法,其特征在于,环境噪声信号包括:工频噪声;去除工频噪声的方法包括:使用50Hz或60Hz的滤波去除;去除眼动与心动伪迹噪声信号的方法包括:使用信号空间投影算法或独立成分分析算法。7. The method according to claim 4, wherein the environmental noise signal comprises: power frequency noise; the method for removing power frequency noise comprises: using 50Hz or 60Hz filtering to remove; removing eye movement and cardiac artifact noise signals Methods include: using a signal space projection algorithm or an independent component analysis algorithm. 8.如权利要求4所述的方法,其特征在于,脑部发出的其他信号包括:α信号和β信号;去除α信号和β信号的方法包括:使用10Hz和20Hz的滤波去除。8 . The method of claim 4 , wherein the other signals emitted by the brain include: alpha signals and beta signals; the method for removing the alpha signals and the beta signals comprises: removing by filtering at 10 Hz and 20 Hz. 9 . 9.一种存储介质,所述存储介质中存储有计算机程序,其中,所述计算机程序被设置为运行时执行权利要求1-8中任一所述方法。9. A storage medium in which a computer program is stored, wherein the computer program is configured to execute the method of any one of claims 1-8 when run. 10.一种电子装置,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器被设置为运行所述计算机程序以执行如权利要求1-8中任一所述方法。10. An electronic device comprising a memory and a processor having a computer program stored in the memory, the processor being arranged to run the computer program to perform the method of any of claims 1-8.
CN202011026442.5A 2020-09-25 2020-09-25 Method and electronic device for detecting single auditory evoked signal in magnetoencephalography Active CN112168167B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011026442.5A CN112168167B (en) 2020-09-25 2020-09-25 Method and electronic device for detecting single auditory evoked signal in magnetoencephalography

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011026442.5A CN112168167B (en) 2020-09-25 2020-09-25 Method and electronic device for detecting single auditory evoked signal in magnetoencephalography

Publications (2)

Publication Number Publication Date
CN112168167A true CN112168167A (en) 2021-01-05
CN112168167B CN112168167B (en) 2022-05-03

Family

ID=73944983

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011026442.5A Active CN112168167B (en) 2020-09-25 2020-09-25 Method and electronic device for detecting single auditory evoked signal in magnetoencephalography

Country Status (1)

Country Link
CN (1) CN112168167B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113988143A (en) * 2021-12-27 2022-01-28 北京航空航天大学杭州创新研究院 Bad track automatic identification method based on optical pump magnetometer measurement brain magnetism
CN117122288A (en) * 2023-09-08 2023-11-28 太原理工大学 Epileptic electroencephalogram signal early warning method and device based on anchoring convolution network

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5224486B1 (en) * 2012-03-08 2013-07-03 学校法人金沢工業大学 Double gradiometer and brain magnetic field measuring device
US20140107464A1 (en) * 2012-06-29 2014-04-17 Commissariat A L'energie Atomique Et Aux Energies Alternatives Method for Locating a Brain Activity
US20140378815A1 (en) * 2011-05-24 2014-12-25 The Regents Of The University Of California Magnetoencephalography source imaging
US20160051162A1 (en) * 2014-07-29 2016-02-25 Commissariat A L'energie Atomique Et Aux Energies Alternatives Method for locating a brain activity associated with a task
US20160157742A1 (en) * 2013-08-05 2016-06-09 The Regents Of The University Of California Magnetoencephalography source imaging for neurological functionality characterizations
CN109009065A (en) * 2018-06-05 2018-12-18 上海理工大学 Brain magnetic information detection system and method based on TMR weak magnetic sensor array
US20190099101A1 (en) * 2017-10-03 2019-04-04 Brain F.I.T. Imaging, LLC Methods and magnetic imaging devices to inventory human brain cortical function
US20190201657A1 (en) * 2016-08-12 2019-07-04 The Board Of Trustees Of The Leland Stanford Junior University Device and method for hearing threshold-adapted acoustic stimulation
US20200218338A1 (en) * 2017-04-07 2020-07-09 Yoonhee Lee System for communicating sensory information with an interactive system and methods thereof
CN111543949A (en) * 2020-05-13 2020-08-18 北京航空航天大学 A diagnostic device for children's ASD based on magnetoencephalography and electroencephalography

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140378815A1 (en) * 2011-05-24 2014-12-25 The Regents Of The University Of California Magnetoencephalography source imaging
JP5224486B1 (en) * 2012-03-08 2013-07-03 学校法人金沢工業大学 Double gradiometer and brain magnetic field measuring device
US20140107464A1 (en) * 2012-06-29 2014-04-17 Commissariat A L'energie Atomique Et Aux Energies Alternatives Method for Locating a Brain Activity
US20160157742A1 (en) * 2013-08-05 2016-06-09 The Regents Of The University Of California Magnetoencephalography source imaging for neurological functionality characterizations
US20160051162A1 (en) * 2014-07-29 2016-02-25 Commissariat A L'energie Atomique Et Aux Energies Alternatives Method for locating a brain activity associated with a task
US20190201657A1 (en) * 2016-08-12 2019-07-04 The Board Of Trustees Of The Leland Stanford Junior University Device and method for hearing threshold-adapted acoustic stimulation
US20200218338A1 (en) * 2017-04-07 2020-07-09 Yoonhee Lee System for communicating sensory information with an interactive system and methods thereof
US20190099101A1 (en) * 2017-10-03 2019-04-04 Brain F.I.T. Imaging, LLC Methods and magnetic imaging devices to inventory human brain cortical function
CN109009065A (en) * 2018-06-05 2018-12-18 上海理工大学 Brain magnetic information detection system and method based on TMR weak magnetic sensor array
CN111543949A (en) * 2020-05-13 2020-08-18 北京航空航天大学 A diagnostic device for children's ASD based on magnetoencephalography and electroencephalography

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
RUI ZHANG ET AL: "Recording brain activities in unshielded Earth’s field with optically pumped atomic magnetometers", 《SCIENCE ADVANCES》 *
SCHUBRING, DAVID,: "Consistency of abnormal sensory gating in first‐admission and chronic schizophrenia across quantification methods.", 《PSYCHOPHYSIOLOGY》 *
蒋玉尔 等: "同类图片命名中信息处理的脑磁图研究", 《中国康复医学杂志》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113988143A (en) * 2021-12-27 2022-01-28 北京航空航天大学杭州创新研究院 Bad track automatic identification method based on optical pump magnetometer measurement brain magnetism
CN117122288A (en) * 2023-09-08 2023-11-28 太原理工大学 Epileptic electroencephalogram signal early warning method and device based on anchoring convolution network

Also Published As

Publication number Publication date
CN112168167B (en) 2022-05-03

Similar Documents

Publication Publication Date Title
CN112656427B (en) Electroencephalogram emotion recognition method based on dimension model
Seeliger et al. Convolutional neural network-based encoding and decoding of visual object recognition in space and time
US4736751A (en) Brain wave source network location scanning method and system
CN113420624A (en) Non-contact fatigue detection method and system
CN113261932B (en) Heart rate measurement method and device based on PPG signal and one-dimensional convolutional neural network
Ebrahimzadeh et al. A novel approach for detection of deception using Smoothed Pseudo Wigner-Ville Distribution (SPWVD)
CN112022153B (en) EEG signal detection method based on convolutional neural network
JP2019516434A (en) Method and magnetic imaging device for cataloging cortical function of human brain
CN114041786A (en) Ballistocardiogram signal detection method, ballistocardiogram signal detection device and ballistocardiogram signal detection equipment
CN111568412A (en) Method and device for reconstructing visual image by utilizing electroencephalogram signal
Qian et al. Decision-level fusion of EEG and pupil features for single-trial visual detection analysis
Wu et al. A real-time method to reduce ballistocardiogram artifacts from EEG during fMRI based on optimal basis sets (OBS)
Raghu et al. Automated biomedical signal quality assessment of electromyograms: Current challenges and future prospects
CN112168167B (en) Method and electronic device for detecting single auditory evoked signal in magnetoencephalography
CN111278353A (en) Method and system for detecting vital sign signal noise
CN112220482A (en) Magnetoencephalography eye movement artifact detection and removal method and electronic device based on neural network
Bathellier et al. Wavelet-based multi-resolution statistics for optical imaging signals: application to automated detection of odour activated glomeruli in the mouse olfactory bulb
CN114569116B (en) Ballistocardiogram ventricular fibrillation auxiliary diagnosis system based on three-channel images and transfer learning
CN113951886A (en) Brain magnetic pattern generation system and lie detection decision system
CN120036730A (en) Neonatal pain classification method and system based on multi-modal evidence
Lee et al. Generating target/non-target images of an RSVP experiment from brain signals in by conditional generative adversarial network
Rathi et al. Detection and characterization of brain tumor using segmentation based on HSOM, wavelet packet feature spaces and ANN
Yang et al. Heart rate estimation from facial videos based on convolutional neural network
KR101413853B1 (en) Method and apparatus for measuring physiological signal usuing infrared image
Zilber ERF and scale-free analyses of source-reconstructed MEG brain signals during a multisensory learning paradigm

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant