US20170164847A1 - Reducing Motion Induced Artifacts in Photoplethysmography (PPG) Signals - Google Patents
Reducing Motion Induced Artifacts in Photoplethysmography (PPG) Signals Download PDFInfo
- Publication number
- US20170164847A1 US20170164847A1 US14/970,434 US201514970434A US2017164847A1 US 20170164847 A1 US20170164847 A1 US 20170164847A1 US 201514970434 A US201514970434 A US 201514970434A US 2017164847 A1 US2017164847 A1 US 2017164847A1
- Authority
- US
- United States
- Prior art keywords
- signal
- ppg
- axis acceleration
- heart rate
- motion compensated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 230000033001 locomotion Effects 0.000 title claims abstract description 127
- 238000013186 photoplethysmography Methods 0.000 title description 4
- 230000001133 acceleration Effects 0.000 claims abstract description 95
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000009532 heart rate measurement Methods 0.000 claims abstract description 12
- 230000003044 adaptive effect Effects 0.000 claims description 50
- MKGHDZIEKZPBCZ-ULQPCXBYSA-N methyl (2s,3s,4r,5r,6r)-4,5,6-trihydroxy-3-methoxyoxane-2-carboxylate Chemical compound CO[C@H]1[C@H](O)[C@@H](O)[C@H](O)O[C@@H]1C(=O)OC MKGHDZIEKZPBCZ-ULQPCXBYSA-N 0.000 description 21
- 238000001228 spectrum Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 230000000875 corresponding effect Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 210000000707 wrist Anatomy 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 210000000624 ear auricle Anatomy 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000001061 forehead Anatomy 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 230000003278 mimic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000005022 packaging material Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1118—Determining activity level
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
- A61B5/721—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Definitions
- Embodiments of the present invention generally relate to heart rate measurement and more specifically relate to reducing motion induced artifacts in photoplethysomography (PPG) signals used for heart rate measurement.
- PPG photoplethysomography
- Photoplethysmography is a well-known optical method for measuring the change in blood flow.
- a light source at a particular wavelength typically, red, infra-red or green
- the change in blood flow can be used to compute vital statistics such as heart rate, respiration rate and other health related statistics.
- vital statistics such as heart rate, respiration rate and other health related statistics.
- Many consumer fitness devices use PPG to estimate the heart rate of the person wearing the device.
- PPG based heart rate estimation during motion can be adversely affected by motion artifacts in the PPG signal.
- the motion artifacts may be caused by, for example, hemodynamic effects, tissue deformation, and sensor movement relative to the skin.
- Embodiments of the present invention relate to methods, apparatus, and computer readable media for reducing motion induced artifacts in photoplethysomography (PPG) signals used for heart rate measurement.
- a method for heart rate measurement in a photoplethysmograph (PPG) heart rate monitor device includes receiving a PPG signal from a PPG sensor of the heart rate monitor device, receiving an X-axis acceleration signal, a Y-axis acceleration signal, and a Z-axis acceleration signal from an accelerometer of the heart rate monitor device, performing motion compensation on the PPG signal, wherein a first motion compensated PPG signal is generated with reference to the X-axis acceleration signal, a second motion compensated PPG signal is generated with reference to the Y-axis acceleration signal, and a third motion compensated PPG signal is generated with reference to the Z-axis acceleration signal, combining the first motion compensated PPG signal, the second motion compensated PPG signal, and the third motion compensated PPG signal to generate a fourth motion compensate
- a photoplethysmograph (PPG) heart rate monitor device includes a PPG sensor component configured to generate an analog PPG signal, an accelerometer configured to generate an X-axis acceleration signal, a Y-axis acceleration signal, and a Z-axis acceleration signal, a motion compensation component coupled to the PPG sensor to receive the PPG signal and coupled to the accelerometer to receive the X-axis acceleration signal, the Y-axis acceleration signal, and the Z-axis acceleration signal, wherein the motion compensation component is configured to generate a first motion compensated PPG signal with reference to the X-axis acceleration signal, a second motion compensated PPG signal with reference to the Y-axis acceleration signal, and a third motion compensated PPG signal with reference to the Z-axis acceleration signal, a combiner configured to combine the first motion compensated PPG signal, the second motion compensated PPG signal, and the third motion compensated PPG signal to generate a fourth motion compensated PPG signal, wherein a first weight is applied to
- FIG. 1 and FIG. 2 are graphs
- FIG. 3 is a block diagram of an example photoplethysmograph (PPG) heart rate monitor device
- FIG. 4 and FIG. 5 are block diagrams illustrating an embodiment of motion compensation of a PPG signal and heart rate tracking that may be performed in a suitably configured heart rate monitor device such as that of FIG. 3 ;
- FIG. 6 is a flow diagram of a method for motion compensation of a PPG signal and heart rate tracking that may be performed in a suitably configured heart rate monitor device such as that of FIG. 3 .
- Some current approaches to improving heart rate estimation when motion artifacts may be present in the PPG signal attempt to remove the motion component in the PPG signal using information from an external signal reference such as an accelerometer.
- an external signal reference such as an accelerometer.
- a spectrum subtraction technique is proposed that removes the spectrum of the acceleration data from that of the PPG signal prior to heart rate estimation.
- FIG. 1 is a graph illustrating this assumption with data from a wrist worn device with a PPG sensor and a three-axis accelerometer sensor in which the hand motion is performed in the direction of the Y-axis of the accelerometer.
- This graph shows that the Y-axis acceleration signal is inversely proportional to the PPG signal when there is motion in the Y-axis.
- the amount of mutual information between a PPG signal and an accelerometer signal can, in a first order approximation, be captured by computing the Pearson correlation coefficient defined as
- FIG. 2 is a graph illustrating the Pearson correlation for the X-axis acceleration signal (accX), the Y-axis acceleration signal (accY) and the z-axis acceleration signal (accZ) of an accelerometer with respect to the PPG signal when the motion is predominantly along the Y-axis.
- the Pearson correlation coefficient associated with the Y-axis acceleration signal is consistently larger than that associated with the X-axis and Z-axis as the motion of the hand is predominantly in the direction of the Y-axis.
- HeartRateTracker computes the heart rate from the power spectrum
- PPG accX is the NLMS filter output with an input of the PPG signal with the accX signal as the reference signal
- PPG accY is the NLMS filter output with an input of the PPG signal with the accY signal as the reference signal
- PPG accZ is the NLMS filter output with an input of the PPG signal with the accZ signal as the reference signal.
- a simple version of HeartRateTracker associates the heart rate with the frequency index of the largest spectrum peak in S(f).
- Embodiments of the disclosure provide for reducing motion induced artifacts in PPG signals with reduced complexity and improved performance over the prior art.
- coherent combining of NLMS filter outputs (or other suitable adaptive filter outputs) is performed followed by a single Fourier Transform (FT).
- FT Fourier Transform
- the motion compensated power spectrum of the PPG signal may be computed as
- a, b, and c are weights applied to the outputs of suitable adaptive filters.
- the adaptive filters are NLMS filters.
- Other adaptive filters such as, for example recursive least squares (RLS) filters.
- NLMS filters are described in S. Haykin, “Adaptive Filter Theory”, pp. 320-324, Pearson Education, Inc., 2002, which is incorporated by reference herein.
- RLS filters are described in Haykin, p. 443, which is incorporated by reference herein. Any suitable technique for determining the weight values of the power spectrum computation may be used. Some examples of suitable techniques are now described.
- the values of the weights are based on the amount of correlation of the PPG signal with the respective acceleration signals.
- the weighting is based on the idea that the higher the correlation of the PPG signal with the corresponding acceleration signal, the more motion compensated the output of the corresponding NLMS filter should be.
- the correlation coefficients between each acceleration signal and the PPG signal are computed in a T second window and the weights are then computed as a function of the correlation coefficients. In one embodiment, the maximum correlation coefficients are used.
- the correlation r 1accX between the accX signal and the PPG signal at time n+m/2 can be computed as
- the maximum of the correlation coefficients in the T second window may be computed as
- r accX max( r 1accX ( n ) . . . r 1accX ( n+m ⁇ 1)).
- the maximum of the correlation coefficients between the accY signal and the PPG signal, r accY , and between the accZ signal and the PPG signal, r accX may be similarly computed.
- the weight values may then be determined as a function of r accX , r accY , and r accZ . For example, in one embodiment the weight values may be determined as
- T is typically chosen to be the size of the data over which the FT is performed and is of the order of a few seconds.
- the values of the weights are based on the predicted motion artifact energy, i.e., interference energy, removed by the NLMS filter over a window of T seconds.
- the NLMS filter output at time n when the sensor reference is the X-axis acceleration signal is given by
- w denotes the NLMS filter coefficients for accX
- the predicted interference energy removed over a window of T seconds (Tsec) by the NLMS filter corresponding to accX may thus be determined as per
- IntEn accX ⁇ Tsec ⁇ ⁇ ⁇ Interf accX ⁇ ( n ) ⁇ 2 .
- the predicted interference energy corresponding to accY, IntEn accY , and accZ, IntEn accZ may be similarly determined.
- the values of the weights may then be determined as a function of IntEn accX , IntEn accY , and IntEn accY .
- the weight values may be determined as
- T is typically chosen to be the size of the data over which the FT is performed and is of the order of a few seconds.
- one of the signals is selected according to some suitable selection criteria.
- the decision rules may be written as
- the decision rules may be written as
- the PPG signal is filtered when the user of the device is in motion and not filtered otherwise. Filtering the PPG signal to remove motion artifacts when the user of the device is not in motion may degrade the PPG signal as motion artifacts are not present.
- Any suitable technique may be used to determine if the user is in motion. As shown in Table 1, one suitable technique that may be used is to compute the variance of the sum of squares of the acceleration signals over a window of T1 seconds having m1 samples and compare the variance to a predetermined threshold Thresh. If the variance exceeds the threshold, the user is determined to be in motion at time instant n; otherwise the user is determined to be stationary. Any suitable threshold value may be used. In some embodiments, the threshold value is determined empirically based on experiments performed on the acceleration signals to determine a good representation of stationary.
- FIG. 3 is a block diagram of an example embedded photoplethysmograph (PPG) heart rate monitor device 300 configured to perform motion compensation on PPG signals as described herein.
- the PPG heart rate monitor device 300 includes a PPG sensor component 304 configured to capture a PPG signal from the person wearing the monitor device 300 and a processing component 302 configured to process the PPG signal to measure the heart rate of a person wearing the monitor device.
- the PPG sensor component 304 may be any suitable PPG sensor.
- the PPG sensor component includes one or more three light emitting diodes (LEDs) and a photodiode arranged to the wavelengths of the light from the one or more LEDs.
- the PPG sensor component 304 includes three light emitting diodes (LEDs): a green light emitting diode (LED), a red LED, and an infrared LED.
- the processing component 302 includes a microcontroller (MCU) 308 , an inertial measurement unit (IMU) 312 , a Bluetooth transceiver 314 , an analog front end (AFE) 310 , and memory 318 .
- the MCU 308 may be any suitable microcontroller, such as, for example, an MSP430 device available from Texas Instruments, Inc.
- the MCU 308 includes memory (one example of a computer readable medium) that can be used to store software instructions needed to perform motion compensation on the PPG signals as described herein and to perform heart rate measurement based on the motion compensated PPG signals. Embodiments of motion compensation that may be executed by the MCU 308 are described in reference to FIGS. 4, 5, and 6 .
- Memory 318 may be any suitable memory device or devices, such as, for example, one or more ferroelectric random access memory (FRAM) devices.
- Memory 318 is coupled to the MCU 308 and may be used, for example, to store data used in the execution of the heart rate monitoring and to store a heart rate history.
- FRAM ferroelectric random access memory
- the IMU 312 is coupled to the MCU 308 via an interface provided by the MCU 308 and provides data regarding the motion of a person wearing the monitor device 300 to the MCU 308 . More specifically, the IMU 312 includes a three axis accelerometer to provide the X-axis acceleration signal (accX), the Y-axis acceleration (accY) signal and the z-axis acceleration (accZ) signal to the MCU 308 .
- the IMU 312 may be any suitable IMU device, such as, for example, an IMU device available from Invensense, Inc.
- the Bluetooth transceiver 314 is coupled to the MCU 308 and may be used, for example, to transfer the computed heart rate measurement of a person wearing the monitor device 300 to another device for display and/or further processing.
- the heart rate measurement may be transmitted to a smart phone or other personal digital assistant, a laptop computer, a desktop computer, a medical monitoring device, etc.
- the Bluetooth transceiver 314 may be any suitable Bluetooth device, such as, for example, a CC25xx Bluetooth system-on-a-chip (SOC) available from Texas Instruments, Inc.
- the USB (universal serial bus) transceiver 316 is coupled to the MCU 308 and may be used, for example, to transfer the computed heart rate measurement of a person wearing the monitor device 300 to another device for display and/or further processing.
- the heart rate measurement may be transmitted to a smart phone or other personal digital assistant, a laptop computer, a desktop computer, a medical monitoring device, etc.
- the associated USB port may also be used to charge a battery (not shown) in the monitor device 300 .
- the AFE 310 is coupled to the MCU 308 and the PPG sensor component 304 and provides a control interface between the MCU 308 and PPG sensor component 304 .
- the AFE 310 includes functionality to receive a PPG signal from the PPG sensor component 304 , convert the signal to a digital signal, apply analog signal conditioning such as changing the gain of the signal, and provide the digital PPG signal to the MCU 308 for heart rate measurement.
- the AFE 310 also includes functionality to drive current to one or more LEDs in the PPG sensor 304 as directed by the MCU 308 . The particular current to be used and the timing of when to apply the current may be controlled by the MCU 308 .
- the AFE 310 may be any suitable AFE device, such as, for example, the AFE4400 provided by Texas Instruments, Inc.
- the PPG heart rate monitor device 300 includes a display component (not shown) coupled to the MCU 308 .
- the display component may be used to display the heart rate of a person wearing the monitor device 300 .
- FIG. 4 and FIG. 5 are block diagrams illustrating an embodiment of motion compensation of a PPG signal and heart rate tracking that may be performed in a suitably configured heart rate monitor device, such as, for example, the heart rate monitor device 300 of FIG. 3 .
- a suitably configured heart rate monitor device such as, for example, the heart rate monitor device 300 of FIG. 3 .
- Embodiments are described in reference to heart rate monitor device 300 for simplicity of explanation.
- embodiments are described assuming the various components are implemented as software executed by a processor such as the MCU 308 .
- the FT component 410 may be an FFT hardware accelerator tightly coupled with the central processing unit of the MCU 308 .
- one or more of the components may be implemented in an application specific integrated circuit (ASIC).
- ASIC application specific integrated circuit
- the bandpass filter 402 receives the PPG signal from the AFE 310 is filtered 402 and outputs a filtered PPG signal.
- the bandpass filter 404 receives the X-axis acceleration signal (accX), the Y-axis acceleration signal (accY) and the z-axis acceleration signal (accZ) from the IMU 312 and output filtered versions of these acceleration signals.
- the bandpass filters 402 , 404 are filters with the same in-band and out-of-band characteristics.
- the motion compensation component (MCC) 406 receives the filtered PPG signal and the filtered acceleration signals and performs motion compensation on the PPG signal based on the acceleration signals. More specifically, the MCC 406 performs motion compensation on the PPG signal with respect to each of the acceleration signals and generates three output signals, PPG accX , PPG accY , and PPG accZ . The motion compensation is described in more detail in reference to FIG. 5 .
- the MCC 406 also outputs weighting data that may be used by the combiner 408 to determine weights for combining the three motion filtered PPG signals, e.g., indications of the correlation of the PPG signal with each acceleration signal and/or the predicted interference energy removed from the PPG signal to generate each of the signals output by the MCC 406 .
- the combiner 408 combines the three signals output by the MCC 406 as per
- the FT component 410 applies an FFT to the final motion compensated PPG signal to convert the signal from the time domain to the frequency domain.
- the resulting signal is then provided to the heart rate tracker 412 to be used for heart rate estimation.
- the block diagram of FIG. 5 illustrates an embodiment of the motion compensation performed by the MCC 406 .
- the correlation indicator component 516 receives the filtered PPG signal and the filtered acceleration signals and generates weighting data for the combiner 408 that is indicative of the correlation between the PPG signal and each of the acceleration signals. That is, the correlation indicator component 516 computes the correlation coefficients between each acceleration signal and the PPG signal in T second window. Computation of correlation coefficients is previously described herein.
- the output of the correlation indicator component 516 is the correlation coefficients for each combination of PPG signal and acceleration signal and the combiner 408 uses the correlation coefficients to determine the weight values used to combine the motion compensated PPG signals PPG accX , PPG accY , and PPG accZ .
- the output of the correlation indicator component 516 is the maximum of the correlation coefficients for each combination of PPG signal and acceleration signal and the combiner 408 uses the maximum coefficients to determine the weight values. Computation of the maximum of the correlation coefficients is previously described herein.
- the interference energy detection component 514 receives the filtered PPG signal, the filtered acceleration signals, and the current filter coefficients from each adaptive filter 502 , 504 , 506 and generates weighting data for the combiner 408 that is indicative of the predicted interference energy removed from the PPG signal by each adaptive filter 502 , 504 , 506 . That is, the interference energy detection component 514 computes the predicted interference energy removed from the PPG signal by each adaptive filter 502 , 504 , 506 in T second window. Computation of predicted interference energy is previously described herein.
- the stationary detector component 518 receives the filtered acceleration signals and generates an indication of whether or not the person wearing the monitor device 300 is stationary, i.e., a stationary indicator.
- the stationary detector component 518 may use any suitable technique to determine if the person is in motion. An example of one such suitable technique is previously described herein.
- the adaptive filters 502 , 504 , 506 each receive the filtered PPG signal and a respective acceleration signal accX, accY, accZ and output respective motion compensated PPG signals PPG accX , PPG accY , and PPG accZ .
- each adaptive filter 502 , 504 , 506 is a normalized least mean squares (NLMS) filter in which the input signal is the PPG signal and the reference signal is the respective acceleration signal.
- NLMS normalized least mean squares
- an NLMS filter is a variant of a least mean squares (LMS) filter which is a class of adaptive filter used to mimic a desired filter by finding the filter coefficients that relate to producing the least mean squares of the error signal, i.e., the difference between the desired signal and the actual signal.
- LMS least mean squares
- other adaptive filters may be used, such as, for example recursive least squares (RLS) filters.
- the adaptive filters 502 , 504 , 506 also receive the stationary indicator from the stationary detector component 518 . If the stationary indicator indicates that the person is stationary, the adaptive filters 502 , 504 , 506 do not filter the input PPG signal and the output of each adaptive filter is the input PPG signal. Further, the combiner 408 performs no combining and passes one of the input PPG signals as the output PPG signal. In some embodiments, if the stationary indicator indicates that the person is stationary, the weights of the adaptive filter are not changed, i.e., the weights remain at their current values. Thus, the filter weights are only updated when the person is in motion. In some embodiments, if the stationary indicator indicates that the person is stationary, the weights of the adaptive filter are reset to their initial values.
- FIG. 5 includes both a correlation indicator component 516 and an interference energy detection component 514 for ease of description.
- a correlation indicator component 516 and an interference energy detection component 514 for ease of description.
- one of these components is not present or is otherwise not used.
- neither component is present, e.g., in embodiments in which the combiner uses a weighting scheme in which the weight values are constants.
- FIG. 6 is a flow diagram of a method for motion compensation of a PPG signal and heart rate tracking that may be performed in a suitably configured heart rate monitor device, such as, for example, the heart rate monitor device 300 of FIG. 3 .
- a suitably configured heart rate monitor device such as, for example, the heart rate monitor device 300 of FIG. 3 .
- Embodiments are described in reference to heart rate monitor device 300 for simplicity of explanation.
- One of ordinary skill in the art will understand embodiments for other hardware architectures.
- the method may be performed continuously while the heart rate monitor device is turned on.
- a PPG signal is received 600 by the MCU 308 from the PPG sensor component 304 via the AFE 310 and is filtered using a bandpass filter.
- Accelerometer signals i.e., an X-axis acceleration signal, a Y-axis acceleration signal, and a Z-axis acceleration signal are also received 602 by the MCU 308 from the IMU 312 and are filtered using a bandpass filter.
- Motion compensation is then performed 604 on the filtered PPG signal using the filtered acceleration signals as reference signals to generate three motion compensated PPG signals.
- the filtered PPG signal is adaptively filtered using each of the X-axis acceleration signal, the Y-axis acceleration signal, and the Z-axis acceleration signal as a reference signal to generate three motion compensated PPG signals.
- the adaptive filter may be, for example, an NLMS filter or an RLS filter.
- the three motion compensated PPG signals are then coherently combined 606 as per
- a single FT is then applied 608 to the final motion compensated PPG signal to generate a frequency domain PPG signal.
- the frequency domain signal is then used to estimate 610 the heart rate.
- PPG heart rate monitor devices housed in form factors suitable for monitoring heart rate in various locations on the body of a person such as, for example, the forehead, an ear lobe, a finger, an upper arm, a wrist, or a chest.
- the stationary detector component is not present and the IMU determines whether the person is in motion and generates the value of the stationary indicator.
- the methods described herein may be implemented in hardware, software, firmware, or any combination thereof. If completely or partially implemented in software, the software may be executed in one or more processors, such as a microcontroller, microprocessor, application specific integrated circuit (ASIC), field programmable gate array (FPGA), or digital signal processor (DSP).
- the software instructions may be initially stored in a computer-readable medium and loaded and executed in the processor. In some cases, the software instructions may also be sold in a computer program product, which includes the computer-readable medium and packaging materials for the computer-readable medium. In some cases, the software instructions may be distributed via removable computer readable media, via a transmission path from computer readable media on another digital system, etc. Examples of computer-readable media include non-writable storage media such as read-only memory devices, writable storage media such as disks, flash memory, memory, or a combination thereof.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Computer Networks & Wireless Communication (AREA)
- Mathematical Physics (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
A method for heart rate measurement in a photoplethysmograph (PPG) heart rate monitor device is provided that includes performing motion compensation on a PPG signal wherein a motion compensated PPG signal PPGaccX is generated with reference to an X-axis acceleration signal, a motion compensated PPG signal PPGaccY is generated with reference to a Y-axis acceleration signal, and a motion compensated PPG signal PPGaccZ is generated with reference to a Z-axis acceleration signal, combining PPGaccX, PPGaccY, and PPGaccZ to generate a final motion compensated PPG signal, wherein a first weight is applied PPGaccX, a second weight is applied to PPGaccY, and a third weight is applied to PPGaccZ, performing a single Fourier Transform (FT) on the final motion compensated PPG signal to generate a frequency domain PPG signal; and estimating a heart rate based on the frequency domain PPG signal.
Description
- Field of the Invention
- Embodiments of the present invention generally relate to heart rate measurement and more specifically relate to reducing motion induced artifacts in photoplethysomography (PPG) signals used for heart rate measurement.
- Description of the Related Art
- Photoplethysmography (PPG) is a well-known optical method for measuring the change in blood flow. In PPG, a light source at a particular wavelength, typically, red, infra-red or green, is kept at a skin surface and the reflected light is measured and related to blood flow. The change in blood flow can be used to compute vital statistics such as heart rate, respiration rate and other health related statistics. Many consumer fitness devices use PPG to estimate the heart rate of the person wearing the device. PPG based heart rate estimation during motion can be adversely affected by motion artifacts in the PPG signal. The motion artifacts may be caused by, for example, hemodynamic effects, tissue deformation, and sensor movement relative to the skin.
- Embodiments of the present invention relate to methods, apparatus, and computer readable media for reducing motion induced artifacts in photoplethysomography (PPG) signals used for heart rate measurement. In one aspect, a method for heart rate measurement in a photoplethysmograph (PPG) heart rate monitor device is provided that includes receiving a PPG signal from a PPG sensor of the heart rate monitor device, receiving an X-axis acceleration signal, a Y-axis acceleration signal, and a Z-axis acceleration signal from an accelerometer of the heart rate monitor device, performing motion compensation on the PPG signal, wherein a first motion compensated PPG signal is generated with reference to the X-axis acceleration signal, a second motion compensated PPG signal is generated with reference to the Y-axis acceleration signal, and a third motion compensated PPG signal is generated with reference to the Z-axis acceleration signal, combining the first motion compensated PPG signal, the second motion compensated PPG signal, and the third motion compensated PPG signal to generate a fourth motion compensated PPG signal, wherein a first weight is applied to the first motion compensated PPG, a second weight is applied to the second motion compensated PPG signal, and a third weight is applied to the third motion compensated PPG signal, performing a single Fourier Transform (FT) on the fourth motion compensated PPG signal to generate a frequency domain PPG signal, and estimating a heart rate based on the frequency domain PPG signal.
- In one aspect, a photoplethysmograph (PPG) heart rate monitor device is provided that includes a PPG sensor component configured to generate an analog PPG signal, an accelerometer configured to generate an X-axis acceleration signal, a Y-axis acceleration signal, and a Z-axis acceleration signal, a motion compensation component coupled to the PPG sensor to receive the PPG signal and coupled to the accelerometer to receive the X-axis acceleration signal, the Y-axis acceleration signal, and the Z-axis acceleration signal, wherein the motion compensation component is configured to generate a first motion compensated PPG signal with reference to the X-axis acceleration signal, a second motion compensated PPG signal with reference to the Y-axis acceleration signal, and a third motion compensated PPG signal with reference to the Z-axis acceleration signal, a combiner configured to combine the first motion compensated PPG signal, the second motion compensated PPG signal, and the third motion compensated PPG signal to generate a fourth motion compensated PPG signal, wherein a first weight is applied to the first motion compensated PPG, a second weight is applied to the second motion compensated PPG signal, and a third weight is applied to the third motion compensated PPG signal, a FT component configured to perform a single Fast Fourier Transform (FFT) on the fourth motion compensated PPG signal to generate a frequency domain PPG signal, and a heart rate tracker component configured to estimate a heart based on the frequency domain PPG signal.
- Particular embodiments will now be described, by way of example only, and with reference to the accompanying drawings:
-
FIG. 1 andFIG. 2 are graphs; -
FIG. 3 is a block diagram of an example photoplethysmograph (PPG) heart rate monitor device; -
FIG. 4 andFIG. 5 are block diagrams illustrating an embodiment of motion compensation of a PPG signal and heart rate tracking that may be performed in a suitably configured heart rate monitor device such as that ofFIG. 3 ; and -
FIG. 6 is a flow diagram of a method for motion compensation of a PPG signal and heart rate tracking that may be performed in a suitably configured heart rate monitor device such as that ofFIG. 3 . - Specific embodiments of the invention will now be described in detail with reference to the accompanying figures. Like elements in the various figures are denoted by like reference numerals for consistency.
- As previously mentioned, motion artifacts in a PPG signal can adversely affect heart rate estimation. Some current approaches to improving heart rate estimation when motion artifacts may be present in the PPG signal attempt to remove the motion component in the PPG signal using information from an external signal reference such as an accelerometer. For example, in H. Fukushima, et al., “Estimating Heart Rate using Wrist-Type Photoplethysmography and Acceleration Sensor While Running,” 34th Annual Internal Conference of the IEEE Engineering in Medicine and Biology Society, pp. 2901-2904, Aug. 28-Sep. 1, 2012, a spectrum subtraction technique is proposed that removes the spectrum of the acceleration data from that of the PPG signal prior to heart rate estimation.
- In another example, in Z. Zhang, et al., “TROIKA: A General Framework for Heart Rate Monitoring using Wrist-type Photoplethysmographic Signals During Intensive Physical Exercise,” IEEE Transactions on Biomedical Engineering, Vol. 62, No. 2, pp. 522-531, February, 2015, an algorithm that uses compressed sensing techniques combined with signal decomposition for de-noising and spectral tracking is proposed. The algorithm was applied to data collected from twelve different individuals carrying out a set of activities. The mean error in beats per minute (BPM) across all data sets was found to be 2.34 BPM. This technique was later improved upon in T. Schack, et al., “A New Method for Heart Rate Monitoring During Physical Exercise using Photoplethysmographic Signals,” 23rd European Signal Processing Conference (EUSIPCO), pp. 2716-2720, Aug. 31-Sep. 4, 2015 (“Schack” herein). In the improved technique, the PPG signal is cleaned up using a normalized least mean squares (NLMS) filter and is non-coherently combined. Applying the algorithm to the twelve data sets resulted in a mean error of 1.77 BPM. As is described in more detail below, this algorithm requires three Fast Fourier Transforms (FFT) which may not be suitable for embedded applications.
- The inherent assumption in accelerometry based cancellation of motion artifacts such the previously described algorithms is that the motion induced artifacts in the PPG signal are highly correlated with the acceleration signal.
FIG. 1 is a graph illustrating this assumption with data from a wrist worn device with a PPG sensor and a three-axis accelerometer sensor in which the hand motion is performed in the direction of the Y-axis of the accelerometer. This graph shows that the Y-axis acceleration signal is inversely proportional to the PPG signal when there is motion in the Y-axis. This observation motivates the use of an interference cancellation procedure such as an NLMS filter to remove the motion component from the PPG signal. - The amount of mutual information between a PPG signal and an accelerometer signal can, in a first order approximation, be captured by computing the Pearson correlation coefficient defined as
-
r=cov(ppg,acc)/σaccσppg - where cov(ppg,acc) denotes the covariance between the PPG signal and the accelerometer signal over a window of time of length T, σacc is the standard deviation of the accelerometer signal over the same window length of time T and, σppg is the standard deviation of the PPG signal over a window of time of length T.
FIG. 2 is a graph illustrating the Pearson correlation for the X-axis acceleration signal (accX), the Y-axis acceleration signal (accY) and the z-axis acceleration signal (accZ) of an accelerometer with respect to the PPG signal when the motion is predominantly along the Y-axis. In this graph, the Pearson correlation coefficient associated with the Y-axis acceleration signal is consistently larger than that associated with the X-axis and Z-axis as the motion of the hand is predominantly in the direction of the Y-axis. - In normal human movement, the motion occurs in all three axes. Thus, the motion interference cancellation should be done with respect to inputs from all three accelerometer axis signals. In Schack, the motion filtered PPG outputs with respect to accX, accY, and accZ are non-coherently combined and the heart rate is computed with respect to the resulting spectrum. The algorithm of Schack can be written as
-
S(f)=|FFT(PPGaccX)|2×|FFT(PPGaccY)|2×|FFT(PPGaccZ)|2 - Heart Rate=HeartRateTracker (S(f))
- where S(f) is the non-coherent power spectrum, HeartRateTracker computes the heart rate from the power spectrum, PPGaccX is the NLMS filter output with an input of the PPG signal with the accX signal as the reference signal, PPGaccY is the NLMS filter output with an input of the PPG signal with the accY signal as the reference signal, and PPGaccZ is the NLMS filter output with an input of the PPG signal with the accZ signal as the reference signal. A simple version of HeartRateTracker associates the heart rate with the frequency index of the largest spectrum peak in S(f).
- Embodiments of the disclosure provide for reducing motion induced artifacts in PPG signals with reduced complexity and improved performance over the prior art. As is described in more detail herein, in some embodiments, instead of requiring three FFTs with non-coherent averaging as in Schack, coherent combining of NLMS filter outputs (or other suitable adaptive filter outputs) is performed followed by a single Fourier Transform (FT). The resulting signal can then be used for heart rate estimation. Experiments have demonstrated a mean error of at least 1.27 BPM.
- In embodiments of the disclosure, the motion compensated power spectrum of the PPG signal may be computed as
-
S(f)=|FFT(a·PPGaccX +b·PPGaccY +c·PPGaccZ)|2 - where a, b, and c are weights applied to the outputs of suitable adaptive filters. Embodiments are described herein assuming the adaptive filters are NLMS filters. One of ordinary skill in the art will understand embodiments in which other adaptive filters are used, such as, for example recursive least squares (RLS) filters. NLMS filters are described in S. Haykin, “Adaptive Filter Theory”, pp. 320-324, Pearson Education, Inc., 2002, which is incorporated by reference herein. RLS filters are described in Haykin, p. 443, which is incorporated by reference herein. Any suitable technique for determining the weight values of the power spectrum computation may be used. Some examples of suitable techniques are now described.
- In some embodiments, the weights are all the same value, i.e., a=b=c=K where K is a constant. This equal weighting may be used, for example, to simplify implementation. The choice of K is implementation dependent. In some embodiments, K=1 such that the three NLMS output signals are added without weighting.
- In some embodiments, the values of the weights are based on the amount of correlation of the PPG signal with the respective acceleration signals. The weighting is based on the idea that the higher the correlation of the PPG signal with the corresponding acceleration signal, the more motion compensated the output of the corresponding NLMS filter should be. To determine the weight values, the correlation coefficients between each acceleration signal and the PPG signal are computed in a T second window and the weights are then computed as a function of the correlation coefficients. In one embodiment, the maximum correlation coefficients are used.
- For example, if there are m samples in a T second window, the correlation r1accX between the accX signal and the PPG signal at time n+m/2 can be computed as
-
- The maximum of the correlation coefficients in the T second window may be computed as
-
r accX=max(r 1accX(n) . . . r 1accX(n+m−1)). - The maximum of the correlation coefficients between the accY signal and the PPG signal, raccY, and between the accZ signal and the PPG signal, raccX, may be similarly computed. The weight values may then be determined as a function of raccX, raccY, and raccZ. For example, in one embodiment the weight values may be determined as
-
- The value of T is typically chosen to be the size of the data over which the FT is performed and is of the order of a few seconds.
- In some embodiments, the values of the weights are based on the predicted motion artifact energy, i.e., interference energy, removed by the NLMS filter over a window of T seconds. The NLMS filter output at time n when the sensor reference is the X-axis acceleration signal is given by
-
PPGaccX(n)=PPG(n)−InterfaccX(n) - where PPG(n) is the PPG signal at time n and InterfaccX(n) is the interference energy removed, which is given by
-
- where w denotes the NLMS filter coefficients for accX, and [accX(n−1), . . . , accX(n−N)] denotes the accX measurements.
- The predicted interference energy removed over a window of T seconds (Tsec) by the NLMS filter corresponding to accX may thus be determined as per
-
- The predicted interference energy corresponding to accY, IntEnaccY, and accZ, IntEnaccZ, may be similarly determined. The values of the weights may then be determined as a function of IntEnaccX, IntEnaccY, and IntEnaccY. For example, in one embodiment, the weight values may be determined as
-
- The value of T is typically chosen to be the size of the data over which the FT is performed and is of the order of a few seconds.
- In some embodiments, rather than combining all the signals output by the NLMS filters to generate the signal used for heart rate estimation, one of the signals is selected according to some suitable selection criteria. One selection criteria may be to select one of the output signals based on the maximum correlation between the PPG signal and the reference acceleration signal. More specifically, the output signal with the maximum correlation may be given a weight of 1 while the other output signals may be given a 0 weight. For example, if raccY is the maximum correlation, then a=0, b=1, and c=0 and
-
S(f)=|FFT(0*PPGaccX+1*PPGaccY+0*PPGaccZ)|2=|FFT(PPGaccY)|2. - The decision rules may be written as
-
- if raccX=max(raccX,raccY,raccZ)
- a=1;b=0;c=0
- if raccY=max(raccX,raccY,raccZ)
- a=0;b=1;c=0
- if raccZ=max(raccX,raccY,raccZ)
- a=0;b=0;c=1
- if raccX=max(raccX,raccY,raccZ)
- Another selection criteria may be to select one of the output signals having the maximum interference energy removed. More specifically, the output signal with the maximum interference energy removed may be given a weight of 1 while the other output signals may be given a 0 weight. For example, if IntEnaccY has the maximum interference energy removed, then a=0, b=1, and c=0 and
-
S(f)=|FFT(0*PPGaccX+1*PPGaccY+0*PPGaccZ)|2=|FFT(PPGaccY)|2. - The decision rules may be written as
- if IntEnaccX=max(IntEnaccX, IntEnaccY, IntEnaccZ)
-
- a=1;b=0;c=0
- if IntEnaccY=max(IntEnaccX, IntEnaccY, IntEnaccZ)
-
- a=0;b=1;c=0
- if IntEnaccZ=max(IntEnaccX, IntEnaccY, IntEnaccZ)
-
- a=0;b=0;c=1
- In some embodiments, rather than constantly filtering the PPG signal to compensate for motion artifacts, the PPG signal is filtered when the user of the device is in motion and not filtered otherwise. Filtering the PPG signal to remove motion artifacts when the user of the device is not in motion may degrade the PPG signal as motion artifacts are not present. Any suitable technique may be used to determine if the user is in motion. As shown in Table 1, one suitable technique that may be used is to compute the variance of the sum of squares of the acceleration signals over a window of T1 seconds having m1 samples and compare the variance to a predetermined threshold Thresh. If the variance exceeds the threshold, the user is determined to be in motion at time instant n; otherwise the user is determined to be stationary. Any suitable threshold value may be used. In some embodiments, the threshold value is determined empirically based on experiments performed on the acceleration signals to determine a good representation of stationary.
-
TABLE 1 Let accsq(n) = accX (n)2 + accY (n)2 + accZ (n)2 Then, if accVar(n) < Thresh −> stationary elseif accVar(n) ≧ Thresh −> motion -
FIG. 3 is a block diagram of an example embedded photoplethysmograph (PPG) heartrate monitor device 300 configured to perform motion compensation on PPG signals as described herein. The PPG heartrate monitor device 300 includes aPPG sensor component 304 configured to capture a PPG signal from the person wearing themonitor device 300 and aprocessing component 302 configured to process the PPG signal to measure the heart rate of a person wearing the monitor device. - The
PPG sensor component 304 may be any suitable PPG sensor. In general, the PPG sensor component includes one or more three light emitting diodes (LEDs) and a photodiode arranged to the wavelengths of the light from the one or more LEDs. In some embodiments, thePPG sensor component 304 includes three light emitting diodes (LEDs): a green light emitting diode (LED), a red LED, and an infrared LED. - The
processing component 302 includes a microcontroller (MCU) 308, an inertial measurement unit (IMU) 312, aBluetooth transceiver 314, an analog front end (AFE) 310, andmemory 318. TheMCU 308 may be any suitable microcontroller, such as, for example, an MSP430 device available from Texas Instruments, Inc. TheMCU 308 includes memory (one example of a computer readable medium) that can be used to store software instructions needed to perform motion compensation on the PPG signals as described herein and to perform heart rate measurement based on the motion compensated PPG signals. Embodiments of motion compensation that may be executed by theMCU 308 are described in reference toFIGS. 4, 5, and 6 .Memory 318 may be any suitable memory device or devices, such as, for example, one or more ferroelectric random access memory (FRAM) devices.Memory 318 is coupled to theMCU 308 and may be used, for example, to store data used in the execution of the heart rate monitoring and to store a heart rate history. - The
IMU 312 is coupled to theMCU 308 via an interface provided by theMCU 308 and provides data regarding the motion of a person wearing themonitor device 300 to theMCU 308. More specifically, theIMU 312 includes a three axis accelerometer to provide the X-axis acceleration signal (accX), the Y-axis acceleration (accY) signal and the z-axis acceleration (accZ) signal to theMCU 308. TheIMU 312 may be any suitable IMU device, such as, for example, an IMU device available from Invensense, Inc. - The
Bluetooth transceiver 314 is coupled to theMCU 308 and may be used, for example, to transfer the computed heart rate measurement of a person wearing themonitor device 300 to another device for display and/or further processing. For example, the heart rate measurement may be transmitted to a smart phone or other personal digital assistant, a laptop computer, a desktop computer, a medical monitoring device, etc. TheBluetooth transceiver 314 may be any suitable Bluetooth device, such as, for example, a CC25xx Bluetooth system-on-a-chip (SOC) available from Texas Instruments, Inc. - The USB (universal serial bus)
transceiver 316 is coupled to theMCU 308 and may be used, for example, to transfer the computed heart rate measurement of a person wearing themonitor device 300 to another device for display and/or further processing. For example, the heart rate measurement may be transmitted to a smart phone or other personal digital assistant, a laptop computer, a desktop computer, a medical monitoring device, etc. The associated USB port may also be used to charge a battery (not shown) in themonitor device 300. - The
AFE 310 is coupled to theMCU 308 and thePPG sensor component 304 and provides a control interface between theMCU 308 andPPG sensor component 304. TheAFE 310 includes functionality to receive a PPG signal from thePPG sensor component 304, convert the signal to a digital signal, apply analog signal conditioning such as changing the gain of the signal, and provide the digital PPG signal to theMCU 308 for heart rate measurement. TheAFE 310 also includes functionality to drive current to one or more LEDs in thePPG sensor 304 as directed by theMCU 308. The particular current to be used and the timing of when to apply the current may be controlled by theMCU 308. TheAFE 310 may be any suitable AFE device, such as, for example, the AFE4400 provided by Texas Instruments, Inc. - In some embodiments, the PPG heart
rate monitor device 300 includes a display component (not shown) coupled to theMCU 308. The display component may be used to display the heart rate of a person wearing themonitor device 300. -
FIG. 4 andFIG. 5 are block diagrams illustrating an embodiment of motion compensation of a PPG signal and heart rate tracking that may be performed in a suitably configured heart rate monitor device, such as, for example, the heartrate monitor device 300 ofFIG. 3 . Embodiments are described in reference to heartrate monitor device 300 for simplicity of explanation. One of ordinary skill in the art will understand embodiments for other hardware architectures. Further, embodiments are described assuming the various components are implemented as software executed by a processor such as theMCU 308. One of ordinary skill in the art will understand embodiments in which one of more of the components may be implemented in hardware. For example, theFT component 410 may be an FFT hardware accelerator tightly coupled with the central processing unit of theMCU 308. In another example, one or more of the components may be implemented in an application specific integrated circuit (ASIC). - Referring first to
FIG. 4 , thebandpass filter 402 receives the PPG signal from theAFE 310 is filtered 402 and outputs a filtered PPG signal. Thebandpass filter 404 receives the X-axis acceleration signal (accX), the Y-axis acceleration signal (accY) and the z-axis acceleration signal (accZ) from theIMU 312 and output filtered versions of these acceleration signals. The bandpass filters 402, 404 are filters with the same in-band and out-of-band characteristics. - The motion compensation component (MCC) 406 receives the filtered PPG signal and the filtered acceleration signals and performs motion compensation on the PPG signal based on the acceleration signals. More specifically, the
MCC 406 performs motion compensation on the PPG signal with respect to each of the acceleration signals and generates three output signals, PPGaccX, PPGaccY, and PPGaccZ. The motion compensation is described in more detail in reference toFIG. 5 . In some embodiments, theMCC 406 also outputs weighting data that may be used by thecombiner 408 to determine weights for combining the three motion filtered PPG signals, e.g., indications of the correlation of the PPG signal with each acceleration signal and/or the predicted interference energy removed from the PPG signal to generate each of the signals output by theMCC 406. - The
combiner 408 combines the three signals output by theMCC 406 as per -
a*PPGaccX+b*PPGaccY+c*PPGaccZ - to generate the final motion compensated PPG signal. Further, the
combiner 408 may determine the values of the weights as per one of the approaches previously described herein based on weighting data from theMCC 406 or the weights may all have the same value, i.e., a=b=c=K where K is a predetermined constant. In some embodiments, K=1. - The
FT component 410 applies an FFT to the final motion compensated PPG signal to convert the signal from the time domain to the frequency domain. The resulting signal is then provided to theheart rate tracker 412 to be used for heart rate estimation. - The block diagram of
FIG. 5 illustrates an embodiment of the motion compensation performed by theMCC 406. Thecorrelation indicator component 516 receives the filtered PPG signal and the filtered acceleration signals and generates weighting data for thecombiner 408 that is indicative of the correlation between the PPG signal and each of the acceleration signals. That is, thecorrelation indicator component 516 computes the correlation coefficients between each acceleration signal and the PPG signal in T second window. Computation of correlation coefficients is previously described herein. - In some embodiments, the output of the
correlation indicator component 516 is the correlation coefficients for each combination of PPG signal and acceleration signal and thecombiner 408 uses the correlation coefficients to determine the weight values used to combine the motion compensated PPG signals PPGaccX, PPGaccY, and PPGaccZ. In some embodiments, the output of thecorrelation indicator component 516 is the maximum of the correlation coefficients for each combination of PPG signal and acceleration signal and thecombiner 408 uses the maximum coefficients to determine the weight values. Computation of the maximum of the correlation coefficients is previously described herein. - The interference
energy detection component 514 receives the filtered PPG signal, the filtered acceleration signals, and the current filter coefficients from eachadaptive filter combiner 408 that is indicative of the predicted interference energy removed from the PPG signal by eachadaptive filter energy detection component 514 computes the predicted interference energy removed from the PPG signal by eachadaptive filter - The
stationary detector component 518 receives the filtered acceleration signals and generates an indication of whether or not the person wearing themonitor device 300 is stationary, i.e., a stationary indicator. Thestationary detector component 518 may use any suitable technique to determine if the person is in motion. An example of one such suitable technique is previously described herein. - The
adaptive filters adaptive filter - The
adaptive filters stationary detector component 518. If the stationary indicator indicates that the person is stationary, theadaptive filters combiner 408 performs no combining and passes one of the input PPG signals as the output PPG signal. In some embodiments, if the stationary indicator indicates that the person is stationary, the weights of the adaptive filter are not changed, i.e., the weights remain at their current values. Thus, the filter weights are only updated when the person is in motion. In some embodiments, if the stationary indicator indicates that the person is stationary, the weights of the adaptive filter are reset to their initial values. -
FIG. 5 includes both acorrelation indicator component 516 and an interferenceenergy detection component 514 for ease of description. One of ordinary skill in the art will understand embodiments in which one of these components is not present or is otherwise not used. Further, one of ordinary skill in the art will understand embodiments in which neither component is present, e.g., in embodiments in which the combiner uses a weighting scheme in which the weight values are constants. -
FIG. 6 is a flow diagram of a method for motion compensation of a PPG signal and heart rate tracking that may be performed in a suitably configured heart rate monitor device, such as, for example, the heartrate monitor device 300 ofFIG. 3 . Embodiments are described in reference to heartrate monitor device 300 for simplicity of explanation. One of ordinary skill in the art will understand embodiments for other hardware architectures. The method may be performed continuously while the heart rate monitor device is turned on. - A PPG signal is received 600 by the
MCU 308 from thePPG sensor component 304 via theAFE 310 and is filtered using a bandpass filter. Accelerometer signals, i.e., an X-axis acceleration signal, a Y-axis acceleration signal, and a Z-axis acceleration signal are also received 602 by theMCU 308 from theIMU 312 and are filtered using a bandpass filter. Motion compensation is then performed 604 on the filtered PPG signal using the filtered acceleration signals as reference signals to generate three motion compensated PPG signals. More specifically, the filtered PPG signal is adaptively filtered using each of the X-axis acceleration signal, the Y-axis acceleration signal, and the Z-axis acceleration signal as a reference signal to generate three motion compensated PPG signals. The adaptive filter may be, for example, an NLMS filter or an RLS filter. - The three motion compensated PPG signals are then coherently combined 606 as per
-
a*PPGaccX+b*PPGaccY+c*PPGaccZ - to generate a final motion compensated PPG signal. The values of the weights may be determined as per one of the approaches previously described herein or the weights may all have the same value, i.e., a=b=c=K where K is a predetermined constant. In some embodiments, K=1.
- A single FT is then applied 608 to the final motion compensated PPG signal to generate a frequency domain PPG signal. The frequency domain signal is then used to estimate 610 the heart rate.
- In some embodiments, a determination is also made as to whether or not the user of the heart rate monitor device is stationary or in motion. Any suitable technique may be used to make this determination. One such suitable technique is previously described herein. If the user is stationary, the PPG signal is not adaptively filtered. In some embodiments, the weights of the adaptive filters remain unchanged if the user is stationary. In some embodiments, the weights of the adaptive filters are reset to initial values if the user is stationary.
- While the invention has been described with respect to a limited number of embodiments, those skilled in the art, having benefit of this disclosure, will appreciate that other embodiments can be devised which do not depart from the scope of the invention as disclosed herein.
- For example, one of ordinary skill in the art will understand that embodiments of may be used in PPG heart rate monitor devices housed in form factors suitable for monitoring heart rate in various locations on the body of a person such as, for example, the forehead, an ear lobe, a finger, an upper arm, a wrist, or a chest.
- In another example, one of ordinary skill in the art will understand embodiments in which the filtered PPG signal and the filtered acceleration signals are downsampled prior to further processing.
- In another example, one of ordinary skill in the art will understand embodiments in which the stationary detector component is not present and the IMU determines whether the person is in motion and generates the value of the stationary indicator.
- In another example, while some embodiments have been described herein assuming an FFT is used, one of ordinary skill in the art will understand embodiments in which other suitable Fourier Transforms are used.
- The methods described herein may be implemented in hardware, software, firmware, or any combination thereof. If completely or partially implemented in software, the software may be executed in one or more processors, such as a microcontroller, microprocessor, application specific integrated circuit (ASIC), field programmable gate array (FPGA), or digital signal processor (DSP). The software instructions may be initially stored in a computer-readable medium and loaded and executed in the processor. In some cases, the software instructions may also be sold in a computer program product, which includes the computer-readable medium and packaging materials for the computer-readable medium. In some cases, the software instructions may be distributed via removable computer readable media, via a transmission path from computer readable media on another digital system, etc. Examples of computer-readable media include non-writable storage media such as read-only memory devices, writable storage media such as disks, flash memory, memory, or a combination thereof.
- Although method steps may be presented and described herein in a sequential fashion, one or more of the steps shown in the figures and described herein may be performed concurrently, may be combined, and/or may be performed in a different order than the order shown in the figures and/or described herein. Accordingly, embodiments should not be considered limited to the specific ordering of steps shown in the figures and/or described herein.
- It is therefore contemplated that the appended claims will cover any such modifications of the embodiments as fall within the true scope of the invention.
Claims (20)
1. A method for heart rate measurement in a photoplethysmograph (PPG) heart rate monitor device, the method comprising:
receiving a PPG signal from a PPG sensor of the heart rate monitor device;
receiving an X-axis acceleration signal, a Y-axis acceleration signal, and a Z-axis acceleration signal from an accelerometer of the heart rate monitor device;
performing motion compensation on the PPG signal, wherein a first motion compensated PPG signal is generated with reference to the X-axis acceleration signal, a second motion compensated PPG signal is generated with reference to the Y-axis acceleration signal, and a third motion compensated PPG signal is generated with reference to the Z-axis acceleration signal;
combining the first motion compensated PPG signal, the second motion compensated PPG signal, and the third motion compensated PPG signal to generate a fourth motion compensated PPG signal, wherein a first weight is applied to the first motion compensated PPG, a second weight is applied to the second motion compensated PPG signal, and a third weight is applied to the third motion compensated PPG signal;
performing a single Fourier Transform (FT) on the fourth motion compensated PPG signal to generate a frequency domain PPG signal; and
estimating a heart rate based on the frequency domain PPG signal.
2. The method of claim 1 , wherein the first motion compensated PPG signal is generated by applying a first adaptive filter to the PPG signal with the X-axis acceleration signal as a reference signal, the second motion compensated PPG signal is generated by applying a second adaptive filter to the PPG signal with the Y-axis acceleration signal as a reference signal, and the third motion compensated PPG signal is generated by applying a third adaptive filter to the PPG signal with the Z-axis acceleration signal as a reference signal.
3. The method of claim 2 , wherein the first adaptive filter, the second adaptive filter, and the third adaptive filter are a same filter type, the filter type selected from one of a normalized least mean squares (NLMS) filter and a recursive least squares (RLS) filter.
4. The method of claim 2 , wherein weights of the first adaptive filter, the second adaptive filter, and the third adaptive filter are updated only when a user of the heart rate monitor device is in motion.
5. The method of claim 2 , wherein weights of the first adaptive filter, the second adaptive filter, and the third adaptive filter are reset to initial values when a user of the heart rate monitor device is stationary.
6. The method of claim 1 , further comprising determining whether or not a user of the heart rate monitor device is stationary.
7. The method of claim 6 , wherein determining whether or not a user of the heart rate monitor device is stationary comprises:
computing a variance of the sum of the squares of the X-axis acceleration signal, the Y-axis acceleration signal, and the Z-acceleration signal over a window of time; and
comparing the variance to a threshold.
8. The method of claim 1 , wherein the first weight, the second weight, and the third weight are equal.
9. The method of claim 1 , wherein a value of the first weight, a value of the second weight, and a value of the third weight are based on an amount of correlation between the PPG signal and the X-axis acceleration signal, an amount of correlation between the PPG signal and Y-axis acceleration signal, and an amount of correlation between the PPG signal and the Z-axis acceleration signal.
10. The method of claim 2 , wherein a value of the first weight, a value of the second weight, and a value of the third weight are based on predicted interference energy removed from the PPG signal by the first adaptive filter, predicted interference energy removed from the PPG signal by the second adaptive filter, and predicted interference energy removed from the PPG signal by the third adaptive filter.
11. A photoplethysmograph (PPG) heart rate monitor device comprising:
a PPG sensor component configured to generate an analog PPG signal;
an accelerometer configured to generate an X-axis acceleration signal, a Y-axis acceleration signal, and a Z-axis acceleration signal;
a motion compensation component coupled to the PPG sensor to receive the PPG signal and coupled to the accelerometer to receive the X-axis acceleration signal, the Y-axis acceleration signal, and the Z-axis acceleration signal, wherein the motion compensation component is configured to generate a first motion compensated PPG signal with reference to the X-axis acceleration signal, a second motion compensated PPG signal with reference to the Y-axis acceleration signal, and a third motion compensated PPG signal with reference to the Z-axis acceleration signal;
a combiner configured to combine the first motion compensated PPG signal, the second motion compensated PPG signal, and the third motion compensated PPG signal to generate a fourth motion compensated PPG signal, wherein a first weight is applied to the first motion compensated PPG, a second weight is applied to the second motion compensated PPG signal, and a third weight is applied to the third motion compensated PPG signal;
a Fourier Transform (FT) component configured to perform a single FT on the fourth motion compensated PPG signal to generate a frequency domain PPG signal; and
a heart rate tracker component configured to estimate a heart based on the frequency domain PPG signal.
12. The PPG heart rate monitor device of claim 11 , wherein the motion compensation component comprises a first adaptive filter configured to generate the first motion compensated PPG signal using the X-axis acceleration signal as a reference signal, a second adaptive filter configured to generate the second motion compensated PPG signal using the Y-axis acceleration signal as a reference signal, and a third adaptive filter configured to generate the third motion compensated PPG signal using the Z-axis acceleration signal as a reference signal.
13. The PPG heart rate monitor device of claim 12 , wherein the first adaptive filter, the second adaptive filter, and the third adaptive filter are a same filter type, the filter type selected from one of a normalized least mean squares (NLMS) filter and a recursive least squares (RLS) filter.
14. The PPG heart rate monitor device of claim 12 , wherein weights of the first adaptive filter, the second adaptive filter, and the third adaptive filter are updated only when a user of the heart rate monitor device is in motion.
15. The PPG heart rate monitor device of claim 12 , wherein weights of the first adaptive filter, the second adaptive filter, and the third adaptive filter are reset to initial values when a user of the heart rate monitor device is stationary.
16. The PPG heart rate monitor device of claim 11 , further comprising a stationary detector component configured to determine whether or not a user of the heart rate monitor device is stationary.
17. The PPG heart rate monitor device of claim 16 , wherein the stationary detector component is coupled to the accelerometer to receive the X-axis acceleration signal, the Y-axis acceleration signal, and the Z-axis acceleration signal and is configured to determine whether or not a user of the heart rate monitor device is stationary by computing a variance of the sum of the squares of the X-axis acceleration signal, the Y-axis acceleration signal, and the Z-acceleration signal over a window of time and comparing the variance to a threshold.
18. The PPG heart rate monitor device of claim 11 , wherein the first weight, the second weight, and the third weight are equal.
19. The PPG heart rate monitor device of claim 11 , wherein a value of the first weight, a value of the second weight, and a value of the third weight are based on an amount of correlation between the PPG signal and the X-axis acceleration signal, an amount of correlation between the PPG signal and Y-axis acceleration signal, and an amount of correlation between the PPG signal and the Z-axis acceleration signal.
20. The PPG heart rate monitor device of claim 12 , wherein a value of the first weight, a value of the second weight, and a value of the third weight are based on predicted interference energy removed from the PPG signal by the first adaptive filter, predicted interference energy removed from the PPG signal by the second adaptive filter, and predicted interference energy removed from the PPG signal by the third adaptive filter.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/970,434 US20170164847A1 (en) | 2015-12-15 | 2015-12-15 | Reducing Motion Induced Artifacts in Photoplethysmography (PPG) Signals |
CN201611138819.XA CN106880351B (en) | 2015-12-15 | 2016-12-12 | Reducing motion-induced artifacts in photoplethysmography (PPG) signals |
US16/271,246 US10893815B2 (en) | 2015-12-15 | 2019-02-08 | Heart rate measurement and a photoplethysmograph (PPG) heart rate monitor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/970,434 US20170164847A1 (en) | 2015-12-15 | 2015-12-15 | Reducing Motion Induced Artifacts in Photoplethysmography (PPG) Signals |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/271,246 Division US10893815B2 (en) | 2015-12-15 | 2019-02-08 | Heart rate measurement and a photoplethysmograph (PPG) heart rate monitor |
Publications (1)
Publication Number | Publication Date |
---|---|
US20170164847A1 true US20170164847A1 (en) | 2017-06-15 |
Family
ID=59018680
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/970,434 Abandoned US20170164847A1 (en) | 2015-12-15 | 2015-12-15 | Reducing Motion Induced Artifacts in Photoplethysmography (PPG) Signals |
US16/271,246 Active 2036-01-08 US10893815B2 (en) | 2015-12-15 | 2019-02-08 | Heart rate measurement and a photoplethysmograph (PPG) heart rate monitor |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/271,246 Active 2036-01-08 US10893815B2 (en) | 2015-12-15 | 2019-02-08 | Heart rate measurement and a photoplethysmograph (PPG) heart rate monitor |
Country Status (2)
Country | Link |
---|---|
US (2) | US20170164847A1 (en) |
CN (1) | CN106880351B (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10667758B2 (en) * | 2016-02-02 | 2020-06-02 | Fujitsu Limited | Sensor information processing apparatus |
CN111616695A (en) * | 2020-06-29 | 2020-09-04 | 歌尔科技有限公司 | Heart rate acquisition method, device, system and medium |
CN113261932A (en) * | 2021-06-28 | 2021-08-17 | 山东大学 | Heart rate measurement method and device based on PPG signal and one-dimensional convolutional neural network |
US11311242B2 (en) * | 2016-05-20 | 2022-04-26 | Sony Corporation | Biological information processing apparatus, biological information processing method, and information processing apparatus |
WO2022139817A1 (en) * | 2020-12-22 | 2022-06-30 | Google Llc | System and method for ppg using a wideband accelerometer with a wearable device |
US11452483B2 (en) * | 2016-09-06 | 2022-09-27 | Soongsil University Research Consortium Techno-Park | Apparatus for eliminating motion artifacts by using PPG signal and method thereof |
US11529096B2 (en) * | 2017-09-29 | 2022-12-20 | Kyocera Corporation | Sleep assessment system, massage system, control method, and electronic device |
EP4140392A1 (en) * | 2021-08-23 | 2023-03-01 | Nokia Technologies Oy | Noise removal in physiological signals |
WO2023159223A1 (en) * | 2022-02-18 | 2023-08-24 | Rajant Health Incorporated | Parallel adaptive motion artifact removal |
CN117084644A (en) * | 2023-10-10 | 2023-11-21 | 荣耀终端有限公司 | Method and device for detecting physiological characteristic index of user and wearable device |
US11850026B2 (en) | 2020-06-24 | 2023-12-26 | The Governing Council Of The University Of Toronto | Remote portable vital signs monitoring |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6723132B2 (en) * | 2016-09-29 | 2020-07-15 | ルネサスエレクトロニクス株式会社 | Pulse measuring device, light intensity control method, and program |
CN108056769B (en) * | 2017-11-14 | 2020-10-16 | 深圳市大耳马科技有限公司 | Vital sign signal analysis processing method and device and vital sign monitoring equipment |
CN108903929B (en) * | 2018-03-30 | 2021-02-26 | 广东思派康电子科技有限公司 | Heart rate detection correction method, device, storage medium and system |
CN110547785B (en) * | 2018-06-04 | 2022-08-05 | 原相科技股份有限公司 | Microcirculation detection system and detection method |
CN109044316B (en) * | 2018-07-11 | 2021-02-02 | 歌尔科技有限公司 | Pure dynamic heart rate signal calculation method and device and intelligent wearable equipment |
CN109222948B (en) * | 2018-09-17 | 2021-07-13 | 歌尔科技有限公司 | Method, device, electronic equipment and storage medium for eliminating motion interference noise |
CN110169764A (en) * | 2019-05-06 | 2019-08-27 | 上海理工大学 | A kind of LMS adaptive-filtering PPG signal heart rate extracting method |
CN113749627A (en) * | 2020-08-20 | 2021-12-07 | 泰科兴业科技(深圳)有限公司 | Novel personnel intelligent social management system and method |
CN111904406A (en) * | 2020-08-25 | 2020-11-10 | 上海交通大学 | Physiological signal motion artifact suppression device and method |
CN112370036A (en) * | 2020-10-20 | 2021-02-19 | 复旦大学 | PPG heart rate extraction device and method based on cascade RLS adaptive filtering |
CN114557685B (en) * | 2020-11-27 | 2023-11-14 | 上海交通大学 | Non-contact type exercise robust heart rate measurement method and measurement device |
CN115886744A (en) * | 2022-11-09 | 2023-04-04 | 广东粤港澳大湾区黄埔材料研究院 | Motion pulse wave denoising method and device |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080152161A1 (en) * | 2006-12-21 | 2008-06-26 | Samsung Electronics Co., Ltd. | System and method for determining application of adaptive filter |
US20120203080A1 (en) * | 2009-11-17 | 2012-08-09 | Min-Joon Kim | Photoplethysmography apparatus |
US20150351646A1 (en) * | 2014-06-09 | 2015-12-10 | Stmicroelectronics S.R.L. | Method for the estimation of the heart-rate and corresponding system |
US20160051157A1 (en) * | 2014-08-22 | 2016-02-25 | Apple Inc. | Frequency domain projection algorithm |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4326486B2 (en) * | 2005-02-28 | 2009-09-09 | 三洋電機株式会社 | Period calculation device and program |
JP4657300B2 (en) * | 2005-09-27 | 2011-03-23 | シチズンホールディングス株式会社 | Heart rate monitor and heart rate detection method |
US8956294B2 (en) * | 2009-05-20 | 2015-02-17 | Sotera Wireless, Inc. | Body-worn system for continuously monitoring a patients BP, HR, SpO2, RR, temperature, and motion; also describes specific monitors for apnea, ASY, VTAC, VFIB, and ‘bed sore’ index |
JP5605269B2 (en) * | 2011-02-28 | 2014-10-15 | セイコーエプソン株式会社 | Beat detector |
BR112015001898A2 (en) * | 2012-08-01 | 2017-07-04 | Koninklijke Philips Nv | system for identifying movement artifacts; and method for identifying movement artifacts |
MX366945B (en) * | 2012-12-14 | 2019-07-31 | Koninklijke Philips Nv | A system and method to detect significant arrhythmic events through a photoplethysmogram (ppg) and accelerometer. |
CN104902809B (en) * | 2012-12-31 | 2020-08-18 | 皇家飞利浦有限公司 | Method and apparatus for reducing motion artifacts in ECG signals |
JP5979604B2 (en) * | 2013-02-06 | 2016-08-24 | カシオ計算機株式会社 | Biological information detection apparatus, biological information detection method, and biological information detection program |
CN105101870B (en) * | 2013-03-29 | 2019-01-22 | 皇家飞利浦有限公司 | Device and method for the removal of ECG motion artifacts |
US9107644B2 (en) * | 2013-07-05 | 2015-08-18 | James Tyler Frix | Continuous transdermal monitoring system and method |
WO2015112266A2 (en) * | 2013-12-06 | 2015-07-30 | The Regents Of The University Of Michigan | Epsilon-tube filter for blunt noise removal |
EP3148404A1 (en) * | 2014-05-28 | 2017-04-05 | Koninklijke Philips N.V. | Motion artifact reduction using multi-channel ppg signals |
US10265024B2 (en) * | 2014-07-26 | 2019-04-23 | Salutron, Inc. | Sensor system for heart rate measurement per axis of shared orientation |
CN104161505A (en) * | 2014-08-13 | 2014-11-26 | 北京邮电大学 | Motion noise interference eliminating method suitable for wearable heart rate monitoring device |
US9901306B2 (en) * | 2015-04-29 | 2018-02-27 | Analog Devices, Inc. | Tracking mechanism for heart rate measurements |
US10758185B2 (en) * | 2015-11-25 | 2020-09-01 | Texas Instruments Incorporated | Heart rate estimation apparatus using digital automatic gain control |
KR101831064B1 (en) * | 2016-09-06 | 2018-02-22 | 숭실대학교산학협력단 | Apparatus for motion artifact removal using ppg signal and method thereof |
-
2015
- 2015-12-15 US US14/970,434 patent/US20170164847A1/en not_active Abandoned
-
2016
- 2016-12-12 CN CN201611138819.XA patent/CN106880351B/en active Active
-
2019
- 2019-02-08 US US16/271,246 patent/US10893815B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080152161A1 (en) * | 2006-12-21 | 2008-06-26 | Samsung Electronics Co., Ltd. | System and method for determining application of adaptive filter |
US20120203080A1 (en) * | 2009-11-17 | 2012-08-09 | Min-Joon Kim | Photoplethysmography apparatus |
US20150351646A1 (en) * | 2014-06-09 | 2015-12-10 | Stmicroelectronics S.R.L. | Method for the estimation of the heart-rate and corresponding system |
US20160051157A1 (en) * | 2014-08-22 | 2016-02-25 | Apple Inc. | Frequency domain projection algorithm |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10667758B2 (en) * | 2016-02-02 | 2020-06-02 | Fujitsu Limited | Sensor information processing apparatus |
US11311242B2 (en) * | 2016-05-20 | 2022-04-26 | Sony Corporation | Biological information processing apparatus, biological information processing method, and information processing apparatus |
US11452483B2 (en) * | 2016-09-06 | 2022-09-27 | Soongsil University Research Consortium Techno-Park | Apparatus for eliminating motion artifacts by using PPG signal and method thereof |
US11529096B2 (en) * | 2017-09-29 | 2022-12-20 | Kyocera Corporation | Sleep assessment system, massage system, control method, and electronic device |
US11850026B2 (en) | 2020-06-24 | 2023-12-26 | The Governing Council Of The University Of Toronto | Remote portable vital signs monitoring |
CN111616695A (en) * | 2020-06-29 | 2020-09-04 | 歌尔科技有限公司 | Heart rate acquisition method, device, system and medium |
WO2022139817A1 (en) * | 2020-12-22 | 2022-06-30 | Google Llc | System and method for ppg using a wideband accelerometer with a wearable device |
CN113261932A (en) * | 2021-06-28 | 2021-08-17 | 山东大学 | Heart rate measurement method and device based on PPG signal and one-dimensional convolutional neural network |
EP4140392A1 (en) * | 2021-08-23 | 2023-03-01 | Nokia Technologies Oy | Noise removal in physiological signals |
WO2023159223A1 (en) * | 2022-02-18 | 2023-08-24 | Rajant Health Incorporated | Parallel adaptive motion artifact removal |
CN117084644A (en) * | 2023-10-10 | 2023-11-21 | 荣耀终端有限公司 | Method and device for detecting physiological characteristic index of user and wearable device |
Also Published As
Publication number | Publication date |
---|---|
CN106880351A (en) | 2017-06-23 |
US20190167129A1 (en) | 2019-06-06 |
CN106880351B (en) | 2021-05-25 |
US10893815B2 (en) | 2021-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10893815B2 (en) | Heart rate measurement and a photoplethysmograph (PPG) heart rate monitor | |
CN107949321B (en) | Temporal interference removal and improved heart rate measurement tracking mechanism | |
US9936886B2 (en) | Method for the estimation of the heart-rate and corresponding system | |
EP3277162B1 (en) | Wearable pulse sensing device signal quality estimation | |
US11672437B2 (en) | Method and device for tissue monitoring and heart rate detection | |
CN109414203B (en) | Online heart rate estimation based on optical measurements | |
Zhang | Heart rate monitoring from wrist-type photoplethysmographic (PPG) signals during intensive physical exercise | |
US20090105556A1 (en) | Measurement of physiological signals | |
Chung et al. | Deep learning for heart rate estimation from reflectance photoplethysmography with acceleration power spectrum and acceleration intensity | |
Arunkumar et al. | Heart rate estimation from photoplethysmography signal for wearable health monitoring devices | |
US20190192079A1 (en) | System and a method for motion artifact reduction in a ppg signal | |
US20180085069A1 (en) | Method and Device for Processing Bio-Signals | |
Islam et al. | PREHEAT: Precision heart rate monitoring from intense motion artifact corrupted PPG signals using constrained RLS and wavelets | |
Periyasamy et al. | Review on heart-rate estimation from photoplethysmography and accelerometer signals during physical exercise | |
Ayesha et al. | Heart rate monitoring using PPG with smartphone camera | |
US20190175037A1 (en) | Motion detection and cancellation using ambient light | |
US11701027B2 (en) | Optical respiration rate detection device | |
WO2020160058A1 (en) | Systems and methods for computationally efficient non-invasive blood quality measurement | |
KR20190011026A (en) | Apparatus and method of blood pressure measurement | |
Arunkumar et al. | Improved heart rate estimation from photoplethysmography during physical exercise using combination of NLMS and RLS adaptive filters | |
US20140371554A1 (en) | Method for improved oxygen saturation estimation in the presence of noise | |
EP3753484B1 (en) | Apparatus and method for calibrating bio-information estimation model, and apparatus for estimating bio-information | |
WO2020012807A1 (en) | Biological information processing device and biological information processing method | |
US20230263474A1 (en) | Parallel adaptive motion artifact removal | |
Bejgam et al. | Signal Processing Method for Enhancement of Pulse Oximeter Signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: TEXAS INSTRUMENTS INCORPORATED, TEXAS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:PANDE, TARKESH;MAGEE, DAVID PATRICK;REEL/FRAME:037299/0368 Effective date: 20151215 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |