Developing an algorithm for pulse oximetry derived respiratory rate (RR(oxi)): a healthy volunteer study

Advances in Photoplethysmography for Personalized Cardiovascular Monitoring

Photoplethysmography (PPG) is garnering substantial interest due to low cost, noninvasiveness, and its potential for diagnosing cardiovascular diseases, such as cardiomyopathy, heart failure, and arrhythmia. The signals obtained through PPG can yield information based on simple analyses, such as heart rate. In contrast, when accompanied by the complex analysis of sophisticated signals, valuable information, such as blood pressure, sympathetic nervous system activity, and heart rate variability, can be obtained. For a complex analysis, a better understanding of the sources of noise, which create limitations in the application of PPG, is needed to get reliable information to assess cardiovascular health. Therefore, this Special Issue handles literature about noises and how they affect the waveform of the PPG caused by individual variations (e.g., skin tone, obesity, age, and gender), physiology (e.g., respiration, venous pulsation, body site of measurement, and body temperature), and external factors (e.g., motion artifact, ambient light, and applied pressure to the skin). It also covers the issues that still need to be considered in each situation.

Integrating Physiological Data Artifacts Detection With Clinical Decision Support Systems: Observational Study

BACKGROUND

Clinical decision support systems (CDSS) have the potential to lower the patient mortality and morbidity rates. However, signal artifacts present in physiological data affect the reliability and accuracy of the CDSS. Moreover, patient monitors and other medical devices generate false alarms while processing physiological data, further leading to alarm fatigue because of increased noise levels, staff disruption, and staff desensitization in busy critical care environments. This adversely affects the quality of care at the patient bedside. Hence, artifact detection (AD) algorithms play a crucial role in assessing the quality of physiological data and mitigating the impact of these artifacts.

OBJECTIVE

The aim of this study is to evaluate a novel AD framework for integrating AD algorithms with CDSS. We designed the framework with features that support real-time implementation within critical care. In this study, we evaluated the framework and its features in a false alarm reduction study. We developed static framework component models, followed by dynamic framework compositions to formulate four CDSS. We evaluated these formulations using neonatal patient data and validated the six framework features: flexibility, reusability, signal quality indicator standardization, scalability, customizability, and real-time implementation support.

METHODS

We developed four exemplar static AD components with standardized requirements and provisions interfaces that facilitate the interoperability of framework components. These AD components were mixed and matched into four different AD compositions to mitigate the artifacts' effects. We developed a novel static clinical event detection component that is integrated with each AD composition to formulate and evaluate a dynamic CDSS for peripheral oxygen saturation (SpO2) alarm generation. This study collected data from 11 patients with diverse pathologies in the neonatal intensive care unit. Collected data streams and corresponding alarms include pulse rate and SpO2 measured from a pulse oximeter (Masimo SET SmartPod) integrated with an Infinity Delta monitor and the heart rate derived from electrocardiography leads attached to a second Infinity Delta monitor.

RESULTS

A total of 119 SpO2 alarms were evaluated. The lowest achievable SpO2 false alarm rate was 39%, with a sensitivity of 80%. This demonstrates the framework's utility in identifying the best possible dynamic composition to serve the clinical need for false SpO2 alarm reduction and subsequent alarm fatigue, given the limitations of a small sample size.

CONCLUSIONS

The framework features, including reusability, signal quality indicator standardization, scalability, and customizability, allow the evaluation and comparison of novel CDSS formulations. The optimal solution for a CDSS can then be hard-coded and integrated within clinical workflows for real-time implementation. The flexibility to serve different clinical needs and standardized component interoperability of the framework supports the potential for a real-time clinical implementation of AD.

Sources of Inaccuracy in Photoplethysmography for Continuous Cardiovascular Monitoring

Photoplethysmography (PPG) is a low-cost, noninvasive optical technique that uses change in light transmission with changes in blood volume within tissue to provide information for cardiovascular health and fitness. As remote health and wearable medical devices become more prevalent, PPG devices are being developed as part of wearable systems to monitor parameters such as heart rate (HR) that do not require complex analysis of the PPG waveform. However, complex analyses of the PPG waveform yield valuable clinical information, such as: blood pressure, respiratory information, sympathetic nervous system activity, and heart rate variability. Systems aiming to derive such complex parameters do not always account for realistic sources of noise, as testing is performed within controlled parameter spaces. A wearable monitoring tool to be used beyond fitness and heart rate must account for noise sources originating from individual patient variations (e.g., skin tone, obesity, age, and gender), physiology (e.g., respiration, venous pulsation, body site of measurement, and body temperature), and external perturbations of the device itself (e.g., motion artifact, ambient light, and applied pressure to the skin). Here, we present a comprehensive review of the literature that aims to summarize these noise sources for future PPG device development for use in health monitoring.

Estimation of Heart Rate and Respiratory Rate from PPG Signal Using Complementary Ensemble Empirical Mode Decomposition with both Independent Component Analysis and Non-Negative Matrix Factorization

This paper proposes a framework combining the complementary ensemble empirical mode decomposition with both the independent component analysis and the non-negative matrix factorization for estimating both the heart rate and the respiratory rate from the photoplethysmography (PPG) signal. After performing the complementary ensemble empirical mode decomposition on the PPG signal, a finite number of intrinsic mode functions are obtained. Then, these intrinsic mode functions are divided into two groups to perform the further analysis via both the independent component analysis and the non-negative matrix factorization. The surrogate cardiac signal related to the heart activity and another surrogate respiratory signal related to the respiratory activity are reconstructed to estimate the heart rate and the respiratory rate, respectively. Finally, different records of signals acquired from the Medical Information Mart for Intensive Care database downloaded from the Physionet Automated Teller Machine (ATM) data bank are employed for demonstrating the outperformance of our proposed method. The results show that our proposed method outperforms both the digital filtering approach and the conventional empirical mode decomposition based methods in terms of reconstructing both the surrogate cardiac signal and the respiratory signal from the PPG signal as well as both achieving the higher accuracy and the higher reliability for estimating both the heart rate and the respiratory rate.

Non-invasive continuous respiratory monitoring using temperature-based sensors

Respiratory rate (RR) is a key vital sign that has been traditionally employed in the clinical assessment of patients and in the prevention of respiratory compromise. Despite its relevance, current practice for monitoring RR in non-intubated patients strongly relies on visual counting, which delivers an intermittent and error-prone assessment of the respiratory status. Here, we present a novel non-invasive respiratory monitor that continuously measures the RR in human subjects. The respiratory activity of the user is inferred by sensing the thermal transfer between the breathing airflow and a temperature sensor placed between the nose and the mouth. The performance of the respiratory monitor is assessed through respiratory experiments performed on healthy subjects. Under spontaneous breathing, the mean RR difference between our respiratory monitor and visual counting was 0.4 breaths per minute (BPM), with a 95% confidence interval equal to [- 0.5, 1.3] BPM. The robustness of the respiratory sensor to the position is assessed by studying the signal-to-noise ratio in different locations on the upper lip, displaying a markedly better performance than traditional thermal sensors used for respiratory airflow measurements.

Recent development of respiratory rate measurement technologies

Respiratory rate (RR) is an important physiological parameter whose abnormality has been regarded as an important indicator of serious illness. In order to make RR monitoring simple to perform, reliable and accurate, many different methods have been proposed for such automatic monitoring. According to the theory of respiratory rate extraction, methods are categorized into three modalities: extracting RR from other physiological signals, RR measurement based on respiratory movements, and RR measurement based on airflow. The merits and limitations of each method are highlighted and discussed. In addition, current works are summarized to suggest key directions for the development of future RR monitoring methodologies.

Acoustic respiration rate and pulse oximetry-derived respiration rate: a clinical comparison study

Respiration rate (RR) is a critical vital sign that provides early detection of respiratory compromise. The acoustic technique of measuring continuous respiration rate (RR_a) interprets the large airway sound envelope to calculate respiratory rate while pulse oximetry-derived respiratory rate (RR_oxi) interprets modulations of the photoplethsymograph in response to hemodynamic changes during the respiratory cycle. The aim of this study was to compare the performance of these technologies to each other and to a capnography-based reference device. Subjects were asked to decrease their RR from 14 to 4 breaths per minute (BPM) and then increase RR from 14 to 24 BPM. The effects of physiological noise, ambient noise, and head movement and shallow breathing on device performance were also evaluated. The test devices were: (1) RR_a, Radical-7 (Masimo Corporation), (2) RR_oxi, Nellcor™ Bedside Respiratory Patient Monitoring System (Medtronic), and (3) reference device, Capnostream20p™ (Medtronic). All devices were configured with their default settings. Twenty-nine healthy adult subjects were included in the study. During abrupt changes in breathing, overall RR_oxi was accurate for longer periods of time than RR_a; specifically, RR_oxi was more accurate during low and normal RR, but not during high RR. RR_oxi also displayed a value for significantly longer time periods than RR_a when the subjects produced physiological sounds and moved their heads, but not during shallow breathing or ambient noise. RR_oxi may be more accurate than RR_a during development of bradypnea. Also, RR_oxi may display a more reliable RR value during routine patient activities.

Incidence of respiratory depression after epidural administration of morphine for cesarean delivery: findings using a continuous respiratory rate monitoring system

BACKGROUND

Epidural morphine is widely used for postoperative analgesia after cesarean delivery. However, respiratory depression can occur after neuraxial administration of morphine. Previous reports describing respiratory depression in obstetric patients have relied on intermittent visual counting of the respiratory rate. In this study, we estimated the incidence of respiratory depression in patients who had received epidural morphine after cesarean delivery, using a continuous respiratory rate monitoring system with a finger sensor.

METHODS

One hundred patients scheduled to undergo elective cesarean delivery and receive intraoperative neuraxial morphine between April and December 2016 were recruited for this single-center, prospective observational study. Postoperatively, all patients received epidural morphine 3 mg and were equipped with the Nellcor respiratory rate monitoring system. Respiratory depression was defined as both bradypnea (respiratory rate ≤10 breaths/min) and oxygen desaturation (mild ≤95%; moderate ≤90%; severe ≤85%) for longer than one minute. The number of patients with respiratory depression between administration of morphine and first ambulation was recorded hourly.

RESULTS

Complete monitoring was obtained for 89 of 100 women. The median duration of monitoring was 19.0 hours. Forty-six patients (52%) developed mild respiratory depression at least once before ambulation, but only one (1%) developed moderate respiratory depression. None required supplemental oxygen or naloxone.

CONCLUSIONS

Approximately half the women experienced mild respiratory depression, but only one developed moderate respiratory depression. Continuous respiratory rate monitoring until ambulation may assist in early identification of respiratory depression after neuraxial administration of morphine.

Extracting Instantaneous Respiratory Rate From Multiple Photoplethysmogram Respiratory-Induced Variations

In this study, we proposed a novel method for extracting the instantaneous respiratory rate (IRR) from the pulse oximeter photoplethysmogram (PPG). The method was performed in three main steps: (1) a time-frequency transform called synchrosqueezing transform (SST) was used to extract the respiratory-induced intensity, amplitude and frequency variation signals from PPG, (2) the second SST was applied to each extracted respiratory-induced variation signal to estimate the corresponding IRR, and (3) the proposed peak-conditioned fusion method then combined the IRR estimates to calculate the final IRR. We validated the implemented method with capnography and nasal/oral airflow as the reference RR using the limits of agreement (LOA) approach. Compared to simple fusion and single respiratory-induced variation estimations, peak-conditioned fusion shows better performance. It provided a bias of 0.28 bpm with the 95% LOAs ranging from −3.62 to 4.17, validated against capnography and a bias of 0.04 bpm with the 95% LOAs ranging from −5.74 to 5.82, validated against nasal/oral airflow. This algorithm would expand the functionality of a conventional pulse oximetry beyond the measurement of heart rate and oxygen saturation to measure the respiratory rate continuously and instantly.

On algorithms for calculating arterial pulse pressure variation during major surgery

OBJECTIVE

Arterial pulse pressure variation (PPV) is widely used for predicting fluid responsiveness and supporting fluid management in the operating room and intensive care unit. Available PPV algorithms have been typically validated for fluid responsiveness during episodes of hemodynamic stability. Yet, little is known about the performance of PPV algorithms during surgery, where fast changes of the blood pressure may affect the robustness of the presented PPV value. This work provides a comprehensive understanding of how various existing algorithmic designs affect the robustness of the presented PPV value during surgery, and proposes additional processing for the pulse pressure signal before calculating PPV.

APPROACH

We recorded arterial blood pressure waveforms from 23 patients undergoing major abdominal surgery. To evaluate the performance, we designed three clinically relevant metrics. Main results and Significance: The results show that all algorithms performed well during episodes of hemodynamic stability. Moreover, it is demonstrated that the proposed processing helps improve the robustness of PPV during the entire course of surgery.

Breathing Rate Estimation From the Electrocardiogram and Photoplethysmogram: A Review

Breathing rate (BR) is a key physiological parameter used in a range of clinical settings. Despite its diagnostic and prognostic value, it is still widely measured by counting breaths manually. A plethora of algorithms have been proposed to estimate BR from the electrocardiogram (ECG) and pulse oximetry (photoplethysmogram, PPG) signals. These BR algorithms provide opportunity for automated, electronic, and unobtrusive measurement of BR in both healthcare and fitness monitoring. This paper presents a review of the literature on BR estimation from the ECG and PPG. First, the structure of BR algorithms and the mathematical techniques used at each stage are described. Second, the experimental methodologies that have been used to assess the performance of BR algorithms are reviewed, and a methodological framework for the assessment of BR algorithms is presented. Third, we outline the most pressing directions for future research, including the steps required to use BR algorithms in wearable sensors, remote video monitoring, and clinical practice.

Video-Based Physiologic Monitoring During an Acute Hypoxic Challenge: Heart Rate, Respiratory Rate, and Oxygen Saturation

BACKGROUND

The physiologic information contained in the video photoplethysmogram is well documented. However, extracting this information during challenging conditions requires new analysis techniques to capture and process the video image streams to extract clinically useful physiologic parameters. We hypothesized that heart rate, respiratory rate, and oxygen saturation trending can be evaluated accurately from video information during acute hypoxia.

METHODS

Video footage was acquired from multiple desaturation episodes during a porcine model of acute hypoxia using a standard visible light camera. A novel in-house algorithm was used to extract photoplethysmographic cardiac pulse and respiratory information from the video image streams and process it to extract a continuously reported video-based heart rate (HRvid), respiratory rate (RRvid), and oxygen saturation (SvidO2). This information was then compared with HR and oxygen saturation references from commercial pulse oximetry and the known rate of respiration from the ventilator.

RESULTS

Eighty-eight minutes of data were acquired during 16 hypoxic episodes in 8 animals. A linear mixed-effects regression showed excellent responses relative to a nonhypoxic reference signal with slopes of 0.976 (95% confidence interval [CI], 0.973-0.979) for HRvid; 1.135 (95% CI, 1.101-1.168) for RRvid, and 0.913 (95% CI, 0.905-0.920) for video-based oxygen saturation. These results were obtained while maintaining continuous uninterrupted vital sign monitoring for the entire study period.

CONCLUSIONS

Video-based monitoring of HR, RR, and oxygen saturation may be performed with reasonable accuracy during acute hypoxic conditions in an anesthetized porcine hypoxia model using standard visible light camera equipment. However, the study was conducted during relatively low motion. A better understanding of the effect of motion and the effect of ambient light on the video photoplethysmogram may help refine this monitoring technology for use in the clinical environment.

Photoplethysmography respiratory rate monitoring in patients receiving procedural sedation and analgesia for upper gastrointestinal endoscopy

The value of capnography during procedural sedation and analgesia (PSA) for the detection of hypoxaemia during upper gastrointestinal (UGI) endoscopic procedures is limited. Photoplethysmography respiratory rate (RRp) monitoring may provide a useful alternative, but the level of agreement with capnography during PSA is unknown. We therefore investigated the level of agreement between the RRp and capnography-based RR (RRc) during PSA for UGI endoscopy. This study included patients undergoing PSA for UGI endoscopy procedures. Pulse oximetry (SpO2) and RRc were recorded in combination with Nellcor 2.0 (RRp) monitoring (Covidien, USA). Bland-Altman analysis was used to evaluate the level of agreement between RRc and RRp. Episodes of apnoea, defined as no detection of exhaled CO2 for minimal 36 s, and hypoxaemia, defined as an SpO2 < 92 %, were registered. A total of 1054 min of data from 26 patients were analysed. Bland-Altman analysis between the RRc and RRp revealed a bias of 2.25 ± 5.41 breath rate per minute (brpm), with limits of agreement from -8.35 to 12.84 brpm for an RR ≥ 4 brpm. A total of 67 apnoea events were detected. In 21 % of all apnoea events, the patient became hypoxaemic. Hypoxaemia occurred 42 times with a median length of 34 (19-141) s, and was preceded in 34 % of the cases by apnoea and in 64 % by an RRc ≥ 8 brpm. In 81 % of all apnoea events, photoplethysmography registered an RRp ≥ 4 brpm. We found a low level of agreement between capnography and the plethysmography respiratory rate during procedural sedation for UGI endoscopy. Moreover, respiratory rate derived from both the capnogram and photoplethysmogram showed a limited ability to provide warning signs for a hypoxaemic event during the sedation procedure.

Multicenter Study Validating Accuracy of a Continuous Respiratory Rate Measurement Derived From Pulse Oximetry: A Comparison With Capnography

Published ahead of print January 17, 2017.

BACKGROUND:

Intermittent measurement of respiratory rate via observation is routine in many patient care settings. This approach has several inherent limitations that diminish the clinical utility of these measurements because it is intermittent, susceptible to human error, and requires clinical resources. As an alternative, a software application that derives continuous respiratory rate measurement from a standard pulse oximeter has been developed. We sought to determine the performance characteristics of this new technology by comparison with clinician-reviewed capnography waveforms in both healthy subjects and hospitalized patients in a low-acuity care setting.

METHODS:

Two independent observational studies were conducted to validate the performance of the Medtronic Nellcor^TM Respiration Rate Software application. One study enrolled 26 healthy volunteer subjects in a clinical laboratory, and a second multicenter study enrolled 53 hospitalized patients. During a 30-minute study period taking place while participants were breathing spontaneously, pulse oximeter and nasal/oral capnography waveforms were collected. Pulse oximeter waveforms were processed to determine respiratory rate via the Medtronic Nellcor Respiration Rate Software. Capnography waveforms reviewed by a clinician were used to determine the reference respiratory rate.

RESULTS:

A total of 23,243 paired observations between the pulse oximeter-derived respiratory rate and the capnography reference method were collected and examined. The mean reference-based respiratory rate was 15.3 ± 4.3 breaths per minute with a range of 4 to 34 breaths per minute. The Pearson correlation coefficient between the Medtronic Nellcor Respiration Rate Software values and the capnography reference respiratory rate is reported as a linear correlation, R, as 0.92 ± 0.02 (P < .001), whereas Lin’s concordance correlation coefficient indicates an overall agreement of 0.85 ± 0.04 (95% confidence interval [CI] +0.76; +0.93) (healthy volunteers: 0.94 ± 0.02 [95% CI +0.91; +0.97]; hospitalized patients: 0.80 ± 0.06 [95% CI +0.68; +0.92]). The mean bias of the Medtronic Nellcor Respiration Rate Software was 0.18 breaths per minute with a precision (SD) of 1.65 breaths per minute (healthy volunteers: 0.37 ± 0.78 [95% limits of agreement: –1.16; +1.90] breaths per minute; hospitalized patients: 0.07 ± 1.99 [95% limits of agreement: –3.84; +3.97] breaths per minute). The root mean square deviation was 1.35 breaths per minute (healthy volunteers: 0.81; hospitalized patients: 1.60).

CONCLUSIONS:

These data demonstrate the performance of the Medtronic Nellcor Respiration Rate Software in healthy subjects and patients hospitalized in a low-acuity care setting when compared with clinician-reviewed capnography. The observed performance of this technology suggests that it may be a useful adjunct to continuous pulse oximetry monitoring by providing continuous respiratory rate measurements. The potential patient safety benefit of using combined continuous pulse oximetry and respiratory rate monitoring warrants assessment.

Ensemble Empirical Mode Decomposition With Principal Component Analysis: A Novel Approach for Extracting Respiratory Rate and Heart Rate From Photoplethysmographic Signal

The photoplethysmographic (PPG) signal measures the local variations of blood volume in tissues, reflecting the peripheral pulse modulated by cardiac activity, respiration, and other physiological effects. Therefore, PPG can be used to extract the vital cardiorespiratory signals like heart rate (HR), and respiratory rate (RR) and this will reduce the number of sensors connected to the patient's body for recording these vital signs. In this paper, we propose an algorithm based on ensemble empirical mode decomposition with principal component analysis (EEMD-PCA) as a novel approach to estimate HR and RR simultaneously from PPG signal. To examine the performance of the proposed algorithm, we used 310 (from 35 subjects) and 632 (from 42 subjects) epochs of simultaneously recorded electrocardiogram, PPG, and respiratory signal extracted from MIMIC (Physionet ATM data bank) and Capnobase database, respectively. Results of EEMD-PCA-based extraction of HR and RR from PPG signal showed that the median RMS error (1st and 3rd quartiles) obtained in MIMIC data set for RR was 0.89 (0, 1.78) breaths/min, for HR was 0.57 (0.30, 0.71) beats/min and in Capnobase data set it was 2.77 (0.50, 5.9) breaths/min and 0.69 (0.54, 1.10) beats/min for RR and HR, respectively. These results illustrated that the proposed EEMD-PCA approach is more accurate in estimating HR and RR than other existing methods. Efficient and reliable extraction of HR and RR from the pulse oximeter's PPG signal will help patients for monitoring HR and RR with low cost and less discomfort.

Respiratory effort from the photoplethysmogram

The potential for a simple, non-invasive measure of respiratory effort based on the pulse oximeter signal - the photoplethysmogram or 'pleth' - was investigated in a pilot study. Several parameters were developed based on a variety of manifestations of respiratory effort in the signal, including modulation changes in amplitude, baseline, frequency and pulse transit times, as well as distinct baseline signal shifts. Thirteen candidate parameters were investigated using data from healthy volunteers. Each volunteer underwent a series of controlled respiratory effort maneuvers at various set flow resistances and respiratory rates. Six oximeter probes were tested at various body sites. In all, over three thousand pleth-based effort-airway pressure (EP) curves were generated across the various airway constrictions, respiratory efforts, respiratory rates, subjects, probe sites, and the candidate parameters considered. Regression analysis was performed to determine the existence of positive monotonic relationships between the respiratory effort parameters and resulting airway pressures. Six of the candidate parameters investigated exhibited a distinct positive relationship (p<0.001 across all probes tested) with increasing upper airway pressure repeatable across the range of respiratory rates and flow constrictions studied. These were: the three fundamental modulations in amplitude (AM-Effort), baseline (BM-Effort) and respiratory sinus arrhythmia (RSA-Effort); two pulse transit time modulations - one using a pulse oximeter probe and an ECG (P2E-Effort) and the other using two pulse oximeter probes placed at different peripheral body sites (P2-Effort); and baseline shifts in heart rate, (BL-HR-Effort). In conclusion, a clear monotonic relationship was found between several pleth-based parameters and imposed respiratory loadings at the mouth across a range of respiratory rates and flow constrictions. The results suggest that the pleth may provide a measure of changing upper airway dynamics indicative of the effort to breathe.

Robust respiration detection from remote photoplethysmography

Continuous monitoring of respiration is essential for early detection of critical illness. Current methods require sensors attached to the body and/or are not robust to subject motion. Alternative camera-based solutions have been presented using motion vectors and remote photoplethysmography. In this work, we present a non-contact camera-based method to detect respiration, which can operate in both visible and dark lighting conditions by detecting the respiratory-induced colour differences of the skin. We make use of the close similarity between skin colour variations caused by the beating of the heart and those caused by respiration, leading to a much improved signal quality compared to single-channel approaches. Essentially, we propose to find the linear combination of colour channels which suppresses the distortions best in a frequency band including pulse rate, and subsequently we use this same linear combination to extract the respiratory signal in a lower frequency band. Evaluation results obtained from recordings on healthy subjects which perform challenging scenarios, including motion, show that respiration can be accurately detected over the entire range of respiratory frequencies, with a correlation coefficient of 0.96 in visible light and 0.98 in infrared, compared to 0.86 with the best-performing non-contact benchmark algorithm. Furthermore, evaluation on a set of videos recorded in a Neonatal Intensive Care Unit (NICU) shows that this technique looks promising as a future alternative to current contact-sensors showing a correlation coefficient of 0.87.

Toward a Robust Estimation of Respiratory Rate From Pulse Oximeters

Goal: Current methods for estimating respiratory rate (RR) from the photoplethysmogram (PPG) typically fail to distinguish between periods of high- and low-quality input data, and fail to perform well on independent “validation” datasets. The lack of robustness of existing methods directly results in a lack of penetration of such systems into clinical practice. The present work proposes an alternative method to improve the robustness of the estimation of RR from the PPG. Methods: The proposed algorithm is based on the use of multiple autoregressive models of different orders for determining the dominant respiratory frequency in the three respiratory-induced variations (frequency, amplitude, and intensity) derived from the PPG. The algorithm was tested on two different datasets comprising 95 eight-minute PPG recordings (in total) acquired from both children and adults in different clinical settings, and its performance using two window sizes (32 and 64 seconds) was compared with that of existing methods in the literature. Results: The proposed method achieved comparable accuracy to existing methods in the literature, with mean absolute errors (median, 25\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{} $\text {th}$\end{document}–75\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\text {th}$\end{document} percentiles for a window size of 32 seconds) of 1.5 (0.3–3.3) and 4.0 (1.8–5.5) breaths per minute (for each dataset respectively), whilst providing RR estimates for a greater proportion of windows (over 90% of the input data are kept). Conclusion: Increased robustness of RR estimation by the proposed method was demonstrated. Significance: This work demonstrates that the use of large publicly available datasets is essential for improving the robustness of wearable-monitoring algorithms for use in clinical practice.

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.

Detection of artifacts on photoplethysmography signals using random distortion testing

In this work, we describe a novel method based on waveform morphology for detecting artifacts in photoplethysmography (PPG) signals and, thus, improve reliability of PPG. By considering inter-individual and measure condition variability, specific parameters are estimated for each record. We introduce a novel metric for comparing pulses, which is the derivative of the correlation coefficient. Then, we propose a detection method based on Random Distortion Testing (RDT), to perform adaptive threasholding for each record. The results show that the proposed method provides pertinent detection of pulses with artifacts. Tested on 104 PPG records, the mean of sensitivity, specificity and accuracy were 84 ± 16%, 83 ± 12% and 83 ± 8%, respectively.

Unobtrusive sleep state measurements in preterm infants - A review

Sleep is important for the development of preterm infants. During sleep, neural connections are formed and the development of brain regions is triggered. In general, various rudimentary sleep states can be identified in the preterm infant, namely active sleep (AS), quiet sleep (QS) and intermediate sleep (IS). As the infant develops, sleep states change in length and organization, with these changes as important indicators of brain development. As a result, several methods have been deployed to distinguish between the different preterm infant sleep states, among which polysomnography (PSG) is the most frequently used. However, this method is limited by the use of adhesive electrodes or patches that are attached to the body by numerous cables that can disturb sleep. Given the importance of sleep, this review explores more unobtrusive methods that can identify sleep states without disturbing the infant. To this end, after a brief introduction to preterm sleep states, an analysis of the physiological characteristics associated with the different sleep states is provided and various methods of measuring these physiological characteristics are explored. Finally, the advantages and disadvantages of each of these methods are evaluated and recommendations for neonatal sleep monitoring proposed.

An assessment of algorithms to estimate respiratory rate from the electrocardiogram and photoplethysmogram

Over 100 algorithms have been proposed to estimate respiratory rate (RR) from the electrocardiogram (ECG) and photoplethysmogram (PPG). As they have never been compared systematically it is unclear which algorithm performs the best. Our primary aim was to determine how closely algorithms agreed with a gold standard RR measure when operating under ideal conditions. Secondary aims were: (i) to compare algorithm performance with IP, the clinical standard for continuous respiratory rate measurement in spontaneously breathing patients; (ii) to compare algorithm performance when using ECG and PPG; and (iii) to provide a toolbox of algorithms and data to allow future researchers to conduct reproducible comparisons of algorithms. Algorithms were divided into three stages: extraction of respiratory signals, estimation of RR, and fusion of estimates. Several interchangeable techniques were implemented for each stage. Algorithms were assembled using all possible combinations of techniques, many of which were novel. After verification on simulated data, algorithms were tested on data from healthy participants. RRs derived from ECG, PPG and IP were compared to reference RRs obtained using a nasal-oral pressure sensor using the limits of agreement (LOA) technique. 314 algorithms were assessed. Of these, 270 could operate on either ECG or PPG, and 44 on only ECG. The best algorithm had 95% LOAs of  -4.7 to 4.7 bpm and a bias of 0.0 bpm when using the ECG, and  -5.1 to 7.2 bpm and 1.0 bpm when using PPG. IP had 95% LOAs of  -5.6 to 5.2 bpm and a bias of  -0.2 bpm. Four algorithms operating on ECG performed better than IP. All high-performing algorithms consisted of novel combinations of time domain RR estimation and modulation fusion techniques. Algorithms performed better when using ECG than PPG. The toolbox of algorithms and data used in this study are publicly available.

Journal of Clinical Monitoring and Computing 2015 end of year summary: respiration

This paper reviews 17 papers or commentaries published in Journal of Clinical Monitoring and Computing in 2015, within the field of respiration. Papers were published covering monitoring and training of breathing, monitoring of gas exchange, hypoxemia and acid-base, and CO2 monitoring.

Respiratory modulations in the photoplethysmogram (DPOP) as a measure of respiratory effort

DPOP is a measure of the strength of respiratory modulations present in the pulse oximetry photoplethysmogram (pleth) waveform. It has been proposed as a non-invasive parameter for the prediction of the response to volume expansion in hypovolemic patients. The effect of resistive breathing on the DPOP parameter was studied to determine whether it may have an adjunct use as a measure of respiratory effort. Healthy volunteers were tasked to breathe at fixed respiratory rates over a range of airway resistances generated by a flow resistor inserted within a mouthpiece. Changes in respiratory efforts, effected by the subjects and measured as airway pressures at the mouth, were compared to DPOP values derived from a finger pulse oximeter probe. It was found that the increased effort to breathe manifests itself as an associated increase in DPOP. Further, a relationship between DPOP and percent modulation of the pleth waveform was observed. A version of the DPOP algorithm that corrects for low perfusion was implemented which resulted in an improved relationship between DPOP and PPV. Although a limited cohort of seven volunteers was used, the results suggest that DPOP may be useful as a respiratory effort parameter, given that the fluid level of the patient is maintained at a constant level over the period of analysis.

A Review of Wavelet Transform Time-Frequency Methods for NIRS-Based Analysis of Cerebral Autoregulation

Near-infrared spectroscopy (NIRS) has been proposed as a suitable technique for the analysis of cerebral autoregulation as it provides a simpler acquisition methodology and more artifact-free signal. A number of sophisticated wavelet transform methods have recently emerged to quantify the cerebral autoregulation mechanism using NIRS and blood pressure signals. These provide an enhanced partitioning of signal information via the time-frequency plane, which facilitates improved extraction of the components of interest. This area is reviewed, and enhancements to this form of analysis are suggested.

Pulse oximetry-derived respiratory rate in general care floor patients

Respiratory rate is recognized as a clinically important parameter for monitoring respiratory status on the general care floor (GCF). Currently, intermittent manual assessment of respiratory rate is the standard of care on the GCF. This technique has several clinically-relevant shortcomings, including the following: (1) it is not a continuous measurement, (2) it is prone to observer error, and (3) it is inefficient for the clinical staff. We report here on an algorithm designed to meet clinical needs by providing respiratory rate through a standard pulse oximeter. Finger photoplethysmograms were collected from a cohort of 63 GCF patients monitored during free breathing over a 25-min period. These were processed using a novel in-house algorithm based on continuous wavelet-transform technology within an infrastructure incorporating confidence-based averaging and logical decision-making processes. The computed oximeter respiratory rates (RRoxi) were compared to an end-tidal CO2 reference rate (RRETCO2). RRETCO2 ranged from a lowest recorded value of 4.7 breaths per minute (brpm) to a highest value of 32.0 brpm. The mean respiratory rate was 16.3 brpm with standard deviation of 4.7 brpm. Excellent agreement was found between RRoxi and RRETCO2, with a mean difference of -0.48 brpm and standard deviation of 1.77 brpm. These data demonstrate that our novel respiratory rate algorithm is a potentially viable method of monitoring respiratory rate in GCF patients. This technology provides the means to facilitate continuous monitoring of respiratory rate, coupled with arterial oxygen saturation and pulse rate, using a single non-invasive sensor in low acuity settings.

Real-time estimation of respiratory rate from a photoplethysmogram using an adaptive lattice notch filter

Background

Many researchers have attempted to acquire respiratory rate (RR) information from a photoplethysmogram (PPG) because respiration affects the waveform of the PPG. However, most of these methods were difficult to operate in real-time because of their complexity or computational requirements. From these needs, we attempted to develop a method to estimate RR from a PPG with a light computational burden.

Methods

To obtain RR information, we adopt a sequential filtering structure and frequency estimation technique, which extracts a dominant frequency from a given signal. In particular, we used an adaptive lattice notch filter (ALNF) to estimate RR from a PPG along with an additional heart rate that is utilized as an adaptation parameter of our method. Furthermore, we designed a sequential infinite impulse response (IIR) notch filtering system (i.e., harmonic IIR notch filter) to eliminate the cardiac component and its harmonics from the PPG. We compared the proposed method with Burg’s AR modeling method, which is widely used to estimate RR from a PPG, using open-source data and measured data.

Results

By using a statistical test, it was determined that our adaptive lattice-type respiratory rate estimator (ALRE) was significantly more accurate than Burg’s AR model method (p <0.0001). Furthermore, the ALRE’s tracking performance was better than that of Burg’s method, and the variances of its estimates were smaller than those of Burg’s method.

Conclusions

In short, our method showed a better performance than Burg’s AR modeling method for real-time applications.

Detection of motion artifact patterns in photoplethysmographic signals based on time and period domain analysis

The presence of motion artifacts in photoplethysmographic (PPG) signals is one of the major obstacles in the extraction of reliable cardiovascular parameters in continuous monitoring applications. In the current paper we present an algorithm for motion artifact detection based on the analysis of the variations in the time and the period domain characteristics of the PPG signal. The extracted features are ranked using a normalized mutual information feature selection algorithm and the best features are used in a support vector machine classification model to distinguish between clean and corrupted sections of the PPG signal. The proposed method has been tested in healthy and cardiovascular diseased volunteers, considering 11 different motion artifact sources. The results achieved by the current algorithm (sensitivity--SE: 84.3%, specificity--SP: 91.5% and accuracy--ACC: 88.5%) show that the current methodology is able to identify both corrupted and clean PPG sections with high accuracy in both healthy (ACC: 87.5%) and cardiovascular diseases (ACC: 89.5%) context.

Increasing signal processing sophistication in the calculation of the respiratory modulation of the photoplethysmogram (DPOP)

DPOP (∆POP or Delta-POP) is a non-invasive parameter which measures the strength of respiratory modulations present in the pulse oximetry photoplethysmogram (pleth) waveform. It has been proposed as a non-invasive surrogate parameter for pulse pressure variation (PPV) used in the prediction of the response to volume expansion in hypovolemic patients. Many groups have reported on the DPOP parameter and its correlation with PPV using various semi-automated algorithmic implementations. The study reported here demonstrates the performance gains made by adding increasingly sophisticated signal processing components to a fully automated DPOP algorithm. A DPOP algorithm was coded and its performance systematically enhanced through a series of code module alterations and additions. Each algorithm iteration was tested on data from 20 mechanically ventilated OR patients. Correlation coefficients and ROC curve statistics were computed at each stage. For the purposes of the analysis we split the data into a manually selected ‘stable’ region subset of the data containing relatively noise free segments and a ‘global’ set incorporating the whole data record. Performance gains were measured in terms of correlation against PPV measurements in OR patients undergoing controlled mechanical ventilation. Through increasingly advanced pre-processing and post-processing enhancements to the algorithm, the correlation coefficient between DPOP and PPV improved from a baseline value of R = 0.347 to R = 0.852 for the stable data set, and, correspondingly, R = 0.225 to R = 0.728 for the more challenging global data set. Marked gains in algorithm performance are achievable for manually selected stable regions of the signals using relatively simple algorithm enhancements. Significant additional algorithm enhancements, including a correction for low perfusion values, were required before similar gains were realised for the more challenging global data set.

Calculation of the Respiratory Modulation of the Photoplethysmogram (DPOP) Incorporating a Correction for Low Perfusion

DPOP quantifies respiratory modulations in the photoplethysmogram. It has been proposed as a noninvasive surrogate for pulse pressure variation (PPV) used in the prediction of the response to volume expansion in hypovolemic patients. The correlation between DPOP and PPV may degrade due to low perfusion effects. We implemented an automated DPOP algorithm with an optional correction for low perfusion. These two algorithm variants (DPOPa and DPOPb) were tested on data from 20 mechanically ventilated OR patients split into a benign "stable region" subset and a whole record "global set." Strong correlation was found between DPOP and PPV for both algorithms when applied to the stable data set: R = 0.83/0.85 for DPOPa/DPOPb. However, a marked improvement was found when applying the low perfusion correction to the global data set: R = 0.47/0.73 for DPOPa/DPOPb. Sensitivities, Specificities, and AUCs were 0.86, 0.70, and 0.88 for DPOPa/stable region; 0.89, 0.82, and 0.92 for DPOPb/stable region; 0.81, 0.61, and 0.73 for DPOPa/global region; 0.83, 0.76, and 0.86 for DPOPb/global region. An improvement was found in all results across both data sets when using the DPOPb algorithm. Further, DPOPb showed marked improvements, both in terms of its values, and correlation with PPV, for signals exhibiting low percent modulations.

Automated prediction of early blood transfusion and mortality in trauma patients

BACKGROUND

Prediction of blood transfusion needs and mortality for trauma patients in near real time is an unrealized goal. We hypothesized that analysis of pulse oximeter signals could predict blood transfusion and mortality as accurately as conventional vital signs (VSs).

METHODS

Continuous VS data were recorded for direct admission trauma patients with abnormal prehospital shock index (SI = heart rate [HR] / systolic blood pressure) greater than 0.62. Predictions of transfusion during the first 24 hours and in-hospital mortality using logistical regression models were compared with DeLong's method for areas under receiver operating characteristic curves (AUROCs) to determine the optimal combinations of prehospital SI and HR, continuous photoplethysmographic (PPG), oxygen saturation (SpO2), and HR-related features.

RESULTS

We enrolled 556 patients; 37 received blood within 24 hours; 7 received more than 4 U of red blood cells in less than 4 hours or "massive transfusion" (MT); and 9 died. The first 15 minutes of VS signals, including prehospital HR plus continuous PPG, and SpO2 HR signal analysis best predicted transfusion at 1 hour to 3 hours, MT, and mortality (AUROC, 0.83; p < 0.03) and no differently (p = 0.32) from a model including blood pressure. Predictions of transfusion based on the first 15 minutes of data were no different using 30 minutes to 60 minutes of data collection. SI plus PPG and SpO2 signal analysis (AUROC, 0.82) predicted 1-hour to 3-hour transfusion, MT, and mortality no differently from pulse oximeter signals alone.

CONCLUSION

Pulse oximeter features collected in the first 15 minutes of our trauma patient resuscitation cohort, without user input, predicted early MT and mortality in the critical first hours of care better than the currently used VS such as combinations of HR and systolic blood pressure or prehospital SI alone.

LEVEL OF EVIDENCE

Therapeutic/prognostic study, level II.

Exploration of key stakeholders' preferences for pre-hospital physiologic monitoring by emergency rescue services

To gather preferences for novel pre-hospital physiologic monitoring technologies from emergency rescue services. Qualitative semi-structured interviews and focus groups were conducted with three groups from UK Search and Rescue (SAR); (1) Extractors (e.g. SAR teams), (2) Transporters (personnel primarily responsible for casualty transport), and (3) Treaters (e.g. Emergency Department doctors). Three themes were defined; SAR casualty management, novel physiologic monitor potential, and physiologic monitor physical properties. Some SAR groups already employed physiologic monitoring but there was no consensus on which monitor(s) to carry or what to monitor and how frequently. Existing monitors also tended to be bulky and heavy and could be unreliable in an unstable environment or if the casualty was cold. Those performing monitoring tended to have only basic first-aid training, and their workload was often high particularly if there was more than one casualty. The potential benefits of employing a novel monitor were strategic and clinical; an opportunity for transmitting data off-scene in order to facilitate monitoring or generate advice (i.e. telemedicine) was also voiced. A range of more intuitive, physical properties was also raised (e.g. small/compact, lightweight). SAR-specific technology should be simple to operate by those with less medical training, which means that clinical data interpretation and presentation should be carefully considered. It would be beneficial if novel monitors carried out a majority of the interpretation, allowing rescuers to proceed with their priority task of removing the casualty to safety.

Multiparameter respiratory rate estimation from the photoplethysmogram

We present a novel method for estimating respiratory rate in real time from the photoplethysmogram (PPG) obtained from pulse oximetry. Three respiratory-induced variations (frequency, intensity, and amplitude) are extracted from the PPG using the Incremental-Merge Segmentation algorithm. Frequency content of each respiratory-induced variation is analyzed using fast Fourier transforms. The proposed Smart Fusion method then combines the results of the three respiratory-induced variations using a transparent mean calculation. It automatically eliminates estimations considered to be unreliable because of detected presence of artifacts in the PPG or disagreement between the different individual respiratory rate estimations. The algorithm has been tested on data obtained from 29 children and 13 adults. Results show that it is important to combine the three respiratory-induced variations for robust estimation of respiratory rate. The Smart Fusion showed trends of improved estimation (mean root mean square error 3.0 breaths/min) compared to the individual estimation methods (5.8, 6.2, and 3.9 breaths/min). The Smart Fusion algorithm is being implemented in a mobile phone pulse oximeter device to facilitate the diagnosis of severe childhood pneumonia in remote areas.

Implementation of artifact detection in critical care: a methodological review

Artifact detection (AD) techniques minimize the impact of artifacts on physiologic data acquired in critical care units (CCU) by assessing quality of data prior to clinical event detection (CED) and parameter derivation (PD). This methodological review introduces unique taxonomies to synthesize over 80 AD algorithms based on these six themes: 1) CCU; 2) physiologic data source; 3) harvested data; 4) data analysis; 5) clinical evaluation; and 6) clinical implementation. Review results show that most published algorithms: a) are designed for one specific type of CCU; b) are validated on data harvested only from one OEM monitor; c) generate signal quality indicators (SQI) that are not yet formalized for useful integration in clinical workflows; d) operate either in standalone mode or coupled with CED or PD applications; e) are rarely evaluated in real-time; and f) are not implemented in clinical practice. In conclusion, it is recommended that AD algorithms conform to generic input and output interfaces with commonly defined data: 1) type; 2) frequency; 3) length; and 4) SQIs. This shall promote: a) reusability of algorithms across different CCU domains; b) evaluation on different OEM monitor data; c) fair comparison through formalized SQIs; d) meaningful integration with other AD, CED and PD algorithms; and e) real-time implementation in clinical workflows.

Identification of apnea during respiratory monitoring using support vector machine classifier: a pilot study

To determine the use of photoplethysmography (PPG) as a reliable marker for identifying respiratory apnea based on time-frequency features with support vector machine (SVM) classifier. The PPG signals were acquired from 40 healthy subjects with the help of a simple, non-invasive experimental setup under normal and induced apnea conditions. Artifact free segments were selected and baseline and amplitude variabilities were derived from each recording. Frequency spectrum analysis was then applied to study the power distribution in the low frequency (0.04-0.15 Hz) and high frequency (0.15-0.40 Hz) bands as a result of respiratory pattern changes. Support vector machine (SVM) learning algorithm was used to distinguish between the normal and apnea waveforms using different time-frequency features. The algorithm was trained and tested (780 and 500 samples respectively) and all the simulations were carried out using linear kernel function. Classification accuracy of 97.22 % was obtained for the combination of power ratio and reflection index features using SVM classifier. The pilot study indicates that PPG can be used as a cost effective diagnostic tool for detecting respiratory apnea using a simple, robust and non-invasive experimental setup. The ease of application and conclusive results has proved that such a system can be further developed for use in real-time monitoring under critical care conditions.

Magnetic resonance thermal imaging combined with SMASH navigators in the presence of motion

This study develops and tests an MR thermometry method combined with SMASH navigators in phantom experiments mimicking human liver motion with the purpose of detecting and correcting motion artifacts in thermal MR images. Experimental data were acquired on a 3T MRI scanner. Motion artifacts of mobile phantoms mimicking human liver motion were detected and corrected using the SMASH navigators and then MR temperature maps were obtained using a proton resonant frequency (PRF) shift method with complex image subtraction. Temperature acquired by MR thermal imaging was compared to that measured via thermocouples. MR thermal imaging combined with the SMASH navigator technique resulted in accurate temperature maps of the mobile phantoms compared to temperatures measured using the thermocouples. The differences between the obtained and measured temperatures varied from 8.2°C to 14.2°C and 2.2°C to 4.9°C without and with motion correction, respectively. Motion correction improved the temperature acquired by MR thermal imaging by >55%. The combination of the MR thermal imaging and SMASH navigator technique will enable monitoring and controlling heat distribution and temperature change in tissues during thermal therapies and will be a very important tool for cancer treatment in mobile organs.

PACS number: 87.57.‐s