Estimating respiratory and heart rates from the correntropy spectral density of the photoplethysmogram

Sleep-phasic heart rate variability predicts stress severity: Building a machine learning-based stress prediction model

We propose a novel approach for predicting stress severity by measuring sleep phasic heart rate variability (HRV) using a smart device. This device can potentially be applied for stress self-screening in large populations. Using a Holter electrocardiogram (ECG) and a Huawei smart device, we conducted 24-h dual recordings of 159 medical workers working regular shifts. Based on photoplethysmography (PPG) and accelerometer signals acquired by the Huawei smart device, we sorted episodes of cyclic alternating pattern (CAP; unstable sleep), non-cyclic alternating pattern (NCAP; stable sleep), wakefulness, and rapid eye movement (REM) sleep based on cardiopulmonary coupling (CPC) algorithms. We further calculated the HRV indices during NCAP, CAP and REM sleep episodes using both the Holter ECG and smart-device PPG signals. We later developed a machine learning model to predict stress severity based only on the smart device data obtained from the participants along with a clinical evaluation of emotion and stress conditions. Sleep phasic HRV indices predict individual stress severity with better performance in CAP or REM sleep than in NCAP. Using the smart device data only, the optimal machine learning-based stress prediction model exhibited accuracy of 80.3 %, sensitivity 87.2 %, and 63.9 % for specificity. Sleep phasic heart rate variability can be accurately evaluated using a smart device and subsequently can be used for stress predication.

A Self-Calibrated Single Wavelength Biosensor for Measuring Oxygen Saturation

Traditional methods for measuring blood oxygen use multiple wavelengths, which produce an intrinsic error due to ratiometric measurements. These methods assume that the absorption changes with the wavelength, but in fact the scattering changes as well and cannot be neglected. We found that if one measures in a specific angle around a cylindrical tissue, called the iso-pathlength (IPL) point, the reemitted light intensity is unaffected by the tissue's scattering. Therefore, the absorption can be isolated from the scattering, which allows the extraction of the subject's oxygen saturation. In this work, we designed an optical biosensor for reading the light intensity reemitted from the tissue, using a single light source and multiple photodetectors (PDs), with one of them in the IPL point's location. Using this bio-device, we developed a methodology to extract the arterial oxygen saturation using a single wavelength light source. We proved this method is not dependent on the light source and is applicable to different measurement locations on the body, with an error of 0.5%. Moreover, we tested thirty-eight males and females with the biosensor under normal conditions. Finally, we show the results of measuring subjects in a hypoxic chamber that simulates extreme conditions with low oxygen.

Breathing Rate Estimation from Head-Worn Photoplethysmography Sensor Data Using Machine Learning

Breathing rate is considered one of the fundamental vital signs and a highly informative indicator of physiological state. Given that the monitoring of heart activity is less complex than the monitoring of breathing, a variety of algorithms have been developed to estimate breathing activity from heart activity. However, estimating breathing rate from heart activity outside of laboratory conditions is still a challenge. The challenge is even greater when new wearable devices with novel sensor placements are being used. In this paper, we present a novel algorithm for breathing rate estimation from photoplethysmography (PPG) data acquired from a head-worn virtual reality mask equipped with a PPG sensor placed on the forehead of a subject. The algorithm is based on advanced signal processing and machine learning techniques and includes a novel quality assessment and motion artifacts removal procedure. The proposed algorithm is evaluated and compared to existing approaches from the related work using two separate datasets that contains data from a total of 37 subjects overall. Numerous experiments show that the proposed algorithm outperforms the compared algorithms, achieving a mean absolute error of 1.38 breaths per minute and a Pearson's correlation coefficient of 0.86. These results indicate that reliable estimation of breathing rate is possible based on PPG data acquired from a head-worn device.

An integrated framework for evaluation on typical ECG-derived respiration waveform extraction and respiration

OBJECTIVE

ECG-derived respiration (EDR) methods have been developed during the past decades to obtain respiration-relevant information. However, it is still necessary to compare the performance of these methods under uniform conditions for reasonable application.

APPROACH

In this paper, the performance of 10 feature-based EDR methods was evaluated comprehensively on three aspects: sampling rate, noise, and window length. The Fantasia database was used in this study, as it contained ECG signals and simultaneously measured respiration signals. The performance was quantified by two parameters: waveform correlation and breathing rate (BR) errors.

MAIN RESULTS

The BR errors of AM_area, AM_QR, AM_R were all below 2 beats per minute (bpm) when the sampling rate was above 150 Hz, while they decreased sharply by about 60% when the sampling rate was below 150 Hz. FM_RR presented stable performance with an error below 2 bpm at different sampling rates. The effect of noise was obviously found in amplitude-based EDR methods, with the maximum decreased by about 40% in waveform correlation. For all EDR methods, significant increase of BR errors occurred with the window shorting from 32 s to 16 s in the frequency-based technique. In addition, about 30%-40% of the window cannot obtain the BR error, calculated based on the time-based technique, within an 8 s window.

SIGNIFICANCE

We proposed a comprehensive and integrated evaluation on typical ECG-derived respiration waveform extraction and respiration rate calculation, providing references for algorithm selection based on different requirements.

Modulation Model of the Photoplethysmography Signal for Vital Sign Extraction

This paper introduces an amplitude and frequency modulation (AM-FM) model to characterize the photoplethysmography (PPG) signal. The model indicates that the PPG signal spectrum contains one dominant frequency component - the heart rate (HR), which is guarded by two weaker frequency components on both sides; the distance from the dominant component to the guard components represents the respiratory rate (RR). Based on this model, an efficient algorithm is proposed to estimate both HR and RR by searching for the dominant frequency component and two guard components. The proposed method is performed in the frequency domain to estimate RR, which is more robust to additive noise than the prior art based on temporal features. Experiments were conducted on two types of PPG signals collected with a contact sensor (an oximeter) and a contactless visible imaging sensor (a color camera), respectively. The PPG signal from the contactless sensor is much noisier than the signal from the contact sensor. The experimental results demonstrate the effectiveness of the proposed algorithm, including under relatively noisy scenarios.

Non-contact breathing rate monitoring in newborns: A review

The neonatal period - the first 4 weeks of life - is the most critical time for a child's survival. Breathing rate is a vital indicator of the health condition and requires continuous monitoring in case of sickness or preterm birth. Breathing movements can be counted by contact and non-contact methods. In the case of newborn infants, the non-contact breathing rate monitoring need is high, as a contact-based approach may interfere while providing care and is subject to interference by non-breathing movements. This review article delivers a factual summary, and describes the methods and processing involved in non-contact based breathing rate monitoring. The article also provides the advantages, limitations, and clinical applications of these methods. Additionally, signal processing, feasibility, and future direction of different non-contact neonatal breathing rate monitoring are discussed.

Magnetomotive Ultrasound Imaging Systems: Basic Principles and First Applications

This review discusses magnetomotive ultrasound, which is an emerging technique that uses superparamagnetic iron oxide nanoparticles as a contrast agent. The key advantage of using nanoparticle-based contrast agents is their ability to reach extravascular targets, whereas commercial contrast agents for ultrasound comprise microbubbles confined to the blood stream. This also extends possibilities for molecular imaging, where the contrast agent is labeled with specific targeting molecules (e.g., antibodies) so that pathologic tissue may be visualized directly. The principle of action is that an external time-varying magnetic field acts to displace the nanoparticles lodged in tissue and thereby their immediate surrounding. This movement is then detected with ultrasound using frequency- or time-domain analysis of echo data. As a contrast agent already approved for magnetic resonance imaging (MRI) by the US Food and Drug Administration, there is a shorter path to clinical translation, although safety studies of magnetomotion are necessary, especially if particle design is altered to affect biodistribution or signal strength. The external modulated magnetic field may be generated by electromagnets, permanent magnets, or a combination of the two. The induced nanoparticle motion may also reveal mechanical material properties of tissue, healthy or diseased, one of several interesting potential future aspects of the technique.

PPG Derived Respiratory Rate Estimation in Daily living Conditions

Respiratory rate (RR) derived from photoplethysmogram (PPG) during daily activities can be corrupted due to movement and other artefacts. We have investigated the use of ensemble empirical mode decomposition (EEMD) based smart fusion approach for improving the RR extraction from PPG. PPG was recorded while subjects performed five different activities: sitting, standing, climbing and descending stairs, walking, and running. RR was obtained using EEMD and smart fusion. The median absolute error (AE) of the proposed method is superior, median AE = 3.05 (range 3.01 to 3.18) breath/min in estimating RR during five different activities. Therefore, the proposed method can be implemented for overcoming the artefact problems when recording continuous RR monitoring during activities of daily living.

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.

FMCW Laser Fuze Multiple Scattering Model and Accurate Fixed-Distance Algorithm in a Smoke Environment

In a smoke environment, suspended particles can scatter and absorb laser photons, making target echo signals extremely weak and difficult to extract and identify, which causes obvious difficulty in fixed-distance of laser fuze. In this paper, the multiple scattering model of frequency-modulated-continuous-wave (FMCW) laser fuze in a smoke environment was established. This model simulates multi-path propagation and multiple scattering of photons. At the same time, we use the correntropy spectral density (CSD) algorithm for accurate fixed-distance of FMCW laser fuze. The absolute error of distance does not exceed 0.15 m in smoke interference environment.

Motion Assessment for Accelerometric and Heart Rate Cycling Data Analysis

Motion analysis is an important topic in the monitoring of physical activities and recognition of neurological disorders. The present paper is devoted to motion assessment using accelerometers inside mobile phones located at selected body positions and the records of changes in the heart rate during cycling, under different body loads. Acquired data include 1293 signal segments recorded by the mobile phone and the Garmin device for uphill and downhill cycling. The proposed method is based upon digital processing of the heart rate and the mean power in different frequency bands of accelerometric data. The classification of the resulting features was performed by the support vector machine, Bayesian methods, k-nearest neighbor method, and neural networks. The proposed criterion is then used to find the best positions for the sensors with the highest discrimination abilities. The results suggest the sensors be positioned on the spine for the classification of uphill and downhill cycling, yielding an accuracy of 96.5% and a cross-validation error of 0.04 evaluated by a two-layer neural network system for features based on the mean power in the frequency bands 〈 3 , 8 〉 and 〈 8 , 15 〉 Hz. This paper shows the possibility of increasing this accuracy to 98.3% by the use of more features and the influence of appropriate sensor positioning for motion monitoring and classification.

A non-contact infection screening system using medical radar and Linux-embedded FPGA: Implementation and preliminary validation

Objectives

In this study, an infection screening system was developed to detect patients suffering from infectious diseases. In addition, the system was also designed to deal with the variability in age and gender, which would affect the accuracy of the detection. Furthermore, to enable a low-cost, non-contact and embedded system, multiple vital signs from a medical radar were measured and all algorithms were implemented on a Field Programmable Gate Array, named PYNQ-Z1.

Methods

The system consisted of two main stages: digital signal processing and data classification. In the former stage, Butterworth filters, with flexible cut-off frequencies depending on age and gender, and a time-domain peak detection algorithm were deployed to compute three vital signs, namely heart rate, respiratory rate, and standard deviation of heart beat-to-beat interval. For the classification problem, two machine learning models, Support Vector Machine and Quadratic Discriminant Analysis, were implemented.

Results

The Student's t-test showed that our proposed digital signal processing algorithms coped well with the variability of human cases in age and gender. Meanwhile, the f₁-score of roughly 98.0% represented the high sensitivity and specificity of our proposed machine learning methods.

Conclusion

This study outlines the implementation of an infection screening system, which achieved competent performance. The system might be beneficial for fast screening of infected patients at public health centers in underdeveloped areas, where people have little access to healthcare.

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.

Assessment of Signal Processing Methods for Measuring the Respiratory Rate in the Neonatal Intensive Care Unit

Knowledge of the pathological instabilities in the breathing pattern can provide valuable insights into the cardiorespiratory status of the critically-ill infant as well as their maturation level. This paper is concerned with the measurement of respiratory rate in premature infants. We compare the rates estimated from the chest impedance pneumogram, the ECG-derived respiratory rhythms, and the PPG-derived respiratory rhythms against those measured in the reference standard of breath detection provided by attending clinical staff during 165 manual breath counts. We demonstrate that accurate RR estimates can be produced from all sources for RR in the 40-80 bpm (breaths per min) range. We also conclude that the use of indirect methods based on the ECG or the PPG poses a fundamental challenge in this population due to their poor behavior at fast breathing rates (upward of 80 bpm).

Estimation of respiratory rate from motion contaminated photoplethysmography signals incorporating accelerometry

Estimation of respiratory rate (RR) from photoplethysmography (PPG) signals has important applications in the healthcare sector, from assisting doctors onwards to monitoring patients in their own homes. The problem is still very challenging, particularly during the motion for large segments of data, where results from different methods often do not agree. The authors aim to propose a new technique which performs motion reduction from PPG signals with the help of simultaneous acceleration signals where the PPG and accelerometer sensors need to be embedded in the same sensor unit. This method also reconstructs motion corrupted PPG signals in the Hilbert domain. An auto-regressive (AR) based technique has been used to estimate the RR from reconstructed PPGs. The proposed method has provided promising results for the estimation of RRs and their variations from PPG signals corrupted with motion artefact. The proposed platform is able to contribute to continuous in-hospital and home-based monitoring of patients using PPG signals under various conditions such as rest and motion states.

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.

Step-to-Step Ankle Inversion/Eversion Torque Modulation Can Reduce Effort Associated with Balance

Below-knee amputation is associated with higher energy expenditure during walking, partially due to difficulty maintaining balance. We previously found that once-per-step push-off work control can reduce balance-related effort, both in simulation and in experiments with human participants. Simulations also suggested that changing ankle inversion/eversion torque on each step, in response to changes in body state, could assist with balance. In this study, we investigated the effects of ankle inversion/eversion torque modulation on balance-related effort among amputees (N = 5) using a multi-actuated ankle-foot prosthesis emulator. In stabilizing conditions, changes in ankle inversion/eversion torque were applied so as to counteract deviations in side-to-side center-of-mass acceleration at the moment of intact-limb toe off; higher acceleration toward the prosthetic limb resulted in a corrective ankle inversion torque during the ensuing stance phase. Destabilizing controllers had the opposite effect, and a zero gain controller made no changes to the nominal inversion/eversion torque. To separate the balance-related effects of step-to-step control from the potential effects of changes in average mechanics, average ankle inversion/eversion torque and prosthesis work were held constant across conditions. High-gain stabilizing control lowered metabolic cost by 13% compared to the zero gain controller (p = 0.05). We then investigated individual responses to subject-specific stabilizing controllers following an enforced exploration period. Four of five participants experienced reduced metabolic rate compared to the zero gain controller (−15, −14, −11, −6, and +4%) an average reduction of 9% (p = 0.05). Average prosthesis mechanics were unchanged across all conditions, suggesting that improvements in energy economy might have come from changes in step-to-step corrections related to balance. Step-to-step modulation of inversion/eversion torque could be used in new, active ankle-foot prostheses to reduce walking effort associated with maintaining balance.

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.

How Nonlinear-Type Time-Frequency Analysis Can Help in Sensing Instantaneous Heart Rate and Instantaneous Respiratory Rate from Photoplethysmography in a Reliable Way

Despite the population of the noninvasive, economic, comfortable, and easy-to-install photoplethysmography (PPG), it is still lacking a mathematically rigorous and stable algorithm which is able to simultaneously extract from a single-channel PPG signal the instantaneous heart rate (IHR) and the instantaneous respiratory rate (IRR). In this paper, a novel algorithm called deppG is provided to tackle this challenge. deppG is composed of two theoretically solid nonlinear-type time-frequency analyses techniques, the de-shape short time Fourier transform and the synchrosqueezing transform, which allows us to extract the instantaneous physiological information from the PPG signal in a reliable way. To test its performance, in addition to validating the algorithm by a simulated signal and discussing the meaning of "instantaneous," the algorithm is applied to two publicly available batch databases, the Capnobase and the ICASSP 2015 signal processing cup. The former contains PPG signals relative to spontaneous or controlled breathing in static patients, and the latter is made up of PPG signals collected from subjects doing intense physical activities. The accuracies of the estimated IHR and IRR are compared with the ones obtained by other methods, and represent the state-of-the-art in this field of research. The results suggest the potential of deppG to extract instantaneous physiological information from a signal acquired from widely available wearable devices, even when a subject carries out intense physical activities.

Human-in-the-loop Bayesian optimization of wearable device parameters

The increasing capabilities of exoskeletons and powered prosthetics for walking assistance have paved the way for more sophisticated and individualized control strategies. In response to this opportunity, recent work on human-in-the-loop optimization has considered the problem of automatically tuning control parameters based on realtime physiological measurements. However, the common use of metabolic cost as a performance metric creates significant experimental challenges due to its long measurement times and low signal-to-noise ratio. We evaluate the use of Bayesian optimization-a family of sample-efficient, noise-tolerant, and global optimization methods-for quickly identifying near-optimal control parameters. To manage experimental complexity and provide comparisons against related work, we consider the task of minimizing metabolic cost by optimizing walking step frequencies in unaided human subjects. Compared to an existing approach based on gradient descent, Bayesian optimization identified a near-optimal step frequency with a faster time to convergence (12 minutes, p < 0.01), smaller inter-subject variability in convergence time (± 2 minutes, p < 0.01), and lower overall energy expenditure (p < 0.01).

Video Synchronization With Bit-Rate Signals and Correntropy Function

We propose an approach for the synchronization of video streams using correntropy. Essentially, the time offset is calculated on the basis of the instantaneous transfer rates of the video streams that are extracted in the form of a univariate signal known as variable bit-rate (VBR). The state-of-the-art approach uses a window segmentation strategy that is based on consensual zero-mean normalized cross-correlation (ZNCC). This strategy has an elevated computational complexity, making its application to synchronizing online data streaming difficult. Hence, our proposal uses a different window strategy that, together with the correntropy function, allows the synchronization to be performed for online applications. This provides equivalent synchronization scores with a rapid offset determination as the streams come into the system. The efficiency of our approach has been verified through experiments that demonstrate its viability with values that are as precise as those obtained by ZNCC. The proposed approach scored 81 % in time reference classification against the equivalent 81 % of the state-of-the-art approach, requiring much less computational power.

An EEMD-PCA approach to extract heart rate, respiratory rate and respiratory activity from PPG signal

The pulse oximeter's photoplethysmographic (PPG) signals, measure 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), respiratory rate (RR) and respiratory activity (RA) and this will reduce the number of sensors connected to the patient's body for recording 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, RR and RA simultaneously from PPG signal. To examine the performance of the proposed algorithm, we used 45 epochs of PPG, electrocardiogram (ECG) and respiratory signal extracted from the MIMIC database (Physionet ATM data bank). The ECG and capnograph based respiratory signal were used as the ground truth and several metrics such as magnitude squared coherence (MSC), correlation coefficients (CC) and root mean square (RMS) error were used to compare the performance of EEMD-PCA algorithm with most of the existing methods in the literature. Results of EEMD-PCA based extraction of HR, RR and RA from PPG signal showed that the median RMS error (quartiles) obtained for RR was 0 (0, 0.89) breaths/min, for HR was 0.62 (0.56, 0.66) beats/min and for RA the average value of MSC and CC was 0.95 and 0.89 respectively. These results illustrated that the proposed EEMD-PCA approach is more accurate in estimating HR, RR and RA than other existing methods.

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.

Photoplethysmography sampling frequency: pilot assessment of how low can we go to analyze pulse rate variability with reliability?

Pulse rate variability (PRV) analysis appears as the first alternative to heart rate variability analysis for wearable devices; however, there is a constraint on computational load and energy consumption for the limited system resources available to the devices. Considering that adjustment of the sampling frequency is one of the strategies for reducing computational load and power consumption, this study aimed to investigate the influence of sampling frequency (f _s) on PRV analysis and to find the minimum sampling frequency while maintaining reliability. We generated 5000, 2500, 1000, 500, 250, 100, 50, 25, 20, 15, 10, 5 Hz down-sampled photoplethysmograms from 10 kHz-sampled PPGs and derived time- and frequency-domain variables of the PRV. These included AVNN, SDNN, SDSD, RMSSD, NN50, pNN50, total power, VLF, LF, HF, LF/HF, nLF and nHF for each down-sampled signal. Derived variables were compared with heart rate variability of the 10 kHz-sampled electrocardiograms, and then statistically investigated using one-way ANOVA test and Bland-Altman analysis. As a result, significant differences (P  <  0.05) were found for SDNN, SDSD, RMSSD, NN50, pNN50, TP, HF, LF/HF, nLF and nHF, but not for AVNN, VLF and LF. Based on the post hoc tests, it was found that the NN50 and pNN50, SDSD and RMSSD, LF/HF and nHF, SDNN, TP and nLF analysis had significant differences at f _s  ⩽  20 Hz, f _s  ⩽  15 Hz, f _s  ⩽10 Hz; f _s  =  5 Hz, respectively. In other words, a significant difference was not seen for any variable if the f _s was greater than 25 Hz. Consequently, our pilot study suggests that analysis of variability in the time and frequency domain from pulse rate obtained through PPG may be potentially as reliable as that derived from the analysis of the electrocardiogram, provided that f _s  ⩾25 Hz sampling frequency is used.

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.

Determination of saturation, heart rate, and respiratory rate at forearm using a Nellcor™ forehead SpO2-saturation sensor

Alterations in arterial blood oxygen saturation, heart rate (HR), and respiratory rate (RR) are strongly associated with intra-hospital cardiac arrests and resuscitations. A wireless, easy-to-use, and comfortable method for monitoring these important clinical signs would be highly useful. We investigated whether the Nellcor™ OxiMask MAX-FAST forehead sensor could provide data for vital sign measurements when located at the distal forearm instead of its intended location at the forehead to provide improved comfortability and easy placement. In a prospective setting, we recruited 30 patients undergoing surgery requiring postoperative care. At the postoperative care unit, patients were monitored for two hours using a standard patient monitor and with a study device equipped with a Nellcor™ Forehead SpO₂ sensor. The readings were electronically recorded and compared in post hoc analysis using Bland-Altman plots, Spearman's correlation, and root-mean-square error (RMSE). Bland-Altman plot showed that saturation (SpO₂) differed by a mean of -0.2 % points (SD, 4.6), with a patient-weighted Spearman's correlation (r) of 0.142, and an RMSE of 4.2 points. For HR measurements, the mean difference was 0.6 bpm (SD, 2.5), r = 0.997, and RMSE = 1.8. For RR, the mean difference was -0.5 1/min (4.1), r = 0.586, and RMSE = 4.0. The SpO₂ readings showed a low mean difference, but also a low correlation and high RMSE, indicating that the Nellcor™ saturation sensor cannot reliably assess oxygen saturation at the forearm when compared to finger PPG measurements.

EMG-Derived Respiration Signal Using the Fixed Sample Entropy during an Inspiratory Load Protocol

Extracting clinical information from one single measurement represents a step forward in the assessment of the respiratory muscle function. This attracting idea entails the reduction of the instrumentation and fosters to develop new medical integrated technologies. We present the use of the fixed sample entropy (fSampEn) as a more direct method to non-invasively derive the breathing activity from the diaphragm electromyographic (EMGdi) signal, and thus to extract the respiratory rate, an important vital sign which is cumbersome and time-consuming to be measured by clinicians. fSampEn is a method to evaluate the EMGdi activity that is less sensitive to the cardiac activity (ECG) and its application has proven to be useful to evaluate the load of the respiratory muscles. The behavior of the proposed method was tested in signals from two subjects that performed an inspiratory load protocol, which consists of increments in the inspiratory mouth pressure (P mouth). Two respiratory signals were derived and compared to the P mouth signal: the ECG-derived respiration (EDR) signal from the lead-I configuration, and the EMG-derived respiration (EMGDR) signal by applying the fSampEn method over the EMGdi signal. The similitude and the lag between signals were calculated through the cross-correlation between each derived respiratory signal and the P mouth. The EMGDR signal showed higher correlation and lower lag values (≥ 0.91 and ≤ 0.70 s, respectively) than the EDR signal (≥ 0.83 and ≤ 0.99 s, respectively). Additionally, the respiratory rate was estimated with the P mouth, EDR and EMGDR signals showing very similar values. The results from this preliminary work suggest that the fSampEn method can be used to derive the respiration waveform from the respiratory muscle electrical activity.

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.

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.

An Ethnographic Observational Study to Evaluate and Optimize the Use of Respiratory Acoustic Monitoring in Children Receiving Postoperative Opioid Infusions

BACKGROUND

Respiratory depression in children receiving postoperative opioid infusions is a significant risk because of the interindividual variability in analgesic requirement. Detection of respiratory depression (or apnea) in these children may be improved with the introduction of automated acoustic respiratory rate (RR) monitoring. However, early detection of adverse events must be balanced with the risk of alarm fatigue. Our objective was to evaluate the use of acoustic RR monitoring in children receiving opioid infusions on a postsurgical ward and identify the causes of false alarm and optimal alarm thresholds.

METHODS

A video ethnographic study was performed using an observational, mixed methods approach. After surgery, an acoustic RR sensor was placed on the participant's neck and attached to a Rad87 monitor. The monitor was networked with paging for alarms. Vital signs data and paging notification logs were obtained from the central monitoring system. Webcam videos of the participant, infusion pump, and Rad87 monitor were recorded, stored on a secure server, and subsequently analyzed by 2 research nurses to identify the cause of the alarm, response, and effectiveness. Alarms occurring within a 90-second window were grouped into a single-alarm response opportunity.

RESULTS

Data from 49 patients (30 females) with median age 14 (range, 4.4-18.8) years were analyzed. The 896 bedside vital sign threshold alarms resulted in 160 alarm response opportunities (44 low RR, 74 high RR, and 42 low SpO2). In 141 periods (88% of total), for which video was available, 65% of alarms were deemed effective (followed by an alarm-related action within 10 minutes). Nurses were the sole responders in 55% of effective alarms and the patient or parent in 20%. Episodes of desaturation (SpO2 < 90%) were observed in 9 patients: At the time of the SpO2 paging trigger, the RR was >10 bpm in 6 of 9 patients. Based on all RR samples observed, the default alarm thresholds, to serve as a starting point for each patient, would be a low RR of 6 (>10 years of age) and 10 (4-9 years of age).

CONCLUSIONS

In this study, the use of RR monitoring did not improve the detection of respiratory depression. An RR threshold, which would have been predictive of desaturations, would have resulted in an unacceptably high false alarm rate. Future research using a combination of variables (e.g., SpO2 and RR), or the measurement of tidal volumes, may be needed to improve patient safety in the postoperative ward.

Estimation of respiratory rate from photoplethysmographic imaging videos compared to pulse oximetry

We present a study evaluating two respiratory rate estimation algorithms using videos obtained from placing a finger on the camera lens of a mobile phone. The two algorithms, based on Smart Fusion and empirical mode decomposition (EMD), consist of previously developed signal processing methods to detect features and extract respiratory induced variations in photoplethysmographic signals to estimate respiratory rate. With custom-built software on an Android phone, photoplethysmographic imaging videos were recorded from 19 healthy adults while breathing spontaneously at respiratory rates between 6 to 32 breaths/min. Signals from two pulse oximeters were simultaneously recorded to compare the algorithms' performance using mobile phone data and clinical data. Capnometry was recorded to obtain reference respiratory rates. Two hundred seventy-two recordings were analyzed. The Smart Fusion algorithm reported 39 recordings with insufficient respiratory information from the photoplethysmographic imaging data. Of the 232 remaining recordings, a root mean square error (RMSE) of 6 breaths/min was obtained. The RMSE for the pulse oximeter data was lower at 2.3 breaths/min. RMSE for the EMD method was higher throughout all data sources as, unlike the Smart Fusion, the EMD method did not screen for inconsistent results. The study showed that it is feasible to estimate respiratory rates by placing a finger on a mobile phone camera, but that it becomes increasingly challenging at respiratory rates greater than 20 breaths/min, independent of data source or algorithm tested.

Respiratory rate assessment from photoplethysmographic imaging

We present a study investigating the suitability of a respiratory rate estimation algorithm applied to photoplethysmographic imaging on a mobile phone. The algorithm consists of a cascade of previously developed signal processing methods to detect features and extract respiratory induced variations in photoplethysmogram signals to estimate respiratory rate. With custom-built software on an Android phone (Camera Oximeter), contact photoplethysmographic imaging videos were recorded using the integrated camera from 19 healthy adults breathing spontaneously at respiratory rates between 6 and 40 breaths/min. Capnometry was simultaneously recorded to obtain reference respiratory rates. Two hundred and ninety-eight Camera Oximeter recordings were available for analysis. The algorithm detected 22 recordings with poor photoplethysmogram quality and 46 recordings with insufficient respiratory information. Of the 232 remaining recordings, a root mean square error of 5.9 breaths/min and a median absolute error of 2.3 breaths/min was obtained. The study showed that it is feasible to estimate respiratory rates by placing a finger on a mobile phone camera, but that it becomes increasingly challenging at respiratory rates higher than 20 breaths/min.

Development of a screening tool for sleep disordered breathing in children using the phone Oximeter™

Background

Sleep disordered breathing (SDB) can lead to daytime sleepiness, growth failure and developmental delay in children. Polysomnography (PSG), the gold standard to diagnose SDB, is a highly resource-intensive test, confined to the sleep laboratory.

Aim

To combine the blood oxygen saturation (SpO₂) characterization and cardiac modulation, quantified by pulse rate variability (PRV), to identify children with SDB using the Phone Oximeter, a device integrating a pulse oximeter with a smartphone.

Methods

Following ethics approval and informed consent, 160 children referred to British Columbia Children's Hospital for overnight PSG were recruited. A second pulse oximeter sensor applied to the finger adjacent to the one used for standard PSG was attached to the Phone Oximeter to record overnight pulse oximetry (SpO₂ and photoplethysmogram (PPG)) alongside the PSG.

Results

We studied 146 children through the analysis of the SpO₂ pattern, and PRV as an estimate of heart rate variability calculated from the PPG. SpO₂ variability and SpO₂ spectral power at low frequency, was significantly higher in children with SDB due to the modulation provoked by airway obstruction during sleep (p-value ). PRV analysis reflected a significant augmentation of sympathetic activity provoked by intermittent hypoxia in SDB children. A linear classifier was trained with the most discriminating features to identify children with SDB. The classifier was validated with internal and external cross-validation, providing a high negative predictive value (92.6%) and a good balance between sensitivity (88.4%) and specificity (83.6%). Combining SpO₂ and PRV analysis improved the classification performance, providing an area under the receiver operating characteristic curve of 88%, beyond the 82% achieved using SpO₂ analysis alone.

Conclusions

These results demonstrate that the implementation of this algorithm in the Phone Oximeter will provide an improved portable, at-home screening tool, with the capability of monitoring patients over multiple nights.

Separating arterial and venous-related components of photoplethysmographic signals for accurate extraction of oxygen saturation and respiratory rate

We propose an algorithm for separating arterial and venous-related signals using second-order statistics of red and infrared signals in a blind source separation technique. The separated arterial signal is used to compute accurate arterial oxygen saturation. We have also introduced an algorithm for extracting the respiratory pattern from the extracted venous-related signal. In addition to real-time monitoring, respiratory rate is also extracted. Our experimental results from multiple subjects show that the proposed separation technique is extremely useful for extracting accurate arterial oxygen saturation and respiratory rate. Specifically, the breathing rate is extracted with average root mean square deviation of 1.89 and average mean difference of -0.69.