WO2002013689A2 - Method and apparatus for reducing contamination of an electrical signal - Google Patents
Method and apparatus for reducing contamination of an electrical signal Download PDFInfo
- Publication number
- WO2002013689A2 WO2002013689A2 PCT/US2001/025480 US0125480W WO0213689A2 WO 2002013689 A2 WO2002013689 A2 WO 2002013689A2 US 0125480 W US0125480 W US 0125480W WO 0213689 A2 WO0213689 A2 WO 0213689A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- signal
- eeg
- electrical
- contaminating
- electrophysiological
- Prior art date
Links
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/30—Input circuits therefor
-
- 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/316—Modalities, i.e. specific diagnostic methods
-
- 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
Definitions
- This invention relates generally to methods and apparatus for signal processing and data collection that are particularly suited for mi- ⁇ imizing artifacts and optimizing signal-to- noise in simultaneous recording of electroencephalographic (EEG) and Magnetic Resonance Imaging (MRI) signals, as well as other environments in which electrical signals are subject to repeated interference.
- EMG electroencephalographic
- MRI Magnetic Resonance Imaging
- the methods of the invention can be applied to other recordings containing repeated electrical interference, including eletromyelographic (EMG), electrocardiographic (ECG) or galvanic skin resistance (GSR) signals recording during fMR and audio recordings or transmissions in the presence of 60 Hz noise or electrical transients.
- EMG eletromyelographic
- ECG electrocardiographic
- GSR galvanic skin resistance
- Electroencephalography (EEG) and functional MRI (fMRI) induce mutual artifacts when recorded concurrently.
- Electroencephalography (EEG) has been a key tool for study of the brain for decades. However, despite its multiple clinical and research uses, such as in epilepsy (Ebersole, 1997), sleep staging (Rechtilles & Kales, 1968) and psychophysiology, little is yet known about the underlying generators of EEG activity in humans.
- Functional MRI (fMRI) recorded in concert with EEG can provide a method for localizing these sources.
- EEG electroencephalography
- fMRI functional MRI
- the invention provides a method of reducing contamination of electrical signals recorded in the presence of repeated interference contamination.
- the method comprises obtaining an electrical signal, wherein the electrical signal was recorded in the presence of a contaminating signal, and detecting a timing signal that occurs at a fixed time point during the electrical signal relative to the onset of the contaminating signal.
- the method further comprises digitizing the electrical signal, wherein the digitizing begins with the timing signal.
- a plurality of digitized electrical signals is then analyzed, wherein the electrical signals are synchronized with respect to the timing signal, to obtain an estimated contaminating signal.
- the estimated contaminating signal is subtracted from the digitized electrical signal, thereby reducing contamination of the electrical signal.
- the analysis to obtain the estimate of the contaminating signal comprises averaging the electrical signals.
- the analysis to obtain the estimate of the contaminating signal comprises calculating a weighted average of the electrical signals.
- the estimate of the contaminating signal can be biased towards recent events, for example, by adding the nth electrical signal to a scalar multiple, w, of the prior estimate of the contaminating signal and dividing this first sum by a second sum obtained by adding the series l+w 2 +w 3 +w 4 ... + ... w n .
- the estimate of the contaminating signal further can be multiplied by a scalar prior to the subtracting step.
- the method is particularly suitable for electrical recordings which comprise an electrophysiological signal, such as an electroencephalographic recording, an electromyelographic recording, an electrocardiographic recording or a measure of galvanic skin resistance.
- an electrophysiological signal such as an electroencephalographic recording, an electromyelographic recording, an electrocardiographic recording or a measure of galvanic skin resistance.
- the method is applicable as well to other types of electrical recordings, including audio recordings.
- the interference comprises interference arising from inductively coupled magnetic fields.
- the interference can also comprise interference arising from alternating current (AC) line noise.
- the digitizing can be performed at a sampling rate below the Nyquist rate for the contaminating signal.
- the electrical signal obtained is passed through a low pass filter prior to the digitizing, at a frequency of approximately one half of the frequency at which the electrical signal is sampled.
- the low pass filter may pass signal frequencies of less than about 200 Hz.
- the method can be performed concurrently with Magnetic Resonance Imaging of the subject.
- the electrical signal comprises an electrophysiological signal and the contaminating signal comprises gradient activity. Examples of a contaminating signal include radio frequency transmitter activity.
- the digitizing is performed at a rate of about 200 to about 5000 samples per second. The digitizing can be performed at rates below 200 and above 5000 samples per second, with representative rates including 100, 250, 500, 1000, 2000, 3000, 4000 and 6000 samples per second.
- the invention additionally provides a method of removing a DC offset from the electrical signal by analog subtraction prior to the digitizing.
- the DC offset is measured and subtracted from the electrical signal using a difference amplifier.
- the DC offset is measured by analog to digital conversion, and averaged over a time period long compared to the lowest frequencies of interest in the electrical signal. An example of such a long time period is approximately 10 times longer than the lowest frequencies of interest in the electrical signal. For example, where lowest frequencies of interest are approximately 3 Hz, the time period is about 30 seconds.
- the analog subtraction comprises converting the averaged signal to an analog voltage and electrically subtracting the averaged signal from the electrical signal through differential amplification.
- the DC offset is measured in an analog integrator having a time constant long compared with lowest frequencies of interest in the signal.
- the method of the invention is useful for electrophysiological recordings, such as in an electroencephalogram that is recorded concurrently with magnetic resonance image acquisition.
- the electrophysiological recordings are used to inform interpretations of magnetic resonance images.
- the electrophysiological recordings can be used in a statistical analysis of change in intensity of the magnetic resonance signal.
- the method can further comprise deterrriining a correlation between change in intensity of the magnetic resonance signal and a feature of the electrophysiological recording.
- the correlation can be used to make statistical images, or image maps, that represent an association between the electrical signals and the intensity of the magnetic resonance signal intensity.
- the feature of the electrophysiological recording comprises a time course of signal intensity change in defined frequency bands contained in the electrophysiological recording.
- the defined frequency bands can be selected to correspond to standard ranges used for clinical interpretations of the electroencephalogram. Representative standard ranges are selected from the group consisting of from 0 to approximately 4 Hz (the Delta band), from approximately 4 to approximately 8 Hz (the Theta band), from approximately 8 to approximately 12 Hz (the Alpha band), from approximately 12 to approximately 30 Hz (the Beta band), and from approximately 30 Hz and greater (the Gamma band). Typically, the frequency bands in this context will not extend beyond 300 Hz.
- the method further comprises convolving the time course of the electrophysiological signal with an estimate of the magnetic resonance hemodynamic impulse response function. In this embodiment, the time course of the electrophysiological signal is suitably conditioned to more accurately reflect the anticipated time course of the magnetic resonance signal change.
- the invention additionally provides a method of reducing magnetic interference during electrophysiological recording from a subject by measuring an electrical potential difference between a pair of electrodes, wherein the pair of electrodes communicate with a differential amplifier via electrical connections, the method comprising twisting the electrical connections together, thereby reducing magnetic interference.
- a method of reducing magnetic interference during electrophysiological recording from a subject by measuring an electrical potential difference between a pair of adjacent electrodes, wherein each electrode comprises two leads, the method comprising twisting each lead together with a lead of an adjacent electrode, thereby reducing magnetic interference.
- the electrophysiological recording comprises an electroencephalographic recording. The method can be performed concurrently with Magnetic Resonance Imaging of the subject.
- the invention further provides an apparatus for processing digitized electrical signals in the presence of a repeated contaminating signal.
- the apparatus comprises a signal processor adapted to receive a recording of an electrical signal; a detector adapted to detect a timing signal that occurs at a fixed time point during an electrical signal relative to the onset of a contaminating signal; a signal accumulator to contain the estimated contaminating signal; and a processor adapted to subtract averaged waveforms from an electrical signal.
- the signal accumulator can be, for example, a signal averager.
- FIG 1A Digital photograph of chained bipolar dual-lead dress. Leads of consecutive electrodes are twisted together to reduce scanner artifact.
- Figure IB Schematic diagram showing electrode connectors 2 attached to the head of a subject 1 using standard electrode gel. Each of the connectors is attached to 2 electrical wires (typically constructed of carbon fiber material, which reduces magnetic artifacts). The wires from adjacent electrodes are twisted tightly together in pairs 3, where each of the two wires from a single electrode is twisted together with a different neighbor. The electrode pairs are then presented to the input of a differential amplifier 4 where the electrical potential difference is amplified to form the electroencephalogram. The inputs to the amplifiers are bridged such that the paired leads enclose a complete loop, thereby minimizing additional differential potentials between amplifiers. The wires are drawn in thick and thin lines for clarity only.
- FIG. 2A Schematic diagram showing how dual lead electrodes allow each bipolar pair 3 to be twisted together for their entire length, sending signal directly to local differential amplifiers 4.
- Figure 2B Schematic diagram showing how twisting of leads leaves only small loops at the head in which e.m.f. can be induced. Current induced in lead twists by motion and gradient switching will be self canceling.
- FIG. 3 Diagram of an EEG data pathway. EEG signal is fed to a local differential amplifier 4, digitized and then sent out of the scanner room 7 via optical fiber for real time display and off-line analysis.
- FIG. 4 Traces illustrating the baUistocardiogram subtraction algorithm, shown using data collected on a normal volunteer.
- Figure 5. EEG of a phantom recorded inside the MR scanner using twisted and untwisted leads, showing recordings both in absence of scanning and during EPI. With or without scanning, EEG recorded using the untwisted leads was significantly noisier.
- FIG. 6 Raw EEG of a normal volunteer recorded inside the MR scanner using untwisted (top) and twisted (bottom) leads. EEG recorded with the untwisted leads was significantly noisier. BaUistocardiogram is visible in the twisted lead data, but is reduced compared to the untwisted leads.
- FIG. 7 EEG recorded on a volunteer during fMRI, before (above) and after (below) post-processing to remove the baUistocardiogram.
- Figure 8 Power spectrum of EEG recorded simultaneously with fMRI, in steps equal to the TR of 2.5 seconds, showing expected increases in the alpha band (8-12 Hz) when the subject's eyes were closed.
- Figure 9 Graph depicting how errors in gradient noise canceUation will occur when the sampling and gradient activity are asynchronous.
- the sampling used to create an error estimate (open circles) has drifted by approximately 200 ⁇ s compared to the current sample (closed circles). Subtracting the error estimate actually increases the residual error (by more than 17%), as indicated by diamonds.
- FIG. 11 Top: Graph representing uncorrected and corrected signals obtained during imaging, using 10 kHz sampling. Bottom: Enlarged (10X) views of the 25 msec periods indicated in dashed lines (uncorrected signal in dotted lines). Note that the artifact suppression varies from cycle to cycle, as a result of phase errors.
- Figure 12 Graph showing simulation of the effects of sampling rate on the efficiency of gradient artifact suppression.
- the trace at the top is actual EEG data recorded at 10 kHz during an echo-planar imaging sequence.
- the three traces below it are the difference between the original signal and a sample lagged at 100 ⁇ s, 200 ⁇ s and 400 ⁇ s as indicated (the worst case errors for 10 kHz, 5 kHz and 2500 Hz sampling).
- the graph below shows a detaU of the period during which echo-planar readout occurs.
- Figure 13 Graph similar to that shown in Figure 4, depicting a simulation that shows the effects of timing errors when sampling at 200 Hz - weU below the Nyquist frequency for the gradients. The waveforms are clearly undersampled and therefore appear quite different from the prior Figure. The magnitudes of the residual artifacts, however, are very similar. The three lower curves show the artifacts that remain after correction with timing errors of 100, 200 and 400 ⁇ s.
- Figure 14 Schematic representation of a low-cost offset nulling differential amplifier circuit for use in fMRI. Power supply connections are omitted for clarity. AU capacitor values are in microfarads.
- Figure 15 MR gradient activity recorded with triggered 200 Hz sampling using a 3s repetition time (TR) and nineteen slices.
- the uncorrected signal is shown at top for a single TR. In the middle is shown the average of 30 TR periods, and at bottom, the difference between the uncorrected and averaged signals ("corrected").
- FIG. 1 Human EEG data coUected during echo-planar functional imaging. The uncorrected data appear at bottom. Above them are corrected records from twenty successive TR periods.
- Figure 17A Graph showing energy as a function of frequency and time, derived from EEG data acquired during scanning.
- FIG 17B Graph showing estimated fMRI activation time course for the EEG data appearing in Figure 9A. For clarity, only the Alpha (solid line) and Theta (dashed line) bands are shown.
- Figure 18A-C Functional MRI statistical maps of signal change correlated with spectral energy at each of five frequency bands (18A left, Delta; 18A right, Theta; 18B left, Alpha; 18B right, Beta; 18C, Gamma), expressed as coefficient of correlation. Note that the color scale for the lower frequencies (18A) is different, as the correlations were overaU lower. AU five images were calculated from the same 4:30 (min:sec) acquisition taken with the subject at rest with eyes open.
- Figure 19 Functional block diagram of MRI-compatible EEG amplifier.
- Figure 20 Alternative differential amplifier circuit for use with fMRI.
- Figure 21 Flow chart illustrating method of signal correction by artifact reduction.
- FIG. 22 Schematic illustration of elements of an apparatus of the invention.
- the invention is based on the discovery that contamination of a digitally encoded electrical signal can be reduced significantly by making use of a timing signal that is associated with the onset of a repeated contamination signal.
- a timing signal can be used to aUgn the digitization of repeated contamination signals for deterrrjining an estimate of the contamination which can then be subtracted from the electrical signal.
- This method is particularly suited for use with electrophysiological signals, such as EEG, ECG, EMG and galvanic skin response (GSR), and for elimination of noise associated with concurrently used methods such as MRI.
- timing signals synchronized with respect to a timing signal means that data corresponding to electrical signals recorded over a period of time are aligned so that a timing signal that occurs within each of the electrical signals recorded over time is superimposed by the alignment.
- the timing signal as weU as any other signal that recurs at a fixed time point relative to the timing signal, will be enhanced relative to non-recurring signals.
- twisted means united by winding, tertwining or coiling. A pair of electrical connections between a pair of electrodes and a differential amplifier is sufficiently twisted if the enclosed magnetic fields are substantiaUy reduced.
- the invention provides a method of reducing contamination of electrical signals recorded in the presence of repeated interference contamination.
- the method comprises obtaining an electrical signal, wherein the electrical signal was recorded in the presence of a contaminating signal, and detecting a timing signal that occurs at a fixed time point during the electrical signal relative to the onset of the contaminating signal.
- the method further comprises digitizing the electrical signal, wherein the digitizing begins with the timing signal.
- a plurakty of digitized electrical signals is then analyzed, wherein the electrical signals are synchronized with respect to the timing signal, to obtain an estimated contaminating signal.
- the estimated contaminating signal is subtracted from the digitized electrical signal, thereby reducing contamination of the electrical signal.
- the analysis to obtain the estimate of the contaminating signal comprises averaging the electrical signals. In some embodiments, the analysis to obtain the estimate of the contaminating signal comprises calculating a weighted average of the electrical signals.
- the use of weighted averages can serve to achieve an adaptive artifact reduction.
- the estimate of the contaminating signal can be biased towards recent events, for example, by adding the nth electrical signal to a scalar multiple, w, of the prior estimate of the contaminating signal and dividing this by the sum of the series l+w 2 +w 3 +w 4 ... + ... w n .
- the estimate of the contaminating signal can be multipked by a scalar prior to the subtracting step.
- Figure 19 is a general functional block diagram of representative analog electronics useful in the method of the invention.
- the method performs optimaUy with an adequately linear electronic signal.
- signal is applied differentiaUy at the input terminals 42.
- contaminating signal from sources such as radio frequency is attenuated before being differentiaUy amplified using standard components 46.
- the differential amplifier 46 is commonly provided with an offset reference input, such that voltages appearing at this terminal are subtracted from the output.
- the amplified DC offset potential, derived from the inputs is first measured, then applied to the differential input amplifier 46.
- an active low pass filter 48 is supplied to the differential input amplifier 46.
- the filtered signal is buffered in the output amplifier 50 before being made available to the digitizer circuit.
- a switching means 54 is provided to the sample and hold circuit to aUow detection of the input DC offset at any time desired.
- the circuit modeled in Figure 19 is an example that is especiaUy suited to the problem of recording electrophysiological signals during Magnetic Resonance Imaging.
- the functional logic of this diagram is shown in somewhat more detail in Figure 20 which, in addition, illustrates a means of transmitting the EEG data out of the MRI suite.
- Figure 20 which, in addition, illustrates a means of transmitting the EEG data out of the MRI suite.
- One skilled in the art wiU see immediately that many different circuit topologies are possible that will accomplish essentiaUy the same function for this or other appHcations.
- the differential ampUfier 46 may be desirable to substitute a bandpass filter for 48.
- the passive filter 44 may be eliminated when any signals that reach the input to the input differential amplifier 46 are of suitably low amplitude to avoid saturation effects.
- EEG signal from the subject 1 is carried via twisted pair leads 3 to the input of a battery powered head amplifier.
- the inputs include RF attenuation via a series inductance and a paraUel capacitor. These inputs are mixed in the inputs to shield driver 64, whose output is applied to a concentric shield 56 surrounding the twisted leads 3 and connected to the subject 1.
- the inputs are coupled to a differential amplifier 58, whose output is applied to an isolation ampUfier, such as the IS0122 from Burr-Brown Corporation.
- This device, and a second simUar device 62 provide electrical isolation to the subject and an added safety factor for EEG recording.
- isolation amp 60 The output from isolation amp 60 is sampled by an analog to digital convertor 66, whose digital output is stored in a latch 68 and converted to an analog voltage by digital to analog convertor 70. This output is applied as a DC correction 63 to differential ampHfier 60 after electrical isolation by isolation ampUfier 62.
- the output from this head ampUfier is presented to a low-pass filter to attenuate signal outside of the desired range of the EEG signal.
- This single channel output may be multiplexed with the outputs of other similar ampUfiers by analog multiplexor 74 clocked by hardware timer 78.
- the output of the multiplexor 74 may be converted to an optical signal by optocoupler 76 and transmitted by .optical fiber, together with the clock signal to a second optocoupler 82 that is located outside of the MR shielded room through a penetration panel 80.
- De-multiplexor 82 is used to separate the signals from multiple ampUfiers and its output may be presented to a differential line driver 88 for transmission of the signals over long distances to an analog to digital convertor for later processing.
- Figure 21 is a flow chart illustrating an appUcation of the method of the invention.
- the raw digitized signal, containing both the desired signal and contaminating artifact are shown as (1) in Figure 21. It is strongly preferred, for the invention to work optimally, that the signal in (1) be faithfully (linearly) recorded and that the digitization be timed with adequate precision to the artifact.
- an estimate of the artifact (2) is simply subtracted from the taw signal.
- the estimate of the artifact may be multipUed by an ampUtude constant, m, to account for differences in the magnitude of the coupling of the artifact to the signal recording system.
- the artifact estimate can be computed as foUows: Each time a new raw sample (1) is avaUable, it is added to the current estimate of the artifact (2) which has been multipUed by a scalar ampUtude constant, w (4). Tliis summed signal is then divided by the sum of the series 1 + w + w 2 + w 3 + . .. resulting in a new representation of the artifact (5), which then replaces the value used for (2). When w is a number less than 1, the process results in a leaky average, where more recent signals have a greater influence on the estimated artifact (2) than do less recent signals. In this way, the system adapts to slow changes in the artifact, if necessary. For the purposes of this disclosure, we have caUed this implementation a "leaky averager.”
- the method is particularly suitable for electrical recordings which comprise an electrophysiological signal, such as an electroencephalographic recording, an electromyelographic recording, an electrocardiographic recording or a measure of galvanic skin resistance.
- the method is appHcable as weU to other types of electrical recordings, including audio recordings.
- the interference comprises interference arising from inductively coupled magnetic fields.
- the interference can also comprise interference arising from other sources, such as alternating current (AC) line noise.
- the digitizing can be performed at a sampling rate below the Nyquist rate for the contaminating signal.
- the electrical signal obtained is passed through a low pass filter prior to the digitizing, at a frequency of approximately one half of the frequency at which the electrical signal is sampled.
- the low pass filter may pass signal frequencies of less than about 200 Hz.
- the method can be performed concurrently with Magnetic Resonance Imaging of the subject.
- the electrical signal comprises an electrophysiological signal and the contaminating signal comprises gradient activity. Examples of a contaminating signal include radio frequency transmitter activity.
- the digitizing is performed at a rate of about 200 to about 5000 samples per second.
- the digitizing can be performed at rates below 200 and above 5000 samples per second, with representative rates including 100, 250, 500, 1000, 2000, 3000, 4000 and 6000 samples per second.
- the invention additionaUy provides a method of removing a DC offset from the electrical signal by analog subtraction prior to the digitizing.
- the DC offset is measured and subtracted from the electrical signal using a difference ampUfier.
- the DC offset is measured by analog to digital conversion, and averaged over a time period long compared to the lowest frequencies of interest in the electrical signal.
- An example of such a long time period is approximately 10 times longer than the lowest frequencies of interest in the electrical signal. For example, where lowest frequencies of interest are approximately 3 Hz, the time period is about 30 seconds.
- the analog subtraction comprises converting the averaged signal to an analog voltage and electricaUy subtracting the averaged signal from the electrical signal through differential ampUfication.
- the DC offset is measured in an analog integrator having a time constant long compared with lowest frequencies of interest in the signal.
- the method of the invention is useful for electrophysiological recordings, such as an electroencephalogram that is recorded concurrently with magnetic resonance image acquisition.
- the electrophysiological recordings are used to inform interpretations of magnetic resonance images.
- the electrophysiological recordings can be used in a statistical analysis of change in intensity of the magnetic resonance signal.
- the method can further comprise deterrriining a correlation between change in intensity of the magnetic resonance signal and a feature of the electrophysiological recording.
- the correlation can be used to make statistical images, or image maps, that represent an association between the electrical signals and the intensity of the magnetic resonance signal intensity.
- the feature of the electrophysiological recording comprises a time course of signal intensity change in defined frequency bands contained in the electrophysiological recording.
- the defined frequency bands can be selected to correspond to standard ranges used for clinical interpretations of the electroencephalogram. Representative standard ranges are selected from the group consisting of from 0 to approximately 4 Hz (the Delta band), from approximately 4 to approximately 8 Hz (the Theta band), from approximately 8 to approximately 12 Hz (the Alpha band), from approximately 12 to approximately 30 Hz (the Beta band), and from approximately 30 Hz and greater (the Gamma band). TypicaUy, the frequency bands for this appUcation will not extend beyond 300 Hz.
- the method further comprises convolving the time course of the electrophysiological signal with an estimate of the magnetic resonance hemodynamic impulse response function. In this embodiment, the time course of the electrophysiological signal is suitably conditioned to more accurately reflect the anticipated time course of the magnetic resonance signal change.
- the invention additionaUy provides a method of reducing magnetic interference during electrophysiological recording from a subject by measuring an electrical potential difference between a pair of electrodes, wherein the pair of electrodes communicate with a differential ampUfier via electrical connections, the method comprising twisting the electrical connections together, thereby reducing magnetic interference. Also provided is a method of reducing magnetic interference during electrophysiological recording from a subject by measuring an electrical potential difference between a pair of adjacent electrodes, wherein each electrode comprises two leads, the method comprising twisting each lead together with a lead of an adjacent electrode, thereby reducing magnetic interference.
- the electrophysiological recording comprises an electroencephalographic recording. The method can be performed concurrently with Magnetic Resonance Imaging of the subject.
- the invention further provides an apparatus for processing digitized electrical signals in the presence of a repeated contaminating signal.
- the apparatus comprises a signal processor 90 adapted to receive a recording of an electrical signal; a detector 92 adapted to detect a timing signal that occurs at a fixed time point during an electrical signal relative to the onset of a contaminating signal; a signal accumulator 94 to contain the estimated contaminating signal; and a processor 96 adapted to subtract averaged waveforms from an electrical signal. Representative variations of the apparatus are described in Figures 3, 14, 19, and 20.
- the signal processor 90 adapted to receive a recording of an electrical signal may be, for example, an electronic circuit (IC) consisting of a integrated differential ampUfier such as a an INA114 from Burr Brown Corp., with a DC offset reference provided through measurement by a sample and hold IC such as the LF298 from National Semiconductor corporation, a low pass active filter and an output buffer aU made using standard operational ampUfier ICs.
- the detector can also include a means of analog to digital conversion, such as a National Instruments NI 603 IE instaUed in a personal computer.
- the detector 92 adapted to detect a timing signal that occurs at a fixed time point during an electrical signal relative to the onset of a contaminating signal may be, for example, an optoisolator IC whose output is conditioned using an IC such an LN555 from National Semiconductor corporation to produces a TTL compatible trigger signal which is then presented to an analog to digital convertor such as a National Instruments NI 603 IE instaUed in a personal computer.
- the signal accumulator 94 may be implemented, for example, in software in the 'C programming language as a vector of numbers, or in the National Instruments LabView prograrriming language as an array of numbers. Either of these may be executed on a personal computer.
- the signal accumulator 94 can be, for example, a signal averager. Other methods, in addition to signal averaging, can be used to generate an estimate of the contaminating signal.
- the processor 96 adapted to subtract averaged waveforms from an electrical signal may be implemented, for example, as a processing routine in the 'C prograx-ming language or the National Instruments LabView programming language running on a personal computer. Those sk ⁇ led in the art will appreciate variations on the above examples of apparatus elements that will serve the same processing, detecting and accumulating functions in accordance with the methods of the invention. Overview of EEG and fMRI in Simultaneous Recording and Mapping
- Electroencephalography is estabUshed firmly as a means to probe changes in electrical signals recorded from the scalp that accompany behavioral tasks, and as a marker for clinical, cognitive or neural states. Determination of the tlxtee-dimensional localization of the EEG signal is ambiguous because the relationship between the actual position of multiple electrical dipoles and the distribution of electrical potentials detected at the scalp has no unique solution.
- Functional Magnetic Resonance Imaging fMRI uses modulations in the magnetic resonance signal that depend on variations in blood oxygenation to distinguish brain regions whose activity is increased or decreased with task demands.
- the foUowing describes a set of solutions to the technical problems in simultaneous recording of fMRI and EEG and shows that data from the two methods may be combined to create tomographic images indicating brain regions whose activity changes as a function of EEG signal intensity in the classicaUy defined spectral bands.
- EEG is now a routine and essential test in clinical neurology. It provides diagnostic information that cannot be gathered through any other commonly obtainable means. Indeed, its ⁇ dispensability derives from this lack of other routine clinical tools to assess a broad region of cerebral neurophysiology with high temporal resolution. EEG depicts moment-to-moment changes in cerebral cortical function, and thus is valuable in any clinical context where such information would guide medical decision making. Such situations are not limited to electrophysiologic abnormaUties, as other pathologic processes often affect neuronal function and, thus, impact the EEG. These include ischemia, metaboUc alterations, mass effect, and infection among others (Markand 1984).
- brain regions active in non-REM sleep include the anterior hypothalamus, dorsal bulbar reticular formation, and nucleus of the soUtary tract (Jones 2000). Regions likely associated with REM activity, and possibly the generation of wakefulness, include the posterior hypothalamus, ventral mesencephaUc pons, basal forebrain, and pontine reticular formation. Such localized brain activity should also be visible in humans through functional neuroimaging. Both PET and SPECT have been used to examine regional activity changes (via cerebral metaboUsm and blood flow respectively) with sleep stage (determined electroencephalographicaUy). These imaging studies have broadly indicated that areas thought to be involved in the active generation of rapid eye movement (REM) sleep are active during this stage of sleep and in non-rapid eye movement (NREM) sleep.
- REM rapid eye movement
- NREM non-rapid eye movement
- the sleep electroencephalogram can be defined by characteristic frequency patterns, and brief electrophysiological phenomena such as k-complexes and sleep spindles. Furthermore, during REM sleep, there is a descending suppression of muscle tone, saccadic eye movements, and the loss of the cortical synchrony that is a hallmark of NREM sleep. Only with temporal resolving power of fMRI wiU it be possible to study activations associated with these transient events.
- One appUcation of the invention is to assess changes in regional brain activity using fMRI and to correlate such activity to classicaUy defined sleep architecture and features. For example, the invention can be used to seek an understanding the brain activity that underlies the general lack of responsiveness to external stimuli, the apparent gating of motor output and dream states.
- fMRI Functional MRI
- the classical or posterior alpha rhythm is found mostly in occipital, parietal, and posterior temporal regions (Adrian and Matthews 1934), and first emerges at about 4 months of age as a 4 Hz osculation, present with eyes closed and blocked with eyes open.
- the frequency of this rhythm increases with age, reaching about 8 Hz by age 3, and by about age 10, it reaches the average adult frequency of 10 Hz (Petersen and Eeg-
- Lopes da S va demonstrated significant thalamocortical coherences in dogs between lateral geniculate nucleus and pulvinar and the cortex (Lopes da SUva, Lierop et al. 1973; Lopes da Silva, Vos et al. 1980).
- Lopes da Silva showed an inverse correlation between EEG alpha power and thalamic metaboUc rate in normal subjects using PET (Lindgren, Larson et al. 1999).
- EEG cardio- electric
- EMG myoelectric
- the noise environment for EEG becomes radicaUy worse when subjects are placed inside an MR imaging system. Almost aU of the noise sources, however, are coupled magneticaUy to the EEG.
- Several of these are non-biological, including: ampUfier noise in the shim and field gradient ampUfiers; large time-varying magnetic fields induced by the field gradients during scanning, and radio-frequency signals generated by the scanner for magnetic resonance induction.
- ampUfier noise in the shim and field gradient ampUfiers large time-varying magnetic fields induced by the field gradients during scanning, and radio-frequency signals generated by the scanner for magnetic resonance induction.
- the magnitude of the voltage induced by these time varying fields is proportional to the first time derivative of the flux, and thus to the ampUtude of the magnetic field, its first time derivative, and the area enclosed by any conducting loop. More specificaUy:
- e.m.f. is the induced electromotive force
- ⁇ is the magnetic flux.
- the MR imaging gradients in state-of-the-art imaging instruments slew at an extremely high rate; the field gradients on a typical scanner operate at 80 T/sec and are thus major sources of noise; the newest generation of MR instruments, with local head gradients, will slew two to three times faster.
- Physiological noise sources are present also. Even smaU motions of the subject are coupled to the EEG as the leads move within the large static magnetic field. Not only minor fidgeting, but also the motion of the whole body with each heartbeat (the baUistocardiogram), produces signal in the ⁇ V to mV range.
- the scalp electrical potentials used in the EEG contain both time varying and static (DC) components. Often the DC signal is much larger than the EEG, but it is seldom of interest for clinical diagnostic purposes (note reference to true DC offsets, not to slowly varying potentials), as it contains essentiaUy no information. However, it does cause trouble for the EEG in several ways. TypicaUy, the DC offset increases the dynamic range needed to digitize the EEG signal. For example, the EEG signal may be only a few ⁇ V, while potentials of a few mV may exist between electrodes, or as a result of the chemical electrode potential. The signal digitization depth will be reduced by the ratio of the EEG to the DC potential.
- the invention further provides a more complete solution that avoids the high pass filter completely.
- the magneticaUy-induced gradient noise is of very large ampUtude (rriilUVolts in a typical scanner) as compared to the EEG, especiaUy in the context of echo-planar imaging.
- One group has implemented a correction scheme for gradient artifacts that is similar to a scheme that group developed for baUistocardiogram removal (AUen, Josephs et al. 2000). Because the fundamental frequencies of the gradient activity are much higher, they developed special recording hardware that aUowed them to use a much higher digital sampling rate of 5 kHz, which they selected as being rapid compared to the nominal Nyquist frequency for the gradient waveshapes. Unfortunately this is not sufficient, as sampling at the Nyquist rate guarantees only against aHasing of the higher frequencies into the pass band, but does not effectively remove the artifact.
- the sampling has drifted with respect to the gradients by 200 ⁇ s (dashed line, open circles).
- the difference signal is the residual artifact, in this case just over 17% greater than the uncorrected signal.
- the AUen group has adopted a sophisticated interpolation scheme that is successful in minimizing the residual contamination.
- the present invention provides a much more effective approach based on an alternative formulation of the digital sampling problem.
- the methods disclosed herein have made it possible to effectively eliminate contamination of the EEG signal by the most severe sources of noise present during MRI scanning in general, and in functional MRI in particular.
- the method has been used in the construction of tomographic maps of brain activity corresponding to the energies in spectraUy-defined components of the EEG.
- Figure 10 shows the timing of a typical echo-planar imaging sequence, as used in typical functional studies.
- the lines for Select, Phase and Readout indicate the ampHtudes of the three orthogonal magnetic field gradients used for imaging.
- the fourth line indicates the timing of the radio frequency channel (only the ampUtude envelope is shown for the RF, as the carrier frequency of 128 MHz is not visible at this resolution.)
- Immediately apparent is the very large high frequency (1400 Hz) osculation of the readout gradient (shown at half the vertical magnification of the other gradients).
- Figure 11 shows raw signal, recorded from an EEG system at a sampling rate of 10 kHz, foUowing analog low-pass filtering at 100 Hz.
- the insets on the bottom of this figure show expanded representations (10X) of the indicated regions of the signals.
- the MRI field gradients induce voltages much higher than the cortical signal. Comparing this and the previous figure reveals that the magnitude of the high frequency components is reduced dramaticaUy by the low pass filter, and that the residual artifacts outlast the gradients themselves.
- the low pass filter provides at least 100-fold reduction in the 1400 Hz oscillations, it does not remove the transients as these gradients turn on and off. These contain energy at very low frequencies, as weU.
- the AC coupled input stage in this traditional ampHfier, is responsible for the extended ring-down of the artifacts (the saturation recovery aUuded to previously), although it is much better than the 0.1 to 1 s ring down reported by others (Lovblad, Thomas et al. 1999), presumably due in large part to the attenuation provided by the differential recording apparatus, which helps to prevent the ampUfiers from going into saturation.
- the maximum phase shift, ⁇ , that can occur at a given sampling rate is equal to 2 ⁇ f 0 /f s , where f s is the sampling frequency, and fo is the frequency of the EPI readout. Comparing the residual artifact during the two scan periods (two expanded frames at bottom of Figure 11) reveals that the canceUation efficiency is unstable as a consequence of the asynchronous timing of the gradient activity and sampling device, which causes the sampling offset, ⁇ , to drift over time. As the sampling rate is increased, the canceUation will become more accurate. Using the approximation that for small ⁇ , sinoc ⁇ and cosoc ⁇ 1, one can see that in this reg ne ⁇ is approximately proportional to ⁇ .
- the signal must be digitized at approximately 100*2* ⁇ times the highest frequency of interest, in this case (the 1400 Hz readout) about 880 kHz/channel, which is impractical for reasons of both cost and overaU data handling.
- the acceptable sampling error is predicted readily by this formula given the low pass filter characteristics and the desired final signal to noise.
- Figure 12 shows the efficiency of artifact subtraction based simply on rapid sampling. On the top is shown the raw artifact. Below it is shown the residual artifact that remains after subtraction if the timing of the sampling and the scanner have drifted from synchrony by 100, 200 and 400 ⁇ sec, corresponding to the worst case errors for sampling at 10 kHz, 5 kHz and the approximate Nyquist rate of 2500 Hz, respectively. It is immediately apparent that the residual artifact, after subtraction, is large even with the smaUest timing offset. In the graph at the bottom of Figure 12, which shows in greater deta only the echo-planar readout segment of that data, one can see that the simple subtraction actuaUy increases the magnitude of the artifact, as predicted in the equation above (and in Figure 9).
- the scanner artifact is about ten times the ampUtude of the EEG. Assuming an EEG signal to noise ratio of 100:1, a thousand-fold suppression is needed in the digital processing, achievable with either 8.8 MHz/channel sampling
- Example 1 Acquiring simultaneous EEG and functional MRI
- the EEG device incorporates numerous hardware modifications to reduce artifact in concurrent EEG/ MRI, and was provided by Telefactor Corporation (W. Conshohocken, PA). Signal is detected from the scalp using sUver chloride plated plastic cup electrodes connected to a compact magnet-compatible local ampUfier (headbox) via 10 foot carbon fiber leads with a resistance of 1 k ⁇ /foot. This design minimized both artifact in the MR images and the induction of RF current loops in the lead wires.
- a lead configuration was devised that minimized unwanted current induction by recording EEG in a hard-wired montage using special dual lead electrodes.
- the lead wires from each electrode pair were twisted together over their entire length, forming a chained bipolar montage for each hemisphere (Jp2- ⁇ , ⁇ -t4, t4-t6, t6-o2, o2-p4,p4-c4, c4-f4, f4-jp2;jp1-f7,f7-t3, t3-t5, t5-o1 , o1p3,p3-c3, c3- ⁇ , ⁇ -jp1).
- the dual leads aUowed each differential pair to be twisted together ( Figure 1A-B and Figure 2A).
- the magnet-compatible headbox contains 32 separate channel inputs, each with a differential ampUfier coupled to an RC filter having a time constant of 0.25 msec.
- EKG electrocardiogram
- a scan trigger channel was also used to receive a pulse from the scanner every TR to aid in postprocessing of the data.
- the signal was filtered in aU channels with a band pass of 0.5-70 Hz to further attenuate high frequency noise.
- a single lead was connected at c to headbox ground as an added patient safety measure, but was not used as a montage reference.
- the signal was fed to a battery powered isocoder containing an A/D converter where it was sampled at 200 Hz., and the digitized signal was carried out of the shielded magnet room via optical fiber to maintain the scanner's electromagnetic isolation.
- the data were routed to a Telefactor Digital EEG (D/EEG) (486 computer).
- D/EEG Telefactor Digital EEG
- the EEG data could then be viewed in real time and sent off-line via a lObaseT Ethernet connection to a post-processing and viewing station for further artifact attenuation.
- Scout scans of the entire brain were first acquired to localize sHce planes paraUel to the AC-PC line through the occipital cortex. To acquire EEG during functional scanning, the scan protocol was then specified to allow windows of readable EEG between gradient bursts.
- Post-processing and viewing was performed on a DeU Inspiron 3000 Pentium PC.
- the EEG data was viewed with Telefactor Twin software, and processed further using home-built software described below to remove remaining artifacts.
- the residual artifact included both noise from e.m.f. induced by the magnetic field gradients (which appeared in the EEG when a sHce is acquired) and baUistocardiogram. The latter occurred in a fairly regular pattern just lagging the Q wave of the EKG, but its morphology and ampUtude differed in each EEG channel.
- the gradient noise was first removed by blanking the EEG and EKG data for the duration of each MR sHce acquisition. FoUowing a trigger every TR, the 90 ms data segments containing scan artifact were replaced with zeros. Therefore, large deflections caused by the MR gradients did not corrupt the EEG data during the averaging and subtraction in further artifact removal.
- FIG. 4B illustrates the averaging and subtraction gorithm performed on the data in each EEG channel. Every trigger-to-trigger section of raw EEG data (A n ) was averaged with aU preceding sections. Because EEG and EKG should be uncorrelated, this method averaged out the EEG signal and left only baUistocardiogram (B n ). Data sections were weighted inversely with their temporal displacement from the current sample to compensate for slow changes in the baUistocardiographic artifact, calculating the baUistocardiogram in each trigger to trigger section using the weighted average
- the EEG data were analyzed to quantify noise reduction due to the lead dress, calculating the loss, in dB, using 30 second data segments as
- the subject was given verbal cues during the scan to keep his eyes open for the first two minutes, closed the next two minutes, then open for the last minute. For the next two tasks, the subject was told to keep his eyes closed.
- the subject was instructed to count backwards by sevens from the four digit number given to him two minutes into the scan. To verify that the task was performed, the subject was asked for the number he ended on after the three minutes of counting.
- the subject was instructed to visualize eating bis favorite meal, again two minutes into the scan.
- the alpha power in each TR was calculated using software developed in house. With this software, the EEG power in user-defined bands with each TR was found using the Fast Fourier Transform. Spectral power in the alpha band was used as a reference function to calculate fMRI signal maps. Results
- FIG. 5 shows the recorded EEG in both twisted and untwisted lead channels.
- the twisted leads reduced random noise power by an average of 5.4 dB.
- gradients caused large artifacts in both the twisted and untwisted lead sets, but this noise power was reduced by an average of 6.3 dB in the twisted leads.
- Figure 6 shows EEG data recorded on a normal volunteer inside the scanner when no scanning was taking place. Twisting the leads (shown here on the left hemisphere vs. the untwisted on the right hemisphere) reduced noise by an average of 7.5 dB across aU channels.
- Figure 7 shows EEG recorded on a volunteer during functional MRI before and after post-processing.
- the post-processing removed significant gradient and RF artifact, as weU as baUistocardiogram.
- this method aUows one to record clean EEG during functional MR scanning. While the EEG is obscured during gradient bursts, it recovers quickly after each sHce acquisition. With this method, then, a trade off must be made between brain coverage in functional scanning and the fraction of useable EEG in the data record. Consideration must also be given to the timing of sHce acquisition so that it does not overlap with desired spectral frequencies.
- the invention provides a potentiaUy powerful tool for localizing the sources of various EEG waveforms. Because the EEG is acquired simultaneously with fMRI instead of serially, it can be used as a direct source for the fMRI reference function. In this way, activation maps could be made of any relevant changes in the EEG: spike and slow wave patterns in epUepsy, spectral changes, or even event related potentials.
- Example 2 Method for removal of artifacts in simultaneous EEG and fMRI
- baUistocardiogram contains substantial energy in the range of most interest to clinical EEG (from 1 to 50 Hz), the other main sources — gradient noise and RF transmit noise — have much higher fundamental frequencies.
- the overwhelming majority of the gradient-related noise is at a fixed frequency of 1400 Hz, with significant energy distributions down to 100 Hz or less.
- the radio frequency energy is at radicaUy higher frequencies weU outside of the interesting pass band for EEG.
- An initial gain stage features a single pole filter to attenuate the large RF signals prior to sampling and provides enough gain to bring the EEG signal into the mV range without the artifacts causing saturation.
- the next stage offers 30 dB/octave attenuation at a corner frequency of 200 Hz (so that gamma range EEG is readily passed.)
- the AC coupling problem is handled in the final ampHfication stage, which is arranged to include a resettable offset-nulling circuit that stores any DC offset across a low-leakage capacitor.
- the offset nulling switch can be, e.g., a mechanical switch, or, so that the nulling can be performed under digital control as needed, presumably when the software detects that the signal is close to digital saturation, a CMOS switch.
- the AC coupling problem is handled in the final ampHfication stage, which is arranged to include a resettable offset-nulling circuit that stores any DC offset across a low-leakage capacitor.
- the offset nulling switch can be, e.g., a mechanical switch, or, so that the nulling can be performed under digital control as needed, presumably when the software detects that the signal is close to digital saturation, a CMOS switch.
- This design has been reduced to printed circuitry and a single channel of this system has been assembled. One can reduce this design to use with a personal computer and assemble a single channel of this system. Under lab test conditions, the circuit is able to hold the DC offset to within 1% for 10 minutes at a time, easUy meeting the system requirements for MR scanning.
- activation mapping data from two subjects were without neurological or radiological abnormaHties as assessed by a brief neurological inventory based on a form developed by the National Institutes of Health NIH and a neurological inventory performed by a board-certified neurologist. Beyond lying in the magnet during scanning with eyes closed, the subjects performed no expHcit cognitive task.
- AU scanning was performed on a General Electric (Waukesha, Wl) 3.0 Tesla Signa® scanner modified by Advanced NMR Systems (Wilmington, MA) for high performance echo-planar imaging (Brady, Cohen et al. 1991; Cohen, KeUey et al. 1996).
- mapping data after image-based shimming (Reese, Davis et al. 1995) and coUection of a scout scan, imaging was performed using a gradient echo EPI scan, as described above. These data, however, were acquired with a longer TR of 4 seconds and only four sUce planes.
- the pulse sequence was modified to include a 5 ⁇ s trigger pulse at the beginning of each TR period, the leading edge of which was used for synchronization of the EEG sampling acquisitions. Digitization of the EEG
- EEG data are acquired using a PCI-1200 (National Instruments, Houston, TX) on a pc- compatible microcomputer.
- LabView National Instruments, Houston, TX
- sampling software was developed that responds to the leading edge of the scanner trigger by acquiring a fixed number of samples at a user-specified rate. SpecificaUy, with a TR of 3 seconds, one can acquire 599 samples at a rate of 200 Hz foUowing each trigger. These data are immediately flushed to a file in a background process; to accommodate this process one sample point is dropped with each TR.
- a frame of EEG data is coUected from each channel with each scan TR, and the frames are averaged separately for each channel to create an accurate representation of the scanner artifact.
- the next step is to manuaUy inspect the EEG data for characteristic artifacts, such as eye blink and facial muscle movement that are recognized easUy by their morphology within and across channels.
- baUistocardiogram suppression procedures have been pubHshed previously by Goldman et al., (Goldman, Stern et al. 2000) and wiU only be summarized here briefly. Because the cardiac-induced motion is nearly repeated with each heartbeat, the resulting artifact is essentiaUy the same and is superimposed onto the desired EEG. To remove this artifact one can detect the heartbeat using the electrocardiographic signal and calculate the average baUistocardiogram over many heartbeats. One can then subtract this average from the EEG signal. Of course, the morphology of the baUistocardiogram differs in each of the EEG channels as the motion of each lead differs sHghtly.
- this averaging and subtraction process is preferably performed separately for each of the EEG lead pairs.
- the method is conceptuaUy like that described by AUen et al., (AUen, Polizzi et al. 1998) but differs in the details by which slow variations in the baUistocardiogram are accommodated.
- scanSTAT available via the internet at the URL for the UCLA Brain Mapping Center, http: / /www.brainmapping.org
- scanSTAT available via the internet at the URL for the UCLA Brain Mapping Center, http: / /www.brainmapping.org
- the convention is to indicate regions of increasing positive correlation in colors from red to yeUow and areas of increasing negative (anti) correlation in colors from blue to cyan.
- regions of high negative correlation are interpreted as indicating decreased blood flow and metaboHc activity.
- Tested was the principle method of triggered sampling by recording the analog waveform that drives the gradient coUs, to determine the efficacy of the method independent of physiological signal fluctuations.
- the uncorrected signal from a single channel appears in Figure 15 (top).
- the averaging process removes noise uncorrelated with the imaging gradients. FinaUy, this averaged signal was subtracted from the uncorrected signal, yielding the corrected signal appearing in the lower portion of Figure 15.
- the large ampUtude gradient activity is removed completely and aU that remains in the signal is the smaU, uncorrelated, noise fluctuation.
- the echo-planar imaging gradient artifact which lasts for 23 ms and goes through 32 sinusoidal osculations in that period, is sampled with only four data points for these data.
- the digitized waveform is thus a very crude record of the actual gradient activity — weU below the Nyquist frequency.
- the morphology of the gradient waveform is grossly different for each of the nineteen sUces that are acquired in this TR period. This is due to the fact that the digitization is phase shifted sHghtly with respect to the readout of each sHce location.
- the subject was scanned while he lay prone in the magnet with his eyes open.
- the raw EEG data were processed as outlined above, first blanking the gradients and then removing the baUistocardiogram.
- the data were then submitted to a Fourier analysis and the energy at each of five pre-defined frequency bands was determined for each TR (that is, for each image time point). For each frequency band, therefore, it was possible to generate a separate time course.
- FIG 17A Shown in Figure 17A is the energy as a function of time in the alpha and theta bands for a four and a half minute period while the subject was at rest. The energy levels in these bands are largely independent for the first minute and a half of this session and seem to co-vary for the latter portion of the recording.
- Figure 17B demonstrates the estimated time course of the BOLD signal response based on the (untested) assumption that the hemodynamic response related to the EEG is similar to that seen in activation studies (e.g., (Cohen 1997; Cohen and DuBois 1999)).
- the alpha rhythm may weU be mtrinsic to the cerebral cortex (though influenced by thalamus) and is often interpreted as a sort of cortical resting state.
- scalp and presumably brain
- magnetic resonance signal intensity presumably indicating increases or decreases in cortical blood flow
- the signal processing methods developed to solve the problems in simultaneous EEG and fMRI have a broad range of additional appHcations.
- these methods can be used, without modification, to study electrical evoked potentials (Bonmassar, Anami et al. 1999) and their localization or, conversely, to yield a better temporal of the "single- trial" evoked MR response (Buckner, Bandettini etal. 1996).
- the clinical appHcations of fuUy integrated fMRI and EEG are substantial, the most immediate harvest will likely be in seizure source localization as an adjunct to other diagnostics used in surgical planning (Engel 2000).
- Interictal spikes are a common finding in the epUeptic brain yet even when recorded from the cortical surface do not provide definitive and reHable source localization. Recently, however, fast ripples have been reported in rat models and in the epUeptic human brain, that are associated closely with seizure foci (Bragin, Wilson et al. 2000). A stiU speculative interpretation is that such activity represents rapid spontaneous action potentials that are ultimately propagated and observable as interictal spikes. If so, this gives added hope to the idea that localization by fMRI of interictal spikes might become a reHable means to plan resective brain surgery, as has been proposed already by several others (e.g., (Warach, Levin et al.
- EEG recorded during MR imaging is particularly noisy, most of the artifacts contended with here are present in conventional EEG, though at lower ampUtude.
- the activity of the heart, particularly the electrocardiogram is a contaminant that can be corrected in the same manner as the minimization of the baUistocardiogram.
- AC line noise is often present as weU, and could be eliminated by using a triggered sampling approach, timing each block of samples to the power line osculations — these could be detected readily with a phase-locked loop.
- Digital sampling based on precise synchronization with weU characterized noise sources has appHcations that extend weU beyond Magnetic Resonance Imaging. It is clear, for example, that a comparable approach could be used to remove noise sources associated with AC power line osciUations; these include the removal of "hum" from digitaUy processed audio, Hght flicker from fluorescent iHumination, and other contaminants of digitaUy sampled signals.
- the scalp electrical potentials used in the EEG contain both time varying and static (DC) components. Often the DC offset is much larger than the EEG, but it is seldom of > interest for clinical diagnostic purposes, as it contains essentiaUy no information. However, it does cause trouble as it increases the dynamic range needed to digitize the EEG signal.
- the EEG may be only a few microVolts ( ⁇ V), while potentials of several tens of milHVolts (mV) may exist between electrodes, or as a result of the chemical electrode potential developed when the electrode is in contact with the scalp.
- the depth of the signal digitization wiU be reduced by the ratio of the EEG to the DC potential: assuming a DC offset of 10 mV, an EEG signal of 10 ⁇ V, and a 12 bit analog to digital converter (ADC), the 4096 different levels representable by the ADC wiU be reduced to only 4 levels for the EEG. Clearly this is unacceptable, as the quantization noise wiU dominate the recordings. For this reason, conventional EEG ampUfiers are equipped with AC-coupled (high pass) inputs', usuaUy a capacitor separating the output of the first stage ampUfiers from the input to the ADC.
- the inputs wiU typicaUy have time constants of several seconds, aUowing fluctuations of 1 Hz or so to pass without significant attenuation.
- One consequence of this AC coupling is that it creates a time constant for signal recovery if the input saturates and, because these filters must pass very low frequencies, the settling time for the analog signal to come back to the center of its nominal range can be quite long. For conventional EEG recordings this is acceptable, as large DC shifts are not frequent.
- the recovery time associated with AC coupling is a particular problem, as the gradient-induced artifacts can be large enough to bring the input stages into saturation.
- the ampUfiers that foUow can be pinned at either the positive or negative supply raUs for several milliseconds.
- the ampUfiers may require significant settling time such that the gradient artifact substantiaUy outlasts the gradient event.
- the invention provides a circuit that behaves as iUustrated in Figure 19.
- the high frequency common mode artifacts typicaUy present in EEG recorded in the MRI environment are first attenuated using passive components before being differentiaUy ampHfied.
- the differential ampHfication further reduces common mode artifacts.
- the data are then presented to a high pass ("anti-aHasing") filter to both attenuate the artifacts and to minimize the possibility of aUasing in analog to digital conversion.
- the final stage of this circuit performs the dual role of ampHfying the signal further and removing any DC bias. This is done by comparing the signal differentiaUy to a DC reference value sampled from the signal itself.
- the sampled signal is stored under a command here shown diagrammaticaUy as a switch connected to ground.
- Figure 14 shows one preferred implementation, in which the first high pass filter consists of a resistor-capacitor (RC) network 20.
- the first differential ampHfier is an integrated instrumentation ampHfier 28.
- the active high pass filter is a 24 dB Butterworth filter.
- the sample/hold circuit 24 stores the DC signal as charge across a low-leakage capacitor 37 with a compensation circuit for the non-zero input offset current found in the ampHfier that drives it.
- the final differential stage is a standard operational ampHfier 34.
- a pair of electrodes is connected to the differential input ampHfier stage 20, bypassed at high frequencies by paraUel capacitance.
- the input ampHfier 28 itself may be an integrated circuit, such as an INA114 from Burr- Brown corporation.
- a pair of matched 5 megohm resistors supply bias current to the integrated circuit.
- the output from the ampHfier is used as an input to a sample and hold stage 24, which utilizes an integrated circuit 36 such as the LF398 from National
- An offset trim adjustment 40 is provided to adjust for small static DC offsets in the integrated circuit 36.
- An input 39 is provided to sample and hold device 36.
- the DC corrected output is presented to the input of a low pass filter 22, that utilizes operational ampHfier integrated circuits 30 and 32, such as the TL072 from Texas Instruments corporation.
- the low- pass filter is implemented in a Chebyshev configuration to provide approximately 30 decibels/octave of attenuation at frequencies above the selected pass band.
- a final gain stage 26 includes a gain trimming adjustment 38 to match the gains across channels.
- the differential recording input stages are supplemented by a shield driver circuit that further reduces the sensitivity to electrical interference from the ongoing scanning activity, or other common mode signals in the leads.
- the output from the initial ampHfiers is coupled through an isolation ampHfier to reduce the possibility of circulating currents from the RF pulses saturating later components of the ampHfiers ( Figure 20).
- the DC offset signal is first detected digitaUy, then held in a digital latch, and returned via a digital to analog convertor as an offset correction. This offset correction is also coupled through an isolation ampHfier.
- the analog output is coupled opticaUy through the RF shielded room to minimize any problems of corruption of the MR signal through factors such as ground loops. Provision is made for the recording of multiple channels of EEG data through the use of multiplexing circuitry.
- Cerebrolysin enhances brain alpha activity and improves cognitive performance in elderly control subjects [In Process Citation]. Journal of Neural Transmission Supplement, 59: 315-28.
- Buckner RL Bandettini PA, O'Craven KM, Savoy RL, Petersen SE, Raichle ME, et al.
Abstract
Description
Claims
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2001283371A AU2001283371A1 (en) | 2000-08-15 | 2001-08-15 | Method and apparatus for reducing contamination of an electrical signal |
US10/344,776 US7286871B2 (en) | 2000-08-15 | 2001-08-15 | Method and apparatus for reducing contamination of an electrical signal |
JP2002518841A JP2004507293A (en) | 2000-08-15 | 2001-08-15 | Method and apparatus for reducing contamination of electrical signals |
EP01962172A EP1355571A2 (en) | 2000-08-15 | 2001-08-15 | Method and apparatus for reducing contamination of an electrical signal |
CA002418478A CA2418478A1 (en) | 2000-08-15 | 2001-08-15 | Method and apparatus for reducing contamination of an electrical signal |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US22538900P | 2000-08-15 | 2000-08-15 | |
US60/225,389 | 2000-08-15 | ||
US26733701P | 2001-02-07 | 2001-02-07 | |
US60/267,337 | 2001-02-07 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2002013689A2 true WO2002013689A2 (en) | 2002-02-21 |
WO2002013689A3 WO2002013689A3 (en) | 2003-08-21 |
Family
ID=26919552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2001/025480 WO2002013689A2 (en) | 2000-08-15 | 2001-08-15 | Method and apparatus for reducing contamination of an electrical signal |
Country Status (6)
Country | Link |
---|---|
US (1) | US7286871B2 (en) |
EP (1) | EP1355571A2 (en) |
JP (1) | JP2004507293A (en) |
AU (1) | AU2001283371A1 (en) |
CA (1) | CA2418478A1 (en) |
WO (1) | WO2002013689A2 (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003100450A1 (en) * | 2002-05-29 | 2003-12-04 | Physio-Tech Co., Ltd. | Simultaneous acquisition of bioelectric signals and magnetic resonance images |
EP1424637A1 (en) * | 2002-11-29 | 2004-06-02 | Instrumentarium Corporation | Artifact removal from an electric signal |
EP1488738A1 (en) * | 2003-06-19 | 2004-12-22 | Instrumentarium Corporation | Patient cable for medical measurements |
WO2005120339A1 (en) * | 2004-06-10 | 2005-12-22 | Unilever Plc | Apparatus and method for reducing interference |
WO2006099919A1 (en) * | 2005-03-03 | 2006-09-28 | Unilever Plc | Apparatus and method for reducing interference |
WO2007073576A1 (en) * | 2005-11-17 | 2007-07-05 | Brain Research Institute Pty Ltd | Apparatus and method for detection and monitoring of electrical activity and motion in the presence of a magnetic field |
US7291123B2 (en) * | 2003-06-04 | 2007-11-06 | Gambro Lundia | Joint for fluid transport lines for medical use |
WO2007140861A2 (en) * | 2006-06-02 | 2007-12-13 | Eberhard-Karls-Universität Tübingen | Medical electrode device |
EP1872715A1 (en) * | 2006-06-30 | 2008-01-02 | General Electric Company | Method and system for multi-channel biosignal processing |
WO2013028680A3 (en) * | 2011-08-25 | 2013-06-06 | Medtronic, Inc. | Method and apparatus for detecting a biomarker in the presence of electrical stimulation |
CN104095632A (en) * | 2013-04-07 | 2014-10-15 | 常州博睿康科技有限公司 | Method for processing electroencephalogram noise under nuclear magnetism |
Families Citing this family (206)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8762065B2 (en) | 1998-08-05 | 2014-06-24 | Cyberonics, Inc. | Closed-loop feedback-driven neuromodulation |
US9375573B2 (en) | 1998-08-05 | 2016-06-28 | Cyberonics, Inc. | Systems and methods for monitoring a patient's neurological disease state |
US9042988B2 (en) | 1998-08-05 | 2015-05-26 | Cyberonics, Inc. | Closed-loop vagus nerve stimulation |
US7747325B2 (en) | 1998-08-05 | 2010-06-29 | Neurovista Corporation | Systems and methods for monitoring a patient's neurological disease state |
US7209787B2 (en) | 1998-08-05 | 2007-04-24 | Bioneuronics Corporation | Apparatus and method for closed-loop intracranial stimulation for optimal control of neurological disease |
US9415222B2 (en) | 1998-08-05 | 2016-08-16 | Cyberonics, Inc. | Monitoring an epilepsy disease state with a supervisory module |
US9320900B2 (en) | 1998-08-05 | 2016-04-26 | Cyberonics, Inc. | Methods and systems for determining subject-specific parameters for a neuromodulation therapy |
WO2003073929A1 (en) * | 2002-02-28 | 2003-09-12 | General Hospital Corporation | Electroencephalograph sensor for use with magnetic resonance imaging and methods using such arrangements |
US7245786B2 (en) * | 2002-05-10 | 2007-07-17 | 976076 Alberta Inc. | Filtering artifact from fMRI data using the stockwell transform |
US7269456B2 (en) * | 2002-05-30 | 2007-09-11 | Collura Thomas F | Repetitive visual stimulation to EEG neurofeedback protocols |
US20060149139A1 (en) * | 2002-11-21 | 2006-07-06 | Giorgio Bonmassar | Apparatus and method for ascertaining and recording electrophysiological signals |
USRE45664E1 (en) * | 2003-03-19 | 2015-09-08 | Church Of Spiritual Technology | System for measuring and indicating changes in the resistance of a living body |
US20050059880A1 (en) * | 2003-09-11 | 2005-03-17 | Mathias Sanjay George | ECG driven image reconstruction for cardiac imaging |
US20050107674A1 (en) * | 2003-09-30 | 2005-05-19 | Jayant Parthasarathy | DC offset cancellation techniques |
US7460737B2 (en) | 2004-02-12 | 2008-12-02 | Hoshiko Llc | Method and apparatus for photograph finding |
US20050228306A1 (en) * | 2004-03-29 | 2005-10-13 | Vivosonic Inc. | System and method for filtering and detecting faint signals in noise |
US7844343B2 (en) * | 2004-03-30 | 2010-11-30 | Medtronic, Inc. | MRI-safe implantable medical device |
US7844344B2 (en) | 2004-03-30 | 2010-11-30 | Medtronic, Inc. | MRI-safe implantable lead |
US20090012387A1 (en) * | 2004-05-25 | 2009-01-08 | Hvidovre Hospital | Encoding and transmission of signals as rf signals for detection using an mr apparatus |
JP2006087526A (en) * | 2004-09-21 | 2006-04-06 | National Institute Of Information & Communication Technology | Biological activity measuring instrument |
US7378844B2 (en) * | 2004-09-30 | 2008-05-27 | General Electric Company | Magnetic resonance system, receiver & method of generating detecting and digitizing analog MR signals solely within the MR shielded environment |
CA2588877A1 (en) * | 2004-10-23 | 2006-05-04 | John O. Bell | Passive monitoring of bioelectical signals and active electrical anesthesia stimulation |
US7545272B2 (en) | 2005-02-08 | 2009-06-09 | Therasense, Inc. | RF tag on test strips, test strip vials and boxes |
US10327701B2 (en) | 2005-05-06 | 2019-06-25 | The General Hospital Corporation | Apparatuses and methods for electrophysiological signal delivery and recording during MRI |
DE102005027438B4 (en) * | 2005-06-14 | 2011-12-22 | Siemens Ag | Method for ECG triggering a measurement sequence of a magnetic resonance device |
DE102005031737B4 (en) * | 2005-07-07 | 2008-04-24 | Mfd Diagnostics Gmbh | Arrangement for the local detection and evaluation of the brain activity of living organisms |
WO2007041766A1 (en) * | 2005-10-10 | 2007-04-19 | Compumedics Limited | Adaptive real-time line noise suppression for electrical or magnetic physiological signals |
CN100423688C (en) | 2005-10-19 | 2008-10-08 | 深圳迈瑞生物医疗电子股份有限公司 | Method and apparatus for inhibiting power frequency common-mode interference |
WO2007102893A2 (en) * | 2005-11-11 | 2007-09-13 | Greatbatch Ltd. | Tank filters placed in series with the lead wires or circuits of active medical devices to enhance mri compatibility |
US20070135727A1 (en) * | 2005-12-12 | 2007-06-14 | Juha Virtanen | Detection of artifacts in bioelectric signals |
US7684856B2 (en) * | 2005-12-12 | 2010-03-23 | General Electric Company | Detection of artifacts in bioelectric signals |
US8868172B2 (en) | 2005-12-28 | 2014-10-21 | Cyberonics, Inc. | Methods and systems for recommending an appropriate action to a patient for managing epilepsy and other neurological disorders |
US8725243B2 (en) | 2005-12-28 | 2014-05-13 | Cyberonics, Inc. | Methods and systems for recommending an appropriate pharmacological treatment to a patient for managing epilepsy and other neurological disorders |
US8406864B2 (en) * | 2006-02-07 | 2013-03-26 | Impulse Dynamics Nv | Assessing cardiac activity |
EP2034885A4 (en) | 2006-06-23 | 2010-12-01 | Neurovista Corp | Minimally invasive monitoring systems and methods |
US7729890B2 (en) * | 2006-08-22 | 2010-06-01 | Analog Devices, Inc. | Method for determining the change of a signal, and an apparatus including a circuit arranged to implement the method |
US8295934B2 (en) | 2006-11-14 | 2012-10-23 | Neurovista Corporation | Systems and methods of reducing artifact in neurological stimulation systems |
US20080159365A1 (en) * | 2006-12-22 | 2008-07-03 | Branislav Dubocanin | Analog Conditioning of Bioelectric Signals |
US20080154116A1 (en) * | 2006-12-22 | 2008-06-26 | Duensing G Randy | Method and apparatus for obtaining electrocardiogram (ECG) signals |
US7847552B2 (en) * | 2007-01-10 | 2010-12-07 | General Electric Company | Exclusion of compromised PET data during simultaneous PET-MR acquisition |
WO2008092119A2 (en) | 2007-01-25 | 2008-07-31 | Neurovista Corporation | Systems and methods for identifying a contra-ictal condition in a subject |
EP2124734A2 (en) | 2007-01-25 | 2009-12-02 | NeuroVista Corporation | Methods and systems for measuring a subject's susceptibility to a seizure |
US9044593B2 (en) | 2007-02-14 | 2015-06-02 | Medtronic, Inc. | Discontinuous conductive filler polymer-matrix composites for electromagnetic shielding |
US20080208008A1 (en) * | 2007-02-28 | 2008-08-28 | Microsoft Corporation | Wireless body sensor with small size background |
US8036736B2 (en) | 2007-03-21 | 2011-10-11 | Neuro Vista Corporation | Implantable systems and methods for identifying a contra-ictal condition in a subject |
KR20080086055A (en) * | 2007-03-21 | 2008-09-25 | 한국과학기술원 | Brain wave measuring method, apparatus and computer readable medium on which program for executing the method is recorded |
CN101711388B (en) | 2007-03-29 | 2016-04-27 | 神经焦点公司 | The effect analysis of marketing and amusement |
US7996076B2 (en) * | 2007-04-02 | 2011-08-09 | The Regents Of The University Of Michigan | Automated polysomnographic assessment for rapid eye movement sleep behavior disorder |
US8483842B2 (en) | 2007-04-25 | 2013-07-09 | Medtronic, Inc. | Lead or lead extension having a conductive body and conductive body contact |
DE102007019962B3 (en) * | 2007-04-27 | 2008-04-03 | Dräger Medical AG & Co. KG | Y-piece for medical respiration systems, has contactless connection element for connection to inspiration branch and/or expiration branch, with breathing tube |
EP2142082A4 (en) | 2007-05-01 | 2015-10-28 | Neurofocus Inc | Neuro-informatics repository system |
WO2008137581A1 (en) | 2007-05-01 | 2008-11-13 | Neurofocus, Inc. | Neuro-feedback based stimulus compression device |
WO2008137161A2 (en) * | 2007-05-07 | 2008-11-13 | The Trustees Of The University Of Pennsylvania | Concurrent electrophysiologic and cerebral blood flow neuroimaging |
US8392253B2 (en) | 2007-05-16 | 2013-03-05 | The Nielsen Company (Us), Llc | Neuro-physiology and neuro-behavioral based stimulus targeting system |
US8694070B2 (en) | 2007-05-22 | 2014-04-08 | Persyst Development Corporation | Electrode applicator for quick press on EEG electrode |
US8494905B2 (en) * | 2007-06-06 | 2013-07-23 | The Nielsen Company (Us), Llc | Audience response analysis using simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) |
US9788744B2 (en) | 2007-07-27 | 2017-10-17 | Cyberonics, Inc. | Systems for monitoring brain activity and patient advisory device |
US8533042B2 (en) | 2007-07-30 | 2013-09-10 | The Nielsen Company (Us), Llc | Neuro-response stimulus and stimulus attribute resonance estimator |
KR20100047865A (en) | 2007-08-28 | 2010-05-10 | 뉴로포커스, 인크. | Consumer experience assessment system |
US8386313B2 (en) | 2007-08-28 | 2013-02-26 | The Nielsen Company (Us), Llc | Stimulus placement system using subject neuro-response measurements |
US8635105B2 (en) | 2007-08-28 | 2014-01-21 | The Nielsen Company (Us), Llc | Consumer experience portrayal effectiveness assessment system |
US8392255B2 (en) | 2007-08-29 | 2013-03-05 | The Nielsen Company (Us), Llc | Content based selection and meta tagging of advertisement breaks |
US20090062680A1 (en) * | 2007-09-04 | 2009-03-05 | Brain Train | Artifact detection and correction system for electroencephalograph neurofeedback training methodology |
US8494610B2 (en) * | 2007-09-20 | 2013-07-23 | The Nielsen Company (Us), Llc | Analysis of marketing and entertainment effectiveness using magnetoencephalography |
US20090083129A1 (en) | 2007-09-20 | 2009-03-26 | Neurofocus, Inc. | Personalized content delivery using neuro-response priming data |
TWI400062B (en) * | 2007-11-08 | 2013-07-01 | 私立中原大學 | Medical devices that record physiological signals |
US20090171168A1 (en) | 2007-12-28 | 2009-07-02 | Leyde Kent W | Systems and Method for Recording Clinical Manifestations of a Seizure |
US9259591B2 (en) | 2007-12-28 | 2016-02-16 | Cyberonics, Inc. | Housing for an implantable medical device |
US9037263B2 (en) | 2008-03-12 | 2015-05-19 | Medtronic, Inc. | System and method for implantable medical device lead shielding |
WO2009117102A2 (en) * | 2008-03-17 | 2009-09-24 | Npe Systems, Inc. | Background light detection system for a flow cytometer |
WO2009134763A1 (en) * | 2008-04-29 | 2009-11-05 | Board of Governors for Higher Education, State of Rhode Island and the Providence Plantations | Biomedical sensors usable on un-prepared contact surfaces |
KR100968505B1 (en) * | 2008-09-08 | 2010-07-07 | 한국과학기술원 | Metal supported solid oxide fuel cell and manufacturing method thereof |
US8059208B2 (en) * | 2008-09-26 | 2011-11-15 | Pelco, Inc. | Line tapping method of video amplification for unshielded twisted pair active receiver |
US8870780B2 (en) | 2008-10-15 | 2014-10-28 | The Board Of Trustees Of The Leland Stanford Junior University | Systems and methods for monitoring heart function |
US8849390B2 (en) | 2008-12-29 | 2014-09-30 | Cyberonics, Inc. | Processing for multi-channel signals |
US8588933B2 (en) | 2009-01-09 | 2013-11-19 | Cyberonics, Inc. | Medical lead termination sleeve for implantable medical devices |
US8464288B2 (en) | 2009-01-21 | 2013-06-11 | The Nielsen Company (Us), Llc | Methods and apparatus for providing personalized media in video |
US9357240B2 (en) | 2009-01-21 | 2016-05-31 | The Nielsen Company (Us), Llc | Methods and apparatus for providing alternate media for video decoders |
US8270814B2 (en) | 2009-01-21 | 2012-09-18 | The Nielsen Company (Us), Llc | Methods and apparatus for providing video with embedded media |
DE102009010501A1 (en) * | 2009-02-25 | 2010-09-09 | Siemens Aktiengesellschaft | Method for noise reduction of CT image data and image processing system |
KR101007965B1 (en) | 2009-02-27 | 2011-01-14 | 고려대학교 산학협력단 | Apparatus for processing of eeg inputted through single channel and processing method using the same |
JP5424676B2 (en) * | 2009-03-13 | 2014-02-26 | キヤノン株式会社 | Image processing device |
US8155736B2 (en) * | 2009-03-16 | 2012-04-10 | Neurosky, Inc. | EEG control of devices using sensory evoked potentials |
US20100250325A1 (en) | 2009-03-24 | 2010-09-30 | Neurofocus, Inc. | Neurological profiles for market matching and stimulus presentation |
US20100249538A1 (en) * | 2009-03-24 | 2010-09-30 | Neurofocus, Inc. | Presentation measure using neurographics |
US20100249635A1 (en) * | 2009-03-26 | 2010-09-30 | Cordial Medical Europe B.V. | Hearing screening system for a subject or a patient, and a method for hearing screening |
US20100249636A1 (en) * | 2009-03-27 | 2010-09-30 | Neurofocus, Inc. | Personalized stimulus placement in video games |
US9186499B2 (en) | 2009-04-30 | 2015-11-17 | Medtronic, Inc. | Grounding of a shield within an implantable medical lead |
US9135731B2 (en) | 2009-05-21 | 2015-09-15 | Navico Holding As | Systems, devices, methods for sensing and processing fishing related data |
US8786624B2 (en) | 2009-06-02 | 2014-07-22 | Cyberonics, Inc. | Processing for multi-channel signals |
WO2010143118A1 (en) * | 2009-06-12 | 2010-12-16 | Koninklijke Philips Electronics N.V. | System and method for biometric identification |
US8655437B2 (en) | 2009-08-21 | 2014-02-18 | The Nielsen Company (Us), Llc | Analysis of the mirror neuron system for evaluation of stimulus |
US10987015B2 (en) | 2009-08-24 | 2021-04-27 | Nielsen Consumer Llc | Dry electrodes for electroencephalography |
US9560984B2 (en) | 2009-10-29 | 2017-02-07 | The Nielsen Company (Us), Llc | Analysis of controlled and automatic attention for introduction of stimulus material |
US20110106750A1 (en) | 2009-10-29 | 2011-05-05 | Neurofocus, Inc. | Generating ratings predictions using neuro-response data |
US8209224B2 (en) | 2009-10-29 | 2012-06-26 | The Nielsen Company (Us), Llc | Intracluster content management using neuro-response priming data |
US8335715B2 (en) | 2009-11-19 | 2012-12-18 | The Nielsen Company (Us), Llc. | Advertisement exchange using neuro-response data |
US8335716B2 (en) | 2009-11-19 | 2012-12-18 | The Nielsen Company (Us), Llc. | Multimedia advertisement exchange |
US9643019B2 (en) | 2010-02-12 | 2017-05-09 | Cyberonics, Inc. | Neurological monitoring and alerts |
US8684742B2 (en) | 2010-04-19 | 2014-04-01 | Innerscope Research, Inc. | Short imagery task (SIT) research method |
US8655428B2 (en) | 2010-05-12 | 2014-02-18 | The Nielsen Company (Us), Llc | Neuro-response data synchronization |
US8040135B1 (en) * | 2010-05-27 | 2011-10-18 | The Board Of Trustees Of The Leland Stanford Junior University | Contrast and resolution enhancement with signal compensation |
US20110295080A1 (en) * | 2010-05-30 | 2011-12-01 | Ralink Technology Corporation | Physiology Condition Detection Device and the System Thereof |
US8494624B2 (en) * | 2010-06-25 | 2013-07-23 | Electrical Geodesics, Inc. | Method and apparatus for reducing noise in brain signal measurements |
US8392250B2 (en) | 2010-08-09 | 2013-03-05 | The Nielsen Company (Us), Llc | Neuro-response evaluated stimulus in virtual reality environments |
US8392251B2 (en) | 2010-08-09 | 2013-03-05 | The Nielsen Company (Us), Llc | Location aware presentation of stimulus material |
US8396744B2 (en) | 2010-08-25 | 2013-03-12 | The Nielsen Company (Us), Llc | Effective virtual reality environments for presentation of marketing materials |
US20130310422A1 (en) | 2010-09-01 | 2013-11-21 | The General Hospital Corporation | Reversal of general anesthesia by administration of methylphenidate, amphetamine, modafinil, amantadine, and/or caffeine |
US8892184B2 (en) | 2010-10-18 | 2014-11-18 | Siemens Medical Solutions Usa, Inc. | Systems and methods for reducing interference in a dual modality imaging system |
US9011346B2 (en) | 2011-01-27 | 2015-04-21 | The Board Of Trustees Of The Leland Stanford Junior University | Systems and methods for monitoring the circulatory system |
US20130331727A1 (en) * | 2011-01-28 | 2013-12-12 | Agency For Science, Technology And Research | Method and system for detecting attention |
US8600492B2 (en) * | 2011-01-31 | 2013-12-03 | Brain Products Gmbh | Method and system for EEG artifact correction in combined EEG and functional MRI recording |
JP2012187166A (en) * | 2011-03-09 | 2012-10-04 | Kiyohara Optics Inc | Cutaneous conductor for mri or meg and method of measuring cutaneous resistance of subject |
US9486332B2 (en) | 2011-04-15 | 2016-11-08 | The Johns Hopkins University | Multi-modal neural interfacing for prosthetic devices |
US9833168B2 (en) | 2011-06-06 | 2017-12-05 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Noise tolerant localization systems and methods |
US9664765B2 (en) * | 2011-06-30 | 2017-05-30 | Hitachi, Ltd. | Magnetic resonance imaging apparatus and gradient magnetic field waveform estimation method |
US20130096440A1 (en) * | 2011-10-14 | 2013-04-18 | Marianna Kiraly | Portable fetal eeg-recording device and method of use |
US9055927B2 (en) | 2011-11-25 | 2015-06-16 | Persyst Development Corporation | User interface for artifact removal in an EEG |
US8666484B2 (en) | 2011-11-25 | 2014-03-04 | Persyst Development Corporation | Method and system for displaying EEG data |
JP6023210B2 (en) | 2011-11-26 | 2016-11-09 | パーシスト ディベロップメント コーポレーション | Method and system for detecting and removing EEG artifacts |
US9569986B2 (en) | 2012-02-27 | 2017-02-14 | The Nielsen Company (Us), Llc | System and method for gathering and analyzing biometric user feedback for use in social media and advertising applications |
US9042623B2 (en) | 2012-03-07 | 2015-05-26 | Vanderbilt University | Automatic detection and correction of magnetic resonance imaging data |
DE102012205629A1 (en) * | 2012-04-05 | 2013-10-10 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Method and magnetic resonance system for functional MR imaging of a predetermined volume portion of a brain of a living examination subject |
WO2013158189A1 (en) | 2012-04-19 | 2013-10-24 | Medtronic, Inc. | Paired medical lead bodies with braided conductive shields having different physical parameter values |
WO2013175349A2 (en) * | 2012-05-24 | 2013-11-28 | Koninklijke Philips N.V. | Reduction of mri interference from the electrocardiogram using lead information |
US9104697B2 (en) | 2012-07-24 | 2015-08-11 | Navico Holdings As | Aquatic geographic information system |
US11083401B2 (en) | 2012-08-09 | 2021-08-10 | Northeastern University | Electric field encephalography: electric field based brain signal detection and monitoring |
US8989835B2 (en) | 2012-08-17 | 2015-03-24 | The Nielsen Company (Us), Llc | Systems and methods to gather and analyze electroencephalographic data |
US20140066739A1 (en) * | 2012-08-29 | 2014-03-06 | Bin He | System and method for quantifying or imaging pain using electrophysiological measurements |
JP6003470B2 (en) * | 2012-09-25 | 2016-10-05 | オムロンヘルスケア株式会社 | Blood pressure measurement device and pulse wave detection method |
US9351653B1 (en) * | 2012-11-29 | 2016-05-31 | Intan Technologies, LLC | Multi-channel reconfigurable systems and methods for sensing biopotential signals |
US9320450B2 (en) | 2013-03-14 | 2016-04-26 | The Nielsen Company (Us), Llc | Methods and apparatus to gather and analyze electroencephalographic data |
WO2014169100A1 (en) * | 2013-04-10 | 2014-10-16 | The General Hospital Corporation | Removing electrophysiologic artifacts from a magnetic resonance imaging system |
JP6503347B2 (en) * | 2013-06-21 | 2019-04-17 | ノースイースタン ユニバーシティ | Sensor assembly for measurement of electrical activity of the brain including electric field electroencephalogram |
US10314503B2 (en) | 2013-06-27 | 2019-06-11 | The General Hospital Corporation | Systems and methods for tracking non-stationary spectral structure and dynamics in physiological data |
US10383574B2 (en) | 2013-06-28 | 2019-08-20 | The General Hospital Corporation | Systems and methods to infer brain state during burst suppression |
US10327663B2 (en) * | 2013-08-31 | 2019-06-25 | Alpha Omega Neuro Technologies Ltd. | Evoked response probe and method of use |
WO2015038969A1 (en) | 2013-09-13 | 2015-03-19 | The General Hospital Corporation | Systems and methods for improved brain monitoring during general anesthesia and sedation |
EP3072266B1 (en) * | 2013-11-19 | 2019-01-09 | The Regents of The University of California | A saturation-tolerant electrophysical recording interface |
US9993638B2 (en) | 2013-12-14 | 2018-06-12 | Medtronic, Inc. | Devices, systems and methods to reduce coupling of a shield and a conductor within an implantable medical lead |
US10929753B1 (en) | 2014-01-20 | 2021-02-23 | Persyst Development Corporation | System and method for generating a probability value for an event |
US9320447B2 (en) | 2014-02-07 | 2016-04-26 | Biosense Webster (Israel) Ltd. | Dynamic cancellation of MRI sequencing noise appearing in an ECG signal |
US20150223712A1 (en) | 2014-02-07 | 2015-08-13 | Biosense Webster (Israel) Ltd. | Analog cancellation of mri sequencing noise appearing in an ecg signal |
US9747421B2 (en) * | 2014-02-10 | 2017-08-29 | Picofemto LLC | Multi-factor brain analysis via medical imaging decision support systems and methods |
US9622702B2 (en) | 2014-04-03 | 2017-04-18 | The Nielsen Company (Us), Llc | Methods and apparatus to gather and analyze electroencephalographic data |
US9568354B2 (en) | 2014-06-12 | 2017-02-14 | PhysioWave, Inc. | Multifunction scale with large-area display |
US10130273B2 (en) | 2014-06-12 | 2018-11-20 | PhysioWave, Inc. | Device and method having automatic user-responsive and user-specific physiological-meter platform |
US9949662B2 (en) | 2014-06-12 | 2018-04-24 | PhysioWave, Inc. | Device and method having automatic user recognition and obtaining impedance-measurement signals |
US9943241B2 (en) | 2014-06-12 | 2018-04-17 | PhysioWave, Inc. | Impedance measurement devices, systems, and methods |
US9546898B2 (en) | 2014-06-12 | 2017-01-17 | PhysioWave, Inc. | Fitness testing scale |
WO2016014427A1 (en) | 2014-07-23 | 2016-01-28 | Medtronic, Inc. | Methods of shielding implantable medical leads and implantable medical lead extensions |
US10155111B2 (en) | 2014-07-24 | 2018-12-18 | Medtronic, Inc. | Methods of shielding implantable medical leads and implantable medical lead extensions |
US9770185B2 (en) * | 2014-08-06 | 2017-09-26 | Verily Life Sciences Llc | Sharing a single electrode between skin resistance and capacitance measurements |
US9693696B2 (en) | 2014-08-07 | 2017-07-04 | PhysioWave, Inc. | System with user-physiological data updates |
US9498137B2 (en) | 2014-08-07 | 2016-11-22 | PhysioWave, Inc. | Multi-function fitness scale with display |
CA2960189C (en) | 2014-09-05 | 2021-11-23 | Hyperfine Research, Inc. | Low field magnetic resonance imaging methods and apparatus |
US9459201B2 (en) | 2014-09-29 | 2016-10-04 | Zyomed Corp. | Systems and methods for noninvasive blood glucose and other analyte detection and measurement using collision computing |
WO2016077417A1 (en) * | 2014-11-11 | 2016-05-19 | Hyperfine Research, Inc. | Low field magnetic resonance methods and apparatus |
US9562840B2 (en) * | 2014-12-03 | 2017-02-07 | Cambridge Viscosity, Inc. | High precision reciprocating bob viscometer |
US20160213277A1 (en) * | 2014-12-23 | 2016-07-28 | Phoenix Children's Hospital, Inc. | Programmable Multimodal Stimulator for Cortical Mapping |
US10010259B2 (en) | 2015-05-01 | 2018-07-03 | Advancer Technologies, Llc | EMG circuit |
US9936250B2 (en) | 2015-05-19 | 2018-04-03 | The Nielsen Company (Us), Llc | Methods and apparatus to adjust content presented to an individual |
US10542961B2 (en) | 2015-06-15 | 2020-01-28 | The Research Foundation For The State University Of New York | System and method for infrasonic cardiac monitoring |
US10945671B2 (en) | 2015-06-23 | 2021-03-16 | PhysioWave, Inc. | Determining physiological parameters using movement detection |
JP6105687B2 (en) * | 2015-07-28 | 2017-03-29 | 株式会社 清原光学 | Ability to judge device |
US11561126B2 (en) | 2015-11-20 | 2023-01-24 | PhysioWave, Inc. | Scale-based user-physiological heuristic systems |
US10980483B2 (en) | 2015-11-20 | 2021-04-20 | PhysioWave, Inc. | Remote physiologic parameter determination methods and platform apparatuses |
US10395055B2 (en) | 2015-11-20 | 2019-08-27 | PhysioWave, Inc. | Scale-based data access control methods and apparatuses |
US10553306B2 (en) | 2015-11-20 | 2020-02-04 | PhysioWave, Inc. | Scaled-based methods and apparatuses for automatically updating patient profiles |
US10923217B2 (en) | 2015-11-20 | 2021-02-16 | PhysioWave, Inc. | Condition or treatment assessment methods and platform apparatuses |
US10436630B2 (en) | 2015-11-20 | 2019-10-08 | PhysioWave, Inc. | Scale-based user-physiological data hierarchy service apparatuses and methods |
US10319144B2 (en) | 2015-12-22 | 2019-06-11 | The Regents Of The University Of California | Computational localization of fibrillation sources |
US9554738B1 (en) | 2016-03-30 | 2017-01-31 | Zyomed Corp. | Spectroscopic tomography systems and methods for noninvasive detection and measurement of analytes using collision computing |
JP6785446B2 (en) * | 2016-04-14 | 2020-11-18 | パナソニックIpマネジメント株式会社 | Biological signal measurement system |
US10390772B1 (en) | 2016-05-04 | 2019-08-27 | PhysioWave, Inc. | Scale-based on-demand care system |
US10365366B2 (en) * | 2016-06-21 | 2019-07-30 | Navico Holding As | Adjustable range viewing of sonar imagery during trip replay |
US9986197B2 (en) | 2016-08-24 | 2018-05-29 | Navico Holding As | Trip replay experience |
CN106443152A (en) * | 2016-09-05 | 2017-02-22 | 东南大学 | High-precision rural network low-voltage transient residual current detection method |
US10215619B1 (en) | 2016-09-06 | 2019-02-26 | PhysioWave, Inc. | Scale-based time synchrony |
KR101870758B1 (en) * | 2016-10-13 | 2018-06-26 | (주)로임시스템 | Bio-signal detection apparatus for bio-signal interference identification |
US10716514B1 (en) | 2017-04-10 | 2020-07-21 | Hrl Laboratories, Llc | System and method for optimized independent component selection for automated signal artifact removal to generate a clean signal |
US11344723B1 (en) | 2016-10-24 | 2022-05-31 | Hrl Laboratories, Llc | System and method for decoding and behaviorally validating memory consolidation during sleep from EEG after waking experience |
US10627464B2 (en) | 2016-11-22 | 2020-04-21 | Hyperfine Research, Inc. | Low-field magnetic resonance imaging methods and apparatus |
US10539637B2 (en) | 2016-11-22 | 2020-01-21 | Hyperfine Research, Inc. | Portable magnetic resonance imaging methods and apparatus |
WO2018102402A1 (en) | 2016-11-29 | 2018-06-07 | The General Hospital Corporation | Systems and methods for analyzing electrophysiological data from patients undergoing medical treatments |
CN110267585B (en) * | 2016-12-08 | 2022-05-31 | 皇家飞利浦有限公司 | Determining corrected timing of stimulation provided to a subject during sleep |
CN107095672B (en) * | 2017-06-10 | 2020-11-03 | 深圳创达云睿智能科技有限公司 | Method for removing artifact signals in electroencephalogram signals |
JP6996135B2 (en) * | 2017-07-03 | 2022-01-17 | 株式会社リコー | Information processing equipment, information processing methods, and programs |
US10555670B2 (en) | 2017-07-10 | 2020-02-11 | International Business Machines Corporation | Adaptive filtration of sweat artifacts during electronic brain monitoring |
WO2019060298A1 (en) | 2017-09-19 | 2019-03-28 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement |
US10581208B2 (en) | 2017-10-16 | 2020-03-03 | Equipement Electroline Inc. | Ethernet connector with electromagnetic filtering |
US11717686B2 (en) | 2017-12-04 | 2023-08-08 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to facilitate learning and performance |
WO2019133997A1 (en) | 2017-12-31 | 2019-07-04 | Neuroenhancement Lab, LLC | System and method for neuroenhancement to enhance emotional response |
WO2019146882A1 (en) * | 2018-01-23 | 2019-08-01 | 인제대학교 산학협력단 | Method for diagnosis and prediction of mental disorders by using heart rate variability and quantitative electroencephalography and electronic apparatus thereof |
US11219416B2 (en) | 2018-03-12 | 2022-01-11 | Persyst Development Corporation | Graphically displaying evoked potentials |
US11457855B2 (en) | 2018-03-12 | 2022-10-04 | Persyst Development Corporation | Method and system for utilizing empirical null hypothesis for a biological time series |
US11364361B2 (en) | 2018-04-20 | 2022-06-21 | Neuroenhancement Lab, LLC | System and method for inducing sleep by transplanting mental states |
CN112639998A (en) | 2018-07-05 | 2021-04-09 | 加利福尼亚大学董事会 | Computational simulation of anatomical structures and body surface electrode positioning |
EP3849410A4 (en) | 2018-09-14 | 2022-11-02 | Neuroenhancement Lab, LLC | System and method of improving sleep |
FR3087332B1 (en) * | 2018-10-23 | 2020-11-06 | Schiller Medical | REAL-TIME CORRECTION PROCESS OF AT LEAST ONE ELECTRO-PHYSIOLOGICAL SIGNAL. |
CN109318207B (en) * | 2018-11-07 | 2020-11-10 | 西安交通大学 | Myoelectric timing lower limb movement preparatory potential detection system and method |
US11786694B2 (en) | 2019-05-24 | 2023-10-17 | NeuroLight, Inc. | Device, method, and app for facilitating sleep |
US11583218B2 (en) | 2019-11-20 | 2023-02-21 | Advancer Technologies, Llc | EMG device |
USD1015545S1 (en) | 2019-11-20 | 2024-02-20 | Advancer Technologies, Llc | Electromyography device |
RU2733111C1 (en) * | 2020-03-26 | 2020-09-29 | федеральное государственное автономное образовательное учреждение высшего образования «Национальный исследовательский Томский политехнический университет» | Method for frequency-time correlation analysis of digital signals |
US11349263B1 (en) * | 2021-05-07 | 2022-05-31 | Rockwell Collins, Inc. | Electromagnetic interference shielding for a coaxial connector using a gasket assembly |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4887609A (en) * | 1987-05-13 | 1989-12-19 | The Methodist Hospital System | Apparatus and method for filtering electrocardiograph signals |
US5217010A (en) * | 1991-05-28 | 1993-06-08 | The Johns Hopkins University | Ecg amplifier and cardiac pacemaker for use during magnetic resonance imaging |
US5230344A (en) * | 1992-07-31 | 1993-07-27 | Intelligent Hearing Systems Corp. | Evoked potential processing system with spectral averaging, adaptive averaging, two dimensional filters, electrode configuration and method therefor |
WO1994027497A1 (en) * | 1993-06-02 | 1994-12-08 | Bang & Olufsen Technology A/S | Apparatus for measuring cardiac signals |
US5513649A (en) * | 1994-03-22 | 1996-05-07 | Sam Technology, Inc. | Adaptive interference canceler for EEG movement and eye artifacts |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4781200A (en) * | 1985-10-04 | 1988-11-01 | Baker Donald A | Ambulatory non-invasive automatic fetal monitoring system |
CA2042271A1 (en) * | 1990-08-06 | 1992-02-07 | Kevin Scott Kreger | Bio-potential signal processor for magnetic resonance imaging |
WO1992015955A1 (en) * | 1991-03-07 | 1992-09-17 | Vital Signals, Inc. | Signal processing apparatus and method |
US5341811A (en) * | 1991-03-26 | 1994-08-30 | Allegheny-Singer Research Institute | Method and apparatus for observation of ventricular late potentials |
WO1993003544A1 (en) * | 1991-08-09 | 1993-02-18 | Siemens Aktiengesellschaft | Process for compensating for quasi-periodic disturbances of measurement signals |
US5445162A (en) * | 1993-08-27 | 1995-08-29 | Beth Israel Hospital Association | Apparatus and method for recording an electroencephalogram during magnetic resonance imaging |
US5713367A (en) * | 1994-01-26 | 1998-02-03 | Cambridge Heart, Inc. | Measuring and assessing cardiac electrical stability |
US5532944A (en) * | 1994-07-28 | 1996-07-02 | Sorrento Electronics, Inc. | Multi-channel analysis system and method using digital signal processing |
US5742900A (en) * | 1994-09-30 | 1998-04-21 | Comsat Corporation | Voltage-variable biased inverting limiter for RFI suppression overview |
US5966684A (en) * | 1995-11-09 | 1999-10-12 | Richardson; Thomas L. | Method and apparatus for cancelling periodic electrical interference |
US5748507A (en) * | 1996-08-06 | 1998-05-05 | Lockheed Martin Corporation | Real-time superresolution signal processing |
US6181961B1 (en) * | 1997-12-16 | 2001-01-30 | Richard L. Prass | Method and apparatus for an automatic setup of a multi-channel nerve integrity monitoring system |
US7572231B2 (en) * | 2000-01-31 | 2009-08-11 | Pearlman Justin D | Method of and system for signal separation during multivariate physiological monitoring |
WO2002002964A1 (en) * | 2000-06-30 | 2002-01-10 | Vistek Inc. | Mechanical signal filter |
EP1417927A1 (en) * | 2002-11-11 | 2004-05-12 | Schiller AG | Method and device for the detection and transmission of electrophysiological signals for use in a MRI system |
EP1424637A1 (en) * | 2002-11-29 | 2004-06-02 | Instrumentarium Corporation | Artifact removal from an electric signal |
US20050277826A1 (en) * | 2004-06-10 | 2005-12-15 | Conopco, Inc. | Apparatus and method for reducing interference |
-
2001
- 2001-08-15 JP JP2002518841A patent/JP2004507293A/en not_active Withdrawn
- 2001-08-15 CA CA002418478A patent/CA2418478A1/en not_active Abandoned
- 2001-08-15 EP EP01962172A patent/EP1355571A2/en not_active Withdrawn
- 2001-08-15 US US10/344,776 patent/US7286871B2/en active Active
- 2001-08-15 WO PCT/US2001/025480 patent/WO2002013689A2/en not_active Application Discontinuation
- 2001-08-15 AU AU2001283371A patent/AU2001283371A1/en not_active Abandoned
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4887609A (en) * | 1987-05-13 | 1989-12-19 | The Methodist Hospital System | Apparatus and method for filtering electrocardiograph signals |
US5217010A (en) * | 1991-05-28 | 1993-06-08 | The Johns Hopkins University | Ecg amplifier and cardiac pacemaker for use during magnetic resonance imaging |
US5230344A (en) * | 1992-07-31 | 1993-07-27 | Intelligent Hearing Systems Corp. | Evoked potential processing system with spectral averaging, adaptive averaging, two dimensional filters, electrode configuration and method therefor |
WO1994027497A1 (en) * | 1993-06-02 | 1994-12-08 | Bang & Olufsen Technology A/S | Apparatus for measuring cardiac signals |
US5513649A (en) * | 1994-03-22 | 1996-05-07 | Sam Technology, Inc. | Adaptive interference canceler for EEG movement and eye artifacts |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003100450A1 (en) * | 2002-05-29 | 2003-12-04 | Physio-Tech Co., Ltd. | Simultaneous acquisition of bioelectric signals and magnetic resonance images |
EP1424637A1 (en) * | 2002-11-29 | 2004-06-02 | Instrumentarium Corporation | Artifact removal from an electric signal |
US7038601B2 (en) | 2002-11-29 | 2006-05-02 | Instrumentatium Corp. | Artifact removal from an electric signal |
US7291123B2 (en) * | 2003-06-04 | 2007-11-06 | Gambro Lundia | Joint for fluid transport lines for medical use |
US7785284B2 (en) | 2003-06-04 | 2010-08-31 | Gambro Lundia Ab | Joint for fluid transport lines for medical use |
US7294785B2 (en) | 2003-06-19 | 2007-11-13 | Ge Healthcare Finland Oy | Patient cable for medical measurements |
EP1488738A1 (en) * | 2003-06-19 | 2004-12-22 | Instrumentarium Corporation | Patient cable for medical measurements |
WO2005120339A1 (en) * | 2004-06-10 | 2005-12-22 | Unilever Plc | Apparatus and method for reducing interference |
AU2005251456B2 (en) * | 2004-06-10 | 2009-06-18 | Unilever Plc | Apparatus and method for reducing interference |
WO2006099919A1 (en) * | 2005-03-03 | 2006-09-28 | Unilever Plc | Apparatus and method for reducing interference |
WO2007073576A1 (en) * | 2005-11-17 | 2007-07-05 | Brain Research Institute Pty Ltd | Apparatus and method for detection and monitoring of electrical activity and motion in the presence of a magnetic field |
WO2007140861A2 (en) * | 2006-06-02 | 2007-12-13 | Eberhard-Karls-Universität Tübingen | Medical electrode device |
WO2007140861A3 (en) * | 2006-06-02 | 2008-03-20 | Univ Eberhard Karls | Medical electrode device |
EP1872715A1 (en) * | 2006-06-30 | 2008-01-02 | General Electric Company | Method and system for multi-channel biosignal processing |
WO2013028680A3 (en) * | 2011-08-25 | 2013-06-06 | Medtronic, Inc. | Method and apparatus for detecting a biomarker in the presence of electrical stimulation |
US9888861B2 (en) | 2011-08-25 | 2018-02-13 | Medtronic, Inc. | Method and apparatus for detecting a biomarker in the presence of electrical stimulation |
CN104095632A (en) * | 2013-04-07 | 2014-10-15 | 常州博睿康科技有限公司 | Method for processing electroencephalogram noise under nuclear magnetism |
Also Published As
Publication number | Publication date |
---|---|
US20040097802A1 (en) | 2004-05-20 |
EP1355571A2 (en) | 2003-10-29 |
US7286871B2 (en) | 2007-10-23 |
JP2004507293A (en) | 2004-03-11 |
AU2001283371A1 (en) | 2002-02-25 |
CA2418478A1 (en) | 2002-02-21 |
WO2002013689A3 (en) | 2003-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7286871B2 (en) | Method and apparatus for reducing contamination of an electrical signal | |
Yuan et al. | Spatiotemporal dynamics of the brain at rest—exploring EEG microstates as electrophysiological signatures of BOLD resting state networks | |
Liebenthal et al. | Simultaneous ERP and fMRI of the auditory cortex in a passive oddball paradigm | |
Srivastava et al. | ICA-based procedures for removing ballistocardiogram artifacts from EEG data acquired in the MRI scanner | |
Salek-Haddadi et al. | Studying spontaneous EEG activity with fMRI | |
Anami et al. | Stepping stone sampling for retrieving artifact-free electroencephalogram during functional magnetic resonance imaging | |
Laufs et al. | EEG-correlated fMRI of human alpha activity | |
Liu et al. | Finding thalamic BOLD correlates to posterior alpha EEG | |
Rosa et al. | Estimating the transfer function from neuronal activity to BOLD using simultaneous EEG-fMRI | |
Elbert et al. | Physical aspects of the EEG in schizophrenics | |
Jacobs et al. | Fast fMRI provides high statistical power in the analysis of epileptic networks | |
Scarff et al. | Simultaneous 3-T fMRI and high-density recording of human auditory evoked potentials | |
Sabri et al. | Attentional modulation in the detection of irrelevant deviance: a simultaneous ERP/fMRI study | |
Formaggio et al. | EEG and FMRI coregistration to investigate the cortical oscillatory activities during finger movement | |
Peters et al. | On the feasibility of concurrent human TMS-EEG-fMRI measurements | |
Ellingson et al. | Ballistocardiogram artifact reduction in the simultaneous acquisition of auditory ERPS and fMRI | |
Ritter et al. | Evaluating gradient artifact correction of EEG data acquired simultaneously with fMRI | |
Gonçalves et al. | Artifact removal in co-registered EEG/fMRI by selective average subtraction | |
Becker et al. | Visual evoked potentials recovered from fMRI scan periods | |
Rusiniak et al. | A modified oddball paradigm for investigation of neural correlates of attention: a simultaneous ERP–fMRI study | |
Storti et al. | Wavelet analysis as a tool for investigating movement-related cortical oscillations in EEG-fMRI coregistration | |
Storti et al. | A multimodal imaging approach to the evaluation of post-traumatic epilepsy | |
Sclocco et al. | EEG-informed fMRI analysis during a hand grip task: estimating the relationship between EEG rhythms and the BOLD signal | |
Salek-Haddadi et al. | EEG quality during simultaneous functional MRI of interictal epileptiform discharges | |
Kotecha et al. | Time, frequency and volumetric differences of high-frequency neuromagnetic oscillation between left and right somatosensory cortices |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A2 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT TZ UA UG US UZ VN YU ZA ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A2 Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
DFPE | Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101) | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2418478 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2002518841 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 10344776 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2001283371 Country of ref document: AU |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2001962172 Country of ref document: EP |
|
REG | Reference to national code |
Ref country code: DE Ref legal event code: 8642 |
|
WWP | Wipo information: published in national office |
Ref document number: 2001962172 Country of ref document: EP |
|
WWW | Wipo information: withdrawn in national office |
Ref document number: 2001962172 Country of ref document: EP |