SCIENCE CHINA Information Sciences, Volume 61 , Issue 2 : 022312(2018) https://doi.org/10.1007/s11432-017-9168-2

A novel approach framework based on statistics for reconstruction and heartrate estimation from PPG with heavy motion artifacts

More info
  • ReceivedMay 2, 2017
  • AcceptedJun 27, 2017
  • PublishedNov 13, 2017



This work was supported by National Natural Science Foundation of China (Grant Nos. 61634006, 61372060, 61335010, 61474107, 81300803), National Key Technologies R&D Program (Grant Nos. 2016YFB0401303, 2016YFB0402405), Basic Research Project of Shanghai Science and Technology Commission (Grant No. 16JC1400101), and Key Research Program of Frontier Science, Chinese Academy of Sciences (Grant No. QYZDY-SSW-JSC004).


[1] Selvaraj N, Jaryal A, Santhosh J, et al. Assessment of heart rate variability derived from finger-tip photoplethysmography as compared to electrocardiography. J Med Eng Technol, 2008, 32: 479--484. Google Scholar

[2] Petersen C L, Chen T P, Ansermino J M, et al. Design and evaluation of a low-cost smartphone pulse oximeter. Sensors, 2013, 13: 16882--16893. Google Scholar

[3] Lee J, Jung W. Design the filter to reject motion artifact of pulse oximetry. Elsevier Comput Sci, 2003, 26: 241--249. Google Scholar

[4] Zhang Z, Pi Z, Liu B. TROIKA: a general framework for heart rate monitoring using wrist-type photoplethysmographic signals during intensive physical exercise. IEEE Trans Biomed Eng, 2015, 62: 522--531. Google Scholar

[5] Murthy N K L, Madhusudana P C, Suresha P, et al. Multiple spectral peak tracking for heart rate monitoring from photoplethysmography signal during intensive physical exercise. IEEE Signal Process Lett, 2015, 22: 2391--2395. Google Scholar

[6] Komaty A, Boudraa A O, Augier B, et al. EMD-based filtering using similarity measure between probability density functions of IMFs. IEEE Trans Instrum Meas, 2014, 63: 27--34. Google Scholar

[7] Yang G, Liu Y, Wang Y, et al. EMD interval thresholding denoising based on similarity measure to select relevant modes. Signal Process, 2015, 109: 95--109. Google Scholar

[8] Sung P, Syed Z, Guttag J. Quantifying morphology changes in time series data with skew. In: Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing, ICASSP 2009, Taipei, 2009. 477--480. Google Scholar

[9] Wu Z, Huang N E. Ensemble empirical mode decomposition: a noise-assisted data analysis method. Adv Adapt Data Anal, 2009, 1: 1--41. Google Scholar

[10] Motin M A, Karmakar C, Palaniswami M. Ensemble Empirical Mode Decomposition with Principal Component Analysis: A Novel Approach for Extracting Respiratory Rate and Heart Rate from Photoplethysmographic Signal.. IEEE J Biomed Health Inform, 2017, : 1-1 CrossRef PubMed Google Scholar

[11] Pang B, Liu M, Zhang X, et al. Advanced EMD method using variance characterization for PPG with motion artifact. In: Proceedings of 2016 IEEE Biomedical Circuits and Systems Conference (BioCAS), Shanghai, 2016. 196--199. Google Scholar

[12] Huang N E, Shen Z, Long S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc Math Phys Eng Sci, 1998, 454: 903--995. Google Scholar

[13] Li P, Liu M, Zhang X, et al. A low-complexity ECG processing algorithm based on the Haar wavelet transform for portable health-care devices. Sci China Inf Sci, 2014, 57: 122303. Google Scholar

[14] Gao Y, Ge G, Sheng Z, et al. Analysis and solution to the mode mixing phenomenon in EMD. Congress Image Signal Process, 2008, 5: 223--227. Google Scholar

[15] Dimitrakopoulos R, Mustapha H, Gloaguen E. High-order statistics of spatial random fields: exploring spatial cumulants for modeling complex non-Gaussian and non-linear phenomena. Math Geosci, 2010, 42: 65--99. Google Scholar

[16] Donoho D L. De-noising by soft-thresholding. IEEE Trans Inf Theory, 1995, 41: 613--627. Google Scholar

[17] Li X, Zhang X, Li P, et al. An 8.12 μW wavelet denoising chip for PPG detection and portable heart rate monitoring in 0.18 $\mu$m CMOS. J Semicond, 2016, 37: 055006. Google Scholar

[18] Zhang Z. Photoplethysmography-based heart rate monitoring in physical activities via joint sparse spectrum reconstruction. IEEE Trans Biomed Eng, 2015, 62: 1902--1910. Google Scholar

  • Figure 1

    (Color online) (a) Clean PPG siganl and PPG with MA corruption; (b) corresponding frequency spectrums.

  • Figure 2

    (Color online) (a) Clean PPG signal and PPG with MA corruption; (b) IMFs extracted from clean signal; protectłinebreak (c) IMFs extracted from PPG with MA corruption.

  • Figure 3

    (Color online) Flowchart for the proposed method. Each processing is shown in module format. Variables passed among modules are shown.

  • Figure 4

    (Color online) AAE varies along ratio $\phi$ changing on 12 samples.

  • Figure 5

    AAE varies along ratio (a) $\phi_{2}$ and (b) $\phi_{1}$ changing on 12 samples.

  • Figure 6

    (Color online) AAE of EEMD-PCA, proposed methods and that without 2-D filtration shown in boxplot view.

  • Figure 7

    (Color online) Pearson correlation between referenced HR and proposed method based PPG derived HR.

  • Figure 8

    (Color online) Effect of 2-D filtration strategy in PPG reconstruction.

  • Table 1   Average AAE and the 1st and 3rd quartiles of the three EMD-based methods
    Method Average AAE $1{\text{st}}$ quartile $3{\text{rd}}$ quartile
    EEMD-PCA[10] 1.78 0.32 2.56
    Proposed method without 2-D filtration 2.56 0.14 3.02
    Proposed method 1.07 0.12 1.41
  • Table 2   AAE (in beat per minute) for each dataset using the proposed method as well as using MISPT, JOSS and TROIKA. Average (SD) is calculated by the test results of each slice
    Strategy Principal methods Dataset Average AAE SD
    TROIKA[4] SSA, SSR, SPT Experimental data 2.34
    MISPT[5] SSA, SSR, SPT Experimental data 1.11 2.33
    JOSS[18] Joint SSR Experimental data 1.28 2.61
    Proposed method EMD, VCS, HWT, 2-D filtration MIMIC II 1.07 1.87

Contact and support