On-demand heart rate estimation based on optical measurements

• US 10,285,651 B2
• Filed: 06/30/2016
• Issued: 05/14/2019
• Est. Priority Date: 06/30/2016
• Status: Active Grant

• Current Assignee: Analog Devices, Inc.
• Original Assignee: Analog Devices, Inc.
• Inventors: Demirtas, Sefa, King, Jason D., Adams, Robert, Akl, Tony Joseph, Bernstein, Jeffrey G.
• Primary Examiner(s): Kahelin, Michael W
• Application Number:
  US15/198,438
• Publication Number:
  US 20180000424A1
• Time in Patent Office: 1,048 Days
• CPC Class Codes: A61B 5/02416 using photoplethysmograph signals, e.g. generated by infrared radiation A61B5/14552 takes precedence
• CPC Class Codes: A61B 5/1118 Determining activity level
• CPC Class Codes: A61B 5/1123 Discriminating type of movement, e.g. walking or running A61B5/1116, A61B5/112 take precedence
• CPC Class Codes: A61B 5/6802 Sensor mounted on worn items
• CPC Class Codes: A61B 5/681 Wristwatch-type devices
• CPC Class Codes: A61B 5/7207 of noise induced by motion artifacts
• CPC Class Codes: A61B 5/7221 Determining signal validity, reliability or quality A61B5/7207 and A61B5/7217 take precedence
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: G16H 50/30 for calculating health indices; for individual health risk assessment
• CPC Class Codes: G16H 50/70 for mining of medical data, e.g. analysing previous cases of other patients

Abstract

Activity monitors and smart watches utilizing optical measurements are becoming widely popular, and users expect to get an increasingly accurate estimate of their heart rate (HR) from these devices. These devices are equipped with a light source and an optical sensor which enable estimation of HR using a technique called photoplethysmography (PPG). One of the main challenges of HR estimation using PPG is the coupling of motion into the optical PPG signal when the user is moving randomly or exercising. The present disclosure describes a computationally feasible and fast HR estimation algorithm to be executed at instances of little or no motion. Resulting HR readings may be useful on their own, or be provided to systems that monitor HR continuously to prevent the problem of such systems being locked on an incorrect HR for long periods of time. Implementing techniques described herein leads to more accurate HR measurements.

Claims

1. A system for determining a heart rate from a heartbeat signal present in a first signal generated by a heartbeat sensor, the system comprising:

means for computing an autocorrelation of data samples of the first signal over a first time period;

means for processing a second signal to evaluate motion of the heartbeat sensor during the first time period, the second signal indicative of a motion of the heartbeat sensor with respect to one or more directions;

means for determining whether the data samples of the first signal over the first time period pass a set of confidence checks, the set of confidence checks comprising at least

determining whether the motion of the heartbeat sensor during the first time period satisfies a motion-related criteria;

means for computing a first estimate of the heart rate from the autocorrelation of the data samples of the first signal over the first time period upon determination that the data samples of the first signal over the first time period pass the set of confidence checks;

means for computing a Discrete Fourier Transform (DFT) for data samples of the first signal over a second time period;

means for computing a second estimate of the heart rate from the DFT; and

means for determining the heart rate based on a comparison of the first estimate and the second estimate.

2. The system according to claim 1, wherein the second time period is greater than and includes the first time period.

3. The system according to claim 2, further comprising means for determining whether the data samples of the first signal over the second time period pass the set of confidence checks, wherein the DFT is computed upon determination that the data samples of the first signal over the second time period pass the set of confidence checks.

4. The system according to claim 3, further comprising means for computing a series of autocorrelations of the data samples of the first signal within the second time period, each autocorrelation computed over a period that is equal to the first time period,

wherein the set of confidence checks further comprises determining whether the outcome of the each autocorrelation satisfies an autocorrelation-related criteria.

5. The system according to claim 4, wherein computing the autocorrelation of the data samples of the first signal over the second time period comprises performing iterations of adding a contribution due to a new data sample to an autocorrelation of a previous iteration and subtracting a contribution due to a former data sample for each new data sample of data samples between a last data sample of the data samples of the first signal over the first time period and a last data sample of the data samples of the first signal over the second time period.

6. The system according to claim 1, wherein the set of confidence checks is a first set of confidence checks, and the system further includes means for determining whether the data samples of the first signal over the first time period pass a second set of confidence checks, wherein the DFT is computed upon determination that the data samples of the first signal over the first time period pass the second set of confidence checks.

7. The system according to claim 6, further comprising means for determining an expected number of local peaks in the autocorrelation of the data samples of the first signal over the first time period based on the first heart rate estimate,

wherein the second set of confidence checks comprises determining whether the expected number of local peaks is equal to a number of local peaks present in the autocorrelation of the data samples of the first signal over the first time period.

8. The system according to claim 1, further comprising means for computing a series of autocorrelations of the data samples of the first signal within the second time period, each autocorrelation computed over a period that is equal to the first time period,

wherein the set of confidence checks further comprises determining whether the outcome of the each autocorrelation satisfies autocorrelation-related criteria or a predetermined number of autocorrelations satisfy autocorrelation-related criteria.

9. The system according to claim 8, wherein determining whether the outcome of the each autocorrelation satisfies the autocorrelation-related criteria is carried out only for a portion of the each autocorrelation for which an autocorrelation lag is between 0.25 seconds and 2 seconds.

10. The system according to claim 9, wherein the autocorrelation-related criteria comprise determining whether an autocorrelation peak that occurs first within the portion of the each autocorrelation is larger than any other autocorrelation value within the portion of the each autocorrelation.

11. The system according to claim 9, wherein the autocorrelation-related criteria comprises determining whether autocorrelation peaks occur within the portion of the each autocorrelation in order of decreasing values.

12. The system according to claim 1, further comprising:

means for using the data samples of the second signal to determine a parameter indicative of the motion of the heartbeat sensor,

wherein the motion-related criteria comprises determining whether the parameter is below a predefined threshold.

13. The system according to claim 1, further comprising means for determining that an absolute value of the DFT has a local peak in a feasible range of heart rates.

14. The system according to claim 1, wherein determining the heart rate based on the comparison of the first estimate and the second estimate comprises:

determining whether the first estimate and the second estimate are within a predefined number of beats per minute (bpm) from one another, and

upon determining that the first estimate and the second estimate are within the predefined number of bpm from one another, determining the heart rate as equal to one of the first estimate, the second estimate, or an estimate computed based on the first estimate and the second estimate.

15. The system according to claim 1, further comprising means for decimating the data samples of the first signal over the second time period prior to computing the DFT and computing the DFT on the decimated data samples.

16. The system according to claim 1, further comprising:

means for processing the first signal and/or the second signal with a pre-processing filter to substantially attenuate signal content outside of a reasonable frequency band of interest corresponding to a range of reasonable heart rates.

17. A system for determining a heart rate from a heartbeat signal present in a signal generated by a heartbeat sensor, the system comprising:

means for computing an autocorrelation of data samples of the signal over a first time period;

means for determining whether the data samples of the signal over the first time period satisfy a motion-related criteria;

means for computing, when the data samples of the signal over the first time period satisfy the motion-related criteria, a first estimate of the heart rate based on the autocorrelation;

means for computing a Fourier Transform (FT) for data samples of the signal over a second time period;

means for computing a second estimate of the heart rate based on the FT; and

means for determining the heart rate based on a comparison of the first estimate and the second estimate.

18. The system according to claim 17, wherein the means for determining whether the data samples of the signal over the first time period satisfy the motion-related criteria include means for determining whether the data samples of the signal over the first time period satisfy the motion-related criteria based on a second signal indicative of a motion of the heartbeat sensor with respect to one or more directions.

19. A system for determining a heart rate from a heartbeat signal present in a first signal generated by a heartbeat sensor, the system comprising:

means for computing an autocorrelation of data samples of the signal over a first time period;

means for determining whether the data samples of the signal over the first time period satisfy a motion-related criteria;

means for computing, when the data samples of the signal over the first time period satisfy the motion-related criteria, a first estimate of the heart rate based on the autocorrelation;

means for computing a Fourier Transform (FT) for data samples of the signal over a second time period, wherein the second time period includes the first time period;

means for computing a second estimate of the heart rate based on the FT; and

means for determining the heart rate based the first estimate and the second estimate.

20. The system according to claim 19, wherein determining the heart rate based on the first estimate and the second estimate includes:

determining whether the first estimate and the second estimate are within a predefined number of beats per minute (bpm) from one another, and

upon determining that the first estimate and the second estimate are within the predefined number of bpm from one another, determining the heart rate as equal to one of the first estimate, the second estimate, or an estimate computed based on the first estimate and the second estimate.