CN112674773B - Magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window - Google Patents
Magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window Download PDFInfo
- Publication number
- CN112674773B CN112674773B CN202011526796.6A CN202011526796A CN112674773B CN 112674773 B CN112674773 B CN 112674773B CN 202011526796 A CN202011526796 A CN 202011526796A CN 112674773 B CN112674773 B CN 112674773B
- Authority
- CN
- China
- Prior art keywords
- magnetoencephalogram
- source
- time window
- signal
- ripple
- 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.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/242—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
- A61B5/245—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetoencephalographic [MEG] signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
- A61B5/4094—Diagnosing or monitoring seizure diseases, e.g. epilepsy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Power Engineering (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
A magnetoencephalography source positioning method and device based on Tucker decomposition and ripple time window are provided, the device comprises: the magnetoencephalogram sensor is used for acquiring a first magnetoencephalogram signal of a user; the ripple detection unit is used for detecting a ripple time window in the first magnetoencephalogram signal as a source-positioned time window through a root mean square method, and obtaining the first magnetoencephalogram signal in the ripple time window as a second magnetoencephalogram signal; a Tucker decomposition unit based on high-order orthogonal iteration, configured to perform Tucker decomposition on an original tensor of the second magnetoencephalogram signal to calculate an estimated value of the original tensor; and the source positioning unit is used for calculating a covariance matrix for the estimated value and calculating the source position corresponding to the second magnetoencephalogram signal by an LCMV inverse problem solving method in a beam forming method. The method eliminates the influence of noise signals, reduces the complexity of calculation, ensures the consistency of each calculation result and improves the accuracy of positioning the seizure area.
Description
Technical Field
The invention relates to a magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window, belonging to the field of biomedicine.
Background
Clinically, up to one third of epileptic patients are refractory epilepsy. While drug therapy with antiepileptic drugs is difficult to control effectively. Epileptic foci can affect cortical excitability, resulting in abnormal discharges, not only directly affecting focal or surrounding tissues, but also affecting distant brain regions. Therefore, these patients need to surgically excise epileptogenic foci (epileptogenic zones) to achieve seizure-free, and whether or not an epileptogenic focus (epileptogenic focus) is postoperatively dependent on the localization of the epileptogenic focus. However, it is difficult to determine the area due to the lack of a tool to directly measure the area.
The Magnetoencephalogram (MEG) is a method for monitoring brain functions in a non-invasive and real-time manner, and in the past decades, several brain source localization methods based on electroencephalogram or electroencephalogram, such as spike-based dipole-fitting methods, have been proposed, but the localization accuracy of these methods still cannot meet the requirement of preoperative accurate localization, the calculation complexity is high, the localization result is susceptible to noise signals, and for the same magnetoencephalogram signal, the results of multiple localizations have deviation and poor consistency.
In view of the above, the present invention is directed to a method and an apparatus for positioning a magnetoencephalogram source based on a Tucker decomposition and a ripple time window to solve one or more of the above technical problems.
Disclosure of Invention
In order to solve one or more technical problems in the prior art, the applicant finds, through long-term research, that accuracy of location of an epileptogenic focus (epileptogenic region) can be improved by extracting a ripple time window of an original magnetoencephalogram signal and combining with Tucker decomposition, and reference and basis can be provided for subsequent operations. Based on the research result, according to an aspect of the present invention, there is provided a magnetoencephalography source positioning device based on a Tucker decomposition and a ripple time window, comprising:
the magnetoencephalogram sensor is used for acquiring a first magnetoencephalogram signal of a user;
a ripple detection unit, configured to detect a ripple time window in the first magnetoencephalogram signal as a source-located time window by a root mean square method, and obtain the first magnetoencephalogram signal in the ripple time window as a second magnetoencephalogram signal, where at least four continuous oscillation signals having a frequency of 80-250Hz and an amplitude higher than a background signal are defined as ripple;
a Tucker decomposition unit based on high-order orthogonal iteration, configured to perform Tucker decomposition on an original tensor of the second magnetoencephalogram signal to calculate an estimated value of the original tensor, where the estimated value calculated by the Tucker decomposition is constrained by the high-order orthogonal iteration, so that the estimated value remains unique, and a noise factor tensor in the original tensor is removed from the estimated value; and
and the source positioning unit is used for calculating a covariance matrix for the estimated value and calculating a source position corresponding to the second magnetoencephalogram signal by an LCMV inverse problem solving method in beam forming.
According to another aspect of the invention, the calculation process of the higher order orthogonal iteration is as follows:
1) calculating factor matrix U of original tensor X by using high-order singular value decomposition algorithm(n)And a core tensor G; let k equal to 0;
2) let k be k +1 and N be 1,2, … N, the following operations are performed:
B(k)←X×1U(1)T…×n-1U(n-1)T×NUNT；
and performs the core tensor B(k)Determining the number Rn of main singular values of the singular value decomposition of the n-mode expansion, wherein B(k)For the core tensor obtained for the kth iteration, then the operations are performed:
U(N)←U(:,1:Rn)；
3) calculating a core tensor B obtained by the k iteration(k)Meter G(k)＝B(k)And judging whether convergence occurs by the following formula:
||G(k)-G(k-1)||F<e; wherein E is a minimum value;
if the convergence condition is met, executing the next step, otherwise, returning to the step 2) to continue iteration until convergence;
4) output factor matrix U(N)And the core tensor G(k)(ii) a According to the factor matrix U(N)And the core tensor G(k)Reconstructing the estimate of the original tensor.
According to another aspect of the present invention, the magnetoencephalogram source positioning device based on the Tucker decomposition and ripple time window further comprises a source display unit for displaying the source position.
According to yet another aspect of the invention, the source location is a seizure region of the brain.
According to another aspect of the present invention, the apparatus for source localization of magnetoencephalogram based on Tucker decomposition and ripple time window further includes a preprocessing unit, configured to perform filtering, linear trend elimination and independent component analysis on the first magnetoencephalogram signal to obtain a preprocessed first magnetoencephalogram signal, and the preprocessed first magnetogram signal is input to the ripple detecting unit for subsequent processing.
According to another aspect of the present invention, the ripple time window is a plurality of windows, and for the second magnetoencephalogram signal corresponding to each ripple time window, the source localization unit calculates a plurality of corresponding source locations, and concentrates a corresponding region according to the plurality of source locations as a localization result.
According to another aspect of the present invention, there is also provided a method for positioning a magnetoencephalogram source based on a Tucker decomposition and a ripple time window, comprising the steps of:
acquiring a first magnetoencephalogram signal of a user through a magnetoencephalogram sensor;
detecting a ripple time window in the first magnetoencephalogram signal as a source-positioned time window through a root mean square method, and obtaining the first magnetoencephalogram signal in the ripple time window as a second magnetoencephalogram signal, wherein at least four continuous oscillation signals with the frequency of 80-250Hz and the amplitude higher than that of a background signal are defined as ripple;
performing Tucker decomposition on an original tensor of the second magnetoencephalogram signal to calculate an estimated value of the original tensor, wherein the estimated value calculated by the Tucker decomposition is constrained by high-order orthogonal iteration to keep the estimated value unique, and a noise factor tensor in the original tensor is removed from the estimated value;
and calculating a covariance matrix for the estimated value, and calculating a source position corresponding to the second magnetoencephalogram signal by an LCMV inverse problem solving method.
According to still another aspect of the present invention, the source location is displayed by a source display unit.
According to yet another aspect of the invention, the source location is taken as a seizure region of the brain.
According to another aspect of the invention, the preprocessing unit filters, trends linearly and analyzes independent components of the first magnetoencephalogram signal to obtain a preprocessed first magnetoencephalogram signal, and the preprocessed first magnetogram signal is used for subsequent ripple time window detection.
According to another aspect of the present invention, the ripple time window is a plurality of windows, for each second magnetoencephalogram signal corresponding to the ripple time window, a plurality of corresponding source positions are calculated, and a region having the first concentration corresponding to the plurality of source positions is used as a positioning result.
Compared with the prior art, the invention has one or more of the following technical effects:
(1) in the magnetoencephalogram source positioning, high-frequency noise and low-frequency noise in magnetoencephalogram signals are effectively removed through a Tucker decomposition method, the influence of noise signals on the positioning result of epileptogenic foci (epileptogenic regions) is reduced, and meanwhile, the complexity of calculation is reduced.
(2) The result of the Tucker decomposition is subjected to uniqueness constraint by using a high-order orthogonal iteration (HOOI) method, and the consistency of the calculation result every time is ensured.
(3) The accuracy of the location of the epileptic focus (epileptogenic region) is further improved by taking ripple as a time window for source location selection and combining with Tucker decomposition, and reference and basis are provided for subsequent operations.
(4) By detecting a plurality of ripple time windows and obtaining a plurality of corresponding source positions, the final positioning result is determined according to the corresponding areas in the source position sets, and the positioning precision is further improved; combining multiple source localization results, taking the brain region with the first concentration (45% -85%) of more than 10 brain regions (source localization regions) as the localization result can further reduce the workload of doctors and eliminate potential interference.
Drawings
So that the manner in which the above recited features of the present invention can be understood in detail, a more particular description of the invention, briefly summarized above, may be had by reference to embodiments. The drawings relate to preferred embodiments of the invention and are described below:
FIG. 1 is a schematic diagram of a magnetoencephalography source positioning device based on a Tucker decomposition and ripple time window according to a preferred embodiment of the present invention;
FIG. 2 is a diagram illustrating the detection result of a ripple time window according to a preferred embodiment of the present invention;
FIG. 3 is a schematic diagram showing the comparison between the positioning accuracy of the tracker decomposition based on the ripple time window and the positioning accuracy of the dipole-fitting based on the spike.
DETAILED DESCRIPTION OF EMBODIMENT (S) OF INVENTION
Reference will now be made in detail to the various embodiments, one or more examples of which are illustrated in each figure. The examples are provided by way of explanation and are not meant as limitations. For example, features illustrated or described as part of one embodiment can be used on or in conjunction with any other embodiment to yield yet a further embodiment. It is intended that the present invention encompass such modifications and variations.
In the following description of the drawings, like reference numerals designate identical or similar structures. Generally, only the differences between the individual embodiments will be described. Descriptions of parts or aspects in one embodiment can also be applied to corresponding parts or aspects in another embodiment, unless explicitly stated otherwise.
Example 1
Referring to fig. 1-3, fig. 1 is a schematic diagram of a positioning apparatus for a magnetoencephalogram source based on a Tucker decomposition and ripple time window according to a preferred embodiment of the present invention; FIG. 2 is a diagram illustrating the detection result of a ripple time window according to a preferred embodiment of the present invention; FIG. 3 is a schematic diagram showing the comparison between the positioning accuracy of the tracker decomposition based on the ripple time window and the positioning accuracy of the dipole-fitting based on the spike. In accordance with a preferred embodiment of the present invention, referring to fig. 1-2, there is provided a magnetoencephalography source localization apparatus based on Tucker decomposition and ripple time window, comprising:
the magnetoencephalogram sensor is used for acquiring a first magnetoencephalogram signal of a user;
a ripple detection unit, configured to detect a ripple time window in the first magnetoencephalogram signal as a source-located time window by a root mean square method, and obtain the first magnetoencephalogram signal in the ripple time window as a second magnetoencephalogram signal, where at least four continuous oscillation signals having a frequency of 80-250Hz and an amplitude higher than a background signal are defined as ripple;
a Tucker decomposition unit based on high-order orthogonal iteration, configured to perform Tucker decomposition on an original tensor of the second magnetoencephalogram signal to calculate an estimated value of the original tensor, where the estimated value calculated by the Tucker decomposition is constrained by the high-order orthogonal iteration, so that the estimated value remains unique, and a noise factor tensor in the original tensor is removed from the estimated value; and
and the source positioning unit is used for calculating a covariance matrix for the estimated value and calculating a source position corresponding to the second magnetoencephalogram signal by an LCMV (linear constrained minimum variance) inverse problem solving method in a beam forming method.
Advantageously, the method improves the detection precision of the source position by the Tucker decomposition based on the high-order orthogonal iteration (HOOI) and taking the ripple time period of the MEG signal as the positioning time window, can be widely applied to deep source and shallow source positioning, and improves the reference value of clinical brain source positioning in preoperative evaluation.
According to another preferred embodiment of the present invention, the calculation process of the higher order orthogonal iteration is as follows:
1) calculating factor matrix U of original tensor X by using high-order singular value decomposition algorithm(n)And a core tensor G; let k equal to 0;
2) let k be k +1 and N be 1,2, … N, the following operations are performed:
B(k)←X×1U(1)T…×n-1U(n-1)T×NUNT；
and performs the core tensor B(k)Determining the number Rn of main singular values of the singular value decomposition of the n-mode expansion, wherein B(k)For the core tensor obtained for the kth iteration, then the operations are performed:
U(N)←U(:,1:Rn)；
3) calculating a core tensor B obtained by the k iteration(k)Meter G(k)＝B(k)And judging whether convergence occurs by the following formula:
||G(k)-G(k-1)||F<e; wherein E is a minimum value;
if the convergence condition is met, executing the next step, otherwise, returning to the step 2) to continue iteration until convergence;
4) output factor matrix U(N)And the core tensor G(k)(ii) a According to the factor matrix U(N)And a coreCardiac tensor G(k)Reconstructing the estimate of the original tensor.
According to another preferred embodiment of the present invention, the apparatus for source localization of magnetoencephalogram based on Tucker decomposition and ripple time window further comprises a source display unit for displaying the source location.
According to yet another preferred embodiment of the invention, the source location is an epileptogenic zone (epileptogenic zone) of the brain.
According to another preferred embodiment of the present invention, the apparatus for source localization of magnetoencephalogram based on Tucker decomposition and ripple time window further comprises a preprocessing unit, configured to perform filtering, linear trend elimination and Independent Component Analysis (ICA) on the first magnetoencephalogram signal to obtain a preprocessed first magnetoencephalogram signal, and the preprocessed first magnetogram signal is input to the ripple detecting unit for subsequent processing.
According to another preferred embodiment of the present invention, the ripple time window is a plurality of windows, and the source localization unit calculates a plurality of corresponding source locations for the second magnetoencephalogram signal corresponding to each ripple time window, and concentrates a corresponding region as a localization result according to the plurality of source locations. Specifically, for example, a region corresponding to at least Y iterations of the Z source positions is output as a final result, and Y is preferably equal to or greater than 3. According to experimental verification, the repeatedly corresponding (indicating) areas for more than 3 times are basically close to the real epileptogenic areas, so that the positioning accuracy is further improved, and reliable reference is provided for later operations of doctors.
Preferably, the ripple time windows are Z, and for the second magnetoencephalogram signal corresponding to each ripple time window, the source localization unit calculates Z source locations corresponding to each ripple time window, determines Z brain regions corresponding to the Z source locations respectively, and takes the brain region with the first concentration in the Z brain regions as a localization result. Research shows that the Z is more than or equal to 10, the first concentration ratio is 45% -85%, accurate positioning of a real seizure area can be achieved, and compared with single source positioning, the method can further reduce workload of doctors and eliminate potential interference.
According to another preferred embodiment of the present invention, there is also provided a method for positioning a magnetoencephalogram source based on a Tucker decomposition and a ripple time window, comprising the steps of:
acquiring a first magnetoencephalogram signal of a user through a magnetoencephalogram sensor;
detecting a ripple time window in the first magnetoencephalogram signal as a source-positioned time window through a root mean square method, and obtaining the first magnetoencephalogram signal in the ripple time window as a second magnetoencephalogram signal, wherein at least four continuous oscillation signals with the frequency of 80-250Hz and the amplitude higher than that of a background signal are defined as ripple;
performing Tucker decomposition on an original tensor of the second magnetoencephalogram signal to calculate an estimated value of the original tensor, wherein the estimated value calculated by the Tucker decomposition is constrained by high-order orthogonal iteration to keep the estimated value unique, and a noise factor tensor in the original tensor is removed from the estimated value;
and calculating a covariance matrix according to the estimated value, and calculating a source position corresponding to the second magnetoencephalogram signal by an LCMV inverse problem solving method in beam forming.
According to still another preferred embodiment of the present invention, the source position is displayed by a source display unit.
According to a further preferred embodiment of the invention, the source location is a seizure region of the brain.
According to another preferred embodiment of the present invention, the pre-processing unit filters, trends linearly, and analyzes the independent components of the first magnetoencephalogram signal to obtain a pre-processed first magnetoencephalogram signal, which is used for subsequent ripple time window detection.
According to another preferred embodiment of the present invention, the ripple time window is a plurality of windows, and for the second magnetoencephalogram signal corresponding to each ripple time window, a plurality of corresponding source positions are calculated, and a corresponding region in the plurality of source positions is determined as a positioning result.
According to another preferred embodiment of the present invention, there is also provided a novel source localization apparatus based on Tucker decomposition and ripple time window, including a ripple detection unit, a tensor decomposition unit, a source localization unit, and a verification unit. Preferably, the Ripple detection unit applies a root mean square method. The MEG (magnetoencephalography) signal adopts Tucker decomposition to calculate the estimation value of an original tensor, the Tucker decomposition result is kept unique based on the mode constraint of high-order orthogonal iteration, and the factor tensor of which the main information is noise is removed. And the source positioning unit calculates a covariance matrix for the estimated value of the tensor decomposed by the Tucker, so that the source position of the epileptic brain magnetic signal is calculated by an LCMV method. The validation unit includes validation on top of the real magnetoencephalography data, wherein the real data-based validation compares the ripple-based Tucker decomposition with spike-based clinical dipole methods. As shown in fig. 3, which shows the comparison result, wherein the present invention greatly improves the positioning accuracy of the seizure-causing region of the brain with respect to the spike-based clinical dipole method.
Preferably, the ripple detection unit defines at least four continuous oscillations of 80-250Hz and amplitude above the background signal as ripple. Preferably, at least four oscillations of 80-250Hz and amplitudes above the first threshold of the background signal are defined as ripple, further improving the positioning accuracy. Preferably, in the present invention, a Root Mean Square (RMS) method is used for detecting ripples in MEG, and the conditions thereof can be specifically selected as follows: the waveform power value should be 3-9 times the background signal and the duration of the waveform should last at least 15ms (corresponding to at least four oscillations). Preferably, false positive ripples can be eliminated manually.
According to another preferred embodiment of the present invention, in the orthogonal iteration process, the original tensor X is input; outputting the factor matrix U after iteration is finished(n)And a core tensor G.
The calculation process of the orthogonal iteration is preferably as follows:
1) calculating a factor matrix U of an original tensor X by using a high-order singular value decomposition (HOSVD) algorithm(n)And a core tensor G. Let k equal to 0;
2) let k be k +1 and N be 1,2, … N, the following operations are performed:
B(k)←X×1U(1)T…×n-1U(n-1)T×NUNT(ii) a Wherein, X is multiplied by nU(n)TRepresenting tensor X and matrix U(n) of Tn-modulo product;
and performs the core tensor B(k)Determining the number Rn of main singular values of the SVD, and then executing the operation:
U(N)←U(:,1:Rn)
3) computing the core tensor B for the kth iteration(k)Meter G(k)＝B(k)And whether to converge is judged by the following convergence formula:
||G(k)-G(k-1)||F<E
e is a minimum value.
If the convergence condition is met, executing the next step, otherwise returning to the step 2 to continue iteration until convergence.
4) Output factor matrix U(N)And the core tensor G(k)(ii) a According to the factor matrix U(N)And the core tensor G(k)Reconstructing the estimate of the original tensor.
It is understood that, among other things, the core tensor G is initialized(0)Is the zero tensor (all elements are 0). I G(k)-G(k-1)||F<E denotes determination of G(k)-G(k-1)Whether a convergence formula of a convergence condition is met.
Preferably, for the estimated value, an inverse problem is solved by using a beamforming (beamforming) analysis method to calculate the covariance matrix and the position of the estimated source, and the principle of beamforming is as follows: the signal is filtered through time domain information and frequency domain information to obtain information in a given direction while attenuating noise interference from other directions.
Compared with the prior art, the invention has one or more of the following technical effects:
(1) in the magnetoencephalogram source positioning, high-frequency noise and low-frequency noise in magnetoencephalogram signals are effectively removed through a Tucker decomposition method, the influence of noise signals on the positioning result of epileptogenic foci (epileptogenic regions) is reduced, and meanwhile, the complexity of calculation is reduced.
(2) The result of the Tucker decomposition is subjected to uniqueness constraint by using a high-order orthogonal iteration (HOOI) method, and the consistency of the calculation result every time is ensured.
(3) The accuracy of the location of the epileptic focus (epileptogenic region) is further improved by taking ripple as a time window for source location selection and combining with Tucker decomposition, and reference and basis are provided for subsequent operations.
(4) By detecting a plurality of ripple time windows and obtaining a plurality of corresponding source positions, the final positioning result is determined according to the corresponding areas in the source position sets, and the positioning precision is further improved; combining multiple source localization results, taking the brain region with the first concentration (45% -85%) of more than 10 brain regions (source localization regions) as the localization result can further reduce the workload of doctors and eliminate potential interference.
While the foregoing is directed to embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.
The above embodiments are merely preferred embodiments of the present invention, which are not intended to limit the present invention, and the features of the embodiments that do not violate each other may be combined with each other. Any modification, equivalent replacement, and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (5)
1. A magnetoencephalography source positioning device based on Tucker decomposition and ripple time window is characterized by comprising:
the magnetoencephalogram sensor is used for acquiring a first magnetoencephalogram signal of a user;
a ripple detection unit, configured to detect a ripple time window in the first magnetoencephalogram signal as a source-located time window by a root mean square method, and obtain the first magnetoencephalogram signal in the ripple time window as a second magnetoencephalogram signal, where at least four continuous oscillation signals having a frequency of 80-250Hz and an amplitude higher than a background signal are defined as ripple;
a Tucker decomposition unit based on high-order orthogonal iteration, configured to perform Tucker decomposition on an original tensor of the second magnetoencephalogram signal to calculate an estimated value of the original tensor, where the estimated value calculated by the Tucker decomposition is constrained by the high-order orthogonal iteration, so that the estimated value remains unique, and a noise factor tensor in the original tensor is removed from the estimated value; and
a source positioning unit, configured to calculate a covariance matrix for the estimated value, and calculate a source position corresponding to the second magnetoencephalogram signal by an LCMV inverse problem solution method in beamforming;
wherein, the calculation process of the high-order orthogonal iteration is as follows:
1) calculating factor matrix of original tensor X by using high-order singular value decomposition algorithmU （n）And a core tensor B（k）(ii) a Let k = 0;
2) let k = k +1 and for N =1,2, … N, perform the following operations:
B（k）←X×1 U （1）T …×n-1 U （n-1）T× N U NT；
and performs the core tensor B（k）Determining the number Rn of main singular values of the singular value decomposition of the n-mode expansion, wherein B（k）For the core tensor obtained for the kth iteration, then the operations are performed:
U (N) ←U(:,1:Rn)；
3) calculating a core tensor B obtained by the k iteration（k）Meter G(k)= B（k）And judging whether convergence occurs by the following formula:
‖G(k) -G(k-1)‖F<e; wherein E is a minimum value;
if the convergence condition is met, executing the next step, otherwise, returning to the step 2) to continue iteration until convergence;
4) output G(k)Sum factor matrix U (N) (ii) a According to G(k)Sum factor matrix U (N) Reconstructing the estimate of the original tensor.
2. The tile decomposition and ripple time window based magnetoencephalogram source localization apparatus of claim 1, further comprising a source display unit for displaying the source location.
3. The magnetocaloric source localization device according to claim 2, wherein the source location is the seizure area of the brain.
4. The tile decomposition and ripple time window based magnetoencephalogram source localization device according to any one of claims 1-3, further comprising a preprocessing unit for performing filtering, linear trend elimination and independent component analysis on the first magnetoencephalogram signal to obtain a preprocessed first magnetoencephalogram signal, and the preprocessed first magnetogram signal is input to the ripple detection unit for subsequent processing.
5. The apparatus according to any of claims 1-3, wherein the ripple time window is a plurality of ripple time windows, and for each second magnetoencephalogram signal corresponding to the ripple time window, the source localization unit calculates a plurality of source locations corresponding to each second magnetoencephalogram signal, and uses a region having a first concentration corresponding to the plurality of source locations as a localization result.
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011526796.6A CN112674773B (en) | 2020-12-22 | 2020-12-22 | Magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window |
| CN202111402038.8A CN113951885B (en) | 2020-12-22 | 2020-12-22 | Magnetoencephalogram source positioning method, storage medium and equipment |
| PCT/CN2021/140340 WO2022135448A1 (en) | 2020-12-22 | 2021-12-22 | Magnetoencephalography source localization method and device based on tucker decomposition and ripple time window |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011526796.6A CN112674773B (en) | 2020-12-22 | 2020-12-22 | Magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111402038.8A Division CN113951885B (en) | 2020-12-22 | 2020-12-22 | Magnetoencephalogram source positioning method, storage medium and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112674773A CN112674773A (en) | 2021-04-20 |
| CN112674773B true CN112674773B (en) | 2021-12-24 |
Family
ID=75450400
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011526796.6A Active CN112674773B (en) | 2020-12-22 | 2020-12-22 | Magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window |
| CN202111402038.8A Active CN113951885B (en) | 2020-12-22 | 2020-12-22 | Magnetoencephalogram source positioning method, storage medium and equipment |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111402038.8A Active CN113951885B (en) | 2020-12-22 | 2020-12-22 | Magnetoencephalogram source positioning method, storage medium and equipment |
Country Status (2)
| Country | Link |
|---|---|
| CN (2) | CN112674773B (en) |
| WO (1) | WO2022135448A1 (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112674773B (en) * | 2020-12-22 | 2021-12-24 | 北京航空航天大学 | Magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window |
| CN113951887B (en) * | 2021-12-22 | 2022-03-25 | 北京航空航天大学杭州创新研究院 | Frequency spectrum matching independent component analysis method and system |
| CN114065825B (en) * | 2022-01-17 | 2022-04-19 | 北京航空航天大学杭州创新研究院 | Brain magnetic MEG source positioning method based on structural similarity |
| CN117503057B (en) * | 2023-11-09 | 2024-05-31 | 电子科技大学 | Epileptic seizure detection device and medium for constructing brain network based on high-order tensor decomposition |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102791332A (en) * | 2009-11-04 | 2012-11-21 | 代理并代表亚利桑那州立大学的亚利桑那董事会 | Devices and methods for modulating brain activity |
| US20130144153A1 (en) * | 2011-12-01 | 2013-06-06 | The Regents Of The University Of California | Functional magnetic resonance imaging apparatus and methods |
| US20160022168A1 (en) * | 2014-07-24 | 2016-01-28 | University Of Lethbridge | Brain state dependent therapy for improved neural training and rehabilitation |
| CN107518898A (en) * | 2017-08-08 | 2017-12-29 | 北京航空航天大学 | Magneticencephalogram source positioner with beam forming is decomposed based on sensor array |
| CN110612057A (en) * | 2017-06-07 | 2019-12-24 | 柯惠有限合伙公司 | System and method for detecting stroke |
| CN111096745A (en) * | 2020-01-02 | 2020-05-05 | 苏州大学 | Steady-state evoked response brain source positioning method based on sparse Bayesian learning |
| CN111209974A (en) * | 2020-01-14 | 2020-05-29 | 华中科技大学鄂州工业技术研究院 | Tensor decomposition-based heterogeneous big data core feature extraction method and system |
| CN111340142A (en) * | 2020-05-14 | 2020-06-26 | 南京慧脑云计算有限公司 | Epilepsia magnetoencephalogram spike automatic detection method and tracing positioning system |
| US20200289005A1 (en) * | 2019-03-14 | 2020-09-17 | Ricoh Company, Ltd. | Multi-sphere head model for dipole localization |
| CN111801046A (en) * | 2017-11-10 | 2020-10-20 | 勒维斯公司 | Efficacy and/or therapeutic parameter recommendation using individual patient data and therapeutic brain network maps |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7461045B1 (en) * | 2003-08-18 | 2008-12-02 | University Of Florida Research Foundation, Inc. | Optimization of spatio-temporal pattern processing for seizure warning and prediction |
| WO2012133185A1 (en) * | 2011-03-31 | 2012-10-04 | 独立行政法人理化学研究所 | Brain wave analysis apparatus, brain wave analysis method, program and recording medium |
| EP2713867A4 (en) * | 2011-05-24 | 2015-01-21 | Univ California | Magnetoencephalography source imaging |
| CN105147288B (en) * | 2015-07-23 | 2018-06-08 | 中国科学院苏州生物医学工程技术研究所 | Brain source strength localization method |
| CN105212895B (en) * | 2015-08-24 | 2019-01-15 | 中国科学院苏州生物医学工程技术研究所 | Dynamic brain source localization method |
| CN105615877B (en) * | 2016-02-22 | 2018-11-27 | 广州视源电子科技股份有限公司 | Positioning system for epilepsia electroencephalogram signal characteristics |
| CN112674773B (en) * | 2020-12-22 | 2021-12-24 | 北京航空航天大学 | Magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window |
-
2020
- 2020-12-22 CN CN202011526796.6A patent/CN112674773B/en active Active
- 2020-12-22 CN CN202111402038.8A patent/CN113951885B/en active Active
-
2021
- 2021-12-22 WO PCT/CN2021/140340 patent/WO2022135448A1/en active Application Filing
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102791332A (en) * | 2009-11-04 | 2012-11-21 | 代理并代表亚利桑那州立大学的亚利桑那董事会 | Devices and methods for modulating brain activity |
| US20190105517A1 (en) * | 2009-11-04 | 2019-04-11 | Arizona Board Of Regents On Behalf Of Arizona State University | Devices and methods for modulating brain activity |
| US20130144153A1 (en) * | 2011-12-01 | 2013-06-06 | The Regents Of The University Of California | Functional magnetic resonance imaging apparatus and methods |
| US20160022168A1 (en) * | 2014-07-24 | 2016-01-28 | University Of Lethbridge | Brain state dependent therapy for improved neural training and rehabilitation |
| CN110612057A (en) * | 2017-06-07 | 2019-12-24 | 柯惠有限合伙公司 | System and method for detecting stroke |
| CN107518898A (en) * | 2017-08-08 | 2017-12-29 | 北京航空航天大学 | Magneticencephalogram source positioner with beam forming is decomposed based on sensor array |
| CN111801046A (en) * | 2017-11-10 | 2020-10-20 | 勒维斯公司 | Efficacy and/or therapeutic parameter recommendation using individual patient data and therapeutic brain network maps |
| US20200289005A1 (en) * | 2019-03-14 | 2020-09-17 | Ricoh Company, Ltd. | Multi-sphere head model for dipole localization |
| CN111096745A (en) * | 2020-01-02 | 2020-05-05 | 苏州大学 | Steady-state evoked response brain source positioning method based on sparse Bayesian learning |
| CN111209974A (en) * | 2020-01-14 | 2020-05-29 | 华中科技大学鄂州工业技术研究院 | Tensor decomposition-based heterogeneous big data core feature extraction method and system |
| CN111340142A (en) * | 2020-05-14 | 2020-06-26 | 南京慧脑云计算有限公司 | Epilepsia magnetoencephalogram spike automatic detection method and tracing positioning system |
Non-Patent Citations (2)
| Title |
|---|
| 《2DPCA在脑磁图棘波信号检测中的应用》;杨宝山;《计算机工程与应用》;20191231;全文 * |
| 《Mechanisms of physiological and epileptic HFO generation》;John G.R. Jefferys;《Progress in Neurobiology》;20121231;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113951885B (en) | 2022-07-05 |
| CN113951885A (en) | 2022-01-21 |
| WO2022135448A1 (en) | 2022-06-30 |
| CN112674773A (en) | 2021-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112674773B (en) | Magnetoencephalogram source positioning method and device based on Tucker decomposition and ripple time window | |
| van der Meer et al. | Carbon-wire loop based artifact correction outperforms post-processing EEG/fMRI corrections—A validation of a real-time simultaneous EEG/fMRI correction method | |
| Roonizi et al. | A signal decomposition model-based Bayesian framework for ECG components separation | |
| Trapp et al. | On the detection of high frequency correlations in resting state fMRI | |
| Oster et al. | Nonlinear Bayesian filtering for denoising of electrocardiograms acquired in a magnetic resonance environment | |
| EP1872715B1 (en) | Method and system for multi-channel biosignal processing | |
| CN103549950A (en) | Improved difference threshold detection algorithm for mobile ECG (electrocardiogram) monitoring | |
| US12023496B2 (en) | Closed-loop adaptive AC stimulation neural network regulation method and system | |
| EP3679860B1 (en) | Control method and system for filtering out working frequency interference signal | |
| US10912470B2 (en) | System and method for tracking cerebral blood flow in fMRI | |
| US20240164695A1 (en) | Systems and methods for enhanced encoded source imaging | |
| CN108509823A (en) | The detection method and device of QRS complex | |
| CN104905786A (en) | Electrocardiographic artifact on-line removal algorithm | |
| CN108289627A (en) | The application of the method and apparatus and this method or this device that quantified to respiratory sinus arrhythmia | |
| CN110327032A (en) | It is a kind of singly to lead the accurate recognizer of electrocardiosignal PQRST wave joint | |
| Calisto et al. | A new approach for investigating intracranial pressure signal: filtering and morphological features extraction from continuous recording | |
| Raatikainen et al. | Combined spatiotemporal ICA (stICA) for continuous and dynamic lag structure analysis of MREG data | |
| Zafar et al. | Initial dip detection based on both HbO and HbR vector-based phase analysis | |
| Hegde et al. | A review on ECG signal processing and HRV analysis | |
| KR102493239B1 (en) | System for analysising sleep stage and method thereof | |
| Sun et al. | Human intracortical responses to varying electrical stimulation conditions are separable in low-dimensional subspaces | |
| Pecher et al. | Phase estimation of visual evoked responses | |
| Fahruzi et al. | Screening of non-overlapping apnea and non-apnea from single lead ECG-apnea recordings using time-frequency approach | |
| Shama et al. | Hamming filter design for ECG signal detection and processing using co-simulation | |
| CN104007823B (en) | Frequency phase hybrid decoding brain-computer interface method based on point-pass filtering and frequency phase hybrid decoding brain-computer interface device based on point-pass filtering |
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 |