Next Article in Journal
Analysis of Plasmonic Sensors Performance Realized by Exploiting Different UV-Cured Optical Adhesives Combined with Plastic Optical Fibers
Next Article in Special Issue
A Vibration Sensing Device Using a Six-Axis IMU and an Optimized Beam Structure for Activity Monitoring
Previous Article in Journal
A Probabilistic Digital Twin for Leak Localization in Water Distribution Networks Using Generative Deep Learning
Previous Article in Special Issue
An Optical Sensory System for Assessment of Residual Cancer Burden in Breast Cancer Patients Undergoing Neoadjuvant Chemotherapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Precision Heart Rate Estimation Using a PPG Sensor Patch Equipped with New Algorithms of Pre-Quality Checking and Hankel Decomposition

Department of Electronics and Electrical Engineering, National Yang Ming Chiao Tung University, Hsinchu 300, Taiwan
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(13), 6180; https://doi.org/10.3390/s23136180
Submission received: 3 June 2023 / Revised: 24 June 2023 / Accepted: 28 June 2023 / Published: 5 July 2023
(This article belongs to the Special Issue Advanced Sensors for Real-Time Monitoring Applications ‖)

Abstract

:
A new method for accurately estimating heart rates based on a single photoplethysmography (PPG) signal and accelerations is proposed in this study, considering motion artifacts due to subjects’ hand motions and walking. The method comprises two sub-algorithms: pre-quality checking and motion artifact removal (MAR) via Hankel decomposition. PPGs and accelerations were collected using a wearable device equipped with a PPG sensor patch and a 3-axis accelerometer. The motion artifacts caused by hand movements and walking were effectively mitigated by the two aforementioned sub-algorithms. The first sub-algorithm utilized a new quality-assessment criterion to identify highly noise-contaminated PPG signals and exclude them from subsequent processing. The second sub-algorithm employed the Hankel matrix and singular value decomposition (SVD) to effectively identify, decompose, and remove motion artifacts. Experimental data collected during hand-moving and walking were considered for evaluation. The performance of the proposed algorithms was assessed using the datasets from the IEEE Signal Processing Cup 2015. The obtained results demonstrated an average error of merely 0.7345 ± 8.1129 beats per minute (bpm) and a mean absolute error of 1.86 bpm for walking, making it the second most accurate method to date that employs a single PPG and a 3-axis accelerometer. The proposed method also achieved the best accuracy of 3.78 bpm in mean absolute errors among all previously reported studies for hand-moving scenarios.

1. Introduction

The first on-market PPG (photoplethysmogram) sensor was the pulse oximeter, which was commercialized around 1980 for hospital use [1]. This PPG sensor device consists of two alternating LEDs in red/infrared and a photodetector (PD) to obtain information about blood vessel pulsations from the output signals of the PD. The devices are now available on the market in the form of wearable smart watches, earphones, etc. [2]. The non-invasiveness, ease of use, and variety in price, design, and uses have made this industry of wearable smart watches flourish within a few years. However, further advances of this PPG technology are nowadays seriously hampered by the unavoidable contamination of PPG signals by noises mainly from motion artifacts (MAs) [2,3,4,5], which are caused by significantly unavoidable changes in the optical power paths from LEDs and PDs during motions, resulting in undesired components to be measured in the PPG signals. These noises become worse when the sensor is dislocated or in loose contact with the body. With MAs, obtaining the correct heart rate (HR), blood oxygen saturation (SpO2), and blood pressure (BP) based on PPGs becomes very difficult and sometimes impossible due to the absence of uncontaminated PPG signals measured from wearable devices.
Many reported studies have been dedicated to detecting, evaluating, mitigating, and/or removing motion artifacts (MAs) from PPG signals. The detection of MAs can be carried out by statistical means. Rajet Krishnan et al. [3] used kurtosis and skewness in both time and frequency domains to distinguish between clean and MA-corrupted data. Some researchers, such as Hanyu et al. [6], have also used the standard deviation and mean error along with selected parameters to detect MAs and remove them. Another method is to consider the spectral analysis of measured raw PPG signals based on ensemble empirical mode decomposition (EEMD) with spectrum subtraction (SS) [7,8]. In this method, corrupted PPG signals are classified into corrupted, moderately corrupted, and clean by the thresholds on the amplitudes and frequencies of dominant peaks. Bashar et al. [9] developed a method of variable frequency complex demodulation (VFCDM) and applied this to a set of 200 subjects for PPGs measured at fingers and wrists as well as with elbow movements. They acquired 449 recordings, with the result that 156 were misclassified and 29 were false-positives. Once detected, the MAs can be removed. Kong et al. [10] also used VFCDM. Some other studies, such as [11,12], conducted motion artifact removal (MAR) with assistance from reference signals from multiple LEDs (blind separation methods, Kalman filtering, etc.) and/or an accelerometer (adaptive filter). With multiple LEDs, Raghuram et al. [11] and Hara et al. [12] showed an accuracy of 0.392 bpm using PPG with finger motions, achieving a RMSE of 6.5 for walking, running, and jumping. On the other hand, using accelerometers, Lin and Ma [13] adopted discrete wavelet transforms for noise reduction. They calculated and tracked heart rates using the Kalman Filter. Mahdi Boloursaz et al. [14] utilized accelerometer signals before applying an LMS filter for MAR. Other techniques include an independent component analysis (ICA), adopted by [15,16,17,18,19]; empirical mode decomposition (EMD), used by [11,20,21,22,23,24,25]; least mean squares (LMS), used by [26,27,28,29]; wavelet transform (WT), used by [30,31,32]; a notch filter, used by [33]; or recently, in 2022, adaptive filters, used by [34]. With advances in technology, some very recent attempts at using machine learning for MAR have also been reported [35,36,37,38,39]. Most recently, in 2023, attention has been turned to lightweight machine learning modules [40] that are implementable in wearable devices.
Different from all the previous studies, this study proposes two sub-algorithms in a series to remove motion artifacts, which were applied to cases of hand movements and walking. The first was a new quality-assessment criterion to disregard highly noise-contaminated PPG signals, while the second employed the Hankel matrix and the associated SVD to remove motion artifacts. The aforementioned pre-screening of noises by the first quality-check sub-algorithm, and the capability of the subsequent Hankel matrices and the associated SVD to identify, decompose, and remove the motion components in measured PPG were expected to deliver a high-precision HR estimation. Two cases of walking and hand movement were considered to demonstrate the effectiveness of the proposed MAR algorithm. The IEEE Signal Processing Cup 2015 dataset was used for performance evaluation. The average error result was 0.7345 ± 8.1129 bpm, with a mean absolute error (MAE) of 1.86 bpm (beats per minute); the second-best of all the reported results. As for hand-moving, this study showed the best accuracy of 3.78 bpm in MAE compared with all the reported studies.
The remainder of this article is organized into four sections. Section 2 describes how the PPG signals were collected. Section 3 details the first algorithm of the quality check on measured PPG signals, while Section 4 describes the algorithm of motion artifact removal using the Hankel matrix and the associated SVD. Section 5 presents the experimental validation of the performance of the proposed method. Finally, Section 6 concludes this study.

2. PPG Signals and Motion Artifacts

2.1. Measuring PPG

A typical PPG signal at the wrist artery, where most wearable devices measure PPG, is shown in Figure 1a, which was measured by the PPG sensor patch developed exclusively by [41], as given in Figure 1b. This typical PPG waveform consists of AC and DC components. The AC component reflects the pulsatile component of blood vessel pulsations synchronized with heart beats, while the DC component is the non-pulsatile component resulting from light absorption in tissue, skin, and bones along the optical paths between LEDs and PD of the PPG devices. To find the best wavelength of LEDs and the best location to measure PPG, a series of experiments with the subject sitting still were conducted with the PPG sensor in Figure 1 emitted at different wavelengths and attached at different locations 1–3, as illustrated in Figure 2a–c. It can be seen from Figure 2a–c that both green and infrared (IR) LEDs lead to larger AC components, while locations 1 and 3 lead to larger ACs in the measured PPGs as opposed to others. This is, in fact, due to the presence of arteries under the locations as shown in Figure 2a, like seen in [42]. For the remainder of this study, to develop the algorithm for removing motion artifacts, the PPG sensor patch in Figure 1b is attached at location 1, as seen in Figure 2b, to obtain the PPG measured from the green LED. Typical PPGs contaminated by motion artifacts, measured by the green LED at location 1 on the wrist artery, are shown in Figure 3a,b, which are those during walking and with hand motion, respectively. Based on the comparison between Figure 3a,b, it can be seen that the PPG waveforms measured with hand movement exhibit much smaller DC drifts than walking, while the AC components of both cases are close to each other. Additionally, the large DC drifts, especially in the case of walking, present a serious, negative effect of motion artifacts on estimating the heart rate, blood oxygen, and/or blood pressure based on measured PPG waveforms.

2.2. Classification of Motion Artifacts

Motion artifacts are understood nowadays as having two different magnitudes, micro and macro motions, as shown in Figure 4. The micro-motion artifacts have been defined as those due to measuring position adjusting, gesture changing, or finger tapping, while the macro-motion artifacts are those having consistent body motion, such as walking, jogging, and running. Figure 4 shows the measured PPGs contaminated by the micro- and macro-motion artifacts due to finger tapping, measuring position adjusting (both leading to micro-MAs), and walking (macro-MAs). It can be seen from this figure that both micro- and macro-motion artifacts cause much more significant fluctuations to the measured PPG signal than sedentary gestures of the subject. All the aforementioned motions could range between 0.1 and 20 Hz, which is right within the frequency range for heart rate (1–4 Hz), causing much trouble for estimating heart rate based on measured PPG signals. Figure 5 shows the frequency response of a contaminated PPG signal with a few components away from the given frequencies. With the presence of micro and macro motions, the estimation accuracy of heart rate could be seriously compromised. To remedy this problem, two algorithms are proposed herein for motion artifact reduction (MAR), as described in the ensuing Section 3 and Section 4.

3. The First Sub-Algorithm of Quality Check on Measured PPG Signals

The algorithm proposed by this study consists of two sub-algorithms to increase the accuracy of heart rate estimation with serious motion artifacts in PPGs measured by the wearable device. The first sub-algorithm is designed to check the quality of real-time raw PPGs measured to rule out those seriously contaminated and keep the information of heart rate still embedded in the measured PPGs. With this first sub-algorithm, the accuracy of estimating heart rate is expected to increase. To this end, a new algorithm for signal quality checking on measured PPG has been developed. Figure 6 elaborates on its computation flow. Note that a qualified PPG signal at the end is defined as the one that is periodic in time domain and has the largest frequency component remaining within a limited range close to that of the legitimate heart rate. Accordingly, specific conditions in a sequence for checking are defined below for subsequent analysis to extract precision heart rate in the next section.
(1)
The number of valleys and peaks of the measured PPG signal is calculated, and then it is checked if
( N u m b e r   o f   V a l l e y s )     ( N u m b e r   o f   P e a k s ) ± 1 ,
to affirm the presence of pulsation in the measured PPG; otherwise, the signal is labeled as “unqualified” and we then return to stage one for motion artifact reduction.
(2)
The ti values are defined as the time intervals between peaks that can be extracted from measured PPGs, as shown in Figure 7. It is next checked if all differences between consecutive ti values,
Δt = ti+1ti,
are each less than or equal to ±0.1 of ti. If this condition is true, the signal is labelled as “qualified,” and then we go to the next step in Figure 6; otherwise, the signal is labelled as “unqualified” and we stop the computation.
(3)
The statistical measures of kurtosis, mean, and standard deviation of the qualified PPG segments are calculated further at each cycle. If all calculated statistical values are within pre-defined thresholds, as shown in Figure 8a, the measured PPG is identified as “qualified;” otherwise, it is determined as as “unqualified” and then the computation is stopped. Note that a similar approach was used by [3] for motion artifact detection.
(4)
A marker SQT (Signal Quality Token) is defined as either ’0’ or ‘1’ to label the measured PPG segment as qualified or unqualified. Figure 8b shows three representative examples of PPG marked with different SQTs. Only the PPG with SQRT = 1 is passed on to calculate heart rate based on the Hankel matrix and SVD decomposition.
Figure 6. Computation flow of the quality-check algorithm on the measured PPG.
Figure 6. Computation flow of the quality-check algorithm on the measured PPG.
Sensors 23 06180 g006
Figure 7. Measured PPG signal in time and frequency domains.
Figure 7. Measured PPG signal in time and frequency domains.
Sensors 23 06180 g007
Figure 8. (a) Mean and kurtosis for qualified and unqualified PPG signals; (b) SQT for the PPG signal.
Figure 8. (a) Mean and kurtosis for qualified and unqualified PPG signals; (b) SQT for the PPG signal.
Sensors 23 06180 g008

4. The Second Sub-Algorithm for Motion Artifact Removal

The PPG signals determined as “qualified” by the first sub-algorithm of quality checking in the previous section are next processed with the motion artifact removal (MAR) algorithm proposed herein.

4.1. The Hankel Matrix and Its SVD

Precision estimation of heart rate is carried out herein by singular value decomposition on the Hankel matrix of the measured PPG and motion via an accelerometer also attached to the PPG sensor, as shown in Figure 1b. The PPG and accelerometer signals are measured and recorded for six seconds and then first organized into Hankel matrices [43], Hppg, Hx, Hy, and Hz, respectively, for PPG and accelerations along x, y, and z. For example, Hppg is in the form of
H ppg = a ( t 1 ) a ( t 2 ) a ( t n ) a ( t 2 ) a ( t n ) a ( t 2 n ) ,
where a(t)’s are PPG data at times of t’s, and there are in total 2n PPG data in the sampled window of six seconds. Next, the single value decompositions (SVDs) [44] are conducted on each of Hppg, Hx, Hy, and Hz, leading to
H ppg = U ppg Σ ppg V T ppg ,
H x = U x Σ x V x T ,
H y = U y Σ y V y T ,
and   H z = U z Σ z V z T
respectively, where U, Σ, and V correspond to orthogonal, diagonal, orthogonal matrices. U contains the orthogonal basis for the column space of H, while V contains the orthogonal basis for the row space of H. Σ contains the eigenvalues of matrix H. The eigenvalues of the diagonal matrices Σ are stored as Λppg, Λx, Λy, and Λz for further processing. In the next steps, the eigenvalues with the associated components highly correlated to motion artifacts are removed. To this end, the correlation matrix of PPG and 3-axis accelerometer signals can be calculated; that is [45],
ρ x , y = corr ( x , y ) = cov ( x , y ) σ x σ y = E ( x - μ x )   E ( y - μ y ) σ x σ y ,
where the correlation value ρx,y gives the similarity index between the two signals, with expected values μx and μy, and standard σx and σy of the two signals. The components in measured PPG that are highly correlated to accelerations are removed. Then, the time-domain PPG signal without motion artifacts can be restored by recovering the frequency response with the non-MR-related components only via the inverse Fourier transform, as seen in Figure 9.

4.2. The Computational Flow of the Proposed MAR Algorithm

A new algorithm of motion artifact removal (MAR) is engineered herein for re-constructing MR-free PPG signals based on the decomposed Hankel matrices in Equations (5)–(7) and their mutual correlation in Equation (8). The associated computation flow is shown in Figure 10, while the pseudo-code is given below. The algorithm (Algorithm 1 in pseudocode) consists of five stages, (a) synthesizing the Hankel matrices Hppg, Hx, Hy, and Hz in the forms of (3); (b) Conducting singular value decomposition (SVD) on H’s; (c) Calculating correlation to remove MR-related components; (d) Conducting discrete Fourier transform on the MR-free PPG signals and then finding the heart rate by taking the maximum peak of the spectrum as the heartbeat component; (e) Finally, the heart rate is estimated again to see if the consecutive estimated heart rates are close to each other to ensure the robustness of the algorithm.
Algorithm 1 In pseudocode: Motion Artifact Algorithm for Walking.
1: Procedure Record PPG signal and accelerometer signal x,y,z for 8 s
2: Initialize HR_est = 78
3: Construct Hankel matrix Hppg for PPG
4: Construct Hankel matrix Hx for x
5: Construct Hankel matrix Hy for y
6: Construct Hankel matrix Hz for z
7: Find SVD of matrix obtained from step 3, 4, 5, 6
8: Construct a correlation matrix between the 3-axis accelerometer and PPG
9: Select eigenvalues
10: Reconstruct using inverse SVD
11: Find DFT of the reconstructed signal
12: Find HR
13: Heart rate estimation
14: End procedure

5. Experimental Results

Experiments were conducted for the cases of hand movement and walking, as shown in Figure 11, to validate the performance of the two sub-algorithms built. A commercially certified oximeter was utilized to provide reference heart rates for validating the performance of the two designed sub-algorithms. All experiments for collecting data were conducted at room temperature, 28~31 °C. The skin temperature of the subject during the experiment is noted to be 35~37 °C. The PPG signals and the accelerometer signals are recorded simultaneously using a printed circuit board (PCB). A Bluetooth connection was established, allowing the subject to freely move his/her arm and walk.

5.1. Hand Movement

The experiment setup with hand movement is shown in Figure 11a,b. A PPG sensor was attached to the wrist artery of a subject, while the fingertip of another hand was clamped with the aforementioned oximeter for reference. The subjects were asked to sit and relax for some time, and then their information, such as skin color, temperature, and reference heart rate from the oximeter, was recorded and noted. Then the subjects were allowed to move the hand up and down with the PPG sensor patch worn at the wrist. The algorithm (consisting of two sub-algorithms) was tested on 10 subjects of three different skin tones (beige, honey, and bronze). The results are shown in Figure 12. The proposed system achieves an accuracy of 0.6525 ± 4.7 bpm with a window of 6 s. The Pearson correlation is 0.6, while the average absolute error (MAE) is 3.78 beats per minute (bpm).

5.2. Walking

The proposed MAR algorithm for walking was tested again on these 10 subjects for walking, as seen in Figure 11c, while the datasets of the IEEE Signal Processing Cup 2015 were used to evaluate the performance of the developed algorithms. The results are shown in Figure 13. Seen in Figure 13a is the Bland–Altman plot of the results predicted by the proposed algorithm for 8-s windows of walking. The resulting average error is as low as 0.7345 ± 8.1129 bpm (beats per minute), with a mean absolute error (MAE) of 1.86 bpm. The associated Pearson correlation is 0.9499, as seen in Figure 13b. Note herein that the MAE of 1.86 bpm for walking being lower than 3.78 bpm for hand-moving is due to much smaller hand-waving amplitudes during walking than intentional hand movement of the subjects.

5.3. Discussion

The performance of the proposed two sub-algorithms is compared herein to the results delivered by other reported past studies. Table 1 [6,13,19] and Table 2 [8,10,22,23,24,25,29,46] show the comparison among the algorithms for hand movement and walking, respectively. The mean absolute error (MAE) is considered the main metric for performance comparison. Table 1 gives the performance comparison with other reported studies on MAR for hand movement, while Table 2 does for walking. It can be seen from Table 1 that the present study shows the best accuracy of 3.78 bpm in mean absolute error (MAE) as compared to all the reported studies. The reasons that the present study renders better results compared to other studies [6,13] are the use of the accelerometer signals as a reference and the much better quality checking offered by this work. In comparison to [19], this work identifies the pulsation from PPG before estimating HR, while the work in [19] relies only on a single accelerometer.
As for walking, the mean absolute error (MAE) offered by our algorithm for walking is 1.86 bpm (beats per minute), the second best but very close to the best [46] of all the reported results. Among all the studies in Table 2, most of them [8,10,22,23,29,46] render better results, with MAEs under 3 bpm, where two-channel PPGs and an accelerometer are employed for HR estimation, the same as in this study, to achieve favorable performance. On top of these studies are those [29,46] that achieve MAEs under 2 bpm for walking. Arunkumar and Bhaskar in 2020 [29] achieved an MAE of 1.89 bpm by developing the recursive least squares (RLS) and normalized least mean squares (NLMS) adaptive filters to remove MAs in the frequency domain. As for the work by Motin et al. in 2019 [46], recursive Wiener filtering was employed. Since Wiener filtering has been considered effective in removing MAs by many studies, the result showed a very favorable low error of HR estimation of 1.85 bpm. Even with very low MAEs [29,46], it could be difficult to remove MAs accurately if the frequency of motion (walking) is close to HR. The methods of Hankel and SVD proposed by this study are based on correlations, which are supposedly more capable of identifying and then removing MAs in the frequencies close to HR. In a nutshell, the favorable precision of HR estimation is due to the pre-screening of noises by the quality-check sub-algorithm and the capability of the subsequent Hankel matrices and the associated SVD to identify, decompose, and remove the components in measured PPGs that are highly correlated to motions.

6. Conclusions

Effort was dedicated to developing new algorithms for motion artifact removal (MAR) to accurately estimate heart rate when the subject is in motion, such as hand movement and walking. The proposed algorithm consists of two sub-algorithms. The first is a new quality-assessment criteria to disregard highly noise-contaminated PPG signals, while the second employs the Hankel matrix and SVD to remove motion artifacts. The second sub-algorithm for MAR (motion artifact removal) is built upon (1) assembling Hankel matrices of PPG and accelerations and (2) singular value decomposition to remove the frequency components related to motions for accurate heart rate estimation. The result shows a mean absolute error (MAE) of 1.86 bpm (beats per minute) for walking, the second best but close to the best among all the reported results. As for hand movement, the algorithm shows the best accuracy of 3.78 bpm in MAE as compared to all the other reported results.

Author Contributions

Conceptualization, P.C.-P.C.; methodology, S.T.; software, S.T.; validation, S.T.; formal analysis, S.T.; investigation, S.T.; resources, P.C.-P.C.; data curation, S.T.; writing—original draft preparation, C.-H.T.; writing—review and editing, P.C.-P.C. and C.-H.T.; visualization, S.T., P.C.-P.C. and C.-H.T.; supervision, P.C.-P.C.; project administration, P.C.-P.C.; funding acquisition, P.C.-P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology, Taiwan grant numbers 111-2223-E-A49-005, 111-2221-E-A49-159-MY3, 112-2811-E-A49-508-MY2 and 112-2223-E-A49-006, and by the Hsinchu and Southern Taiwan Science Park Bureaus, Ministry of Science and Technology, Taiwan grant numbers 108A31B, 110CE-2-02 and 112AO28B.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of National Chiao Tung University (protocol code NCTU-REC-108-111E and date of approval: 7 July 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study is available on request from the corresponding author. The data is not publicly available due to privacy and ethical restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zenian, A. Technology Spotlight: Taking the Pulse of the Pulse Oximetry Market. Available online: https://24x7mag.com/medical-equipment/patient-care-equipment/technology-spotlight-taking-pulse-pulse-oximetry-market/ (accessed on 17 March 2021).
  2. Castaneda, D.; Esparza, A.; Ghamari, M.; Soltanpur, C.; Nazeran, H. A review on wearable photoplethysmography sensors and their potential future applications in health care. Int. J. Biosens. Bioelectron. 2018, 4, 195. [Google Scholar] [PubMed] [Green Version]
  3. Krishnan, R.; Natarajan, B.; Warren, S. Analysis and detection of motion artifact in photoplethysmographic data using higher order statistics. In Proceedings of the 2008 IEEE International Conference on Acoustics, Speech and Signal Processing, Las Vegas, NV, USA, 31 March–4 April 2008; pp. 613–616. [Google Scholar]
  4. Lee, Y.; Shin, H.; Jo, J.; Lee, Y.-K. Development of a wristwatch-type PPG array sensor module. In Proceedings of the 2011 IEEE International Conference on Consumer Electronics-Berlin (ICCE-Berlin), Berlin, Germany, 6–8 September 2011; pp. 168–171. [Google Scholar]
  5. Thomas, S.S.; Nathan, V.; Zong, C.; Soundarapandian, K.; Shi, X.; Jafari, R. BioWatch: A noninvasive wrist-based blood pressure monitor that incorporates training techniques for posture and subject variability. IEEE J. Biomed. Health Inform. 2015, 20, 1291–1300. [Google Scholar] [CrossRef]
  6. Shao, H.Y.; Chen, X.H. Motion artifact detection and reduction in PPG signals based on statistics analysis. In Proceedings of the 2017 29th Chinese control and decision conference (CCDC), Chongqing, China, 28–30 May 2017; pp. 3114–3119. [Google Scholar]
  7. Zhang, Y.; Liu, B.; Zhang, Z. Combining ensemble empirical mode decomposition with spectrum subtraction technique for heart rate monitoring using wrist-type photoplethysmography. Biomed. Signal Process. Control 2015, 21, 119–125. [Google Scholar] [CrossRef]
  8. Friman, S.; Vehkaoja, A.; Perez-Macias, J.M. The use of wrist EMG increases the PPG Heart Rate accuracy in smartwatches. IEEE Sens. J. 2022, 22, 24197–24204. [Google Scholar] [CrossRef]
  9. Bashar, S.K.; Han, D.; Soni, A.; McManus, D.D.; Chon, K.H. Developing a novel noise artifact detection algorithm for smartphone PPG signals: Preliminary results. In Proceedings of the 2018 IEEE EMBS International Conference on Biomedical & Health Informatics (BHI), Las Vegas, NV, USA, 4–7 March 2018; pp. 79–82. [Google Scholar]
  10. Kong, Y.; Chon, K.H. Heart rate tracking using a wearable photoplethysmographic sensor during treadmill exercise. IEEE Access 2019, 7, 152421–152428. [Google Scholar] [CrossRef]
  11. Raghuram, M.; Sivani, K.; Reddy, K.A. Use of complex EMD generated noise reference for adaptive reduction of motion artifacts from PPG signals. In Proceedings of the 2016 International Conference on Electrical, Electronics, and Optimization Techniques (ICEEOT), Chennai, India, 3–5 March 2016; pp. 1816–1820. [Google Scholar]
  12. Hara, S.; Shimazaki, T.; Okuhata, H.; Nakamura, H.; Kawabata, T.; Cai, K.; Takubo, T. Parameter optimization of motion artifact canceling PPG-based heart rate sensor by means of cross validation. In Proceedings of the 2017 11th International Symposium on Medical Information and Communication Technology (ISMICT), Lisbon, Portugal, 6–8 February 2017; pp. 73–76. [Google Scholar]
  13. Lin, W.-J.; Ma, H.-P. A physiological information extraction method based on wearable PPG sensors with motion artifact removal. In Proceedings of the 2016 IEEE International Conference on Communications (ICC), Kuala Lumpur, Malaysia, 22–27 May 2016; pp. 1–6. [Google Scholar]
  14. Mashhadi, M.B.; Asadi, E.; Eskandari, M.; Kiani, S.; Marvasti, F. Heart rate tracking using wrist-type photoplethysmographic (PPG) signals during physical exercise with simultaneous accelerometry. IEEE Signal Process. Lett. 2015, 23, 227–231. [Google Scholar] [CrossRef] [Green Version]
  15. Lee, J.; Kim, M.; Park, H.-K.; Kim, I.Y. Motion artifact reduction in wearable photoplethysmography based on multi-channel sensors with multiple wavelengths. Sensors 2020, 20, 1493. [Google Scholar] [CrossRef] [Green Version]
  16. Kim, B.S.; Yoo, S.K. Motion artifact reduction in photoplethysmography using independent component analysis. IEEE Trans. Biomed. Eng. 2006, 53, 566–568. [Google Scholar] [CrossRef] [PubMed]
  17. Lee, H.; Chung, H.; Ko, H.; Lee, J. Wearable multichannel photoplethysmography framework for heart rate monitoring during intensive exercise. IEEE Sens. J. 2018, 18, 2983–2993. [Google Scholar] [CrossRef]
  18. Xiong, J.; Cai, L.; Jiang, D.; Song, H.; He, X. Spectral matrix decomposition-based motion artifacts removal in multi-channel PPG sensor signals. IEEE Access 2016, 4, 3076–3086. [Google Scholar] [CrossRef]
  19. Zhao, C.; Zeng, W.; Hu, D.; Liu, H. Robust heart rate monitoring by a single wrist-worn accelerometer based on signal decomposition. IEEE Sens. J. 2021, 21, 15962–15971. [Google Scholar] [CrossRef]
  20. Yousef, Q.; Reaz, M.; Ali, M.A.M. The analysis of PPG morphology: Investigating the effects of aging on arterial compliance. Meas. Sci. Rev. 2012, 12, 266–271. [Google Scholar] [CrossRef] [Green Version]
  21. Motin, M.A.; Karmakar, C.K.; 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, 22, 766–774. [Google Scholar] [CrossRef] [PubMed]
  22. Koneshloo, A.; Du, D. A novel motion artifact removal method via joint basis pursuit linear Program to accurately monitor heart rate. IEEE Sens. J. 2019, 19, 9945–9952. [Google Scholar] [CrossRef]
  23. He, W.; Ye, Y.; Lu, L.; Cheng, Y.; Li, Y.; Wang, Z. Robust heart rate monitoring for quasi-periodic motions by wrist-type PPG signals. IEEE J. Biomed. Health Inform. 2019, 24, 636–648. [Google Scholar] [CrossRef] [PubMed]
  24. Savaskan, D.A.; Mahanoglu, A.; Soner, B.; Kholmatov, A. Heart Rate Measurement from Motion Compensated Photoplethysmographic Signals. In Proceedings of the 2020 28th Signal Processing and Communications Applications Conference (SIU), Gaziantep, Turkey, 5–7 October 2020; pp. 1–4. [Google Scholar]
  25. Huang, N.; Selvaraj, N. Robust ppg-based ambulatory heart rate tracking algorithm. In Proceedings of the 2020 42nd Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC), Montreal, QC, Canada, 20–24 July 2020; pp. 5929–5934. [Google Scholar]
  26. Temko, A. Accurate heart rate monitoring during physical exercises using PPG. IEEE Trans. Biomed. Eng. 2017, 64, 2016–2024. [Google Scholar] [CrossRef] [Green Version]
  27. Tanweer, K.T.; Hasan, S.R.; Kamboh, A.M. Motion artifact reduction from PPG signals during intense exercise using filtered X-LMS. In Proceedings of the 2017 IEEE International Symposium on Circuits and Systems (ISCAS), Baltimore, MD, USA, 28–31 May 2017; pp. 1–4. [Google Scholar]
  28. Wu, C.-C.; Chen, I.-W.; Fang, W.-C. An implementation of motion artifacts elimination for PPG signal processing based on recursive least squares adaptive filter. In Proceedings of the 2017 IEEE Biomedical Circuits and Systems Conference (BioCAS), Turin, Italy, 19–21 October 2017; pp. 1–4. [Google Scholar]
  29. Arunkumar, K.; Bhaskar, M. Robust de-noising technique for accurate heart rate estimation using wrist-type PPG signals. IEEE Sens. J. 2020, 20, 7980–7987. [Google Scholar] [CrossRef]
  30. Zhang, Y.; Song, S.; Vullings, R.; Biswas, D.; Simões-Capela, N.; Van Helleputte, N.; Van Hoof, C.; Groenendaal, W. Motion artifact reduction for wrist-worn photoplethysmograph sensors based on different wavelengths. Sensons 2019, 19, 673. [Google Scholar] [CrossRef] [Green Version]
  31. Bousefsaf, F.; Maaoui, C.; Pruski, A. Continuous wavelet filtering on webcam photoplethysmographic signals to remotely assess the instantaneous heart rate. Biomed. Signal Process. Control 2013, 8, 568–574. [Google Scholar] [CrossRef]
  32. Teng, X.; Zhang, Y. Continuous and noninvasive estimation of arterial blood pressure using a photoplethysmographic approach. In Proceedings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (IEEE Cat. No. 03CH37439), Cancun, Mexico, 17–21 September 2003; Volume 4, pp. 3153–3156. [Google Scholar]
  33. Zhang, Q.; Xie, Q.; Wang, M.; Wang, G. Motion artifact removal for PPG signals based on accurate fundamental frequency estimation and notch filtering. In Proceedings of the 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Honolulu, HI, USA, 18–21 July 2018; pp. 2965–2968. [Google Scholar]
  34. Guo, J.; Chen, X.; Zhao, J.; Zhang, X.; Chen, X. An Effective Photoplethysmography Heart Rate Estimation Framework Integrating Two-Level Denoising Method and Heart Rate Tracking Algorithm Guided by Finite State Machine. IEEE J. Biomed. Health Inform. 2022, 26, 3731–3742. [Google Scholar] [CrossRef]
  35. Shuzan, M.N.I.; Chowdhury, M.H.; Hossain, M.S.; Chowdhury, M.E.; Reaz, M.B.I.; Uddin, M.M.; Khandakar, A.; Mahbub, Z.B.; Ali, S.H.M. A novel non-invasive estimation of respiration rate from motion corrupted photoplethysmograph signal using machine learning model. IEEE Access 2021, 9, 96775–96790. [Google Scholar] [CrossRef]
  36. Biswas, D.; Everson, L.; Liu, M.; Panwar, M.; Verhoef, B.-E.; Patki, S.; Kim, C.H.; Acharyya, A.; Van Hoof, C.; Konijnenburg, M.; et al. CorNET: Deep learning framework for PPG-based heart rate estimation and biometric identification in ambulant environment. IEEE Trans. Biomed. Circuits Syst. 2019, 13, 282–291. [Google Scholar] [CrossRef]
  37. Shin, H.C.; Roth, H.R.; Gao, M.; Lu, L.; Xu, Z.; Nogues, I.; Yao, J.; Mollura, D.; Summers, R.M. Deep convolutional neural networks for computer-aided detection: CNN architectures, dataset characteristics and transfer learning. IEEE Trans. Med. Imaging 2016, 35, 1285–1298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Schäck, T.; Muma, M.; Zoubir, A.M. Computationally efficient heart rate estimation during physical exercise using photoplethysmographic signals. In Proceedings of the 2017 25th European Signal Processing Conference (EUSIPCO), Kos, Greece, 28 August–2 September 2017; pp. 2478–2481. [Google Scholar]
  39. Chowdhury, S.S.; Hyder, R.; Hafiz, M.S.B.; Haque, M.A. Real-time robust heart rate estimation from wrist-type PPG signals using multiple reference adaptive noise cancellation. IEEE J. Biomed. Health Inform. 2016, 22, 450–459. [Google Scholar] [CrossRef]
  40. Yu, S.-N.; Wang, C.-S.; Chang, Y.P. Heart Rate Estimation from Remote Photoplethysmography Based on Light-weight U-Net and Attention Modules. IEEE Access 2023, 11, 54058–54069. [Google Scholar] [CrossRef]
  41. Kao, Y.-H.; Chao, P.C.-P.; Wey, C.-L. Design and validation of a new PPG module to acquire high-quality physiological signals for high-accuracy biomedical sensing. IEEE J. Sel. Top. Quant. 2018, 25, 69000210. [Google Scholar] [CrossRef]
  42. Abreu, S. The Pulse in the Wrist Is Felt Over Which Artery? Available online: https://socratic.org/questions/the-pulse-in-the-wrist-is-felt-over-which-artery (accessed on 1 August 2017).
  43. Peller, V. An Introduction to Hankel Operators. In Hankel Operators and Their Applications; Springer: Berlin/Heidelberg, Germany, 2003; pp. 1–59. [Google Scholar]
  44. Banerjee, S.; Roy, A. Linear Algebra and Matrix Analysis for Statistics; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  45. Rodgers, J.L.; Nicewander, W.A. Thirteen ways to look at the correlation coefficient. Am. Stat. 1988, 42, 59–66. [Google Scholar] [CrossRef]
  46. Motin, M.A.; Karmakar, C.K.; Palaniswami, M. PPG derived heart rate estimation during intensive physical exercise. IEEE Access 2019, 7, 56062–56069. [Google Scholar] [CrossRef]
Figure 1. (a) AC and DC components of a typical PPG signal; (b) PPG sensor patch developed by [41].
Figure 1. (a) AC and DC components of a typical PPG signal; (b) PPG sensor patch developed by [41].
Sensors 23 06180 g001
Figure 2. (a) Human arteries at wrist; (b) three different locations on the wrist to measure PPG; (c) measured PPGs using green, red, and infrared LEDs at three different locations.
Figure 2. (a) Human arteries at wrist; (b) three different locations on the wrist to measure PPG; (c) measured PPGs using green, red, and infrared LEDs at three different locations.
Sensors 23 06180 g002
Figure 3. PPGs measured at a wrist artery contaminated by two different motion artifacts: (a) walking; (b) hand movement.
Figure 3. PPGs measured at a wrist artery contaminated by two different motion artifacts: (a) walking; (b) hand movement.
Sensors 23 06180 g003
Figure 4. PPG contaminated by motion and mispositioning.
Figure 4. PPG contaminated by motion and mispositioning.
Sensors 23 06180 g004
Figure 5. Frequency response of the measured PPG signal.
Figure 5. Frequency response of the measured PPG signal.
Sensors 23 06180 g005
Figure 9. Estimating heart rate for hand movement using the proposed algorithms.
Figure 9. Estimating heart rate for hand movement using the proposed algorithms.
Sensors 23 06180 g009
Figure 10. Computation flow of the proposed second sub-algorithm.
Figure 10. Computation flow of the proposed second sub-algorithm.
Sensors 23 06180 g010
Figure 11. (a) Experiment setup for (b) hand movement and (c) walking.
Figure 11. (a) Experiment setup for (b) hand movement and (c) walking.
Sensors 23 06180 g011
Figure 12. (a) Bland–Atman plot for heart rate estimation. (b) Correlation plot of estimated HRs vs. ground-truth HRs for the case of hand movement.
Figure 12. (a) Bland–Atman plot for heart rate estimation. (b) Correlation plot of estimated HRs vs. ground-truth HRs for the case of hand movement.
Sensors 23 06180 g012
Figure 13. (a) Bland-Atman plot for heart rate estimation; (b) Correlation plot of estimated HRs vs ground-truth HRs for the case of walking.
Figure 13. (a) Bland-Atman plot for heart rate estimation; (b) Correlation plot of estimated HRs vs ground-truth HRs for the case of walking.
Sensors 23 06180 g013
Table 1. Comparison among various techniques for MAR on hand movement.
Table 1. Comparison among various techniques for MAR on hand movement.
YearTechniqueSensors (Database)Reference SignalMovementMean Absolute Error in (bpm) Mean Error (bpm)Measurement Location
This study2023Quality Check and Notch Filtering with peak selection and current and gain tuningOne-channel PPG and 3-axis accelerometer signals recorded in the labAccelerometerWaving the hand3.78
95% of HR estimation within ± 9.3 bpm
0.6525Wrist
Lin and Ma [13]2016DWTPPG signalsNoneWaving the hand6.87NANA
Hanyu and Xiao hui [6] 2017Statistical
Evaluation
PPG signalsNoneFinger tapping or hand swinging7.85NANA
Chao Zhao et al. [19]2021ICA, VMD, WSST, SSA, and Kalman SmoothingA three-axis acceleration signalsNoneFinger tapping or hand swinging95% of HR estimation within ±8.86 bpmNAWrist
Table 2. Performance comparison among various techniques for MAR on walking.
Table 2. Performance comparison among various techniques for MAR on walking.
YearTechniqueSensors (Database)Reference SignalMovementMean Absolute Error (bpm) Mean Error
(bpm)
Measurement Location
This Study2023Hankel Matrix, SVD and Spectral AnalysisTwo-channel PPG signals, three-axis accelerationsAccelerometer and a single PPGWalking 1.860.7345Wrist
Amirhossein Koneshloo et al. [22]2019Joint Basis Pursuit Linear ProgramTwo-channel PPG signals, three-axis accelerationsAccelerometer and PPG signal.Walking and running2.61NAWrist
Mohammod Abdul Motin et al. [46]2019Recursive Wiener FilteringTwo-channel PPG signals, three-axis accelerationsAccelerometer and PPG signalWalking and running1.85NAWrist
Wenwen He et al. [23]2020Motion tracking, Sparse Representation-based MA elimination, and Spectral Peak Tracking for HRPPG signals with 3-axis accelerometer signalAccelerometerQuasi-periodic motions.2.40NAWrist
Deniz Alp Savaskan et al. [24]2020SPECMAR, TROIKA and JOSS methods along with pre and post processingTwo-channel PPG signals, three-axis acceleration signals for 12 samplesAccelerometer and PPG signalWalking and running4.19NAWrist
Youngsun Kong et al. [10]2019VFCDM approach, Cubic SplineTwo-channel PPG signals, three-axis acceleration signalsAccelerometer and PPG signalWalking and running2.94NAWrist
Two-channel PPG signals, three-axis acceleration signals (lab)Accelerometer and PPG signalWalking and runningNAForehead
Nicholas Huang et al. [25]2020TAPIR MethodTwo-channel PPG signals, three-axis acceleration signalsAccelerometer
and PPG signal
Walking9.21NAWrist
K.R. Arunkumar et al. [29]2020Recursive Least Squares (RLS) and Normalized Least Mean Squares (NLMS)Two-channel PPG signals, three-axis acceleration signals recorded for 23 samplesAccelerometer and PPG signal.Walking and running1.89NAWrist
S. Friman et al. [8]2022Electromyogram (EMG) and accelerometer (ACC)PPG signals with three-axis acceleration signals Accelerometer, EMG and PPG signalWalking and running2.83NAWrist
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Thakur, S.; Chao, P.C.-P.; Tsai, C.-H. Precision Heart Rate Estimation Using a PPG Sensor Patch Equipped with New Algorithms of Pre-Quality Checking and Hankel Decomposition. Sensors 2023, 23, 6180. https://doi.org/10.3390/s23136180

AMA Style

Thakur S, Chao PC-P, Tsai C-H. Precision Heart Rate Estimation Using a PPG Sensor Patch Equipped with New Algorithms of Pre-Quality Checking and Hankel Decomposition. Sensors. 2023; 23(13):6180. https://doi.org/10.3390/s23136180

Chicago/Turabian Style

Thakur, Smriti, Paul C.-P. Chao, and Cheng-Han Tsai. 2023. "Precision Heart Rate Estimation Using a PPG Sensor Patch Equipped with New Algorithms of Pre-Quality Checking and Hankel Decomposition" Sensors 23, no. 13: 6180. https://doi.org/10.3390/s23136180

APA Style

Thakur, S., Chao, P. C. -P., & Tsai, C. -H. (2023). Precision Heart Rate Estimation Using a PPG Sensor Patch Equipped with New Algorithms of Pre-Quality Checking and Hankel Decomposition. Sensors, 23(13), 6180. https://doi.org/10.3390/s23136180

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop