Removing motion-related artifacts in heart rate measurement systems using iterative mask estimation in frequency-domain

Patent metadata

Field	Value
Publication number	US-2016354038-A1
Application number	US-201615170373-A
Country	US
Kind code	A1
Filing date	Jun 1, 2016
Priority date	Jun 3, 2015
Publication date	Dec 8, 2016
Grant date	—

Abstract

Heart rate monitors are plagued by noisy photoplethysmography (PPG) data, which makes it difficult for the monitors to output a consistently accurate heart rate reading. Noise is often caused by motion. Using known methods for processing accelerometer readings that measure movement to filter out some of this noise may help, but not always. The present disclosure describes an improved filtering approach, referred to herein as an iterative frequency-domain mask estimation technique, based on using frequency-domain representation (e.g. STFT) of PPG data and accelerometer data for each accelerometer channel to generate filters for filtering the PPG signal from motion-related artifacts prior to tracking frequency of the heartbeat (heart rate). Implementing this technique leads to more accurate heart rate measurements.

First claim

What is claimed is:

1 . A computer-implemented method for assisting separation of a heartbeat signal present in a first signal generated by a heartbeat sensor, the method comprising:
processing the first signal to compute a time-frequency-domain representation (p obs (f,t)) of the first signal; processing a second signal to compute a time-frequency-domain representation (p accx (f,t)) of the second signal, the second signal indicative of a motion of the heartbeat sensor with respect to a first direction (x); initializing a first source model for a first source (s 1 ), the first source representing a source of the heartbeat signal; initializing a second source model for a second source (s 2 ), the second source representing a source of a contribution to the first signal due to the motion of the heartbeat sensor with respect to the first direction, wherein the second source model is initialized based on the time-frequency-domain representation (p accx (f,t)) of the second signal; performing a plurality of iterations of modifying one or more parameters of the first source model and/or one or more parameters of the second source model based on the time-frequency-domain representation of the first signal; and following the plurality of iterations, computing an estimate of the heartbeat signal based on the one or more parameters of the first source model and/or the one or more parameters of the second source model.

2 . The method according to claim 1 , wherein each iteration of the plurality of iterations comprises:
using the first source model to generate an estimate of a time-frequency-domain representation (q(f,t|s 1 )) for the first source, using the second source model to generate an estimate of a time-frequency-domain representation (q(f,t|s 2 )) for the second source, and based on the estimate of the time-frequency-domain representation for the first source and the estimate of the time-frequency-domain representation for the second source, computing a first mask function q(s 1 |f,t) for separating the heartbeat signal from the first signal and computing a second mask function q(s 2 |f,t) for separating, from the first signal, the contribution due to motion of the heartbeat sensor with respect to the first direction.

3 . The method according to claim 2 , wherein said each iteration of the plurality of iterations further comprises:
applying the first mask function to the time-frequency-domain representation of the first signal to generate an updated estimate of a time-frequency-domain representation (p(f,t|s 1 )) for the first source, applying the second mask function to the time-frequency-domain representation of the first signal to generate an updated estimate of a frequency-domain representation for the second source (s2), applying the updated estimate of the time-frequency-domain representation for the first source to update one or more parameters of the first source model, and applying the updated estimate of the frequency-domain representation for the second source to update one or more parameters of the second source model.

4 . The method according to claim 3 , wherein:
applying the first mask function comprises performing a point-wise multiplication of the first mask function and the time-frequency-domain representation of the first signal, and applying the second mask function comprises performing a point-wise multiplication of the second mask function and the time-frequency-domain representation of the first signal.

5 . The method according to claim 1 , wherein:
the first source model comprises a Non-negative Tensor Factorization (NTF) source model or a one-peak source model, and/or the second source model comprises the NTF source model or a constant source model.

6 . The method according to claim 1 , further comprising:
processing a third signal to compute a time-frequency-domain representation (p accy (f,t)) of the third signal, the third signal indicative of a motion of the heartbeat sensor with respect to a second direction (y); and initializing a third source model for a third source (s3), the third source representing a source of a contribution to the first signal due to the motion of the heartbeat sensor with respect to the second direction, wherein the third source model is initialized based on the time-frequency-domain representation (p accy (f,t)) of the third signal, wherein the estimate of the heartbeat signal is computed based on the one or more parameters of the first source model and/or the one or more parameters of the second source model and/or the one or more parameters of the third source model.

7 . The method according to claim 6 , wherein the plurality of iterations further include modifying one or more parameters of the third source model based on the frequency-domain representation of the first signal.

8 . The method according to claim 1 , wherein:
the heartbeat sensor comprises an optical sensor configured to generate the first signal by detecting light at a first frequency or range of frequencies, and the method further comprises performing processing applied to the first signal to a further signal generated by a further heartbeat sensor, the further heartbeat sensor comprising an optical sensor configured to generate the further signal by detecting light at a second frequency or range of frequencies, and the estimate of the heartbeat signal is computed further based on the further signal.

9 . The method according to claim 1 , wherein the plurality of iterations are performed until a predefined maximum number of iterations is reached or until a divergence between the estimate of the time-frequency-domain representation for the first source (i.e. q(f,t|s1)) generated by the first source model and an updated estimate of the frequency-domain representation for the first source (i.e. p(f,t|s1)) satisfies a predefined criteria.

10 . The method according to claim 9 , wherein the predefined criteria is based on a Kullback-Leibler (KL) divergence.

11 . The method according to claim 1 , wherein performing the plurality of iterations comprises carrying out an iterative algorithm based on a probabilistic inference approach.

12 . The method according to claim 1 , further comprising:
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 rate frequencies.

13 . The method according to claim 12 , wherein:
the pre-processing filter is a low-pass filter or a band-pass filter; and the reasonable frequency band of interest comprises frequencies between 0.5 Hertz and 4 Hertz.

14 . The method according to claim 1 , further comprising:
applying a filter or mask to or removing a portion of the first signal indicative of a saturation condition of the heartbeat sensor.

15 . The method according to claim 1 , further comprising:
generating a time-frequency-domain representation of the estimate of the heartbeat signal computed based on the one or more parameters of the first source model and/or the one or more parameters of the second source model; and tracking one or more contours present in the time-frequency-domain representation to track the heartbeat signal.

16 . The method according to claim 1 , wherein the heartbeat sensor comprises one or more of the following: optical sensor, audio sensor, capacitive sensor, magnetic sensor, chemical sensor, humidity sensor, moisture sensor, pressure sensor, and biosensor.

17 . The method accord

Assignees

• Analog Devices Inc

Inventors

• Demirtas Sefa
• Adams Robert
• Bernstein Jeffrey G
• Stein Noah Daniel
• Wingate David

Classifications

• A61B2562/0219

  Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches · CPC title

• A61B5/725

  using specific filters therefor, e.g. Kalman or adaptive filters (specific diagnostics methods using using bioelectric or biomagnetic signals A61B5/316) · CPC title

• A61B5/721 Primary

  using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured · CPC title

• A61B5/02438

  with portable devices, e.g. worn by the patient · CPC title

• A61B5/02416

  using photoplethysmograph signals, e.g. generated by infrared radiation (A61B5/14552 takes precedence) · CPC title