Estimating instantaneous respiratory rate from the photoplethysmogram

The photoplethysmogram (PPG) obtained from pulse oximetry shows the local changes of blood volume in tissues. Respiration induces variation in the PPG baseline due to the variation in venous blood return during each breathing cycle. We have proposed an algorithm based on the synchrosqueezing transform (SST) to estimate instantaneous respiratory rate (IRR) from the PPG. The SST is a combination of wavelet analysis and a reallocation method which aims to sharpen the time-frequency representation of the signal and can provide an accurate estimation of instantaneous frequency. In this application, the SST was applied to the PPG and IRR was detected as the predominant ridge in the respiratory band (0.1 Hz - 1 Hz) in the SST plane. The algorithm was tested against the Capnobase benchmark dataset that contains PPG, capnography, and expert labelled reference respiratory rate from 42 subjects. The IRR estimation accuracy was assessed using the root mean square (RMS) error and Bland-Altman plot. The median RMS error was 0.39 breaths/min for all subjects which ranged from the lowest error of 0.18 breaths/min to the highest error of 13.86 breaths/min. A Bland-Altman plot showed an agreement between the IRR obtained from PPG and reference respiratory rate with a bias of -0.32 and limits agreement of -7.72 to 7.07. Extracting IRR from PPG expands the functionality of pulse oximeters and provides additional diagnostic power to this non-invasive monitoring tool.

[1]  A. Awad,et al.  The Use of Joint Time Frequency Analysis to Quantify the Effect of Ventilation on the Pulse Oximeter Waveform , 2006, Journal of Clinical Monitoring and Computing.

[2]  Hau-Tieng Wu,et al.  The Synchrosqueezing algorithm for time-varying spectral analysis: Robustness properties and new paleoclimate applications , 2011, Signal Process..

[3]  A. Guyton,et al.  Textbook of Medical Physiology , 1961 .

[4]  J. N. Watson,et al.  Developing an algorithm for pulse oximetry derived respiratory rate (RRoxi): a healthy volunteer study , 2012, Journal of Clinical Monitoring and Computing.

[5]  Guy A. Dumont,et al.  CapnoBase: Signal database and tools to collect, share and annotate respiratory signals , 2010 .

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

[7]  David R. Burton,et al.  Ridge extraction algorithms for one-dimensional continuous wavelet transform: a comparison , 2007 .

[8]  Walter Karlen,et al.  Multiparameter Respiratory Rate Estimation From the Photoplethysmogram , 2013, IEEE Transactions on Biomedical Engineering.

[9]  I. Daubechies,et al.  Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool , 2011 .

[10]  Ingrid Daubechies,et al.  A Nonlinear Squeezing of the Continuous Wavelet Transform Based on Auditory Nerve Models , 2017 .

[11]  A. Buda,et al.  Effect of intrathoracic pressure on left ventricular performance. , 1979, The New England journal of medicine.

[12]  G Y Zou,et al.  Confidence interval estimation for the Bland–Altman limits of agreement with multiple observations per individual , 2013, Statistical methods in medical research.

[13]  J. C. Dorlas,et al.  Photo-electric plethysmography as a monitoring device in anaesthesia. Application and interpretation. , 1985, British journal of anaesthesia.

[14]  John G. Webster,et al.  Design of Pulse Oximeters , 1997 .

[15]  Walter Karlen,et al.  Empirical mode decomposition for respiratory and heart rate estimation from the photoplethysmogram , 2013, Computing in Cardiology 2013.

[16]  L. Tarassenko,et al.  Photoplethysmographic derivation of respiratory rate: a review of relevant physiology , 2012, Journal of medical engineering & technology.

[17]  W. Marsden I and J , 2012 .

[18]  J. N. Watson,et al.  Measurement Of Respiratory Rate From the Photoplethysmogram In Chest Clinic Patients , 2007, Journal of clinical monitoring and computing.