The present invention relates broadly to a device and method for sleep monitoring, in particular to a device and method for determining time-to-sleep and wake periods during sleep and to a device and method for determining rapid eye movement (REM) sleep and non REM (NREM) sleep.
Having a good sleep at night is a key to perform the best during the day and to keep the fitness and well-being.
Many research studies show that there is a major link between sleep problems and a variety of serious health conditions including depression, heart disease, obesity and shorter life expectancy. Only loss of one hour of sleep in several nights can cause significantly negative effect on performance, learning skill, mood and safety. Long night sleepers who sleep more than 9 hrs or more also show risk of coronary heart disease and risk of stroke.
A personal device for sleep monitoring is thus desirable.
To track sleep condition, important parameters as referred to herein are Time-to-sleep, Total time-in-bed, Total sleep time, Time awake during sleep, Sleep efficiency and sleep quality (architecture/stages).
Time-to-sleep is also called Sleep Latency or Sleep Onset. Normal people without significant sleep deprivation typically take more than 20 minutes to fall asleep. Time-to-sleep is also correlated to sleep deprivation by referring to the MSLT (Multiple Sleep Latency Test) table shown in Table 1 below. MSLT provides a sleepiness of a subject and the severity of their sleep debt from the time taken to fall asleep.
Total time-in-bed is the recorded time that the user has spent in-bed in total when they have entered and exited the sleep monitoring mode.
Total sleep time is the total sleep time recorded, which is the difference between the total time-in-bed and the time awake during sleep.
Time awake during sleep are the periods of wakefulness/restlessness identified during sleep and the record of the number of times awake and their duration.
Sleep efficiency is determined by the ratio of total sleep time over total time-in-bed.
Sleep quality can be determined by one or more of total sleep time, amount of REM, NREM sleep and Sleep stages, amount of movement and wakefulness and sleep diary, i.e., record of daily sleep hours & the feeling of the next day to know how much sleep necessary for individual.
REM is also sometimes referred to as “dream” sleep. NREM includes 3 stages called N1, N2 and N3.
Many of a user's physiological functions such as brain wave activity, breathing, and heart rate are quite variable during REM sleep, but are extremely regular in NREM sleep.
It has been found that during REM sleep, the brain is restored and captures memories which allow learning to take place, etc. Heart rate, blood pressure and body temperature will typically increase. Generally, 20-25% of total sleep time is REM sleep. N1 is a transition between wakefulness and sleep. N2 is during light sleep, in which the heart rate is slower. Generally, 50-55% of total sleep time is N2 sleep. N3 is during deep sleep to restore the physical body, during which the body temperature and blood pressure will typically decrease.
Sleep Cycle consists of consecutive REM and NREM sleep stages. An average duration for each cycle is about 90 to 110 minutes and there are about 4 to 6 cycles for normal sleeping hours over the night (compare
There are several devices in the market to monitor sleep efficiency or quality. Polysomnography (PSG) is the current gold standard for sleep study to diagnose sleep disorder. PSG includes the monitoring of many different physiological signals such as heart rate variability (HRV), respiration, electroencephalography (EEG), eletromyography (EMG), electrooculagram (EOG); and it needs to be performed in a sleep laboratory under the supervision of sleep experts. Although PSG is an important tool for sleep diagnosis, it is an uncomfortable and costly procedure, especially when multiple nights of observation are required. Some wearable devices have also been developed to ease these inconveniences. However, those devices are not generally able to measure sleep quality nor sleep efficiency accurately.
Embodiments of the present invention provide at least an alternative device and method for sleep monitoring.
According to a first aspect of the present invention there is provided a method for determining time-to-sleep and wake periods during sleep, the method comprising the steps of obtaining motion data representative of motion of a user; detecting the time-to-sleep from the motion data based on a first time-above-threshold (TAT) threshold and a first proportional integration method (PIM) threshold; and detecting the wake periods during sleep from the motion data based on a second TAT threshold and a second PIM threshold.
According to a second aspect of the present invention there is provided a device for determining time-to-sleep and wake periods during sleep, the device comprising a sensor for obtaining motion data representative of motion of a user; and a processor for detecting the time-to-sleep from the motion data based on a first time-above-threshold (TAT) threshold and a first proportional integration device (PIM) threshold; and for detecting the wake periods during sleep from the motion data based on a second TAT threshold and a second PIM threshold.
According to a third aspect of the present invention there is provided a method for determining rapid eye movement (REM) sleep and non REM (NREM) sleep, the method comprising the steps of obtaining physiological signal data of a user; splitting the physiological signal data into respective data subsets; and detecting REM sleep and non REM (NREM) sleep in each data subset based on one or more heart rate variability (HRV) features extracted from each data subset based on adaptive thresholds for each HRV feature.
According to a fourth aspect of the present invention there is provided a device for determining rapid eye movement (REM) sleep and non REM (NREM) sleep, the device comprising a sensor for obtaining physiological signal data of a user; and a processor for splitting the physiological signal data into respective data subsets; and for detecting REM sleep and non REM (NREM) sleep in each data subset based on one or more heart rate variability (HRV) features extracted from each data subset based on adaptive thresholds for each HRV feature.
Embodiments of the invention will be better understood and readily apparent to one of ordinary skill in the art from the following written description, by way of example only, and in conjunction with the drawings, in which:
Embodiments of the present invention provide a device and method for sleep monitoring, in particular for determining sleep condition, especially sleep stages (REM, NREM), and/or sleep and wake states.
In the example embodiments described, the sleep stages are determined based on heart rate variability (HRV) and through adaptive thresholds derived from the average of a data-subset to determine sleep cycles. Sleep and wake states are identified based on acceleration magnitude and a combination of TAT (Time-above-threshold) and PIM (Proportional Integration Method) thresholds.
Advantageously, embodiments of the present invention can measure sleep stages accurately and effectively with power consumption efficiency, thus reducing battery spent for the wearable device.
Also, the example embodiments described advantageously provide for accurate detection of sleep-onset latency (time taken to fall asleep) through usage of different threshold levels to further differentiate motions during and before sleep. Stringent TAT and PIM thresholds in each level are applied to differentiate motions related to wakefulness and sleep.
In one embodiment, three stages of wake during sleep, REM sleep, and NREM sleep are calculated simultaneously by using motion data, e.g. an acceleration signal measured by an accelerometer (ACC) sensor or a gyroscope, and physiological signal data, for example a photoplethysmography (PPG) signal measured by a PPG sensor.
Example embodiments use stringent TAT and PIM thresholds obtained from experimental data to differentiate motions related to wakefulness and sleep when both conditions are satisfied. Accurate detection of sleep-onset latency is preferably enabled through the combined usage of high/low threshold levels that are set to further differentiate motions during sleep and when trying to fall asleep. Advantageously, high sensitivity thresholds are sensitive to motions when trying to fall asleep, and low sensitivity thresholds are sensitive to motions during sleeping.
The present specification also discloses an apparatus, which may be internal and/or external to the wearable device in example embodiments, for performing the operations of the methods. Such apparatus may be specially constructed for the required purposes, or may comprise a general purpose computer or other device selectively activated or reconfigured by a computer program stored in the computer. The algorithms and displays presented herein are not inherently related to any particular computer or other apparatus. Various general purpose machines may be used with programs in accordance with the teachings herein. Alternatively, the construction of more specialized apparatus to perform the required method steps may be appropriate. The structure of a conventional general purpose computer will appear from the description below. In addition, the present specification also implicitly discloses a computer program, in that it would be apparent to the person skilled in the art that the individual steps of the method described herein may be put into effect by computer code. The computer program is not intended to be limited to any particular programming language and implementation thereof. It will be appreciated that a variety of programming languages and coding thereof may be used to implement the teachings of the disclosure contained herein. Moreover, the computer program is not intended to be limited to any particular control flow. There are many other variants of the computer program, which can use different control flows without departing from the spirit or scope of the invention.
Furthermore, one or more of the steps of the computer program may be performed in parallel rather than sequentially. Such a computer program may be stored on any computer readable medium. The computer readable medium may include storage devices such as magnetic or optical disks, memory chips, or other storage devices suitable for interfacing with a general purpose computer. The computer readable medium may also include a hard-wired medium such as exemplified in the Internet system, or wireless medium such as exemplified in the GSM mobile telephone system. The computer program when loaded and executed on such a general-purpose computer effectively results in an apparatus that implements the steps of the preferred method.
The invention may also be implemented as hardware modules. More particularly, in the hardware sense, a module is a functional hardware unit designed for use with other components or modules. For example, a module may be implemented using discrete electronic components, or it can form a portion of an entire electronic circuit such as an Application Specific Integrated Circuit (ASIC). Numerous other possibilities exist. Those skilled in the art will appreciate that the system can also be implemented as a combination of hardware and software modules.
The described embodiments of the invention described herein relate to a wearable device and a method for sleep monitoring based on motion signals acquired from the user with a motion sensor such as an ACC and/or a gyroscope and based on physiological signals acquired from the user with sensor such as PPG sensor.
In one embodiment, the device can be worn on any location of the user with sufficient skin area to allow the light emitting diode-photo detector (LED-PD) arrangement to acquire the PPG signal and allows the tri-axial ACC to acquire motion signals.
The device 200 according to an example embodiment shown in
With reference to
The total sleep data is split into subsets corresponding to each estimated sleep cycle duration (step 406), and thresholds are set (step 408). For example, a sleep cycle is estimated as 1 hr. Thresholds are set based on the average during each subset/estimated sleep cycle, in the example embodiment.
REM sleep is determined when the HRV features are greater than the thresholds, otherwise determined as NREM sleep (step 410-412). If REM sleep is determined for data falling within the initial period of the whole night data (step 414), for example within the first 45 minutes, the determination is changed to NREM (step 412), otherwise the REM determination is maintained (step 416). The combination of the REM sleep and NREM sleep determinations is used to generate a first or intermediate result for the sleep stages (step 418), where S(i) indicates the sleep stage result in each 3 min measurement interval in the example embodiment. For example, S(i)=3 for REM stage result, and S(i)=2 for NREM stage result.
Smoothing by checking nearest neighbor of sleep stage results for a 3 min measurement interval (step 420) is performed, to remove false condition prior to output of the final sleep stages result (step 422). Details of the nearest neighbor checking method in the example embodiment are shown in the flowchart 500 in
The Sleep cycle is estimated as 1 hr in the example embodiment, and experimental results show a close relation compared to PSG reference. By estimating the sleep cycle as 1 hr, we made the calculation process simple and effective.
As illustrated in
By using both TAT and PIM, the results the example embodiment advantageously reflect substantially all the important factors of movement, including duration, frequency, acceleration and amplitude.
On the other hand, the inventors have unexpectedly found that the ZCM (Zero Crossing Mode) parameter, which is often used in existing techniques, does not fully describe the motion and provides less information related to jerk or toss movement. This is illustrated in Table 2 below.
Large movements (i.e. toss) during sleep are assumed very uncommon during light and deep sleep due to the body slowing down for physical restoration. It is however possible to have sudden muscle jerks, but these are unrelated to wakefulness.
In the example embodiment, lower sensitivity level thresholds for detecting wake-during-sleep are set to 90% of the “Slow Toss (1×)” during sleep values for the TAT and PIM scores, noting again that the example embodiment does deliberately not use ZCM scores for the reasons explained above and illustrated in Table 2.
For very small movements (i.e. jerk), the values of TAT and PIM are very low. In the example embodiment, the thresholds are set at a higher sensitivity level based on the values for “Jerk (1×)” to identify small movements. Because small movements are unlikely related to movements made when awake, these higher sensitivity level thresholds are used in the example embodiment to identify time-to-sleep.
As mentioned above, for larger movements (i.e. toss), the values of TAT and PIM are much higher. The thresholds are set at this lower sensitivity level to identify larger movements that are better correlated to restlessness/wakefulness during sleep in the example embodiments to identify wakefulness/restlessness during sleep, also referred to herein as wake period during sleep, or wake-during-sleep.
In one embodiment, the higher sensitivity threshold levels for TAT & PIM are set to be 1 and 10 respectively, and the lower sensitivity threshold levels for TAT & PIM are set to be 100 and 62 respectively. It is noted again that to identify wake-during-sleep status and Time-to-sleep, both criteria derived from TAT and PIT score need to be satisfied in the example embodiment, to advantageously make the result more accurate. No ZCM scores are used in this example embodiment.
Returning to
Six sleep parameters can be calculated in an example embodiment. The six parameters are Time-to-sleep, Number of awakenings, Total time awake during actual sleep period, Total sleep time, Total time-in-bed, and Sleep efficiency.
Time-to-sleep (sleep-onset latency) is identified based on high sensitivity thresholds (steps 1010 and 1012). If both TAT and PIM scores are lower than the high sensitivity thresholds, the 60 seconds window is classified as quiet period, and quiet period must satisfy consecutive ‘N’ windows, i.e. N windows of little or no movements. N can be about 5-20, preferably about 8-15 in example embodiments. Wake periods during sleep are identified when TAT and PIM scores exceed pre-determined low sensitivity thresholds (steps 1010 and 1014). If both TAT and PIM are higher than the low sensitivity thresholds, the 60 seconds window is classified as wake period during sleep.
Sleep Efficiency is determined by calculating Total sleep time/Total time in bed. Actionable feedback can be provided for sleep efficiency, MSLT score, Sleep debt and optimal alarm function. If sleep efficiency is greater than about 85%, it can be considered normal according to current understanding. The MSLT score can be used to show how serious the user's sleep deprivation is. Sleep debt shows whether the user is getting enough sleep hours. Optimal alarm function can be set and vibration can be used in an example embodiment.
One non-limiting example of a preferred accelerometer that can be adapted for use in the device is a triple-axis accelerometer MMA8652FC available from Freescale Semiconductor, Inc. This accelerometer can provide the advantage of measuring acceleration in all three directions with a single package. Alternatively, several single-axis accelerometers oriented to provide three-axis sensing can be used in different embodiments.
The device 1601 also includes a second sensing module 1603, such as an LED-PD module, for obtaining a physiological signal of the user. The device 1601 also includes a data processing and computational module 1604, such as a processor, which is arranged to receive and process the acceleration information from the signal sensing module 1602 and the physiological signal from the measurement module 1603. The device 1601 also includes a display unit 1606 for displaying a result to a user of the device 1601 and for receiving user input via touch screen technology. The device 1601 in this embodiment further includes a wireless transmission module 1608 arranged to communicate wirelessly with a telecommunications device 1610 of the assembly 1600. The telecommunication device 1610 includes a wireless receiver module 1612 for receiving signals from the wearable device 1601, a display unit 1614 for displaying a result to a user of the telecommunication device 1610 and for receiving user input via touch screen technology.
The first and second TAT thresholds may be different. The first TAT threshold may be lower than the second TAT threshold.
The first and second PIM thresholds may be different. The first PIM threshold may be lower than the second PIM threshold.
Detecting the time-to-sleep from the motion data may comprise dividing the motion data into time windows; determining TAT and PIM scores for each time window, and identifying windows in which the TAT and PIM scores are below the first TAT threshold and the first PIM threshold.
Detecting the wake periods during sleep from the motion data may comprise dividing the motion data into time windows; determining TAT and PIM scores for each time window, and identifying windows in which the TAT and PIM scores exceed the second TAT threshold and the second PIM threshold.
The motion data may comprise multi-axis motion signals. The method may further comprise calculating a resultant magnitude of the multi-axis motion signals using bandpass filtering and root-mean-square (RMS) calculation.
The first and second TAT thresholds may be respective number of times a magnitude derived from the motion data is above an acceleration threshold. The acceleration threshold may be in a range from 0.1 to 0.2 G, and preferably about 0.15 G. The first and second PIM thresholds may be respective areas under a magnitude curve derived from the motion data. The respective areas may be estimated using a trapezoid rule.
The determining of the time-to-sleep and wake periods during sleep may not be based on zero-crossing-mode detection based on the motion data.
The first and second TAT thresholds may be different. The first TAT threshold may be lower than the second TAT threshold.
The first and second PIM thresholds may be different. The first PIM threshold may be lower than the second PIM threshold.
Detecting the time-to-sleep from the motion data may comprise dividing the motion data into time windows; determining TAT and PIM scores for each time window, and identifying windows in which the TAT and PIM scores are below the first TAT threshold and the first PIM threshold.
Detecting the wake periods during sleep from the motion data may comprise dividing the motion data into time windows; determining TAT and PIM scores for each time window, and identifying windows in which the TAT and PIM scores exceed the second TAT threshold and the second PIM threshold.
The motion data may comprise multi-axis motion signals. The processor may further be configured for calculating a resultant magnitude of the multi-axis motion signals using bandpass filtering and root-mean-square (RMS) calculation.
The first and second TAT thresholds may be respective number of times a magnitude derived from the motion data is above an acceleration threshold. The acceleration threshold may be in a range from 0.1 to 0.2 G, and preferably about 0.15 G.
The first and second PIM thresholds may be respective areas under a magnitude curve derived from the motion data. The respective areas may be estimated using a trapezoid rule.
The determining of the time-to-sleep and wake periods during sleep may not be based on zero-crossing-mode detection based on the motion data.
Detecting REM sleep and NREM sleep may comprise detecting REM sleep in respective time windows within the data subset. The time windows may correspond to on-stages of a detector for the physiological signal data, the detector operating in an on/off operation mode. The on-stage may be about 3 minutes and the detector may operate in an about 50% on/off operation mode.
In each time window REM sleep and NREM may be detected based on the adaptive thresholds.
The method may further comprise changing a detected REM sleep to a detected NREM sleep if the detected REM sleep is within an initial time period of the obtained physiological signal data. The initial time period may be about 45 minutes.
The method may further comprise comparing a REM sleep and NREM sleep detection result for one time window with respective results for its nearest neighboring time windows. The method may comprise maintaining the detection result in said one window if said detection result is similar to the respective results for the nearest neighboring time windows, and changing said detection result otherwise.
The HRV features may comprise a mean heart rate (meanHR) and a low frequency/high frequency (LF/HF) ratio derived from the physiological signal data.
A first adaptive threshold may be an average of a first HRV feature in each data subset. A second adaptive threshold may be an average of a second HRV feature in each data subset. REM sleep may be detected if the first HRV feature is above a first adaptive threshold and the second HRV feature is above a second adaptive threshold, and NREM sleep may be detected otherwise.
Detecting REM sleep and NREM sleep may comprise detecting REM sleep in respective time windows within the data subset. The time windows may correspond to on-stages of a detector for the physiological signal data, the detector operating in an on/off operation mode. The on-stage may be about 3 minutes and the detector may operate in an about 50% on/off operation mode.
In each time window REM sleep and NREM may be detected based on the adaptive thresholds.
The processor 2104 may further be configured for changing a detected REM sleep to a detected NREM sleep if the detected REM sleep is within an initial time period of the obtained physiological signal data. The initial time period may be about 45 minutes. The processor 2104 may further be configured for comparing a REM sleep and NREM sleep detection result for one time window with respective results for its nearest neighboring time windows. The processor 2104 may be configured for maintaining the detection result in said one window if said detection result is similar to the respective results for the nearest neighboring time windows, and for changing said detection result otherwise.
The HRV features may comprise a mean heart rate (meanHR) and a low frequency/high frequency (LF/HF) ratio derived from the physiological signal data.
A first adaptive threshold may be an average of a first HRV feature in each data subset. A second adaptive threshold may be an average of a second HRV feature in each data subset. REM sleep may be detected if the first HRV feature is above a first adaptive threshold and the second HRV feature is above a second adaptive threshold, and NREM sleep may be detected otherwise.
It will be appreciated by a person skilled in the art that numerous variations and/or modifications may be made to the present invention as shown in the specific embodiments without departing from the spirit or scope of the invention as broadly described. The present embodiments are, therefore, to be considered in all respects to be illustrative and not restrictive. Also, the invention includes any combination of features, in particular any combination of features in the patent claims, even if the feature or combination of features is not explicitly specified in the patent claims or the present embodiments.
For example, while a wrist-worn device is described in some embodiments, the device may be worn on any part of the arms, hip, waist or foot of the user.
Also, according to the human sleep behavior as currently understood, reduction in heart rate and blood pressure occur during NREM sleep. In REM sleep, there is more variation in cardiovascular activity which can cause overall increases in blood pressure and heart rate. The example embodiments described employ mean HR and LF/HF as HRV features. However, it will be appreciated that different HRV features (e.g Standard Deviation of Heart Rate for period of interest (SDHR), Percent of NN intervals>50 ms different from previous (NN) for period of interest (PNN50), Root mean square of successive differences of NN interval for period of interest (RMSSD) and blood flow features (e.g mean Pulse Pressure for the period of interest (mean PP), Average standard deviation of pulse pressure for the period of interest (ASDPP)) can additionally or alternatively be applied in different embodiments, to improve performance.
This application is a continuation of U.S. application Ser. No. 15/540,670, filed Jun. 29, 2017, which is the U.S. National Stage of International Application No. PCT/SG2014/000624, filed on Dec. 30, 2014, published in English. The entire teachings of the above application(s) are incorporated herein by reference
Number | Date | Country | |
---|---|---|---|
Parent | 15540670 | Jun 2017 | US |
Child | 17451938 | US |