[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Extraction of Rice Phenological Differences under Heavy Metal Stress Using EVI Time-Series from HJ-1A/B Data
Next Article in Special Issue
Indoor Trajectory Tracking Scheme Based on Delaunay Triangulation and Heuristic Information in Wireless Sensor Networks
Previous Article in Journal
The Lateral Tracking Control for the Intelligent Vehicle Based on Adaptive PID Neural Network
Previous Article in Special Issue
An Improved Compressive Sensing and Received Signal Strength-Based Target Localization Algorithm with Unknown Target Population for Wireless Local Area Networks
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vital Sign Monitoring and Mobile Phone Usage Detection Using IR-UWB Radar for Intended Use in Car Crash Prevention

Department of Electronics and Computer Engineering, Hanyang University, 222 Wangsimini-ro, Seongdong-gu, Seoul 04763, Korea
*
Author to whom correspondence should be addressed.
Co-first authors, these authors contributed equally to this work.
Sensors 2017, 17(6), 1240; https://doi.org/10.3390/s17061240
Submission received: 23 March 2017 / Revised: 20 May 2017 / Accepted: 25 May 2017 / Published: 30 May 2017
Figure 1
<p>Noisy vital signal obtained by conventional algorithm during motion period.</p> ">
Figure 2
<p>“Vital signal” vs slow time and fast time.</p> ">
Figure 3
<p>Process block diagram for extraction of driver vital signs.</p> ">
Figure 4
<p>Data fitting algorithm results (<b>a</b>) fitting frequency higher than the respiration frequencies range (<b>b</b>) fitting frequency lower than the respiration frequencies range (<b>c</b>) fitting magnitude less than vital signal threshold (<b>d</b>) fitting signal for signal having information of vital signs.</p> ">
Figure 4 Cont.
<p>Data fitting algorithm results (<b>a</b>) fitting frequency higher than the respiration frequencies range (<b>b</b>) fitting frequency lower than the respiration frequencies range (<b>c</b>) fitting magnitude less than vital signal threshold (<b>d</b>) fitting signal for signal having information of vital signs.</p> ">
Figure 5
<p>R square values while (<b>a</b>) stationary (<b>b</b>) moving hand near abdominal region (<b>c</b>) moving upper parts of body slightly.</p> ">
Figure 5 Cont.
<p>R square values while (<b>a</b>) stationary (<b>b</b>) moving hand near abdominal region (<b>c</b>) moving upper parts of body slightly.</p> ">
Figure 6
<p>A best fit signal which has the corresponding r-square value of 0.46.</p> ">
Figure 7
<p>The vital signal construction from the candidate sub-signal components of the best fit signals.</p> ">
Figure 8
<p>The first sub-signal component of vital_constructed_signal.</p> ">
Figure 9
<p>The constructed vital signal after second iteration.</p> ">
Figure 10
<p>The final constructed vital signal.</p> ">
Figure 11
<p>State diagram of dual-mode background subtraction algorithm for phone usage detection.</p> ">
Figure 12
<p>Signal characteristics when using a mobile phone (<b>a</b>) received signal over fast time (<b>b</b>) channel impulse response (<b>c</b>) the maximum value of the received signal over slow time.</p> ">
Figure 13
<p>Signal characteristics when there is a background change (<b>a</b>) received signal over fast time (<b>b</b>) channel impulse response.</p> ">
Figure 14
<p>(<b>a</b>) Experimental setup inside the car (<b>b</b>) IR-UWB radar module.</p> ">
Figure 15
<p>(<b>a</b>) Measurement setup for the ECG signal acquisition; (<b>b</b>) The signal obtained from the ECG sensor module PSL-iECG2 (measurement reference device).</p> ">
Figure 16
<p>(<b>a</b>) Comparison of the vital signal obtained from the conventional vs proposed algorithm while the human is slightly moving (<b>b</b>) R-square and normalized variance values at each fast-time index during motion state.</p> ">
Figure 16 Cont.
<p>(<b>a</b>) Comparison of the vital signal obtained from the conventional vs proposed algorithm while the human is slightly moving (<b>b</b>) R-square and normalized variance values at each fast-time index during motion state.</p> ">
Figure 17
<p>Spectrum obtained by applying FFT algorithm to the signal obtained by proposed algorithm in <a href="#sensors-17-01240-f016" class="html-fig">Figure 16</a>a.</p> ">
Figure 18
<p>The R-square values in case of huge body movements (all the values are below the threshold value of 0.3).</p> ">
Figure 19
<p>The maximum value of the received signal over slow time when using a mobile phone (<b>a</b>) texting (<b>b</b>) scrolling, touching (<b>c</b>) viewing.</p> ">
Figure 20
<p>The maximum value of the received signal over slow time with momentary movements. (<b>a</b>) moving hand (<b>b</b>) moving water bottle (<b>c</b>) moving cell phone.</p> ">
Figure 20 Cont.
<p>The maximum value of the received signal over slow time with momentary movements. (<b>a</b>) moving hand (<b>b</b>) moving water bottle (<b>c</b>) moving cell phone.</p> ">
Figure 21
<p>The maximum value of the received over slow time when there is a background change (<b>a</b>) water bottle (appearing) (<b>b</b>) water bottle (disappearing)</p> ">
Versions Notes

Abstract

:
In order to avoid car crashes, active safety systems are becoming more and more important. Many crashes are caused due to driver drowsiness or mobile phone usage. Detecting the drowsiness of the driver is very important for the safety of a car. Monitoring of vital signs such as respiration rate and heart rate is important to determine the occurrence of driver drowsiness. In this paper, robust vital signs monitoring through impulse radio ultra-wideband (IR-UWB) radar is discussed. We propose a new algorithm that can estimate the vital signs even if there is motion caused by the driving activities. We analyzed the whole fast time vital detection region and found the signals at those fast time locations that have useful information related to the vital signals. We segmented those signals into sub-signals and then constructed the desired vital signal using the correlation method. In this way, the vital signs of the driver can be monitored noninvasively, which can be used by researchers to detect the drowsiness of the driver which is related to the vital signs i.e., respiration and heart rate. In addition, texting on a mobile phone during driving may cause visual, manual or cognitive distraction of the driver. In order to reduce accidents caused by a distracted driver, we proposed an algorithm that can detect perfectly a driver's mobile phone usage even if there are various motions of the driver in the car or changes in background objects. These novel techniques, which monitor vital signs associated with drowsiness and detect phone usage before a driver makes a mistake, may be very helpful in developing techniques for preventing a car crash.

1. Introduction

Every year, according to the statistics, more than a million people die on the world’s roads due to car crashes and the cost of dealing with the consequences of these crashes runs to billions [1]. Researchers are trying to reduce these numbers by employing different techniques [2,3,4,5,6,7,8]. Many people are saved during collisions by air bags, seat belts, and the car frames that can absorb a huge amount of energy during a car crash [9]. However, the ultimate solution is to keep cars from smashing into each other in the first place. Adaptive cruise control (ACC) systems use laser beams or radar to measure the distance between the cars and their relative speed [10]. All of the ACC systems available today are built around sensors that detect the vehicle ahead with either radar or light detecting and ranging (LIDAR), the laser-based analog to radar [11]. These systems do detect and prevent some types of driver error. However, technologies cannot detect and prevent all of the mistakes that drivers make and can result in collisions. The major driver mistakes are caused by drowsiness [12] and mobile phone usage during driving.
According to a recent study, 20% of crashes occurred due to drowsiness of the driver [13]. Continuous monotonous driving for a long time causes drowsiness, which may result in fatal road accidents. For long haul/distance drivers, fatigue and drowsiness can hardly be avoided [13]. Thus, research works are carried out to investigate approaches to monitor driver condition. Various methods are used to detect if a driver is tired or falling asleep. Researchers have attempted to determine driver drowsiness using the following measures: (1) vehicle-based measures [14,15] (2) behavioral measures [16,17,18] and (3) physiological measures [19,20,21,22,23]. Mercedes-Benz pioneered one of the first techniques, which uses a computer algorithm that compares the steering behavior of a driver to the trained behavior at the start of the trip [9]. Other systems monitor the location of car within its lane of travel, looking for erratic maneuvers indicative of inattention. Some researchers also track the driver eye movement with an in-car camera, noting rapid or prolonged eye blinks. In [14,15], the researchers used metrics like lane position, steering wheel position, driving speed and car yaw angle for detecting the drowsiness of drivers. These metrics are useful, but this method is an alarm for the mistake after the happening, so it cannot prevent in advance the potential mistakes caused by drowsiness. An eye-tracking based scheme is proposed for detection of driver drowsiness using an unscented Kalman filter [16]. In reference [17], the researchers used facial image sequences to detect the driver fatigue based on multiscale dynamic features, but camera-based techniques can’t perform in dark environments. Biomedical variables like respiration rate related to the nervous system may provide direct information of driver physiological conditions. Hence, they are extremely important to determine the drowsiness cycle and anticipate risky situations during driving [22,23,24]. In [25], the tracking and localization of cyclists using UWB technology is discussed, however, the researchers did not consider the crashes due to cars and other big vehicles and it has no method to detect the mobile phone usage or drowsiness detection of the drivers. Electroencephalogram (EEG) was used by researchers [18] to detect driver drowsiness, but in case of EEG measurement electrodes are used for the measurements which require physical contact with the person. The researchers in reference [19] used EEG, electrocochleogram (EOG), and electrocardiogram (ECG) signals for detection of driver drowsiness based on the physiological conditions of the driver, but this solution is really uncomfortable as the driver has to physically wear the devices. The ANSWatch device was used for detection of heart rate in reference [20], but it requires the driver to always wear the watch. Researchers in [21,22,23] also used EEG and/or ECG for physiological condition monitoring which is cumbersome for the driver. Vital sign measurement for drivers using radar can solve all the above problems at the same time. The radar operates robustly in the dark, and by measuring the vital signs in a noncontact manner without any inconvenience to the driver, an alarm can be generated in advance before the driver makes a mistake due to drowsiness.
In our work, we monitored the vital signs, i.e., respiration and heart rate of the driver inside a real car using IR-UWB radar. The UWB technology offer several benefits such as accurate ranging, robustness to multipath interference, better penetration and rejection of interferences [26,27]. Although there are algorithms in the literature [28,29,30,31] that deal with vital sign measurement, however, the conventional algorithms do not estimate the vital signs during the motion state. In our paper, we present a novel algorithm that can extract the vital signal from the radar signals reflected from the human body which are slightly distorted by the driver motions due to steering use or electronic and mechanical devices operation. The main concept of our vital sign estimation algorithm is that the vital signals are reflected from many points of the chest and abdominal area of the human sitting in front of the radar. Although the driving activity decreases the signal to noise ratio (SNR) of some part of the vital signal reflected from some area of the body, but still it doesn’t decrease the SNR of all the signals reflected from the body. Our goal is to search for those segments of the signals which have the best vital signs information and which are least distorted by the motion caused by driver activity. To this end, we first determined those signals which have better sinusoidal shapes and hence better vital signal information by using the sinusoidal fitting algorithm. Then we segment all those sinusoidal signals into sub-signal components and select those sub-signal components which have better correlation with the vital signal. In this way, we construct our desired vital signal from many signals reflected from the driver body. Then FFT algorithm is applied to the constructed vital signal to find the respiration and heart rate values.
Another main cause of collisions is mobile phone use during driving. Distracted driving caused due to the usage of cell phone (texting or speaking), is associated with 26% of all crashes and is increasing in frequency [32]. The mobile phone usage during driving results in restriction of sight; limiting drivers’ ability to monitor the road since their line of vision is focused on the handset [33]. It also reduces the concentration level and situational awareness [34,35,36]. Motivated by its impact on public safety and property, several state and federal governments have enacted regulations that prohibit driver mobile phone usage while driving [35]. In [36] a computer vision-based method for detection of driver cell phone usage by using a near infrared (NIR) camera system inside the car was presented. Smith et al. [37] have presented an algorithm for the detection of handset usage by the driver by analyzing images taken inside the car. Xu et al. [38] have proposed a machine-learning-based method for detecting driver cell phone usage using a camera system directed at the vehicle’s front windshield. These methods, however, are dependent on a camera, which may not provide good performance in dark environments. In our work, we proposed an algorithm to detect the use of mobile phones while driving by using IR-UWB radar that is unaffected by the light conditions. The proposed algorithm can accurately detect the use of mobile phones even in environments where driver's various movements or change of background objects occur inside the car. We propose a dual mode background subtraction algorithm using two different clutter signals for optimum detection, which is explained in Section 3.2.
The main contribution of our work is that we propose an algorithm using a single IR-UWB radar for monitoring the vital signs and detecting the mobile phone usage of the driver that can be used to develop technologies to prevent accidents caused by drowsiness driving and cell phone use. Our main purpose of the monitoring of driver’s vital signs is to provide a method for future researchers to use it for drowsiness detection so that the accidents due to drowsiness may be reduced. We propose an algorithm to accurately measure vital signs during the motion caused by the driving activities. We also suggest a way to detect the use of a cell phone that takes the driver's attention away. The proposed algorithm suggests a new method to optimally detect drivers’ cell phone usage considering various situations that can occur in the car. The technologies presented in this paper may become one of the technologies to save human life by being utilized in application technologies for preventing drowsiness driving and cell phone use during driving that are the root cause of many road crashes.
In Section 2 of the paper, the problems related to the previous methods for vital measurement and mobile phone detection are discussed. The signal preprocessing such as clutter removal is also discussed in Section 2. Section 3 is about our proposed algorithms for vital sign measurement during various driving conditions and phone detection in a specific area inside a car. Section 4 is related to the experimental results and in Section 5 the paper is concluded.

2. Problem Description and Signal Preprocessing

2.1. Problem Statement

In this paper, we cover two topics to prevent car crashes. One is to measure the driver's vital signs to prevent drowsy driving and the other is to detect the use of mobile phones to prevent accidents caused by their use while driving.

2.1.1. Monitoring of Vital Signs

Although there exist algorithms in the literature to measure the vital signs of a human [28,29,30,31], vital sign measurement during driving is a totally different problem as it involves different hand and body movements. The vital signs obtained during motion are deteriorated by the motion and vital signs can’t be measured from the signal. An example of such a signal is shown in Figure 1. The detection of the vital signs is difficult due to the movement of the driver as there is no trace of any periodicity in the received signal, which is a characteristic of vital signs. Researchers in [28] have discussed the problem of vital sign measurement of a non-stationary human, but the approach is very simple as they detect the motion of the body and lock the vital sign measurement until the body motion is stopped. In this work, however, we not only detect the body motion but also measure the vital signs during the motion period. The conventional method of extraction of vital signs at a fast time index of maximum variance doesn’t always hold true when random body movement is involved, therefore, we have used the data fitting method (see Section 3.1.1) to find those signals located at fast time indexes which gives better sinusoid fits. The vital signals are constructed from those better fit signals which exist even in presence of motion of some parts of the human body. After the vital signal construction, the conventional Fourier transform technique is applied to find the breathing and the heart rate.

2.1.2. Detection of Phone Usage

A cellular phone are composed of electronic parts, PCB, and antenna. Their main components are metals such as copper, so cell phones have a very large radar cross-section (RCS). Although a large RCS makes it easy to detect a cell phone, there are some difficulties in detecting the use of a cell phone in the car. First, in the mobile phone sensing area of a car, various objects besides a mobile phone can appear. For example, several movements such as hand movements to press a button or grab beverage bottles to drink can occur, and the radar can detect it. In fact, these short time movements may be necessary to drive and those are not dangerous actions that distract the driver, therefore any detection algorithm should ignore them.
Secondly, in a really dangerous situation, like when the driver is staring at his cell phone for a while, there is little movement of the cell phone. If only the conventional background removal method [39] in Section 2.2 is used to detect the mobile phone in this situation, the mobile phone may be recognized as a background object and may not be detected. Considering the above two situations, an algorithm for detecting the use of mobile phones in a car should ignore very short time actions that are not dangerous to the driver and generate an alarm for any type of mobile phone use which is dangerous to the driver.
For a momentary movement that does not interfere with the driving, the alarm can be disabled by measuring the duration of the detection. In other words, we can ignore the detection of movement that appears for less than a certain time. In addition, to solve the problem that the signal reflected from the mobile phone is removed as background noise when there is almost no movement of the mobile phone, such as watching the screen quietly, we proposed a dual-mode background subtraction algorithm in Section 3.2. In this algorithm, background removal is performed in two ways, and mobile phone usage is detected by using two received signals obtained by each background removal method.
In the dual-mode background subtraction algorithm, immediately after the handset is detected, the clutter signal is not updated and the background is removed by using the clutter immediately before sensing the cell phone. In this way, it is possible to detect mobile phones with minimal movement. However, when the mobile phone is detected in this way, the mobile phone detection alarm will continue to be sounded when a new background object is placed in the detection area. To solve this problem, we have noted the difference between the received signal from the new background object and the received signal from usage of the mobile phone. There is a slight movement, even if mobile phone is held static like when just looking at a screen. Therefore, when the background is removed while updating the clutter signal, the received signal is small but is detected. On the other hand, the new background object has no motion, so if the background is removed while updating the clutter signal, the received signal is not detected at all. By using these differences, it is possible to distinguish a new background object from a mobile phone.

2.2. Signal Pre-Processing

The loopback filter is applied to remove the clutter from the raw signal [39]. The loopback filter is represented by the following equations:
c k ( t ) = c k 1 ( t ) + ( 1 ) r k ( t )
y k ( t ) = r k ( t ) c k ( t )
In Equations (1) and (2), the symbol is a constant used for weighting. The   value in Equation (2) used for background subtraction was set to 0.97 for vital signs monitoring whereas for the mobile phone detection it was kept 0.8. A low alpha value means that the background removal rate is fast but vulnerable to noise, and if the alpha value is large, the background removal rate is slow but robust against noise. We have chosen the alpha value which is optimal for the corresponding application through several experiments. The symbol r k ( t ) show the received signal and symbol c k ( t ) shows the clutter signal, which is made until the k-th received sample. The background-subtracted signal is represented by the symbol   y k ( t ) . We need to store each filtered signal waveform and combine them into matrix Wmn of size “ m × n . The “m” represents the slow time length whereas the “n” represents the fast time length of the matrix. The value of “n” depends on the observation distance and has a value of 256 for our experiments (which means 1 meter observation distance) while “m” is a user choice and it depends on how long is the observation time. The concept of the slow time and fast time is shown in Figure 2.
Moreover, we only use specified areas for the detection of vital signs and phone detection by extracting only the fast time data of a specific area. The sensing area for phone detection has 25 cm range whereas sensing area for vital signs monitoring has a range of around 50 cm. This helps ignore environmental noise that comes from outside the sensing area. In other words, all the human activities that are out of this range are not taken into account during application of the vital sign measurement and the phone detection algorithm.

3. Proposed Algorithm

3.1. Vital Sign Measurement

The process block diagram for extracting the vital signs from the radar scans is shown in Figure 3. During driving, the hands continuously move because the driver has to control the steering, brakes and other devices inside the car. The conventional algorithm for vital sign measurement of a stationary human supposes that the biggest movement caused by our body is due to the vital signals of breathing and human heart, which is true for the case of stationary humans, but in case of body motions this is not always correct, so conventional algorithms for finding the breathing and heart rate will not provide accurate measurements in such conditions.
We observed that during driving, the major part of the motion exists in the upper part of the driver body and the lower part (abdominal area) is relatively stationary. Therefore, we first find which part of the body has motion due to the vital signals and which part of the body has motion due to random hand gestures, steering control or brake application. To this end, we use the sinusoidal data fitting algorithm, explained below.

3.1.1. Extracting Receive Signals Which Have Useful Vital Sign Information

After we combine the waveforms into a matrix, we need to find which fast-time domain samples represent the sinusoidal motion, which is caused by the breathing and human heart motion. We use sinusoidal fitting algorithm (Algorithm 1).to show how much the received data fit into the sinusoid [40,41]. The input signal is s [ n ] and we want to estimate the frequency, amplitude and phase shift of this signal.
Algorithm 1: Sinusoidal data fitting algorithm for parameters estimation
  • The signal s [ n ]   in its general form can be represented as sum of sinusoids as follows:
    s [ n ] = l = 1 p ( α l sin ( ω l n + φ l ) + C l )
    In Equation (3), α l   is real valued constant describing the amplitude and ω l   represent angular frequency, φ l   is the initial phase. The constant C l   shows the mean value other than zero. It is convenient to write the Equation (3) as follows.
    s [ n ] = l = 1 p ( A l cos ( ω l n ) + B l sin ( ω l n ) + C l )
    It is convenient to use a vector representation of Equation (4). Let’s stack the samples in column vector as follows:
    s = l = 1 p H l θ l
    In the above Equation (5);
    H l = [ cos ( ω l )     sin ( ω l )   1 cos ( 2 ω l )   sin ( 2 ω l )   1 . . . cos ( N ω l )   sin ( N ω l )   1 ]   θ l = [ A l B l C l ]
  • In our case, the model given by Equation (4) is simplified as; p = 1   and C = 0 . Consider a signal s[n] with unknown parameters A , B , C and ω i.e., Equation (4) with   p = 1 . The measured signal is then given by s [ n ] deteriorated by additive noise   w [ n ] as shown in Equation (7):
    x [ n ] = s [ n ] + w [ n ]   n = 1 , 2 , 3 ,   ..   ,   N
    The noise is assumed zero mean white Gaussian noise. The estimation problem is to estimate the signal parameters using the measured samples of the input data x = [ x ( 1 ) ,   x ( 2 ) ,   .   x ( N ) ] .  
  • The probability density function (pdf)
    p ( x ; θ , ω ) = 1 ( 2 π σ 2 ) exp [ 1 2 σ 2 ( x H θ ) T ( x H θ ) ]
    The above equation describes the probability per infinitesimal volume of receiving the data samples x given a set of parameters   { θ , ω } .
  • The maximum likelihood estimator (MLE) tries to maximize the pdf with respect to unknown parameters for given values of x and use those parameters as estimates i.e.,
    [ θ ^ , ω ^ ] = arg   θ , ω m a x p ( x ; θ , ω )
  • Finally, the estimate of θ is given by least-squares solution as follows
    θ ^ = ( H T H ) 1 H T x  
  • By using the least-squares solution, the MLE criterion function can be concentrated to one parameter as follows
    g ( w ) = x T H ( H T H ) 1 H T x
  • The frequency estimate is then obtained from maximizing the function g(w), that is
    ω ^ = a r g ω m a x g ( w )
    The Equation (12) can be solved either by using iterative step method i.e., Gauss-Newton iteration [42] or by a non-linear search.
The input data used for sinusoidal fitting is the slow time domain data at each fast time index. The R-square value is used for finding the fit of the signal, which is defined as follows:
R 2 = 1 i = 1 n ( y i y ^ i ) 2 i = 1 n ( y i y ¯ ) 2
In Equation (13), “ y ^ ” represent the estimated values of “y” by the fitting algorithm whereas “ y ¯ ” shows the mean of “y” [43]. R square is a statistic that give information about the goodness of fit of a model. In general, a model fits the data well if the differences between the observed values and the model’s predicted values are small and unbiased. It determines how well the regression line or curve approximates the real data points. An R square value of 1 indicates that the regression line or curve perfectly fits the real data points. R-square can take on any value between 0 and 1, with a value closer to 1 indicating that a greater proportion of variance is accounted for by the model. For example, an R-square value of 0.73 means that the fit explains 73% of the total variation in the data about the average. The higher value of R square means that the sinusoidal regression model is more accurate and hence the motion is more sinusoidal and hence caused by the respiration and heart rate signals whereas the lower value of R square means non-sinusoidal motion caused by any random body/hand motion. We observed the R square values for all the fast time indexes. We have set the threshold value of 0.3 for finding the best fit sinusoids that contain useful information about the respiration and heart rate. If the value is kept too small then it may also select the signals that doesn’t have good vital signal information because they would not match sinusoidal data and if it is kept too high then we may not get any signal above the threshold because the vital signal is not perfectly sinusoidal and moreover, the movement of the body causes the signal to noise ratio of the vital signal to be very low. Therefore, we chose the optimal threshold value of 0.3 for finding the best fit signals. The detailed algorithm for finding the best sinusoid fits is given in Algorithm 2:
Algorithm 2: Selection of the fast time indexes that gives the best fit signals
  • Find the sinusoid fit as shown in Algorithm 1, which gives three parameters i.e., magnitude, frequency and phase shift of the sinusoid that best fits the signal.
  • Find the R square values for the signals at each fast time index.
  • If the fitting frequency of a signal is not in the range of the respiration frequencies i.e., 0.15 Hz to 0.5 Hz, as shown in Figure 4a,b then replace its R square value by zero.
  • If the fitting sinusoid magnitude is less than a certain threshold value as shown in Figure 4c then replace its R square value by zero. Figure 4d shows the fitting for the signal that have the vital signs information as it has relatively higher magnitude and its frequency lies in the breathing frequencies range i.e., 0.15 Hz to 0.5 Hz.
  • Apply the moving averaging filter to the resulting R square values as shown in Figure 5.
  • Extract received signals with R values above threshold value.
Figure 5a shows the result of the fitting algorithm for the case when the person sitting in front of the radar is stationary. Figure 5b shows the R square values for the case when we move hand in front of the abdominal region, which results in lower R square values at the location of the abdominal region while the chest region signal gives better R square values. In experimental scenario for the Figure 5c, the upper parts of the body are moving due to hand gesture and/or steering motion exerted by the hands. From the comparison of these three figures, we know that the R square values decrease at the points (abdominal area in Figure 5b and chest area in Figure 5c due to the motion of the body parts, even though the variance of the signal increases. If we use the conventional algorithm in which the signal at the highest variance point is considered as vital signal, then it will result in wrong measurements when applied in a condition when the maximum variance is caused by the random body motion instead of the vital signs.
After finding the R square values of the signals in fast time index, we chose signals which have the good R-square values and ignore signals which have small R-square values. Figure 6 shows an example of the selected best fit signals which contain useful information about the vital sign signal. The reason that we expressed the slow time in term of index values rather than seconds is that Algorithm 3 uses the slow time index information for the construction of vital signal.

3.1.2. Estimate of Breathing and Heart Rate through the Vital Signal Reconstruction

Now we have to construct the vital signal from the best-fit signals. First, we divide each signal into sub-signals based on the zero crossings as shown in Figure 6. A sub-signal component of a signal is defined as the signal samples between two consecutive zero crossing values. After dividing the best-fit signals into segments (sub-signal components), we constructed whole vital signal recursively from those segments using correlation concept as explained in Algorithm 3. For the sake of clarity, we have inserted some examples enclosed in the rectangular boxes inside the Algorithm 3.
We constructed the vital signal from the segmented good-fit signals by combining each sub-signal component to the previously constructed vital signal recursively and finding its autocorrelation. We chose the optimum sub-signal component that shows the highest autocorrelation. We have determined the autocorrelation of the signal after combining the new sub-signal-component with the previously constructed vital signal. If these signals have good correlation then the width of the autocorrelation ( w i d t h _ c o r r e l a t i o n ) will be wider, otherwise the width of the autocorrelation will be narrower [28].
In Figure 7, the big picture of the Algorithm 3 is explained by a diagram. The algorithm explains how to construct the vital signal from multiple best fit signals reflected from the body during motion of the driver.
The step by step explanation of the vital signal construction is explained by the following algorithm (Algorithm 3).
Algorithm 3: Step by step construction of the vital signal from the best fit signals
(1) Initialization
Find “k” best-fit signals whose R-square values are above certain set threshold, say R-square_min = 0.3.The size of each best-fit signal is 256 samples. So, we have a matrix of b e s t _ f i t m a t r i x   ( k × 256 ) . In our example, k = 9.
(2) for i t e r a t i o n s = 1 : l a s t _ z e r o _ c r o s s i n g
If (iteration == 1)
Then Find the first zero-crossing fast time index ( z e r o _ c r o s s i n g s ) of all the “k” best-fit signals. And find the sub-signal component, which has the maximum value of z e r o _ c r o s s i n g s
For our example,   f i r s t _ z e r o _ c r o s s i n g s are found to be as follows:
first   z e r o _ c r o s s i n g s = [ 234   231   223   231   230   235   224   220 ] ;
m a x _ z e r o _ i n d e x = m a x   ( z e r o _ c r o s s i n g s )   =   235 ;
k _ m a x _ c r o s s i n g = a r g m a x k ( z e r o _ c r o s s i n g s ) = 6
Construct the vital signal according to the following expression, the result is shown in Figure 8.
c o n s t r u c t e d _ v i t a l _ s i g n a l ( 1 : max _ z e r o _ i n d e x ) = b e s t _ f i t _ m a t r i x ( k _ m a x _ c r o s s i n g , 1 : m a x _ z e r o _ i n d e x )
For our example; Figure 8 show the construction of first sub-signal component of the vital signal:
c o n s t r u c t e d _ v i t a l _ s i g n a l ( 1 : 235 ) = b e s t _ f i t _ m a t r i x ( 6 , 1 : 235 ) ;
Else if (iteration > 1) Find the next zero-crossing fast time index ( z e r o _ c r o s s i n g s ) of all the “k” best-fit signals
for i = 1 : k
Append the i t h sub-signal component to the previously constructed vital signal.
Find the correlation of the resulting signal
c o r r e l a t i o n ( i ) = a u t o c o r r ( [ constructed _ vital _ signal ;   i th   sub _ signal _ component ] ) ;
end
max ( w i d t h _ c o r r e l a t i o n ( i ) ) ;   i = 1 : k
Append that sub-signal component to the previously constructed vital signal, which has maximum correlation as found by the above expression.For our example, when iteration = 2:
s e c o n d   z e r o _ c r o s s i n g s   =   [ 354       359       351       344       347       360       354       351 ] ;
After we appended each sub-component to the previously constructed vital signal and found the correlation of each signal, it was found to be:
c o r r e l a t i o n   = [ 17.1     13.3     18.3     9.3     16.2     23.8     18.2     22.8 ] ;
k m a x c o r r = a r g m a x k ( c o r r e l a t i o n ) = 6
max _ c o r r _ i n d e x = a r g m a x z e r o _ c r o s s i n g s   ( c o r r e l a t i o n ) = 360
Here the 6th signal sub-component shown the maximum correlation with the previously constructed signal. Now we append the next sub-signal component of the 6th signal to the previously constructed vital signal to find the constructed_vital_signal as follows in Figure 9.
c o n s t r u c t e d _ v i t a l _ s i g n a l ( e n d : max _ c o r r _ i n d e x ) = b e s t _ f i t ( k _ m a x _ c o r r , e n d : max _ c o r r _ i n d e x )
where e n d = l e n g t h ( p r e v i o u s   c o n s t r u c t e d   v i t a l   s i g n a l ) + 1
For our example:
c o n s t r u c t e d _ v i t a l _ s i g n a l   ( 236 : 360 )   =   b e s t _ f i t _ m a t r i x   ( 6 ,   236 : 360 ) ;
(3) Check, if: there are further zero crossings, then go to step 2
Else; stop the loop;
Save; the constructed_vital_signalThe final vital signal obtained is given in the Figure 10.
end
In Algorithm 3, the vital signal is constructed from the sub-signal components of best-fit signals. Now we have to extract the breathing and heart values from the constructed vital signal using the Fast Fourier transform (FFT). The size of FFT used in the evaluation of the frequency domain signal was 215 and the slow time sampling frequency was around 110 samples/s. The frequency resolution thus obtained was 0.2014/min. The respiration rate is extracted from the FFT signal by finding the location of the maximum peak of the spectrum as discussed in detail in the reference [28].

3.2. Dual-Mode Background Subtraction Algorithm for Phone Detection

We proposed a dual-mode background subtraction algorithm to detect the driver’s use of the mobile phone accurately and robustly. Figure 11 shows a diagram for the proposed algorithm. State 1 indicates that the magnitude of the reflected signal exceeding the threshold is not detected in the radar detection area. The received signal is obtained by using the background subtraction method in Section 2.2. If the peak value of the received signal exceeds threshold 1, a transition from state 1 to state 2 occurs. In Equation (14), c k _ s t o p p e d ( t ) is the clutter signal immediately before state 1 → state 2 transition. And k t r a n s i e n t is the slow time index at transition from state 1 to state 2. Parameter k m a r g i n is the slow time margin for transition:
c k _ s t o p p e d ( t ) = c k t r a n s i e n t k m a r g i n 1 ( t ) + ( 1 ) r k t r a n s i e n t k m a r g i n ( t )
State 2 indicates a state in which a signal exceeding a threshold value 1 is detected by an instantly high noise or a moving object such as a hand gesture. In state 2, two received signals are obtained by a dual mode background subtraction method using different clutters. The first is the signal y k ( t ) obtained by using the background substation method in Section 2.2 and the second is the signal y ˘ k ( t ) obtained by using c k _ s t o p p e d ( t ) clutter as shown in the following Equation (15):
y ˘ k ( t ) = r k ( t ) c k _ s t o p p e d ( t )
If it is an instantaneous moving object, the maximum value of y ˘ k ( t ) will be smaller than threshold 1, which means state 2 → state 1 transition. If a new background object appears or an existing background object disappears, the maximum value of y ˘ k ( t ) is still larger than the threshold 1, but the state 2 → state 4 → state 1 transition will occur because the maximum value of y k ( t ) in the T y slow time interval becomes smaller than the threshold 2. If the driver's mobile phone is detected, the transition from state 2 to state 3 occurs because the maximum value of y k ( t ) in the T y slow time interval is greater than threshold 2 and the maximum value of y ˘ k ( t ) is greater than threshold 1 for a certain period of time ( T h o l d ) . State 3 indicates a state in which the use of the mobile phone is detected. If the maximum value of y ˘ k ( t ) becomes smaller than threshold 1, a transition occurs to state 3 → state 1, which means that the use of mobile phone is stopped. Or when the maximum value of y k ( t ) in the T y slow time interval becomes smaller than threshold 2, a transition occurs to state 3 → state 4→ state1, which means that a new background object appears or that the existing background object disappears.
Figure 12a below shows the magnitude of the received signal over fast time when the driver uses the mobile phone. It is shown that the magnitude of the reflected signal is much larger when using the proposed background removal method y ˘ k ( t ) than the conventional background removal method   y k ( t ) . This is because, in the conventional method, the mobile phone is regarded as the background, and the signal reflected from the mobile phone is attenuated.
Figure 12b shows the channel impulse response (CIR). The CIR was estimated using the clean algorithm in [44]. This also shows that the CIR of y ˘ k ( t ) is much larger than the CIR of y k ( t ) for the above reasons. Figure 12c shows the maximum value of the received signal over slow time in the same situation. At about 0.7 s, the mobile phone entered the sensing area, which caused both y k ( t ) and y ˘ k ( t ) to increase in magnitude. At this point, the algorithm transitions from state1 to state2. At about 1 s, mobile phone almost becomes stationary except some minor motions due to texting or touching the screen, and after that point, the magnitude of y ˘ k ( t ) still remains a large value. y k ( t ) represents a relatively small value because the mobile phone is recognized as the background. The mobile phone is continuously used in the detection area, and therefore, the maximum value of y ˘ k ( t ) and the maximum value of y k ( t ) in the T y slow time interval are both larger than the respective threshold values, and thus the state transitions from state 2 to state 3 (mobile phone detection) occurs at about 1.9 s. In addition to the signal characteristics for cell phone use described here, the signal characteristics for various types of driver behavior are specified in the results section. Figure 13a,b below shows the magnitude of the reflected signal over fast time and its CIR when there is a background change. y k ( t ) represents the magnitude of a noise level while the magnitude of y ˘ k ( t ) represent a large value. This characteristics of y k ( t ) can be used to distinguish background changes from using mobile phones with minor motion.

4. Results and Discussion

4.1. Experimental Setup and Reference Data Measurements

In our experiments, we used the commercially available single-chip IR-UWB radar transceiver NVA6201 made by NOVELDA (Novelda AS, Kviteseid, Norway). The radar has a center frequency of 6.8 GHz, a bandwidth of 2.3 GHz, and a transmission output power of −53 dBm/Hz. The pulse repetition frequency is 100 MHz and the slow time sampling frequency (measurement rate) is 110 samples/s. The radar transceiver supports staggered pulse repetition frequency (PRF), a transmission process in which the time between each coherent pulse transmission is patterned and slightly changed. This function extends the maximum unambiguous range (MUR) by making it possible to clearly distinguish the return pulses corresponding to the transmission pulses among many return pulses. In practice, the radar transceiver has a range of almost 10 m and the maximum range is determined by the SNR. The operating temperature of the radar transceiver chip is −40 to 80 degrees, which ensures stable operation under any environment.
The experimental setup of the radar inside the car is shown in Figure 14a. The sensing area is divided into two parts. The region which is near to the radar module is the phone detection region because in this region a driver usually positions his/her phone while reading messages or other mobile phone related activity which can cause distraction of the driver from his primary task of monitoring the road. The region right after the mobile sensing region is the vital sign detection region. It covers the whole front body area of the driver. Figure 14b shows the structure of the IR-UWB radar module. The radar module consists of a radar transceiver, a transmit antenna and a receive antenna. A patch antenna is used, and the antenna and transceiver board are connected by a SubMiniature version A (SMA) connector.
In order to verify the results of our proposed algorithm for the heart rate measurement, we have to compare it some benchmark device such as ECG sensor. We used the ECG sensor module PSL-iECG2 which require only 5 V input voltage (Vcc) and the current consumption is below 50 mA. The amplification is 750 V/V. It uses a selectable 50 or 60 Hz enabled by the notch filter switch. The measurement setup for the connection of the electrodes to the human hands is illustrated in Figure 15a [45].
The sample output of the ECG signal is shown in Figure 15b. We also used some reference measurement for the validity of our respiration rate value by the proposed algorithm. In many studies, respiration rate is measured by wearing a band on the belly or chest [46], but in case of our experiments, if we cover the chest or belly with a band, the measurements are affected because the reflected signal from the belly or chest is distorted by the band. Therefore, we used nasal breath sound recordings from a smartphone as a reference for respiration rate measurements. As in [47], the reference respiration rate can be measured accurately with the errors less than 1 %   for all breathing ranges.

4.2. Respiration & Heart Rate Results

First, we give an example to show the performance of our algorithm as compared to a conventional algorithm such as the one described in reference [28] during the motion period. After the example, we verify the proposed algorithm for different types of motions related to the driving activities and measure the vital signs during the motion period and the results are given in Table 1 and Table 2. In order to show that our algorithm is repeatable, we tested it for five different human subjects and calculated the difference between the reference results and the results from our proposed algorithm and summarized the results in Table 3.
In order to show that the proposed algorithm extract the vital signal even if there is some motion of the body, we made an experiment as follows. In the experiment, slight motion of the upper part of the body is made. The vital signals that are obtained by the conventional algorithm as well as by the proposed algorithm are shown in Figure 16.
As is shown in Figure 16, in a state of slight motion, conventional algorithm looks for the signal at fast time index having maximum variance as the vital signal, however, the motion is not due to the respiration and heartbeat motion but it is due to the random body motion as shown by the values of the R-square in blue-color line in Figure 16b. The vital signal obtained by the conventional algorithm in Figure 16a does not have any trace of periodicity, whereas the respiration signal obtained by the proposed algorithm in this example has good periodic pattern and it captures the vital signal information. In Figure 17, the frequency domain signal is shown. The highest peak i.e., the peak at 21 cycles/minutes represent the breathing frequency whereas the second highest peak in the heart frequency range i.e., at 62 cycles per minute represent the heart rate of the human.
Driving involve different kinds of motion of hand, head and body such as looking at the side mirror to check the road behind, or picking a glass in front of the driver or hand motion associated with the steering operation. We also taken into account, the external effects like turning the car, applying the brakes due to external factor and accelerating the car. The conventional algorithms for vital signs measurement didn’t specifically deal with vital signs measurement during such motions. We measured the vital signs during the motion activities and summarized the results in Table 1 and Table 2 as follows. The values of the measurements are rounded to the nearest integer values.
Table 1 show the results for the estimated respiration rate for different body motion states and it also shows the number of best fit signals reflected from body for each motion case. In case of body motion related to driving such as in the Table 1, our proposed algorithm can find some good-fit signals which can be used to extract the vital signal and hence it can give good measurement results.
We also measured the heart rate during different body movements which may be made by a driver. The results are summarized in Table 2 and compared with the reference heart rate measurements. The values of the measurements are rounded off to the nearest integer values.
The above Table 2 shows the results of the heart rate measurements during the body movements of a driver. The results are related to the general movements of human as well as specific motions related to driving activities. The results in Table 2 show that the algorithm works well under various motion conditions. During huge body movements, which rarely happen during driving activities such as turning the car speedily and passing at high speed over a speed bump, there is no signal which can fit a sinusoidal signal (the R-square values of the signals obtained during huge motion of the car are shown in Figure 18), which clearly indicates that all the R-square values are below the threshold value for the best fit signals i.e., 0.3. Usually, during driving the most of motions are made by hands/arms, shoulder and chest so vital signs measurements by our proposed algorithm is suitable for the driving conditions.
In order to show the validation of the proposed algorithm, we used different human subjects for our experiments. The average error (difference between the reference value and the estimated value) results of the vital signs for five humans are shown below in Table 3. The people involved in these experiments were healthy and aged between 24–38 years old.
The results in Table 3 show that the algorithm works correctly and robustly regardless of the person.

4.3. Mobile Phone Detection Results

In order to verify the proposed cell phone detection algorithm, several experiments were conducted. As shown in the diagram in Figure 11, the signals received in the sensing area are processed in two ways to remove the background and two post-processed signals are obtained. The use of the mobile phone is detected while comparing the maximum values of the processed signals with the threshold value. Threshold 1 is 3.0 and threshold 2 is 0.8. T h o l d is 1.2 s and T y is 1.0 s. These configuration parameters have been optimized so that the alarms are generated and released naturally and the detection error is minimized through several experiments.
The first is an experiment to determine whether the proposed algorithm correctly recognizes that the mobile phone is in use when the driver uses the mobile phone in various ways. Figure 19 shows the characteristics of the received signal over time for each case. Both y k ( t ) and y ˘ k ( t ) are consistently higher than thresholds. In this case, the state transition in the algorithm occurs in order state 1 (no detection) → state 2 (no detection or moving object) → state 3 (mobile phone detection).
Next is whether the algorithm will detect moving objects as a mobile phone or not when there is a moving object in the sensing area for a while. Figure 20 shows the characteristics of the received signal over time for each case. It can be seen that the magnitudes of the signals y k ( t ) and y ˘ k ( t ) instantaneously become larger than the thresholds and the magnitude becomes smaller as the object moves out of the sensing area. In this case, state transitions occur as state 1 (no detection) →state 2 (no detection or moving object) → state 1 (no detection).
Next is, if there is a change of the background object in the sensing area, whether the proposed algorithm detects this change as the use of the mobile phone or not. Figure 21 shows the characteristics of the reflected signal over time for each case.
The magnitude of the signal y ˘ k ( t ) is consistently greater than the threshold 1, but the magnitude of y k ( t ) is only increased when the object enters the sensing area, and soon it become less than the threshold 2. In this case, the state transition occurs as state 1 (no detection) → state 2 (no detection or moving object) → state 4 (change of background object) → state 1 (no detection) or as state 1 (no detection) → state 2 (no detection or moving object) → state 3 (mobile phone detection) → state 4 (change of background object) → state 1 (no detection). Table 4 shows the detection results when 5 people repeated 50 repetitive actions for each experimental case. There is no miss detection at all and it also shows that there is no false alarm when there is no cell phone use. Moreover, the proposed algorithm does not generate false alarms for instantaneous hand gestures or some moving objects and even if a new background object appears or disappears.

5. Conclusions

We have presented techniques based on IR-UWB radar that can be helpful in developing methods for preventing car crashes from happening. The monitoring of vital signs related to drowsiness driving and the detection of mobile phone usage by the driver are the two main purposes of this paper. Firstly, we presented how to extract vital signs signal from the measurements during motion of the body due to driving activity. We found the fast time locations which have better information about the vital signal of the driver. We divided the signals at those locations into segments and constructed the vital signal based on the correlation concept. After construction of the vital signal, an FFT algorithm was applied and the respiration and heart rate were found. Another objective of this paper was to detect the use of mobile phones while driving. The proposed algorithm distinguishes and detects the driver’s cell phone use from various other actions or changes inside the car using the dual mode background subtraction method. Experimental results shown that the proposed mobile phone detection algorithm works perfectly in most of the scenarios that can occur in a car. If these technologies are combined, then it may be very useful for avoiding the car crashes due to drowsiness or mobile phone usage of drivers.

Acknowledgments

This research was supported by the Research Grant from Hanyang University through the Korea Agency for Infrastructure Technology Advancement funded by the Ministry of Land, Infrastructure and Transport of the Korean government (17RTRP-B067918-05).

Author Contributions

Seong Kyu Leem and Faheem Khan designed and carried out the experiments, computer simulations and wrote the paper under the supervision of Sung Ho Cho.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Violence and Injury Prevention. Global Status Report on Road Safety 2013: Supporting a Decade of Action; World Health Organization: Geneva, Switzerland, 2013. [Google Scholar]
  2. Horswill, M.S.; McKenna, F.P. Drivers’ hazard perception ability: Situation awareness on the road. In A Cognitive Approach to Situation Awareness: Theory and Application; Banbury, S., Tremblay, S., Eds.; Ashgate: Aldershot, UK, 2004; pp. 155–175. [Google Scholar]
  3. Klauer, S.G.; Dingus, T.A.; Neale, V.L.; Sudweeks, J.D.; Ramsey, D.J. The Impact of Driver Inattention on Near-Crash/Crash Risk: An Analysis Using the 100-Car Naturalistic Driving Study Data; National Highway Traffic Safety Administration: Washington, DC, USA, 2006. [Google Scholar]
  4. Romin, L.A.; John, A.; Bissonette, J.A. Deer: Vehicle collisions: Status of state monitoring activities and mitigation efforts. Wildl. Soc. Bull. 1996, 24, 276–283. [Google Scholar]
  5. Retting, R.A.; Ferguson, S.A.; McCartt, A.T. A review of evidence-based traffic engineering measures designed to reduce pedestrian–motor vehicle crashes. Am. J. Public Health 2003, 93, 1456–1463. [Google Scholar] [CrossRef] [PubMed]
  6. Wanvik, P.O. Effects of road lighting: An analysis based on Dutch accident statistics 1987–2006. Accid. Anal. Prev. 2009, 41, 123–128. [Google Scholar] [CrossRef] [PubMed]
  7. Pérez, J.; Seco, F.; Milanés, V.; Jiménez, A.; Díaz, J.C.; De Pedro, T. An RFID-based intelligent vehicle speed controller using active traffic signals. Sensors 2010, 10, 5872–5887. [Google Scholar] [CrossRef] [PubMed]
  8. Jones, W.D. Keeping cars from crashing. IEEE Spectr. 2001, 38, 40–45. [Google Scholar] [CrossRef]
  9. Vahidi, A.; Eskandarian, A. Research advances in intelligent collision avoidance and adaptive cruise control. IEEE Trans. Intell. Transp. Syst. 2003, 4, 143–153. [Google Scholar] [CrossRef]
  10. Widmann, G.R.; Daniels, M.K.; Hamilton, L.; Humm, L.; Riley, B.; Schiffmann, J.K.; Schnelker, D.E.; Wishon, W.H. Comparison of Lidar-Based and Radar-Based Adaptive Cruise Control Systems; SAE Technical Paper No. 2000-01-0345; SAE International: Warrendale, PA, USA, 2000. [Google Scholar] [CrossRef]
  11. Ji, Q.; Zhu, Z.; Lan, P. Real-time nonintrusive monitoring and prediction of driver fatigue. IEEE Trans. Veh. Technol. 2004, 53, 1052–1068. [Google Scholar] [CrossRef]
  12. Rodríguez-Ibáñez, N.M.; García-González, A.; Fernández-Chimeno, M.; Ramos-Castro, J. Drowsiness detection by thoracic effort signal analysis in real driving environments. In Proceedings of the 2011 Annual International Conference IEEE Engineering in Medicine and Biology Society (EMBS), Boston, MA, USA, 20 August–3 September 2011. [Google Scholar]
  13. Liu, C.C.; Hosking, S.G.; Lenné, M.G. Predicting driver drowsiness using vehicle measures: Recent insights and future challenges. J. Saf. Res. 2009, 40, 239–245. [Google Scholar] [CrossRef] [PubMed]
  14. Forsman, P.M.; Vila, B.J.; Short, R.A.; Mott, C.G.; Van Dongen, H.P.A. Efficient driver drowsiness detection at moderate levels of drowsiness. Accid. Anal. Prev. 2013, 50, 341–350. [Google Scholar] [CrossRef] [PubMed]
  15. Fan, X.; Yin, B.-C.; Sun, Y.-F. Yawning detection based on Gabor wavelets and LDA. J. Beijing Univ. Technol. 2009, 35, 409–413. [Google Scholar]
  16. Zhang, Z.; Zhang, J. A new real-time eye tracking based on nonlinear unscented Kalman filter for monitoring driver fatigue. J. Control Theory Appl. 2010, 8, 181–188. [Google Scholar] [CrossRef]
  17. Yin, B.-C.; Fan, X.; Sun, Y.-F. Multiscale dynamic features based driver fatigue detection. Int. J. Pattern Recognit. Artif. Intell. 2009, 23, 575–589. [Google Scholar] [CrossRef]
  18. Akin, M.; Kurt, M.B.; Sezgin, N.; Bayram, M. Estimating vigilance level by using EEG and EMG signals. Neural Comput. Appl. 2008, 17, 227–236. [Google Scholar] [CrossRef]
  19. Khushaba, R.N.; Kodagoda, S.; Lal, S.; Dissanayake, G. Driver drowsiness classification using fuzzy wavelet-packet-based feature-extraction algorithm. IEEE Trans. Biomed. Eng. 2011, 58, 121–131. [Google Scholar] [CrossRef] [PubMed]
  20. Liang, W.C.; Yuan, J.; Sun, D.C.; Lin, M.H. Changes in physiological parameters induced by indoor simulated driving: Effect of lower body exercise at mid-term break. Sensors 2009, 9, 6913–6933. [Google Scholar] [CrossRef] [PubMed]
  21. Yang, G.; Lin, Y.; Bhattacharya, P. A driver fatigue recognition model based on information fusion and dynamic Bayesian network. Inf. Sci. 2010, 180, 1942–1954. [Google Scholar] [CrossRef]
  22. Kokonozi, A.K.; Michail, E.M.; Chouvarda, I.C.; Maglaveras, N.M. A study of heart rate and brain system complexity and their interaction in sleep-deprived subjects. In Proceedings of the 35th Annual Computers in Cardiology Conference, Bologna, Italy, 14–17 September 2008; pp. 969–971. [Google Scholar]
  23. Reyner, L.A.; Horne, J.A. Falling asleep whilst driving: Are drivers aware of prior sleepiness? Int. J. Leg. Med. 1998, 111, 120–123. [Google Scholar] [CrossRef]
  24. Mittal, A.; Kumar, K.; Dhamija, S.; Kaur, M. Head movement-based driver drowsiness detection: A review of state-of-art techniques. In Proceedings of the 2016 IEEE International Conference on Engineering and Technology (ICETECH), Coimbatore, India, 17–18 March 2016. [Google Scholar]
  25. Davide, D.; Decarli, N.; Guerra, A.; Al-Rimawi, A.; Puchades, V.M.; Prati, G.; De Angelis, M.; Fraboni, F.; Pietrantoni, L. High-Accuracy Tracking Using Ultrawideband Signals for Enhanced Safety of Cyclists. Mob. Inf. Syst. 2017, 2017, 8149348. [Google Scholar]
  26. Khan, F.; Leem, S.L.; Cho, S.H. Hand-Based Gesture Recognition for Vehicular Applications Using IR-UWB Radar. Sensors 2017, 17, 833. [Google Scholar] [CrossRef] [PubMed]
  27. Win, M.Z.; Dardari, D.; Molisch, A.F.; Wiesbeck, W.; Zhang, J. History and applications of UWB. Proc. IEEE 2009, 97, 198–204. [Google Scholar] [CrossRef]
  28. Khan, F.; Cho, S.H. A Detailed Algorithm for Vital Sign Monitoring of a Stationary/Non-Stationary Human through IR-UWB Radar. Sensors 2017, 17, 290. [Google Scholar] [CrossRef] [PubMed]
  29. Khan, F.; Choi, J.W.; Cho, S.H. Vital sign monitoring of a non-stationary human through IR-UWB radar. In Proceedings of the 4th IEEE International Conference on Network Infrastructure and Digital Content (IC-NIDC), Beijing, China, 19–21 September 2014. [Google Scholar]
  30. Lazaro, A.; Girbau, D.; Villarino, R. Techniques for clutter suppression in the presence of body movements during the detection of respiratory activity through UWB radars. Sensors 2014, 14, 2595–2618. [Google Scholar] [CrossRef] [PubMed]
  31. Hu, X.; Jin, T. Short-Range Vital Signs Sensing Based on EEMD and CWT Using IR-UWB Radar. Sensors 2016, 16, 2025. [Google Scholar] [CrossRef] [PubMed]
  32. Engelberg, J.K.; Hill, L.L.; Rybar, J.; Styer, T. Distracted driving behaviors related to cell phone use among middle-aged adults. J. Transp. Health 2015, 2, 434–440. [Google Scholar] [CrossRef]
  33. Odero, W.; Garner, P.; Zwi, A. Road traffic injuries in developing countries: A comprehensive review of epidemiological studies. Trop. Med. Int. Health 1997, 2, 445–460. [Google Scholar] [CrossRef] [PubMed]
  34. Nasar, J.L.; Troyer, D. Pedestrian injuries due to mobile phone use in public places. Accid. Anal. Prev. 2013, 57, 91–95. [Google Scholar] [CrossRef] [PubMed]
  35. Violanti, J.M.; Marshall, J.R. Cellular phones and traffic accidents: An epidemiological approach. Accid. Anal. Prev. 1996, 28, 265–270. [Google Scholar] [CrossRef]
  36. Yusuf, A.; Bulan, O.; Loce, R.P.; Paul, P. Driver cell phone usage detection from HOV/HOT NIR images. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition Workshops, Columbus, OH, USA, 23–28 June 2014. [Google Scholar]
  37. Smith, M.J.; Stephens, D.R. Detecting Use of a Mobile Device by a Driver of a Vehicle, Such as an Automobile. U.S. Patent Application No. 13/290,126, 23 August 2012. [Google Scholar]
  38. Xu, B.; Loce, R.P. A machine learning approach for detecting cell phone usage. In Proceedings of the IS&T/SPIE Electronic Imaging, International Society for Optics and Photonics, San Francisco, CA, USA, 4 March 2015. [Google Scholar]
  39. Yim, D.H.; Cho, S.H. An Equidistance Multi-Human Detection Algorithm Based on Noise Level Using Mono-static IR-UWB Radar System. In Proceedings of the 2014 International Conference on Future Communication, Information and Computer Science (FCICS 2014), Beijing, China, 22–23 May 2014. [Google Scholar]
  40. Tomas, A. Parameter Estimation and Waveform Fitting for Narrowband Signals. Doctoral Thesis, KTH Royal Institute of Technology, Stockholm, Sweden, 2005. [Google Scholar]
  41. Rife, D.; Boorstyn, R. Single tone parameter estimation from discrete-time observations. IEEE Trans. Inf. Theory 1974, 20, 591–598. [Google Scholar] [CrossRef]
  42. Petre, S.; Moses, R.L.; Friedlander, B.; Soderstrom, T. Maximum likelihood estimation of the parameters of multiple sinusoids from noisy measurements. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 378–392. [Google Scholar]
  43. Steel, R.G.D.; Torrie, J.H. Principles and Procedures of Statistics with Special Reference to the Biological Sciences; Mcgraw-Hill Book Company: New York, NY, USA, 1960. [Google Scholar]
  44. Yano, S.M. Investigating the ultra-wideband indoor wireless channel. In Proceedings of the IEEE 55th Vehicular Technology Conference, VTC Spring 2002, Birmingham, AL, USA, 6–9 May 2002. [Google Scholar]
  45. PSL-iECG2 (Mini-Size 2ch ECG Sensor Module with Isolation). Available online: http://physiolab.en.ec21.com/PSL-iECG2_Mini-size_2ch_ECG_sensor--9063127_9816937.html (accessed on 29 May 2017).
  46. AL-Khalidi, F.Q.; Saatchi, R.; Burke, D.; Elphick, H.; Tan, S. Respiration rate monitoring methods: A review. Pediatr. Pulmonol. 2011, 46, 523–529. [Google Scholar] [CrossRef] [PubMed]
  47. Nam, Y.; Reyes, B.A.; Chon, K.H. Estimation of Respiratory Rates Using the Built-in Microphone of a Smartphone or Headset. IEEE J. Biomed. Health Inform. 2016, 26, 1493–1501. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Noisy vital signal obtained by conventional algorithm during motion period.
Figure 1. Noisy vital signal obtained by conventional algorithm during motion period.
Sensors 17 01240 g001
Figure 2. “Vital signal” vs slow time and fast time.
Figure 2. “Vital signal” vs slow time and fast time.
Sensors 17 01240 g002
Figure 3. Process block diagram for extraction of driver vital signs.
Figure 3. Process block diagram for extraction of driver vital signs.
Sensors 17 01240 g003
Figure 4. Data fitting algorithm results (a) fitting frequency higher than the respiration frequencies range (b) fitting frequency lower than the respiration frequencies range (c) fitting magnitude less than vital signal threshold (d) fitting signal for signal having information of vital signs.
Figure 4. Data fitting algorithm results (a) fitting frequency higher than the respiration frequencies range (b) fitting frequency lower than the respiration frequencies range (c) fitting magnitude less than vital signal threshold (d) fitting signal for signal having information of vital signs.
Sensors 17 01240 g004aSensors 17 01240 g004b
Figure 5. R square values while (a) stationary (b) moving hand near abdominal region (c) moving upper parts of body slightly.
Figure 5. R square values while (a) stationary (b) moving hand near abdominal region (c) moving upper parts of body slightly.
Sensors 17 01240 g005aSensors 17 01240 g005b
Figure 6. A best fit signal which has the corresponding r-square value of 0.46.
Figure 6. A best fit signal which has the corresponding r-square value of 0.46.
Sensors 17 01240 g006
Figure 7. The vital signal construction from the candidate sub-signal components of the best fit signals.
Figure 7. The vital signal construction from the candidate sub-signal components of the best fit signals.
Sensors 17 01240 g007
Figure 8. The first sub-signal component of vital_constructed_signal.
Figure 8. The first sub-signal component of vital_constructed_signal.
Sensors 17 01240 g008
Figure 9. The constructed vital signal after second iteration.
Figure 9. The constructed vital signal after second iteration.
Sensors 17 01240 g009
Figure 10. The final constructed vital signal.
Figure 10. The final constructed vital signal.
Sensors 17 01240 g010
Figure 11. State diagram of dual-mode background subtraction algorithm for phone usage detection.
Figure 11. State diagram of dual-mode background subtraction algorithm for phone usage detection.
Sensors 17 01240 g011
Figure 12. Signal characteristics when using a mobile phone (a) received signal over fast time (b) channel impulse response (c) the maximum value of the received signal over slow time.
Figure 12. Signal characteristics when using a mobile phone (a) received signal over fast time (b) channel impulse response (c) the maximum value of the received signal over slow time.
Sensors 17 01240 g012
Figure 13. Signal characteristics when there is a background change (a) received signal over fast time (b) channel impulse response.
Figure 13. Signal characteristics when there is a background change (a) received signal over fast time (b) channel impulse response.
Sensors 17 01240 g013
Figure 14. (a) Experimental setup inside the car (b) IR-UWB radar module.
Figure 14. (a) Experimental setup inside the car (b) IR-UWB radar module.
Sensors 17 01240 g014
Figure 15. (a) Measurement setup for the ECG signal acquisition; (b) The signal obtained from the ECG sensor module PSL-iECG2 (measurement reference device).
Figure 15. (a) Measurement setup for the ECG signal acquisition; (b) The signal obtained from the ECG sensor module PSL-iECG2 (measurement reference device).
Sensors 17 01240 g015
Figure 16. (a) Comparison of the vital signal obtained from the conventional vs proposed algorithm while the human is slightly moving (b) R-square and normalized variance values at each fast-time index during motion state.
Figure 16. (a) Comparison of the vital signal obtained from the conventional vs proposed algorithm while the human is slightly moving (b) R-square and normalized variance values at each fast-time index during motion state.
Sensors 17 01240 g016aSensors 17 01240 g016b
Figure 17. Spectrum obtained by applying FFT algorithm to the signal obtained by proposed algorithm in Figure 16a.
Figure 17. Spectrum obtained by applying FFT algorithm to the signal obtained by proposed algorithm in Figure 16a.
Sensors 17 01240 g017
Figure 18. The R-square values in case of huge body movements (all the values are below the threshold value of 0.3).
Figure 18. The R-square values in case of huge body movements (all the values are below the threshold value of 0.3).
Sensors 17 01240 g018
Figure 19. The maximum value of the received signal over slow time when using a mobile phone (a) texting (b) scrolling, touching (c) viewing.
Figure 19. The maximum value of the received signal over slow time when using a mobile phone (a) texting (b) scrolling, touching (c) viewing.
Sensors 17 01240 g019
Figure 20. The maximum value of the received signal over slow time with momentary movements. (a) moving hand (b) moving water bottle (c) moving cell phone.
Figure 20. The maximum value of the received signal over slow time with momentary movements. (a) moving hand (b) moving water bottle (c) moving cell phone.
Sensors 17 01240 g020aSensors 17 01240 g020b
Figure 21. The maximum value of the received over slow time when there is a background change (a) water bottle (appearing) (b) water bottle (disappearing)
Figure 21. The maximum value of the received over slow time when there is a background change (a) water bottle (appearing) (b) water bottle (disappearing)
Sensors 17 01240 g021
Table 1. Experimental results for breathing rate for different motion conditions.
Table 1. Experimental results for breathing rate for different motion conditions.
Movement CaseNo. of Best-Fit Columns above R-Square Threshold (0.3)Reference Measurement (Breathing Rate)Estimated Value (Breathing Rate)
stationary191818
general movementshand motion near body81819
shoulder motion for steering41816
gesture made by hand away from body111818
lips motion during speaking141818
specific driving movementshead movement to watch the mirror of the car71818
turning the car41820
applying the brakes71817
accelerating the car101818
Table 2. Experimental results for heart rate measurements for different motion conditions.
Table 2. Experimental results for heart rate measurements for different motion conditions.
Movement CaseNo. of Best-Fit Columns above R-Square Threshold (0.3)Reference Measurement (Heart Rate)Estimated Value (Heart Rate)
stationary177272
general movementshand motion near body97271
gesture made by hand away from body117372
lips motion during speaking127273
specific driving movementshead movement to watch the mirror of the car87274
turning the car57275
applying the brakes77471
accelerating the car127272
Table 3. Average error for respiration and heart rate of different human subjects (five humans).
Table 3. Average error for respiration and heart rate of different human subjects (five humans).
Movement CaseAverage Error (Bpm) (Respiration Rate)Average Error (Bpm) (Heart Rate)
stationaryalmost zeroalmost zero
general movementshand motion near body0.71.1
gesture made by hand away from body0.30.9
lips motion during speaking0.20.6
specific driving movementslead movement to watch the mirror of the car0.32.0
turning the car1.22.5
applying the brakes0.51.7
accelerating the car0.41.3
Table 4. Mobile phone detection results.
Table 4. Mobile phone detection results.
ObjectDetail ActionDetection RateFalse Alarm Rate
mobile phonetexting100% (1)0% (2)
scrolling, touching100% (1)
viewing100% (1)
moving object for a momentmoving handN/A0% (3)
moving water bottle0% (3)
moving cell phone0% (3)
change of backgroundwater bottle (appearing)0% (4)
water bottle (disappearing)0% (4)
(1) The probability of being recognized by a mobile phone when using a real mobile phone. (2) The probability of recognizing the use of a mobile phone even though there is nothing in the detection area. (3) The probability of detecting an instantly moving object as the use of mobile phone. (4) The probability of detecting change of background as the use of mobile phone.

Share and Cite

MDPI and ACS Style

Leem, S.K.; Khan, F.; Cho, S.H. Vital Sign Monitoring and Mobile Phone Usage Detection Using IR-UWB Radar for Intended Use in Car Crash Prevention. Sensors 2017, 17, 1240. https://doi.org/10.3390/s17061240

AMA Style

Leem SK, Khan F, Cho SH. Vital Sign Monitoring and Mobile Phone Usage Detection Using IR-UWB Radar for Intended Use in Car Crash Prevention. Sensors. 2017; 17(6):1240. https://doi.org/10.3390/s17061240

Chicago/Turabian Style

Leem, Seong Kyu, Faheem Khan, and Sung Ho Cho. 2017. "Vital Sign Monitoring and Mobile Phone Usage Detection Using IR-UWB Radar for Intended Use in Car Crash Prevention" Sensors 17, no. 6: 1240. https://doi.org/10.3390/s17061240

APA Style

Leem, S. K., Khan, F., & Cho, S. H. (2017). Vital Sign Monitoring and Mobile Phone Usage Detection Using IR-UWB Radar for Intended Use in Car Crash Prevention. Sensors, 17(6), 1240. https://doi.org/10.3390/s17061240

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop