Method and system for recognizing workload of a person

0001The invention relates to a method and system for recognizing and/or detecting workload of a person, especially of a driver of a vehicle, by detecting and evaluating head movements of the person, and for recognizing and/or detecting the resulting impairment and/or inattention of the person due to a related distraction of the person caused by the workload.

0002The vast majority of traffic accidents - about 90% according to an American study (Treat, 1979) - occurs due to driver error. Driver inattention is the most prevalent primary cause of collisions, accounting for an estimated 25-55% of crashes (according to an NHTSA study of the CDS database representing over 2,5 million tow-away crashes per year - Wang, et al, 1996). Substantially, at least one of the following three components contribute to driver inattention: drowsiness, distraction and/or workload. Common crash types and locations caused by inattention are rear-end crashes, crashes at intersections, at lane change/merge locations, at road departures, as well as single vehicle crashes or crashes on low speed limit roadways.

0003More in detail, collisions caused by drowsiness are a very serious road safety problem. 57% of drivers who have been asked, reported that they have driven while being drowsy and 23% said that they had actually fallen asleep while driving (John-son, 1998). Drivers cannot reliably tell when they are impaired by fatigue to the point of having a serious vigilance lapse or uncontrolled sleep attack (Wylie, et al, 1996; Brown, 1997). Low vigilance and later stages of drowsiness impairs alertness and judgment, and hence driving performance (Williamson et al, 2001). Fatigue has been estimated to be involved in between 2% and 23% of all crashes (Knipling & Wang, 1995). These crashes occur more often late at night or in mid-afternoon. They typically involve a single vehicle drifting out of the lane or road at high speeds with no indications that the driver attempted to avoid this (NHTSA, 1998). The injuries in case of a crash are often very serious. Sleepy drivers are also often involved in rear-end and head-on collisions (op. cit.). These statistics are most likely underestimates of the real problem since drivers tend to be unwilling to recognize or admit that sleepiness was a factor in causing their crash (Stutts, Wilkens & Vaughn, 1999).

0004Furthermore, drivers often attend to things that are not related to driving. These things are often called secondary tasks and are potential reasons for distraction. Unfortunately, the driver is often unaware of the effects that distraction has on the driving task (Boase, et al, 1988; Stern, et al, 1996). As a result, distraction (including visual distraction and mental distraction, i.e. looking without seeing) is found to account for at least 22.9% of crashes (Wang, et al, 1996). Since 1996 there has been a rapid increase and market penetration of new information systems and functionality that may be used within the vehicle, e.g. mobile telephones, navigation aids, internet, e-services. These distracters may capture a driver's attention excessively and thus potentially increase collisions based on distraction causes.

0005Generally, a secondary task becomes a distraction (including visual-, auditory-, mental-, cognitive-, and biomechanical distraction) when the driver's attention is captured by information that is irrelevant to the driving situation to the degree that: <ol id="ol0001" compact="compact"><li>(a) insufficient attention is left for the primary control task of driving, and/or</li><li>(b) driving performance (e.g. lane keeping or speed control) and safety are compromised.</li></ol>

0006Finally, workload refers to how busy a driver is and the amount of effort the driver needs to perform the tasks at hand. Having a lot to do ("high workload") creates much to attend to (i.e. a high perceptual demand) and much to act upon (i.e. a high action demand). For example, driver workload increases both with the primary driving task demand (caused e.g. by road and environment conditions, traffic, driving style, etc.) as well as with the execution of secondary tasks (distraction) as mentioned above. Thus, high workload situations increase the probability of drivers becoming distracted.

0007<patcit id="pcit0001" dnum="WO0160254A1"><text>WO 01/60254 A1</text></patcit> discloses a method and means for monitoring driver alertness by imparting an impulse to the steering or some other part of the vehicle with which the driver actively interacts, whilst underway, by means of one or more actuators. The spontaneous reaction of the driver to this impulse is detected by one or more sensors and is compared with the impulse itself in order to draw conclusions regarding the driver's presence of mind. However, imparting an impulse to any part of the vehicle is considered as detrimental. Furthermore, the driver's alertness cannot be monitored continuously but only at such a moment in which the impulse is imparted.

0008Another approach to the determination of driver alertness is a "Driver alertness indication system" (DAISY) of Spericon Ltd. as disclosed in the article "Principles of operation and system feasibility" dated June 2002. By this system, the accumulated action of external forces exerted on the vehicle and disturbing the movement of the vehicle along its path is detected, and the action of the driver who acts to maintain the vehicle on the road is analyzed. The analysis of the driving pattern vis-à-vis the pattern of the disturbances allows a continuous determination of the driver level of alertness. However, the necessary separation between the movements initiated by the driver (which originate at the steering wheel) and the movements due to the external disturbances (which originate at the road wheel) has to be evaluated in consideration of the free play and the time lag in the steering system of the vehicle so that an intensive use of mathematical and engineering techniques as well as the application of complex algorithms is required.

0009<patcit id="pcit0002" dnum="USPS5691693A"><text>US-PS 5,691,693</text></patcit> discloses an "Impaired transportation vehicle operator system" and a related method in which the head position and head movement of the driver is sensed by means of three capacitive coupling sensors. The output signals of the sensors are analyzed by a microprocessor which records, tracks and compares the operator head movement data with stored profiles of normal head movement and profiles of impaired head movement to detect any head movement indicative of operator impairment. However, this system and method is not able to detect and quantify what types of head movement are part of the normal and are part of the impaired profiles.

0010It is an object of the invention to provide a method and system for detecting and/or recognizing workload of a person when monitoring a process, e.g. a vehicle driver when driving a vehicle, which workload might lead to reduced safety of driving.

0011It is another object of the invention to provide a method and system for quantifying the level of workload of a person, especially a vehicle driver.

0012These objects are solved by a method according to claims 1 and by a system according to claim 6.

0013The dependent claims disclose advantageous embodiments of the methods and system according to claims 1 and 6, respectively.

0014Further details, features and advantages of the invention are evident from the following description of exemplary embodiments of the invention in connection with the drawings, in which schematically shows: <dl id="dl0001" compact="compact"><dt>Fig. 1</dt><dd>a schematic view of normal regions and outside normal regions of head movements;</dd><dt>Fig. 2</dt><dd>a schematic view of sub-regions within the outer normal regions of head movements;</dd><dt>Fig. 3</dt><dd>a plot of the density of horizontal and vertical head rotations of alert driv- ers;</dd><dt>Fig. 4</dt><dd>a plot of the density of horizontal and vertical head rotations of drowsy drivers;</dd><dt>Fig. 5</dt><dd>a first power spectral density plot of vertical head rotations of alert drivers in comparison to drowsy drivers;</dd><dt>Fig. 6</dt><dd>a second power spectral density plot of vertical head positions of alert drivers in comparison to drowsy drivers;</dd><dt>Fig. 7</dt><dd>a third power spectral density plot of the horizontal (depth of) head posi- tions of alert drivers in comparison to drowsy drivers;</dd><dt>Fig. 8</dt><dd>a schematic flow chart of a first embodiment of a method according to the invention; and</dd><dt>Fig. 9</dt><dd>a schematic flow chart of a second embodiment of a method according to the invention.</dd></dl>

0015In the following, the invention shall first be explained in the form of preferred methods and systems for recognizing drowsiness.

0016The basis of this inventive methods and systems is that muscles fall into a relaxed state when a person becomes drowsy. Such drowsiness can be caused e.g. by fatigue and/or alcohol and/or drugs and/or medicine and/or intoxication etc.. As a background, skeletal muscle tone results entirely from nerve impulses coming from the spinal cord. These in turn are controlled partly by impulses transmitted from the brain and partly by impulses that originate in muscle spindles located in the muscle itself. The muscle spindle impulses provide the nervous system with continuous feedback of information from each muscle, giving the status of each muscle at each instant, which is length, tension, and how rapidly length or tension of the muscle is changing. Thus, there is a general relaxing of muscle tone with drowsiness. As a result of loss of muscle tone especially in the neck and of sluggishness in reaction to disturbances, the head starts to move around more and with larger amplitudes of movements. This change of head movements can be quantified by the methods and systems according to the invention.

0017The following measures of drowsiness may all be used as absolute or relative values. As absolute values, the measures are taken to indicate a level of drowsiness on a scale from alert to asleep. As relative values, the measures are used to indicate the difference when comparing to a driving segment in which the driver is alert. Since head movements may differ from person to person, the method can have an initiation phase where normal or alert behavior is recorded and/or analyzed and quantified and used for comparison. The degree of deviation from this normal or alert behavior then indicates drowsiness.

0018Generally, it has to be distinguished between detecting drowsiness from head movements alone (first embodiment of the invention) and between detecting drowsi-ness by comparing head movements with a disturbance (second embodiment of the invention).

0019The first embodiment mentioned in the preceding paragraph shall now be explained in more details.

0020The head movement is generally described by three head rotation components (Pitch, Yaw and Roll), and three head translation components according to a Cartesian coordinate system with the axis x, y and z (Pos_x, Pos_y, Pos_z). These six components can be combined to a head movement signal <i>HMOVE</i> that completely describes the head movement.

0021Preferably, each of the six signal components is pre-processed to remove measurement noise and long term postural changes of the driver. One way to achieve this is to use a high-pass filter to remove the long term postural changes, which by nature are of low frequency, and a low-pass filter to suppress or at least attenuate the signal noise. Both can of course be achieved by an appropriate band-pass filter as well.

0022After this pre-processing, the global measurement of the total head movement <i>HMOVE</i> can e.g. be expressed and calculated as a function of the six degrees of freedom and their first- and second-degree time-dependent derivatives according to equation (1.1): <maths id="math0001" num="(1.1)"><math display="block"><mi mathvariant="italic">HMOVE</mi><mo>=</mo><mo>ℑ</mo><mfenced><msub><mi mathvariant="italic">Pos</mi><mi>X</mi></msub><mo>,</mo><msub><mi mathvariant="italic">Pos</mi><mi>Y</mi></msub><mo>,</mo><msub><mi mathvariant="italic">Pos</mi><mi>Z</mi></msub><mo>,</mo><mi mathvariant="italic">Pitch</mi><mo mathvariant="italic">,</mo><mi mathvariant="italic">Yaw</mi><mo mathvariant="italic">,</mo><mi mathvariant="italic">Roll</mi><mo>,</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>⁢</mo><msub><mi mathvariant="italic">Pos</mi><mi>X</mi></msub><mo>,</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>⁢</mo><msub><mi mathvariant="italic">Pos</mi><mi>Y</mi></msub><mo>,</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>⁢</mo><msub><mi mathvariant="italic">Pos</mi><mi>Z</mi></msub><mo>,</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>⁢</mo><mi mathvariant="italic">Pitch</mi><mo>,</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>⁢</mo><mi mathvariant="italic">Yaw</mi><mo mathvariant="italic">,</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>⁢</mo><mi mathvariant="italic">Roll</mi><mo>,</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>⁢</mo><msub><mi mathvariant="italic">Pos</mi><mi>X</mi></msub><mo>,</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>⁢</mo><msub><mi mathvariant="italic">Pos</mi><mi>Y</mi></msub><mo>,</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>⁢</mo><msub><mi mathvariant="italic">Pos</mi><mi>Z</mi></msub><mo>,</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>⁢</mo><mi mathvariant="italic">Pitch</mi><mo>,</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>⁢</mo><mi mathvariant="italic">Yaw</mi><mo mathvariant="italic">,</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>⁢</mo><mi mathvariant="italic">Roll</mi></mfenced></math><img file="EP2072010A2_D0001.tif" /></maths>

0023This general function or model can be chosen and tuned differently for different applications. For example, a global head movement measure could be defined according to equation (1.2): <maths id="math0002" num="(1.2)"><math display="block"><mi mathvariant="italic">HMove</mi><mo>=</mo><msqrt><mi>A</mi><mo>*</mo><msup><mi mathvariant="italic">pitch</mi><mn>2</mn></msup><mo>+</mo><mi>B</mi><mo>*</mo><msup><mi mathvariant="italic">yaw</mi><mn>2</mn></msup><mo>+</mo><mi>C</mi><mo>*</mo><msup><mi mathvariant="italic">roll</mi><mn>2</mn></msup><mo>+</mo><mi>D</mi><mo>*</mo><msup><msub><mi mathvariant="italic">Pos</mi><mi>x</mi></msub><mn>2</mn></msup><mo>+</mo><mi>E</mi><mo>*</mo><msup><msub><mi mathvariant="italic">Pos</mi><mi>y</mi></msub><mn>2</mn></msup><mo>+</mo><mi>F</mi><mo>*</mo><msup><msub><mi mathvariant="italic">Pos</mi><mi>z</mi></msub><mn>2</mn></msup></msqrt></math><img file="EP2072010A2_D0002.tif" /></maths> wherein A to F are weighting factors which determine the sensitivity for different types of head movements and thus the sensitivity of the measurement.

0024For detecting drowsiness from this head movement signal <i>HMOVE,</i> the inventive methods and systems are based on and make use of a (first) reference calculation of a normal region center point (NRCP) X or of a normal region (NR). Reference is made to <figref idref="f0001">Figure 1</figref> which shows such a NRCP which is surrounded by a NR. Furthermore, a first, a second and a third outer normal region ONR1, ONR2, ONR3 are indicated in the form of concentric circles surrounding the NR. (Although <figref idref="f0001">Figure 1</figref> shows a two-dimensional area (e.g. pitch and yaw movements), this can be applied and extended to e.g. all six dimensions as mentioned above).

0025At first, the NRCP as shown in <figref idref="f0001">Figure 1</figref> is calculated for each of the head rotation and head translation signal components measured as mentioned above during the first few minutes (e.g. three minutes) of driving or the first few minutes (e.g. three minutes) of driving in a certain environment. Exemplary results of such calculations are shown in <figref idref="f0002">Figures 3 and 4</figref>.

0026<figref idref="f0002">Figure 3</figref> shows a plot of the resulting density of signal components related to horizontal and vertical head rotations (as an example) of a plurality of alert drivers, whereas <figref idref="f0002">Figure 4</figref> shows this density for a plurality of drowsy drivers. On the x-axis of both these Figures, signal components of vertical rotations (pitch movements of the head in degree) and on the y-axis signal components of horizontal rotations (yaw movements of the head in degree) are indicated, whereas on the z-axis the density of these movements in percent per bin or funnel is indicated. The NRCP according to <figref idref="f0001">Figures 1 and 2</figref> corresponds with the point at values x=0 and y=0.

0027<figref idref="f0002">Figure 3</figref> shows the result of such a first reference calculation for head rotation signals of 12 alert drivers during preferably the first three minutes of their respective driving. In comparison to this, <figref idref="f0002">Figure 4</figref> shows the result of second calculations (as explained below) for 12 drivers in case that they are drowsy.

0028For both alert and drowsy drivers there is a distinct peak of the density of the signal components which is considered as the NRCP (X). By comparing <figref idref="f0002">Figures 3 and 4</figref> it is obvious that finding the NRCP does not rely on the driver being alert because the NRCP results independently of the impairment and/or inattention of the driver.

0029At second, it is evaluated how a change in the distribution of the density of the signal components between the different regions in <figref idref="f0002">Figures 3 and 4</figref> could be classified as drowsiness. In order to compare a current distribution with another distribution, a basic trough for each individual driver is defined. The basic trough is established on the basis of head movements during the first few minutes of driving (or the first few minutes of driving in a particular environment). Thus, for a relative measure on the basis of the basic trough it has to be assumed that the driver is alert and not under the influence of drugs during these first few minutes of calibration. For an absolute measure, a threshold is pre-set and no calibration period is needed.

0030Generally, the basic trough can be defined either on the basis of a plurality of alert drivers so that it can be used as a driver-independently standard for any driver, or the basic trough is established individually for each specific (alert) driver who is monitored by the related system

0031The sensor used for these measurements is e.g. a head- and eye-tracker based on one or two cameras and a computer software for processing the output signals of the camera(s). The sensor especially measures head position and direction in three-dimensional coordinates. Of course, other measurement techniques or other sensors (e.g. infrared sensors) can be used as well.

0032In order to calculate the distribution of head movements within each region (NRCP and each NR and ONR e.g. according to <figref idref="f0001">Figure 1</figref>), each signal component (in each head movement signal) is counted in its corresponding region. The more signal components (e.g. Pos_x and Pos_y) are inside one region, the higher the value will be for this region and the lower for the other regions (all regions sum up to 100 % for each head movement signal). If the regions of <figref idref="f0001">Figure 1</figref> are laid into the x/y-plane of <figref idref="f0002">Figures 3 and 4</figref> they will indicate the amount of horizontal and vertical head movements within each of the regions (most movements are in the normal region NR) for these two signal components. From <figref idref="f0002">Figures 3 and 4</figref> it is obvious as well that the distribution between the regions differs considerably for alert and drowsy drivers.

0033The regions have the form of a grid (rectangular, circular according to <figref idref="f0002">Figures 3 and 4</figref>, or they have another shape) of bins/funnels and the distribution is a normalized histogram for the signal components based on the sizes of the bins/funnels. For example it is assumed that for a particular signal component a first region between 0 and 1, a second region between 1 and 2 and a third region between 2 and 3 is defined and five samples of the signal component having values of 0,5; 1,1; 2,4; 0,9 and 0,3 are considered. The resulting distribution would be 60 % of the values within region 1 and each 20 % in region 2 and 3, respectively. However, according to the invention there are up to six-dimensional histograms/distributions which are evaluated in this way and it is assumed that certain "shapes" of these distributions indicate drowsiness.

0034For each head rotation and head translation signal component, each sample is added to a related (rotation- and translation-) funnel wherein a funnel size of e.g. 0.98 by 0.98 degrees (128 x 128 for +/-30 degrees from straight ahead, or the zero point) can be used. Then, the mode of each funnel (samples with the largest frequency in a funnel) is set as a normal region center point (NRCP) X for each rotation and translation signal components, respectively, and stored.

0035Subsequently, the above first (reference) calculation (by which the NRCP or the NR has been calculated) is repeated preferably periodically with a predetermined time distance as a plurality of second calculations. Each result of the subsequent second calculations is compared with the first calculation. If the difference between both exceeds a predetermined value, it is supposed that the driver is drowsy and a signal is output or an alarm device is activated. It is noted that it is the distribution which is continuously calculated. For this purpose a sliding time window is used that allows to focus e.g. on the past minute or another time interval of driving and disregard what happened before this time interval. The continuously calculated distribution is then compared to the reference distribution (or a pre-defined "alert"-distribution). When the difference exceeds a predefined threshold value or when the continuously calculated distribution is deviating from a predefined reference shape or has a certain predefined shape it is assumed that the driver is drowsy.

0036In this disclosure, a "sliding time window" is a window which is sliding along a signal. As a new sample is added, the oldest sample in the time window will be deleted. All measurements are calculated based on the samples within the sliding time window.

0037As shown in <figref idref="f0001">Fig. 2</figref>, the outer normal regions ONR1, ONR2, ONR3 can be further divided into sub-regions 2 to 9, 10 to 17, and 18 to 25, respectively (region 1 is considered as the normal region NR). By this, the number of funnels/bins and their shape and size can be different from the division shown in <figref idref="f0001">Figure 1</figref>. This might be advantageous for evaluating and comparing certain distributions with respect to certain signal components (e.g. in the case of <figref idref="f0002">Figures 3 and 4</figref> with respect to the horizontal and vertical head rotations).

0038For evaluating a percentage normal region (PNR) measure, at first each sample within a sliding time window (e.g. a one minute time window with a 60Hz update frequency) is classified as being either a "1", i.e. lying within a normal region (NR) or a "0", i.e. lying outside the normal region (ONR) on the basis of a normal region criteria. The normal region criteria is, for example, calculated by taking the distance from normal region center point X and setting a cutoff threshold, e.g. 5 degrees for rotational position of the head and/or 5 cm for translational position of the head as a radius around the normal region center point X, indicating e.g. a first normal region NR as indicated in <figref idref="f0001">Figure 1</figref>.

0039Of course several normal region criteria with different thresholds can be established leading to a second normal region NR2, a third normal region NR3, a fourth normal region NR4, etc., each having a radius which corresponds to the related threshold value.

0040The cutoff threshold(s) could for example also be defined as data within one standard deviation from normal region center point X, or else wise identified by clustering algorithms, thereby allowing normal region(s) NR to be non-circular as well.

0041Again, those samples that fall within the cutoff threshold are classified as "within the related normal region (NR)" and those that fall outside are classified as "outside the related normal region (ONR)".

0042In the next step, a calculation of percentage is made by dividing the number of data-points within the related normal region NR by the total number of data-points within the time window and e.g. multiplying the product by 100 so that a percentage normal region measure (PNR) is achieved.

0043Besides the percentage normal region measure (PNR), a Standard Deviation of distance to normal region center point (SD_DNRC) from mode road center (SD-MRC) - can be calculated in the same time window as mentioned above, according to the standard deviation formula, with the exception that the mean value is replaced by mode:

0044More in detail, the standard deviation (SD) of the distance between the NRCP and the current signal value (called DNRC) is calculated. However, the usual mean value used in the SD-formula has to be replaced with the mode i.e. the NRCP. The calculation is performed on the limited amount of data contained within a sliding time window. Thus it becomes a RMS (Root Mean Square) of the DNRC. A threshold at e.g. two times the SD of an alert driver is set, and when the threshold is reached or exceeded, it its assumed that the driver is drowsy: <maths id="math0003" num="(1.3)"><math display="block"><mi mathvariant="italic">DNRC</mi><mo>=</mo><msqrt><mrow><mo>(</mo><mi mathvariant="italic">HMove</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">NRCP</mi><mo>⁢</mo><msup><mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mspace width="7em" /><mi mathvariant="italic">HMOVE</mi><mo mathvariant="italic">,</mo><mi mathvariant="italic">NRCP</mi><mo>∈</mo><msup><mi>R</mi><mi>N</mi></msup></math><img file="EP2072010A2_D0003.tif" /></maths><maths id="math0004" num="(1.4)"><math display="block"><mi mathvariant="italic">SD_DNRC</mi><mo>=</mo><msqrt><mfrac><mstyle displaystyle="false"><munder><mo>∑</mo><mspace width="1em" /></munder><msup><mi mathvariant="italic">DNRC</mi><mn>2</mn></msup></mstyle><mi mathvariant="italic">WindowLenght</mi></mfrac></msqrt><mspace width="4em" /><mi mathvariant="italic">DNRC</mi><mo>∈</mo><mi>R</mi></math><img file="EP2072010A2_D0004.tif" /></maths> R^N in formula (1.3) are real numbers (not complex numbers) of dimension N, wherein N = 1, 2, 3, 4, 5, 6, 7, 8,......depending on how many signal components were used to form <i>HMOVE</i> and <i>NRCP.</i> Consequently, R (or R¹) in formula (1.4) are real numbers of one dimension.

0045In equation (1.3), <i>HMOVE</i> could e.g. be defined according to formula (1.5): <maths id="math0005" num="(1.5)"><math display="block"><mi mathvariant="italic">HMOVE</mi><mo>=</mo><mfenced><mi>A</mi><mo>*</mo><msub><mi mathvariant="italic">Pos</mi><mi>X</mi></msub><mo>,</mo><mi>B</mi><mo>*</mo><msub><mi mathvariant="italic">Pos</mi><mi>Y</mi></msub><mo>,</mo><mi>C</mi><mo>*</mo><msub><mi mathvariant="italic">Pos</mi><mi>Z</mi></msub><mo>,</mo><mi>D</mi><mo>*</mo><mi mathvariant="italic">Pitch</mi><mo>,</mo><mi>E</mi><mo>*</mo><mi mathvariant="italic">Yaw</mi><mo>,</mo><mi>F</mi><mo>*</mo><mi mathvariant="italic">Roll</mi></mfenced><mo>⁢</mo><mi mathvariant="italic">HMOVE</mi><mo>∈</mo><msup><mi>r</mi><mn>6</mn></msup></math><img file="EP2072010A2_D0005.tif" /></maths> wherein A to F are weighting factors.

0046Drowsiness it then calculated based on pattern recognition (e.g. neural net or statistical) of the distribution of head movements between the normal regions NR; NR1, NR2,...and/or outer normal regions ONR1, ONR2,... and the segments 2, 3,... 25. An alert driver will have a different distribution in the regions according to <figref idref="f0001">Figure 1 or 2</figref> compared to a drowsy driver.

0047In more details, for each instant in time and for each driver there will be a unique distribution of head movements in the different regions NR, ONR. Each distribution indicates how the driver has moved his head during the past X minutes (depending on the length of the sliding time window). It is assumed that the distributions of drowsy drivers differ considerably from alert drivers (again, see <figref idref="f0002">Figures 3 and 4</figref>). In order to classify a driver as having the distribution of a drowsy driver, pattern recognition methods are used. According to <figref idref="f0002">Figure 4</figref> the drowsy drivers have a much more smeared out distribution of the density of the signal components compared to the alert drivers according to <figref idref="f0002">Figure 3</figref>.

0048In a very simple example, drowsiness could be defined as having a predefined value X % of the head movements outside a predefined circle (which in <figref idref="f0002">Figures 3 and 4</figref> is placed in the x/y-plane). However, <figref idref="f0002">Figures 3 and 4</figref> only show two dimensions, but the distribution and recognition can be done in N-dimensions (signal components), analyzing e.g. both head position and head rotation, as well as velocity or angular velocity of the head movement.

0049<figref idref="f0003">Figure 5</figref> shows a plot of the power spectral density (indicated on the y-axis) of vertical head rotations (pitch movements) of drowsy drivers (curve B) in comparison to alert drivers (curve A) in terms of the frequency in Hertz (indicated on the x-axis).

0050<figref idref="f0003">Figure 6</figref> shows such a plot of the power spectral density for vertical head positions (translation movement) of drowsy drivers (curve B) in comparison to alert drivers (curve A) whereas <figref idref="f0004">Figure 7</figref> shows such a plot of the power spectral density for horizontal (depth of) head positions (translation movement), preferably comprising movement-components in a back-forward direction, again of drowsy drivers (curve B) in comparison to alert drivers (curve A). The y-axis in <figref idref="f0003">Figures 6</figref> and <figref idref="f0004">7</figref> is again denoted with the power spectral density and the x-axis of <figref idref="f0003">Figures 5 and 6</figref> is denoted with the frequency in Hertz.

0051<figref idref="f0003 f0004">Figures 5 to 7</figref> show a general increase in head movements of drowsy drivers. More in detail, curve A represents the frequency content of head movements of 12 alert drivers driving for 20 minutes each. In <figref idref="f0003 f0004">Figures 5 to 7</figref> the alert drivers show less head movements across all frequencies.

0052The drowsy driver data have been measured in blocks of 20 minutes each (e.g. drowsy block (DB) 1 represents the first 20 min of all 12 drowsy drivers, drowsy block (DB) 2 represents the second 20 min, etc). The trend in data from such DB 1 to a DB 4 (which in this example begins 40 minutes after the end of DB 1 and lasts for 20 minutes) shows that the drowsy drivers are effected by time-on-task. It has been measured that there is a general increase in head movements from DB1 to DB3 because drivers are driving accordingly longer. This trend, however, is reversed sometimes in that the DB4 has a level slightly below DB3, but yet above DB2. The reason for this could be that the drivers expect that the trip would be ending and thus became slightly more alert towards the end of the trip. The curves B in <figref idref="f0003 f0004">Figures 5 to 7</figref> represent the power spectral density curves as a sum of all such drowsiness blocks 1 to 4 (DB1 to DB4). The curves of these drowsiness blocks DB1 to DB4 are indicated in <figref idref="f0003 f0004">Figures 5 to 7</figref> with dashed lines.

0053Presenting the data in DBs provides an illustration of the trend of increasing head movements with increasing drowsiness. The indicated curves B show the power spectral density content corresponding with the same data as shown in <figref idref="f0002">Figures 3 and 4</figref>.

0054These power spectral density data provide a strong evidence showing that there is indeed an increase in head movements as drivers become increasingly drowsy.

0055Another alternative is to measure the power spectral density data on a real-time basis so that a real-time measurement of drowsiness can be conducted. For this purpose, the power spectral density data could be calculated in real-time either across the frequencies e.g. shown in <figref idref="f0003 f0004">Figures 5 to 7</figref> or a subset of these frequencies.

0056The second embodiment of a method and system according to the invention is proposed for detecting drowsiness by comparing head movement reactions in case of a disturbance and shall now be explained in more details.

0057Generally, the drowsier a driver is the less control he has of his corrective head movement, that means in this case head movements which try to counteract a disturbance. The disturbance could be side wind, rough pavement or simply movements induced by the vehicle itself (common in trucks). A basic principle is to measure the disturbance and the corrective head movement of the driver. The difference between these signals (the disturbance and the reaction) which is the head movement reaction (HMR) measure is an estimate of how well the driver compensates for the disturbance. An alert driver has faster, more precise, and cushioned reactions whereas a drowsy driver exhibits slower, less precise and sloppy reactions with higher amplitudes of dislocations ("rubberneck"). Consequently, the difference between disturbance and reaction is greater for drowsy drivers.

0058First, the disturbance itself has to be measured. This could be done in a number of ways (or combinations of these) like e.g. the following: <ul id="ul0001" list-style="dash" compact="compact"><li>Measuring dynamic properties of objects in a camera view, for example body movements or loose objects relative to the background;</li><li>Measuring the force acting on the driver with an acceleration sensor (for example placed in the seat);</li><li>Measuring the cabin air-suspension and tilt angle compared to the rest of the truck;</li><li>Measuring steering wheel movements.</li></ul>

0059The reactions to the disturbances are contained in head movements (head movements can be either stimulus-driven [reactions] or goal-driven [intentional]). Head movement reactions (HMR) to disturbances can be identified and filtered out in a number of ways as for example: <ul id="ul0002" list-style="dash" compact="compact"><li>A threshold trigger can be set on the disturbance signal. Reactions are only calculated for disturbances which are larger than a certain threshold. Thus, the measurement is probe based (the probe being the disturbance) and not a continuous signal;</li><li>Head movements could only be used if the vehicle is traveling fairly straight for a predefined period of time and above a predefined speed threshold, i.e. not while turning onto a new road, going around a roundabout or while in a city;</li><li>Head movements should be ignored when signaling to change lanes.</li></ul>

0060Finally, the combination of GPS- and map-matching data and/or statistical pattern recognition data can be used to define the environment and to evaluate whether the vehicle is traveling or is going to travel fairly straight and when disturbances can be used or have to be used and when not.

0061As mentioned above, the HMR measure is a function of disturbance and reaction: <maths id="math0006" num="(2.1)"><math display="block"><mi mathvariant="italic">HMR</mi><mo>=</mo><mo>ℑ</mo><mfenced><mi mathvariant="italic">Distubance</mi><mo>⁢</mo><mi mathvariant="italic">Reaction</mi></mfenced></math><img file="EP2072010A2_D0006.tif" /></maths>

0062The disturbance is a force acting on the driver and the reaction is the response reaction of the driver with respect to his head movements. The HMR function could for example be defined to be: <ul id="ul0003" list-style="dash" compact="compact"><li>The difference (e.g. with respect to displacement) between the reaction and the normal region center point (NRCP) X. The difference is larger for drowsy drivers than for alert drivers;</li><li>The RMS (root mean square) value of the difference between the reaction and the normal region center point X. Again, a low RMS value indicates that head movement is being controlled well and the driver is alert. A high RMS value indicates a drowsy driver;</li><li>Pattern recognition of HMR;</li><li>The error between a reference model for a head movement in response to a disturbance and the measured head movement in response to the same disturbance. The error represents the amount of drowsiness. The alert driver head movement response to disturbances could for example be a physical model of a particle mass on top oftwo rods connected with an upper joint and a lower joint at the bottom, which is a simple model of the head and spine. The lower joint represents the seat contact with the driver, the lower rod represents the body trunk, the upper joint represents the neck, and the particle mass represents the head. Disturbances in the driver environment (truck cab) are measured and fed into the model which result in modeled head movements, which are then compared with measured head movements. The model is tuned to react as the driver when she/he is alert. Models such as those of crash test dummies could be used.</li></ul>

0063According to a third embodiment of a method and system according to the invention, distraction of the driver can be recognized by applying the above methods in the following way: <ul id="ul0004" list-style="none" compact="compact"><li>Especially in a truck, a bus, a train, but even in a smaller car, interaction with integrated systems within the cab or cabin sometimes demands large eye movements (back and forth between road center and system) which are often accompanied by smaller head movements back and forth in the same direction. This is also true for side mirror checks and over the shoulder checks. These head movements could be identified and classified as distractive tasks and detected and evaluated as explained above.</li></ul>

0064According to a fourth embodiment of a method and system according to the invention, workload of the driver can be recognized by applying the above methods considering the fact that workload can be quantified in terms of the amount of head movement. The following two equations exemplify this: <maths id="math0007" num="(4.1)"><math display="block"><mi mathvariant="italic">SD</mi><mfenced><mi mathvariant="italic">HMOVE</mi></mfenced></math><img file="EP2072010A2_D0007.tif" /></maths><maths id="math0008" num="(4.2)"><math display="block"><mi mathvariant="italic">RMS</mi><mfenced><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mfenced><mi mathvariant="italic">HMOVE</mi></mfenced></mfenced></math><img file="EP2072010A2_D0008.tif" /></maths><i>SD</i> stands for standard deviation, <i>RMS</i> for the root mean square and <i>HMOVE</i> can be defined as exemplified in equation (1.2).

0065Equations (4.1) and (4.2) are calculated over a sliding time window of a few (e.g. 4) seconds. These measures produce a continuous quantification of workload. A threshold can be used to obtain preferably a Boolean value, i.e. high or low workload, and a high-to-low workload delay is used to prevent flickering of the value. The use of high weighting factors for the absolute values (equation 1.1) makes the estimation more sensitive to simultaneous changes in several of the terms that are used to describe head movement, for example directional movement, towards the center stack, which consists of both pitch and yaw movements, and less sensitive to head movements caused by the road curvature, which are mainly yaw movements. The dynamics of the signal can be tuned by adjusting the length of the time window and/or the magnitude of the weighting factors. Instead of one threshold and a Boolean value, a plurality of thresholds could be used to obtain several levels of workload accordingly.

0066The time-dependent derivatives of head position and head rotation could also be used in a similar way to detect workload. For example, fast head movements would indicate a higher workload than slow head movements.

0067The invention is not only applicable to vehicle drivers but also to airplane pilots, air traffic controllers, engine drivers or any persons who have to be aware of not falling asleep and sit fairly still and look most of her time in the same direction.

0068All of the measures described in this disclosure can be individualized by saving individual specific data (e.g. distributions of densities of certain head movement signal components for an alert condition) on a driver identification chip, for example a smart card.

0069<figref idref="f0005">Figure 8</figref> shows a schematic flow chart of the first embodiment of the method according to the invention for detecting drowsiness simply by comparing an alert distribution of signal-densities (or a pre-stored distribution of such densities) with a current distribution of the related signal-densities of head movements.

0070The components in <figref idref="f0005">Figure 8</figref> are denoted as follows: <ul id="ul0005" list-style="none" compact="compact"><li>Arrow H: Input of driver head movement signals (all six components as mentioned above); <ul id="ul0006" list-style="none" compact="compact"><li>1: Pre-processing stage including noise suppression and removal of postural effects;</li><li>2: Initialization and adaptation block for initializing the system parameters that are specific to the driver (this may include a continuous adaptation as well);</li><li>3: Head movement distribution calculation for the alert initialization period (or adaptation of the pre-defined head movement distribution of the alert driver). This step is paused after a predefined period of time and serves then as a the reference distribution for an alert driver;</li></ul></li><li>Arrow E: Enables the drowsiness calculation when initialization has been finished; <ul id="ul0007" list-style="none" compact="compact"><li>4: Continuous calculation of head movement signal distributions based on the pre-processed head movement signals of the driver delivered by component 1 for the last predefined X minutes of driving ("sliding time window");</li><li>5: Comparing the current head movement signal distribution with the head movement signal distribution of the alert driver. This component contains the criteria and / or thresholds for when the driver is to be considered drowsy; Arrow O: Output of continuous or discrete drowsiness values.</li></ul></li></ul>

0071<figref idref="f0005">Figure 9</figref> shows a schematic flow chart of the second embodiment of the method according to the invention for detecting drowsiness by comparing head movement reactions in case of a disturbance.

0072The components in <figref idref="f0005">Figure 9</figref> are denoted as follows: <ul id="ul0008" list-style="none" compact="compact"><li>Arrow H: Input of driver head movement signals (all six components as mentioned above); <ul id="ul0009" list-style="none" compact="compact"><li>1: Pre-processing stage including noise suppression and removal of postural effects;</li><li>2: Initialization and adaptation block for initializing the system parameters that are specific to the driver (this may include a continuous adaptation as well);</li></ul></li><li>Arrow E: Enables the drowsiness calculation when initialization has been finished;</li><li>Arrow D: Input of a disturbance signal (e.g. accelerometer signal); <ul id="ul0010" list-style="none" compact="compact"><li>3: Probe trigger and/or simple reference model of the driver. It is decided which disturbances are large enough (and exceed a predefined value) to be used and which are feed through to the comparison component 5. The component 3 includes a simple reference model of the driver (e.g. the mechanical model described above) as well. Disturbances D are fed into the model and the output (simulated head movements) is then fed through to the comparison component 5;</li></ul></li><li>Arrow P: Trigger signal, indicates when to compare and/or calculate measures; <ul id="ul0011" list-style="none" compact="compact"><li>4: Continuous calculation of head movement measures (including head movement distribution) for the last predefined X minutes of driving;</li><li>5: Comparing the current head movement measures with the modeled head movement distribution. This component contains the criteria and / or thresholds for when the driver is to be considered drowsy based on this comparison and the values of head movement- and disturbance measures;</li></ul></li><li>Arrow O: Output of continuous or discrete drowsiness values.</li></ul>

0073Summarizing the above, the following methods and systems are provided: <ul id="ul0012" list-style="none" compact="compact"><li>Generally, a method and system for recognizing and/or detecting impairment and/or inattention of a person by detecting and evaluating head movements of the person on the basis of at least one of three head rotation components and/or three head translation components is provided, comprising the following steps: <ul id="ul0013" list-style="dash" compact="compact"><li>conducting a first calculation for evaluating a reference density pattern of the at least one component during a first time interval in which the person has a reference level of impairment and/or inattention;</li><li>conducting at least one second calculation for evaluating an actual density pattern of the at least one component during at least one second time interval; and</li><li>comparing the actual density pattern with the reference density pattern, evaluating the difference between both and outputting a signal if the difference exceeds a predefined value.</li></ul></li></ul>

0074The reference density pattern is preferably evaluated in the form of a normal region or a normal region center point including the greatest density of samples of the measured component.

0075Preferably one or a plurality of thresholds is predetermined for defining one or a plurality of normal regions, respectively.

0076The thresholds are preferably defined with respect to a standard deviation from a normal region center point.

0077The comparing of the actual density pattern with the reference density pattern preferably includes evaluating the difference between both in different normal regions and/or segments of at least one normal region.

0078Generally, another method and system is provided for recognizing and/or detecting impairment and/or inattention of a person, especially a driver of a vehicle, who is actingly or reactingly connected with an object, especially a vehicle, wherein the object or vehicle is influenced or affected by an external disturbance which causes a head movement reaction of the person, comprising the following steps: <ul id="ul0014" list-style="dash" compact="compact"><li>detecting the disturbance exerted onto the object or vehicle by external forces;</li><li>detecting head movement reaction of the person upon the disturbance on the basis of at least one of three head rotation components and/or three head translation components;</li><li>evaluating the difference between the head movement reaction and a predefined reference head movement reaction or a stored previous head movement reaction upon a previous disturbance and</li><li>outputting a signal if the difference exceeds a predefined value.</li></ul>

0079Preferably, the object is a vehicle and the detection of the disturbance includes measuring of at least one of the following values: <ul id="ul0015" list-style="dash" compact="compact"><li>dynamic properties of objects within the vehicle;</li><li>an acceleration force acting on the person;</li><li>the cabin air-suspension and tilt angle compared to the rest of the vehicle, especially in case of a truck; and</li><li>steering wheel movements.</li></ul>

0080Preferably, detecting the head movement reaction includes at least one of the following measures and evaluations, respectively: <ul id="ul0016" list-style="dash" compact="compact"><li>predefining a threshold value for the disturbance signal and evaluating head movement reactions if the disturbance exceeds the threshold value;</li><li>evaluating head movements if the vehicle is traveling substantially in a straight direction for substantially a predefined period of time and above a predefined speed threshold value; and</li><li>ignoring head movements if the vehicle is not traveling substantially in a straight direction for substantially the predefined period of time and/or above the predefined speed threshold value.</li></ul>

0081The traveling of the vehicle substantially in a straight direction and/or the speed of the vehicle is preferably evaluated on the basis of GPS- and map-matching data and/or similar navigation/localization systems and data.

0082Preferably the detecting of the head movement reaction is activated and/or deactivated on the basis of GPS- and map-matching data and/or similar navigation/localization systems and data.

0083The evaluation of the difference preferably includes at least one of the following measures: <ul id="ul0017" list-style="dash" compact="compact"><li>the difference between the reaction and the normal region center point;</li><li>the RMS (root mean square) value of the difference between the reaction and the normal region center point; and</li><li>the error between a model for a head movement in response to a disturbance and the measured head movement in response to the same disturbance.</li></ul>

0084The evaluation of the head movements is preferably conducted for detecting drowsiness and/or distraction and/or workload of a person.

0085In this case, recognizing and/or detecting workload of a person is preferably conducted by detecting and evaluating head movements of the person on the basis of the amount and/or the time-dependent derivative of at least one of three head rotation components and/or three head translation components.

0086The above methods are preferably conducted by means of a computer program comprising computer program code means adapted to perform the method when said program is run on a programmable microcomputer.

0087Such a computer program can be adapted to be downloaded to the above systems or one of their components when run on a computer which is connected to the internet.

0088Finally, a computer program product stored on a computer usable medium can be provided comprising the above computer program code means.