US20040181162A1 - Method and system for detecting seizures using electroencephalograms - Google Patents
Method and system for detecting seizures using electroencephalograms Download PDFInfo
- Publication number
- US20040181162A1 US20040181162A1 US10/810,949 US81094904A US2004181162A1 US 20040181162 A1 US20040181162 A1 US 20040181162A1 US 81094904 A US81094904 A US 81094904A US 2004181162 A1 US2004181162 A1 US 2004181162A1
- Authority
- US
- United States
- Prior art keywords
- computer readable
- seizure
- program code
- code means
- readable program
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- 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
- A61B5/369—Electroencephalography [EEG]
-
- 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
- A61B5/369—Electroencephalography [EEG]
- A61B5/372—Analysis of electroencephalograms
-
- 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S128/00—Surgery
- Y10S128/92—Computer assisted medical diagnostics
- Y10S128/925—Neural network
Definitions
- Such simultaneous recording of voltage waveforms across different pairs of electrodes is commonly referred to as a montage, a typical example of which is illustrated in FIG. 2.
- the voltage waveform across a given pair of electrodes in the montage of an EEG recording is commonly referred to as a channel.
- a seizure is typically manifested by a highly rhythmic pattern of voltage waveforms on an EEG recording.
- an onset of seizure is sometimes not readily discernable by a human reader from a montage of an EEG recording.
- a seizure may manifest itself as a random waveform pattern across a montage of an EEG recording.
- recording errors may occur in one or more channels of a montage of an EEG recording.
- an onset of seizure is not shown on an EEG recording as a rhythmic pattern of waveforms, but rather as an abnormal change from the background waveform pattern.
- a human reader of an EEG recording may need to go through hours or even days of recorded waveforms to determine the onset, duration, and type of seizures that may have occurred during that time.
- the human reader may miss an occurrence of a seizure, which is referred to as a false negative, or may mark a segment of the waveforms as a seizure event, which is referred to as a false positive.
- the present invention provides a method of detecting a seizure, generally comprising the steps of:
- the method in an embodiment according to the present invention which utilizes advanced numerical analysis techniques including matching pursuit, neural network rules and connected-object clustering is capable of improving seizure detection with a high degree of accuracy and reliability.
- FIG. 1 is a diagram illustrating a typical arrangement of electrodes positioned on the scalp of an epilepsy patient
- FIG. 2 shows a typical montage of an EEG recording with eight channels of voltage waveforms across eight pairs of electrodes
- FIG. 3 is a flow chart illustrating an embodiment of a method for detecting seizures according to the present invention
- FIG. 5 is a block diagram illustrating a computer system for processing EEG recordings to identify seizure events in an embodiment according to the present invention.
- FIG. 3 is a flow chart illustrating a method of seizure detection using advanced numerical algorithms in an embodiment according to the present invention.
- digitized waveforms of an electroencephalogram (EEG) recording indicated by block 100 is received by a computer which processes these waveforms to detect seizure events.
- a singular waveform is defined as one of a plurality of waveforms of an EEG recording measured across a given pair of electrodes as illustrated in FIG. 1.
- An EEG recording typically includes a plurality of voltage waveforms recorded simultaneously across a plurality of predetermined pairs of electrodes. A simultaneous recording of such multiple waveforms across multiple pairs of electrodes is commonly referred to as a montage.
- the waveform measured across each designated pair of electrodes for indicating activity in a given region of the brain is commonly referred to as a channel.
- the computer determines whether unprocessed EEG data is available at the time T as indicated by block 104 and causes a time delay or a sleep of a predetermined duration, for example, one second, if unprocessed EEG data is unavailable at the time T, as indicated by block 106 .
- each digitized waveform of the EEG recording is entered for processing in an iterative manner until the processing of all waveforms of the montage of the EEG recording is complete, as indicated by block 108 .
- the processing of the first waveform 10 is followed by the processing of the second waveform 12 and so on, until the last waveform 24 is processed.
- the method according to an embodiment of the present invention utilizes three algorithmic techniques, namely, matching pursuit, small neural network (NN) rules and connected-object hierarchical clustering algorithms as indicated by blocks 110 , 112 and 114 , respectively, for enhanced seizure detection.
- an EEG recording which typically contains EEG data over a period of many hours or even days is divided into relatively small segments of time to facilitate data processing.
- Typical characteristics of seizures allow for partitioning of time to increase computational efficiency while detecting seizures with a high degree of accuracy.
- a short seizure lasts for about two seconds.
- the matching pursuit algorithm has been developed as a result of the desire to improve the time and frequency resolution of traditional fast Fourier transform (FFT) algorithms. It is well known in the art that wavelets typically give high-frequency bands better time resolution with corresponding loss in frequency resolution. It is also well known in the art that wavelet packets typically allow the time-frequency resolution to fit to the signal for different frequency bands, thereby resulting in a trade-off between the frequency resolution and the time resolution within a given frequency band. For example, with conventional wavelet packet processing, each frequency band may be processed with a high frequency resolution and a low time resolution, or alternatively, with a low frequency resolution and a high time resolution.
- the matching pursuit algorithm which is an improvement over conventional wavelet packet algorithms, affords complete flexibility to different parts of signals, even within the same frequency band, in a trade-off between frequency resolution and time resolution.
- t time
- s>0 the scale
- u the translation
- v the frequency
- the Gabor dictionary i.e., all values of the index I
- NN neural network
- Small neural networks can be combined with domain expertise to create NN rules, which are similar to expert system rules that propagate the strength or uncertainty in the inputs to the NN.
- a mathematical description of such NN rules is described in Wilson et al., Spike detection: reduction of model complexity via small neural networks constrained by domain expertise, Clinical Neurophysiology , submitted Mar. 5, 2002.
- Tendency rules are used to encode weak expert knowledge and describe the general shape of the NN's transfer function, i.e., the graph of the NN inputs vs. output.
- Inputs are displayed on the x and y axes, and the output is displayed on the z axis. If the NN has more than two inputs, then the other inputs are held constant at their average value, and pairs of inputs are examined.
- the tendency rule SzRhythm/BkRhythm and SzAmp/BkAmp->Perception encompasses the general notion that seizures should be more rhythmic and probably at higher voltages compared to the background activity. Any. NN-rule which has a transfer function that contradicts this knowledge is not accepted.
- Rules are traditionally written in the form if x and Y then Z.
- Boolean logic describes how to combine these types of rules if X, Y and Z are variables taking the possible values of true and false.
- Expert and fuzzy systems are developed from rules with the understanding that many interesting problems can best be expressed using non-Boolean variables. They extend Boolean logic to allow rules that propagate the uncertainty and/or strength of the premise to the rule's conclusion. If the rule has only a single conclusion, then the general form of these rules is a function that takes one or more inputs and has a single output, which is the essence of a neural network. If the NN is sufficiently simple, it may be possible to describe it linguistically as one or more rules.
- the rules that are developed for detection of seizures in one embodiment are very small NNs consisting of two to five inputs and a single output. Expert knowledge is used to choose the inputs and the related output. After the NN is trained, it is only accepted if the response surface, i.e., what the rule says, is consistent with the expert knowledge. The requirement that the NN be consistent with expert knowledge, rather than simply a “black box,” is important to this embodiment. It is preferred that small neural networks that can be analyzed via their response surface be implemented to allow for penetration of the conventional “black box” shell. In contrast, the response surface of a poorly defined network which displays discontinuities or contradicts domain knowledge is undesirable in these applications.
- an EEG recording includes EEG data recorded over a relatively long period of time, from several hours to several days, for example.
- the EEG data covering the total time period of many hours or days usually need not be analyzed at once.
- connected-object clustering algorithm is applied to a plurality of epochs in a sliding window of a pre-determined duration, for example, 60 seconds, as indicated by block 114 in the flow chart of FIG. 3.
- EEG data that are older than sixty seconds are discarded.
- One of the complicating aspects of seizure detection is the need to perform two tasks at roughly the same time: to identify the event and to determine its approximate endpoints.
- a novel type of clustering is implemented in an embodiment to identify the onset and approximate duration of an event of brain activity, and then to decide whether the event is a seizure event.
- Clustering methods fall into two general categories: hierarchical and partitional (Jain, A. K. and Dubes, R. C., Algorithms for clustering data , Englewood Cliffs, N.J., Prentice Hall, 1988). Partitional methods require as input the appropriate number of groups, and the groups found are dependent on the order that the objects are submitted.
- the computational complexity i.e., the approximate number of operations
- Hierarchical methods generate a nested sequence that allows the appropriate number of groups to be identified interactively, using a dendrogram, for example, or via a heuristic method.
- Step 2 Choose index k to maximize the difference between the two groups. Stop if the desired number of groups is found. Stop this branch of the algorithm if the difference between the two groups is too small.
- Step 2 requires a function that rates the difference between two groups.
- the Davies-Bouldin index (Jain, A. K. and Dubes, R. C., Algorithms for clustering data , Englewood Cliffs, N.J., Prentice Hall, 1988, p. 186) is used, which was developed as a way to determine the “correct” number of clusters.
- the nested structure of the clustering is captured by storing the indexes k in a tree at each step.
- DB(K) can also be used to choose overall partitioning of the tree, because for strongly clustered data it will be minimized at the “correct” number of clusters.
- Klatchko et al. (Klatchko, A., Raviv, G., Webber, W. R., and Lesser, R. P., Enhancing the detection of seizures with a clustering algorithm , Electroencephalogr. Clin. Neurophysiol 106: 52-63, 1998) apply clustering (via directed graph) to seizure detection to enhance an existing detector.
- Candidate seizure epochs are grouped together if they independently have a probability that exceeds a cutoff value, are in the same or neighboring channels, and occur simultaneously or sequentially. This method does not work when there are gaps, i.e., non-seizure epochs, in the sequence of candidate epochs. Because the grouping condition only considers proximal pairs of epochs, the same result can be achieved without clustering.
- FIG. 4 is a flow chart illustrating an example of a connected-object clustering algorithm with two branches in an embodiment according to the present invention.
- N which is the total number of groups in the nested structure, is initially set to one, as indicated by block 202 .
- k is selected such that the difference Diff[k] between the two groups of k and N ⁇ k objects is maximized, and M is incremented by 1, as indicated by block 214 in FIG. 4.
- a seizure point (SzPt) is established from the clustering result and described with neural network rules, as indicated by block 116 .
- a determination is then made as to whether the seizure point is valid and better than a pre-existing seizure point, if any, as indicated by block 118 . If the seizure point derived from the current waveform is better than a pre-existing seizure point established from a previous waveform, the current seizure point is then regarded as the “best” seizure point, as indicated by block 120 .
- the pre-existing seizure point remains as the “best” seizure point.
- the matching pursuit, neural network rules and connected-object clustering algorithms as indicated by blocks 110 , 112 , and 114 are applied to other digitized waveforms in the montage of the EEG recording until all of the waveforms within the sixty-second sliding window for all channels of the montage are processed.
- one or more seizure events are established by analyzing a plurality of proximal seizure points from the waveforms of the montage, as indicated by block 122 in FIG. 3. After a seizure event is established with the seizure points that are temporally close to each other from the waveforms of the montage, the algorithm determines whether such a seizure event (SzEvt) is valid, as indicated by block 124 in FIG. 3.
- the sliding window is shifted by a small increment of time, for example, one second, until the whole EEG recording is processed, as indicated by blocks 130 and 132 . If the EEG recording has been completely processed, the program ends as indicated by block 134 . If the processing of the EEG recording has not been completed, one second is added to the time T and the sliding window is shifted by one second as indicated by block 132 . The program then loops back to block 104 to process new EEG data within the newly shifted sixty-second window.
- FIG. 5 is a simplified block diagram illustrating a computer system or an article of manufacture for processing EEG data according to embodiments of the present invention.
- the computer system includes a central processing unit (CPU) 300 , a computer readable medium such as a memory 302 , an input device such as a keyboard 304 , and output devices such as a display monitor 306 and a printer 308 .
- the digitized EEG recording 310 is read into the memory 302 for processing by the CPU 300 , which is capable of running a computer program code comprising instructions for performing process steps that include matching pursuit, neural network rules, and connected-object clustering algorithms according to embodiments of the present invention.
Abstract
A method and a system are provided for detection of seizures by applying advanced numerical analysis techniques to digitized waveforms of an electroencephalogram (EEG) recording. In an embodiment, the advanced numerical analysis techniques implemented in the method and system for detecting seizures include a combination of matching pursuit, neural network rules, and connected-object clustering algorithms.
Description
- 1. Field of the Invention
- The present invention relates to detection of seizures, and more particularly, to detection of seizures with advanced algorithms.
- 2. Background Art
- Epilepsy is a condition characterized by recurrent seizures which are the outward manifestation of excessive or hypersynchronous abnormal electrical activity of neurons in the cerebral cortex of the brain. A seizure patient may suffer from several different types of seizures, or any combination thereof. For example, a common type of epilepsy is called the grand mal seizure, which is manifested by symptoms of convulsions with tonic-clonic contractions of muscles. Another type of epilepsy is called the absence seizure, which is characterized by brief and sudden loss of consciousness. Other types of seizures include complex partial seizure, which is characterized by a complete loss of consciousness, and psychomotor seizure, which is characterized by clouding of consciousness for one to two minutes. Some types of seizures may involve the entire brain, while other types of seizures may affect only a local portion of the brain.
- Electroencephalograms (EEGs) have been employed to record electrical signals generated by different parts of the brain. In a typical EEG, a plurality of electrodes are placed across the scalp of a seizure patient with predetermined spacings. FIG. 1 shows a diagram illustrating a typical arrangement of electrodes positioned on the scalp of an epilepsy patient along standard lines of measurements. Voltages measured across predetermined pairs of electrodes, for example, C3-F3, F3-F7, F7-T3, T3-T5, etc., are measured simultaneously and recorded as waveforms over a long period of time. Such simultaneous recording of voltage waveforms across different pairs of electrodes is commonly referred to as a montage, a typical example of which is illustrated in FIG. 2. The voltage waveform across a given pair of electrodes in the montage of an EEG recording is commonly referred to as a channel. A seizure is typically manifested by a highly rhythmic pattern of voltage waveforms on an EEG recording. However, depending upon the individual patient, different types of seizures, and various other factors, an onset of seizure is sometimes not readily discernable by a human reader from a montage of an EEG recording. For example, sometimes a seizure may manifest itself as a random waveform pattern across a montage of an EEG recording. Sometimes recording errors may occur in one or more channels of a montage of an EEG recording. Sometimes an onset of seizure is not shown on an EEG recording as a rhythmic pattern of waveforms, but rather as an abnormal change from the background waveform pattern.
- A human reader of an EEG recording may need to go through hours or even days of recorded waveforms to determine the onset, duration, and type of seizures that may have occurred during that time. The human reader may miss an occurrence of a seizure, which is referred to as a false negative, or may mark a segment of the waveforms as a seizure event, which is referred to as a false positive.
- Conventional algorithms have been developed to assist a human reader in detecting seizures using traditional fast Fourier transform (FFT) or other spectral analysis techniques. While these traditional techniques are usually effective in detecting highly rhythmic patterns of waveforms in order to identify a seizure event, some types of seizures which are not manifested by such highly rhythmic patterns may still be missed. Signal processing using conventional spectral analysis techniques may also sometimes return false positives, depending upon the perimeters set by the algorithm.
- Therefore, there is a need for an improved system and method for detecting seizures from EEG recordings with a high degree of reliability.
- The present invention provides a method of detecting a seizure, generally comprising the steps of:
- (a) dividing a digitized waveform of an electroencephalogram (EEG) recording into a plurality of epochs each having a first predetermined duration;
- (b) computing matching pursuit for a given one of the epochs to obtain a plurality of seizure atoms;
- (c) describing the seizure atoms and the given epoch with at least one neural network (NN) rule;
- (d) applying connected-object clustering to the epochs in a sliding window of a second predetermined duration to obtain a clustering result; and
- (e) establishing a seizure point from the clustering result.
- Advantageously, the method in an embodiment according to the present invention which utilizes advanced numerical analysis techniques including matching pursuit, neural network rules and connected-object clustering is capable of improving seizure detection with a high degree of accuracy and reliability.
- The present invention will be described with particular embodiments thereof, and references will be made to the drawings in which:
- FIG. 1 is a diagram illustrating a typical arrangement of electrodes positioned on the scalp of an epilepsy patient;
- FIG. 2 shows a typical montage of an EEG recording with eight channels of voltage waveforms across eight pairs of electrodes;
- FIG. 3 is a flow chart illustrating an embodiment of a method for detecting seizures according to the present invention;
- FIG. 4 is a flow chart illustrating an embodiment of connected-object clustering for establishing a seizure point according to the present invention; and
- FIG. 5 is a block diagram illustrating a computer system for processing EEG recordings to identify seizure events in an embodiment according to the present invention.
- FIG. 3 is a flow chart illustrating a method of seizure detection using advanced numerical algorithms in an embodiment according to the present invention. Initially, digitized waveforms of an electroencephalogram (EEG) recording indicated by block100 is received by a computer which processes these waveforms to detect seizure events. In this patent application, a singular waveform is defined as one of a plurality of waveforms of an EEG recording measured across a given pair of electrodes as illustrated in FIG. 1. An EEG recording typically includes a plurality of voltage waveforms recorded simultaneously across a plurality of predetermined pairs of electrodes. A simultaneous recording of such multiple waveforms across multiple pairs of electrodes is commonly referred to as a montage. The waveform measured across each designated pair of electrodes for indicating activity in a given region of the brain is commonly referred to as a channel.
- FIG. 2 illustrates a typical montage of an EEG recording with eight channels of simultaneous waveforms recorded over a given period of time. For example, a
first waveform 10 in FIG. 2 is a voltage waveform measured across a pair of electrodes C3 and F3 in FIG. 1, while asecond waveform 12 in FIG. 2 is a voltage waveform measured across another pair of electrodes F3 and F7 in FIG. 1. Theremaining waveforms - Referring to FIG. 3, after the digitized waveforms' of the EEG recording is received by the computer, the algorithm is initialized and the start of the EEG recording is set at a time T=0 as indicated by
block 102. The computer then determines whether unprocessed EEG data is available at the time T as indicated byblock 104 and causes a time delay or a sleep of a predetermined duration, for example, one second, if unprocessed EEG data is unavailable at the time T, as indicated byblock 106. - If unprocessed EEG data is available at the time T, each digitized waveform of the EEG recording is entered for processing in an iterative manner until the processing of all waveforms of the montage of the EEG recording is complete, as indicated by
block 108. For example, referring to FIG. 2, the processing of thefirst waveform 10 is followed by the processing of thesecond waveform 12 and so on, until thelast waveform 24 is processed. For each waveform of the digitized EEG recording, the method according to an embodiment of the present invention utilizes three algorithmic techniques, namely, matching pursuit, small neural network (NN) rules and connected-object hierarchical clustering algorithms as indicated byblocks - Before each digitized EEG waveform is processed with matching pursuit, small neural network rules and connected-object clustering algorithms, an EEG recording which typically contains EEG data over a period of many hours or even days is divided into relatively small segments of time to facilitate data processing. Typical characteristics of seizures allow for partitioning of time to increase computational efficiency while detecting seizures with a high degree of accuracy. A short seizure lasts for about two seconds.
- In an embodiment, a digitized EEG waveform is divided into two-second intervals, each interval commonly referred to as an epoch. To enhance the probability of detection of short-duration seizures, the processing of EEG data is applied to each two-second epoch with a time shift of one second. For example, referring to FIG. 2, matching pursuit is computed for a given waveform for a first epoch of two seconds between times T=0 and T=2, followed by a second epoch between times T=1 and T=3, followed by a third epoch between times T=2 and T=4, and so on. Such an overlap of one second for each successive two-second epoch allows for a greater probability of detection of short-duration seizures.
- The matching pursuit algorithm that is applied to each two-second epoch as indicated by
block 110 in the flow chart of FIG. 3 is described mathematically in Wickerhauser, Adapted Wavelet Analysis From Theory to Software, A. K. Peters, Wellesley, Mass., 1994, and Mallat, A Wavelet Tour of Signal Processing, Academic Press, San Diego, 1999, both of which are incorporated herein by reference. - The matching pursuit algorithm has been developed as a result of the desire to improve the time and frequency resolution of traditional fast Fourier transform (FFT) algorithms. It is well known in the art that wavelets typically give high-frequency bands better time resolution with corresponding loss in frequency resolution. It is also well known in the art that wavelet packets typically allow the time-frequency resolution to fit to the signal for different frequency bands, thereby resulting in a trade-off between the frequency resolution and the time resolution within a given frequency band. For example, with conventional wavelet packet processing, each frequency band may be processed with a high frequency resolution and a low time resolution, or alternatively, with a low frequency resolution and a high time resolution. The matching pursuit algorithm, which is an improvement over conventional wavelet packet algorithms, affords complete flexibility to different parts of signals, even within the same frequency band, in a trade-off between frequency resolution and time resolution.
-
- where t is time, s>0 is the scale, u is the translation and v is the frequency. The index I=(s, v, u) describes the chosen set of parameters. If g(t)=e−πt
2 , then gI(t) is called a Gabor atom, which is a sinusoid of frequency v with an exponentially localized time window centered on u. From the set of all possible atoms, called the Gabor dictionary (i.e., all values of the index I), the best single atom is found in each step and then subtracted from the waveform. This continues for a fixed number of iterations or until the error of the residual waveform reaches an acceptable size. - The use of Gabor atoms as described above offers precise time-frequency resolution, but it is relatively slow—about ˜22 ms per atom for each two-second epoch of 32 channels sampled at 32 Hz on a 866 MHz processor coded in C++. A faster version of matching pursuit is achieved by using a wavelet packet dictionary in place of the Gabor dictionary. (Like the Gabor dictionary, the wavelet packet is a redundant division of the time-frequency plane into various sized rectangles.) This results in a speed improvement of about eighteen times with some loss in resolution, which is usually acceptable in seizure detection applications.
- After the matching pursuit algorithm is applied to each two-second epoch (SzEpoch) by decomposing the epochs of the EEG waveform into atoms (SzAtoms), the seizure atoms and the seizure epoch are described with one or more neural network (NN) rules, as indicated by
block 112 in the flow chart of FIG. 3. Small neural networks can be combined with domain expertise to create NN rules, which are similar to expert system rules that propagate the strength or uncertainty in the inputs to the NN. A mathematical description of such NN rules is described in Wilson et al., Spike detection: reduction of model complexity via small neural networks constrained by domain expertise, Clinical Neurophysiology, submitted Mar. 5, 2002. Using this methodology, the NNs are not over-trained and the number of required training samples can be reduced by reducing the model complexity. Tendency rules are used to encode weak expert knowledge and describe the general shape of the NN's transfer function, i.e., the graph of the NN inputs vs. output. - Inputs are displayed on the x and y axes, and the output is displayed on the z axis. If the NN has more than two inputs, then the other inputs are held constant at their average value, and pairs of inputs are examined. For example, the tendency rule SzRhythm/BkRhythm and SzAmp/BkAmp->Perception encompasses the general notion that seizures should be more rhythmic and probably at higher voltages compared to the background activity. Any. NN-rule which has a transfer function that contradicts this knowledge is not accepted.
- Rules are traditionally written in the form if x and Y then Z. Boolean logic describes how to combine these types of rules if X, Y and Z are variables taking the possible values of true and false. Expert and fuzzy systems are developed from rules with the understanding that many interesting problems can best be expressed using non-Boolean variables. They extend Boolean logic to allow rules that propagate the uncertainty and/or strength of the premise to the rule's conclusion. If the rule has only a single conclusion, then the general form of these rules is a function that takes one or more inputs and has a single output, which is the essence of a neural network. If the NN is sufficiently simple, it may be possible to describe it linguistically as one or more rules.
- The rules that are developed for detection of seizures in one embodiment are very small NNs consisting of two to five inputs and a single output. Expert knowledge is used to choose the inputs and the related output. After the NN is trained, it is only accepted if the response surface, i.e., what the rule says, is consistent with the expert knowledge. The requirement that the NN be consistent with expert knowledge, rather than simply a “black box,” is important to this embodiment. It is preferred that small neural networks that can be analyzed via their response surface be implemented to allow for penetration of the conventional “black box” shell. In contrast, the response surface of a poorly defined network which displays discontinuities or contradicts domain knowledge is undesirable in these applications.
- It is typical that an EEG recording includes EEG data recorded over a relatively long period of time, from several hours to several days, for example. The EEG data covering the total time period of many hours or days usually need not be analyzed at once. In an embodiment, connected-object clustering algorithm is applied to a plurality of epochs in a sliding window of a pre-determined duration, for example, 60 seconds, as indicated by block114 in the flow chart of FIG. 3.
- As illustrated in FIG. 2, a sixty-second sliding window includes thirty epochs of two seconds each and initially starts with T=0. The initial sixty-second sliding window covers a time period between T=0 and T=60. After the EEG data for this sixty-second window is processed, the sliding window shifts to the right by one second covering a time period between T=1 and T=61. As the sixty-second sliding window shifts progressively in one-second intervals, EEG data that are older than sixty seconds are discarded. By applying the sliding window to the EEG waveform, connected-object clustering can be performed on the epochs to obtain a clustering result with improved computational speed and efficiency.
- One of the complicating aspects of seizure detection is the need to perform two tasks at roughly the same time: to identify the event and to determine its approximate endpoints. A novel type of clustering is implemented in an embodiment to identify the onset and approximate duration of an event of brain activity, and then to decide whether the event is a seizure event.
- Clustering methods fall into two general categories: hierarchical and partitional (Jain, A. K. and Dubes, R. C.,Algorithms for clustering data, Englewood Cliffs, N.J., Prentice Hall, 1988). Partitional methods require as input the appropriate number of groups, and the groups found are dependent on the order that the objects are submitted. The computational complexity (i.e., the approximate number of operations) for these algorithms is O(N*K*R), where N is the number of objects, K is the number of groups, and R is the number of iterations. Hierarchical methods generate a nested sequence that allows the appropriate number of groups to be identified interactively, using a dendrogram, for example, or via a heuristic method. The high computational complexity [O(N*N*N)] of hierarchical algorithms makes them undesirable for seizure detection applications. For example, an eight-hour record includes N=28,800 one-second epochs. These requirements can be reduced to O(N*log(M)), where M is the maximum number of groups in the nested structure, by taking advantage of the connected nature of the epochs. The following are the steps of an embodiment of the connected-object clustering algorithm:
-
Step 1. Given an ordered list of N connected objects, break the list into two groups of k and N−k objects for k=1 . . . N−1. If N=1, then stop this branch of the algorithm. -
Step 2. Choose index k to maximize the difference between the two groups. Stop if the desired number of groups is found. Stop this branch of the algorithm if the difference between the two groups is too small. - Step 3. Apply the algorithm to the two resulting groups.
-
Step 2 requires a function that rates the difference between two groups. The Davies-Bouldin index (Jain, A. K. and Dubes, R. C., Algorithms for clustering data, Englewood Cliffs, N.J., Prentice Hall, 1988, p. 186) is used, which was developed as a way to determine the “correct” number of clusters. For K clusters, the index is given by - for K>1.
-
- where nj, μj and σj are the number of objects, the mean and the standard deviation of cluster j. The splitting index k is chosen to maximize DB(K).
- The nested structure of the clustering is captured by storing the indexes k in a tree at each step. DB(K) can also be used to choose overall partitioning of the tree, because for strongly clustered data it will be minimized at the “correct” number of clusters.
- A desirable feature of this algorithm, not easily supported by other clustering methods, is that new objects (e.g., new EEG epochs as they are acquired) can be added to the end of the list without having to start over. Komzak (Komzak, Jiri,Hierarchical clustering speed up using position lists and data position hierarchy, www.kmi.open.ac.uk/papers/kmi-tr-115.pdf, 2001) has developed a more general, but slower, O(N*N), method for speeding up clustering by using sorted objects.
- Klatchko et al. (Klatchko, A., Raviv, G., Webber, W. R., and Lesser, R. P.,Enhancing the detection of seizures with a clustering algorithm, Electroencephalogr. Clin. Neurophysiol 106: 52-63, 1998) apply clustering (via directed graph) to seizure detection to enhance an existing detector. Candidate seizure epochs are grouped together if they independently have a probability that exceeds a cutoff value, are in the same or neighboring channels, and occur simultaneously or sequentially. This method does not work when there are gaps, i.e., non-seizure epochs, in the sequence of candidate epochs. Because the grouping condition only considers proximal pairs of epochs, the same result can be achieved without clustering.
- FIG. 4 is a flow chart illustrating an example of a connected-object clustering algorithm with two branches in an embodiment according to the present invention. Initially, an ordered list of N connected objects is supplied to the algorithm as indicated by
block 200. M, which is the total number of groups in the nested structure, is initially set to one, as indicated by block 202. The algorithm then determines whether N=1 as indicated byblock 204, and the current branch of the algorithm is stopped if N=1, as indicated byblock 206. If N is greater than 1, the algorithm enters an iterative loop as indicated byblocks - Subsequently, k is selected such that the difference Diff[k] between the two groups of k and N−k objects is maximized, and M is incremented by 1, as indicated by
block 214 in FIG. 4. The connected-object clustering algorithm determines whether M is equal to a desired count as indicated byblock 216, and stops all branches of the algorithm as indicated byblock 218 if M is equal to the desired count. If M is not equal to the desired count, the algorithm is branched via recursion as indicated by block 220. One branch of the algorithm sets N=k as indicated byblock 222 while another branch sets N=N−k as indicated byblock 224. These two branches of the algorithm are then looped back to block 204 as illustrated in FIG. 4. - Referring back to FIG. 3, after the connected-object clustering algorithm is applied to the epochs in the sixty-second sliding window as indicated by block114, a seizure point (SzPt) is established from the clustering result and described with neural network rules, as indicated by block 116. A determination is then made as to whether the seizure point is valid and better than a pre-existing seizure point, if any, as indicated by
block 118. If the seizure point derived from the current waveform is better than a pre-existing seizure point established from a previous waveform, the current seizure point is then regarded as the “best” seizure point, as indicated byblock 120. Otherwise, the pre-existing seizure point remains as the “best” seizure point. Subsequently, the matching pursuit, neural network rules and connected-object clustering algorithms as indicated byblocks - After all of the waveforms in the sixty-second sliding window of the montage are processed, one or more seizure events (SzEvts) are established by analyzing a plurality of proximal seizure points from the waveforms of the montage, as indicated by
block 122 in FIG. 3. After a seizure event is established with the seizure points that are temporally close to each other from the waveforms of the montage, the algorithm determines whether such a seizure event (SzEvt) is valid, as indicated byblock 124 in FIG. 3. - If the seizure event is valid, the computer notifies the user that a seizure has been identified, as indicated by block126 in FIG. 3. The time of onset of the seizure and a description of the seizure based upon the EEG data, including, for example, the duration, type, and severity of the seizure, may be saved in the computer in an embodiment as indicated by
block 128 in FIG. 3. - After all waveforms of the montage within a given sixty-second sliding window is processed, the sliding window is shifted by a small increment of time, for example, one second, until the whole EEG recording is processed, as indicated by
blocks 130 and 132. If the EEG recording has been completely processed, the program ends as indicated by block 134. If the processing of the EEG recording has not been completed, one second is added to the time T and the sliding window is shifted by one second as indicated by block 132. The program then loops back to block 104 to process new EEG data within the newly shifted sixty-second window. - FIG. 5 is a simplified block diagram illustrating a computer system or an article of manufacture for processing EEG data according to embodiments of the present invention. As illustrated in FIG. 5, the computer system includes a central processing unit (CPU)300, a computer readable medium such as a
memory 302, an input device such as akeyboard 304, and output devices such as adisplay monitor 306 and aprinter 308. Thedigitized EEG recording 310 is read into thememory 302 for processing by theCPU 300, which is capable of running a computer program code comprising instructions for performing process steps that include matching pursuit, neural network rules, and connected-object clustering algorithms according to embodiments of the present invention. - From the above description of embodiments of the invention, it is manifest that various equivalents can be used to implement the concepts of the present invention without departing from its scope. Moreover, while the invention has been described with specific reference to certain embodiments, a person of ordinary skill in the art would recognize that changes can be made in form and detail without departing from the spirit and the scope of the invention. The described embodiments are to be considered in all respects as illustrative and not restrictive. It should also be understood that the invention is not limited to the particular embodiments described herein, but is capable of many equivalents, rearrangements, modifications, and substitutions without departing from the scope of the invention as set forth in the claims.
Claims (45)
1. A method of detecting a seizure, comprising the steps of:
a) dividing a digitized waveform of an electroencephalogram (EEG) recording into a plurality of epochs each having a first predetermined duration;
(b) computing matching pursuit for a given one of the epochs to obtain a plurality of seizure atoms;
(c) describing the seizure atoms and the given epoch with at least one neural network (NN) rule;
(d) applying connected-object clustering to the epochs in a sliding window of a second predetermined duration to obtain a clustering result; and
(e) establishing a seizure point from the clustering result.
2. The method of claim 1 , further comprising the step of determining whether the seizure point is valid.
3. The method of claim 1 , further comprising the steps of repeating steps (b)-(e) and determining whether a succeeding seizure point is better than a preceding seizure point.
4. The method of claim 1 , wherein the first predetermined duration is less than the second predetermined duration.
5. The method of claim 4 , wherein the first predetermined duration is about 2 seconds.
6. The method of claim 4 , wherein the second predetermined duration is about 60 seconds.
7. The method of claim 1 , further comprising the step of setting an initial time T=0 prior to the step of computing matching pursuit for the given epoch.
8. The method of claim 7 , further comprising the steps of determining whether unprocessed EEG data is available at the time T, and causing a time delay of a third predetermined duration if unprocessed EEG data is unavailable at the time T.
9. The method of claim 8 , wherein the digital waveform is one of a plurality of waveforms forming a montage of the EEG recording, the method further comprising the steps of obtaining a subsequent waveform of the EEG recording and repeating steps (b)-(e) if unprocessed EEG data is available at the time T.
10. The method of claim 9 , further comprising the step of establishing a seizure event with a plurality of proximal seizure points from the waveforms of the montage.
11. The method of claim 10 , further comprising the step of determining whether the seizure event is valid.
12. The method of claim 11 , further comprising the step of notifying a user that a seizure is identified if the seizure event is valid.
13. The method of claim 12 , further comprising the step of saving the time T and a description of the seizure.
14. The method of claim 11 , further comprising the steps of adding a time of the third predetermined duration to the time T and shifting the sliding window by the time of the third predetermined duration until the EEG recording is finished.
15. The method of claim 14 , wherein the third predetermined duration is about 1 second.
16. A computer readable medium comprising a plurality of instructions, which when executed by a computer, cause the computer to perform the steps of:
(a) dividing a digitized waveform of an electroencephalogram (EEG) recording into a plurality of epochs each having a first predetermined duration;
(b) computing matching pursuit for a given one of the epochs to obtain a plurality of seizure atoms;
(c) describing the seizure atoms and the given epoch with at least one neural network (NN) rule;
(d) applying connected-object clustering to the epochs in a sliding window of a second predetermined duration to obtain a clustering result; and
(e) establishing a seizure point from the clustering result.
17. The computer readable medium as set forth in claim 16 , further comprising instructions which cause the computer to perform the step of determining whether the seizure point is valid.
18. The computer readable medium as set forth in claim 16 , further comprising instructions which cause the computer to perform the steps of:
repeating steps (b)-(e); and
determining whether a succeeding seizure point is better than a preceding seizure point.
19. The computer readable medium as set forth in claim 16 , wherein the first predetermined duration is less than the second predetermined duration.
20. The computer readable medium as set forth in claim 19 , wherein the first predetermined duration is about 2 seconds.
21. The computer readable medium as set forth in claim 19 , wherein the second predetermined duration is about 60 seconds.
22. The computer readable medium as set forth in claim 16 , further comprising instructions which cause the computer to perform the step of setting an initial time T=0 prior to the step of computing matching pursuit for the given epoch.
23. The computer readable medium as set forth in claim 22 , further comprising instructions which cause the computer to perform the steps of:
determining whether unprocessed EEG data is available at the time T; and causing a time delay of a third predetermined duration if unprocessed EEG data is unavailable at the time T.
24. The computer readable medium as set forth in claim 23 , wherein the digitized waveform is one of a plurality of waveforms of the EEG recording which includes a montage of multiple channels of simultaneous waveforms, the computer readable medium further comprising instructions which cause the computer to perform the steps of:
obtaining a subsequent waveform of the EEG recording; and
repeating steps (b)-(e) if unprocessed EEG data is available at the time T.
25. The computer readable medium as set forth in claim 24 , further comprising instructions which cause the computer to perform the step of establishing a seizure event with a plurality of proximal seizure points from the waveforms of the montage.
26. The computer readable medium as set forth in claim 25 , further comprising instructions which cause the computer to perform the step of determining whether the seizure event is valid.
27. The computer readable medium as set forth in claim 26 , further comprising instructions which cause the computer to perform the step of notifying a user that a seizure is identified if the seizure event is valid.
28. The computer readable medium as set forth in claim 27 , further comprising instructions which cause the computer to perform the step of saving the time T and a description of the seizure.
29. The computer readable medium as set forth in claim 26 , further comprising instructions which cause the computer to perform the steps of:
adding a time of the third predetermined duration to the time T; and
shifting the sliding window by the time of the third predetermined duration until the EEG recording is finished.
30. The computer readable medium as set forth in claim 29 , wherein the third predetermined duration is about 1 second.
31. An article of manufacture, including a computer readable medium having computer readable program code means embodied therein for detecting a seizure, the computer readable program code means in the article of manufacture comprising:
(a) computer readable program code means for dividing a digitized waveform of an electroencephalogram (EEG) recording into a plurality of epochs each having a first predetermined duration;
(b) computer readable program code means for computing matching pursuit for a given one of the epochs to obtain a plurality of seizure atoms;
(c) computer readable program code means for describing the seizure atoms and the given epoch with at least one neural network (NN) rule;
(d) computer readable program code means for applying connected-object clustering to the epochs in a sliding window of a second predetermined duration to obtain a clustering result; and
(e) computer readable program code means for establishing a seizure point from the clustering result.
32. The article of manufacture as set forth in claim 31 , wherein the computer readable program code means further comprises computer readable program code means for determining whether the seizure point is valid.
33. The article of manufacture as set forth in claim 31 , wherein the computer readable program code means further comprises:
computer readable program code means for repeating steps performed by computer readable program code means (b)-(e); and
computer readable program code means for determining whether a succeeding seizure point is better than a preceding seizure point.
34. The article of manufacture as set forth in claim 31 , wherein the first predetermined duration is less than the second predetermined duration.
35. The article of manufacture as set forth in claim 34 , wherein the first predetermined duration is about 2 seconds.
36. The article of manufacture as set forth in claim 34 , wherein the second predetermined duration is about 60 seconds.
37. The article of manufacture as set forth in claim 31 , wherein the computer readable program code means further comprises computer readable program code means for setting an initial time T=0 prior to computing matching pursuit for the given epoch.
38. The article of manufacture as set forth in claim 37 , wherein the computer readable program code means further comprises:
computer readable program code means for determining whether unprocessed EEG data is available at the time T; and
computer readable program code means for causing a time delay of a third predetermined duration if unprocessed EEG data is unavailable at the time T.
39. The article of manufacture as set forth in claim 38 , wherein the digitized waveform is one of a plurality of waveforms of the EEG recording which includes a montage of multiple channels of simultaneous waveforms, and wherein the computer readable program code means further comprises:
computer readable program code means for obtaining a subsequent waveform of the EEG recording; and
computer readable program code means for repeating steps performed by computer readable program code means (b)-(e) if unprocessed EEG data is available at the time T.
40. The article of manufacture as set forth in claim 39 , wherein the computer readable program code means further comprises computer readable program code means for establishing a seizure event with a plurality of proximal seizure points from the waveforms of the montage.
41. The article of manufacture as set forth in claim 40 , wherein the computer readable program code means further comprises computer readable program code means for determining whether the seizure event is valid.
42. The article of manufacture as set forth in claim 41 , wherein the computer readable program code means further comprises computer readable program code means for notifying a user that a seizure is identified if the seizure event is valid.
43. The article of manufacture as set forth in claim 42 , wherein the computer readable program code means further comprises computer readable program code means for saving the time T and a description of the seizure.
44. The article of manufacture as set forth in claim 41 , wherein the computer readable program code means further comprises:
computer readable program code means for adding the time of the third predetermined duration to the time T; and
computer readable program code means for shifting the sliding window by the time of the third predetermined duration until the EEG recording is finished.
45. The article of manufacture as set forth in claim 44 , wherein the third predetermined duration is about 1 second.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/810,949 US20040181162A1 (en) | 2002-04-15 | 2004-03-25 | Method and system for detecting seizures using electroencephalograms |
US10/830,177 US20040199482A1 (en) | 2002-04-15 | 2004-04-21 | Systems and methods for automatic and incremental learning of patient states from biomedical signals |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/123,834 US6735467B2 (en) | 2002-04-15 | 2002-04-15 | Method and system for detecting seizures using electroencephalograms |
US10/810,949 US20040181162A1 (en) | 2002-04-15 | 2004-03-25 | Method and system for detecting seizures using electroencephalograms |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/123,834 Continuation US6735467B2 (en) | 2002-04-15 | 2002-04-15 | Method and system for detecting seizures using electroencephalograms |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/830,177 Continuation-In-Part US20040199482A1 (en) | 2002-04-15 | 2004-04-21 | Systems and methods for automatic and incremental learning of patient states from biomedical signals |
Publications (1)
Publication Number | Publication Date |
---|---|
US20040181162A1 true US20040181162A1 (en) | 2004-09-16 |
Family
ID=28790828
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/123,834 Expired - Fee Related US6735467B2 (en) | 2002-04-15 | 2002-04-15 | Method and system for detecting seizures using electroencephalograms |
US10/810,949 Abandoned US20040181162A1 (en) | 2002-04-15 | 2004-03-25 | Method and system for detecting seizures using electroencephalograms |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/123,834 Expired - Fee Related US6735467B2 (en) | 2002-04-15 | 2002-04-15 | Method and system for detecting seizures using electroencephalograms |
Country Status (1)
Country | Link |
---|---|
US (2) | US6735467B2 (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050080828A1 (en) * | 2003-10-08 | 2005-04-14 | Johnson Mark W. | Method for displaying spectral trends in complex signals |
WO2006122201A3 (en) * | 2005-05-10 | 2007-12-21 | Salk Inst For Biological Studi | Automated detection of sleep and waking states |
US20080102236A1 (en) * | 2006-10-27 | 2008-05-01 | Fish Robert B | Pipes containing nanoclays and method for their manufacture |
CN106821376A (en) * | 2017-03-28 | 2017-06-13 | 南京医科大学 | A kind of epileptic attack early warning system and method based on deep learning algorithm |
US10803145B2 (en) | 2016-02-05 | 2020-10-13 | The Intellectual Property Network, Inc. | Triggered responses based on real-time electroencephalography |
US11219405B2 (en) | 2018-05-01 | 2022-01-11 | International Business Machines Corporation | Epilepsy seizure detection and prediction using techniques such as deep learning methods |
US11273283B2 (en) | 2017-12-31 | 2022-03-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
US11364361B2 (en) | 2018-04-20 | 2022-06-21 | Neuroenhancement Lab, LLC | System and method for inducing sleep by transplanting mental states |
US11452839B2 (en) | 2018-09-14 | 2022-09-27 | Neuroenhancement Lab, LLC | System and method of improving sleep |
US11717686B2 (en) | 2017-12-04 | 2023-08-08 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to facilitate learning and performance |
US11723579B2 (en) | 2017-09-19 | 2023-08-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement |
US11786694B2 (en) | 2019-05-24 | 2023-10-17 | NeuroLight, Inc. | Device, method, and app for facilitating sleep |
Families Citing this family (54)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
US7747325B2 (en) | 1998-08-05 | 2010-06-29 | Neurovista Corporation | Systems and methods for monitoring a patient's neurological disease state |
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 |
US8762065B2 (en) | 1998-08-05 | 2014-06-24 | Cyberonics, Inc. | Closed-loop feedback-driven neuromodulation |
US7853329B2 (en) | 1998-08-05 | 2010-12-14 | Neurovista Corporation | Monitoring efficacy of neural modulation therapy |
US7596535B2 (en) * | 2003-09-29 | 2009-09-29 | Biotronik Gmbh & Co. Kg | Apparatus for the classification of physiological events |
US7276183B2 (en) * | 2005-03-25 | 2007-10-02 | Sarnoff Corporation | Metal silicate-silica-based polymorphous phosphors and lighting devices |
WO2007075477A2 (en) * | 2005-12-19 | 2007-07-05 | University Of Florida | Closed-loop state-dependent seizure prevention systems |
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 |
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 |
GB0611872D0 (en) * | 2006-06-15 | 2006-07-26 | Hypo Safe As | Analysis of EEG signals to detect hypoglycaemia |
EP2034885A4 (en) | 2006-06-23 | 2010-12-01 | Neurovista Corp | Minimally invasive monitoring systems and methods |
US7765170B2 (en) * | 2006-07-11 | 2010-07-27 | Samsung Electronics Co., Ltd. | Data segmentation algorithm using a decrementing sliding walker chunking approach |
US20080146958A1 (en) * | 2006-10-12 | 2008-06-19 | Kenneth Shane Guillory | Self-contained seizure monitor and method |
US20080091089A1 (en) * | 2006-10-12 | 2008-04-17 | Kenneth Shane Guillory | Single use, self-contained surface physiological monitor |
US20080091090A1 (en) * | 2006-10-12 | 2008-04-17 | Kenneth Shane Guillory | Self-contained surface physiological monitor with adhesive attachment |
US8295934B2 (en) | 2006-11-14 | 2012-10-23 | Neurovista Corporation | Systems and methods of reducing artifact in neurological stimulation systems |
WO2008058343A1 (en) * | 2006-11-15 | 2008-05-22 | The University Of Queensland | A method for detecting eeg seizures in a newborn or a young child |
US9898656B2 (en) | 2007-01-25 | 2018-02-20 | Cyberonics, Inc. | Systems and methods for identifying a contra-ictal condition in a subject |
US20080183097A1 (en) | 2007-01-25 | 2008-07-31 | Leyde Kent W | Methods and Systems for Measuring a Subject's Susceptibility to a Seizure |
US8036736B2 (en) | 2007-03-21 | 2011-10-11 | Neuro Vista Corporation | Implantable systems and methods for identifying a contra-ictal condition in a subject |
US9554721B1 (en) | 2007-04-23 | 2017-01-31 | Neurowave Systems Inc. | Seizure detector, brain dysfunction monitor and method |
FR2919406B1 (en) * | 2007-07-23 | 2009-10-23 | Commissariat Energie Atomique | METHOD AND DEVICE FOR RECOGNIZING THE POSITION OR MOVEMENT OF A DEVICE OR LIVING. |
US9788744B2 (en) | 2007-07-27 | 2017-10-17 | Cyberonics, Inc. | Systems for monitoring brain activity and patient advisory device |
EP2211708A1 (en) * | 2007-10-23 | 2010-08-04 | Optima Neuroscience, Inc. | System for seizure monitoring and detection |
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 |
US20110082381A1 (en) * | 2008-06-25 | 2011-04-07 | Basim M Uthman | Time Frequency Transformation Analysis for Detection and Quantification of Epileptiform Activity Load in Generalized Epilepsies |
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 |
GB0906029D0 (en) | 2009-04-07 | 2009-05-20 | Nat Univ Ireland Cork | A method of analysing an electroencephalogram (EEG) signal |
US8786624B2 (en) | 2009-06-02 | 2014-07-22 | Cyberonics, Inc. | Processing for multi-channel signals |
US9643019B2 (en) | 2010-02-12 | 2017-05-09 | Cyberonics, Inc. | Neurological monitoring and alerts |
US9558762B1 (en) * | 2011-07-03 | 2017-01-31 | Reality Analytics, Inc. | System and method for distinguishing source from unconstrained acoustic signals emitted thereby in context agnostic manner |
US9814426B2 (en) | 2012-06-14 | 2017-11-14 | Medibotics Llc | Mobile wearable electromagnetic brain activity monitor |
US20140275840A1 (en) * | 2013-03-15 | 2014-09-18 | Flint Hills Scientific, L.L.C. | Pathological state detection using dynamically determined body data variability range values |
US10929753B1 (en) | 2014-01-20 | 2021-02-23 | Persyst Development Corporation | System and method for generating a probability value for an event |
WO2017106490A1 (en) | 2015-12-15 | 2017-06-22 | Mayo Foundation For Medical Education And Research | Systems and methods for linear-time clustering for bounded, repeatable, rare events in physiological signals |
US11020035B2 (en) | 2016-02-01 | 2021-06-01 | Epitel, Inc. | Self-contained EEG recording system |
US10772527B2 (en) * | 2016-03-18 | 2020-09-15 | Robert Kocher | Brain matching |
US11152086B2 (en) | 2017-11-02 | 2021-10-19 | Bank Of America Corporation | Electroencephalogram triggered experience modification system |
US10567961B2 (en) | 2017-11-02 | 2020-02-18 | Bank Of America Corporation | System for electroencephalogram patterning recognition for authentication |
US10456054B2 (en) | 2017-11-02 | 2019-10-29 | Bank Of America Corporation | Electroencephalogram triggered resource distribution query system |
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 |
CN108852350B (en) * | 2018-05-18 | 2021-06-29 | 中山大学 | Modeling method for recognizing and positioning scalp electroencephalogram seizure area based on deep learning algorithm |
TWI680746B (en) * | 2018-11-21 | 2020-01-01 | 晶神醫創股份有限公司 | Method and system for neurological event detection |
CN110292379A (en) * | 2019-07-02 | 2019-10-01 | 重庆大学 | The recognizer of epilepsy non-invasive diagnosis based on convolutional neural networks |
US10743809B1 (en) | 2019-09-20 | 2020-08-18 | CeriBell, Inc. | Systems and methods for seizure prediction and detection |
US20210307672A1 (en) | 2020-04-05 | 2021-10-07 | Epitel, Inc. | Eeg recording and analysis |
CN114947885B (en) * | 2022-07-05 | 2023-07-04 | 华南师范大学 | Method for evaluating postoperative effect of radio frequency thermal coagulation damage epileptic focus |
US11857330B1 (en) | 2022-10-19 | 2024-01-02 | Epitel, Inc. | Systems and methods for electroencephalogram monitoring |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5311876A (en) * | 1992-11-18 | 1994-05-17 | The Johns Hopkins University | Automatic detection of seizures using electroencephalographic signals |
US5559929A (en) * | 1994-07-29 | 1996-09-24 | Unisys Corporation | Method of enhancing the selection of a training set for use in training of a neural network |
US5857179A (en) * | 1996-09-09 | 1999-01-05 | Digital Equipment Corporation | Computer method and apparatus for clustering documents and automatic generation of cluster keywords |
US5943661A (en) * | 1991-07-11 | 1999-08-24 | Texas Instruments Incorporated | Hybrid neural network classifier, systems and methods |
US6092059A (en) * | 1996-12-27 | 2000-07-18 | Cognex Corporation | Automatic classifier for real time inspection and classification |
US6157921A (en) * | 1998-05-01 | 2000-12-05 | Barnhill Technologies, Llc | Enhancing knowledge discovery using support vector machines in a distributed network environment |
US20010014776A1 (en) * | 1994-09-21 | 2001-08-16 | Nancy E. Oriol | Fetal data processing system and method employing a time-frequency representation |
US6282305B1 (en) * | 1998-06-05 | 2001-08-28 | Arch Development Corporation | Method and system for the computerized assessment of breast cancer risk |
US6324532B1 (en) * | 1997-02-07 | 2001-11-27 | Sarnoff Corporation | Method and apparatus for training a neural network to detect objects in an image |
US6336119B1 (en) * | 1997-11-20 | 2002-01-01 | International Business Machines Corporation | Method and system for applying cluster-based group multicast to content-based publish-subscribe system |
US6337997B1 (en) * | 1998-04-30 | 2002-01-08 | Medtronic, Inc. | Implantable seizure warning system |
US20020007237A1 (en) * | 2000-06-14 | 2002-01-17 | Phung Tam A. | Method and system for the diagnosis of vehicles |
US6442421B1 (en) * | 2000-04-27 | 2002-08-27 | Centre National De La Recherche Scientifique | Method for the medical monitoring in real time of a patient from the analysis of electroencephalograms to characterize and differentiate between physiological or pathological conditions, and a method for anticipating epileptic seizures |
US6549804B1 (en) * | 1996-01-23 | 2003-04-15 | University Of Kansas | System for the prediction, rapid detection, warning, prevention or control of changes in activity states in the brain of a subject |
US6594524B2 (en) * | 2000-12-12 | 2003-07-15 | The Trustees Of The University Of Pennsylvania | Adaptive method and apparatus for forecasting and controlling neurological disturbances under a multi-level control |
US6658287B1 (en) * | 1998-08-24 | 2003-12-02 | Georgia Tech Research Corporation | Method and apparatus for predicting the onset of seizures based on features derived from signals indicative of brain activity |
-
2002
- 2002-04-15 US US10/123,834 patent/US6735467B2/en not_active Expired - Fee Related
-
2004
- 2004-03-25 US US10/810,949 patent/US20040181162A1/en not_active Abandoned
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5943661A (en) * | 1991-07-11 | 1999-08-24 | Texas Instruments Incorporated | Hybrid neural network classifier, systems and methods |
US5311876A (en) * | 1992-11-18 | 1994-05-17 | The Johns Hopkins University | Automatic detection of seizures using electroencephalographic signals |
US5559929A (en) * | 1994-07-29 | 1996-09-24 | Unisys Corporation | Method of enhancing the selection of a training set for use in training of a neural network |
US20010014776A1 (en) * | 1994-09-21 | 2001-08-16 | Nancy E. Oriol | Fetal data processing system and method employing a time-frequency representation |
US6549804B1 (en) * | 1996-01-23 | 2003-04-15 | University Of Kansas | System for the prediction, rapid detection, warning, prevention or control of changes in activity states in the brain of a subject |
US5857179A (en) * | 1996-09-09 | 1999-01-05 | Digital Equipment Corporation | Computer method and apparatus for clustering documents and automatic generation of cluster keywords |
US6092059A (en) * | 1996-12-27 | 2000-07-18 | Cognex Corporation | Automatic classifier for real time inspection and classification |
US6324532B1 (en) * | 1997-02-07 | 2001-11-27 | Sarnoff Corporation | Method and apparatus for training a neural network to detect objects in an image |
US6336119B1 (en) * | 1997-11-20 | 2002-01-01 | International Business Machines Corporation | Method and system for applying cluster-based group multicast to content-based publish-subscribe system |
US6337997B1 (en) * | 1998-04-30 | 2002-01-08 | Medtronic, Inc. | Implantable seizure warning system |
US6157921A (en) * | 1998-05-01 | 2000-12-05 | Barnhill Technologies, Llc | Enhancing knowledge discovery using support vector machines in a distributed network environment |
US6282305B1 (en) * | 1998-06-05 | 2001-08-28 | Arch Development Corporation | Method and system for the computerized assessment of breast cancer risk |
US6658287B1 (en) * | 1998-08-24 | 2003-12-02 | Georgia Tech Research Corporation | Method and apparatus for predicting the onset of seizures based on features derived from signals indicative of brain activity |
US6442421B1 (en) * | 2000-04-27 | 2002-08-27 | Centre National De La Recherche Scientifique | Method for the medical monitoring in real time of a patient from the analysis of electroencephalograms to characterize and differentiate between physiological or pathological conditions, and a method for anticipating epileptic seizures |
US20020007237A1 (en) * | 2000-06-14 | 2002-01-17 | Phung Tam A. | Method and system for the diagnosis of vehicles |
US6594524B2 (en) * | 2000-12-12 | 2003-07-15 | The Trustees Of The University Of Pennsylvania | Adaptive method and apparatus for forecasting and controlling neurological disturbances under a multi-level control |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050080828A1 (en) * | 2003-10-08 | 2005-04-14 | Johnson Mark W. | Method for displaying spectral trends in complex signals |
US7395292B2 (en) | 2003-10-08 | 2008-07-01 | Johnson Mark W | Method for displaying spectral trends in complex signals |
US20090094305A1 (en) * | 2003-10-08 | 2009-04-09 | Johnson Mark W | System for Displaying Spectral Trends in Complex Signals |
WO2006122201A3 (en) * | 2005-05-10 | 2007-12-21 | Salk Inst For Biological Studi | Automated detection of sleep and waking states |
KR101395197B1 (en) * | 2005-05-10 | 2014-05-15 | 더 설크 인스티튜트 포 바이오로지컬 스터디즈 | Automated detection of sleep and waking states |
US11690557B2 (en) | 2005-05-10 | 2023-07-04 | The Salk Institute For Biological Studies Intellectual Property And Technology Transfer | Automated detection of sleep and waking states |
US20080102236A1 (en) * | 2006-10-27 | 2008-05-01 | Fish Robert B | Pipes containing nanoclays and method for their manufacture |
US10803145B2 (en) | 2016-02-05 | 2020-10-13 | The Intellectual Property Network, Inc. | Triggered responses based on real-time electroencephalography |
CN106821376A (en) * | 2017-03-28 | 2017-06-13 | 南京医科大学 | A kind of epileptic attack early warning system and method based on deep learning algorithm |
US11723579B2 (en) | 2017-09-19 | 2023-08-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement |
US11717686B2 (en) | 2017-12-04 | 2023-08-08 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to facilitate learning and performance |
US11273283B2 (en) | 2017-12-31 | 2022-03-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
US11318277B2 (en) | 2017-12-31 | 2022-05-03 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
US11478603B2 (en) | 2017-12-31 | 2022-10-25 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
US11364361B2 (en) | 2018-04-20 | 2022-06-21 | Neuroenhancement Lab, LLC | System and method for inducing sleep by transplanting mental states |
US11219405B2 (en) | 2018-05-01 | 2022-01-11 | International Business Machines Corporation | Epilepsy seizure detection and prediction using techniques such as deep learning methods |
US11452839B2 (en) | 2018-09-14 | 2022-09-27 | Neuroenhancement Lab, LLC | System and method of improving sleep |
US11786694B2 (en) | 2019-05-24 | 2023-10-17 | NeuroLight, Inc. | Device, method, and app for facilitating sleep |
Also Published As
Publication number | Publication date |
---|---|
US20030195429A1 (en) | 2003-10-16 |
US6735467B2 (en) | 2004-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6735467B2 (en) | Method and system for detecting seizures using electroencephalograms | |
Wilson et al. | Spike detection: a review and comparison of algorithms | |
Wilson et al. | Seizure detection: evaluation of the Reveal algorithm | |
Tzallas et al. | Epileptic seizure detection in EEGs using time–frequency analysis | |
Iasemidis et al. | Quadratic binary programming and dynamical system approach to determine the predictability of epileptic seizures | |
Sharanreddy et al. | EEG signal classification for epilepsy seizure detection using improved approximate entropy | |
Pincus et al. | Heart rate control in normal and aborted-SIDS infants | |
EP0898460B1 (en) | Systems for prediction, rapid detection, warning, prevention or control of changes in activity states in the brain | |
James et al. | Detection of epileptiform discharges in the EEG by a hybrid system comprising mimetic, self-organized artificial neural network, and fuzzy logic stages | |
Tarassenko et al. | Identification of inter-ictal spikes in the EEG using neural network analysis | |
Tzallas et al. | The use of time-frequency distributions for epileptic seizure detection in EEG recordings | |
Sameer et al. | ROC Analysis for detection of Epileptical Seizures using Haralick features of Gamma band | |
Slimen et al. | Epileptic seizure prediction based on EEG spikes detection of ictal-preictal states | |
De Lucia et al. | A novel method for automated classification of epileptiform activity in the human electroencephalogram-based on independent component analysis | |
Yadav et al. | Morphology-based automatic seizure detector for intracerebral EEG recordings | |
Cranstoun et al. | Time-frequency spectral estimation of multichannel EEG using the auto-SLEX method | |
US20010020137A1 (en) | Method and computer program product for assessing neurological conditions and treatments using evoked response potentials | |
Li et al. | An improved sparse representation over learned dictionary method for seizure detection | |
Tajmirriahi et al. | Modeling of seizure and seizure-free EEG signals based on stochastic differential equations | |
Rasoulzadeh et al. | A comparative stationarity analysis of EEG signals | |
Oliveira et al. | Multitaper-based method for automatic k-complex detection in human sleep EEG | |
Park et al. | Detection of epileptiform activities in the EEG using neural network and expert system | |
CN113974653B (en) | Method and device for detecting optimized spike based on about log index, storage medium and terminal | |
Shukla et al. | An Improved Decision Support System for Identification of Abnormal EEG Signals Using a 1D Convolutional Neural Network and Savitzky-Golay Filtering | |
Park et al. | Detection of epileptiform activity using wavelet and neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |