Generalized Deep Neural Network Model for Cuffless Blood Pressure Estimation with Photoplethysmogram Signal Only
Abstract
:1. Introduction
2. Literature Review
3. Materials and Methods
- Data preprocessing: This part comprises signal smoothing of raw PPG data and the removal of abnormal data following standard procedures suggested by [23]. Next, we partition the preprocessed PPG into an approximately 2.17 million heart cycles.
- Feature extractor: Features from the preprocessed data are further extracted and selected as the input set.
- Deep neural network predictor: We feed the feature set into a deep neural network predictor, which consists of five fully connected layers, and each layer contains 2000+ units of fully connected perceptrons, responsible for predicting BP in each heart cycle from 32 extracted physiological parameters.
3.1. Data Source
3.2. PPG Raw Data Preprocessing
- Noise removal: Fast Fourier transform (FFT) is applied to every PPG data segment to convert it from its time domain into the frequency domain. Let , represent the PPG, and the FFT of x[n] is denoted as . We remove the frequency components that are lower than 0 Hz or higher than 8 Hz by turning off those frequency components, as follows
- Normalization and 1st and 2nd derivative of PPG calculation (denoted as “dPPG” and “sdPPG”): All the raw values of PPG are positive, so min–max normalization is applied to every PPG data segment. The equation of min–max normalization can be represented as (5):
- Feature point detection: Before feature extraction, a few points should be marked and detected in every cycle of the heartbeat for every signal (PPG, dPPG and sdPPG) for cycle segmentation and alignment. Firstly, the systolic peaks of PPG waves of each heart cycle are marked by taking advantage of an algorithm mentioned in [24]. The correctness and validity of the systolic peak detection algorithm is of vital importance because the rest of the feature point detection algorithm is based on it. Secondly, the onset and offset valley points of PPG are detected by finding the minimum between two consecutive systolic peaks. Thirdly, with the valley points of PPG found, the location with the maximal and minimal slope values of PPG and dPPG can easily be derived by computing their gradients. Fourthly, the dicrotic notch points of PPG are detected by finding the secondary peaks of the sdPPG contour [20]. An example set of waveforms is shown in Figure 4.
- Partitioning and abnormal cycle removal: After feature points are located, each PPG data segment and its corresponding dPPG and sdPPG waves are partitioned into fragments by reserving each PPG data segment from one valley point of PPG to the next consecutive valley point of PPG. Abnormal heart cycles are also removed following the criteria mentioned in [23]. After abnormal cycle removal is done, the histograms of distribution of SBP and DBP are plotted, as seen in Figure 5, and approximately 2.17 million PPG, dPPG and sdPPG data fragments of heart cycles are obtained.
3.3. Feature Extraction and Selection Index
3.3.1. Feature Extraction
3.3.2. Selection Index
- f = probability mass function of standardized target feature, its estimated precision is down to k decimal places.
- (u) = probability density function of standard normal distribution.
- = .
- = 0, the mean of the standardized target feature.
- = C, where is the standard deviation of the target feature (in the case of standardized features, and are equal to 0 and 1) and C is an integer. For evaluation, the definition of the values of features ranging from – C to + C is used.
3.4. Deep Neural Network Predictor
3.4.1. Introduction to Fully Connected Neural Network
3.4.2. Neural Network Selection
4. Experiments and Results
4.1. Feature Point Detection and Abnormal Cycle Removal
4.2. Characteristic Features of Cardiac Cycles
4.3. Model of Deep Neural Network Predictor
4.4. Performance of Proposed Model
4.4.1. Performance Evaluation by AAMI Standards
4.4.2. Performance Evaluation by BHS Standards
4.4.3. Pearson’s Correlation and Bland–Altman Analysis
4.4.4. Comparison with Other Works
5. Conclusions and Future Works
Author Contributions
Funding
Conflicts of Interest
References
- WHO. Cardiovascular Diseases. 2017. Available online: https://www.who.int/news-room/fact-sheets/detail/cardiovascular-diseases-(cvds)#:~:text=CVDs%20are%20the%20number%201,to%20heart%20attack%20and%20stroke (accessed on 1 May 2019).
- How Accurate Are Drugstore Blood Pressure Machines? Available online: https://www.health.harvard.edu/blood-pressure/how-accurate-are-drugstore-blood-pressure-machines (accessed on 27 July 2020).
- O’brien, E.; Waeber, B.; Parati, G.; Staessen, J.A.; Myers, M.G. Blood pressure measuring devices: Recommendations of the European Society of Hypertension. BMJ 2001, 322, 531–536. [Google Scholar] [CrossRef] [Green Version]
- O’Brien, E.; Petrie, J.; Littler, W.; De Swiet, M.; Padfield, P.L.; O’Malley, K.; Jamieson, M.; Altman, D.; Bland, M.; Atkins, N. The British Hypertension Society protocol for the evaluation of automated and semi-automated blood pressure measuring devices with special reference to ambulatory systems. J. Hypertens. 1990, 8, 607–619. [Google Scholar] [CrossRef] [Green Version]
- Westerhof, N.; Stergiopulos, N.; Noble, M.I.M. Snapshots of Hemodynamics, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar] [CrossRef]
- Bramwell, J.C.; Hill, A.V. The velocity of pulse wave in man. Proc. R. Soc. Lond. Ser. B Contain. Pap. Biol. Character 1922, 93, 298–306. [Google Scholar]
- Tanveer, S.; Hasan, K. Cuffless blood pressure estimation from electrocardiogram and photoplethysmogram using waveform based ANN-LSTM network. Biomed. Signal. Process. Control. 2019, 51, 382–392. [Google Scholar] [CrossRef] [Green Version]
- Hughes, D.J.; Babbs, C.F.; Geddes, L.A.; Bourland, J.D. Measurements of Young’s modulus of elasticity of the canine aorta with ultrasound. Ultrason. Imaging 1979, 1, 356–367. [Google Scholar] [CrossRef] [Green Version]
- Kachuee, M.; Kiani, M.M.; Mohammadzade, H.; Shabany, M. Cuff-less high-accuracy calibration-free blood pressure estimation using pulse transit time. In Proceedings of the International Symposium on Circuits and Systems (ISCAS), Lisbon, Portugal, 24–27 May 2015; pp. 1006–1009. [Google Scholar]
- Kachuee, M.; Kiani, M.M.; Mohammadzade, H.; Shabany, M. Cuffless Blood Pressure Estimation Algorithms for Continuous Health-Care Monitoring. IEEE Trans. Biomed. Eng. 2017, 64, 859–869. [Google Scholar] [CrossRef]
- Zhang, Y.; Poon, C.C.Y.; Chan, C.; Tsang, M.W.W.; Wu, K. A Health-Shirt using e-Textile Materials for the Continuous and Cuffless Monitoring of Arterial Blood Pressure. In Proceedings of the 3rd IEEE/EMBS International Summer School on Medical Devices and Biosensors, Cambridge, MA, USA, 4–6 September 2006; pp. 86–89. [Google Scholar]
- Isakadze, N.; Martin, S. How useful is the smartwatch ECG? Trends Cardiovasc. Med. 2019, 30, 442–448. [Google Scholar] [CrossRef]
- De Moraes, J.L.; Rocha, M.X.; Vasconcelos, G.G.; Filho, J.E.D.V.; De Albuquerque, V.H.C.; De Alexandria, A.R. Advances in Photopletysmography Signal Analysis for Biomedical Applications. Sensors 2018, 18, 1894. [Google Scholar] [CrossRef] [Green Version]
- Khalid, S.G.; Zhang, J.; Chen, F.; Zheng, D. Blood Pressure Estimation Using Photoplethysmography Only: Comparison between Different Machine Learning Approaches. J. Health Eng. 2018, 2018, 1–13. [Google Scholar] [CrossRef] [Green Version]
- Ibtehaz, N.; Rahman, M.S. PPG2ABP: Translating Photoplethysmogram (PPG) Signals to Arterial Blood Pressure (ABP) Waveforms using Fully Convolutional Neural Networks. arXiv 2020, arXiv:2005.01669. [Google Scholar]
- Lin, W.-H.; Li, X.; Li, Y.; Li, G.; Chen, F. Investigating the physiological mechanisms of the photoplethysmogram features for blood pressure estimation. Physiol. Meas. 2020, 41, 044003. [Google Scholar] [CrossRef] [PubMed]
- Wang, C.; Yang, F.; Yuan, X.; Zhang, Y.; Chang, K.; Li, Z. An End-to-End Neural Network Model for Blood Pressure Estimation Using PPG Signal. In Artificial Intelligence in China; Springer: Singapore, 2020; pp. 262–272. [Google Scholar] [CrossRef]
- Baldoumas, G.; Peschos, D.; Tatsis, G.; Chronopoulos, S.K.; Christofilakis, V.; Kostarakis, P.; Varotsos, P.; Sarlis, N.V.; Skordas, E.S.; Bechlioulis, A.; et al. A Prototype Photoplethysmography Electronic Device that Distinguishes Congestive Heart Failure from Healthy Individuals by Applying Natural Time Analysis. Electronics 2019, 8, 1288. [Google Scholar] [CrossRef] [Green Version]
- Sun, Y.; Thakor, N.V. Photoplethysmography Revisited: From Contact to Noncontact, From Point to Imaging. IEEE Trans. Biomed. Eng. 2015, 63, 463–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Elgendi, M. On the analysis of fingertip photoplethysmogram signals. Curr. Cardiol. Rev. 2012, 8, 14–25. [Google Scholar] [CrossRef]
- Elgendi, M.; Liang, Y.; Ward, R. Toward Generating More Diagnostic Features from Photoplethysmogram Waveforms. Diseases 2018, 6, 20. [Google Scholar] [CrossRef] [Green Version]
- Wang, L.; Zhou, W.; Xing, Y.; Zhou, X.-G. A Novel Neural Network Model for Blood Pressure Estimation Using Photoplethesmography without Electrocardiogram. J. Health Eng. 2018, 2018, 1–9. [Google Scholar] [CrossRef]
- Sun, J.X.; Reisner, A.T.; Mark, R.G. A signal abnormality index for arterial blood pressure waveforms. In Proceedings of the 2006 Computers in Cardiology, Valencia, Spain, 17–20 September 2006; pp. 13–16. [Google Scholar]
- Gent, P.V.; Farah, H.; Nes, N.; van Arem, B.V. Heart Rate Analysis for Human Factors: Development and Validation of an Open Source Toolkit for Noisy Naturalistic Heart Rate Data. In Proceedings of the 6th Humanist Conference, The Hague, The Netherlands, 13–14 June 2018; pp. 170–175. [Google Scholar]
- Li, Y.; Wang, Z.; Zhang, L.; Yang, X.; Song, J. Characters available in photoplethysmogram for blood pressure estimation: Beyond the pulse transit time. Australas. Phys. Eng. Sci. Med. 2014, 37, 367–376. [Google Scholar] [CrossRef]
- Ding, X.-R.; Zhang, Y.-T.; Liu, J.; Dai, W.-X.; Tsang, H.K. Continuous Cuffless Blood Pressure Estimation Using Pulse Transit Time and Photoplethysmogram Intensity Ratio. IEEE Trans. Biomed. Eng. 2015, 63, 964–972. [Google Scholar] [CrossRef]
- Fukushima, H.; Kawanaka, H.; Bhuiyan, M.S.; Oguri, K. Cuffless blood pressure estimation using only photoplethysmography based on cardiovascular parameters. In Proceedings of the 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Osaka, Japan, 3–7 July 2013; pp. 2132–2135. [Google Scholar]
- Ding, X.-R.; Yan, B.P.; Zhang, Y.; Liu, J.; Su, P.; Zhao, N. Feature Exploration for Knowledge-guided and Data-driven Approach Based Cuffless Blood Pressure Measurement. arXiv 2019, arXiv:1908.10245. [Google Scholar]
- Su, P.; Ding, X.; Zhang, Y.; Liu, J.; Miao, F.; Zhao, N. Long-term blood pressure prediction with deep recurrent neural networks. In Proceedings of the EMBS International Conference on Biomedical & Health Informatics (BHI), Las Vegas, NV, USA, 4–7 March 2018; pp. 323–328. [Google Scholar]
- Kurylyak, Y.; Lamonaca, F.; Grimaldi, D. A Neural Network-based method for continuous blood pressure estimation from a PPG signal. In Proceedings of the International Instrumentation and Measurement Technology Conference (I2MTC), Minneapolis, MN, USA, 6–9 May 2013; pp. 280–283. [Google Scholar]
- Kingma, D.; Ba, J. Adam: A Method for Stochastic Optimization. In Proceedings of the International Conference on Learning Representations, Banff, AB, Canada, 14–16 April 2014. [Google Scholar]
- Keras, F.C.O. Available online: https://github.com/fchollet/keras (accessed on 5 July 2020).
- Oliphant, T. Numpy: A Guide to Numpy. Available online: http://www.numpy.org/ (accessed on 29 July 2020).
- Botchkarev, A. A New Typology Design of Performance Metrics to Measure Errors in Machine Learning Regression Algorithms. Interdiscip. J. Inf. Knowl. Manag. 2019, 14, 45–76. [Google Scholar] [CrossRef] [Green Version]
- Association for the Advancement of Medical Instrumentation. American National Standard Manual, Electronic or Automated Sphygmonanometers; Association for the Advancement of Medical Instrumentation: Arlington, VA, USA, 2003; Volume AASI/AAMI SP 10:2002; Available online: https://webstore.ansi.org/standards/aami/ansiaamisp102002a12003 (accessed on 4 October 2020).
- Bewick, V.; Cheek, L.; Ball, J. Statistics review 7: Correlation and regression. Crit. Care 2003, 7, 451–459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Bland, D.G.A.M. Measurement in Medicine: The Analysis of Method Comparison Studies. J. R. Stat. Soc. Ser. D 1983, 32, 307. [Google Scholar] [CrossRef]
- Mousavi, S.S.; Charmi, M.; Firouzmand, M.; Hemmati, M.; Moghadam, M.; Ghorbani, Y. Cuff-Less Blood Pressure Estimation Using Only the Photoplethysmography Signal by A Frequency Whole-based Method. In Proceedings of the 2018 8th International Conference on Computer and Knowledge Engineering (ICCKE), Mashhad, Iran, 25–26 October 2018; pp. 147–152. [Google Scholar]
- Slapničar, G.; Mlakar, N.; Luštrek, M. Blood Pressure Estimation from Photoplethysmogram Using a Spectro-Temporal Deep Neural Network. Sensors 2019, 19, 3420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Denotation | Feature | Definition of Feature | |
---|---|---|---|
hr | Heart rate | # 0.01281 * | |
t1 | Time interval of S1 (as seen in Figure 4) | - | |
t2 | Time interval of S2 (as seen in Figure 4) | 0.01653 | |
t3 | Time interval of S3 (as seen in Figure 4) | 0.01637 | |
t4 | Time interval of S4 (as seen in Figure 4) | 0.01604 | |
t5 | Time interval of dAA (as seen in Figure 4) | 0.01634 | |
t6 | Time interval of sdDA (as seen in Figure 4) | 0.01628 | |
t7 | Time interval of sdAA (as seen in Figure 4) | 0.01673 | |
t8 | Time interval of sdDA (as seen in Figure 4) | 0.01628 | |
AS | Ascending slope of PPG (slope from onset point to maximum peak) | # 0.01276 * | |
dAS | Ascending slope of dPPG | 0.01455 | |
sdAS | Ascending slope of sdPPG | # 0.01405 | |
DS | Descending slope of PPG (slope from maximum peak to offset point) | # 0.01405 * | |
dDS | Descending slope of dPPG | 0.01646 | |
sdDS | Descending slope of sdPPG | # 0.01298 | |
S1 | Area under PPG curve between onset point and maximum slope point (as seen in Figure 4) | 0.01556 * | |
S2 | Area under PPG curve between maximum slope point and maximum peak (as seen in Figure 4) | 0.01411 * | |
AA | Ascending area of PPG (as seen in Figure 4) | # 0.01381 * | |
AA | Ascending area of PPG (as seen in Figure 4) | # 0.01381 * | |
dAA | Ascending area of dPPG (as seen in Figure 4) | # 0.01255 * | |
sdAA | Ascending area of sdPPG (as seen in Figure 4) | # 0.01298 * | |
DA | Descending area of PPG (as seen in Figure 4) | # 0.01232 * | |
dDA | Descending area of dPPG (as seen in Figure 4) | # 0.01228 * | |
sdDA | Descending area of sdPPG (as seen in Figure 4) | # 0.01265 * | |
RAAD | Ratio of ascending area to descending area, AA/DA | - | |
dRAAD | dAA/dDA | - | |
sdRAAD | sdAA/sdDA | - | |
PI | Peak intensity of PPG | # 0.01261 * | |
dPI | Peak intensity of dPPG | # 0.01313 * | |
sdPI | Peak intensity of sdPPG | # 0.01305 * | |
dVI | Valley intensity of dPPG | # 0.01296 * | |
sdVI | Valley intensity of sdPPG | # 0.01299 * | |
AID | Ascending intensity difference of PPG, intensity difference between maximum peak and onset point | # 0.01324 * | |
dAID | Ascending intensity difference of dPPG, intensity difference between maximum peak and onset point | # 0.01311 * | |
sdAID | Ascending intensity difference of sdPPG, intensity difference between maximum peak and onset point | # 0.01305 * | |
dDID | Descending intensity difference of dPPG, intensity difference between offset point and maximum peak | # 0.01322 * | |
sdDID | Descending intensity difference of sdPPG, intensity difference between offset point and maximum peak | # 0.01310 * | |
PIR | Peak intensity ratio of PPG, ratio of maximum peak intensity to onset intensity | - | |
dPIR | Peak intensity ratio of dPPG, ratio of maximum peak intensity to onset intensity | - | |
sdPIR | Peak intensity ratio of sdPPG, ratio of maximum peak intensity to onset intensity | - | |
dRIPV | Ratio of maximum peak intensity to valley intensity of dPPG | # 0.01305 * | |
sdRIPV | Ratio of maximum peak intensity to valley intensity of sdPPG | # 0.01350 * | |
AT | Ascending time interval of PPG | # 0.01348 * | |
dAT | Ascending time interval of sPPG | 0.01634 | |
sdAT | Ascending time interval of sdPPG | 0.01673 | |
DT | Descending time interval of PPG | 0.01490 | |
dDT | Descending time interval of dPPG | 0.01628 | |
sdDT | Descending time interval of sdPPG | 0.01628 | |
dTVO | Time interval between valley point and offset point of dPPG | 0.01569 | |
sdTVO | Time interval between valley point and offset point of sdPPG | 0.01438 | |
Slope_a | Slope from maximum peak to dicrotic notch of PPG | # 0.01308 * | |
S3 | Area under PPG curve between maximum peak and dicrotic notch (as seen in Figure 4) | # 0.01333 * | |
S4 | Area under PPG curve between dicrotic notch and offset point (as seen in Figure 4) | # 0.01323 * | |
RtArea | Ratio of systolic area to diastolic area, (S1 + S2 + S3)/S4 (as seen in Figure 4) | - | |
NI | Dicrotic notch intensity | # 0.01230 * | |
AI | Augmentation index = NI/PI | # 0.01277 * | |
AI1 | Augmentation index 1 = (PI − NI)/PI | # 0.01274 * | |
RSD | Ratio of systolic duration to diastolic duration, (t1 + t2 + t3)/t4 | # 0.01405 * | |
RSC | Ratio of diastolic duration to cardiac cycle, t4/(t1 + t2 + t3 +t4) | # 0.01286 * | |
RDC | Ratio of systolic duration to cardiac cycle, (t1 + t2 + t3)/(t1 + t2 + t3 + t4) | 0.01611 |
Error ≤ 5 mmHg | Error ≤ 10 mmHg | Error ≤ 15 mmHg | ||
---|---|---|---|---|
BHS [4] | Grade A | 60% | 85% | 95% |
Grade B | 50% | 75% | 90% | |
Grade C | 40% | 65% | 85% | |
Our Model | SBP | 80.63% | 95.86% | 98.78% |
DBP | 90.19% | 98.29% | 99.59% |
Researchers | Dataset | Input | Performance |
---|---|---|---|
Mousavi et al. [38] | MIMIC II (441 subjects) | PPG | BHS standard: Grade A for DBP and Grade D for SBP AAMI: only the results of DBP satisfy the standards MAE, RMSE: not mentioned in the paper |
Slapnivcar et al. [39] | MIMIC II (510 subjects) | PPG | MAE: DBP = 9.43 mmHg, SBP = 6.88 mmHg RMSE: not mentioned in the paper |
Ibtehaz and Rahman [15] | MIMIC II (942 subjects) | PPG | BHS standard: Grade A for DBP and Grade B for SBP AAMI: the results of both DBP and SBP satisfy the standards MAE: DBP = 3.45 mmHg, SBP = 5.73 mmHg RMSE: not mentioned in the paper |
Our proposed model | MIMIC II (2,176,188 records of BP in total) | PPG | BHS standard: Grade A for both DBP and SBP AAMI: the results of both DBP and SBP satisfy the standards MAE: DBP = 2.23 mmHg, SBP = 3.21 mmHg RMSE: DBP = 3.21 mmHg, SBP = 4.63 mmHg |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Hsu, Y.-C.; Li, Y.-H.; Chang, C.-C.; Harfiya, L.N. Generalized Deep Neural Network Model for Cuffless Blood Pressure Estimation with Photoplethysmogram Signal Only. Sensors 2020, 20, 5668. https://doi.org/10.3390/s20195668
Hsu Y-C, Li Y-H, Chang C-C, Harfiya LN. Generalized Deep Neural Network Model for Cuffless Blood Pressure Estimation with Photoplethysmogram Signal Only. Sensors. 2020; 20(19):5668. https://doi.org/10.3390/s20195668
Chicago/Turabian StyleHsu, Yan-Cheng, Yung-Hui Li, Ching-Chun Chang, and Latifa Nabila Harfiya. 2020. "Generalized Deep Neural Network Model for Cuffless Blood Pressure Estimation with Photoplethysmogram Signal Only" Sensors 20, no. 19: 5668. https://doi.org/10.3390/s20195668
APA StyleHsu, Y. -C., Li, Y. -H., Chang, C. -C., & Harfiya, L. N. (2020). Generalized Deep Neural Network Model for Cuffless Blood Pressure Estimation with Photoplethysmogram Signal Only. Sensors, 20(19), 5668. https://doi.org/10.3390/s20195668