ALERTNESS AND DROWSINESS DETECTION AND TRACKING SYSTEM

Abstract: A drowsiness detection system (100) constructed according to the invention provides improved performance by preserving, and analyzing newly discovered information contained in rhythmic signal components above 30 Hz which in prior art has universally ignored or discarded as noise. In the first (700), and second (800) embodiments of the invention, one or more analog EEG signals are collected from a subject, appropriately filtered (616), converted into digital form (212), and subjected to frequency analysis. Selected signal components from the frequencies above 30 Hz, which have a high correlation with subject drowsiness, are isolated, and their amplitude, energy, or power contribution to the EEG signal is used to produce a continuous output measure reflecting the subject's alertness or drowsiness (140). The output measure (140) may be compared with a threshold to provide an indication of whether the subject is excessively drowsy.

Description: ALERTNESS AND DROWSINESS DETECTION AND TRACKING SYSTEM BACKGROUND GF THE INVENTION This invention relates to systems for determining a subject s state of , alertness and drowsiness, and more particularly to apparatus and methods for analyzing EEG signals acquireα from a subject to definitively determine the ό subject ' s state of alertness, drowsiness, or his or her sleep stage, and to determine whether the subject ' s performance is impaired thereby THE NEED FOR ALERTNESS MONITORING. DROWSINESS DETECTION. AND SLEEP STAGING : J Sleep deprivation has become one of the most significant causes of ;, error and accident throughout our society. The United States Department of Transportation estimates that 200.000 tr fic accidents each year may be 12 fatigue- or sleep-related. In transportation alone, sleep-reiated accidents annually n claim over 5,000 lives, and cause hundreds of thousands of injuries, with an is accumulated cost for health care, death, lost productivity, and damage to property d in the billions of dollars. U.S. Department of Health and Human Services (1992). j. ' Pilots say their schedules often force them to snooze in the cockpit in :_ order to get enough sleep. Industry insiders report that flight attendants need to ιι periodically check to ensure that the crew is awake. The National Traπsportalion :o Safety Board (NTSB) cited pilot fatigue as either the cause or a contributing 2 i factor in 69 airplane accidents from 1983 through 1986 Stanford Sleep Disorders Clinic and Research Center.(1991). Tj Recent analyses of spectacular accidents and catastrophes suggest that :•» sleepiness may have played an imDortant rcle n sucn events including the ?5 Three Mile Island nuclear disaster and the E' oπ Valdez oil SDIII Mitler et al :-, (1938) Such accidents endangered large segments of the population and the 2 * environment .3 The presidential Commission on the Space Shuttle Challenger accident ■9 ruled that ground crew fatigue was a contπbut.ng cause of the 1986 disaster in , the near catastrophic launch of the shuttle Columbia only three weeks before, operator fatigue was reported as one or the major factors contributing to this ,: incident Stanford Sleep Disorders Clinic and Pesearch Center (1991) Thus 23 the importance of human vigilance or attention is critical to the performance of i individuals in various types of occupations Each human being requires a sDecifir amount of sleep in each 24-hour JO period to maintain a functional level of alertness If an individual obtains less 27 sleep, he/she will be less alert the following day Moreover sleep loss accumulates s from one night to the next as a " sleep debt ' Therefore only a modest loss of j9 sleep per night may produce a serious sleep debt when sustained over several JO nights The more sleep lost each day the greater the sleep debt and the larger the 4i impairment Because individuals often do not recognize that they are sleepy, they 42 seldom guard against involuntary sleep episodes Much like intoxicated drivers 2 sleepy drivers do not realize that they are incapable of adequate performance, and 4 may therefore deny drowsiness and impairment U S Department of Health and 5 Human Services (1992). o The effects of sleep loss can be amplified by the bi-modal circadian rhythm 7 Evidence of this can be found in the temporal patterns of accidents attributed to .3 " falling asleep " or even to mere lapses in operator attention Studies of 49 vehicle truck accidents in Israel. Texas, and New York ail reveal two distinct peaks ■ύ in the time of day when these accidents occurred. La.'ieetal (1986). Langlois et si al. (1985); G V. Duff (unpublished observations) One peak occurs in the early 5 2 morning hours from l a m. to 7 a.m. and a lower peak occurs during the mid- 53 afternoon from 1 p m to 4 p m 54 Another factor which raises risk of accidents is the increasing level of 55 automation For example, drivers using cruise control and pilots using au omatic 5d flight control systems are more susceptible to drowsiness due to the removal 5/ of stimulating influences The Exxon Valdez was on automatic pilot during id the critical minutes leading to its grounding as it hit Bligh Reef at 12:04 a.m. 5 Stanford Sleep Disorders Clinic and Research Center (1991) The NTSB ' s 00 investigation of the accident indicated that the third mate was asleep on his 01 feet and failed to respond to the warning light and alarm identifying the reef 6? U.S. Department, of Health and Human Services (1992). Although automation ,j has provided tremendous benefits, it tends to limit operator activity to vigilant o monitoring of the system. Over a period of time, this can reduce the awareness 65 level of the operators and impair their ability to react properly to an external stimulus. In addition, vigilance is further degraded by sleep loss and fatigue. ? Thus, it would be highly desirable to produce an automated real-time όs system to track the changes in levels of alertness, such as the transition from 69 alertness to drowsiness, or the onset of sleep. In addition, there are a number of other applications in which an automated system for measuring an individual's alertness drowsiness cr stage of sleeo would be highlv useful For example, sleep stagin — i e , the identifica ion of a subject s stage cr condition of sleep based on physiological indicators — is used clinically fzr cagnosmg and treating sleep disorders Sleep staging is also of interest m medicai research. Normally, sleep staging is performed by a highly trained physician or technician by studying voluminous EEG records collected while a subject sleecs A totally automated system for sleep staging could improve consistency and reduce research and treatment costs. Although the sleep scoring field is weii established, the greatest disagreement among sleep scorers analyzing identical segments of sleep data occurs when scoring the transition from ' stage W " (a state of wakefulness) to " stage 1 " sleep (an initial stage of sleep sometimes referred to as the sleep onset s: period) 3 No system is currently available which can effectively use the EEG signal 34 for continuous drowsiness tracking and detection A recent report to the United 5 States Department of Transportation (DOT) surveying methods of drowsiness n detection acknowledged that automated processing of the EEG signal has proved 37 very difficult to implement. Wierwille (1994) Various phases or stages of sleep 3 are identifiable using automated methods. However, drowsiness and the onset 9 of sleep are much less distinguishable in the EEG waveform, and therefore, i much more difficult to identify using automated methoGs. than are sleep stages. 1 Research surveyed in the DOT report suggests using a manual method of analyzing EEG and EOG signals. Wierwille (1994), c/r/πjjPlaπque (1991). 9 It should be noted that in many real-worid applications it is insufficient j4 to detect sleep, as normally understood because it is often essential to provide ' ■5 a warning before an individual s performance is impaired In particular for -,n critical applications in which a lack of vigilance could affect health and safety it •v is necessary to detect extreme sleepiness ' Extreme sleepiness ' is used herein >s to refer to the state during which sleep is perceived as difficult to resist, the 15 individual struggles against sleep performance lapses occur and sleep will 100 eventually ensue but has not yet occurred By detecting the onset of extreme lϋl sleepiness in an individual, the individual may be alerted or disqualified from ιo: service, before they reach a state in which they are incaoable of safely performing a task PRIOR-ART APPROACHES TO ALERTNESS MONITORING /us A variety of methods have previously been used or pro posed for detec ing W6 or evaluating sleep or drowsiness in a subject Although some prior-art methods ιo7 have been reasonably successful at automated detection of actual sleep, there toe has heretofore been no automated system capable of consistently and definitively :o9 detecting the onset of extreme sleepiness before an individual becomes unable to no safely perform a task. ALERTNESS MONITORING USING EXTERNAL MANIFESTATIONS OF u2 DROWSINESS n Several prior-art approaches to the automated sleep detection problem have in relied on externally observable manifestations of sleep. For example, Chiu U.S. us Patent No.4,875,030 discloses a system which observes the state of a subject ' s .l eyelids while performing a task sucn as driving If the subject s eveύds remain closed for a period greater tπaπ a normal blιnκ interval the suoiect is determined . to be succumbing to slee p and an ^arm is gιve π This s /ste could be extended n to sound an alarm when some othe blink-related characteristic such as blink ac duty cycle or blink frequency deviates from established norms Kishi U S. Patent ι:ι No 5,311,877 discloses a system for estimating a waking degree (which might be a measure ofalertness) using an indi idual s eve blink frecue n cv or the time 122 required for the individual to react to a visual stimulus Estrada U S Patent No 3953,831, discloses a system which attempts to monitor the attitude of the as suoject s head, if the head is observed to droop then the subject is determined i.', to be succumbing to sleep and an alarm is given Slansky U S Patent No z.v 4,617559, discloses a fatigue alarm system which employs a pressure-operated i23 switch disposed in a wrapper for a steering wheel or the like when the grip of the i29 user becomes relaxed, the switch operates and the alarm is given i2o These prior-art approaches suffer from a number of important ij disadvantages, especially in real-time applications involving health or safety A U2 primary disadvantage is that these approaches oo not detect drowsiness or lack i23 of vigilance early enough If a subject closes his or her eyes, or allows his or her i3 head to droop, and that behavior is the result of drowsiness, then the subject us may already be at the far end of the drowsiness spectrum, and performance may 136 already be impaired In certain applications, eye closure or head drooping is i37 meaningless in a drowsiness detection system, because such behavior is neither u3 necessary nor sufficient for a conclusion of drowsiness or sleep For example, there ay are some environments in which the subject is permittee to rest his or her eyes for :J3 brief intervals, provided that he or she remains awake and vigilant Furthermore. :n detection of eye closure or head droopin produces a binary output which has :i no further sensitivity once the targeted behavior has been detected Thus, in the aforementioned example, once the subject closes his or her eyes, the system . • 4 cannot distinguish among wakefulness, drowsiness or sleep. .-45 Another disadvantage of these approaches is that it is relatively difficult O in practice to usefully monitor a subject ' s physical activity, such as the state of 1 ,' a subject ' s eyes or head position This disadvantage is further compounded by ,J3 the difficulty in using information about physical activity to distinguish fatigue 4 from transient, but normal, variations in behavior Estrada, for example, discloses no using a mercury switch to monitor the position of the subject ' s head Such a i5 switch provides a discrete-valued output, is difficult to adjust, and may produce is: false indications during normal subject and vehicle movement. Although imaging is systems have been proposed for capturing and interpreting an image of the i54 subject ' s eye or eyes, such systems require expensive image and signal processing (55 components. Further, any imaging system must contend with: various eyewear .-so and clothing which the subject may employ (e.g., prescription glasses, sunglasses, ιs contact lenses, hats); large variations in normal lighting conditions; contamination 153 of the vehicle environment by high-amplitude spurious lighting (e.g., illumination ιs9 by the headlights of another vehicle); and normal movement by the subject. too Seko et al. U.S. Patents Nos. 4,564.833 and 4.604.611 disclose drowsiness i5i detection systems for motor vehicles which detect the onset of sleep in a driver by .j observing a change in the number, rate or amplitude c ' certain steering inputs o> When a vehicle is being driven along a linear path the civer typically makes 0 4 frequent steering corrections bv performing a series or Ξ ~ a 11 displacements of .5 the steering wheel in either direction because these displacements typically so are small in magnitude and frequently are in the direction opposite that of the o previous displacement, they are sometimes referred to as ' micro-reversals." As oi a driver falls asleep the resolution of the driver s stee r ι"g control degrades This (o is detectable as a change in the frequency and amplituce of steering inputs. As i -J a driver becomes sleepy and loses attentiveness or if the driver falls asleep for a / * ι brief instant and then awakens the driver will have failed to provide appropriate i - steering inputs during that interval Upon awa ening or regaining attentiveness. i,j the driver attempts to rαpidlv supply steering inputs which correct the entire i,-4 steering error that accumulated during the sleep interval Thus, changes in the i5 pattern of steering reversals may indicate that the driver has fallen asleep or is 1,-6 about to do so i77 This method of detecting rowsiness or sleep also has a number of (-5 disadvantages It is difficult to distinguish between abnormal changes in steering i patterns caused by the onset of sleep and normal steering pattern changes 130 required by road or traffic conditions which may mimic crowsiness- or sleep- isi induced changes. As a result, systems which rely on this method are prone to false i32 alarms. The method also fails to take into account characteristics of individual i3j drivers. In addition, because this method relies on measurement of the subject ' s i3 actual task performance, it does not detect the onset of sleep until the subject ' s S performance is already noticeably impaired Such detection may not provide sufficient warning to avoid an accident Furtnεrmore. manv tasks which require vigilance do not employ any user inputs under normal conditions, other tasks as require some user inputs but such inputs may not form recognizable patterns .'59 from which abnormal user behavior may be distinguished. :90 ALERTNESS MONITORING USING INTERNAL MANIFESTATIONS OF ai DROWSINESS Other prior-art automated sleep detection approaches have attempted to measure directly one or more of a subject ' s internal physiological characteristics which may indicate alertness or drowsiness. Yoshimi et αl U.S. Patent No as 4.928.090 discloses a system for judging " arousal level " based on a measurement uo of skin potential level A disadvantage of this system is that skin potential level . , • may be affected by many factors other than arousal or drowsiness, and therefore 198 it is difficult to distinguish drowsiness-related changes in skin potential level and (59 changes caused by other factors 2oo EEG-BASED APPROACHES oι The electroencephalogram (EEG) is a recording or " the low-voltage 0 electrical activity produced in specific regions of the brain. The EEG provides 3 a powerful tool for studying both normal and abnormal brain function, and has 04 been commonly used to measure and define wakefulness and sleep. There is 05 considerable evidence that physiological sleepiness is directly related to the 26 rapidity of the onset of EEG-defined sleep. Torsvall et al. (1987); Torsvall et al. 20 (1989); Akerstedt etal (1990), Akerstedt et al (1991). Wierwille et. al. (1992); 08 Dingus et al (1937) 309 Sleep is regarded as an active and complex state and has been 2io characterized in terms of various stages and cycles The term sleep architecture " :tι is used to describe these stages and cycles The stages and cycles of sleep may be .'.•j defined using both external manifestations and internal physiological processes u which can be externally measured such as EEG and electrooculogram signals i4 The electrooculogram (EOG) is a recording of the low-voltage electrical 3is activity associated with eye movement It has Peen observed that the transition to 3i6 sleep is also frequently accompanied by slow rolling eye movements (SEMs) that 3i can be detected in EOG signals 13 EEG signals include both periodically recurring or rhythmic features 21 (waves), and transient features such as spindles which do not recur on a 320 periodic basis and are highly localized in time One method of summarizing 32i and evaluating the content of EEG records is through analysis of frequency i2 components contained in the EEG signal Traditional EEG doctrine states that 223 the information content of an EEG signal is band-limited between roughly 05 23 and 30 Hz. Within this range, standardized names have been given to individual 325 frequency bands. (See Table 1) TABLE 1 TRADITIONAL EEG SIGNAL FRE: 'UENCY BANDS 337 Band Lowei r Upper 2 3 Name Limit (Hz) Limit, Hz) 339 Delta 05 2 230 Theta 3 7 321 Alpha 8 12 333 Beta 13 30 33 A substantial amount of research has been conducted in an attempt to j characterize the relationships between a subject ' s EEG .vaveforms the subject ' s j75 state of alertness, drowsiness or sleep, and the subject s ability to perform a 2 jo task The relationship between performance degradation (slower reaction rates 237 and attention lapses) and increased sleepiness has been established by several j3 researchers. Wilkinson and Houghton (1975), O ' Hanlon and Kelley (1977), 239 Diπges (1988); Molodofsky (1992), Trejo and Shensa (1993), Makeig and Inlow 2 0 (1993). Similarly, a strong correlation has been noticed between performance 4i degradation and particular patterns in the EEG waveform Horvath etal. (1976), 42 O'Hanlon and Beatty (1977); O ' Hanlon and Kelley (1977), Makeig and Inlow 3 (1993). In turn, the EEG waveform has also been correlated with vigilance and 44 sleepiness of subjects in various studies Gale (1977), Daniel (1967), Fruhstorfer 45 etal. (1977); Santamaria and Chiappa (1987) Kishi U S. Patent No. 5,311,877 46 purports to employ " brain waves " a'πd performance measurement in A system 3 for estimating a " waking degree " Although the meaning of the term " -waking 313 degree " is unclear it is treated as equivalent to react. c, time (see Fig. 11), and it .19 might be a measure of alertness Kishi discloses the use of " brain waves ' analyzed 350 by a " brain wave processor " to supply inputs to the waK πg degree estimation .-si unit However, the particular analysis performed is unclear, and it is likewise .'5.' unclear what relationship may exist between a subject s cram waves and his or her waking degree, cr how any sucn relationship might De exploited ;5 Sleep researchers rely on the EEG in the classificαt.on of various phases or 55 stages of sleep However compared to sleep stages drowsiness and the onset .'bo of sleep a re not as easily distinguishable in the EEG waveform In particular, 57 when drowsiness or sleep onset occur, the changes directly apparent in the EEG .'53 waveform are less profound, nd m y be masked by even s or processes which 59 affect the EEG but which are not directly related to drowsiness 60 ALERTNESS MONITORING USING EEG FREQUENCY ANALYSIS oι One commonly applied method of analyzing EEG signals has been to study 63 the frequency bands in which significant or predominant components of the 63 signals reside. The term " predominant " is used herein to refer to the frequency 64 bands or components which contain most of the energy in EEG signals. The 65 research to date indicates that drowsiness is associated with a re-distribution of OO- energy in the traditional frequency bands toward lower frequencies However, it 6 is difficult to use this general guideline in a drowsiness detection system because 63 of the varied characteristics of EEG signals among subjects For example, subject 69 behavior and physiology unrelated to drowsiness can produce changes in the EEG JO which are similar to those correlated with drowsiness Additionally whether the eyes are open or closed can make a substantial difference in the level of activity m particular frequency bands of interest in the 27 EEG Based on sleep stage classification studies Smith claims that the EEG can be used to identify states of severe sleepiness if the subject ' s eyes are open Smith 375 (1987) In alert subjects with their eyes open the predominant energy in the .v. EEG signal appears in the beta frequency band (13-30 Hz) A shift of energy .' into the alpha band (8-12 Hz) occurs as the subjects become drowsy O ' Hanlon j-3 and Beatty (1977) For an individual with his or her eyes open studies clearly 39 indicate that increases in alpha and theta activity in the EEG may correlate with 3so sleepiness as well as reduced performance Fruhstorfer et al (1977), O ' Hanlon 3t and Kelley (1977), Daniel (1967) Horvath et αl (1976), O ' Hanlon and Beatty 32 (1977) Thus, the appearance of alpha activity in the EEG may be an indicator of 23 drowsiness (an incipient indicator of a loss of vigilance), if the subject has his or 3 her eyes open. Santamaria and Chiappa (1987), Makeig and Inlow (1993) 85 If the subject ' s eyes are closed, however, it is much more difficult to 2 s6 differentiate between sleepy and alert states. Smith (1987). in individuals with 237 eyes closed, EEG signal energy is predominantly located in the alpha frequency ?33 band even though they are wide awake. As sleepiness develops, a subject whose 39 eyes are closed generally experiences a reduction in alpha band energy and an 290 increase in theta band energy (and possibly delta band energy) present in the 21 occipital channel of an EEG. Smith observed a pattern of SEMs during sleepiness 292 with "open eyes" and reported that the studied subjects alternated between 3 open and almost closed eyes Smith (1987) This may be an explanation for the 394 increased alpha activity observed d ring drowsiness Thus, when the EEG of an 2 5 initially-alert individual with eyes open displays a s ift in predominant energy from 296 the beta band to the alpha band, that shift might indicate that the individual has 29 become drowsy, but it also might simply indicate that the individual has closed his 398 or her eyes but otherwise remains alert. O ' Hanlon and Beatty (1977); Makeig and .'9 Inlow (1993), Santamaria and Chiappa (1987) 3oo A number of known techniques seek to characterize the subject ' s state 2oι of alertness or sleep by observing the re-distπbutioπ of energy in the subject ' s JO: EEG signal among the traditional frequency bands. Some of these refinements 303 are directed to examining ratios of the energies present in two or more frequency 304 bands, or examining the ratio of the energy contained in a predefined frequency os band to the total energy in the 0-30 Hz band. 3 06 Other techniques of EEG signal analysis, which do not principally rely on 307 the frequency domain, have also been tried. Sleep spindles and K-complexes are J O3 perturbations in the EEG signal which are associated with stage 2 sleep and are 309 generally recognizable in the time domain. But the presence of sleep spindles and JZO K-complexes occurs too late to be useful in detecting drowsiness. Researchers in have also tried to use SEMs. which are detectable in EOG signals, as an indicator i 2 of drowsiness. However, none of these have produced a reliable indicator of 3i3 drowsiness. 3i4 Thus, although it is clear that there exists some information in EEG signals 3i5 which indicates drowsiness, and although some prior-art methods have been 3iά reasona bly successful a t a u toma ted detection of actua l sleep, heretofore there has been no automated system ca pa ble of consisten tly πc definitively detecting the :. ■ ? onset of extreme sleepiness (which functions as a precursor to sleep). 21 9 INADEQUACIES IN PRIOR-A RT METHODS OF ACQ UIRING. RECORDING 130 AND ANALYZING EEG SIGNALS Although the acq uisition a nd recordi ng of EEG signals from h uman subjects has been practiced for many years conven tion a; acquisition a nd analysis techniques have not resulted in strong a nd consistent correlations between the 324 EEG signal and alertness. Fig 1 is a simplified block digram of a system 60 typical of those used for recording and analyzing EEG signals for the primary purpose of -36 medical diagnostics. At least one signal electrode 64 a nd a t least one reference electrode 66 333 are coupled to a su bject 62 from whom it is desired to receive EEG signals. τ9 Although only one signal electrode 64 is shown , it is common in clinical and 230 research applications to place signal electrodes at several ' standard locations 33i on the subject ' s head to obtain EEG informa tion which correlates with various 33 physiological a nd or cognitive processes tha t occur in different regions of the 33 brain. The electrodes 64, 66 are typically mecha nically secured to the subject ' s j scalp to provide a relatively low impeda nce electrical con nection therewith; often, jj5 a conductive gel or paste is used to en hance the conductivity of the connection. 336 Although one electrode 66 is designated as the " reference " electrode, such jj7 designation is arbitrary because the EEG signal is measured differentially. The j s electrodes 64 a nd 66 and certain related signal acquisition components (to be 339 further described) may represent a plurality of duplicate sets of electrodes and 1 related signal acquisition components 24i In both clinical and researcn applications the trsc tioπal means of 242 displaying and recording EEG signals has been to suopi. the amplified signals to : 3 a "pen recorder " 74 as shown in Fig 1 A pen recorder t/pically has a strip of 74 paper or chart medium 76 which moves at a predetermined speed in a longitudinal ; ? direction corresponding to the time axis of the chart One or more pens 88 are ji mounted to trace one or more continuous lines on the cr.art as the chart moves 247 longitudinally A suitable mechanism causes transverse c.solacement of the pen j4 responsive to a corresponding electrical signal input so tnat the pen traces a j record of the amplitude of the signal input over time 350 With respect to EEG signals the human body is a relatively high- jsi impedance source which produces very low voltage signals (in the microvolt 5? range) at very low currents In order to provide sufficient signal amplitude to drive J 5J the pen recorder 74, the EEG signals must be amplified The electrodes 64, 66 35 are connected to an EEG amplifier 70 using suitable leads 68, which are selected J55 to minimize noise contamination Leads 68 may be provided with a shield 58 5s as an additional noise reduction step The EEG amplifier 70 is typically a high- 3S quality high-input-impedance linear amplifier having a several user-selectable 53 gam settings in the range of lO-^-lO 5 . A variety of EEG amplifiers are available 359 as commercial products and in various configurations, one EEG amplifier popular 3oo in clinical applications is a component of a commercially available product sold άi under the name " Grass instruments Model 12C Neurodata Acquisition System. " 2 Traditional EEG doctrine states that the information content of an EEG 263 signal is band-limited between roughly 05 and 30 Hz. n addition pen recorders jo necessarily have limited bandwidth because they incorporate electro-mechanical JO5 transducers Accordingly, conventional EEG acquisition systems include a suitable 366 filter (shown in simplified form as filter 70 of Fig. 1) as 03rt of the amplifier 3=7 electronics Typically, a first-order Butterworth type filter with a 50 °O amplitude :=7 response at approximately 30 Hz is used Fig 2 is a graph 88 showing the 1- 2o9 500 Hz frequency response 90 experimentally derived from a commercial EEG 2 filter of the first-order Butterworth type, which is available as a component 7,-ι of the aforementioned "Grass Instruments Model 12C fJeurodata Acquisition 37 2 System. " As best seen in Fig. 2. such a filter has a -3 dβ point at approximately 3.-3 21-22 Hz, and a filter roll-off of approximately 20 dB per decade Also, in many 3 74 environments, there may be a large amount of electrical noise at the commercial 3 7s power-line frequency (60 or 50 Hz). Substantial contamination of EEG signals 3 76 can occur at these frequencies, and therefore, most commercial EEG equipment 3 77 provides additional filtering to minimize the effect of this contamination 73 The analysis of EEGs generally involves the opinions of expert clinicians in 379 conjunction with an off-line data analysis procedure. In order to display and record 330 the signal, the amplified and filtered EEG signal 78 is typically supplied to the pen 33i recorder 74 In some applications, it may also be desirable to record or process 32 the EEG signal in electronic form The amplified and filtered EEG signal 78 may 333 be supplied to an optional post-processing system 80 for storing and processing 23 the signal. The post-processing system 80 may employ any of a wide variety of 735 storage means 82. such as instrument tape recorders and digital storage systems. 336 Any suitable post-processing functions 84 may be applied to the EEG data stored 23 in storage means 82. The results of processing the EEG data may be provided on a 333 lead or data path 86. 3 In some applications, the signal 78 may be converted using conventional 70 techniques into digital form for storage on any suitable digital storage medium. :?( When a band-limited analog signal is periodically sampled for conversion into 3)2 digital form, according to yquist ' s theorem the sampling frequency must be 97 at least two times the highest frequency in the analog signal Accordingly, in 94 applications where EEG signals have been recorded and processed digitally, 7 5 sampling rates have generally been in the range of 90-275 Hz EEG examin tions 6 may employ one or several electrodes Because the examinations are conducted 39 7 over periods of five minutes to several hours, a large amount of data must be 98 collected, stored, and processed. Accordingly, it is highly desirable to minimize J the sampling rate to the extent possible consistent with EEG signal bandwidth in 0 0 order to minimize the amount of data storage and post-processing required. 01 Fig.3 is a graph 92 showing an excerpt of a sampled EEG waveform 94 4 02 acquired from a human subject using conventional EEG recording techniques as joj generally shown in Fig. 1. This sample waveform was obtained from the occipital 4 0 EEG channel of an awake subject during an alertness test in which the subject was os presented with a time-separated series of visual stimuli and scored on their ability 4 06 to promptly respond to each stimulus. This sampled waveform corresponds to 07 a successful response. This waveform shows the amplitude of the EEG signal in 403 microvolts over a 2-second interval and was acquired a: a sampling rate of 256 Hz , 4,9 resulting in 512 data samples io Fig 4 is a graph 9 showing the power spectral density 93 of the sampled HI waveform depicted in Fig 3 over the frequency range -128 Hz The energy 412 within the 0-30 Hz range accounts for more than 99% cf the total spectral energy. 4( Although approximately 1% of the total spectral energy :s contained in the 4 frequency range above 30 Hz it is not visible on the grapπ of Fig 4 which is 4(5 scaled for viewing the predominant frequency content i.n the 0-30 Hz range. i6 An ιnterestιn feature of graph 96 ιs the pαιrof promιπent peaks56aπd 4(7 58 appearing in the power spectrum in the r nge of 9- II Hz, which is within the (3 alpha frequency band These peaks account for approximately 505% of the total ( energy in the spectrum Established EEG doctrine asserts that predominant 420 alpha band energy in the EEG of a subject whose eyes are open is an indication 2i of extreme drowsiness and is often associated with poor task performance. 422 Surprisingly, although alpha band energy was present (and even predominant) 3 in this sample, performance was not only acceptable, but corresponds to one 4. of the fastest reaction times this subject achieved. (The fact that the subject 42s responded to the visual stimulus, and corresponding data from a simultaneously 26 acquired vertical electrooculogram (EOG) signal, verify that the subject ' s 427 eyes were open.) This type of contradiction between a subject ' s alertness, as .s demonstrated by task performance, and that predicted from the subject ' s EEG 429 signals according to traditional doctrine can be seen in other test samples. o Therefore, an alertness monitoring system that relies heavily on the presence or 471 shift of the predominant energy in rhythmic EEG s.'gna 5 among the frequency bands that have traditional!/ been of interest in ΞΞC -essarcr. is prone to error in the form of disagreement bet.veen predicted and a::_a behavior It is taught throughout the EEG literature απc assumed in the design of 475 prior-art alertness monitoring systems that all of t _.=ful informa ion in the EEG signal is contained in the 0-30 Hz frequency bα-c For example. Kishi J7 U.S Patent No 5.311877 purports to employ bra - .-.αves ano performance measurement in a system for estimating a ' waking az'te " However, Kishi ' s brain wave processor provides band-pass filters to se'ec: only the traditional EEG frequency bands, and is thus not directed to analysis c: EEG signal components above 30 Hz. 44.' It is believed that all methods, analyses, and systems built on the 447 hypothesis that the useful information in the EEG signal is contained in the 0-30 Hz band share the following characteristics: 445 1. the analysis focuses on the frequency banc of the EEG signal containing the predominant energy; 447 2. the analysis is directed to rhythmic signal activity at frequencies 443 below approximately 30 Hz (i.e.. exclusive y in the standard frequency 449 bands according to conventional EEG literature), and 450 3. the analysis treats as noise, or otherwise c .cards signal components 451 above appro<imately 30 Hz (e.g., by low-cass filtering). 45.T We have observed that higher-frequency components of EEG rhythms, 45J although not the predominant components, nonetheless contain information 54 useful for automated monitoring of alertness and drowsiness 455 Although some researcπ n the EEG fie'd has employed -wider niter o bandwidths and higher sampling rates to preserve frequency components in 57 the EEG signals above 30 Hz sucn research has generally been directed to - • 3 the detection of physiological conditions other than drowsiness or sleep or has 5; primarily relied on analysis of transient or non-periodic events There are several 4oo well-known transient events during which higher frequency components (above oi 30 Hz) are present in thr- EEG signal Common sources cf high frequency energy 4o.' are most often associated .vith burst activity sucn as the neurological "spike " 4 ) activity observed during an epileptic seizure, and muscle ' " movement artifact. o4 Muscle/movement artifact is an extremely common source of high-frequency 4ό5 contamination found in nearly every EEG record Sleep spindles are brief bursts 4o of 12-14 Hz activity in the EEG signal which have been associated with Stage 4o7 2 sleep. Sleep spindles are events which are highly localized in time and are o3 detectable primarily using time-domain analysis Because they occur only late 4 0 in a subject ' s sleep cycle they are not generally useful in monitoring alertness or 0 drowsiness. i Researchers trying to identify spike activity in the EEG (commonly required 2 in epileptic research) have increased the bandwidth of their low pass (anti- 47j aliasing) filters and have used faster sampling rates during digitization Because 47 the hig_h frequency components in spike or burst activity are highly localized and 45 sharply defined in the time domain, researchers use increased filter bandwidth 7ό and sampling rates to capture more of the frequency content of the spike, thereby 47 improving identification, analysis, and reconstruction in the time domain 4 a Researchers interested in transient time-domain events generally have not 49 employed frequency analysis techniques to study them For example Olsen et al. 4,70 U.S. Patent No 5,311,876 discloses a system for automatic seizure detection 81 using EEG signals, and uses a sampling rate of 200 samples per second However, 43-τ Olsen is directed to detecting seizures not drowsiness, and attempts to classify 431 events or features in the signals using time-domain analysis techniques such as 434 counting the number of signal extrema which occur in an analysis epoch. 75 The transient time-domain events to which the aforementioned EEG 4 ό research has been addressed should not be confused with the continuous 4i7 appearance of high-frequency energy or " rhythms " We are unaware of any prior 433 art establishing the presence of sustained high frequency rhythms in EEG signals 439 or relating such rhythms to alertness, drowsiness, or sleep. 90 The current literature discourages those skilled in the art from exploring 1 the usefulness of high frequency rhythmic activity in EEG signals. For example, .' Gaillard refers to frequencies above the beta band as ' high frequency noise " . 493 Gaillard (1987), 9-11. Pritchard refers to the higher frequency signals as white 494 or near-white noise. Pritchard (1995), 378. Carskadon and Rechtswchaffen state 5 that an upper filter cut-off frequency in the range of 30 to 35 Hz will generally 6 pass through the " essential " wave forms, while minimizing high frequency 497 interference. Carskadon and Rechtswchaffen (1987), 668. O ' Hanlon and 9β Beatty also refer to " noise " greater than 30 Hz. O ' Hanlon and Beatty (1977), 9 195. Thus, all of the prior art relating to EEG-based sleep detection has either 5oo ignored evidence of higher frequency rhythmic activity in the EEG. or has failed : ι to recognize that such activity may oe usefully correlated with alertness and J O drowsiness in a practical monitoring or detection system 5 3 OBJECTS AND SUMM RY OF THE INVENTION 0 It is therefore an object of the present invention to provide an automated 5 05 alertness and drowsiness monitoring system which avoids the disadvantages of -co prior-art systems 50 7 It is another object of the invention to provide an automated alertness and 50 3 drowsiness monitoring system which reliably provides an indication when a subject 50 becomes excessively drowsy ,o It is a further object of the invention to provide an automated alertness and i d drowsiness monitoring system which reliably provides an indication of drowsiness 5 j before a subject becomes incapable of safely performing a task 5 ij It is another object of the invention to provide an automated alertness and 5 i drowsiness monitoring system which reliably provides as an output a continuous 5 i 5 measure representing the alertness or drowsiness of--a subject. i It is a further object of the invention to provide an automated alertness and 5 i drowsiness monitoring system which provides an indication of a subject ' s alertness sis or drowsiness based on information in a subject ' s EEG signal, including frequency 5.-9 components above 30 Hz. 5: 0 It is another object of the invention to provide an automated alertness and 52 i drowsiness monitoring system which provides an indication of a subject's alertness 5 .τ or drowsiness based on information in a subject ' s EEG signal, including non- 577 predominant components thereof 534 It is a further object of the invention to provide n automated alertness and 75 drowsiness monitoring system which provides a reliable indication of a subject ' s 52r > alertness or drowsiness from a subject s EEG signal anc which js suitable for use 7 in real-time applications 523 A drowsiness detection system constructed according to the present ■ 2) invention avoids the aforementioned disadvantages of fe oπor art by preserving 570 and analyzing newly discovered rhythmic signal components in selected frequency 51 bands which the prior art has universally ignored or discarded as noise " 577 An EEG-based system for monitoring or detecting alertness, drowsiness. 533 and sleep is provided which exhibits improved performance over prior-art systems 54 in detecting the onset of drowsiness in a human subject be fore the subject actually 5J5 succumbs to sleep or suffers a performance failure The system is referred to 536 herein as a drowsiness monitoring or drowsiness detection system, although 5 J the invention may also find application in alertness monitoring, sleep staging, 5J3 state-of-coπsciousness monitoring, anesthesia monitoring, and other related 5J applications. According to the invention, subject drowsiness is strongly correlated 50 with the energy present in certain rhythmic components of the subject ' s EEG S4i signal at frequencies above 30 Hz. 5 4 2 An EEG-based drowsiness monitoring system constructed according to the 5 4 present invention includes' acquisition components, signal analysis components, 5 44 artifact detection components, and threshold components The signal acquisition 5 4 5 components sense the subject ' s EEG signal amplify the signal for further analysis, 54ό and filter certain signal components which apparently do not contain useful 5 information and which degrade further processing and analysis steps The signal 543 analysis component receives the amplified ano filtered EΞG signal, determines 543 the amplitudes or energies of the components located in several predefined i5 frequency ranges, including at least a portion of the frecuency range 30-500 Hz, i and determines an output measure signal which represents the drowsiness of the 57 subject. The artifact detection components examine the subject s EEG signal, 557 and optionally examine secondary indicators of the subject ' s physical activity, 55 and determine therefrom whether the subject s EEG signal acquired during a 555 particular time interval is likely to be contaminated by artifact, and therefore 55o should not be used in drowsiness detection The threshold component establishes 557 a threshold against which the output measure of the signal analysis component 553 may be compared to determine whether the output measure indicates that the 559 subject is excessively drowsy The threshold may be determined for the subject 5όo under examination using EEG signals collected from the subject in a known state sάi of alertness, or may be a universal threshold applicable to the population as a 56 whole. 5sj In first and second embodiments of an EEG signal analysis component 54 constructed according to the invention, primarily digital signal processing 565 techniques are used. The analog EEG signal is sensed, amplified, and low-pass ■ filtered for anti-aliasing. Next, the analog signal is converted to a digital signal 5 by an analog to digital converter system. The digital signal may be analyzed oπ- •8 line in real time, or may be stored for off-line processing and analysis. During 5-9 analysis, a series of small overlapping win dews or b a tcr.es of data corresponding to bri f time intervals of the signal are sequen ially selected Frequ ncy analysis rast r-cι r a n s r o r m i is us d to conv rt irr.e-comain sig onal into a frequency domain output vector B cause the ΞΞG signal ' is non -sta ionary, a windowing function is applied prior to frequency analysis The power spectral 5 density (or " power spectrum " ) of the FFT output vector is determined, in order to 5-5 obtain the power in eacn frequency componen The power spectrum components 5 c are grouped into a small number of pre-seiectec spectral bins whicr, correspond -:- tc predefined free nancy ranges The spectral bin components are aggregated 73 to determine the total energy in each bin, and each of the resulting spectral bin 5. '9 energy values is inverted Weights are applied to the respective inverted spectral 530 energy values An " output measure ' signal is determined as the sum ot the 53i weighted inverted energy values. The output measure is a continuous-valued 53.T signal indicating the drowsiness of the subject The output measure may be 537 compared with the aforementioned threshold to produce a simplified output signal 53 indicating, for example, that the subject is approaching extreme drowsiness, or is 535 asleep, or is likely to be incapable of safely performing a task. 536 In a third embodiment of an EEG signal analysis component constructed 537 according to the invention, primarily analog signal processing techniques are used. 533 The signal processing components are broadly analogous to the digital signal 539 processing components provided in the first and second embodiments 5 BRIEF DESCRIPTION OF THE DRAWINGS 591 These and other features of this invention will be best understood by 59? reference to the following detailed description of a ρre'V r ed of the invention r.aκen in conjunction with tne accompanying drawings .n whicn Fig 1 ,s a bloc^ O'agram of an EEG recording s stem typ cai of those used according to the prior art 5.0 Fig. 2 's a graph snowing the frequency response cf a low-pass filter used in the prior-art ΞEG recording system of Fig 1 Fig 2.s a graph showing an excerpt of an EEG waveform acq ired using an EEG recording system according to prior-art techniques as shown in Fig 1 and employing a ιθw-pass filter having the response shown in rig 2 Fig is graph showing the power spectral density calculated over the frequency range 0-128 Hz from the EEG waveform whicn is depicted in Fig 3 and σ07 which was acquired using the prior-art EEG recording system and low-pass filter of o04 Figs. 1-2, o05 Fig. 5 is a simplified block diagram showing the general structure of first oOo and second embodiments of an EEG-based drowsiness monitoring system which is o07 constructed according to the present invention and which is adapted for use with 608 digital signal processing and related techniques; o09 Fig.6a is a simplified schematic diagram showing a low-pass filter which 610 may be used for anti-aliasing in the inventive EEG-baseo drowsiness monitoring oil system of Fig. 5 612 Fig. 6b is a graph showing the frequency response of the low-pass filter of olJ Fig. 6a, for use in the inventive EEG-based alertness monitoring system of Fig.5 Fig 7 is a graph showing an excerpt of an EEG waveform acquired using an !5 experimental embodiment of tne inventive alertness m onitoring system of Fi ' σ. 5 - o Fig 8 is a grapn snowing t~e power spectral c nsit/ calculated over the - • frequency range 0-123 Hz from the EEG wavercrm ..r en is depicted in Fig 7 and which 'was acquired using the ι /entive alertness mc tcrmg system of Fig. 5; ->) Fig 9 is a graph showing the power spectral de n sity calculated over the - ; frequency range 31-475 Hz from the EEG waveform ..men is depicted in Fig 7 and wnich was acquired using a expeπme tal emoccimeπt of the inventive alertness monitoring system of F'g 5, - 7 Fig 10 is a graph showing the power spectral c = n sιty (PSD) calculated , by averaging the PSDs over the frequency range !00--~5 Hz from 25 EEG data 0-5 segments collected during experimental trials in whicn a subject successfully -Ό responded to a presented stimulus o37 Fig. 11 is a graph showing the power spectral density (PSD) calculated 0:3 by averaging the PSDs over the frequency range 100--75 Hz from 20 EEG data o29 segments collected during experimental trials in which a subject failed to respond oo to a presented stimulus, o i Fig. 12a is a block diagram showing the configuration of the subject OJ interface portion ofa first embodiment of the inventive drowsiness monitoring o3 system shown generally in Fig. 5, the embodiment being adapted for collecting 074 drowsiness-related EEG data in a clinical or research environment- ojs Fig. 12b is a block diagram showing the configuration of the information oj-o processing and control portion of a first embodiment of the inventive drowsiness cj monitoring system shown generally in Fig. 5 the embodiment being adapted for o7 collecting drowsiπess-re'ated EEG data in a chπicai or researcn environment n Fig 13a is a oiocκ diagram snowin the ccr • u'αtic of the subject ? interface portion or a second e oociment cr he ι / π t, /e d rows mess monitoring o~ι system shown gene r ally ,n Fig 5 the embedment ee g adapted for collecting o .τ drowsiness-related EEG data in a stand-alone or tasΛ-based environment- o 7 Fig 13b is a block diagram shewing the con:,guration of the information processing a no con t r oι portion of a second e oocim n of the inventive 045 drowsiness monitoring system shewn generall / in F g Ξ the e ocdiment being r4 adapted for collecting d r owsιness-related EEG data <n a stane-a'one or task-based o4 environment, o4.ι Fig 14 is a blocκ diagr m of a first embodiment of an EEG signal processor 049 component which may be used with the inventive d'owsiness detection system 05 of Fig. 5, the signal processor component being arranged to use digital signal 5i processing techniques and in conjunction with the clinical or research environment 652 of Figs. 12a-12b, 053 Fig. 15 is a block diagram of a second embodiment of an EEG signal 054 processor component which may be used with the inventive drowsiness detection 055 system of Fig. 5, the signal processor component being arranged to use digital ose signal processing techniques and in conjunction with the stand-alone or task- o57 based environment of Figs 13a— 13b, o Fig. 16 is a data-flow diagram depicting, in simplified form, the processing o5 of EEG signal information acquired from a subject to produce a useful output oo measure indicating the drowsiness of the subject as that information is operated upon by the processin means of the embodiments cf Figs 5 ana 1--I5 OO Fig. 17 is a block diagram of a third empooimentcf aπ Ee: --eased - i drowsiness monitoring system wnicr, is const'-cted according to fe present oo invention a d w ich is adapted for use with a ~aiog signal process ng and related COO techniques; ooό Fig. 13 is a block diagram of a method for use π conjunction -with the o inventive drowsiness detection system cf Fig Ξ for detecting artifacts in EEG oo signal data collected from a subject, o Fig 19 is a block diagram cf a method ' ' or use in conjunction with the o o inventive drowsiness detection system of Fig Ξ for determining a drowsiness 7i threshold against which a drowsiness measure cfa subject may be comp red to o 2 determine when the subject s performance is iikely to be impaired 37 Fig. 20a is a graph depicting the output measure produced by an 674 experimental embodiment of the inventive drowsiness monitoring system during 675 two separate periods of an examination ofa human subject, during the first of o7ό which the subject was alert and performing a visual test, and during the second of 677 which the subject was in bed going to sleep; and 678 Figs. 20b and 20c are graphs depicting the output measure produced by an o7 experimental embodiment of the inventive drowsiness monitoring system during c o an examination ofa human subject, in which the output measure accurately o3i predicted the subject ' s performance failures. 682 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS o3j According to one aspect of the present invention, an EEG-based alertness 634 and drowsiness monitoring and detection system is prcioed which exhibits cos improved performance over prior-art systems in detecting the onset of drowsiness oio in a human subject before the subject actually succumes tc sleep cr suffers a 37 performance failure The inventive system lso provides improved performance o3s in a variety of ether applications reouiππg monitoring eetectioπ or continuous o39 tracking of a subject s state of alertness drowsiness cr sleep. The inventive o o system also may provide improvements in applications relating to other aspects of i a subject ' s alertness which may not necessarily relate to sleep, such as monitoring o aspects of a su eject s state of consciousness or response to anesthesia Thus, o97 although the system disclosed herein will generally be referred to henceforth as a 04 " drowsiness monitoring system " cr ' drowsiness detection system. " the scope of 0 5 the present invention shall not be limited by the use of this shorthand terminology. oό Broadly defined, an EEG-based drowsiness monitoring system constructed 6 according to the present invention comprises- means for acquiring one or more 69β EEG signals from a subject; means for filtering the acquired signals to eliminate o 9 signal components which are unimportant or spurious, while retaining signal 7oo components having a usable correla ion with the subject s state of alertness or 7oi drowsiness; means for processing the signals to select particular components 702 which are relevant to alertness and drowsiness and to measure characteristics 73 of the selected signals; means for eliminating certain signal information which, 704 based on analysis of the EEG signals or on external information, appears to be 05 contaminated by events unrelated to the subject ' s alertness or drowsiness; and 70s means for determining an output measure which indicates the subject ' s alertness 37 or drowsiness. Optionally, the system may also comprise means for comparing the 03 output measure to a threshold to produce a simplified output signal indicating, 27 for example, that the subject is approaching extreme drowsiness, cr is asleep, or :o is likely to be incapable of safely performing a task. ιι An EEG-based drowsiness monitoring system constructed according to a the present invention may take several different forms ~ herefore, this patent application discloses three preferred embodiments or the invention, from which i4 an appropriate one may be selected depending on the particular environment in which the invention is to be applied and on the cost and availability of 16 implementing technologies. 17 Fig. 5 is a simplified block diagram showing the general structure 100 i3 of a drowsiness detection system constructed according to an aspect of the i9 present invention. Fig. 5 is a generic drawing which is applicable to all three 2 specific embodiments with the understanding that everything outside block 2i 124 is optional and may or may not be present in the specific embodiments. 2 Figs. 12a, 12b, 14, and 16 show the structure of a first preferred embodiment 2 700 of the invention. Figs. 13a, 13b, 15, and 16 show the structure ofa second 4 preferred embodiment 800 of the invention. The first and second preferred s embodiments share many similarities in their methods of acquiring, processing, 6 and analyzing EEG information, and may be primarily implemented using digital 7 signal processing and related technologies. Thus, these two embodiments will 8 often be discussed together in this application. Fig. 17 is a simplified block 9 diagram showing the structure of the third preferred embodiment 400. The third -7o preferred embodiment of the invention may be pπmaπi ; implemented using analog signal processing and related technologies Although these three embodiments mav be implemented using different technologies according to one aspect of this invention they all share the features -7 of including particular EEG signal components which have heretofore been -;5 discarded or ignored in accord with the teachings of the prior art and exploiting -a the information contained in these components to provide a reliable measure ofa subject s alertness or drowsiness The analog signal processing functions of the third embodiment are analogous to those performed using digital signal 9 processing in the first two embodiments Accordingly, the first and second 40 embodiments will be discussed first, and the third embcciment will be discussed 74( subsequently with reference to analogous signal processing functions of the first 74? and second embodiments 74j The first embodiment, which is shown in greater cetail in Figs 12a, 12b, 74 and 14, may be most appropriate for use in a sleep laboratory, clinical sleep - 5 analysis, or other medical or laboratory applications in which it may be desirable 746 to acquire large amounts of EEG data for post-processing and/or on-line analysis 7 which is not necessarily limited to the recognition of alertness, drowsiness or sleep. 3 The second embodiment, which is shown in greater detail in Figs. 13a, 13b, and 4 15, may be most appropriate in a stand-alone alertness/drowsiness monitoring 750 and alarm application For example, the second embodiment might be used as a 75i self-contained, real-time device for on-line monitoring of the alertness/drowsiness 75? ofa sonar or radar operator, a nuclear power plant or industrial process control- 753 room operator cr a vehicle operator 54 In general the variations between these two embodiments relate to the -55 scale and structure of the data acquisition and processing components These 7S variations are the result cf differences in the amount cf information which may be 757 feasibly collected and the uses to which the raw data and analyzed results are put 753 in the respective application environments Fig 15 is a simplified data-flow-type " 59 diagram showing the transformation of signal information at various stages, and is -o0 common to both embodiments 7oi The simolihed block diagram of Fig 5 presents a basic platform for 7o? realizing an EEG-based drowsiness detection system 100 which is constructed 7o3 in accordance with an aspect of the present invention The system 100 of 764 Fig 5 may primarily employ either digital or analog signal processing and 7≤5 related technologies Fig 5 is applicable to all three preferred embodiments. 66 The three embodiments discussed herein differ primarily in their respective 767 application environments, and in the technology used to implement an EEG signal 763 processor component 124 The signal processing components outside block 124, 69 particularly the threshold means and the artifact detection means are optional 770 and may or may not be present in the specific embodiments. Further, according i to another aspect of the invention, a drowsiness analysis component 136 could 772 be used independently to analyze EEG signals acquired and stored by other EEG 7 examination systems. 77 The system 100 comprises suitable means (such as electrodes 112 and 75 114) for obtaining one or more EEG signals or "channels " 116 from a subject 110, 776 amplifier means 118 for receiving and amplifying the EEG signals, signal analysis means 120 for receiving the amplified EEG signals on lead 130 and producingan 773 output signal 142 indicating that the subject is excessively drowsy (or another 779 suitable result of the analysis), and means 122 for presenting a display or alarm 730 indicating the result of the analysis. 711 The signal analysis means 120 may be conceptually divided into three main 732 functional components EEG signal processor means 124 receives analog EEG 737 signals on lead 130 from amplifier means 118, filters and converts the signals into a plurality of digital samples representin the EEG signals analyzes the 735 digital samples, and responsively produces an output measure 138 indicating the subject ' s state of alertness, drowsiness, or sleep An optional artifact detection 73! means 128 receives the digital samples from EEG signal processor means 124, 783 and analyzes the samples, and possibly other information, to determine whether 739 the samples are apparently contaminated by artifacts, and therefore cannot be 790 reliably used. An optional threshold means 126 generates a suitable threshold 791 value (on lead 146) which is compared with the output measure 138 produced 792 by the EEG signal processor means 124 The result of the co parison may be 793 a simplified output signal 142 indicating that the subject has reached or passed 794 a threshold stage of drowsiness or sleep. The threshold means 126 may use as 795 threshold value 146 a universal value believed to be generally applicable to the 796 human population, or may optionally determine the threshold value 146 for an 797 individual using baseline EEG signals acquired therefrom. 793 The conceptual division of signal processing and analysis means 120 into 799 several functional components is useful in describing its operation but when 3co implemented, signal analysis means 120 need not have separate corresponding oi physical components For example, depending on application recuirements, all 30 of the functions of signal analysis means 120 could be implemented using one sc3 or more high-performance computer-based systems cr a single- board signal so processing system or even a single-chip digital signal processor 305 EEG signais may be ootained from a human subject 110 using any 3ύo appropriate means, of which several methods are well known in the art The most so- commonly used method oi obtaining EEG signals is to apply suitable electrodes. 803 such as electrodes 112 and 114 (Fig 5), at various locations on the subject ' s so-) scalp in order to detect particular rhythmic EEG signals known to carry desired 3io information sii The " International 10-20 System " of electrode placement has become the aa standard instrument in research and clinical neurophysioiogy The 10-20 System 3i3 determines electrode locations based on the size of the subject ' s head and is 3i4 therefore specific to individuals Use of the 10-20 System of electrode placement 3i5 assures accurate and repeatable placement of electrodes for an individual across sio multiple occasions and examination facilities, and allows the comparison of 3i EEG signals between subjects. The detailed head measurements and electrode sis application techniques will not be described here, but are disclosed in many 3i9 handbooks on EEG technology. Some commonly used location names are: frontal 820 (F), central (C), parietal (P), occipital (O), and anterior (A). It is believed that 321 among the standardized scalp locations, the location pairs O^-A 1 and O^-A^ are 333 good for use in detecting drowsiness 332 Typical EEG eiectrode connections in clinical applications may nave an 324 impedance in the range cf 5-10 K ohms It .s generally desirable to minimize 335 the impedance of the connections between tne electrode and tre sub j ect. 330 Laboratory-grade electrodes may be used to form connections with impedances 3 7 below 2 K ohms A conductive caste or gel may be applied to the electrode to further improve the conductivity and mecnaπical st bility of t e connection 829 Commercially available active electrodes which provide an amplifier on or near 830 the electrode may also be used Needle electrodes may be applied subcutaneously 831 in laboratory applications 333 In some other applications the mechanical configuration cr stability 333 of the electrodes, or the ease of applying the electrodes to the subject, or the 834 compatibility of the electrodes with the subject ' s mobility, are high priorities 835 Capacitively coupled electrodes may also be used Additionally a lightweight EEG 836 sensor could be advantageously provided in a suitable headset (not shown) having 837 a sufficient number of electrodes 833 Because the electrical currents of the EEG signals produce associated 839 magnetic fields, methods have been proposed for sensing these magnetic fields to 840 acquire signals equivalent to conventional EEG measurements without electrical 841 contact with the subject. Such signals are referred to as magnetoencephalogram 842 (MEG) signals. One of ordinary skill in the art will appreciate that the present 843 invention could be used with MEG signals with little or no modification. 844 To maximize clarity, the block diagram of Fig.5 depicts collection and 345 processing of a single EEG cnannel However cne skilled in the art will appreciate 86 that an inventive rowsiness detection system eo iα ac y an ageo us /collect and 84- process multiple EEG cnanne's 1 ' desιreα collection are processing of multiple 343 EEG channels may be acccmpnsπed by simply replicating all cf the components 8 Q shown on a per-cπannel basis Alterna ively one could use multiplexing, computer 550 m ulti- tasking and other teehmques which are known in the art for processing o5i multiple signals with a smalier numDer of processing means aπc signal patns 35? The numoer of EEG signals to be processed depends on the environment ssj in which tne drowsiness monitoring system is to oe used In laboratory or 35 clinical environments, it may be desirable and relatively convenient to collect 355 and process EEG signals obtained at several locations on the subject s scalp 35o Skilled technicians are avail ble to apply e electrodes a d porta ility of the 357 signal acquisition and analysis equipment is not a high priority In stand-alone 353 environments, where a subject is to be monitored while performing a real-world 35 task, technicians may not be available to apply the electrodes, and subject 360 mobility requirements may limit the size and weight of the signal acquisition 8όi and analysis equipment In such environments, it may be feasible to collect 862 and process only a few channels, or only a single channel, of EEG data When 863 processed and analyzed according to the present invention a single channel of 864 EEG data is sufficient to reliably track or detect the onset of extreme drowsiness 865 in a human subject. see The number of electrodes required depends on the number of EEG 867 signals to be processed At least two electrodes (one signal electrode, such as 363 electrode 112 (Fig 5) and one reference electrode such as electrode 11*1) are dό recuirεd to obtain one EEG cnanne! The number of e-ectroces needed to obtain 37 greater numoers of EEG signal channels depends on wnether separate reference a-.- electrcdes or a single reference electrode is used. s τ ? As noted previously the EEG signals as measured are very low voltage 3-7 signals (in the microvolt range) In order to provide sufficient voltage levels for 8 further processing, the EEG signals must ce amplified -πy suitable cabies or 375 wires 116 may be used to connect the electrodes 112 anc 114 to an appropriate 8 o EEG amplifier means 118 Care should be taken to minimize interference from 37 electrical noise sources For example, cable 116 may include a shield conductor 373 174 37a EEG amplifier means 113 may be implemented using a suitable high-quality 380 amplifier having a high input impedance and sufficient ga to amplify the EEG 3 i signals for input to the signal analysis means 120 As wiil be discussed further, 88 the signal analysis means 120 is nominally configured to receive an output signal asj on lead 130 in the ±2.5 volt range from amplifier means 118. However, another 33 appropriate voltage range could be selected for the amplified signal, provided 335 that both the output of amplifier means 118 and the input of the signal analysis ssό means 120 are compatible High linearity, low distortion fiat frequency response. 887 and good common-mode rejection are desirable characteristics for amplifier ass means 118. EEG signal levels available at the input of amplifier means 118 may 889 vary, depending on the subject, the type of electrode (or other probe device) ago used, and the qualityof the connections obtained Preferably, amplifier means 39i 118 provides several operator-seiectabie gain settings to accommodate such 39? variations Because the amplifier me ns 113 will be electrically connected to a 397 person, amplifier means 113 must be designed and constructed consistent with 504 applicable safety standards for such equipment. 35 Several commercially available amplifiers which have been specifically 55 developed for use in collecting EEG signals and meet these requirements, could be 397 used One such amplifier popular in clinical EEG applications is a component ofa 398 commercially available product sold under the name " Grass Instruments Model 3? 12C Neurodata Acquisition System " In stand-alone applications where low oo cost and subject mobility are high priorities, commercially available operational oi amplifiers ( " op-amps " ) in single-chip or module form might be used in this 9 ? application One skilled in the art would appreciate how to select a suitable op- oj amp (and support components) or other suitable means, such as a portable EEG 0 system, for use in this application 05 As best seen in Fig. 5. an amplified EEG signal is provided on lead 130 to 06 an EEG signal processor means 124 which is a component of the signal analysis 907 means 120. The signal on lead 130 is an " analog " signal The first and second 90s preferred embodiments of this invention, are preferably implemented using digital 99 signal processing and related technologies. Accordingly, in those embodiments, 9io the analog EEG signal must first be digitized— that is, converted into a series of 9ii digital samples which represent the original EEG signal with sufficient accuracy. 9i3 EEG signal processor means 124 comprises an EEG signal acquisition means 9i 132 (including a filter 164 and an analog-to-digital (A/D) converter 212) and 9i4 a drowsiness analysis means 136 The filter 164 filters tne EEG signal and the 9.5 converter 212 samples ano converts tne EEG signai into cigita form and may i store the converted samples i he drowsiness analysis means 136 analyzes the 5i samples and produces therefrom a output signai on lead 133 which serves as a 9iε measure ' of a subject s state cf alertness drowsiness or sleep 9!-> Prior-art attempts at automated drowsiness detection using EEG signals 7c have systematical! /discarded or ignored EEG signal components at frequencies 3.- ' greater than approximately 30 he a factor which we believe nas been a significant 933 contributor to the failure of such systems According to one aspect of the ?3 present invention reliable and accurate information concerning a subject ' s ?4 state of alertness drowsiness or sleep is extracted from EEG signal components 9?5 in frequency ranges well above tnose considered relevant by traditional EEG ?o doctrine Prior-art EEG recording systems have filtered out these frequencies 92 as " noise. " According to this invention, these frequencies are preserved and 9s analyzed. 929 We have determined experimentally that particularly userul information 30 resides in the frequency range of approximately 80-475 Hz Accordingly, the 93i frequency range of 80-475 Hz is discussed herein merely as an example of 9j 2 frequencies above 30 Hz which are of interest in the analysis of EEG signals 933 for information concerning a subject s state of alertness, drowsiness, or sleep. 934 Drowsiness detection systems which use this exemplary frequency range may 935 be constructed consistent with the preferred embodiments of the invention 96 disclosed herein, and would provide significantly improved performance over pπor- 93 art systems -owever it is believed that at least some of the advantages cf the 973 present in vent, en may be cc tamed using any reasoπac'v selected subranges of frequencies between approximately 30 Hz and ΞC0 H: In additicn frequencies 9JO above 30 Hz - hich are outside the exemplary frequency range cf 80-475 Hz may 9 i also be fouπc serul upon further study, and if used might enable drowsiness 94? detection systems to provide satisfactory or improved performance. Therefore. 94j embodiments of the invention may be designed to a ecu ire and analyze an 9 4 expanded range of frequencies, even above 475 Hz 95 In order to analyze such high frequency information, the EEG signal must 6 be sampled at a sufficiently high rate during digitization, and a properly designed 7 anti-aliasing filter must be used to condition the signal prior to sampling. It 9 4 3 is known that a band-limited signal can be uniquely represented by a plurality of regularly spaced samples if the rate f s at which samples are taken is at least 950 twice the highest frequency f max in the original sample For example, in order 9 5 1 to sample a signal containing components at frequencies up to an exemplary 95 ? maximum freouency f ma χ of 475 Hz. the minimum sampling rate f s would be 5j 950 samples per second (SPS), but any sampling rate greater than 950 SPS 5 could be used. In a prototype embodiment of the invention, a sampling rate f s 955 of 950 Hz was used to collect data (discussed below In detail) concerning the 5ό performance of the invention in detecting drowsiness. Although the sampling rate 957 may be increased to the limits of available technology, faster sampling increases 953 the cost of equipment and the amount of signal data which must be processed. 59 Inexpensive commercially available equipment could be advantageously employed 960 in preferred embodiments cf the invention to allow exploitation cf EEG signal oi information at higher frequencies o? Drowsiness monitoring system ICO preferably comprises suitable anti- 9-J aliasing low-pass filter means 164 for attenuating signal components above 9o the maximum frequencies desired for crowsiπess detection thereby providing 9όϊ a band-limited signal for digitization Filter means 164 is shown in Fig.5 as an 9oό integral part of EEG signai acquisition means 122. but it may be located elsewhere 97 provided that it operates on the EEG signal before sampling, and could be, for 63 example integrated with the EEG signal amplifier 113 It is conventional to 969 employ a low pass filter 70 (Fig 1) in an EEG recording system, and a variety of 970 such filters are commercially available However, the commercially available filters 97i designed for conventional EEG recording applications are not preferred for anti- 9? aliasing because they employ a low-order design (such as first-order Butterworth) 93 resulting in a shallow filter roll-off (see Fig. 2) 9 If the previously mentioned exemplary frequency band of 80-475 Hz were 975 selected for use with the first and second preferred embodiments discussed herein, 976 a minimum sampling rate f s of 950 Hz could be used. A low-pass (anti-aliasing) 977 filter providing 50 dB of attenuation in the stop band would be sufficient. A sixth 978 order Butterworth filter having a -3 dB point selected at 182.1 Hz provides at least 9 50 dB of attenuation at 475 Hz and above and can be affordably implemented 980 with commercially available components. This type of filter has a more rapid 931 transition from the pass band to the stop band than the shallower first-order 982 response of filters typically implemented in EEG equipment. 53? Fig. 6a is a schematic diagram of a low-pass anti-aliasing filter means 164a 53 which has tnese characteristics and which is suitable for use in conjunction with 92- the drowsiness monitoring system 100 (Figs 5 - 15) of the present invention. 98o Fig 6b is a graph 644 shewing the frequency response 646 of the low-pass filter of 95 Fig. 6a As best seen in Fig cc the attenuation of the filter 164a begins within 98S the frequency band of interest However subsecuent processing may be used to 989 compensate for this attenuation (see filter compensation means 230. Figs 14 and =0 15), or a higher sampling rate could be used to allow the filter s -3 dB point to be 9i moved beyond 132 1 Hz thereby extending the pass band of the filter 9? As best seen in Fig 6a. anti-aliasing filter means 164a comprises a suitable 99? power supply 612 (which may be shared with other signal acquisition equipment), 99 an optional input buffer amplifier 614, an optional output buffer amplifier 620. 995 a filter module or circuit 616, and an oscillator module 618 for determining 996 the cut-off frequency of the filter module 616. Power supply 612 may be any 997 suitable regulated DC power supply providing at least a first output " V- " 626 993 at approximately -5 volts, and a second output " V— " 628 at approximately -j-5 999 volts, referenced to a common ground 654 Power supply 612 is preferably of the woo "linear" type which generally provides reduced noise at the power supply outputs toot (compared to " switching " power supplies). ιao2 A suitable power supply 612 is commercially available from Power-One, looj Inc., 740 Colte Piano, Camaπllo, California 93012, under the designation ιoo4 "International Series, Model HCC5-6/OVP-A, " but other power supplies could 1005 also be used. The power supply outputs 626, 628, and ground reference 654 loco may be supplied to other data accuisitioπ components such as analos-to- ιoc digital conversion means 212 (see Figs 14 and 15) to allow the analog-to-digital loos conversion means 212 access to a low-noise source of power without providing ιoo9 an additional power supply The power supply 612 may receive power from any wio suitable source, such as leads 622 and 624 connected to a domestic AC power ion source In a vehicular application a suitable DC-DC power supply could be used, ioa and leads 622 and 624 could be connected to the vehicle electrical system In an aa application where portability is desired, power supply 612 could be replaced with a 7014 suitable battery. ioi5 Any suitable high-performance low-pass filter circuit 616 providing i sufficient attenuation could be used to implement the filtering function of ioa filter means 164a In a preferred embodiment of the invention, a sixth-order ιoi8 Butterworth-type switched capacitor low-pass filter may be used as the filter ιoi9 circuit 616. Such a filter is commercially available, packaged in a single integrated ισ?σ circuit (IC), from National Semiconductor Corporation, Santa Clara, CA under ιo i the designation " LMF60CIN-100 " However, other filter circuits could also be ιo2 used. ιo3 The LMF60 switched capacitor filter circuit advantageously provides a ιo high-performance filter and requires few external components. The filter circuit ιo5 616 uses an externally derived clock signal supplied on lead 634 to determine its ιo26 cut-off frequency. The cut-off frequency is obtained by dividing the clock signal by ιo27 1D0. The LMF60 ' s sixth-order Butterworth design provides a minimum stop-band io8 attenuation rate of approximately 36 dB per octave. In an embodiment which 1C29 employs a sampling rate f 5 of 950 Hz, the resulting sampling system upper limit 1020 frequency f. 2ιs< 5h ' z A cut-cff frequency of 2C0 Hz is suitable for use with lc?. ' this filter to provide at least 50 d ' B attenuation at 475 Hz Accordingly, a 20 KHz icj-r filter clock signal may be provided en lead 634 to seiect a 200 Hz cut-off One of 1022 skill in the art will appreciate that filters having different cut-off frequencies and lc?-- stop-band attenuation rates could also be used as cesired provided that sufficient 175 attenuation is provided at the sampling system upper limit frequency f s /2. 1030 Filter means 164a comprises an appropriate oscillator circuit 618 which 10J7 provides the clock signal for filter circuit 616 Any suitable 20 KHz oscillator iojs could be used In a preferred embodiment of the invention, a crystal-controlled ιo ' i oscillator having an integrated externally-programmable divider chain may be ιo4o used to provide the 20 KHz clock signal An integrated crystal oscillator and 10 ! divider circuit 618 which may be used in the filter means 164 of the present io2 invention is commercially available from Epson America, Inc. under the ioj designation " SFG-8651B " As best seen in Fig 6a, oscillator 618 has a plurality 04 of control inputs 632 by which operation of the internal divider chain may be 105 programmed. In order to produce an output of 20 KHz with the SPG-8651B ιo46 integrated circuit, control inputs Cl and C3 must be connected to V-f- via lead ιo47 628, and control inputs C2 and C4-C6 must be connected to ground via lead ιo48 654. The output signal from oscillator 618 is supplied to filter 616 on lead 634. o9 The use ofa programmable divider advantageously allows a variety of different ιo5o clock frequencies to be selected as desired for a particular application. Other wsi suitable oscillator/divider circuits which can produce different selections of output frequencies are commercially available With reference once again to Fig 5 the EEG Signal from EEG amplifier 115 1J5- is provided to filter means 164 on lead 121 The anti-a.iased EEG output sisπal 2C55 from anti-alias filter means 164. is supplied to signal sampling and conversion I?:; means 212 (Figs 5 14 15 I on lead 310 As best seen π Fig. όa. the EEG signal 105- _ιs preferably buffered by an input buffer amplifier 61- ~ea 130 is connected to 1053 the input of the input buffer amplifier 61- The output of the Duffer amplifier 614 15? is connected to the input of the filter means via lead 632 The input impedance loό of the filter circuit 615 may be undesirably low The input buffer amplifier 614 loάi may be used to present α higher input imoedance to the EEG signal supplied on ιoo? lead 130. Any suitable amplifier could be used for input buffer amplifier 614 The ιo6 LMF60 filter IC includes two uncommitted operational amplifiers, of which one ιo6 could be used as a buffer amplifier 614 However, another commercially available 1065 operational amplifier could also be used Alternately, lead 130 may be connected 106 directly to the input of filter circuit 616 ιoό7 An output buffer amplifier 620 is preferably provided to isolate the filter ιo63 circuit 616 from the filter means output lead 310. The filter circuit output lead ιo69 636 is connected to the input of the output buffer amplifier 620: and the output ιoo of amplifier 620 is connected to the filter means output lead 310. Any suitable ioi amplifier could be used for output buffer amplifier 620. The LMF60 filter IC io3 includes two uncommitted operational amplifiers, of which one could be used 703 as an output buffer amplifier 620. However, another commercially available ιo4 operational amplifier could also be used. Alternately, filter means output lead 310 1075 may be connected directly to tne output of the filter circuit 616. ιo75 The anti-aliased EEG signal output from filter means 164 is provided to io7 signal sampling and conversion means 212 (Figs 5.14. and 15) on lead 310. 107$ Signal sampling and conversion means 212 may be implemented using any io5 suitable analog-to-digital converter system capable cf sampling at the desired rate 1080 and providing sufficient resolution (or precision) Analcg-to-digital conversion at 1 3 a resolution of 12 bits (i e 4096 discrete values) is adequate, but conversion at 1083 higher resolution (such as 16 bits) may result in greater overall performance of the io33 drowsiness detection system 103 In clinical or laboratory applications, such as the environment in which the 1085 first preferred embodiment 700 (see Figs 12a. 12b) of the invention may be used, 1086 a general purpose computer system may be employed as an acuisitioπ control and 1037 processing means 216a (Fig 12b) for EEG data collection The computer system waβ is preferably equipped with an interface card or accessory containing one or more ιo3 analog-to-digital converters and suitable control software to drive the converters 1090 and save the converted samples The EEG signal is converted into digital samples 1091 for subsequent analysis using digital signal processing techniques. For example, 1093 many small "personal computers " may be used. ιo93 A suitable analog-to-digital interface card for use in a general-purpose ιo94 computer is commercially available from Keithley MetraByte, Inc..440 Myles 1095 Standish Boulevard, Taunton, MA, 02780, under the designation " DAS- 1096 1802HR." As noted above when processed and analyzed according to the present ιo7 invention, a single channel of EEG data is sufficient to reliably detect the onset ιo93 of extreme drowsiness in a human subject In addition a samoiing rate of 950 105'- samples per second appears to ce adequate to accuirε an exemplary range of HOC frequencies of interest in drowsiness detection However it may oe advantageous uoi to collect EEG signal data obtained from a plurality of electrode locations, or no? to collect data at higher rates The aforementioned interface card is capable of no: samoiing and converting up to 8 fully differential analog channels to digital data iio4 at a resolution of 16 bits at an aggregate rate of 100000 samples per second If al ιι 5 eight channels are used the interface card can accommodate a sampling rate of iioo approximately 12500 samples per second which allows signals up to 6.250 Hz to ιιo7 be sampled. Other analog-to-digital interface cards could also be used nos In a stand-alone application such as the environment in which the second 1709 preferred embodiment 800 (see Figs 13a, 13b) of the invention may be used, mo it may be a priority to implement a drowsiness detection system 840 in an iiu inexpensive and highly portable package. Signal sampling and conversion means a 212 may be implemented using a small general-purpose computer which may be aa equipped with an analog-to-digital interface card as described above ιιi As an alternative, signal sampling and conversion means 212 could be uis implemented using any appropriate single-board computer which is equipped me for signal processing applications, of which several are commercially available ιu7 For example, Texas Instruments produces several board-level products for ma digital signal processing applications which may be used A suitable single- ιιi9 board computer for use in this application would provide, at least one channel of mo analog-to-digital conversion capable of a sampling rate of at least 950 samples U31 per second a cigital signal processor or other high-performance microorocessor .3 to control data acquisition sufficient memory for storing centre, 1 software and U32 acquired data and means for presenting the acquired data for analysis. u 4 Another alternative would be to implement signal sampling and conversion 17.-5 means 212 using a single-chiD microcomputer having suitable on-chip facilities u3o for aπalog-to-oigital conversion For example a smgie-chip microcomputer U37 commercially a/ailable from Motorola, Ine Austin TX 85284 under the l s designation " UC68HC11A8 " includes an analog-to-digital converter which can U29 sample and convert up to four analog signal inputs at approximately 15,000 mo samples per second, producing 8- bit results Although higher resolution A-D 11 3 1 conversion results would be desirable.8 bits of A-D resolution may be adequate iij? for some drowsiness detection applications Π J Once EEG signal data is acquired and converted by the signal sampling and U 3 conversion means 212 (Figs.5.14, 15), the data is provided to the drowsiness l as analysis means 136 (Figs.5.14, 15). which analyzes the data and produces iij6 therefrom an output signal on lead 138 which serves as a "measure " of a subject s U37 state of alertness, drowsiness, or sleep The operation of the drowsiness analysis ins means is described in greater detail in connection with Figs.14-16. Any suitable U 39 general-purpose computer, or a special-purpose digital signal processing system, mo may be used to implement the drowsiness analysis means 136 A general purpose mi computer may be used to control EEG signal acquisition means 132 and to U 2 perform all of the other processing functions of signal analysis means 120, ιi43 including those of the drowsiness analysis means 136. threshold means 126. and 1144 artifact detection means 123 ιi 5 If a general-purpose computer is used capabilities required of the computer 114o W vill depend on the sampι,~g rate the number of EEG Signal channels to be IN; nalyzed, and whether th r esults of the analysis are required in an on-line mode ma In some applications sucr as in clinical diagnostics, it may be appropriate to m collect a large batch of Ed 3 data for subsequent anaivsis In other applications, m and particularly in staπo-a.one environments in which a drowsiness detection us.- system is used to provide a warning when a subject performing a task becomes U S3 unable to perform the task safely, it is essential to analyze the data and produce a 1153 real-time result For a smg'e EEG signal sampled at approximately 1000 samples 1154 per second, a medium performance personal computer navingan Intel ' i486 " ιi 55 or equivalent central processor is suitable for implementing the functions of the ιi56 drowsiness analysis means 136 in real time, lower performance computers, such ii5 as a single-chip microcomputer, may also be usable. For a plurality of EEG signals 1153 sampled at higher rates (e.g , 4 signals sampled at 12.000 samples per second), a 1159 mid-range workstation or h'gh-performaπce personal computer may be suitable, 1160 such as personal computer models having Intel Pentium brand processors. 1161 As best seen in Fig. 5. the signal analysis means 120 preferably also 1162 comprises a threshold means 126 responsive to the value of an output measure 1163 138 produced by the EEG signal processor means 124 to produce an output signal 1164 142. This output signal 142 would indicate the level of alertness, drowsiness, or πό5 sleep of the subject relative to the threshold As best seen in Figs.20a, 20b. and 1166 20c (discussed further in greater detail), the output measure 138 is preferably a cont i nuous q uant i ty represent i ng the instantaneous drowsiness of t h e sub j ect Λn exemplary met h od of determining the output measure is d iscussed further i n Ilo9 greater deta i l ,n connect i on w i th F.gs 14-16 but ,t may be brie f l y c h aracter i zed as a we i ghted sum of the inverse of the energy of a su bj ect s EEG s i smal m ir selected frequency bands ir: We have determ i ned exper i mentall y that for any particu l ar su bj ect the 11-3 output measure 138 ,s useful as a relative indicator of drowsiness The measure π -J i ncreases i n magn i tude i n response to the increased drowsiness o f the sub j ect la. Accord i ng to one aspect of the present invention EEG data may oe collected 1176 from an alert i nd i v i dual and used to establish a threshold such that values of u the measure wh i c h exceed the threshold reliably indicate excessive drows i ness 1173 i n that i nd i v i dual As best seen ,n F,g 5 threshold value generator means 144 1179 preferably rece i ves on lead or data path 134 certain baseline EEG data wh i ch 1180 was collected by EEG s i gnal acquisition means 132 and provides a drows i ness 1131 output measure threshold value or signal on lead 146 A process 510 for use i n 1182 connect i on w i th threshold value generator means 144 for determining a sub j ect s 1133 drows i ness threshold from baseline EEG data acquired from the sub j ect ,s shown 1184 i nF.g 19 , and ,s d i scussed further i n greater detail However it may not always 1135 be des i rable or conven i ent to obtain baseline EEG data from the sub j ect It 1136 i s bel.eved that a predetermined universal threshold value which has been 1187 determ i ned exper i mentally may be used with most of the population to prov i de ass acceptable drowsiness detection results Threshold means 126 preferably includes comparison means 140 for U90 comparing the output measure 138 produced by the EEG signal processor means U91 124 with the threshold vaiue on lead 146 to produce a simplified output signal ii92 on lead 142 indicating that the subject has reached or passed a threshold of U93 drowsiness or sleep Comparison means 140 preferably also receives a signal u on lead 150 from artifact detection means 128 indicating that the EEG signal U95 currently being analyzed appears to be contaminated by artifact and one may U96 choose not to use this data in determining the subject ' s drowsiness. ji97 Any appropriate comparison device or process could be used to implement U98 comparison means 140 depending on the form of the signals or data to be U99 compared and on the desired form of the output In an embodiment which ce employs digital signal processing, the values of the output measure, the i oι threshold value, and the artifact detection signal may be stored in locations 1202 which are updated from time to time. In that case, the comparison means 140 103 may be implemented using a software comparison routine. However, in other 120 implementations, the output measure, the threshold signal, and the artifact 1205 detection signal, could appear as individual digital signals, or as three separate 1206 analog signals, and any suitable comparator hardware, such as magnitude 1207 comparators or analog comparators, could be used to implement comparison nos means 140. The output signal on lead 142 may be used to drive an appropriate 1209 display or alarm 122 device to indicate that the subject being monitored by the i 2 io drowsiness detector is excessively drowsy. i2iι As best seen in Fig.5, the signal analysis means 120 preferably comprises i2i2 artifact detection means 128 for determining whether the EEG samples acquired 1213 from the subject are apparently contaminated by artifact Artifact detection 121 means 128 preferably provides an output signal en lead 150 to comparison means ais 140. In this ease comparison means 140 may responsively inhibit drowsiness ι ? io indications produced from artifact-contaminated data This process can i 3 i7 effectively ignore or discard the EEG data acquired during intervals when artifact- 1218 contaminated data is present i2i9 Artifact detection means 128 preferably comprises means 148 for analyzing 122 0 EEG data, anc optionally, information obtained from non-EEG sources which 121 may indicate tnat artifact may be present A process 540 for use in con j unction 122 3 with artifact analysis means 148 for evaluating whether artifact is likely to be 1 present from the actual EEG data acquired from a subject is shown in Fig 18, 1 2 and is discussed further in greater detail The artifact analysis means 148 may 122 be implemented using any suitable hardware and /or software components. 122 6 In the first and second preferred embodiments, which employ digital signal 122 7 processing techniques, and in which sampled EEG data is readily available, the 12 8 artifact analysis means 148 may be implemented using the same general purpose 122 9 computer used to implement the drowsiness analysis means 136 and may operate i 2 3o concurrently with the drowsiness analysis components thereof. i 2 3i However according to an aspect of the present invention, any artifact- 12 3 2 related information which can be conveniently acquired and analyzed may be 123 3 used to assist in determining that artifact is present in the EEG signals. Thus, 12 34 artifact detection may also be performed by analyzing signals (shown generally ι 5 as secondary artifact indicators 824.176; Figs. 5.13a, and 13b) representing 1236 external observations of the subiect s movement or behavior such as signals az- produced by an external motion detector 1238 In some applications it may be unnecessary to explicitly provide artifact i::9 detection means 128 including means for implementing the artifact detection i 2 o method 540 of Fig 18 In general many types cf artifact produce a broad¬ ,2 * 1 band increase in the energy of a subject s EEG signal '- prior-art EEG-based 1242 drowsiness or sleep monitoring systems which rely on increases in low frequency 12 energy as an indication of drowsiness the presence of artifact may be perceived 124 as an indication of drowsiness so that artifact tends to produce a false positive ι?5 alarm In the present invention which relies on a change in energy at higher i24 frequencies as an indication of drowsiness the presence cf artifact produces a i?47 response opposite that of drowsiness and therefore the inventive system provides ι 4s substantial immunity to many types of artifact ι 2 4 Figs 12a and 12b are block diagrams showing the configuration ofa first use embodiment 700 of a drowsiness detection system whicn is constructed according i25i to the present invention and adapted for use in a clinical or research environment us: such as a sleep laboratory The first embodiment 700 is an implementation 1253 based on the general structure of the drowsiness detection system 100 shown in 1254 Fig 5, and employs primarily digital signal processing and related techniques ass Fig 12a is a block diagram showing the configuration of a portion 710 of the i25ό first embodiment directed to the function of interfacing with a subject to be 1257 monitored Fig 12b is a block diagram showing the configuration of a portion ass 712 of the first embodiment directed to acquiring EEG signals and other signals ι 2 59 relating to examination of the s bject processing the data, dispiavingand ι 6o recording results, and performing certain control tasks Figs 12a and 12b are i26i shown separately because the interface portion 710 and processing portion 712 162 may be functionally distinct It may be convenient to locate the processing and 12 63 analysis equipment in a location sufficiently remote from the sub j ect so that 1 264 its routine operation does not disturb or distract the subject However, some 1 265 components of the subject interface portion 710 and the processing portion 712 i 2 6o could be physically integrated if desired 1267 As best seen in Fig 12a. an exemplary application environment for 12 68 the first embodiment 700 may be a clinical or research laboratory in which a 1 2 69 subject ' s alertness, drowsiness or sleep stage is being monitored The application i 2 7o environment shown in Figs 12a and 12b is similar to one which was used in the i 2 i development of the present invention for the detection of the onset of extreme 12 72 drowsiness A brief description of the experimental work may aid in understanding i 2 3 the con figuration of the first embodiment and how the invention could be 12 74 advantageously used in a research or clinical setting The experiments were 12 5 directed to obtaining EEG signal information from a subject while performing an 12 76 assigned vigilance task, and subsequently while resting and falling asleep 12 77 Initially, " baseline " EEG signal information was recorded from a wakeful 12 78 subject over a period of several minutes Next, a performance test was conducted. 1279 The subject was seated in front of a stimulus display and input device The 12 80 stimulus display had several diffuse stimuli which appeared temporarily as dimly i28i illuminated port ions of an otherwise blank display field The subject was provided 1282 with an input device to respond wneπ the stimuli were ooserved A computer 1233 was used to present stimuli at random intervals and the subject was instructed 1234 to activate the input device when each stimulus appeared The stimuli were 735 physically spaced on the display field so that a subject monitoring the stimul i i28o would have to visually scan the display field to observe each stimulus If the 7237 subject failed to respond to a stimulus subsequent stimuli were presented more ass frequently until the subject successfully responded Video cameras , EOG data , s9 and other external indicators were used to determine whether a performance i29o failure was cue to drowsiness or sleep or non-sleep-related behavior Subsequent 121 to the performance testing phase the subject was allowed to rest on a bed and to 1292 fall asleep 123 EEG EOG and other information regarding the sub j ect was recorded 1294 throughout the testing, and an experimental embodiment of the present invention 1295 was used to determine the output measure indicating the sub j ect's drowsiness 1296 Results ofthe trials on human subjects using the experimental embodiment 1297 are discussed further in connection with Figs 7-11 The performance test was 198 intended to simulate a task typical of those in which drowsiness is a particular 1299 problem such as monitoring a sonar or radar screen, or controlling the operation i3oo of a process Such tasks are characterized by limited physical activity and external i3oι stimuli, and requires a high level of vigilance and concentration Successful U02 operation of the experimental embodiment in monitoring the drowsiness of a 1303 subject who participated in the aforementioned performance testing is believed to 1304 be representative of the performance of the invention in a real-world application 13C5 The interface portion 710 may be located in a suitable room ''not shown j lZCo which is preferaeiv isolated from spurious noise lighting and otner subject 130 7 distractions - cnair 730 is pro iαed for tne subject to use while performing the nos assigned tas* - oed 726 is preferably provided for use by the subject when resting ,309 or sleeping as part of the test Appropriate sensors 773 780 ana 782 collect EEG 1310 and EOG signals from the subject and transmit them via a group 732 of signal 1311 leads to an arra y of EEG amplifiers 118 (Fig 12b) (or other suitable amplifiers) 1312 Sensors 778 780 and 782 may be any suitable electrodes transducers or other 1313 sensors, such as conventional EEG and EOG electrodes When tne subject is using 1314 bed 726, sensors 778 780 and 782 may be connected to the group 732 of signal 1315 leads via an alternate connection 728 1316 A display means 742 facing the subject when in a seated position and an 1317 input device 734 are preferably provided for use by the subject during performance 1318 testing The display board 742 and input device 734 are connected to a signal 1319 conditioning unit 768 (Fig 12b) by means of leads 744 and 736 respectively The 1320 signal conditioning unit 768 is connected to a monitoring and control computer 1321 772 using leads 770 The monitoring and control computer 772 is primarily 1322 responsible for interacting with the subject in the context of an experiment or 1323 test, such as the aforementioned subject performance test As best seen in Fig 1324 12b, the examination monitoring and control computer 772 and an acquisition 1325 control and processing means 216a may be provided as separate units However 1326 a single computer could also be used to perform these tasks Display board 742 1327 may be any suitable display which the subject can observe during a test 1 323 In an experimental e bodiment c f the in yen tic ° the c. splay means 742 1 3 2 9 was constructed as a p I a π a r surface ha /mg five spaced indiαa 746 which could i33o be illuminated under control of the monitoring and control computer 772 When 1331 not illuminated the indicia 746 were essentially invisible to the subject The I JJ? input device 734 may be any device act atable by the subject in response to an 1 333 external stimulus including the indicia 746 of display board 742 For example 133 the input device 734 may be a push-button switch Other devices could be used 1335 to implement the display means 742 indicia 7-16 and the input device 734 The i33o signal conditioning unit 76S may be any suitable I/O interface for connecting a I3J7 computer to real-world inputs and outputs For example the signal conditioning ijjs unit 768 may be implemented using a model CIO-DIO24 digital input/output 1 339 module coupled using suitable interconnection means to model DCO5-B DC ijo output isolation modules and model DC15-B DC input isolation modules, all of i 4i which are commercially available from Omega Engineering, Iπc , P O Box 4047, 142 Stamford, CT, 06907 Other interface equipment could also be used 1 343 In order to observe and record the subject s movements suitable camera 134 means 714 and 722 may be trained on the chair 730 and the bed 726 respectively 135 Camera means 714 and 722 may be any suitable cameras providing a recordable 136 output; in a preferred embodiment of the invention, camera means 714 and 722 1347 are CCD video cameras having high sensitivity to visible and infra-red light, and 133 provide a standard NTSC video output on leads 716 ano 724, respectively Leads 1349 716 and 724 are preferably connected to a suitable monitor 750 for viewing by a jj-50 remote person, such as a nurse or a person supervising the experiment or test. An 357 intercom 7^0 in the testing room may be connected to the monitor 750 via lead 7352 748 to allow the subject and the remote person to communicate crallv. 1353 The monitor 750 preferably includes means to allow an operator to 135 select the camera from which an image is displayed, and provides a video signal ss corresponding to the signal from selected camera on lead 752 to a recorder 754. 1356 Recorder 754 may be any suitable image and sound recording means such as a 135 conventional VHS-format " video ' cassette recorder The monitor 750 provides a ass suitable audio output signal on lead 760 to an audio mixer pane! 762 The output 1359 of the audio mixer panel is provided to recorder 754 A monitor 758, which may nee be a conventional television receiver or monitor is connected to the recorder to 1361 present the recorded video and audio to a user An infra-red light source 718 1362 supplied with power on lead 720 illuminates the testing room to enable operation 1363 of the cameras even when room lighting is subdued or off 1364 The image information supplied by the cameras 714 and 722 may be IJ65 useful in the evaluation of the cause of some artifacts in the EEG record, and to 1366 determine whether a failure by the subject to respond to a stimulus was due to 7367 some cause other than drowsiness or sleep. For example, the subject may be fully 1368 alert, but may not be looking at the display board 742. Accordingly, information 1369 obtained from the cameras may be supplied to the drowsiness detection system i37o 100 as an optional secondary input to the artifact detection means 148 (see lead 132 176, Fig.5). Systems (not shown) which detect motion or change in a video 1372 signal are known in the art and could also be used to produce a signal indicating 1373 the presence of movement at particular times in the subject s EEG record. Such 1374 a signal could be used by art i fact detection means 148 (Fig 5) either alone or in 1375 conjunction with other artifact indicators (such as tnose derived from the EEG 1376 data i tself and from EOG signals) for artifact detection as previously discussed 1377 As best seen i n F i g 12b , EEG and related signals may be provided to an 13 78 array of EEG ampl i fiers 118 on signal leads or bus 732 As described prev i ously 1379 in connect i on w i th F i g 5 any suitable EEG amplifiers may be used Suitable use anti-aliasing filters 164 are preferably provided for each channel of EEG or related i38i data , and may either be integrated with the EEG amplifiers (as depicted i n Fig 2JS2 12b), provided as self-contained filters, or provided as a component of the EEG 1383 signal acquisition means 132a A suitable anti-aliasing filter 164a for a s i ngle 133 channel has been described previously in connection with Fig 6a , but other f i lter ass designs could also be used EEG amplifier 118 may be connected to suitable 1386 test equipment, such as oscilloscope 786 and signal generator 788 for use by the 138 7 testing supervisor or other personnel in verifying the operation of the amplifier and use in calibration. 1389 Amplified EEG and related signals are provided on signal leads or bus 776 to i39o a mon i tor i ng and control computer 772 and a acquisition control and process i ng i39i means 216a. Each computer preferably cooperates with suitable means for 392 interacting with a testing supervisor, such as displays 790, 792 , respectively The 13 93 acquisition control and processing means 216a samples EEG signals , converts the 139 4 sampled signals into representative digital data, and stores the digital signals so 1395 that they may be analyzed in real-time or in subsequent processing An additional 1396 computer system 218a may be provided to analyze the digitized EEG signals 1397 and present results Analysis control and processing means 218a preferably also 1398 cooperates wιt n suitable means for interacting with a testing supervisor, such as 1399 display 798 i4oo Acquisition control and processing means 216a provides the digital data 7402 representing the EEG signals to the signal analysis computer 218a on a suitable 140 computer-to-comouter interconnection medium 794 The interconnection 403 medium 794 may be implemented using any suitable medium, such as a serial data 1404 link or a local area network If it is desired that the signal analysis computer 218a 1405 analyze EEG data in real time, interconnection medium 794 must be capable of 2406 relatively high-s p eed data transfers For example a 10 Mbps Ethernet network, or 70 another network providing suitable performance could be used nos As best seen in Fig 12b, the first preferred embodiment 700 may comprise 1 409 separate computers (control and processing means) 772, 216a, and 218a for i4io examination control and monitoring, data acquisition, and signal analysis i4u functions However, all of these functions could also be implemented using a 2422 single computer in the configuration of Fig 12b computer 216a is used to i4i3 implement the functions of the EEG signal acquisition means 212 of Fig. 5, and i i4 computer 218a is used to implement the functions of the drowsiness analysis i4i5 means 136 of Fig 5 Computer 218a may also be used to implement the other i4i6 functions of signal analysis means 120, including those of threshold means i4i7 126, and artifact detection means 128. The drowsiness analysis means 136a, i4i8 is described further in greater detail in connection with Figs. 14, 15, and 16. i i9 In the first preferred embodiment 700, in which several (e.g., 4-8) channels of i 2 o EEG data may be simultaneouslyacαuired at a sampling rate of 950 samples ,4?ι per second per channel medium performance personal computers having Intel 122 "i486 or equivalent central processors are suitable for use as computers 772 and 123 216a and 218a however other computers could also be used and the minimum .2 performance required of computers 772 216a and 218a has not been determined 1 5 The inventive drowsiness detection system may analyze acquired data and 1420 produce a result on-line in a real-time mode or may perform these functions 4?7 off-line In a clinical or research laboratory environment it is typically desirable 1423 to store a lasting record of the EEG and related signals accuired from a sub j ect 1429 for further review and analysis even if analysis is initially performed on-line If i 3o sampled EEG data is to be stored for later review large amounts of storage may i 3i be required in computer 216a If the data is not compressed acouisition of 8 1 4 2 channels of EEG data at 95012-bit samples per second per channel requires 1433 approximately 534 megabytes of storage per hour If the sampled EEG data is 1434 not stored for later review only a small amount of data storage is required As 35 discussed further in greater detail, the drowsiness analysis means 136a operates 1 436 on successive, overlapping segments of data or " windows representing the 137 subject s EEG signal during short time intervals As a minimum, storage must 1 438 be provided in either the data acquisition means 132a (computer 216a) or in 1 439 drowsiness analysis means 136a (computer 218a) for at least one complete 1440 "current" window of data, plus approximately one half wind o of new data which 1 44 2 may arrive during or after analysis of the " current" window Additional storage 1 442 may be required for other functions, such as threshold determination and artifact 143 detection 1444 A synchronization device 764 is provided to allow correlation of the 145 contents of the video and audio record made by recording means 754 with subject i4o performance test results recorded by computer 772 and EEG data acquired 2447 and recorded by computer 216a Synchronization device 764 may be activated we automatically (by monitoring and control computer 772 via signal conditioning 2 9 unit 768) to supply a synchronization signal on leads 766 to recording means 754 i45o via audio mixer panel 762 A functionally equivalent signal may be supplied to 2451 computer 216a The synchronization signal allows a particular instant in time 45? to be identified in the separately acquired subject performance video and audio, 1453 and EEG records The synchronization signal may be any suitable signal usable 15 and recording means 754 For example, synchronization signal may be an audio- 255 frequency signal generated by a conventional oscillator or signai generator. Other 2456 synchronization signals and methods, such as time stamps derived from a master 2457 clock, could also be used 1458 Fi s. 13a and 13bareblockdιagramsshowιngtheconfιguratιon ofa second 1459 embodiment 800 of a drowsiness detection system which is constructed according 2460 to the present invention and adapted for use in a real-world alertness/drowsiness 46 monitoring and alarm application For example, the second embodiment might be 262 used to monitor the drowsiness ofa person assigned to the operating position of 1463 a sonar or radar system, a power or chemical process plant, or a vehicle In order 2 o to be most useful in such environments, the second embodiment 800 is preferably 2465 adapted for implementation in a relatively small, portable, configuration 840 i 6o which mav be assembled into a single enclosure free of external components 16 7 except for sensors le g electrodes) for EEG and related signals 1468 The secon d embodiment 800 is an implementation base d on the General 2 69 structure of the drows i ness detection system 100 shown in Fig 5 and employs i o pr i mar i l y d i g i tal s i gnal process i ng and related technioues Fig 13a is a block i 47 i diagram snowing tne configuration ofa portion 810 of the secon d embod i ment i 7i d i rected to the funct i on of i nterfac i ng w i th a subject to be monitored Fι * 13b 2 4 73 is a block diagram showing the configuration ofa portion 812 of the second 1 embodiment directed to acquiring EEG signals and other signals relating to 1 5 monitoring of the subject processing the data displaying and recording results , 2476 and performing certain control tasks Figs 13a and 13b are shown separately 1477 because the interface portion 810 and processing portion 812 may be functionally 1 78 d i st i nct Preferably process i ng portion 812 is integrated in a small equipment 24 package 840 which may be carried on the subject s person or installed in an i48o operator s console or in a vehicle, at a location near the sub j ect The processing 2 4 3 portion 812 may be conceptually divided into an analog signal processing section 1482 842 and a control and processing means 216b which may be conceptually viewed 4 S as a digital signal processing section The conceptual division of processing 1484 portion 812 into two functional sections is useful in describing its operation , 1485 but when implemented the processing portion 812 need not have separate i486 corresponding physical components and could be implemented using a single 1487 pπnted-circuit module or even a single integrated circuit 188 As best seen i n Fig 13a an exemplary application for the second 2439 embodiment 800 may be to monitor for drowsiness a subject 110 assigned to i o the operating position 866 of a sonar system snown schematically as 818 One i49i of skill in the a r t would ap p reciate now the second embodiment 800 could also s? be used in man/ other applications with minor or no modification The sonar 193 system 818 includes a display 820 and operating controls 822 with which the 149 subject is expected to interact Display 820 is driven bv signals on lead 868 from 195 sonar processing equipment (not snown) The ooeratmg controls 822 produce i4 signals on lead 870 for use bv the sonar processing equipment A signal on lead 79 824 may be oeπ » ed from one or both of the signals on leads 868 and 870 for use 1 98 as a task-related indication of subject activity and may be supplied to artifact 1499 detection means 128 (Figs 13b and 5) (discussed further in greater detail) Other i5o means 872 such as a video camera or a motion detector may be used to provide a i5oι secondary subject activity signal to artifact detection means 128 on lead 176 1 502 At least one EEG signal must be collected from the subject Additional 2503 EEG signals, and other physiologically relevant signals such as EOG signals, 50 may also be collected and used In the stand-alone environment of the second 1505 embodiment 800 because of equipment size, weight cost, and operator mobility 2506 constraints, it may be feasible to collect and process only a few channels, or only 1507 a single channel of EEG data When processed and analyzed according to the 503 present invention a single channel of EEG data is sufficient to reliably detect 250 the onset of extreme drowsiness in a human subject Accordingly the second aw embodiment 800 will be described hereafter as acquiring and processing a single 2522 EEG signal channel However, one skilled in the art will appreciate that the second 1512 embodiment 800 could be extended to collect and process multiple EEG signal 1512 channels as previously described 1514 Appropriate sensors 814 collect an EEG signal from the subject ana 1515 transmit it on leads 816 to a suitable amplifier 118b Sensors 814 may oe any 151o suitable electrodes transducers or other sensors such as conventional EEG and 151 " EOG electrodes Amplιfιe r 118b provides the amplified EEG signal to a suitable 152; anti-aliasing filter 164 Tne amplified and filtered EEG signal is then provided 2529 to data acquisition conversion and storage means 212 Amplifier 118b may 1520 be implemented using a special purpose portable EEG amplifier Alternatively 1521 amplifier 118b may be im p lemented using any suitable operational amplifiers 1522 arranged to provide a high input impedance and the desired amplification 1523 which may be in the range of 10 3 -10 5 Alternatively amplifier 118b may be 1524 implemented using a high- gam amplifier module A variety of appropriate 1525 operational amplifiers and high-gain amplifier modules are commercially available 1526 and may be used in the second embodiment 800 1527 Due to differences among subjects and subject-electrode connections the 1528 amplitude range of the EEG signal acquired from a subject is only predictable to 1529 within a range Because of this variability in order that the digitized EEG signal 1530 data has sufficient resolution for analysis amplifier 118b preferably provides a 1531 user-adjustable gain setting Incorrect gain settings may cause high-amplitude 1532 signals to be " clipped " or low-amplitude signals to occupy too small a portion 1533 of the analog-to-digital converter data range, resulting in reduced signal quality 1534 Amplifier 118b may provide a user-adjustable gam setting which may have a 1535 plurality of selectable discrete gam settings or a continuous range of stain settings 25J However, in the stand-alone environment to which second embo d iment 800 i s 2537 directed, it may be inconvenient to perform such adjustments in the field , as 1533 might be required to accommodate multiple subjects at various working t i mes or 2539 shifts. 50 A suitable automatic gam control (AGC) (not shown) may be prov i ded 151 in amplifier 118b to ensure that the amplified output signal remains with i n the 25 dynamic range of downstream data acquisition, conversion, and storage means 254 212b. The automatic gain control (AGC) preferably provides an output signal 1544 indicating the current gain setting which may be recorded with the associated 255 EEG data: when the digitized samples are analyzed, they may be normalized by i54o multiplying them by a corresponding scaling factor so that the numerical value of 257 all samples will have a consistent physical meaning regardless of the amplifier gain 1548 used to acquire them. 254 An alternative to providing an AGC or manual gain adjustment on the 2550 amplifier 118b is to select a fixed amplifier gain such that the maximum s i gnal i55i amplitude never exceeds the full-scale value of the data acquisition means 212b 1552 and to provide a greater resolution for the analog-to-digital converter in the 1553 data acquisition means 212b. If the analog-to-digital converter has sufficient 2554 resolution, even signals at the low end of the range of expected signal amplitudes ass will be represented with sufficient resolution that the drowsiness analysis means i55ό can operate effectively. The resolution required in excess of that needed for the 1557 analysis means to operate with full-scale signals depends on the difference in 2553 amplitudes between the largest and smallest expected EEG signals 2559 Filter means 164 receives the amplified EEG signal from amplifier 118b i5oo Any suitable anti-aliasing filter designed to provide sufficient attenuation and a i5oi desired cut-off freouency selected according to the frequency band of interest and 2562 the sampling rate of data acquisition, conversion, and storage means 212b, may 2563 be used to implement filter means 164 The sixth-order Butterworth switched- 156 capacitor filter 164 of Figs 6a and 6b may be a suitable choice 265 Data acquisition, conversion, and storage means 212b implements the 5 o functions of signal sampling and conversion The data acquisition means 156 212b receives the amplified and filtered EEG signal from filter means 164, and i5os samples the data converts the data to digital samples, and stores a sufficient 25 quantity of the resulting samples for use by drowsiness analysis means 136b. As 250 discussed below in greater detail, the drowsiness analysis means 136b operates on 52 successive, overlapping segments of data or " windows " representing the sub j ect ' s 157 EEG signal during corresponding time intervals For example, the window time 1573 interval might be approximately 2 seconds, and the drowsiness analysis means 1574 136b may be configured to determine the output measure approximately once per 157 second, using approximately one second ' s worth of new data, and one second's 1576 worth of "old " data (i e overlapping a portion of the data used in the previous 2577 determination). The data acquisition means 212b must provide sufficient storage 1578 for the data required by the drowsiness analysis means 136b. For example, if the 579 drowsiness analysis means 136b operates on successive overlapping segments of 2530 data, then the data acquisition means 852 should provide sufficient storage for at 15S2 least one complete " current " window plus approximately one half window of new 253 data which may arrive during or after analysis of the current ' window Storage 253? may be reouired for other functions such as tnreshold determination and artifact 158 detection, and would certainly be required if a permanent record of the acquired 535 EEG signals is desired 1586 The EEG signal data acquired by the data acouisition means 212b is 158 provided via leads 134 to drowsiness analysis means 136b, threshold value 2533 generator means 144b, and artifact detection means 128b on output leads 134, 2539 which operate as discussed previously in connection with Fig 5 Drowsiness ago analysis means 136b will be discussed in greater detail in connection with Figs. 151 15 and 16 A process 510 for use in conjunction with threshold value generator 2592 means 144b for determining an appropriate drowsiness detection threshold for 159 an individual subject based on EEG signal data acquired from that subject is 254 shown in Fig. 19, and is discussed below in greater detail A process 540 for use 1595 in conjunction with artifact detection means 128b for evaluating whether artifact 259ό is likely to be present from the actual EEG data acquired from a subject is shown 1597 in Fig. 18, and is discussed below in greater detail 2595 Drowsiness analysis means provides an output signal on lead 138 which 1599 functions as a continuous measure of the subject ' s state of alertness, drowsiness 26oo or sleep. The outputs from drowsiness analysis means 136b, threshold value 602 generator means 144, and artifact detection means 128b are provided to 1602 comparison means 140 which compares the output measure with the threshold 260 value 146 to produce a simplified output signal indicating that the subject has 1 604 reached or passed a threshold stage of drowsiness or sleep Suitable amplifiers, 1605 line drivers, or relays (not shown) may be provided to convert the output signal 606 into a visual alarm output on lead 832 and an audible aiarm output on lead 836 for 1 607 driving a visual alarm 830 and an audible alarm 834 respectively Alarms 830 and 1 608 834 advantagecusly provide a warning to the subject cr to others in the vicinity 2609 that the drowsiness detection system 100 has determined that the subject is low becoming excessively drowsy before the subject becomes unable to perform their 2 6 2 2 assigned task as intended. The subject cr others observing or hearing the alarms i6i2 can then take appropriate action. i6i Any suitable digital signal processing means may be used to implement 1 6 2 4 the functions of the control and processing means 216b (including drowsiness io analysis means 136b, threshold value generator means 144b, artifact detection i6i6 means 128b, and comparison means 140b) For example, as best seen in Figs. i6i7 13b and 15, the control and processing means 216b digital signal processing 1 618 section may employ a general-purpose central processing unit (CPU) 846 i6i9 incorporating program and data storage means, input/output means, and 1 6 2 0 associated support components The CPU 846 may be implemented using a 2 622 single-board digital signal processing computer, a single-board control computer, 16 2 2 or even a high-performance single-chip microcomputer. The central processing 1 6 2 3 unit 846 preferably controls the data acquisition means 212b, which may be 1 6 2 4 incorporated with CPU 846 on the same card or integrated circuit. A single-board 1 6 2 5 control computer having an Intel " i486" or equivalent central processor, and 1 66 approximately 1 megabyte of storage for program code and data (RAM, ROM, 162 EEPROM, or equivalent) is suitable for use as CPU 846 Such computers are 1623 commercially available from numerous sources 26?9 Fig 14 is a block diagram showing a first embodiment of an EEG iό3o signal processor means 12^ including EEG Signal acquisition means 132a. lόJi and drowsiness analysis means 136a and certain preceding signal acquisition 1633 components, configured for use in conjunction with the first embodiment 700 1633 of a drowsiness detection system 100 constructed according to the present 2634 invention Fig 15 is a block diagram showing a second embodiment of an 2635 EEG signal processor means 124 including signal acquisition means 132b 1636 and drowsiness analysis means 136b, and certain preceding signal acquisition 163 components, configured for use in conjunction with the second embodiment 633 800 of a drowsiness detection system 100 constructed according to the present 1639 invention Fig. 16 is a data-flow diagram depicting, in simplified form, the 1640 processing of EEG signal information acquired from a subject, as that information i64i is operated upon by the processing means of the embodiments of Figs.5 and 14- 6? 15. Interconnections between the components of EEG signal processor means 124 264 are designated herein as signal leads, or in the case of digital data, " data paths." 2644 For each such signal lead or data path in Figs. 14-15, there is a corresponding 2645 signal or data element shown in Fig 16 The elements of Fig 16 represent the i6ό form of the signal ordata being transferred over the corresponding signal lead or 267 data path. 163 The principal difference between the two embodiments of EEG signal 2649 processor means 124 shown in Figs 14 and 15 is that in the embodiment of Fig. 2650 14 separate control and process i ng means 216a 218a are provided for the data 1652 acqu i s i t i on and analys i s respect i vely while in the emoociment of Fig 15 a single 2652 control an d process i ng means 216b serves Doth functions Accordingly , Fι<*s 14- 265J 16 w i ll be d i scussed together and unless otherwise specified the embodiments of 65 Figs 14 and 15 should be considered to operate similarly ; °" As best seen i n F i gs 14 and 15 and discussed in detail previously at least 265o one channel of EEG signal information is acouired from the subject 110 us i n^ 165 7 electrodes 112 and 114 or other suitable sensors and is supplied to amplif i er 118 loss on leads 116 As best seen in Fig 16 the signal on leads 116 is an analog , time- 2 59 doma i n s i gnal X(t) Ampl i f i er 118 (Figs 14 and 15) provides an amplified version 26o o of the s i gnal on lead 130 (Figs 14, 15 and 16) Anti-aliasing filter 164 (Figs 14 i66i and 15) supplies a low-pass-filtered version of the amplified EEG signal on lead 166 310 (F i gs 14-16) Accord i ng to an aspect of the present invention, the cut-off 26 frequency and other characteristics of the filter 164 are selected such that the 266 filtered output signal includes frequency components substantially above 30 Hz 665 According to traditional EEG doctrine signal components above 30 Hz have been 2666 considered unusable or noise ' and have therefore been discarded or ignored 66 Analog-to-digital conversion means 212 (Figs 14 and 15) samples the 1668 filtered and amplified EEG signal at a selected sampling rate and converts each 1669 sample to a word " or un i t of digital data As discussed above, the sampling i67o rate f s must be at least twice the maximum frequency f ma of the signal being i67i converted For example if the maximum frequency f max is 475 Hz , then the 1672 sampling rate f s must be at least 950 samples per second The resolution of the 167 analog-to-digital conversion means 212— hat is the length, in bits, of the unit ό74 of digital data produced by the ana log-to- digital conversion means 212 for each 5 sample-determines the smallest difference between analog samples which can be i67ό represented by distinct output codes 1 677 It is believed that, in the current state of the A-D converter art, word 1678 lengths in the 12- bit range provide adequate precision at moderate cost. A iό79 shorter word length might provide insufficient precision but could be usable if i68o predominant signal components (such as those from approximately 0-30 Hz lost and a power line noise spike at 50 or 60 Hz), which mav not be of interest in 1632 the inventive drowsiness detection system 100, are substantially attenuated by 1683 appropriate filtering in the analog domain before the EEG signal is presented to 1 6 3 the analog-to-digital conversion means 1685 As best seen in Figs. 14-15, control and processing means 216 is connected 1 686 to analog-to-digital conversion means 212 and to storage means 214 Control and 2637 processing means 216 controls the A-D conversion process, routes the resulting 1 638 digital data corresponding to converted digital samples to the storage means 214, 1 689 and forwards the sampled data to the drowsiness analysis means 136. 760 As best seen in Fig. 14, in a first embodiment 700 of the invention, two i69i separate control and processing means 216a and 218a may be provided to serve 1692 the signal acquisition means 132a and drowsiness analysis means 136 respectively. 69J As best seen in Fig. 15, in a second embodiment 800 of the invention, a single 1694 control and processing means 216b may be provided to serve both signal 2 695 acquisition means 132b and drowsiness analysis means 136b. Each of the control 2696 and processing means 216a 218a. and 216b could provide additional functions 2 = 7 as required in the application, such as artifact detection and thresholding. In 2693 the first and second preferred embodiments 700, 800 of the invention, control 2 69 and processing means 216 and 218 are preferably implemented using a suitable ace programmable computer having a central processor, such as a commerαally- i7oι available general-purpose microprocessor or digital signal processor, coupled 102 to appropriate storage means 214 and input/output support components (not 1703 shown), or a single-chip microcomputer with storage means 214 and input/output no support disposed internally or externally For example, a personal computer or 2 05 single-board computer having an Intel " i486 " or equivalent central processor 2706 could be used. 1 707 If separate control and processing means 216 and 218 are provided for data 1708 acquisition and analysis, these two systems are preferably interconnected using 20 any suitable computer interconnection medium 794 In the preferred embodiment 270 800, an Ethernet local area network is used as the interconnection medium, but 2722 others, such as a high-speed serial data link, or an IEEE-488-type instrumentation 2722 bus, could also be used. 723 The data from analog-to-digital conversion means 212 corresponding to i7i4 the sampled EEG signal is supplied via data path (or leads) 314 (Figs. 14-16) to 225 storage means 214. As discussed previously, the drowsiness analysis means 136 i7i operates on successive, overlapping segments of data or " windows " representing aa the subject ' s EEG signal during corresponding time intervals. Storage means 273 214 should provide sufficient storage for these operations. Storage means 214 i7i9 may be implemented using any suitable read-write storage or memory For i o example, storage means 214 may be implemented using semiconductor RAM" 2722 " EEPROM, " or " FLASH ' memory rotating magnetic media, or any other 1 7 22 appropriate storage 17 In the first and second preferred embodiments 700, 800, the functions of 272 drowsiness analysis means 136 are preferably imolemented using known digital 2725 signal processing techniques Thus, as best seen in Fig 14, in the first preferred 1726 embodiment 700, all of the functional blocks or components within the confines 272? of drowsiness analysis means 136. may be implemented using appropriate software 178 components running in analysis control and processing means 218a As best seen 129 in Fig 15, in the second preferred embodiment 800, all of the functional blocks or i 3o ' components of drowsiness analysis means 136 are preferably implemented using 2732 appropriate software components running in the control and processing means 1732 216a. However, each functional block of drowsiness analysis means 136 could also 1733 be implemented using appropriate digital signal processing hardware components. 274 Interconnections between the components of drowsiness analysis means 135 136 are designated herein as " data paths " In an embodiment in which the 2 6 functions of drowsiness analysis means 136 are implemented using software 177 running in a programmable general-purpose digital computer, the data paths 2 JS discussed herein will typically be regions in memory which are shared by or 2 9 accessible to the relevant software components In an embodiment in which the 1740 functions of drowsiness analysis means 136 are implemented using digital signal i74i processing hardware, the data paths discussed herein may be any suitable data 142 interconnection between the components such as a serial data link, a data bus 274? or shared memory 274 Stored EEG signal data is supplied from storage means 214 to data 75 selection means 220 via a data path 316 The digitized EEG data is stored as a 4o plurality of EEG signal samples corresponding the subject s EEG signal acquired 1747 over a time interval of interest, and may be arranged as a contiguous group of the 273 most recently acquired samples Fig 7 is a graph 152 showing the time evolution 179 of an exemplary EEG signal 154 which was acquired from a human subject usιπ°- 1750 an experimental embodiment similar to the first preferred embodiment 700 2752 disclosed herein The graph 152 was plotted from data which had been sampled 2752 as disclosed herein and stored in a form corresponding to data element 316 of 2753 Fig 16 (see below) The de p endent variable in Fig 7 is the amplitude of the EEG 1754 signal in microvolts The sampling rate at which this signal was acquired was ass 950 samples per second. Data corresponding to 2048 samples are plotted. The 756 independent variable in Fig 7 is an ordinal number corresponding to the position 275 of each data point in the sequence in which the data wasacquired Because 2753 the sampling rate was constant, this sample number is equivalent to time. The 1759 resulting time interval represented in the graph is approximately 2.16 seconds. A i76o first-order Butterworth filter having a half-amplitude frequency of approximately i76i 100 Hz was used as an anti-aliasing filter in the experimental embodiment. 1762 Although acceptable results were achieved using this filter, a higher performance 2763 filter, such as filter 610 (Fig. 6a), may provide improved system performance a-s 276 it does provide improved anti-aliasing performance. Comparing graph 152 with 75 graph 92 (see Fig 3) which shows an equivalent signal 94 acquired according i7όo to prior-art methods it is apparent that significant high-frequency information is 1 6 present in signal 154 due to acouisition according to the present invention i7θ A sufficient αuantity of data must be available for processing by drowsiness 1 769 analysis means 136 In a drowsiness analysis means 136 that operates on i77o successive data batches comprising overlapping segments of data or windows ' i77i representing the subject s EEG signal during a brief ccrresponding time interval 1 77 2 the window defines the quantity of EEG signal data on which the drowsiness 1 775 analysis means 136 operates at a particular time Each successive operation of 1 77 drowsiness analysis means 136 may operate on a segment of newly obtained 1 77 5 data and a segment of previously analyzed data so that each analysis window 17 7 6 overlaps the previous analysis window by a selected amount In the first and 1 777 second preferred embodiments 700 800 of the invention, each window overlaps 1 77 8 the previous window by one halfof the width ofthe analysis window, but other 1 779 amounts of overlap could also be used 730 The size and overlap of analysis windows may be selected according to 1781 application requirements The window size should be large enough to provide 1782 sufficient resolution for downstream analysis (discussed further in greater detail), 17 83 but small enough that the analysis can be computed acceptably quickly using 734 affordable computing resources The window overlap should also be selected to 2 S 5 allow the output measure to be frεαuently determined, so that adequate warning 2 7 3 6 is provided if a subject becomes excessively drowsy 1 37 Data selection means 220 extracts a quantity of data corresponding to ass one analysis window, and presents it to windowing 2 " 22 ~ on data 172) path 320 The process of extracting tne cata may be re: ' erred to as " rectanzular az windowing ' The extraction, cr selection process is ecu: /aient to multiplying ail . - 9. of the data by a function navmg the value 1 during the window period and zero at all other times The windowing function means 222 multiplies the data in 1737 the selected window by a window function which reduces the a plitude of the samples in the neig borhood cr the begin nmg and end of the win ow 179: A number of windowing functions such as recta ne iar. Hamming 175o Hannmg, and Blackman. are known in the signal processing art In the preferred 179 " embodiments 700.800. α Hannmg window function, which results in consicierably 2 3 smaller side-lobes than the rectangular function (at the expense of a broader I .-?.; main lobe), may be used If analysis warrants, other windowin g functions may 1300 be employed (such as the Dolph-Chebyshev function) to provide improved i3oi characteristics over the Hannmg function at the expense of some increased 2 307 complexity. One subsequent processing step performed by drowsiness analysis 1303 means 136 is spectral or frequency analysis. The data window type (rectangular. 30 Hannmg, Hamming, or Blackman, etc.) influences the resolution of the 1305 downstream frequency analysis by " blurring " the results in the frequency domain. 1306 The output of the windowing function means 222 is provided to frequency 130 analysis means 224 via data path 322 as a vector of time-domain digital samples. 2303 The length of the vector (i.e.. the number of samples) is equivalent to the window 1309 width. aa Frequency analysis means 224 employs any suitable frequency or spectral is.-.- analysis technique to determine tne amplitude and phase of components cf the aa sampled EEG signal at various freaueπcies thrcugncu: a defined frequency raπ^ aa cf interest. For example in a preferred err oodiner.t cf the invention tne Fast 15 4 Fourier Transform (FFT) may be used i3is The FFT is a frequency analysis technique which is -well known in the signal 16 processing art, and which converts a vector cf TJ regularly-spaced time-domain isi • samples of a signal into a vector of data representing tne signai in the frequency 113 domain Each element of the FFT vector is a complex number representing the i i amplitude and phase of the signai component at one cf the IM regularly-spaced 1 30 frequencies The FFT is a good candidate for use in the preferred embodiments i3?ι of the present invention because its properties are well understood and there is?? exist high-speed implementations of the algorithm which are commercially and i3 2 3 otherwise available in software and hardware Other suitable frequency analysis 132 techniques could also be used to implement frequency analysis means 224 i325 The output of the frequency analysis means 224 is provided to power 132 6 spectral density determination means 226 via data path 324 The frequency 132 analysis output 324 is an output vector in which each element is a complex 1 823 number representing the amplitude and phase of the signal component at one of 32 the N regularly-spaced frequencies. l3 3o The power spectral density determination means 226 determines the 232 power present in the sampled EEG signal at each frequency point by squaring the 1332 magnitude of each component of the FFT output vector. The use of frequency or 1333 spectral analysis allows the isolation of EEG signals in specific frequency ranges 1334 wr hich carry information related to alertness drowsiness and sleep Proper selection cf anti-aliasing filter characteristics and sampling rates is essent i al in 1 330 order to preserve signal components in these frequency ranges In previously is ? - discussed tnals using an experimental embodiment of the invention a sampling z 3 rate of 950 samples per second with a suitable anti-aliasing filter was used to 2339 capture signal components up to 475 Hz Frequency components in the ranσ c from 80 to 475 Hz were found particularly useful in detecting the onset of extreme 2 3 1 drowsiness , and provided improved drowsiness detection performance compared 34? to that reported for prior-art systems However it is believed that otner frequency 13 4 ? ranges above 30 Hz may also produce useful results and perhaps may provide 1344 even better drowsiness detection performance 'θ 5 Fig 8 is a graph 156 of the power spectral density 158 from 0-123 Hz , 36 determ i ned from the sampled signal 154 shown in Fig 7 Almost all of the 24 7 energy in this signal is below 30 Hz, despite the relatively high sampling rate 133 and an appropriately selected anti-aliasing filter which were used to acquire the 23 data from which Figs 7-8 were produced A tiny spike appears at 60 Hz , and 2350 represents signals attributable to commercial power line noise acquired with the last signal. Fig. 9 is a graph 160 of the power spectral density 162 from 31-475 Hz. 2352 determined from the same sampled signal 154 used to produce graph 156 of F 35 8. Elimination of the frequencies from 0-30 Hz allows the graph to be re-scaled so 1354 that the power present at frequencies above 30 Hz is now visible 1SSS Figs. 10-11 are graphs depicting two power spectrums in the frequency 2356 range 100-475 Hz, derived from EEG signal data obtained from a human subject 1357 using an exper i mental embod i ment of the present invention Table 2. represents 1353 i n numer i cally analyzed form relevant aspects cf the cαtα sro-wπ graphically 1 5? m Figs 10-11 -- comparison of the data in Table 2. arc the erapπs cf Fισ S 15 0 10-11 , w i ll prov i de i ns i ght into the drowsiness-related information contained i n ae i higher-frequency components of a subject ' s EEG signal In Fig 10. the pcwer 1362 spectrum corresponds to EEG data obtained during hits ' - that is intervals i3o dur i ng wh i ch the subject successfully responded to a visual stimulus presented 1364 accord i ng to the exper i ment In Fig 11, the power spectrum corresponds to EEG 1305 data obtained dur i ng misses " - that is, intervals during which the subject fa i led 1366 to respond to a v i sual stimulus Table 2 shows, for several selected frequency ranges and subranges the amount of energy m the EEG signals attributable to the 130 frequencies i n each range, for hits and misses The " Change in Energy" column 230 presents a factor for each selected frequency range indicating the relative amount 1370 of energy in that range for EEG signals acquired during "hit" events , compared to 137! the energy in the range for signals acquired during " miss " events If the entry for 1372 a frequency range in the " Change in Energy column is positive then it represents 13 the quotient formed by the energy in the " miss " column divided by the energy in 2374 the " hits" column. If the entry for a frequency range in the "Change in Energy 35 column is negative , then it represents the quotient formed by energy in the "hits" 2376 column divided by the energy in the "miss " column 237 ' E I PΞSUL " ;Ξ " S U; 137? Frequency -nct,c r • rjci.cr, '_π jnse 135 Ranse Energy z- Iota; ■r ' Ts'.ji ir 1 fcπersy 13S1 delta 105-1 200ie-G2 -> 545\ i 7 Ξ - - •! I 41601 -ffJβT 1332 thet J i 371 44 ό -01 ;) rj : 1 15rr-.r 2 Ξ 51 N -255 ; 133 alprκι i'3 12, 9370^-01 _2 ό5N - 13^-1 " . :• 180% 2385 1334 beta C13 301 2303 -01 5253°o 224 -C1 - 72% .1028 335 0 to 30h ' z 4083e — 02 1-10.91% J <~-~ ^ U? <IOI- , n , 336 0 to 40Hz 4135e-02 -01Q7° U J 46ve-G2 99141% -103' 133 0 to 475hz 4293---02 100000 , 150,:- • • _ '.'30000% ■ 1 ().; 33 31 to l00Hz 1547e-01 73521° 0 90 ^-2^ 91672 2593 2339 100 to 200Hz 4239e-r00 20140% - ύl- Ol 7224% -9017 180 100 to 300Hz 5.249e÷00 24946% 5225e-01 8029% -10046 2391 100 to475Hz 5.571e-00 26479% 5419e-01 8328% -10232 1392 200 to 475Hz 1333e τ 00 3333% 7 180e-02 1103% -18559 1393 300 to 475Hz 3.225e-0i 1533% 194 -02 0298% -16613 13 31 to475Hz 2104e-01 100000% 5507e-00 100000% -3234 i89s To produce the graph of Fig. 10, the EEG data corresponding to 25 hits was 1896 selected according to the following criteria: 1397 1. the subject successfully responded to a visual stimulus: 898 . 2. no muscle or movement artifact was detected in the record; 2399 3. the event appeared within a string of at least eight consecutive hits; and 2500 4 the event did not occur at a hit/ miss transition or boundarv 1701 The power spectrum was calculated fcr the EEG data corresooncir ϊ to each 190 selected hit At each discrete frequency in the powe r spectrum the mean of the 1903 ppwer spectrum value from the 25 events was determined and this mean was 1904 used to produce the graph Thus tne graph of Fig 10 represents the mean of the 105 power spectrum calculated from the 25 hits r-cc To p roduce the grapπ of Fig 11 the EEG da to corresponding to 20 misses 2907 was selected according to the following criteria i03 1 the subject failed to respond to a visual stimulus and the failure was not 1909 attributable to external causes i9iϋ 2 no muscle or movement artifact was detected in the record 112 3 the event appeared within a string of at least eight consecutive misses, i9i2 4 the event did not occur at a hit/miss transition or boundary and 9 5 the event was not one of the first few misses among a string of misses 292 The mean power spectrum was calculated for the misses in a manner similar to 2925 that described for the hits i9i6 Activity at 60 Hz and odd harmonics thereof has been eliminated by noise 927 removal means 228 (described further in detail) from the power spectrum data 913 from which the graphs 164, 170 of Figs.10 and 11 were produced 229 The EEG data shown in Fig.10 (i e , the hit events) corresponds on 190 average, to a subject in an alert state commensurate with good performance. The i92i EEG data shown in Fig. 11 (i.e., the miss events) corresponds, on average, to a 1922 subject in a state of extreme drowsiness commensurate with failed performance. 2923 The scale of the vertical a/is cf graph 164 of F'g 10 is 10 times that of graoh 170 192 of Fist 11 Comparing the ςrapns of igs 10-11 it can ce ccserved that tne tota 1525 eπergv represented in tne cow?' spectrum 164 fπg 10 i cc-esocnding to hits is 2? s approximately 10 times greate r than the total e n ergy represented in the power 1927 spectrum 170 (Fig 11) corresponding to misse. Thus in this experimental data 1 3 there is a significant ccrrela ic, between a subject s state cf alertness drowsiness ? :-. o'sleep and the energy in tne EEG signai in certam frequeic, oaπαs The graphs i 9 jo of Fies 10-11 illustrate this pneπomeπon for the fre uencv ran e of 100-475 Hz i 9 τ. Easec on these and other experimental results it is believed that this correlation 1 32 is a σeneral characteristic of the human population that extends across a wide 1933 variety of different tasks It is believed that (T ) the recognition of a strong correlation between a 195 subject ' s sta e of alertness, drowsiness, or sleep, and the energy or amplitude 1 36 of the EEG signal at frequencies above 30 Hz, and (2) the exploitation of this 1937 correlation in a system for monitoring a subject ' s state of alertness, drowsiness 1933 or sleep, are novel and nonobvious aspects of the present invention As best seen 2939 in Figs 4 and 8, almost all of the energy in the subject s EEG signal is contained 1940 in components at frequencies below 30 Hz, regardless of the subject ' s state of 1 42 alertness, drowsiness, or sleep In Fig. 8, the only significant component visible 194 above approximately 30 Hz is a small peak at 60 Hz, which corresponds to noise 1943 derived from the commercial electrical power supply system As best seen in Fig. 1944 9, although the amplitude of the 60 Hz noise peak is minuscule compared to that 2945 of the predominant 12-13 Hz peak of Fig. 8, it is many times largerthan any 19J6 of the h i gher freαuencv s i gna l com p onents disovered to pe useful ,n drows i ness n-" detection Because or the l arge disoaπt / between t^e amp' tu d es c the 3-30 Hz 194) components an d any com p onents which ma / occ.r at -,gner rreαuences when 295 , di splayed at normal scales and resolutions the C-30 h: components effect i vely a,, mask the presence of h.gher-freαuenc/ components T us neitner the presence cf userul EEG s i gncl er^g, at r^nαes over 30 nz nor tne strong correlat i on 2 5 between sucn s i gnal energy and a su p ject s drowsiness „ould have been apparent us i ng normal data v i sual i zat i on and analvsis techniques Further accord i ng to 2955 trad i t i onal EEG d octr i ne rhythmic components ,n EEG signals at frequences n o aoove approx i matel y 30 Hz have been consiαereπ to be noise and ,n the pr i or art such components have been systematically discarded by fi l ter i ng or at best 1953 i gnored Thus the pr i or art discouraged practitioners from acquiring or explo i t i ng 2959 rhythm i c EEG s i gnal components above approximately 30 Hz in drows i ness mon i tor i ng and other related applications effectively teaching away from the 2962 present invention Components 226-244 are directed to extracting and exploiting i nformat i on 1963 wh i ch has been determ i ned to be useful m drowsiness detection from the acqu i red 1964 EEG s i gnal data that has been transformed using frequency analysis techn i ques The output of the power spectral density determinate means 226 i s 1900 suppl i ed to no i se removal means 228 via data path 326 The output of the power 1967 spectral dens i ty determ i nat i on means 226 is a vector m which each element i s a 1968 scalar represent i ng the power ,n the sampled EEG signal at each frequency po i nt 1569 produced in the frequencv analysis .)") The noise removal ans 223 replaces elements cr the power spectral 1 '1 density vector at cr adjacent to known noise frequences with πe'ghboπnσ values It is very difficult to avoid infiltration of noise at the commercial power line frequency (and its harmonics) in many real-world ap p lications due to the presence !9"4 of lighting and electrical eαuipment in the region in wh ch the data is collecte d i)-5 Even in clinical and research settings, it is nearly impossible to avoid acquiring 197 ^ such noise Accordingly fcr harmonics of the commercial power supply frequency :)7- (typically 60 Hz in North America) the noise removal means 223 replaces the 1975 power spectral density entries corresponding to that frequency and its upper and 179 lower nearest neighbors with the value cf one cf the nearest neighbors It may also isso be necessary to filter this power line noise in the analog domain depending on the 2931 amplitude of the noise which is present Other known sources of noise could be 1 982 similarly eliminated. 1 983 The output of the noise removal means 228 is provided to filter 2 9-* compensation means 230 via data path 328 The output 328 of the noise removal 935 means is essentially the same as the power spectral density (PSD) vector 326, 1986 with entries in the neighborhood of noise frequencies adjusted. 198 As discussed previously, the attenuation or cut-off portion of the frequency 1938 response of the anti-aliasing filter means 164 may intrude into the frequency range 1989 of interest in drowsiness detection, thereby attenuating signal components which i 9o may be useful The filter compensation means 230 compensates for the effect of 2S92 the anti-aliasing filter means 164 by multiplying the entries of the power spectral 29 9 2 density vector 326 by a suitable compensating function The comoensat i n σ 1) 2 funct i on depends on the characteristics of the particular filter used to implement 1 . filter means 164 For most f i lter designs of known characteristics such as the 2 95 s i xth order Butterworth filter 610 of Fig 6a, the compensating funct i on i s 199 relatively s i mple. The PSD vector entries corresponding to frequences in the 1 99 attenuation region are multiplied by a compensating factor computed from the 1992 filter character i st i cs The output cf the niter compensation means 230 i s prov i ded 2 to spectral bin grouping means 232 via data path 330 2 000 Accord i ng to an aspect of the present invention components 232-244 201 determine , respons i ve to the compensated PSD vector 330 , an output measure 2002 s i gnal on lead 138 ind i cating the subject ' s state of alertness , drowsiness , or sleep 2 003 In tr i als of an exper i mental embodiment of the invention using human sub j ects , a 20 04 high correlation has been observed between values of the output measure which 200 5 exceed a predefined threshold and the subjects' failures to respond to visual 2006 stimuli. 2007 Spectral bin grouping means 232 defines a plurality of selected frequency 008 ranges for which the energy in each range will be determined. For example , in a 2 009 preferred embodiment, five spectral bins may be defined as: 80-120 Hz; 120- 2 oιo 178 Hz; 181-240 Hz; 240-299 Hz: and 301-420 Hz. Although these spectral bins 2022 are used in preferred embodiments of the invention, other spectral bin choices 2oi2 could also be used. In fact, it is believed that at least some of the advantages 2 o i3 of the present invention may be obtained using one or more reasonably selected 20 subranges of frequencies between approximately 30 Hz and 500 Hz. Further, 2 0 1 5 the spectral b i n grouping means could also include a substantially wider ran * 2oi o of frequences than that represented in the spectral bins descπoed above For 2 o i7 example the spectra! b i n grouping means 232 could use some frequences below 023 30 Hz provided that a substantial contribution to the drowsiness output measure 2oi (discussed further i n detail) is based on signal information from frequences 22 above 20 Hz Although the value to a drowsiness detection system of EEG 2 C 21 signal components at requenc s above 300 H: has net been established , 2022 further investigation may reveal that components at frequencies above 500 Hz 2 0 2 3 may enable drowsiness detection systems to provide even better performance 2 0 2 The output 332 of the spectral bin grouping means 232 appears as groups 2075 correspond i ng to each bin containing the power-spectral density values for the 2 0 2 o frequences in that bin The output of the spectral bin grouping means 232 is 027 provided to bin energy determination means 234 via data path 332 ?0 » The b i n energy determination means 234 determines the energy in each 2029 spectral bin as the arithmetic sum of the power values in each bin The output 2 0 2 0 of the bin energy determination means 234 is a vector containing an entry for 2 0 i each bin representing the energy in the frequency range associated with that bin. 03 2 The output of the bin energy determination means 234 is provided to bin energy 2 0 3 inversion means 236 via data path 334. The bin energy inversion means 236 2 034 determines the inverse of the energy in each spectral bin. The use of the inverse in 2 03 5 determining the drowsiness output measure was chosen, in part , because it causes 06 a desired diminution of the output measure in response to increased energy in the 03 frequency bands of interest. Because increases in high-frequency energy tend to 038 indicate wakefulness the inverse operation produces an output measure havinσ a 2 9 low value for a wakeful alert supject and higher values lor drowsy suojects The 2 00 use of the i nverse also provides a non-linearity in the cutout measure which is 2 04 1 desirable because very low values of high-frequency energy correspond to extreme ?04 drows i ness and the non-linearity provides increased sensitivity at such low energy 2042 values However tne des i red output s i gnal response could be achieved using other 0, - funct i ons In some appl i cat i ons d i fferent types of output signals may be desired , 2 ,5 eliminat i ng the need for the inverse or the energy values in the frequency bands of 20 t interest could be used directly without additional processing 20 7 The output of the bin energy inversion means 236 is provided to weighting 2 0 S means 238 v i a data path 336 A weighting function generator means 240 suppl i es 20 49 a weighting function for use by the weighting means 23S in determining how to 2 050 weight the inverse energy value for each spectral bin in the final output measure 2 05 2 In a currently preferred embodiment of the invention each bin may be equally 2 05 weighted-e g the function provided by weighting function generator means 240 053 may be 1 , so that the weighting means simply multiplies each inverse energy value 205 by 1 Alternat i vely each b i n could be we i ghted by a different scalar value and 2055 these values could be generated by a function Improveα weighting functions may 056 be discovered through additional experimentation 57 The weighted inverted spectral bin energy values are provided by weighting 053 means 238 to summ i ng means 242 via a data path 338 The summing means 2 0 5 9 242 determines the arithmetic sum of the weighted inverse energy values and 3060 optionally presents it to interpolation means 244 via data path 340 2061 Interpolat i on means 244 may optionally be provided to generate an 2063 output s i gnal 138 (wh i ch may be an analog signal a discrete-valued s i gnal or 263 values i n memory) The i nterpolation means 244 thus continuously presents an 2064 output measure s i gnal 138 wh i ch represents the state cf the sub j ect ' s alertness , 3065 drows i ness , or sleep , dur i ng the period in which the sub j ect's EEG s i gnal was 20όo acqu i red. If the i nterpolat i on means 244 ,s not used, then the output of summ i ng 3067 means 242 may be used i n a similar manner to produce the drows i ness output 20o3 measure by connect i ng lead 340 to output signal lead 133 (as ind i cated bv 306 broken l i ne 138a) If the drowsiness analysis process is performed essent i ally i n 3070 " real-t i me ," then the output measure signal 138 will lag the sub j ect's actual 3071 drows i ness by an i ns i gn i f i cant amount Thus the output measure s i gnal 138 may 3073 be effect i vely used to produce an alarm when the output measure signai exceeds a 3073 threshold i nd i cat i ng that the subject has become excessively drowsy and may be 2074 unable to safely perform a task As best seen in Fig 5 , the output measure s i gnal 2075 138 may be cond i t i oned by comparison with a threshold and by indications that 2076 artifact may be present i n the subject ' s EEG signal, before it is used to produce a 2077 display or alarm. 2078 Figs. 20a-20c are graphs of the output measure produced in an 2079 exper i mental embod i ment of the i nvention using EEG data acquired during 2080 trials us i ng human subjects. The output measures shown are equivalent to and 2081 determ i ned i n substant i ally the same manner as the output measure suppl i ed on 2082 lead 133 from the first and second embodiments 700 and 800 (Figs. 12-16) of 20S3 the i nvent i on The dependent variable in each graph is the value of the output 2084 measure Each po i nt on the plots corresponds to the output value determ i ned 20S 1 from a 2048-poιnt data analvs i s window of the subject s EEG d ata -sα « cr i D « d 208o prev i ously i n connect i on w i th data selection means 220 adjacent w i n d ows of data 2087 are offset Dy 1024 samples so that each wmαow includes 1024 samp l es of new ' 20SS data and 1024 samples from the previous window The independent var i able i n 2 OS 9 each of F i gs 20a-20c t i me i n seconds over which the measure was determ i ned 209C The graphs of F i gs 20a-20c each include upper and lower thresnoid 2091 markers 916 and 918 which extend horizontally at 17 and 12 units on the 2093 output measure scale respectively Output measure values increase w i th 2093 sub j ect drows i ness Output measure values which exceed the upper t h reshold 201 916 i nd i cate that the subject was excessively drowsy (or asleep) dur i ng the 2095 correspond i ng t i me i nterval Output measure values which are below tne 209o lower threshold 918 i nd i cate that the subject was wakeful or alert dur i ng the 2097 correspond i ng t i me i nterval Output measure values which lie between the 2098 thresholds 916 and 918 do not definitely indicate either alertness or d rows i ness 2099 However , the slope of the output measure can be used to interpret behav i oral 2100 trends For example , i f the output measure is between the thresholds and the 2101 slope i s pos i t i ve the subject is approaching an excessively drowsy cond i t i on If the 2102 slope is negative , the subject is approaching a wakeful condition 2103 The thresholds 916 , 918 were selected based on results from ex p erimental 2104 exam i nat i ons of several subjects, and are believed to be generally app li cable 2105 for drows i ness detect i on As best seen ,n Fig 19 (described further m deta i l) , 2106 a su i table method for determining an appropriate drowsiness thresholc wh i ch 3io" is tailored to a particular individual using EEG signal data acquired from that individual could also ce used in conjunction with the drowsiness detection svstem 2ic9 of the present invention 2110 Fig 20a is a graph 910 showing two superimposed output measure plots 912 and 914 which were determined using EEG signals acquired from a human i subject during two different examination periods Each of the segments shown represents approximately P 62 minutes A first plot 912 was derived from data 722 4 acquired while the subject was observed to be in a state of wakefulness. based on 3ia the subject s physical behavior while performing a visual task A second plot 914 3iio was derived from data acquired later, while the subject was resting on a bed and 3ii7 allowed to fall asleep in the first clot 912 the output measure never exceeds the 3U3 lower threshold 918 In the second plot 914 the output measure exhibits several 2229 initial episodes in which the measure crosses the upper threshold followed by a 22 0 period during which the measure continuously exceeds the upper threshold 916 3i i and tends to increase over time Plot 912 is believed to be representative of the i33 behavior of the output measure for wakeful human subjects Plot 914 is believed 2123 to be representative of the behavior of the output measure for a human subject 224 who becomes drowsy and falls asleep The output measure is believed to be a 2125 sensitive indicator of the drowsiness of a human subject 2i 6 Fig. 20b is a graph 920 showing an output measure plot 922 determined 212 7 using EEG signals acquired from a human subject during approximately 9.29 212 8 minutes of an examination period The EEG signals were acquired during subject 21 9 performance trials conducted as previously described in connection with an 2130 exper i mental embod i ment of the invention The suoject was instructed to respond to a ser i es of visual st i muli presented at random intervals If the suPject fa i led to respond to a st i mulus w i thin a predetermined period of t i me , additional 21Z3 stimul i were presented at regular time intervals until the subject responded to 2134 two consecut i ve st i mul i Each stimulus is represented en the graph 920 as a solid 2i35 vert i cal l i ne (such as l i ne 924) ,f the subject responded successfully or as a broken vert i cal l i ne (such as l i ne 926) if the suoject failed to respond Graph 920 reveals excellent agreement between the subject s actual drowsiness , as i nd i cated by test 2i3 s performance and that pred i cted by the output measure 922 In general adequate 3i39 subject performance wh i ch is indicated by solid vertical lines , ,s accompan i ed by 7240 values of the output measure below the lower threshold 918 Response fa i lures 928 and 930 i nclud i ng a group 928 of 15 contiguous response failures attr i buted to a sleep ep i sode last i ng approximately 90 seconds are predicted by excurs i ons of 2i43 the output measure above the upper threshold 916 prior to the actual fa i lure 3,44 Fι ε 20c ιs a S ra P h 32 show i ng an output measure plot 934 determined 2245 us i ng EEG s i gnals acqu i red from a human subject during approx i mately 883 214o m i nutes of an exam i nat i on period EEG signals were acquired dur i ng subject 2147 performance tr i als conducted as previously described in connect i on w i th an 2243 experimental embod i ment of the invention Graph 932 also shows excellent 249 agreement between the subject ' s actual drowsiness as indicated by test performance , and that predicted by the value of the output measure 934 Graph 2252 932 suggests that even i f the output measure does not exceed the upper threshold 3i53 916, performance fa i lures due to αrowsmess may be predicted pr i or to a failure 3i53 when the va l ue of the out p ut measure is both above the lower threshold and i54 i ncreasing (i e tne output measure has positive slope ) , 5 Graph 932 a i so i ncludes two p erformance failure episodes wh i ch were not the „ result of excess i ve drows i ness i nstead these failures occurred because the sub j ect was not pay i ng attention to the visual stimulus display 725 F i g 15 i s a flow d i agram showing the steps of a method 540 for 2255 determ i n i ng wnether a segment or window of sampled EEG d ata l i kely conta i ns 2160 art i fact Accord i ng to an as p ect of the present invention the method 540 may ?lol be used ,n con j unct i on w i th artifact detection means 148 of d rows i ness detect i on 22o2 system 500 , i nclud i ng the f i rst or second preferred embodiments 700 , 800 The 2263 art i fact detect i on method 540 ma y be implemented using data acqu i s i t i on and process i ng har d ware and software in common with that used to i mplement the 2165 f i rst or secon d p referred embod i ments 700 800 Alternatively separate data ?16o acqu i s i t i on anc process i ng hardware and software could also be used The art i fact 2167 detect i on met h od may also be used ,n other applications ,n which ,t ,s des i red to 31 S analyze only artifact-free EEG data 3169 In step 5 <1 2. appropr i ate EEG data ,s acquired from a sub j ect , and a w i ndow of data of appropr i ate length ,s selected therefrom If the artifact detect i on method ,s operat i ng i n conjunction with one of the first or second preferred embod i ments 700 800 of the invention discussed previously the selected data 21 3 may , for example be obta i ned from tne output of the data select i on means 220 2174 of Figs ϋ and 15 Thus ste p 542 may represent all of the data acqu.s i t.on and related funct i ons and steps Performed up to and including data selecfon means 2276 220 of Figs 14 and 15 However i f i t i s desired to use the artifact detection 7177 method 540 i n conjunct i on with an embodiment which does not itself requ i re 21 7 3 these data acqu i s i t i on and selection functions then means for acquir i ng and 2279 selecting the data may be separately provided - ,J5 ° ln ste P 558 lf ar N samples contained in the selected data window appear 2i ει to be saturated the data w i ndow is rejected A sample may be considered to be 37 saturated if its value is the largest or smallest value which can be produced by 3i s 3 the analog-to-d i g i tal conversion means 212, or if the analog-to-digital convers i on 21S means 212 provides an overflow or similar out-of-range indication is5 In step 544 a suitable transform or process is applied to the selected 213 c data to convert i t from a time-domain representation to a frequency-doma i n 13 7 representation A πy suitable frequency analysis method may be used Several 2i3s usable frequency analysis methods are discussed in connection with the frequency 2139 analys i s means 224 of F i gs 14-15 and step 544 may be implemented using the 2i9 o same frequency analysis technique as used to implement frequency analysis means 2191 224 2i9 In step 546 the power spectral density of the selected data segment is 2i93 determined using the results of the frequency analysis step The power spectral 22 density may be determined as described previously in connection with power 2i95 spectral density determination means 226 of Fig 14-15 i96 In step 548 one or more spectral bins which correspond to selected i 7 frequency ranges , are def i ned For each spectral bin the energy contained in that 2198 b i n is determ i ned using the power spectral density of all frequency components 2199 within that b i n In a preferred embodiment of the invention a single spectral b i n 32 oo correspond i ng to the freαuency range 30-60 H: was definec This frequency 2 0 2 range appears to have adequate sensitivity for use in an artifact detection 2 30 algor i thm However a d i fferent number of spectral bins, covering other frequency 2203 ranges could also be used If multiple spectral bins are defined the energy in all 220 spectral bins may be combined to produce an aggregate energy value 2205 In step 550 the mean value of the energy αetermmed in step 548 over all 22 0 windows in the current segment of EEG data is determined This mean value 3307 is then multipl i ed by an artifact threshold factor to obtain an artifact threshold 208 value The art i fact threshold factor may be experimentally determined In an 3 2 09 experimental embod i ment of the invention an artifact threshold factor of 20 was 2 io used 2 i ι In step 552 , the aggregate energy values which were determined in step 2222 548, are compared to the artifact threshold value For all of the windows i n a 2i 3 particular EEG data segment if the aggregate energy in any window exceeds the 33i 4 artifact threshold value determined in step 550 that window is considered to 222 5 be contaminated by artifact and is indicated accordingly in the artifact output 2226 measure on lead 150 2i7 In step 556 , wh i ch is optional, inputs 164 representing certain secondary 2218 artifact detection factors or indicators may be considered in determining whether i9 EEG data i s contam i nated by artifact For example, signals or other informat i on 2 0 may be derived from a video camera or a motion detector whicn could be useful 2221 for artifact detection Thus , the artifact detection method 540 may produce an 3223 artifact present signal 150 if artifact is detected in either the comparison step 552 222z or the secondary artifact consideration step 556 222 In some applications it may De unnecessary to explicitly provide artifact :25 detection means 128 including means for implementing the artifact detection 222o method 540 of Fig 18 Implementations which produce an output measure as 2227 described in connection with the first and second preferred embodiments 700 7225 800 are tolerant of artifact, in that tne presence of artifact is accompanied by a 322° significant increase in EEG signal energy across all frequencies which causes a 230 decrease in the value of the output measure Artifact resulting from sustained 3231 movement or other activity indicating wakefulness changes the output measure 3333 in the direction of wakefulness and does not detract frcm the performance of the 22J3 drowsiness detection system 223 Fig 19 is a flow diagram showing the steps of a method 510 for 2235 determining a drowsiness threshold for an individual based on EEG signals 226 acquired from the individual while alert According to an aspect of the present 727 invention, the method 510 may be used in conjunction with the threshold value 2238 generator means 144 of drowsiness detection system 500 including the first or 2239 second preferred embodiments 700 800 The individual threshold determination 240 method 510 may be implemented using data acquisition and processing hardware 24i and software in common with that used to implement the first or second preferred 2 embodiments 700, 800 Alternatively, separate data acquisition and processing 2 3 hardware and software could also be used 2244 In step 512 EEG data is acquired for use in the method from an alert or 2245 wakeful subject To obta i n sufficient EEG signal information from an i nd i v i dual to 246 determ i ne an appropr i ate threshold data >s preferably acquired over a thresnold 22 7 determ i nat i on per i od of a p proximately 10 minutes although adequate results 223 may be obta i ned even i f z substantial deviation from this period i s used The EEG s i gnal i nformat i on needed for use with the threshold determinat i on method 3250 510 may be acqu i red anc stored for subsecueπt processing , or may be used i n 2 1 real t i me If the threshol d determination method 510 is operating in conjunct i on w i th one of the f i rst or second preferred emoodiments 700, 800 of the invent i on 3353 discussed prev i ously the αata ma y for example be obtained from the output of 7254 storage means 214 of F i gs 14 and 15 Thus, step 512 may represent all of the 2255 data acqu i s i t i on and related functions and steps performed up to and i nclud i ng 7250 storage means 214 of Figs 14 and 15 2 7 In step 520 , the drowsiness output measure is determined using the data 228 collected during the thresnold determination period The output measure may 2259 be determ i ned as prev i ously disclosed in connection with the drowsiness analys i s 2250 means 856 , 126 In some appl i cat i ons one of the drowsiness analysis means 136 , 2261 756 , 856 of F i gs 14 or 15 may be operated during the threshold determ i nat i on 2262 per i od , and the output measure produced tnerefrom may be recorded for 26 subsequent use In step 522. an artifact indication signal is preferably determined 64 over the same threshold determination period For example , the artifact detect i on 25 means 148 may be operated in conjunction with the artifact determinat i on 726o method 540 of F i g 18 dur i ng the thresnold determination period , and the artifact 2 67 i nd i cat i on signal producec therefrom may be recorded for subsequent use 2263 In step 514 the results of artifact detection step 522 are used to condition 69 the results of the output measure determination step 520, so that only output 22 0 measure information produced during artifact-free periods is retained In step 3 7i 516 the mean of the output measure over the artifact-free portion of the 2272 threshold determination period is determined The mean output measure value 2 7 3 determined in step 516 establishes an alert baseline value of the output measure 2274 for the individual subject being examined Because artifact in EEG data tends to 2275 reduce the output measure removal of output measure values determined during 7 76 artifact-contaminated periods avoids an artificial reduction of the Paseime value. 7277 thereby reducing the potential for false positive drowsiness indications 22 78 In step 518 the subject ' s alert baseline output measure value determined in 22 79 step 516 is multipl i ed by an appropriate threshold factor to determ i ne a threshold 2230 value for that indiv i dual The threshold value may be supplied on lead 146 for 22 81 use by comparison means 140 (Fig 5) When the output measure produced 72 3 2 by the drowsiness analysis means 136, 756, 856 exceeds the threshold value 83 determined in step 518 , the individual may be considered to be excessively drowsy 22 3 For an experimental embodiment of the invention similar to the first preferred 2 3 5 embodiment 700 disclosed herein, a threshold factor of 3.0 was found appropriate 2286 for advance detection of extreme drowsiness likely to affect performance 237 Other threshold factors may be appropriate for different applicat i ons For 2233 example , a smaller threshold factor could be used to increase sensit i vity to detect 2239 even mild drowsiness Larger threshold factors may be appropriate for use i n detecting particular sleep stages or for monitoring aspects of a subject s state of 2291 consciousness or response to anesthesia 393 F, δ 17 ιs a sιm P'ιπed block diagram showing the general structure of a 293 th i rd embo di ment 400 of an EEG signal processor means 124 for use i n an EEG- 229 based drows i ness mon i toring system 100 which is constructed accord i ng to the 2295 present i nvent i on an d wh i ch is adapted for use with analog signal process i ng and 229o related techn i ques Analog signal processing techniques can often be eff i c i ently i mplemented at very low cost As best seen in F,g 17 the third embodiment 400 2298 rece i ves ampl i f i ed analog EEG signals as input data on lead 130 analyzes the 2299 s i gnals , and respons i vely produces an output measure on lead 138 ind i cat i n σ the 23oo subject ' s state of alertness drowsiness or sleep 2 oι The t hi rd embod i ment 400 provides functions similar and largely 307 analogous, to those of the first and second embodiments 700 , 800 , but 2303 employs a d i fferent conf i guration of functional components suitable for analog 2304 s i gnal processing Like the first and second embodiments 700 , 800 , the th i rd 2305 embod i ment compr i ses means for selecting particular EEG signal components 730o wh i ch have heretofore p een discarded or ignored in accord with the teachings 37 of the pr i or art , and means for exploiting the information contained in these 308 components to prov i de a reliable measure of a sub j ect s alertness or drows i ness 2309 In general , the th i rd embodiment 400 defines several frequency ranges of interest 21 (denoted a-m) and i ndependently processes the selected frequency ranges i n 3U respect i ve analog s i gnal processing channels 408a-408m The reference numbers 23i2 referr i ng to the components of each channel are denoted with a suff i x letter 2313 i nd i cat i ng the channel to which they belong like reference numbers (exclud i ng the 324 suffix letter) generally denote eouivalent components in each channel A referen ce 2 3i5 lacking a suffix letter made to a component common to all channels i s i ntended 2 3 i o to refer to the respect i ve equivalent components in all channels Any suitable 3 3 1 7 number of frequency ranges and associated separate processing channels 408 , 2318 may be used (such as the five frequency ranges defined in the first and second .- 2 embodiments 700 800) For this discussion and the diagram of Fι* 17, the 7220 choice of subscripts i s not meant to imply the presence of a particular number of 232i frequency ranges or processing channels In each channel 408a-408m amplified analog EEG data received on 2 J 2 3 lead 130 is presented to frequency selection means 410 For each channel , the 23 4 frequency selection means 410 allows signals within tne channel s predef i ned 2325 frequency range for that channel to pass while substantially attenuat i ng s i gnals 2 3 2 6 outside of the channel s frequency range Frequency selection means 410 may be 2327 i mplemented using any appropriate filter design, of which many suitable ones are 32 8 well known in the signal processing art The selection of the frequency range of 7J 2 each channel is preferably done cooperatively with the design of the correspond i ng 330 frequency selection means Depending on cost availability, and other practical 233i aspects of filter design , the frequency ranges to which respective channels are 7332 directed may overlap or be mutually exclusive 22Z3 The output of frequency selection means 410 which appears as a band-pass 2 334 filtered version of the original analog EEG input signal is provided via lead 412 to 2335 energy determ i nation means 414 For each channel, energy determination means . "3o 414 integrates the i nstantaneous power of the bandpass-filtered signal over a time 2J37 interval I preced i ng the current time r J/_ τ χ2, (r) d - Energy determ i nat i on 233S m< eans 414 may be i mplemented using anv suitable integrator , designs for which 2339 are well known in the art 2340 The output of energy determination means 414 is supplied v i a lead 416 to 2341 energy value i nversion means 418 Inversion means 418 determines the i nverse of the energy i n the band-pass filtered signal (for each channel) As d i scussed above 3343 in connect i on w i th bin energy inversion means 236 of Figs. 14-15 , the use of an 2344 invers i on i n produc i ng the drowsiness output measure was selected because i t 345 produces an output measure having a particular desired response character i st i c. 24 However , the des i red output signal response could be achieved us i ng other 347 funct i ons. In some applications, different types of output signals may be des i red. 7343 el i minat i ng the need for the inversion, or the energy values in the frequency bands 7 9 of interest could be used directly without additional processing. 2350 The output of the energy value inversion means 418 is provided v i a lead 420 to we i ght i ng funct i on means 422. For each channel, the weighting funct i on 2352 means 422 appl i es a respective predetermined weighting function to the i nverted 2353 energy value. In a preferred embodiment of the invention , the energy from each 2354 channel may be equally weighted using the weighting function f ch ( u ) = ku. In th i s 2355 case, the weight i ng function is multiplication by the same constant value k In a 2356 preferred embod i ment , k may be set to 0.2. However, each bin could be we i ghted 7357 by a different scalar value, and more generalized weighting functions could also be 2358 used The we i ght i ng function means may be implemented using a commerc i ally 259 available operational amplifier or a multiplier. 2360 The outputs of the we i ghting function means 422a-422m for all channleels 23c2 are preferably prov i ded via leads 424a-424m to a single summation means 426 2362 wh i ch determ i nes the sum of the weighted inverse energy values over all channels 2363 The summat i on means 426 may be implemented using a commercially ava i lable 2364 analog summer or adder or a network of operational amplifiers The result i ng 23o5 sum is provided as the drowsiness output measure signai on lead 138 2Joo The above-oescπbed embod i ments of the invention are merely examples of 2367 ways i n wh i ch the i nvent i on may be carried out Other ways may also be poss i ble , 23oS and are w i th i n the scope of the following claims defining the invention 2369 BIBLIOGRAPHY 2370 Akerstedt T and M G i liberg (1990) " Subjective and Objective Sleepiness in th< 23 i Active Individual " Int J of Neuroscence, Vol 52 , pp 29-37 732 Akerstedt , T , G Kecklund and A Knutsson (1991). " Manifest Sleepiness and 233 the Spectral Content of the EEG during Shift Work ' Sleep Vol 14 No 3 2374 pp. 221-225 3375 Carskadon, M A A Rechtschaffen (1987) Chapter 73 Monitoring and Staging 3376 Human Sleep -from Principles and Practice of Sleep Medicine - section 3377 7 methodology, by W B Saunders, Philadelphia PA 2378 Chiu U S Patent No 4875030 3379 Daniel R S (1967) " Alpha and Theta EEG in Vigilance ' D ercept Mot Skills 7330 Vol 25, pp 697-703 232 Dmges, D F (1988), " The Nature of Sleepiness Causes Contexts and 332 Consequences ," Perspect i ves in Behavioral Medicine A Baum and A 23S Stunkard, eds , New Jersey, Erlbaum 233 Dingus, T A , H L Hardee, and W W Wierwille (1987) ' Development of Models 235 for On-Board Detection of Driver Impairment, " Accident Analysis and 2336 Prevention, Vol 19. No 4 2337 Duff , G W , unpubl i shed observat i ons (referenced in Mitler et al 1988) 388 Estrada U.S. Patent No 3.953.831 2389 Fruhstorfer. H P Laπganke , K Meinzer, J H Peter, and U Pfaff (1977) 2390 " Neurooπys i olog i cal V i gilance Indicators and Operational Analys i s of a 2391 Tra i n V i g i lance Mon i toring Device A Laboratorvand Field Study , " i n 2392 V i g i lance ed i ted by R R Mackie Plenum Press New York , pp 147-162 Ga i llard. J M (1987) , ■ EEG Fitting A New Method for Numerical Analys i s of 2 EEG, " Neuropsychobiology, Vol. 17 , pp 9-14 2395 Gale A (1977) Some EEG Correlates of Sustained Attention " ,n V i g i lance , 2396 ed i ted by R R Mack i e , Plenum Press, New York pp 263-283 2397 Horvath M E Frant i k K Kopπva and J Meissner (1976) "EEG Theta 2398 Act i v i ty Increase Co i nc i ding with Performance Decrement in a Monotonou; 2399 Task ," Act i v Nerv Sup (Praha) Vol 18, pp 207-210 700 Kishi U S Patent No 5.311,877 402 Langlo.s. P H . M H Smolensky. B P Hsi, and F W We.r (1985) , Temporal 240 Patterns of Reported Single- Vehicle Car and Truck Accidents in Texas , 2403 U.S. A .dur i ng 1980-1983, " Chronobiology Int , Vol 2. pp. 131-140 3404 Lav i e , P , M Wollman , and I Pollack (1986), " Frequency of Sleep Related Traff i c 2405 Acc i dents and Hour of the Day ," Sleep Research. Vol 15, p 275 2406 Make i g , S. and M Inlow (1993), " Lapses in Alertness Coherence of Fluctuations 2407 i n Performance and EEG Spectrum, " Electroencephalography and Clin i cal 203 Neurophysioiogy, 86(1993)23-35 2409 Mitler , M M , M A Carskadon, C A Czeisler, W C Dement. D F D i nges , and 2410 R.C. Graeber (1988) , " Catastrophes, Sleep, and Public Policy Consensus ™* Report " Sleep, Vol 11, No 1, pp. 100-109 2422 Molodofsky , H. (1992). ' Evaluation of Daytime Sleepmess." Clinics in Chest < i2 Medicine, Vol. 13 No 3, pp 417-425. 2414 O ' Hanlon, J F and j Beatty (1977) , " Concurrence of Electroencephalograph i c 2425 and Performance Changes During a Simulated Radar Watch and Some 2416 Impl i cat i ons for the Arousal Theory of Vigilance," ln Vigilance , ed i ted by 2417 R.R Mack i e , Plenum Press. New York pp. 189-202. 2418 O ' Hanlon. J F and G R Kelley (1977), " Comparison of Performance and 2419 Phys i olog i cal Changes Between Drivers Who Perform Well and Poorly 2420 Dur i ng Prolonged Veh i cular Operation ,n Vigilance , ed i ted by R.R - i2i Mackie, Plenum Press, New York, pp. 87-110. 42 Olsen etal U S Patent No 5,311,876 2422 Planque , S. , D Chaput , C Pet i t , and C Tarπere (1991) , "Analys i s of EOG and EEG Signals to Detect Lapses of Alertness in Car Simulation Driv i ng , " 2425 Presented at the 13th ESV Conference, Paris France, November 4-7.1991 2426 (referenced in Wierwille (1994)) 2427 Pritchard, W.S. (1995) , " Measuring ' Chaos ' m the Brain. A Tutor i al Rev i ew of 2428 EEG Dimens i on Est i mat i on, " Brain and CogniMon , Vol. 27 , pp. 353-397. 2429 Santamaria, J. and K.H. Ch i appa (1987), " The EEG of Drows i ness i n Normal 2430 Adults ," J. of Clinical Neurophysioiogy, Vol 4 , pp.327-382 23i Sekoetal. U.S. Patent No 4.564,833. 232 Sekoet al. U.S. Patent No 4,604,611. 43 3 Slansky U.S. Patent No 4,617,559. 2434 Smith. J.R. (1987) , " Computer Analysis of Sle«p Data," in Clinical Applications 23 of Computer Analys i s of EEG and other Neurophysiological Signals , 2436 Handbook of EEG and Clinical Neurophysioiogy (revised series. Vol. 2<37 2 ). F - H - L °P es da Silva . W . Storm van Leeuwen, and A. Remond , Eds.. 438 Amsterdam: Elsevier, pp. 131-145. 2439 Stanford Sleep Disorders Clinic and Research Center (1991) , "Why Should We 2440 Care About Sleep ? The Toll of Daytime Sleepiness and Sleep D i sorders on 24 Society. " Torsvall , L . and T Akerstedt (1987), " Sleepiness on the Job: Continuously 7443 Measured EEG Changes in Train Drivers," Electroen. and Clinical 44 Neurophy., Vol.66, pp.502-511 . 2445 Torsvall , L , T. Akerstedt, K. Gillander, and A Knutsson (1989) , "Sleep on 3446 the Night Shift: 24-Hour EEG Monitoring of Spontaneous Sleep/Wake 2 Behavior, " Psychophysiology, Vol.26, No.3 , pp.353-358 . 448 Trejo , L.J. , and MJ. Shensa (1993), " Linear and Neural Network Models for 24 Predicting Human Signal Detection Performance from Event-Related 450 Potentials: A Comparison of the Wavelet Transform with Other Feature 451 Extraction Methods ," Navy Personnel Research and Development Center 2452 (San Diego , CA) , and Naval Command Control and Ocean Surveillance 2i5 Center. RDT&E Division (San Diego , CA) . 245 U.S. Department of Health and Human Services (1992) , Wake Up America: A 45 National Sleep Alert, Vol. 1: Executive Summary and Executive Report , 5 Report of the National Commission on Sleep Disorders Research. 2457 Wierwille , W.W. , S.S. Wreg i tt, and M.W. Mitchell (1992), "Research on Vehicle- 245 Based Dr i ver Status/Performance Monitoring, " , First Semiannual Research 25 Report, DTNH22-91-Y-07266. 24ό o Wierwille, W. W. , S. S. Wreggit, C. L. Kirn, L. A. Ellsworth, R. J. Fairbanks 461 (1994) , " Research on Vehicle-Based Driver Status/Performance 3463 Monitor i ng: Development, Validation, and Refinement of Algorithms for 63 Detection of Dr i ver Drowsiness, " Final Report, Report DOT HS 808247 2464 Wilkinson , R T. and D. Houghton (1975), " Portable Four-Choice React i on Time 465 Test w i th Magnet i c Memory, " Behav Res. Meth. Instru. Comp., Vol.7, p. 2466 441. 2 4 6 7 Yoshimiet al U.S. Patent No. 4,928,090.

Claims: WHAT WE CLAIM IS: 1. A system adapted to measure in a human subject a sterτe"a " KTftg a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia comprising: means for acquiring a brain-wave signal having rhythmic components from the subject; means for selecting components of said brain-wave signal lying in a predetermined range of frequencies including frequencies above 30 Hz: means responsive to said selecting means for determining a contribution to said brainwave signal due to said components lying in said range; and means responsive to said contribution-determining means for producing an output- measure signal indicative of said state of the subject. 2. A system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia comprising: means for acquiring a brain-wave signal having rhythmic components from the subject; means for defining a range of frequencies including above 30 Hz; means for analyzing frequency content of said acquired brain-wave signal responsive to said defined range of frequencies for defining a plurality of frequencies or frequency bands within said range, and for determining for each of said defined frequencies or frequency bands a contribution due to a component of said acquired brain-wave signal at such frequency or within such frequency band; and output-measure-determining means responsive to said frequency analysis means for producing an output signal indicative of said state of the subject. 3. A system adapted to measure in a human subject a state along a continuum relating to aleπness, drowsiness, sleep, unconsciousness, or anesthesia comprising: means for acquiring a brain-wave signal having rhythmic components from the subject; means for defining a range of frequencies including frequencies above 30 Hz; means for analyzing frequency content of said acquired brain-wave signal responsive to said defined range of frequencies for defining a plurality of frequencies or frequency bands within said range, and for determining for each of said defined frequencies or frequency bands an amplitude or distribution of amplitudes corresponding to a component or components of said acquired brain-wave signal at such frequency or within such frequency band; power determination means for determining for each of said defined frequencies the power contained in said corresponding component or for each of said defined frequency bands the distribution of power contained in said corresponding components; and spectral-grouping means for defining a limited number of spectral groups, each of said spectral groups defining a subrange of said range and including at least one of said defined frequencies or frequency bands; group-energy-determining means responsive to said power determination means and said spectral-grouping means for determining for each of said spectral groups an aggregate energy represented in a plurality of said components corresponding to the defined frequencies or frequency bands included in such group: and output-measure-determining means responsive to said group-energy-determining means for producing an output signal indicative of said state of the subject. 4. A system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia comprising: means for acquiring a brain-wave signal having rhythmic components from the subject; anti-aliasing filter means for substantially removing from said brain-wave signal components outside a first predefined frequency range; means responsive to said brain-wave signal-acquiring means for converting said brainwave signal into a sequential plurality of digital data elements representing the amplitude of said brain-wave signal at a corresponding plurality of discrete times: means for selecting for analysis a limited number of said data elements; means for defining a second range of frequencies including frequencies above 30 Hz: means for analyzing frequency content of said acquired brain-wave signal, said means being operative to define a plurality of frequencies or frequency bands within said second range, and said means also being responsive thereto and to said selected data elements to determine for each such defined frequency or frequency bands an amplitude or distribution of amplitudes of a component or components of said acquired brain-wave signal at that frequency or within that frequency band: power determination means for determining for each of said defined frequencies the power contained in said corresponding component or for each of said defined frequency bands the distribution of power contained in said corresponding components; spectral-grouping means for defining a limited number of spectral groups, each of said spectral groups defining a subrange of said second range and including at least one of said defined frequencies or frequency bands; group-energy-determining means responsive to said power determination means for determining for each of said spectral groups an aggregate energy represented in a plurality of said components corresponding to the defined frequencies or frequency bands included in such group; and output-measure-determining means responsive to said group-energy-determining means for producing an output signal indicative of said state of the subject. 5 The system of claim 4 wherein said power-determining means further comprises noise-removal means for identify ing frequencies or frequency bands at which noise is expected, and for substituting tor the power value or power distribution at each such identified frequency or frequency band a power value or power distribution corresponding to a neighboring frequency or frequency band 6 The system of claim 4 further comprising means for causing said output- measure determination to occur at repeated predefined intervals, and wherein said output- measure-determining means further comprises interpolation means for continuously presenting an output signal representing a current value of said output measure until a next-determined value of said output measure becomes available 7 A system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia comprising means for obtaining a bram-v ave signal from the subject. subject state analysis means responsive to said brain-wave signal obtained from the subject for (a) selecting components of said brain-wave signal lying in a first predetermined range of frequencies including lrequencies above 30 Hz, (b) determining a contribution to said brain-wave signal due to said components lying in said range and (c) producing, responsive to at least said determined contribution, an output- measure signal indicative of said state of the subject 8 The system ot claim 7 further comprising threshold means responsive to said output-measure signal for producing a non-continuous-valued output signal indicating that the subject is drowsy 9 The system of claim 8 wherein said threshold means comprises means for determining a threshold value relating to said output-measure signal whereby said threshold value defines a range of output-measure-signal values corresponding to a defined physiological state of the subject, and means for comparing said output-measure signal with said threshold value to produce a non-continuous-valued output signal indicating that the subject is in said defined physiological state 10 The system of claim 9 wherein said threshold-value-determining means employs a preselected threshold value believed to be generally applicable to human subjects. 1 1 The system of claim 10 wherein said threshold-value-determining means employs a threshold determined for the subject responsive to a brain-wave signal obtained from the subject while in a known state along said continuum. 12. The system of claim 8 wherein said threshold means comprises: means for determining a threshold value relating to said output-measure signai whereby said threshold value defines a range of output-measure-signal values indicating that the subject is drowsy; and means for comparing said output-measure signal with said threshold value to produce a non-continuous-valued output signal indicating that the subject is drowsy. 13. The system of claim 7 further comprising threshold means responsive to said output-measure signal for producing a non-continuous-valued output signal indicating said state of the subject. 14. The system of claim 7 further comprising artifact detection means for determining whether said brain-wave signal obtained from the subject is apparently contaminated by artifact. 15. The system of claim 14 further comprising means responsive to said artifact detection means for eliminating from said output signal an effect of said brain-wave signal obtained when artifact is determined to be present. ! 6. The system of claim 14 wherein said artifact detection means is responsive to said brain-wave signal acquired from the subject to determine whether said brain-wave signal obtained from the subject is apparently contaminated by artifact. 17. The system of claim 14 wherein said artifact detection means is responsive to information other than said brain-wave signal obtained from the subject to determine whether said brain-wave signal obtained from the subject is apparently contaminated by artifact. 18. The system of claim 7 further comprising means responsive to an output signal of said system to produce a warning indicating that the subject is not fully alert. 19. The system of claim 18 wherein said warning is an audible alarm. 20. The system of claim 18 wherein said warning is a visual display. 21. The system of claim 7 wherein said subject state analysis means comprises: at least one brain-wave signal-processing channel; each such brain-wave signal-processing channel being directed to a corresponding range of frequencies and supplying an output representing a contribution to said brain-wave signal obtained from said subject due to components thereof lying substantially within said corresponding range of frequencies; and means responsive to said outputs supplied by a plurality of such brain-wave signal- processing channels to produce said output-measure signal indicative of said state of the subject. power determination means for determining for each ofp j ibJKid frequencies a power represented in said corresponding components or tor each ot said defined frequency bands a power distribution represented in said corresponding components, and spectral-grouping means for defining a limited number of spectral groups, each of said spectral groups defining of frequencies within said first predefined range of frequencies and including at least one of said defined frequencies or frequency bands, group-energy-determining means responsive to said power determination means and said spectral-grouping means for determining for each of said spectral groups an aggregate energy represented in a plurality of said components corresponding to the defined frequencies or frequency bands included in such spectral group, and output-measure-determining means responsive to said group-energy -determining means for producing said output-measure signal indicative of said state of the subject 28 The system of claim 27 wherein said brain-wave signal acquisition means includes a first digital control and processing means and said subject state analysis means includes a second digital control and processing means 29 The system of claim 27 wherein said system includes a digital control and processing means operatively connected to said brain-wave signal acquisition means and to said subject state analysis means, whereby said digital control and processing means performs at least some processing functions tor each of said brain-wave signal acquisition means and said subject state analysis means 30 The system ot claim 7 further comprising brain-wave signal acquisition means, said brain-wave signal acquisition means including filter means tor receiving a representation of said brain-wave signal obtained from the subject, and for substantially attenuating signal components outside a second predefined frequency range and substantially preserving signals in said range above 30 Hz; and conversion means responsive to said brain-wave signal obtained from said subject for converting said brain-wave signal into a sequential plurality of digital data elements representing the amplitude of said brain-wave signal at a corresponding plurality of discrete times. 31 The system of claim 30 wherein said subject state analysis means further comprises. means for selecting for analysis a limited number of said data elements, means for analyzing frequency content of said acquired brain-wave signal, said means being operative to define a plurality of frequencies or frequency bands within said first predefined range, said means also being responsive thereto and to said selected data elements to 22. The system of claim 21 wherein each brain-wave signal-processing channel comprises: frequency selection means responsive to said obtained brain-wave signal for passing substantially exclusively brain-wave signal components within said corresponding range of frequencies; and energy-determination means responsive to an output of said frequency selection means for determining energy contained in said brain-wave signal components within said corresponding range of frequencies during a predefined time interval. 23. The system of claim 22 wherein each brain-wave signal processing channel further comprises: means responsive to an output of said energy-determination means for determining an inverse of said energy contained in said brain-wave signai components within said corresponding range of frequencies. 24. The system of claim 23 wherein each brain-wave signal processing channel further comprises: weighting-function means responsive to an output of said inverse-energy-determining means for applying a predetermined weight to said inverse of said energy. 25. The system of claim 24 wherein said means responsive to said outputs supplied by a plurality of such brain-wave signal-processing channels to produce said output-measure signai comprises a summation means responsive to an output from said weighting function means of each brain-wave signal-processing channel. 26. The system of claim 7 further comprising brain-wave signal acquisition means, said brain-wave signal acquisition means including: filter means for receiving a representation of said brain-wave signal obtained from the subject, and for substantially attenuating signal components above a predefined cutoff frequency and substantially preserving signals in said range above 30 Hz; and conversion means responsive to said brain-wave signal obtained from the subject for producing a digital representation of said brain- wave signal. 27. The system of claim 26 wherein said subject state analysis means comprises: means for analyzing frequency content of said acquired brain- wave signal responsive to said digital representation of said brain-wave signal and said first predefined range of frequencies for defining a plurality of frequencies or frequency bands within said range, and for determining for each of said defined frequencies or frequency bands an amplitude or distribution of amplitudes corresponding to a component or components of said brain-wave signal at such frequency or within such frequency band: determine for each such defined frequency or frequency band an^rrTjJRuTde'or distribution of amplitudes corresponding to a component or components of said acquired brain-wave signal at that frequency or within that frequency band; power determination means for determining for each of said defined frequencies a power represented in said corresponding component or for each of said defined frequency- bands a distribution of power represented in said corresponding components: spectral-grouping means for defining a limited number of spectral groups, each of said spectral groups defining a group of frequencies within said first predefined range and including at least one of said defined frequencies or frequency bands; group-energy-determining means responsive to said power determination means and said spectral-grouping means for determining for each of said spectral groups an aggregate energy represented in a plurality of said components corresponding to the defined frequencies or frequency bands included in such spectral group; and output-measure-determining means responsive to said group-energy-determining means for producing an output signal indicative of said state of the subject. 32. A system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia comprising: means for obtaining a brain-wave signal having rhythmic components from the subject over at least a first predetermined range of frequencies including a second predefined range of frequencies including frequencies above 30 Hz; brain-wave signal acquisition means including: (a) filter means for receiving a representation of said brain-wave signal obtained from said subject, and for substantially attenuating signal components outside a third predefined frequency range and substantially preserving signals within said second predefined subrange and (b) conversion means responsive to a representation of said brain- wave signal obtained from the subject for producing a digital representation of said brain-wave signal; and subject state analysis means responsive to said digital representation of said brain-wave signal for: (a) selecting components of said brain- wave signal lying in said first predetermined range of frequencies including said second predefined range (b) determining a contribution to said brain-wave signal due to said components lying in said second predefined range and (c) producing, responsive to at least said determined contribution, an output- measure signal indicative of said state of the subject. 33. The system of claim 32 wherein said subject state analysis means comprises: means for analyzing frequency content of said acquired brain-wave signai responsive to said first predefined range of frequencies for defining a plurality of frequencies or frequency bands within said range, and for determining for each of said defined frequencies or frequency bands an amplitude or distribution of amplitudes corresponding to said component or components of said acquired brain-wave signai at such frequency or within such frequency band; power determination means for determining for each of said defined frequencies a power represented in said corresponding components or for each of said defined frequency bands a distribution of power represented in said corresponding components; and spectral-grouping means for defining a limited number of spectral groups, each of said spectral groups defining a group of frequencies within said first predefined range and including at least one of said defined frequencies or frequency bands: group-energy-determining means responsive to said power determination means and said spectral-grouping means for determining for each of said spectral groups an aggregate energy represented in a plurality of said components corresponding to the frequencies or frequency bands included in such spectral group; and output-measure-determining means responsive to said group-energy determining means for producing said output-measure signal indicative of said state of the subject. 34. The system of claim 33 wherein said brain-wave signal acquisition means includes a first digital control and processing means and said subject state analysis means includes a second digital control and processing means. 35. The system of claim 33 wherein said system includes a digital control and processing means operatively connected to said brain-wave signal acquisition means and to said subject state analysis means, whereby said digital control and processing means performs at least some processing functions for each of said brain- wave signal acquisition means and said subject state analysis means. 36. The system of claim 32 wherein said conversion means converts said brain-wave signal into a sequential plurality of digital data elements representing amplitude of said brainwave signal at a corresponding plurality of discrete times. 37. The system of claim 36 wherein said subject state analysis means further comprises: means for analyzing frequency content of said acquired brain-wave signal, said means being operative to define a plurality of frequencies or frequency bands within said first predefined range, and said means also being responsive thereto and to said selected data elements to determine for each such defined frequency or frequency band an amplitude or distribution of amplitudes corresponding to a component or components of said acquired brainwave signal at such frequency or within such frequency band; power determination means for determining for each of said frequencies the power represented in said corresponding component or for each of said frequency bands the distribution of power represented in said corresponding components; spectral-grouping means for defining a limited number of spectral groups, each of said spectral groups defining a group of frequencies within said first predefined range and including at least one of said frequencies or frequency bands; group-energy-determining means responsive to said power determination means and said spectral-grouping means for determining for each of said spectral groups an aggregate energy represented in a plurality of said components corresponding to the frequencies or frequency bands included in such spectral group; and output-measure-determining means responsive to said group-energy-determining means for producing an output signal indicative of said state of the subject. 38. For use with a system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia, said system having means for obtaining a brain-wave signal including rhythmic components from the subject, and conversion means responsive to said brain-wave signal for converting the subject's brain-wave signal into a sequential plurality of digital data elements representing amplitude of said brain-wave signal at a corresponding plurality of discrete time intervals, and producing therefrom an output-measure signal indicating said state of the subject; the method of determining a threshold value which defines a range of output-measure signal values indicating that the subject is not fully alert, comprising the steps of: (a) obtaining a sequential plurality of digital data elements representing amplitude of said brain- wave signal of the subject over a period of time, said signal including rhythmic components (b) selecting a threshold-determination segment of said digital data elements obtained in step (a) corresponding to a period during which the subject is in a state of wakefulness; (c) determining a value of said output-measure signal at each of said plurality of discrete time intervals represented in said threshold-determination segment; (d) determining, for each of said plurality of discrete time intervals represented in said threshold-determination segment, whether said output-measure signal is contaminated by artifact during said discrete time interval; (e) discarding from said output-measure-signal values substantially all output- measure-signal values corresponding to said discrete time intervals during which said output-measure-signal value is contaminated by artifact; (f) determining a known-state baseline output-measure value for the subject as a combined value of substantially all undiscarded output-measure-signal values corresponding to said discrete time intervals represented in said threshold-determination segment; (g) multiplying said known-state baseline output-measure value for the subject by a predetermined threshold factor to obtain an output-measure threshold value for the subject corresponding to said selected physiological state. 39. For use with a system for using rhythmic brain-wave signals to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia, said system having means for obtaining a brain-wave signal from said subject, and conversion means responsive to said brain-wave signal for converting said brain-wave signal into a sequential plurality of digital data elements representing amplitude of said brain-wave signal at a corresponding plurality of discrete time intervals, and producing therefrom an output-measure signal indicating said state of said subject; the method of determining a threshold value which defines a range of output-measure-signal values corresponding to a selected physiological state of said subject, comprising the steps of: (a) obtaining a sequential plurality of digital data elements representing amplitude of said brain-wave signal of said subject over a period of time: (b) selecting a threshold-determination segment of said digital data elements obtained in step (a) corresponding to a period during which said subject is in a known state along said continuum; (c) determining a value of said output-measure signal at each of said plurality of discrete time intervals represented in said threshold-determination segment; (d) determining, for each of said plurality of discrete time intervals represented in said threshold-determination segment, whether said output-measure (e) discarding from said output-measure-signal values substantially all output- measure-signal values corresponding to said discrete time intervals during which said output-measure-signal value is contaminated by artifact; (f) determining a known-state baseline output-measure value for said subject as a combined value of substantially all undiscarded output-measure-signal values corresponding to said discrete time intervals represented in said threshold-determination segment; (g) multiplying said known-state baseline output-measure value for said subject by a predetermined threshold factor to obtain an output-measure threshold value for sai l ? subject corresponding to said selected physiological state. 40 For use in analyzing brain-wave signal information acquired from a human subject, said brain-wave signal information being presented as a sequential plurality of digital data elements representing amplitude of the subject ' s brain-wave signal at a corresponding plurality of discrete time intervals, and producing therefrom an output-measure signal indicating the subject ' s state along a continuum relating to alertness, dro siness, sleep, unconsciousness, or anesthesia; the method of determining a threshold value which defines a range of output-measure-signal values indicating that the subject is not fully alert, comprising the steps of (a) a sequential plurality ot digital data elements representing amplitude of a brain-wave signal of the subject over a period of time, said signal including rhythmic components (b) selecting a threshold-determination segment ot said digital data elements obtained in step (a) corresponding to a period during which the subject is in a state of wakefulness; (c) determining a value of said output-measure signal at each of said plurality of discrete time intervals represented in said threshold-determination segment; (d) determining, for each of said plurality of discrete time intervals represented in said threshold-determination segment, whether said output-measure signal is contaminated by artifact during said discrete time interval: (e) discarding from said output-measure-signal values substantially all output- measure-signal values corresponding to said discrete time intervals during which said output-measure-signal value is contaminated by artifact; (f) determining a known-state baseline output-measure value for the subject as a combined value of substantially all undiscarded output-measure-signal values corresponding to said discrete time intervals represented in said threshold-determination segment: and (g) multiplying said known-state baseline output-measure value for the subject by a predetermined threshold factor to obtain a drowsiness output-measure threshold value for the subject. 41. For use in analyzing brain-wave signal information acquired from a human subject, said brain-wave signal information being presented as a sequential plurality of digital data elements representing amplitude of the subject ' s brain-wave signal at a corresponding plurality of discrete time intervals, and producing therefrom an output-measure signal indicating the subject's state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia; the method of determining a threshold value which defines a range of output-measure-signal values corresponding to a selected physiological state of the subject, comprising the steps of: (a) obtaining a sequential plurality of digital data elements representing amplitude of a brain-wave signal of the subject over a period of time, said signal including rhythmic components (b) selecting a threshold-determination segment of said digital data elements obtained in step (a) corresponding to a period during which the subject is in a known state along said continuum; (c) determining a value of said output-measure signal at each of said plurality of discrete time intervals represented in said threshold-determination segment; (d) determining, for each of said plurality of discrete time intervals represented in said threshold-determination segment, whether said output-measure signal is contaminated by artifact during said discrete time interval; (e) discarding from said output-measure-signal values substantially all output- measure-signal values corresponding to said discrete time intervals during which said output-measure-signal value is contaminated by artifact; (f) determining a known-state baseline output-measure value for the subject as a combined value of substantially all undiscarded output-measure-signal values corresponding to said discrete time intervals represented in said threshold-determination segment; and (g) multiplying said known-state baseline output-measure value for the subject by a predetermined threshold factor to obtain an output-measure threshold value for the subject corresponding to said selected physiological state. 42. For use with a system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia, said system having means for obtaining a brain-wave signal including rhythmic components from the subject, and conversion means responsive to said brainwave signal for converting said brainwave signal into a sequential plurality of digital data elements representing amplitude of said brain- wave signal at a corresponding plurality of discrete time intervals, and producing therefrom an output-measure signal indicating said state of the subject; the method of detecting whether said brain-wave signal is contaminated by artifact during a selected time period, comprising the steps of: (a) obtaining an artifact determination segment of data comprising a sequential plurality of digital data elements representing amplitude of selected rhythmic components of said brain-wave signal of the subject over a period of time, said signal including rhythmic components (b) selecting in turn for further processing a plurality of data windows, each comprising a subset of said artifact determination segment corresponding to an adjacent predefined window interval within said period of time; (c) rejecting as saturated and in its entirety any data window containing a digital data element having at least one value equal to the maximum value which may be produced by said conversion means; (d) analyzing frequency content of said brain-wave signal using said data elements of said data window, thereby determining for each of a plurality of defined frequencies or frequency bands an amplitude or distribution of amplitudes corresponding to a component or components of said brain-wave signal at such defined frequency or within such frequency band; (e) determining for each of said defined frequencies power represented in said corresponding components or for each of said defined frequency bands a distribution of power represented in said corresponding components; (0 defining at least one spectral group, each of said spectral groups comprising a group of frequencies which may be present in said brain- wave signal and including at least one of said defined frequencies or frequency bands; (g) determining for each of said spectral groups an aggregate energy represented in a plurality of said components corresponding to the frequencies included in such group; (h) determining a combined value of said aggregate energy present in a plurality of said spectral groups over a plurality of said data windows in said artifact determination segment; (i) multiplying said combined value of said aggregate energy by a predetermined artifact threshold factor to produce an artifact threshold value; and (j) comparing said aggregate energy for each data window to said artifact threshold value, and identifying as contaminated by artifact each window interval in which the corresponding aggregate energy exceeds said artifact threshold value. 43. The method of claim 42 further comprising the step of: (k) inhibiting use of said output-measure signal corresponding to each window interval identified as contaminated by artifact in step (j). 44. The method of claim 42 further comprising the step of: (k) further identifying as contaminated by artifact each window interval in which said corresponding data window was rejected in step (c). 45. The method of claim 42 further comprising the steps of: (k) receiving at least one additional signal, other than said brain-wave signal, providing information concerning whether said brain-wave signal is contaminated by artifact; and (I) further identifying as contaminated by artifact each window interval in which said additional signal indicates said brain-wave signal is contaminated by artifact. 46. For use in analyzing brain-wave signal information acquired from a human subject, said brain-wave signal information being presented as a sequential plurality of digital data elements representing amplitude of the subject ' s brain-wave signal at a corresponding plurality of discrete time intervals, said digital data elements being derived from said brainwave signal through a conversion process, and producing therefrom an output-measure signal indicating the subject ' s state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia: the method of detecting whether said brain-wave signal data corresponding to a selected time period is contaminated by artifact, comprising the steps of: (a) obtaining an artifact determination segment of data comprising a sequential plurality of digital data elements representing amplitude of said brainwave signal of the subject over a period of time, said signal including rhythmic components; (b) selecting in turn for further processing a plurality of data windows, each comprising a subset of said artifact determination segment corresponding to an adjacent predefined window interval within said period of time; (c) rejecting as saturated and in its entirety any data window containing a digital data element having at least one value equal to a maximum value which may be produced in said conversion process; (d) analyzing frequency content of said brain-wave signal using said data elements of said data window, thereby determining for each of a plurality of defined frequencies or frequency bands an amplitude or distribution of amplitudes corresponding to a component or components of said brain-wave signal at such defined frequency or within such frequency band; (e) determining for each of said defined frequencies power represented in said corresponding components or for each of said defined frequency bands a distribution of power represented in said corresponding components; (f) defining at least one spectral group, each of said spectral groups comprising a group of frequencies which may be present in said brain-wave signal and including at least one of said defined frequencies or frequency bands; (g) determining for each of said spectral groups an aggregate energy represented in a plurality of said components corresponding to the frequencies included in such group; (h) determining a combined value of said aggregate energy present in a pluralitv of said spectral groups over a plurality of said data windows in said artifact determination segment; (i) multiplying said combined value of said aggregate energy by a predetermined artifact threshold factor to produce an artifact threshold value; and (j) comparing said aggregate energy for each data window to said artifact threshold value, and identifying as contaminated by artifact each window interval in which the corresponding aggregate energy exceeds said artifact threshold value. 47. For use with a system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia, said system having means for obtaining a brain-wave signal including rhythmic components from the subject, and conversion means responsive to said brain-wave signal for converting said brainwave signal into a sequential plurality of digital data elements representing amplitude of said brain-wave signal at a corresponding plurality of discrete time intervals, and producing therefrom an output-measure signal indicating said state of the subject; the method of producing said output-measure signal, comprising the steps of: (a) obtaining selected rhythmic components of said brain-wave signal from the subject, said signal including rhythmic components (b) defining at least one range of frequencies including frequencies above 30 Hz; (c) for each such defined range, determining a contribution to said brainwave signal due to components of said signal at frequencies lying within that defined range; and (d) determining said output-measure signal responsive to said contributions to said brain- wave signal determined for a plurality of said defined ranges in combination. 48. The method of claim 47 wherein step (b) further comprises the step of defining at least one range including frequencies above 30 Hz. 49. The method of claim 47 wherein step (a) further comprises the steps of: (al) acquiring said brain-wave signal as an electrical signal; (a2) amplifying said electrical signal; (a3) filtering said electrical signal such that signal components outside a second predefined frequency range are substantially attenuated each 1 df said defined ranges are substantially preserved. 50. The method of claim 47 wherein step (a) further comprises the steps of: (a) acquiring said brain-wave signal as an electrical signal; and (b) converting said electrical brain-wave signal into a sequential plurality of digital data elements representing amplitude of said brain-wave signal at a corresponding plurality of discrete times. 51. The method of claim 51 wherein step (c4) further comprises the steps of: (cl ) defining a plurality of data windows, each comprising a subset of said digital data elements corresponding to a predefined window interval: (c2) selecting each of said windows in turn for further processing; (c2) applying a windowing function to adjust values of said digital data elements of said selected window; (c3) analyzing frequency content of said brain-wave signal using said digital data elements of said selected window, thereby determining for each o a plurality of defined frequencies or frequency bands an amplitude or amplitude distribution corresponding to a component or components of said brain-wave signal at such frequency or within such frequency band; (c4) for each of said plurality of said defined frequencies or frequency bands determining the power or distribution of power contributed by said corresponding component or components components of said brain-wave signal; (c5) defining at least one spectral bin and assigning to each of said defined spectral bins at least one of said defined frequencies or frequency bands; and (c6) for each of said spectral bins, determining a total energy contributed by said components of said EEG signal corresponding to a plurality of said defined frequencies or frequency bands assigned to such spectral bin. 52. The method of claim 51 wherein step (c4) further comprises the steps of: (c4i) defining at least one of said defined frequencies or frequency bands as an undesired frequency or frequency band; and (c4ii) for each of said undesired frequencies or frequency bands, whenever further processing would require use of said determined contribution due to a component or components corresponding to such frequency or frequency band, substituting therefor a value derived from said determined contribution due to a component or components corresponding to at least one nearby frequency or frequency band. 53. The method of claim 51 wherein: step (a) further comprises the step of (al) applying an anti-aliasing filter to said brain-wave signal obtained from said subject; and step (c4) further comprises the step of (c4i) adjusting said determined contribution to compensate for an effect of said filtering step in at least one of said subranges. 54. The method of claim 47 wherein step (d) further comprises the steps of: (dl) for each such defined range, determining an inverse of energy contained in said brain-wave signal by signal components corresponding to frequencies within that defined range and (d2) for each such defined range, applying a respective predefined weight to said corresponding inverse-energy value; and (d3) determining a value for said output-measure signal by arithmetically combining said weighted inverse-energy values of a plurality of said defined ranges. 55. For use in analyzing brain-wave signal information acquired from a human subject, said brain-wave signal information being presented as a sequential plurality of digital data elements representing amplitude of said subject ' s brain-wave signal at a corresponding plurality of discrete time intervals, said signal including rhythmic components, said digital data elements being derived from said brain-wave signal through a conversion process, and producing therefrom an output-measure signal indicating the subject's state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia; the method of producing said output-measure signal, comprising the steps of: (a) analyzing frequency content of said brain- wave signal including rhythmic components using said digital data elements, thereby determining for each of a plurality of defined frequencies or frequency bands an amplitude or amplitude distribution corresponding to a component or components of said brain-wave signal at such frequency or within such frequency band; (b) defining at least one range of frequencies including frequencies above 30 Hz; (c) for each such defined range, determining a contribution to said brainwave signal due to components of said signal at frequencies lying within that defined range and (d) determining said output-measure signal responsive to said contributions to said brain-wave signal determined for a plurality of said defined ranges in combination. 56. The method of claim 55 wherein step (b) further comprises the step of defining at least one range including frequencies above 30 Hz. 57. A system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia comprising: means for acquiring an electrical brain-wave signal including rhythmic components from the subject; at least one signal-processing channel means; each such signal-processing channel means being associated with a respective predefined range of frequencies at least one of said signal-processing channel means being associated with a range of frequencies including frequencies above 30 Hz; each such signal-processing channel means having: means for receiving said acquired brain-wave signal; means for substantially limiting response of said signal-processing channel to its associated predefined range of frequencies and means for determining a contribution to said brain-wave signal due to components thereof lying within said associated predefined range of frequencies; and means for receiving from a plurality of such signal-processing channel means respective signals representing said determined contributions and for responsively determining an output- measure signal indicative of said state of the subject. 58. The system of claim 57 wherein said means for substantially limiting response of said signal-processing channel comprises a band-pass filter. 59. The system of claim 57 wherein each such signal-processing channel means further comprises: means for determining power contributed to said brain-wave signal due to components lying within said associated predefined range of frequencies. 60. The system of claim 59 wherein each such signal-processing channel means further comprises: means responsive to said power-determining means for determining energy contributed to said brain-wave signal over a predefined time interval due to components lying within said associated predefined range of frequencies; means responsive to said energy-determining means for determining an inverse value of said energy contribution; and weighting means responsive to said inverse-value-determining means for applying a weighting function to said determined inverse value. 61. The system of claim 60 wherein said energy-determining means comprises an integrator. 62. The system of claim 60 wherein said weighting means comprises means for establishing a weighting parameter, means for receiving a signal representing said inverse- energy value, and multiplier means for multiplying said signal by said inverse-energy value. 63. The system of claim 57 wherein each such signal-processing channel means further comprises: means for determining energy contributed to said brain-wave signal over a predefined time interval due to components lying within said associated predefined range of frequencies. 64. The system of claim 57 wherein said output-measure-determining means comprises a summer. 65. A system adapted to measure in a human subject a state along a continuum relating to alertness, drowsiness, sleep, unconsciousness, or anesthesia comprising: means for acquiring a brain-wave signal including rhythmic components from the subject; means for selecting components of said brain-wave signal lying in a predetermined range of frequencies including frequencies above 30 Hz; and means responsive to components of said signal corresponding to frequencies lying in said range for producing an output-measure signal indicative of said state of the subject; said means for producing said output-measure signal being adapted to modify its response when artifact is present in said brain-wave signal. 66. The system of claim 65 wherein: said means for producing said output-measure signal causes said output-measure signal to change in a first direction when said state of the subject approaches alertness; and said means for producing said output-measure signal causes said output measure signal to change in said first direction responsive to presence of artifact in said brain- wave signal. 67. The system of claim 65 wherein: said means for producing said output-measure signal causes said output- measure signal to change in a first direction when said state of the subject deviates from alertness; and said means for producing said output-measure signal causes said output-measure signal to change in a direction opposite from said first direction responsive to presence of artifact in said brain-wave signal.

Inventors: KAPLAN FREDERIC RICHARD, LOPARO KENNETH ALAN, KAPLAN, FREDERIC, RICHARD, LOPARO, KENNETH, ALAN
Applicants: CONSOLIDATED RESEARCH OF RICHM, CONSOLIDATED RESEARCH OF RICHMOND, INC.
Application year: 1998
Main classification: A61N5
IPC codes: A61B5