US5020538A - Low noise magnetoencephalogram system and method - Google Patents
Low noise magnetoencephalogram system and method Download PDFInfo
- Publication number
- US5020538A US5020538A US07/435,754 US43575489A US5020538A US 5020538 A US5020538 A US 5020538A US 43575489 A US43575489 A US 43575489A US 5020538 A US5020538 A US 5020538A
- Authority
- US
- United States
- Prior art keywords
- noise
- brain
- channels
- meg
- signals
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/022—Measuring gradient
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S128/00—Surgery
- Y10S128/901—Suppression of noise in electric signal
Definitions
- the present invention relates to medical instruments and more particularly to a low noise magnetoencephalogram (MEG) system.
- MEG magnetoencephalogram
- the human brain produces both electrical and magnetic signals. It is conventional to detect both ongoing spontaneous electrical brain activity and evoked electrical brain activity (EP) by non-invasive electrodes connected to the scalp of the patient using an electroencephalograph (EEG). For example, evoked brain activity (EP) may be stimulated by an external stimulus such as a flashed light.
- EEG electroencephalograph
- the weak magnetic fields at the scalp may be detected by a non-invasive magnetoencephalogram (MEG) whose detecting coils are positioned close to the scalp.
- MEG non-invasive magnetoencephalogram
- the magnetoencephalogram may be used to locate the origin of seizures within the human brain.
- the earliest MEG devices were constructed inside a magnetically shielded room, since the magnetic fields sought to be detected are much weaker than environmental fields.
- the environmental field which is "noise" in the system, changes character over time.
- MEG magnetic shielded room
- SQUID superconducting (S) quantum (QU) interference (I) device
- the SQUID operates at cryogenic temperature at which metal loses its electrical resistivity.
- a metal film is deposited on a cylindrical quartz cylinder having a narrow bridge "weak link".
- a magnetic field will increase the energy, making the superconducitng state unstable so that the weak link converts to its normal state (non-superconducting), allowing one flux quantum to enter.
- the conduction of the ring is monitored by a radio frequency circuit whose sensitivity is increased with a feedback current circuit.
- a “flux transporter” in which a primary coil (“detection coil”), of superconductive wire, is connected to a secondary coil (“input coil”) contained in the superconducting chamber with the SQUID.
- the chamber is a "dewar” vacuum chamber filled with liquid helium.
- five primary coils are placed closely about the patient's head to detect the brain' s magnetic fluxes. Those fluxes create responses in the detection coils which are communicated to the SQUIDS, which produces a voltage proportional to the net magnetic flux on each detection coil.
- the flux transporter and the SQUID constitute a "magnetometer". When used to detect brain activity, the system is called a “magnetoencephalograph” or "MEG”.
- a MEG system may be used in a magnetically unshielded room by forming the detection coil with oppositely directed loops (a "gradiometer"), so that external magnetic flux will be detected in both loops and their effects will be canceled.
- the loops may be arranged, for relative insensitivity to spacially uniform fields and gradients, in a "second derivative gradiometer".
- the brain's spontaneous neuromagnetic activity is weak (10 -12 Tesla) and the neurpmagnetic evoked response is even weaker (10 -13 Tesla). Consequently, it is difficult to detect such activity in the presence of external magnetic fields which produce noise in the system.
- the signal/noise ratio may be improved by using an averaging computer which passes signals which are in synchronism with the stimulus. Since the noise is random, the cumulative average should represent the signals as the random noise cancels itself out. Improvements of signal/noise ratio of 10:1 have been reported. However, that technique is applicable only to evoked response (EP) and not to spontaneous brain activity.
- Efficient and effective noise cancellation is critical to the development of the magnetoencephalogram (MEG) as a clinical tool.
- This noise contamination from environmental magnetic fields is a major factor, even in the latest generation, multichannel, second-order gradiometer of DC SQUID (superconducting quantum interference device) systems.
- SQUID superconducting quantum interference device
- removal of this environmental magnetic field contamination requires using the information from four SQUIDS (which act as environmental noise reference sensors) to filter the corresponding noise out of the seven channels which sense both bonafide brain signals and environmental noise.
- the channels used only to detect environmental noise will be called the noise reference signals or channels, and the channels which detect both neuromagnetic brain signals and environmental noise will be called the brain signals or channels.
- the noise reference channels monitor the three orthogonal components of the ambient (environmental) magnetic field and, usually, the gradient of the field along the axis of the dewar.
- these noise reference signals are scaled through a manual procedure which involves static setting of eight independent parameters, four for the noise amplitude and four for the time derivatives of the noise, for each of the seven brain signal channels. This is a total of 56 manual adjustments.
- Such manual adjustments necessarily are suboptimal and cannot compensate for rapid variations in the relationship between the noise reference and brain signal channels, for example, due to changes in the environment or movement of the dewar from one position to another.
- the present invention is a magnetoencephalogram (MEG) system and method which effectively removes the environmental noise from MEG recordings without manual adjustments.
- the system utilizes adaptive noise cancellation, which continuously and automatically adjust the weights of the reference noise channels to remove the effect of time-varying environmental sources of contamination.
- Another application of the present invention is the removal of other types of contaminants in MEG recordings, such as those arising from magnetic fields generated by muscle artifact, such as movement of the eyeballs.
- MEG magnetoencephalogram
- FIG. 1 is a graph of one MEG brain signal channel of the prior art showing the presence of environmental noise
- FIG. 2 is a graph of one MEG brain signal channel illustrating the effect of the filter system and method of the present invention, the signal channel being connected to a detector in the presence of a subject;
- FIG. 3 is a graph similar to the graph of FIG. 2, but not related to a subject;
- FIG. 4 is a schematic block diagram of the MEG system of the present invention.
- Channel lines 13 and 18 denote N s +L channels and channel lines 19 denote N s channels.
- FIG. 5 is a schematic block diagram of another embodiment of the MEG system of the present invention illustrating the case where noise reduction is not separately applied to frequency bands but is applied to full spectrum signals.
- FIG. 6 is a schematic diagram of yet another embodiment of the MEG system of the present invention illustrating the case where digital bandpass filters are used instead of analog filters.
- Channel lines 13' denote N s +L channels; channel lines 19' denote N s channels.
- FIG. 7 is a diagram schematizing the signal recording situation referenced by the present invention description.
- Three SQUIDs are shown placed about a subject's head; these are the “signal” channels.
- Three other SQUIDs are shown place remote from the subject's head; these are the “reference” channels.
- MEG data are treated using multidimensional adaptive noise cancellation methods.
- most adaptive filter literature has discussed noise cancellation in other contexts, from a single time series, the present invention generalizes to the multichannel case.
- noise contributions to each MEG brain signal channel are estimated with a linear combination of past and present sample values from each of the MEG noise reference channels.
- the problem is formulated as the least squares identification of the multidimensional transfer function between the noise pickup at the MEG noise reference channels and the noise pickup at the MEG brain signal channels.
- covariance between the MEG noise reference channels and the MEG brain signal (plus noise) channels is used to optimally (in the least squares sense) estimate and remove the noise component from the brain signal channels.
- a block computation uses a frame of data to form a matrix for the least squares solution. This is solved using the Householder transformation, explained below, yielding the best linear combination of noise reference channels to estimate the noise in each contaminated brain signal channel.
- k is the time index
- L is the number of reference channels
- M is the prediction memory (i.e., how many samples into the past a current value is regressed on past reference values)
- n 1 (k) is the k th data value for the first noise reference channel
- n 2 (k) is the k th data value for the second noise reference channel.
- N is the data matrix representing a system of equations in which the coefficients are the time series from the noise reference channels
- a is the noise coefficient vector, i.e., the unknown weights
- s is the brain signal vector
- e is the error vector
- the error length is minimized when a is chosen to be the solution of the transformed equations.
- the resulting error length is equal to the square of the lower right-hand element.
- the noise coefficient vector (a) may be used to filter the reference noise channels to best approximate the noise in each signal channel.
- each element of the coefficient vector is itself a vector, but otherwise the basic computation is the same.
- the best translation between the noise reference channels and the brain signal channels is found by a calibration procedure using the above-described least squares method.
- the optimal noise reference weightings in fact, vary significantly over the time of an MEG recording. This means that the least squares calculations must be done during the actual recordings and recalculated frequently, preferably every 0.25 seconds.
- SNR signal-to-noise ratio
- bandpass filters 10 separate the noise reference channels 12 (channels 1 through L) and signal channels 11 (channels 1 through N s ) into contiguous frequency bands 13, e.g., 1-16 Hz ("Low frequency”), 16-56 Hz (“Mid1 frequency”), 56-64 Hz (“Mid2 frequency”), and 64-120 Hz ("High frequency”).
- Filter 10 operates on both the signal and reference channels.
- A/D converters 14a-14d For each signal and noise channel and for each frequency band of these channels, convert to digital data in A/D (analog to digital) converters 14a-14d. Break these data into 50% overlapped time windows. For example, assuming a 0.5 second window, the first time window would span from time 0 to time 0.5 seconds; the second time window would span from time 0.25 to 0.75 seconds, etc..
- the Householder transformation translates the over-determined system of equations (as described in (1) where the system is "over-determined" since the total number of time analyzed, K, is greater than the number of times, M, used for modeling noise) into an equivalent equation. For example, the equation which has the smallest sum of squared errors will give the best solution. This happens to be an upper triangular matrix which is easily solved. Its solution gives the weights of the reference time series which yield the beat approximation to the signal time series.
- the bandlimited noise time series will be weighted by these values and summed to yield an estimate of the signal time series. This estimate is then subtracted from the original time series, using calculation filter system 17. Essentially, this subtracts out the portion of the signal time series which is correlated with the noise time series.
- Householder Transformation Matrix is explained in Chapters 2, 3 and 10 of Lawson and Hanson, SOLVING LEAST-SQUARES PROBLEMS (Prentice-Hall 1974), incorporated by reference herein.
- LS least squares
- e; where ⁇ +1 if v 1 ⁇ 0 and -1 if v 1 ⁇ 0; ##EQU4##
- the subroutine H12 given in Appendix C of the Lawson text implements the steps 1-11 above on a digital computer using the FORTRAN programming language.
- the output of the Householder block is the best estimate of the bandlimited noise component in the original signal.
- the bandlimited noise estimates are then passed to the Calculation Filter System block where they are summed to produce full spectrum estimates for the noise in each signal channel.
- Full spectrum noise estimates for each signal channel for the first half of each time window are then subtracted from corresponding time windows of digitized data from each of the signal channels. For example, using the 0.5 second window example above, the initial window of "clean" signal produced would span times 0 to 0.25 seconds, the second would span times 0.25 to 0.5 seconds, etc.
- FIGS. 1, 2 and 3 are, respectively, samples of auditory evoked (EP) magnetic field and ambient magnetic field data.
- the EP data consist of 5 brain signal channels connected to 5 SQUID brain signal sensor means and 4 noise reference channels connected to 4 SQUID environmental magnetic activity sensor means, 100 recordings are taken each of 2-second long trials from a healthy, right-handed adult male subject hearing a 500 msec auditory tone burst beginning 1.024 seconds into each trial, with the dewar positioned over the right temporal lobe to span both the maxima and minima of the evoked field.
- the data was sampled at 250 HZ and the data was averaged.
- the 9 channels of data were lowpass filtered at 100 HZ, highpass filtered at 0.2 HZ (48 dB/octave at each corner frequency), and filtered every 4 milliseconds. Consequently, components below 0.2 HZ and above 100 HZ were removed. Also prior to analysis, the current standard analog noise cancellation was performed after manual adjustments were made to subtract out noise; line-locked comb filtering was also performed to reduce components at 60 HZ and its harmonics.
- FIGS. 2 and 3 show the dramatic improvement in one MEG channel resulting from application of the filtering method and system of the present invention.
- the improved evoked MEG field response was obtained from the average of 100 trials that had been noise-filtered using 0.5 second windows with a 50% overlap, and 4 consecutive noise removal steps over differing (0-16, 16-56, 56-64, 64-128 HZ) frequency bands. Remanent line-related noise was removed in an effectively perfect notch without any hand-tuning or special notch filters. Low frequency noise was also removed.
- the filtering of the present invention eliminates most of the low- and high-frequency noise.
- the evoked brain magnetic field peak at approximately 95 msec is quite evident.
- FIG. 3 shows the average of 50 segments of MEG recorded with the same apparatus but without a human subject.
- the noise removal has subtracted almost all of the noise.
- the remanent is likely due to the intrinsic residual noise of the SQUID detectors themselves, and not noise from external contaminants.
- FIG. 3 shows the average MEG waveforms for 50 epochs of 3 seconds, without a human subject in the apparatus. Prior to filtering, the data were analog filtered to remove components below 0.2 HZ and above 50 HZ using 48 dB/octave filters. Sampling was at 125 HZ.
- FIG. 3 shows the original waveform from one channel of a 5-channel MEG recording along with the same waveform after noise filtering. The filtering used 0.5 second analysis windows with 50% overlap and four separate frequency bands. The separation into frequency bands tend to separate out the different sources of noise which have different optimal reduction solutions. This is due to the facts that simple least square methods will be dominated by the largest frequency component and the sources of noise will tend to have both radically different spectra and direction of original.
- the embodiment shown in FIG. 4 includes digitization of the outputs of analog bandpass filters as well as the wideband channels.
- special-purpose dedicated digital hardware including digital bandpass filter 10A, is used to perform the noise cancellation.
- software filters have been designed to work within the constraints of a general-purpose laboratory computing environment.
- the signal channels 11, 11', and 11", respectively, are channels 1 through N s ; the environmental noise channels 12, 12', and 12", respectively are channels 1 through L.
- the channels 11 and 12 are connected to bandpass filters 10 whose outputs, on four multi-channel lines 13 (low freq, med 1 freq, and 2 freq, high freq) are to respective analogue-digital converters (A/D) 14a-14d.
- A/D 14a-14d is connected to a Householder least square means 16a-16d in digital computer system 15, and each of the four least square means 16a-16d is connected to calculation filter system (subtractor) 17, to which A/D 14e is also connected.
- A/Ds 14a', 14b' are connected to HLS (Householder Least Squares) 16a', 16b', respectively and A/Ds 14c', 14d' are connected to HLS 16a', 16b', respectively on digital computer system 15', which also includes calculation filter system (subtractor) 17'.
- HLS Householder Least Squares
- A/Ds 14c', 14d' are connected to HLS 16a', 16b', respectively on digital computer system 15', which also includes calculation filter system (subtractor) 17'.
- three SQUIDS 20-22 are connected to brain signal channels (channels 11 of FIG. 4) and three SQUIDS 23-25 are connected to environmental noise reference channels (channels 12 of FIG. 4).
- the system uses one to fourteen magnetic sensor means adjacent the head of the subject to detect neuromagnetic brain signals and three to eight magnetic sensor means to detect the environment.
Abstract
A magnetoencaphalogram (MEG) detects the neuromagnetic brain signals of a subject by bringing a set of magnetic sensors, preferably SQUIDs (superconducting quantum interference device), close to the scalp of the subject. The adverse effect of environmental magnetic noise is reduced by using a second set of magnetic sensors which are placed near the subject. The system includes a bandpass filter, to divide the signals from the various sensors into predetermined frequency bands, and a computer to provide a least squares noise component estimate of the noise from each brain wave sensor and to calculate therefrom a filtered brain wave signal having reduced noise.
Description
This application is a continuation-in-part application based upon application Ser. No. 07/231,385, filed Aug. 12, 1988, now abandoned, and Ser. No. 027,264, filed Mar. 18, 1987, now abandoned.
1. Field of the Invention
The present invention relates to medical instruments and more particularly to a low noise magnetoencephalogram (MEG) system.
2. Description of the Related Art
The human brain produces both electrical and magnetic signals. It is conventional to detect both ongoing spontaneous electrical brain activity and evoked electrical brain activity (EP) by non-invasive electrodes connected to the scalp of the patient using an electroencephalograph (EEG). For example, evoked brain activity (EP) may be stimulated by an external stimulus such as a flashed light.
The weak magnetic fields at the scalp, produced by the flow of electrical current within the brain, may be detected by a non-invasive magnetoencephalogram (MEG) whose detecting coils are positioned close to the scalp. There are indications, for example, that the magnetoencephalogram may be used to locate the origin of seizures within the human brain. The earliest MEG devices were constructed inside a magnetically shielded room, since the magnetic fields sought to be detected are much weaker than environmental fields. The environmental field, which is "noise" in the system, changes character over time. Even when the magnetic detecting instrument is brought up close to the patient's scalp, the brain's magnetic signals are so faint that they may be drowned out, even after compensating for the earth's magnetic field, by the stronger magnetic field of a truck passing by outside the building. Unfortunately, a properly magnetically shielded room is expensive. The room must be large enough for the patient, the operators and the sensitive equipment and may cost over two hundred thousand dollars.
An alternative to the use of a magnetic shielded room is to improve the signal sensitivity and noise suppression of the magnetic field detector. The preferred type of MEG uses a "SQUID" (superconducting (S) quantum (QU) interference (I) device (D)). The SQUID operates at cryogenic temperature at which metal loses its electrical resistivity. In one type of SQUID a metal film is deposited on a cylindrical quartz cylinder having a narrow bridge "weak link". A magnetic field will increase the energy, making the superconducitng state unstable so that the weak link converts to its normal state (non-superconducting), allowing one flux quantum to enter. The conduction of the ring is monitored by a radio frequency circuit whose sensitivity is increased with a feedback current circuit. Generally a "flux transporter" is used in which a primary coil ("detection coil"), of superconductive wire, is connected to a secondary coil ("input coil") contained in the superconducting chamber with the SQUID. For example, the chamber is a "dewar" vacuum chamber filled with liquid helium. For example, five primary coils are placed closely about the patient's head to detect the brain' s magnetic fluxes. Those fluxes create responses in the detection coils which are communicated to the SQUIDS, which produces a voltage proportional to the net magnetic flux on each detection coil. The flux transporter and the SQUID constitute a "magnetometer". When used to detect brain activity, the system is called a "magnetoencephalograph" or "MEG".
It has been suggested that a MEG system may be used in a magnetically unshielded room by forming the detection coil with oppositely directed loops (a "gradiometer"), so that external magnetic flux will be detected in both loops and their effects will be canceled. The loops may be arranged, for relative insensitivity to spacially uniform fields and gradients, in a "second derivative gradiometer".
However, the brain's spontaneous neuromagnetic activity is weak (10-12 Tesla) and the neurpmagnetic evoked response is even weaker (10-13 Tesla). Consequently, it is difficult to detect such activity in the presence of external magnetic fields which produce noise in the system. It has been suggested that the signal/noise ratio may be improved by using an averaging computer which passes signals which are in synchronism with the stimulus. Since the noise is random, the cumulative average should represent the signals as the random noise cancels itself out. Improvements of signal/noise ratio of 10:1 have been reported. However, that technique is applicable only to evoked response (EP) and not to spontaneous brain activity.
Efficient and effective noise cancellation is critical to the development of the magnetoencephalogram (MEG) as a clinical tool. This noise contamination from environmental magnetic fields is a major factor, even in the latest generation, multichannel, second-order gradiometer of DC SQUID (superconducting quantum interference device) systems. For example, for the Biomagnetic Technology Inc. Model 607 magnetometer, removal of this environmental magnetic field contamination requires using the information from four SQUIDS (which act as environmental noise reference sensors) to filter the corresponding noise out of the seven channels which sense both bonafide brain signals and environmental noise. The channels used only to detect environmental noise will be called the noise reference signals or channels, and the channels which detect both neuromagnetic brain signals and environmental noise will be called the brain signals or channels. The noise reference channels monitor the three orthogonal components of the ambient (environmental) magnetic field and, usually, the gradient of the field along the axis of the dewar. In the current commercial systems, these noise reference signals are scaled through a manual procedure which involves static setting of eight independent parameters, four for the noise amplitude and four for the time derivatives of the noise, for each of the seven brain signal channels. This is a total of 56 manual adjustments. Such manual adjustments necessarily are suboptimal and cannot compensate for rapid variations in the relationship between the noise reference and brain signal channels, for example, due to changes in the environment or movement of the dewar from one position to another.
The book, H. Weinberg et al, BIOMAGNETISM: APPLICATIONS AND THEORY (Pergamon Press 1985) contains a chapter, Williamson et al, FIVE CHANNEL SQUID INSTALLATION FOR UNSHIELDED NEUROMAGNETIC MEASUREMENTS (pgs. 46-51, incorporated by reference herein). In Williamson nine SQUIDS are used in an unshielded environment. Five SQUIDS, having second order gradiometers, monitor brain activity (5-channel magnetometer) and four SQUIDS monitor three components of the ambient field and one component of the gradient. The four SQUIDS attempt to cancel background (environmental) noise by an analog noise cancellation system requiring 40 static adjustments for the 5 brain signal channel magnetometer. The results from this cumbersome procedure are that MEG recordings are often unusable because of overwhelming residual contamination by environmental noise sources. In practice, users of MEG devices often must purchase magnetically shielded recording rooms costing about as much as the MEG instrument itself.
The present invention is a magnetoencephalogram (MEG) system and method which effectively removes the environmental noise from MEG recordings without manual adjustments. The system utilizes adaptive noise cancellation, which continuously and automatically adjust the weights of the reference noise channels to remove the effect of time-varying environmental sources of contamination. Another application of the present invention is the removal of other types of contaminants in MEG recordings, such as those arising from magnetic fields generated by muscle artifact, such as movement of the eyeballs.
It is an objective of the present invention to provide a low noise magnetoencephalogram (MEG) system and method which sufficiently cancels environmental magnetic noise so that the MEG system may be utilized in a nonmagnetically shielded room.
It is a further objective of the present invention to provide such an MEG system which will automatically and rapidly compensates for changes in the magnetic environment so that the manual adjustments are not necessary.
It is a further objective of the present invention to provide such an MEG system which will not greatly add to the costs or complexity of the MEG system, but instead may utilize the existing number of channels in commercially available systems and may utilize available digital computer systems.
It is a further objective of the present invention to provide such an MEG system which will provide a set of "clean", i.e., filtered, MEG signals whose signal/noise ratio is improved relative to alternative systems and methods.
Other objectives of the present invention will be apparent from the following detailed description of the invention, taken in conjunction with the accompanying drawings.
In the drawings:
FIG. 1 is a graph of one MEG brain signal channel of the prior art showing the presence of environmental noise;
FIG. 2 is a graph of one MEG brain signal channel illustrating the effect of the filter system and method of the present invention, the signal channel being connected to a detector in the presence of a subject;
FIG. 3 is a graph similar to the graph of FIG. 2, but not related to a subject;
FIG. 4 is a schematic block diagram of the MEG system of the present invention. Channel lines 13 and 18 denote Ns +L channels and channel lines 19 denote Ns channels.
FIG. 5 is a schematic block diagram of another embodiment of the MEG system of the present invention illustrating the case where noise reduction is not separately applied to frequency bands but is applied to full spectrum signals.
FIG. 6 is a schematic diagram of yet another embodiment of the MEG system of the present invention illustrating the case where digital bandpass filters are used instead of analog filters. Channel lines 13' denote Ns +L channels; channel lines 19' denote Ns channels.
FIG. 7 is a diagram schematizing the signal recording situation referenced by the present invention description. Three SQUIDs are shown placed about a subject's head; these are the "signal" channels. Three other SQUIDs are shown place remote from the subject's head; these are the "reference" channels.
In accordance with the present invention, MEG data are treated using multidimensional adaptive noise cancellation methods. Although most adaptive filter literature has discussed noise cancellation in other contexts, from a single time series, the present invention generalizes to the multichannel case. In the present multidimensional adaptive filter, noise contributions to each MEG brain signal channel are estimated with a linear combination of past and present sample values from each of the MEG noise reference channels.
The problem is formulated as the least squares identification of the multidimensional transfer function between the noise pickup at the MEG noise reference channels and the noise pickup at the MEG brain signal channels. In this method, covariance between the MEG noise reference channels and the MEG brain signal (plus noise) channels is used to optimally (in the least squares sense) estimate and remove the noise component from the brain signal channels. A block computation uses a frame of data to form a matrix for the least squares solution. This is solved using the Householder transformation, explained below, yielding the best linear combination of noise reference channels to estimate the noise in each contaminated brain signal channel.
Given a system of prediction equations (one equation for each time point in the frame) as follows: ##EQU1## where: k is the time index; L is the number of reference channels; M is the prediction memory (i.e., how many samples into the past a current value is regressed on past reference values); n1 (k) is the kth data value for the first noise reference channel; and n2 (k) is the kth data value for the second noise reference channel. Or, expressed as a matrix equation
Na-s=e (2)
where N is the data matrix representing a system of equations in which the coefficients are the time series from the noise reference channels, a is the noise coefficient vector, i.e., the unknown weights, s is the brain signal vector, and e is the error vector.
The least squares criterion is met by that choice of a which minimizes the Euclidean length of e, written as
∥e∥=∥Na-s∥ (3)
If Q is an orthogonal matrix, then
∥QNa=Qs∥=∥Na-s∥ (4)
If Q is further chosen to reduce N|s--the matrix N augmented by column vector s--to upper triangular form, the error length is minimized when a is chosen to be the solution of the transformed equations. The resulting error length is equal to the square of the lower right-hand element.
After such a computation, the noise coefficient vector (a) may be used to filter the reference noise channels to best approximate the noise in each signal channel. For the multidimensional case, each element of the coefficient vector is itself a vector, but otherwise the basic computation is the same.
In principle, the best translation between the noise reference channels and the brain signal channels is found by a calibration procedure using the above-described least squares method. However, the optimal noise reference weightings, in fact, vary significantly over the time of an MEG recording. This means that the least squares calculations must be done during the actual recordings and recalculated frequently, preferably every 0.25 seconds. Further, since the actual signal-to-noise ratio (SNR) in MEG recordings is both frequency dependent and unknown, and since the nature of the measured noise is also unknown in each frequency band (i.e., how much is independent instrumental noise, how much is dependent instrumental noise, how much is correlated with the reference channels, etc.), consecutive passes through the data over discrete frequency bands give the best overall noise cancellation.
The basic steps, in sequence, of the noise reduction method, as illustrated in FIG. 4, are:
(1) Using the bandpass filters 10, separate the noise reference channels 12 (channels 1 through L) and signal channels 11 (channels 1 through Ns) into contiguous frequency bands 13, e.g., 1-16 Hz ("Low frequency"), 16-56 Hz ("Mid1 frequency"), 56-64 Hz ("Mid2 frequency"), and 64-120 Hz ("High frequency"). Filter 10 operates on both the signal and reference channels.
(2) For each signal and noise channel and for each frequency band of these channels, convert to digital data in A/D (analog to digital) converters 14a-14d. Break these data into 50% overlapped time windows. For example, assuming a 0.5 second window, the first time window would span from time 0 to time 0.5 seconds; the second time window would span from time 0.25 to 0.75 seconds, etc..
(3) Input each time window to the Householder blocks for analysis. In the computer system 15, the best least-squares translation function between the filtered noise reference channels 12 and the filtered signal channels is coupled consecutively for each time window using the multidimensional Householder blocks 16a-16d.
The Householder transformation translates the over-determined system of equations (as described in (1) where the system is "over-determined" since the total number of time analyzed, K, is greater than the number of times, M, used for modeling noise) into an equivalent equation. For example, the equation which has the smallest sum of squared errors will give the best solution. This happens to be an upper triangular matrix which is easily solved. Its solution gives the weights of the reference time series which yield the beat approximation to the signal time series.
Once the weights described above are calculated, the bandlimited noise time series will be weighted by these values and summed to yield an estimate of the signal time series. This estimate is then subtracted from the original time series, using calculation filter system 17. Essentially, this subtracts out the portion of the signal time series which is correlated with the noise time series.
The Householder Transformation Matrix is explained in Chapters 2, 3 and 10 of Lawson and Hanson, SOLVING LEAST-SQUARES PROBLEMS (Prentice-Hall 1974), incorporated by reference herein.
In general, the least squares ("LS") problem is solved based on a certain decomposition called an "orthogonal decomposition of A". A is A=HRKT where H is an mXm orthogonal matrix, R is an mXm matrix of the form ##EQU2## R11 is a kXk submatrix of rank k and k is an nXn orthogonal matrix.
The Householder transformtion matrix for a non-zero vector V is ##EQU3## where Q is an orthogonal matrix, Im is a matrix; u=v+σ|v|e; where σ=+1 if v1 ≧0 and -1 if v1 <0; ##EQU4##
One suitable series of steps is as follows, where p, l, m, and w are integers; v is a vector (m-vector) if w>0 then array contains m-vectors, cj, j=1, . . . , w. The following steps compute vector u, number b, vector y=Qv, and if w>0 the vectors cj =Qcj ; where j=1, . . . , W.
______________________________________ Step Description ______________________________________ ##STR1## 2 If v.sub.p > 0, set s: - s 3 Set h: = v.sub.p - s, v.sub.p : = s. 4 Comment: The construction of the transformation is complete. At Step 5 the application of the transformation to the vectors c.sub.j begins. 5 Set b: = v.sub.p h. 6 If b = 0 or W = 0, go to Step 11 7 for j: = 1, . . . , w, do Steps 8-10 8 ##STR2## 9 Set c.sub.pj : = c.sub.pj + sh. 10 For i: l, . . . , m, set c.sub.ij : = c.sub.ij + sV.sub.i 11 completed. ______________________________________
The above algorithm is used to reduce the N matrix augmented by the column vector s of the matrix equation Na-s=e. The subroutine H12 given in Appendix C of the Lawson text implements the steps 1-11 above on a digital computer using the FORTRAN programming language.
The output of the Householder block is the best estimate of the bandlimited noise component in the original signal.
(4) The bandlimited noise estimates are then passed to the Calculation Filter System block where they are summed to produce full spectrum estimates for the noise in each signal channel. Full spectrum noise estimates for each signal channel for the first half of each time window are then subtracted from corresponding time windows of digitized data from each of the signal channels. For example, using the 0.5 second window example above, the initial window of "clean" signal produced would span times 0 to 0.25 seconds, the second would span times 0.25 to 0.5 seconds, etc.
Evidence of the efficacy of the filter is shown in FIGS. 1, 2 and 3. These are, respectively, samples of auditory evoked (EP) magnetic field and ambient magnetic field data. The EP data consist of 5 brain signal channels connected to 5 SQUID brain signal sensor means and 4 noise reference channels connected to 4 SQUID environmental magnetic activity sensor means, 100 recordings are taken each of 2-second long trials from a healthy, right-handed adult male subject hearing a 500 msec auditory tone burst beginning 1.024 seconds into each trial, with the dewar positioned over the right temporal lobe to span both the maxima and minima of the evoked field. The data was sampled at 250 HZ and the data was averaged. Prior to analysis, the 9 channels of data were lowpass filtered at 100 HZ, highpass filtered at 0.2 HZ (48 dB/octave at each corner frequency), and filtered every 4 milliseconds. Consequently, components below 0.2 HZ and above 100 HZ were removed. Also prior to analysis, the current standard analog noise cancellation was performed after manual adjustments were made to subtract out noise; line-locked comb filtering was also performed to reduce components at 60 HZ and its harmonics.
In FIGS. 1-3, on the X-axis, 0.2=200 msec post-stimulus. The Y-axis is in arbitrary units proportional to femtotesla where 1 femtotesla=10-15 tesla. Each figure shows one channel of a 5-channel MEG recording. FIG. 1 is with averaging; but without the noise reduction system and method of the present invention. FIGS. 2 and 3 show the dramatic improvement in one MEG channel resulting from application of the filtering method and system of the present invention. The improved evoked MEG field response was obtained from the average of 100 trials that had been noise-filtered using 0.5 second windows with a 50% overlap, and 4 consecutive noise removal steps over differing (0-16, 16-56, 56-64, 64-128 HZ) frequency bands. Remanent line-related noise was removed in an effectively perfect notch without any hand-tuning or special notch filters. Low frequency noise was also removed.
As seen in FIG. 2, the filtering of the present invention eliminates most of the low- and high-frequency noise. The evoked brain magnetic field peak at approximately 95 msec is quite evident.
FIG. 3 shows the average of 50 segments of MEG recorded with the same apparatus but without a human subject. The noise removal has subtracted almost all of the noise. The remanent is likely due to the intrinsic residual noise of the SQUID detectors themselves, and not noise from external contaminants.
FIG. 3 shows the average MEG waveforms for 50 epochs of 3 seconds, without a human subject in the apparatus. Prior to filtering, the data were analog filtered to remove components below 0.2 HZ and above 50 HZ using 48 dB/octave filters. Sampling was at 125 HZ. FIG. 3 shows the original waveform from one channel of a 5-channel MEG recording along with the same waveform after noise filtering. The filtering used 0.5 second analysis windows with 50% overlap and four separate frequency bands. The separation into frequency bands tend to separate out the different sources of noise which have different optimal reduction solutions. This is due to the facts that simple least square methods will be dominated by the largest frequency component and the sources of noise will tend to have both radically different spectra and direction of original.
Since the frequency bandsplitting is a significant fraction of the computation, the embodiment shown in FIG. 4 includes digitization of the outputs of analog bandpass filters as well as the wideband channels. In another embodiment, shown in FIG. 5, special-purpose dedicated digital hardware, including digital bandpass filter 10A, is used to perform the noise cancellation. In a third embodiment, shown in FIG. 6, software filters have been designed to work within the constraints of a general-purpose laboratory computing environment.
In FIGS. 4, 5, and 6 the signal channels 11, 11', and 11", respectively, are channels 1 through Ns ; the environmental noise channels 12, 12', and 12", respectively are channels 1 through L. In FIG. 4, the channels 11 and 12 are connected to bandpass filters 10 whose outputs, on four multi-channel lines 13 (low freq, med 1 freq, and 2 freq, high freq) are to respective analogue-digital converters (A/D) 14a-14d. Each A/D 14a-14d is connected to a Householder least square means 16a-16d in digital computer system 15, and each of the four least square means 16a-16d is connected to calculation filter system (subtractor) 17, to which A/D 14e is also connected. In FIG. 5 the A/Ds 14a', 14b' are connected to HLS (Householder Least Squares) 16a', 16b', respectively and A/Ds 14c', 14d' are connected to HLS 16a', 16b', respectively on digital computer system 15', which also includes calculation filter system (subtractor) 17'.
In the embodiment of FIG. 7 three SQUIDS 20-22 are connected to brain signal channels (channels 11 of FIG. 4) and three SQUIDS 23-25 are connected to environmental noise reference channels (channels 12 of FIG. 4).
Preferably the system uses one to fourteen magnetic sensor means adjacent the head of the subject to detect neuromagnetic brain signals and three to eight magnetic sensor means to detect the environment.
Claims (10)
1. A low noise magnetoencephalogram (MEG) system in an environment having magnetic activity, the system being for the detection of brain neuromagnetic activity of a human subject, including:
a first plurality of magnetic sensor means at least part of each sensor means being adapted to be positioned adjacent the head of a subject to detect neuromagnetic brain signals; a plurality of brain signal channels, each brain signal channel connected to one of said first magnetic sensor means;
a second plurality of magnetic sensor means to detect signals from environmental magnetic activity; a plurality of environmental noise reference channels, each reference channel being connected to one of said second magnetic sensor means;
a bandpass filter means, connected to the brain signal channels and the reference channels to separate said brain signals and environmental signals into a plurality of predetermined frequency bands and reject said signals not within said predetermined frequency bands;
computer transformation means which act separately on each of said frequency bands to provide a least squares noise component estimate of the noise component present in each brain signal channel;
computer calculation means to algebraically subtract said least square noise component estimate from the signals in each brain wave channel thereby providing a filtered brain signal having reduced noise.
2. A magnetoencephalogram (MEG) system as in claim 1 wherein the transformation means uses a Householder transformation matrix to provide the least square estimate.
3. A magnetoencephalogram (MEG) system as in claim 1 the system comprises at least 5 brain signal channels and 4 reference channels.
4. A magnetoencephalogram (MEG) system as in claim 1 wherein said bandpass filter means include at least 4 bandpass filters.
5. A magnetoencephalogram (MEG) system as in claim 1 wherein each of the sensor means is a SQUID (superconducting quantum interference device).
6. A noise reduction method in a magnetoencephalogram (MEG) for the detection of human brain neuromagnetic activity of a subject in an environment having noise magnetic activity including:
positioning a plurality of magnetic sensor means adjacent to the head of a subject to detect neuromagnetic brain signals; connecting a brain signal channel to each of said brain signal magnetic sensor means;
positioning a plurality of magnetic sensor means to detect the environmental magnetic activity; connecting a reference channel to each of the environmental magnetic senor means;
separating the signals from each of the channels into a plurality of frequency bands using a bandpass filter means connected to the brain signal channels and the reference channels;
automatically providing a least squares estimate of the noise component present in each brain signal channel and applying said least squares estimate to each of said frequency bands;
automatically algebraically subtracting said noise component estimate from the signals in each brain wave channel thereby providing a filtered brain signal having reduced noise.
7. A noise reduction method as in claim 6 wherein 3 brain signal channels are connected to 3 brain signal magnetic sensor means and 3 reference channels are connected to 3 environment magnetic sensor means.
8. A noise reduction method as in claim 6 and using a Householder transformation matrix to provide the least squares estimate.
9. A noise reduction method as in claim 6 wherein the frequencies within the bandpass range are divided into at least 4 frequency bands.
10. A noise reduction method as in claim 6 and including the step of locating the sensor means in a non-magnetically shielded room.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US07/435,754 US5020538A (en) | 1988-08-12 | 1989-11-14 | Low noise magnetoencephalogram system and method |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US23138588A | 1988-08-12 | 1988-08-12 | |
US07/435,754 US5020538A (en) | 1988-08-12 | 1989-11-14 | Low noise magnetoencephalogram system and method |
Related Parent Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US07027264 Continuation-In-Part | 1987-03-18 | ||
US23138588A Continuation-In-Part | 1988-08-12 | 1988-08-12 |
Publications (1)
Publication Number | Publication Date |
---|---|
US5020538A true US5020538A (en) | 1991-06-04 |
Family
ID=26925077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US07/435,754 Expired - Fee Related US5020538A (en) | 1988-08-12 | 1989-11-14 | Low noise magnetoencephalogram system and method |
Country Status (1)
Country | Link |
---|---|
US (1) | US5020538A (en) |
Cited By (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1994012100A1 (en) * | 1992-11-30 | 1994-06-09 | Risto Ilmoniemi | Method and apparatus for separating the different components of evoked response and spontaneous activity brain signals as well as of signals measured from the heart |
US5440232A (en) * | 1993-12-06 | 1995-08-08 | The United States Of America As Represented By The Secretary Of The Navy | System for monitoring and analyzing field energy exposure |
US5458142A (en) * | 1993-03-19 | 1995-10-17 | Farmer; Edward J. | Device for monitoring a magnetic field emanating from an organism |
WO1996014796A1 (en) * | 1994-11-14 | 1996-05-23 | Cambridge Heart, Inc. | Using related signals to reduce ecg noise |
US5557199A (en) * | 1994-04-29 | 1996-09-17 | The United States Of America As Represented By The Department Of Health And Human Services | Magnetic resonance monitor |
US5574369A (en) * | 1994-08-19 | 1996-11-12 | Hibbs; Andrew D. | Detection and communications device employing stochastic resonance |
US5633583A (en) * | 1995-06-07 | 1997-05-27 | Gas Research Institute | Magnetic telescope with enhanced noise suppression |
EP0710849A3 (en) * | 1994-11-07 | 1997-11-05 | Philips Patentverwaltung GmbH | Method to determine spatial field distribution |
US5697958A (en) * | 1995-06-07 | 1997-12-16 | Intermedics, Inc. | Electromagnetic noise detector for implantable medical devices |
US5737433A (en) * | 1996-01-16 | 1998-04-07 | Gardner; William A. | Sound environment control apparatus |
WO1998040757A1 (en) * | 1997-03-10 | 1998-09-17 | The Secretary Of State For Defence | Magnetic gradiometer |
GB2337128A (en) * | 1997-03-10 | 1999-11-10 | Secr Defence | Magnetic gradiometer |
EP0982597A2 (en) * | 1998-08-28 | 2000-03-01 | Neuromag Oy | Method and apparatus for eliminating background interference signals from multichannel detector arrays |
US6795724B2 (en) | 2002-02-19 | 2004-09-21 | Mark Bradford Hogan | Color-based neurofeedback |
US6853186B2 (en) | 2002-01-15 | 2005-02-08 | National University Of Singapore | Variable permeability magnetic field sensor and method |
US20050146326A1 (en) * | 2002-01-15 | 2005-07-07 | Xiaoping Li | Variable permeability magnetic field sensor and method |
US20060031038A1 (en) * | 2003-03-14 | 2006-02-09 | Elekta Neuromag Oy | Method and system for processing a multi-channel measurement of magnetic fields |
US20060100525A1 (en) * | 2004-10-21 | 2006-05-11 | Takusige Katura | Biological optical measurement system |
KR100583767B1 (en) * | 2003-10-21 | 2006-05-26 | 한국표준과학연구원 | Elimination of Mimic Signals from the Signal-processing Multichannel Stimulus-evoked Magnetoencephalogram |
US20060119794A1 (en) * | 2004-12-03 | 2006-06-08 | Hillis W D | Vision modification with reflected image |
US20060122530A1 (en) * | 2004-12-03 | 2006-06-08 | Goodall Eleanor V | Adjustable lens system with neural-based control |
US20060119793A1 (en) * | 2004-12-03 | 2006-06-08 | Hillis W D | Temporal vision modification |
US20060122531A1 (en) * | 2004-12-03 | 2006-06-08 | Goodall Eleanor V | Method and system for adaptive vision modification |
US20060146281A1 (en) * | 2004-12-03 | 2006-07-06 | Goodall Eleanor V | Method and system for vision enhancement |
EP1727072A1 (en) * | 2005-05-25 | 2006-11-29 | The Babraham Institute | Signal processing, transmission, data storage and representation |
US20070010757A1 (en) * | 2004-12-03 | 2007-01-11 | Searete Llc | Method and system for adaptive vision modification |
US20070019272A1 (en) * | 2004-12-03 | 2007-01-25 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Vision modification with reflected image |
US20070019279A1 (en) * | 2004-12-03 | 2007-01-25 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Adjustable lens system with neural-based control |
US20070019157A1 (en) * | 2004-12-03 | 2007-01-25 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Vision modification with reflected image |
US20070028931A1 (en) * | 2004-12-03 | 2007-02-08 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Temporal vision modification |
US20070060058A1 (en) * | 1993-07-23 | 2007-03-15 | Shattil Steven J | Cancellation Systems for Multicarrier Transceiver Arrays |
US20080137031A1 (en) * | 2004-12-03 | 2008-06-12 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Vision modification with reflected image |
US20090093709A1 (en) * | 2007-05-18 | 2009-04-09 | Beth Israel Deaconess Medical Center, Inc. | Noise reduction system and methods for magnetic resonance imaging |
US20090195751A1 (en) * | 2004-12-03 | 2009-08-06 | Searete Llc | Vision modification with reflected image |
US20100085837A1 (en) * | 2008-10-03 | 2010-04-08 | Kambiz Iranpour | Harmonic attenuation using multiple sweep rates |
US20110150253A1 (en) * | 2009-12-22 | 2011-06-23 | Siemens Medical Instruments Pte. Ltd. | Method and device for setting a hearing device by detecting listening effort |
US9155483B2 (en) | 2004-12-03 | 2015-10-13 | The Invention Science Fund I, Llc | Vision modification with reflected image |
US20180333062A1 (en) * | 2017-05-19 | 2018-11-22 | Taishi WATANABE | Biomagnetic measurement method, biomagnetic measuring device, and biomagnetic measuring system |
JP2018192236A (en) * | 2017-05-19 | 2018-12-06 | 株式会社リコー | Biomagnetic measurement method, biomagnetic measuring device, and biomagnetic measuring system |
US11273283B2 (en) | 2017-12-31 | 2022-03-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
US11364361B2 (en) | 2018-04-20 | 2022-06-21 | Neuroenhancement Lab, LLC | System and method for inducing sleep by transplanting mental states |
US11452839B2 (en) | 2018-09-14 | 2022-09-27 | Neuroenhancement Lab, LLC | System and method of improving sleep |
US11717686B2 (en) | 2017-12-04 | 2023-08-08 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to facilitate learning and performance |
US11723579B2 (en) | 2017-09-19 | 2023-08-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement |
US11786694B2 (en) | 2019-05-24 | 2023-10-17 | NeuroLight, Inc. | Device, method, and app for facilitating sleep |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3697870A (en) * | 1971-09-24 | 1972-10-10 | Burroughs Corp | Digitally nulled magnetic detection system |
US3980076A (en) * | 1974-10-02 | 1976-09-14 | The Board Of Trustees Of Leland Stanford Junior University | Method for measuring externally of the human body magnetic susceptibility changes |
US4324255A (en) * | 1980-03-07 | 1982-04-13 | Barach John P | Method and apparatus for measuring magnetic fields and electrical currents in biological and other systems |
US4386318A (en) * | 1980-09-26 | 1983-05-31 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence | Method and apparatus to compensate a gradiometer having first and second unwanted terms |
US4389612A (en) * | 1980-06-17 | 1983-06-21 | S.H.E. Corporation | Apparatus for reducing low frequency noise in dc biased SQUIDS |
US4613817A (en) * | 1983-07-05 | 1986-09-23 | Siemens Aktiengesellschaft | Superconducting gradiometer coil system for an apparatus for the multi-channel measurement of weak nonstationary magnetic fields |
US4675606A (en) * | 1985-07-29 | 1987-06-23 | Suman Ganguly | Magnetometers for detecting metallic objects in earth's magnetic field |
-
1989
- 1989-11-14 US US07/435,754 patent/US5020538A/en not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3697870A (en) * | 1971-09-24 | 1972-10-10 | Burroughs Corp | Digitally nulled magnetic detection system |
US3980076A (en) * | 1974-10-02 | 1976-09-14 | The Board Of Trustees Of Leland Stanford Junior University | Method for measuring externally of the human body magnetic susceptibility changes |
US4324255A (en) * | 1980-03-07 | 1982-04-13 | Barach John P | Method and apparatus for measuring magnetic fields and electrical currents in biological and other systems |
US4389612A (en) * | 1980-06-17 | 1983-06-21 | S.H.E. Corporation | Apparatus for reducing low frequency noise in dc biased SQUIDS |
US4386318A (en) * | 1980-09-26 | 1983-05-31 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence | Method and apparatus to compensate a gradiometer having first and second unwanted terms |
US4613817A (en) * | 1983-07-05 | 1986-09-23 | Siemens Aktiengesellschaft | Superconducting gradiometer coil system for an apparatus for the multi-channel measurement of weak nonstationary magnetic fields |
US4675606A (en) * | 1985-07-29 | 1987-06-23 | Suman Ganguly | Magnetometers for detecting metallic objects in earth's magnetic field |
Cited By (81)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1994012100A1 (en) * | 1992-11-30 | 1994-06-09 | Risto Ilmoniemi | Method and apparatus for separating the different components of evoked response and spontaneous activity brain signals as well as of signals measured from the heart |
US5458142A (en) * | 1993-03-19 | 1995-10-17 | Farmer; Edward J. | Device for monitoring a magnetic field emanating from an organism |
US20070060058A1 (en) * | 1993-07-23 | 2007-03-15 | Shattil Steven J | Cancellation Systems for Multicarrier Transceiver Arrays |
US5440232A (en) * | 1993-12-06 | 1995-08-08 | The United States Of America As Represented By The Secretary Of The Navy | System for monitoring and analyzing field energy exposure |
US5557199A (en) * | 1994-04-29 | 1996-09-17 | The United States Of America As Represented By The Department Of Health And Human Services | Magnetic resonance monitor |
US5574369A (en) * | 1994-08-19 | 1996-11-12 | Hibbs; Andrew D. | Detection and communications device employing stochastic resonance |
EP0710849A3 (en) * | 1994-11-07 | 1997-11-05 | Philips Patentverwaltung GmbH | Method to determine spatial field distribution |
US5704365A (en) * | 1994-11-14 | 1998-01-06 | Cambridge Heart, Inc. | Using related signals to reduce ECG noise |
WO1996014796A1 (en) * | 1994-11-14 | 1996-05-23 | Cambridge Heart, Inc. | Using related signals to reduce ecg noise |
US5697958A (en) * | 1995-06-07 | 1997-12-16 | Intermedics, Inc. | Electromagnetic noise detector for implantable medical devices |
US5845500A (en) * | 1995-06-07 | 1998-12-08 | Podney; Walter N. | Magnetic telescope with enhanced noise suppression |
US5633583A (en) * | 1995-06-07 | 1997-05-27 | Gas Research Institute | Magnetic telescope with enhanced noise suppression |
US5737433A (en) * | 1996-01-16 | 1998-04-07 | Gardner; William A. | Sound environment control apparatus |
WO1998040757A1 (en) * | 1997-03-10 | 1998-09-17 | The Secretary Of State For Defence | Magnetic gradiometer |
GB2337128A (en) * | 1997-03-10 | 1999-11-10 | Secr Defence | Magnetic gradiometer |
GB2337128B (en) * | 1997-03-10 | 2001-12-05 | Secr Defence | Magnetic gradiometer |
US6339328B1 (en) | 1997-03-10 | 2002-01-15 | The Secretary Of State For Defence | Magnetic gradiometer incorporating global feedback |
EP0982597A2 (en) * | 1998-08-28 | 2000-03-01 | Neuromag Oy | Method and apparatus for eliminating background interference signals from multichannel detector arrays |
EP0982597A3 (en) * | 1998-08-28 | 2009-01-28 | ELEKTA AB (publ.) | Method and apparatus for eliminating background interference signals from multichannel detector arrays |
US6538436B1 (en) | 1998-08-28 | 2003-03-25 | Neuromag Oy | Method and apparatus for eliminating background interference signals from multichannel signal measurements |
US20050146326A1 (en) * | 2002-01-15 | 2005-07-07 | Xiaoping Li | Variable permeability magnetic field sensor and method |
US6853186B2 (en) | 2002-01-15 | 2005-02-08 | National University Of Singapore | Variable permeability magnetic field sensor and method |
US7196514B2 (en) | 2002-01-15 | 2007-03-27 | National University Of Singapore | Multi-conductive ferromagnetic core, variable permeability field sensor and method |
US6795724B2 (en) | 2002-02-19 | 2004-09-21 | Mark Bradford Hogan | Color-based neurofeedback |
US7035686B2 (en) | 2002-02-19 | 2006-04-25 | Mark Bradford Hogan | Color-based neurofeedback |
US20040210156A1 (en) * | 2002-02-19 | 2004-10-21 | Hogan Mark Bradford | Color-based neurofeedback |
US20060031038A1 (en) * | 2003-03-14 | 2006-02-09 | Elekta Neuromag Oy | Method and system for processing a multi-channel measurement of magnetic fields |
US7463024B2 (en) * | 2003-03-14 | 2008-12-09 | Elekta Neuromag Oy | Method and device for processing a multi-channel measurement of magnetic fields |
KR100583767B1 (en) * | 2003-10-21 | 2006-05-26 | 한국표준과학연구원 | Elimination of Mimic Signals from the Signal-processing Multichannel Stimulus-evoked Magnetoencephalogram |
US20060100525A1 (en) * | 2004-10-21 | 2006-05-11 | Takusige Katura | Biological optical measurement system |
US7565192B2 (en) * | 2004-10-21 | 2009-07-21 | Hitachi, Ltd. | Biological optical measurement system |
US7334894B2 (en) | 2004-12-03 | 2008-02-26 | Searete, Llc | Temporal vision modification |
US20100103373A1 (en) * | 2004-12-03 | 2010-04-29 | Searete Llc | Vision modification with reflected image |
US20070019272A1 (en) * | 2004-12-03 | 2007-01-25 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Vision modification with reflected image |
US20070019279A1 (en) * | 2004-12-03 | 2007-01-25 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Adjustable lens system with neural-based control |
US20070019157A1 (en) * | 2004-12-03 | 2007-01-25 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Vision modification with reflected image |
US20070028931A1 (en) * | 2004-12-03 | 2007-02-08 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Temporal vision modification |
US10409088B2 (en) | 2004-12-03 | 2019-09-10 | Gearbox, Llc | Adjustable contact lens system |
US20060146281A1 (en) * | 2004-12-03 | 2006-07-06 | Goodall Eleanor V | Method and system for vision enhancement |
US7334892B2 (en) | 2004-12-03 | 2008-02-26 | Searete Llc | Method and system for vision enhancement |
US20060122531A1 (en) * | 2004-12-03 | 2006-06-08 | Goodall Eleanor V | Method and system for adaptive vision modification |
US7344244B2 (en) | 2004-12-03 | 2008-03-18 | Searete, Llc | Adjustable lens system with neural-based control |
US7350919B2 (en) | 2004-12-03 | 2008-04-01 | Searete Llc | Vision modification with reflected image |
US20080137031A1 (en) * | 2004-12-03 | 2008-06-12 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Vision modification with reflected image |
US7390088B2 (en) | 2004-12-03 | 2008-06-24 | Searete Llc | Adjustable lens system with neural-based control |
US20080161673A1 (en) * | 2004-12-03 | 2008-07-03 | Searete Llc. | Method and system for adaptive vision modification |
US20060119793A1 (en) * | 2004-12-03 | 2006-06-08 | Hillis W D | Temporal vision modification |
US7470027B2 (en) | 2004-12-03 | 2008-12-30 | Searete Llc | Temporal vision modification |
US20060122530A1 (en) * | 2004-12-03 | 2006-06-08 | Goodall Eleanor V | Adjustable lens system with neural-based control |
US7486988B2 (en) | 2004-12-03 | 2009-02-03 | Searete Llc | Method and system for adaptive vision modification |
US9459470B2 (en) | 2004-12-03 | 2016-10-04 | Gearbox, Llc | Vision modification with reflected image |
US20060119794A1 (en) * | 2004-12-03 | 2006-06-08 | Hillis W D | Vision modification with reflected image |
US20090195751A1 (en) * | 2004-12-03 | 2009-08-06 | Searete Llc | Vision modification with reflected image |
US7594727B2 (en) | 2004-12-03 | 2009-09-29 | Searete Llc | Vision modification with reflected image |
US7656569B2 (en) | 2004-12-03 | 2010-02-02 | Searete Llc | Vision modification with reflected image |
US9155483B2 (en) | 2004-12-03 | 2015-10-13 | The Invention Science Fund I, Llc | Vision modification with reflected image |
US20070010757A1 (en) * | 2004-12-03 | 2007-01-11 | Searete Llc | Method and system for adaptive vision modification |
US20100177279A1 (en) * | 2004-12-03 | 2010-07-15 | Searete Llc | Vision modification with reflected image |
US7931373B2 (en) | 2004-12-03 | 2011-04-26 | The Invention Science Fund I, Llc | Vision modification with reflected image |
US8562540B2 (en) | 2004-12-03 | 2013-10-22 | The Invention Science Fund I, Llc | Method and system for adaptive vision modification |
US8104892B2 (en) | 2004-12-03 | 2012-01-31 | The Invention Science Fund I, Llc | Vision modification with reflected image |
US8109632B2 (en) | 2004-12-03 | 2012-02-07 | The Invention Science Fund I, Llc | Vision modification with reflected image |
US8244342B2 (en) | 2004-12-03 | 2012-08-14 | The Invention Science Fund I, Llc | Method and system for adaptive vision modification |
US8282212B2 (en) | 2004-12-03 | 2012-10-09 | The Invention Science Fund I, Llc | Vision modification with reflected image |
EP1727072A1 (en) * | 2005-05-25 | 2006-11-29 | The Babraham Institute | Signal processing, transmission, data storage and representation |
US20090093709A1 (en) * | 2007-05-18 | 2009-04-09 | Beth Israel Deaconess Medical Center, Inc. | Noise reduction system and methods for magnetic resonance imaging |
US8947976B2 (en) * | 2008-10-03 | 2015-02-03 | Westerngeco L.L.C. | Harmonic attenuation using multiple sweep rates |
US20100085837A1 (en) * | 2008-10-03 | 2010-04-08 | Kambiz Iranpour | Harmonic attenuation using multiple sweep rates |
US20110150253A1 (en) * | 2009-12-22 | 2011-06-23 | Siemens Medical Instruments Pte. Ltd. | Method and device for setting a hearing device by detecting listening effort |
US8559645B2 (en) * | 2009-12-22 | 2013-10-15 | Siemens Medical Instruments Pte. Ltd. | Method and device for setting a hearing device by detecting listening effort |
US11051737B2 (en) * | 2017-05-19 | 2021-07-06 | Ricoh Company, Ltd. | Biomagnetic measurement method, biomagnetic measuring device, and biomagnetic measuring system |
JP2018192236A (en) * | 2017-05-19 | 2018-12-06 | 株式会社リコー | Biomagnetic measurement method, biomagnetic measuring device, and biomagnetic measuring system |
US20180333062A1 (en) * | 2017-05-19 | 2018-11-22 | Taishi WATANABE | Biomagnetic measurement method, biomagnetic measuring device, and biomagnetic measuring system |
US11723579B2 (en) | 2017-09-19 | 2023-08-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement |
US11717686B2 (en) | 2017-12-04 | 2023-08-08 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to facilitate learning and performance |
US11273283B2 (en) | 2017-12-31 | 2022-03-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
US11318277B2 (en) | 2017-12-31 | 2022-05-03 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
US11478603B2 (en) | 2017-12-31 | 2022-10-25 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
US11364361B2 (en) | 2018-04-20 | 2022-06-21 | Neuroenhancement Lab, LLC | System and method for inducing sleep by transplanting mental states |
US11452839B2 (en) | 2018-09-14 | 2022-09-27 | Neuroenhancement Lab, LLC | System and method of improving sleep |
US11786694B2 (en) | 2019-05-24 | 2023-10-17 | NeuroLight, Inc. | Device, method, and app for facilitating sleep |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5020538A (en) | Low noise magnetoencephalogram system and method | |
CA2149915C (en) | Method and apparatus for separating the different components of evoked response and spontaneous activity brain signals as well as of signals measured from the heart | |
US5436564A (en) | Bio-potential signal processor for magnetic resonance imaging | |
Vrba et al. | Signal processing in magnetoencephalography | |
Masterton et al. | Measurement and reduction of motion and ballistocardiogram artefacts from simultaneous EEG and fMRI recordings | |
JP4875696B2 (en) | Method and apparatus for suppressing interference in electromagnetic multi-channel measurements | |
US5269325A (en) | Analysis of biological signals using data from arrays of sensors | |
EP0982597B1 (en) | Method and apparatus for eliminating background interference signals from multichannel detector arrays | |
US10307105B2 (en) | Method and device for recognizing and removing undesired artifacts in multichannel magnetic field or electric potential measurements | |
JPH1052409A (en) | Device and method for wave filter to eliminate emi | |
Reite et al. | The human magnetoencephalogram: some EEG and related correlations | |
Taulu et al. | Novel noise reduction methods | |
JP6890484B2 (en) | Magnetic field measuring device and measured magnetic field display method | |
de Pinto | Filters for the reduction of baseline wander and muscle artifact in the ECG | |
EP0638869A1 (en) | Procedure for a high-resolution spectral analysis for multichannel observations | |
EP0443069A1 (en) | Method of measuring the electric or magnetic field pattern with a sensor arrangement | |
EA039153B1 (en) | Noise removal in magnetometer for medical use | |
parimita Swain et al. | A feasibility study to measure magnetocardiography (MCG) in unshielded environment using first order gradiometer | |
Reermann et al. | Real-time biomagnetic signal processing for uncooled magnetometers in cardiology | |
WO1997005823A1 (en) | Artifact removal from physiological signals | |
Vrba et al. | Noise cancellation by a whole-cortex SQUID MEG system | |
Fife et al. | Synthetic gradiometer systems for MEG | |
Sternickel et al. | Nonlinear noise reduction using reference data | |
McKay et al. | Implementation of a multi-channel biomagnetic measurement system using DSP technology | |
Kawakatsu | Application of ICA to MEG noise reduction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
REMI | Maintenance fee reminder mailed | ||
LAPS | Lapse for failure to pay maintenance fees | ||
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20030604 |