A Novel Time-Varying Spectral Filtering Algorithm for Reconstruction of Motion Artifact Corrupted Heart Rate Signals During Intense Physical Activities Using a Wearable Photoplethysmogram Sensor

Accurate estimation of heart rates from photoplethysmogram (PPG) signals during intense physical activity is a very challenging problem. This is because strenuous and high intensity exercise can result in severe motion artifacts in PPG signals, making accurate heart rate (HR) estimation difficult. In this study we investigated a novel technique to accurately reconstruct motion-corrupted PPG signals and HR based on time-varying spectral analysis. The algorithm is called Spectral filter algorithm for Motion Artifacts and heart rate reconstruction (SpaMA). The idea is to calculate the power spectral density of both PPG and accelerometer signals for each time shift of a windowed data segment. By comparing time-varying spectra of PPG and accelerometer data, those frequency peaks resulting from motion artifacts can be distinguished from the PPG spectrum. The SpaMA approach was applied to three different datasets and four types of activities: (1) training datasets from the 2015 IEEE Signal Process. Cup Database recorded from 12 subjects while performing treadmill exercise from 1 km/h to 15 km/h; (2) test datasets from the 2015 IEEE Signal Process. Cup Database recorded from 11 subjects while performing forearm and upper arm exercise. (3) Chon Lab dataset including 10 min recordings from 10 subjects during treadmill exercise. The ECG signals from all three datasets provided the reference HRs which were used to determine the accuracy of our SpaMA algorithm. The performance of the SpaMA approach was calculated by computing the mean absolute error between the estimated HR from the PPG and the reference HR from the ECG. The average estimation errors using our method on the first, second and third datasets are 0.89, 1.93 and 1.38 beats/min respectively, while the overall error on all 33 subjects is 1.86 beats/min and the performance on only treadmill experiment datasets (22 subjects) is 1.11 beats/min. Moreover, it was found that dynamics of heart rate variability can be accurately captured using the algorithm where the mean Pearson’s correlation coefficient between the power spectral densities of the reference and the reconstructed heart rate time series was found to be 0.98. These results show that the SpaMA method has a potential for PPG-based HR monitoring in wearable devices for fitness tracking and health monitoring during intense physical activities.

[1]  Paulo S. R. Diniz,et al.  Adaptive Filtering: Algorithms and Practical Implementation , 1997 .

[2]  Giuseppe Baselli,et al.  The influence of exercise intensity on the power spectrum of heart rate variability , 2004, European Journal of Applied Physiology and Occupational Physiology.

[3]  Boreom Lee,et al.  Improved elimination of motion artifacts from a photoplethysmographic signal using a Kalman smoother with simultaneous accelerometry , 2010, Physiological measurement.

[4]  Pierre Comon,et al.  Independent component analysis, A new concept? , 1994, Signal Process..

[5]  Heng Tao Shen,et al.  Principal Component Analysis , 2009, Encyclopedia of Biometrics.

[6]  Haruki Kawanaka,et al.  Estimating heart rate using wrist-type Photoplethysmography and acceleration sensor while running , 2012, 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[7]  Jianchu Yao,et al.  A Short Study to Assess the Potential of Independent Component Analysis for Motion Artifact Separation in Wearable Pulse Oximeter Signals , 2005, 2005 IEEE Engineering in Medicine and Biology 27th Annual Conference.

[8]  Kirk H. Shelley,et al.  Estimation of Respiratory Rate From ECG, Photoplethysmogram, and Piezoelectric Pulse Transducer Signals: A Comparative Study of Time–Frequency Methods , 2010, IEEE Transactions on Biomedical Engineering.

[9]  C. Collet,et al.  Autonomic nervous system responses correlate with mental rehearsal in volleyball training , 1998, European Journal of Applied Physiology and Occupational Physiology.

[10]  Zhilin Zhang,et al.  Photoplethysmography-Based Heart Rate Monitoring in Physical Activities via Joint Sparse Spectrum Reconstruction , 2015, IEEE Transactions on Biomedical Engineering.

[11]  Marzo Edir Da Silva-Grigoletto,et al.  Heart rate variability during high-intensity exercise , 2013, J. Syst. Sci. Complex..

[12]  Steve Warren,et al.  Two-Stage Approach for Detection and Reduction of Motion Artifacts in Photoplethysmographic Data , 2010, IEEE Transactions on Biomedical Engineering.

[13]  B. Hofmann-Wellenhof,et al.  Introduction to spectral analysis , 1986 .

[14]  F. Iellamo Neural mechanisms of cardiovascular regulation during exercise , 2001, Autonomic Neuroscience.

[15]  Willis J. Tompkins,et al.  A Real-Time QRS Detection Algorithm , 1985, IEEE Transactions on Biomedical Engineering.

[16]  Yang Hao,et al.  Detecting Vital Signs with Wearable Wireless Sensors , 2010, Sensors.

[17]  J. Lee,et al.  Respiratory Rate Extraction Via an Autoregressive Model Using the Optimal Parameter Search Criterion , 2010, Annals of Biomedical Engineering.

[18]  S. Seyedtabaii,et al.  Kalman Filter Based Adaptive Reduction of Motion Artifact from Photoplethysmographic Signal , 2008 .

[19]  A. Jeukendrup,et al.  Heart Rate Monitoring , 2003, Sports medicine.

[20]  Zhilin Zhang,et al.  TROIKA: A General Framework for Heart Rate Monitoring Using Wrist-Type Photoplethysmographic Signals During Intensive Physical Exercise , 2014, IEEE Transactions on Biomedical Engineering.

[21]  [Autonomic nervous system adaptation during physical exercise]. , 2001, Italian heart journal. Supplement : official journal of the Italian Federation of Cardiology.

[22]  J. Elsner,et al.  Singular Spectrum Analysis: A New Tool in Time Series Analysis , 1996 .

[23]  Vasilis Z. Marmarelis,et al.  Nonlinear analysis of renal autoregulation under broadband forcing conditions , 1993, Annals of Biomedical Engineering.

[24]  Jo Woon Chong,et al.  Arrhythmia Discrimination Using a Smart Phone , 2013, IEEE Journal of Biomedical and Health Informatics.

[25]  Andriy Temko,et al.  Estimation of heart rate from photoplethysmography during physical exercise using Wiener filtering and the phase vocoder , 2015, 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[26]  N.V. Thakor,et al.  Applications of adaptive filtering to ECG analysis: noise cancellation and arrhythmia detection , 1991, IEEE Transactions on Biomedical Engineering.

[27]  K. Tremper,et al.  Pulse oximetry. , 1989, Anesthesiology.

[28]  L. Sinoway,et al.  Cardiovascular response to exercise. , 1999, The American journal of physiology.

[29]  Sun K. Yoo,et al.  Motion artifact reduction in photoplethysmography using independent component analysis , 2006, IEEE Transactions on Biomedical Engineering.

[30]  Domenico Prattichizzo,et al.  Application of Kalman Filter to Remove TMS-Induced Artifacts from EEG Recordings , 2008, IEEE Transactions on Control Systems Technology.

[31]  B. Thompson Canonical Correlation Analysis: Uses and Interpretation , 1984 .

[32]  A. Aubert,et al.  Heart Rate Variability in Athletes , 2003, Sports medicine.

[33]  Ki H. Chon,et al.  Estimation of Respiratory Rate From Photoplethysmogram Data Using Time–Frequency Spectral Estimation , 2009, IEEE Transactions on Biomedical Engineering.

[34]  Jo Woon Chong,et al.  Photoplethysmograph Signal Reconstruction based on a Novel Motion Artifact Detection-Reduction Approach. Part II: Motion and Noise Artifact Removal , 2014, Annals of Biomedical Engineering.

[35]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .

[36]  Toshiyo Tamura,et al.  Relationship Between Measurement Site and Motion Artifacts in Wearable Reflected Photoplethysmography , 2011, Journal of Medical Systems.

[37]  Ki H. Chon,et al.  Time-Varying Autoregressive Model-Based Multiple Modes Particle Filtering Algorithm for Respiratory Rate Extraction From Pulse Oximeter , 2011, IEEE Transactions on Biomedical Engineering.