Abstract

Measuring blood oxygen saturation (SpO2) is crucial in a triage process for identifying patients with respiratory distress or shock, since low SpO2 levels indicate compromised hemostability and the need for priority treatment. This paper explores the use of wearable mechanical deflection sensors based on laser-induced graphene (LIG) for SpO2 estimation. The LIG sensors are attached to a subject’s chest for real-time monitoring of respiratory signals. We have developed a novel database of the respiratory signals, with corresponding SpO2 values ranging from 86% to 100%. The database is used to develop an artificial neural network model for SpO2 estimation. The neural network performance is promising, with regression metrics mean squared error = 0.184, mean absolute error = 0.301, root mean squared error = 0.429, and R-squared = 0.804. The use of mechanical respiration sensors in combination with neural networks in biosensing opens new possibilities for noninvasive SpO2 monitoring and other innovative applications.

1. Introduction

Real-time health monitoring has the potential to profoundly alter the strategies employed in the prevention and management of diseases, particularly in triage situations with a high number of casualties, by facilitating prompt identification of the most severely injured individuals. In particular, the ability to monitor vital parameters, such as heartbeat, respiration, blood oxygen saturation (SpO2), and blood pressure in real-time can save lives in those critical situations. Timely detection of anomalies in these parameters can provide early warning of medical emergencies, enabling healthcare professionals to take swift action and prevent unwanted outcomes [1].

Measuring SpO2 is an important part of the triage process, as it is used to establish patient hemostability and to prioritize patient treatment based on the severity of each patient’s condition [2]. The main reason for measuring this vital parameter is detecting victims or patients with severe respiratory distress or shock: patients in shock may have decreased tissue perfusion and low SpO2 levels. The common way of measuring SpO2 is by utilizing a photoplethysmography (PPG) sensor, i.e., by analyzing the light absorption characteristics of oxygenated and deoxygenated hemoglobin. However, there are some drawbacks to using PPG for SpO2 measurements. One limitation is that PPG-based SpO2 measurements may not be as accurate as those obtained from arterial blood gas analysis or pulse oximetry using dedicated sensors. PPG signals can also be affected by motion artifacts, ambient light interference, and variations in skin pigmentation or thickness [3].

Wearable mechanical deflection sensors have been shown to be an effective and reliable method for real-time monitoring of respiration, with potential applications in sleep apnea monitoring, exercise physiology, and respiratory disease management [46]. Laser-induced graphene (LIG), an emerging material recently used for mechanical deflection sensors, is an excellent candidate for wearable respiration monitoring [710]. LIG is a type of graphene that is created by irradiating a polymer or other organic material with a laser [11]. The laser energy causes the material to carbonize and transform into a graphene-like structure. Graphene, a single layer of carbon atoms arranged in a hexagonal lattice, has unique properties that make it attractive for sensor applications [12]. Its high electrical conductivity, large surface area, and biocompatibility have led to research exploring the use of graphene-based sensors for various biomedical purposes [13, 14]. In our prior publication [15], we established a correlation between LIG-obtained signals and heartbeat parameters using the HeartPy toolkit implemented in Python [16].

In this paper, we investigate the use of wearable mechanical deflection sensors made of LIG, used as respiration monitors, for SpO2 estimation as an alternative for oximeters with PPG signals, with a vision of using such sensors as components of integrated wearable patch-like multisensor devices, since there is a relationship between the magnitude of respiratory peaks and the oxygen concentration [17]. Our main goal is to create an integrated patch-like sensor capable of delivering the four crucial hemostability parameters (heart rate, respiratory rate, blood pressure, and SpO2). These parameters are vital for monitoring the medical condition of patients or injured individuals in emergency situations [18]. This patch-like sensor serves as the sole device required for continuous monitoring of the biosignals, eliminating the need for at least three separate devices (oximeter, blood pressure device, ECG/PPG sensors).

The LIG sensors employed in this paper detect changes in mechanical deflection, such as chest or abdominal movements, which are used to measure respiration rate and volume. The sensors are noninvasive and can be easily integrated into wearable devices such as belts or patches, making them suitable for continuous monitoring of respiratory parameters, as well as other vital parameters [15]. An artificial neural network (ANN) is used to extract SpO2 from respiration measurements, making use of the fact that oxygenation and respiration are closely interconnected physiological processes. A direct link between real-time respiration signals and SpO2, such as the one that we make, has not been demonstrated before. The use of LIG in biosensors for the real-time SpO2 estimation and prediction opens the door to novel applications that are not available with other types of sensing devices. This paper also elaborates on the creation of a novel database consisting solely of LIG-obtained respiratory signals with corresponding SpO2 references in the range of 86%–100%.

In the past, researchers have demonstrated the use of graphene in prototype flexible photodetectors for measuring heart rate, arterial oxygen saturation (SpO2), and respiratory rate [19]. Such prototypes are of good use for wearable respiration monitoring using optical technology; however, they suffer from the same disadvantages as traditional PPG technology. Similarly, there is a plethora of research focused on SpO2 estimation using machine learning and deep learning models, as well as calculations based on PPG signals [20, 21], but none based on real-time respiration signals. Our approach overcomes some of the mentioned drawbacks of using PPG signals, as well as the need for preprocessing of PPG signals when used in ML/ANN/DNN models [22].

The paper is organized as follows: in Section 2, we describe the production of the LIG sensor and material characterization (Section 2.1). We also describe the novel LIG database (Section 2.2) and the used ANN model for SpO2 estimation (Section 2.3). In Section 3, we elaborate on the obtained results, presenting a discussion in Section 4. Section 5 presents concluding remarks and several ideas for future work.

2. Materials and Methods

2.1. Graphene Sensor Production, Characterization, and Respiration Monitoring

LIG was produced by scanning a CO2 laser beam across the surface of polyimide tape, as illustrated in Figure 1 [23]. The laser used was DBK FL-350, with a maximum power of 60 W. The power was set to 18%, the scanning speed was 400 mm/s, and the scanning resolution was 800 DPI.

The devices were formed by laser-writing LIG in the shape of rectangles with dimensions 1 × 3 cm, as in Figure 2. The graphene was transferred to double-sided adhesive medical tape (Duplomed 8411, Lohmann, GmbH, Neuwied, Germany). Raman spectra of the samples were recorded with a DXR Raman microscope (Thermo Fisher Scientific, Waltham, MA, USA). The samples were excited with a diode laser at a wavelength of 532 nm and a power of 10 mW, focused on a 2.1 μm spot on the surface. Spectra were obtained as averages of three measurements from different positions on the sample (10 exposures, 10 s each per position). The recorded Raman spectra were treated to correct the fluorescence background using a 5th-order polynomial baseline correction method built-in the Omnic software (OMNIC for Dispersive Raman 9.2.41.). SEM with EDS was done with a PhenomProX scanning electron microscope (Phenom, Thermo Fisher Scientific, Waltham, MA, USA). XRD analysis was carried out on a D8 Advance Eco diffractometer with Cu K radiation ( = 0.15,406 nm, Bruker, Germany). Electrical contacts were made to the LIG by attaching conductive copper tapes at the device’s ends. Wires were soldered to copper tapes and connected to the measuring device. Such contacting provides ohmic contact to the graphene. The sensor was attached to the body by adhesion of the medical tape to the skin in the chest area.

The experimental setup consisted of a one-channel source-measure unit (SMU) Keithley 2450, a desktop computer that controls the measuring device and acquires data using a Windows application written in C#, a graphene sensor attached to the body using medical tape and connected with cables to the SMU Keithley 2450, and a commercial Onyx II oximeter (Nonin Medical, Plymouth, MN, USA) placed on a finger of a volunteer subject (Figure 3). The specified precision of the oximeter of +/−2% is sufficient to evaluate the concept of applying machine learning in combination with mechanical deflection sensors to predict SpO2 from respiration data.

Measurements were performed in the constant current mode with the current set to 0.1 mA, and the voltage was measured over a period of 30 s. The placement of the graphene sensor on the torso allows for the capture of respiratory signals, which provide valuable information regarding the breathing pattern (Figure 4). The measurement shown in Figure 4 starts with exhalation, followed by a sequence of seven breaths, as labeled on the graph. Although we measure voltage while keeping the current fixed, it is more intuitive to plot the resistance than the voltage. Hence, in Figure 4, we convert the measured voltage to resistance using the simple linear relation V = IR, and we plot resistance over time during breathing.

2.2. LIG Database

Respiration data were collected from three healthy volunteer subjects, as our primal interest is the triage process of healthy injured victims of mass victim incidents. The number of subjects relates to already published studies—one and three subjects, respectively [24, 25]. During our measuring, the subjects sat in a chair and breathed normally while wearing the LIG sensor. Simultaneously, reference SpO2 measurements were obtained using the oximeter. The signals were recorded and conducted with approval from the Ethical Committee. We created a novel database consisting of 89 measurements of respiratory signals. The majority of measurements in the novel database exhibit normal oxygen saturation (SpO2 95%), but we also recorded some values below 95%, even down to 86%, as depicted in Figure 5. To record data with SpO2 less than 90%, i.e., the critical level [22], we applied the breath-holding technique [26, 27]. The subjects held their breath for 30, 45 s, 1, or 2 min, resuming breathing after that. The first 30–40 s of breathing after holding breath were recorded. Including SpO2 values lower than 95% in the database is important for applications in clinical and casualty situations in which patients may exhibit low SpO2. We have qualitatively observed that, although intense motion may cause signal artifacts, the artifacts do not survive the postprocessing treatment. Moreover, these sensors are expected to be used in cases when wounded or sick individuals are lying still and not moving.

As some commercial SpO2 sensors provide readings every 10 s [28], we aimed at similarly fast measurements by utilizing 20–30 s of the LIG sensor signals in our analysis. By utilizing an approach that allows for signals of different temporal durations, we aim to expand on the range of opportunities for applying our neural network algorithm to respiration signals.

Unlike PPG signals typically used for SpO2 estimation, the LIG-obtained signals did not require preprocessing such as normalization, Butterworth band-pass filtering, or Hampel decision filtering. The LIG-obtained respiratory signals were used raw and unprocessed for the estimation of SpO2.

In the process of generating an ANN learning model for SpO2 prediction, a crucial step is to generate a list of features from the respiratory signals. These features serve as input vectors for the model and play a vital role in training the network to identify significant patterns related to SpO2 estimation. By extracting relevant features from the respiratory signals, the ANN can learn to recognize patterns and establish relationships between the input signals and SpO2 levels. This feature extraction step is essential for capturing the relevant information present in the respiratory signals and facilitating accurate predictions of SpO2. These extracted features can then be analyzed and utilized to characterize the respiratory pattern [7, 8].

Using empirical analysis, a set of 10 common features was selected to describe the LIG respiratory signals. These specific features were chosen based on their significance in capturing essential characteristics of the respiratory signals. The NumPy [29] and SciPy [30] Python libraries were used to support the features calculation.

Table 1 presents the elaborated respiratory features in a structured way, together with a short explanation. The comprehensive definitions and explanations for all features are available in the Appendix section.

2.3. ANN Model

This research employed a supervised learning approach by building an ANN to estimate the SpO2 value as a regression problem. The ANN is designed to output a single value representing the predicted SpO2 value, while the input is a matrix containing the 10 respiratory features extracted from the input signal.

The dataset was divided into a training set, consisting of 70% of the data (62 measurements), and a testing set containing the remaining 30% (27 measurements). Prior to splitting the dataset, a standardization process was applied to ensure consistent scaling and facilitate effective training and evaluation of the model.

The implementation of the ANN involved the use of the Keras API written in Python, which runs on top of the TensorFlow machine learning platform. The architecture of the employed ANN in this study is illustrated in Figure 6.

The Keras library’s sequential module was employed to build a sequence of ANN layers that are stacked in consecutive order. To specify the number of neurons, the Dense Keras module is used to define each layer. As depicted in Figure 6, the learning neural network is fully connected, comprising several hidden layers with a particular number of neurons and one output layer with only one neuron that predicts the SpO2 value.

To perform calculations within each neuron, we employed the rectified linear unit (ReLU) function as the activation function, which is the most widely used activation function [31]. The ReLU function produces an output of zero if the input value is less than zero; otherwise, it is equal to the provided input value. To compute the loss, we used mean squared error (MSE) as it is the most commonly used loss function for regression [32].

To determine the optimal accuracy, we conducted a tuning process on the ANN model, exploring different combinations of the “epoch” and “batch_size” values. Specifically, we used grid search cross-validation to test various values for these hyperparameters. The “batch_size” parameter refers to the number of training examples used in a single forward/backward pass, while ’epochs’ indicates the number of times the learning algorithm runs through the entire training dataset.

To assess the accuracy of the model, we used standard regression metrics: mean absolute error (MAE), MSE, root mean squared error (RMSE), root mean log squared error (RMLSE), and R-squared (R2).

MAE represents the average absolute difference between the predicted and actual values, providing a measure of the model’s overall accuracy. MSE implements the average squared difference between the predicted and actual values, emphasizing larger errors. RMSE is the square root of MSE, which provides a more interpretable metric in the original scale of the target variable. R2, also known as the coefficient of determination, measures the proportion of the variance in the target variable that can be explained by the model. A higher R2 value indicates a better fit of the model to the data. By considering these metrics, we can assess the model’s accuracy, precision, and ability to explain the variability in the target variable, aiding in the evaluation and comparison of different ML/ANN models.

3. Results

The Raman spectrum of LIG (Figure 7(a)) exhibits characteristic features of graphene. Prominent D and G bands appear at 1,338 and 1,590 cm, respectively. The G peak indicates the presence of sp carbon atoms, while the D peak signals defects in graphene. The 2D peak at 2,700 cm is attributed to second-order zone boundary phonons [33, 34].

The XRD pattern of LIG (Figure 7(b)) shows a significant peak at , suggesting the presence of a graphitic phase with an inter-planar distance of 3.4 Å corresponding to the (002) plane [11]. The interplanar distance is calculated using Bragg’s law for the first-order diffraction, X-ray wavelength of 1.54 Å, and an incident beam angle of . A characteristic peak at indicates the production of graphite oxide, while the peak at is formed due to the formation of reduced graphene oxide, caused by oxygen contamination on the polyimide films during the laser induction process [35].

To examine both the morphology and elemental composition of LIG, SEM-EDS analysis was conducted (Figure 7(c)). The produced LIG film displays a distinctive foam-like appearance with porous structures, a result of the rapid release of gaseous products during the laser-induction process [36]. The sample was further used to assess the elemental composition of LIG (Figure 7(d)). EDS analysis reveals the presence of carbon, oxygen, silicon, and chlorine, while the atomic percentage of carbon equals 95.25%. These findings indicate that negligible oxidation occurred during the graphene formation process.

Figure 8 illustrates the training loss during the training process. The decreasing trend observed in the training loss curve indicates an improvement in the model’s performance over time. This suggests that the model is effectively learning and adapting to the training data without overfitting. This outcome is desirable as it indicates the model’s ability to generalize well to unseen data and make accurate predictions beyond the training set. The decreasing trend in the training loss curve provides confidence in the model’s learning process and its potential effectiveness in predicting SpO2 values. We tested the optimal number of epochs by performing experiments with different numbers of epochs. The validation loss ceases to improve as the number of epochs increases beyond 90, which was the final number of epochs that we used.

Table 2 presents the actual and predicted SpO2 values from the testing set containing 27 measurements.

Table 3 shows results from the model’s performance evaluation metrics in predicting SpO2 from respiratory signals.

The ANN model’s performance in predicting SpO2 from respiratory signals yields promising results. The MSE of 0.184 suggests that, on average, the model’s predictions deviate from the actual SpO2 values by a relatively small squared difference. The MAE value of 0.301 indicates that the model’s average absolute difference between predicted and actual SpO2 values is also reasonably low. The RMSE value of 0.429, derived from the MSE, provides a more interpretable metric in the original scale of SpO2, showing that the model’s predictions align well with the true values. Furthermore, the R2 value of 0.804 demonstrates that approximately 80.4% of the variance in SpO2 can be explained by the model, indicating a good fit to the underlying patterns in the data. These results collectively demonstrate the effectiveness and accuracy of the presented ANN learning model in predicting SpO2 from LIG-obtained respiratory signals.

4. Discussion

Our research aimed to explore the feasibility of utilizing LIG sensors as an alternative SpO2 estimator by leveraging the nature of LIG-obtained signals as respiratory signals. The impact of the demonstrated results is twofold. First, the results demonstrate that LIG sensors can be used to measure respiration in real time, which opens doors for the use of low-profile wearable LIG-based patches for medical and wellness applications. Second, we have shown that temporal respiration data obtained from signals of chest expansion and contraction during breathing can be used as input to an ANN that outputs SpO2 with high precision and correctness.

The underlying physiological relationship between SpO2 and respiration action is rooted in the process of oxygen transfer and circulation within the body. When we breathe, oxygen enters the respiratory system and is transported to the bloodstream, where it binds to hemoglobin in red blood cells. This oxygenated blood is then distributed throughout the body, ensuring the delivery of oxygen to various tissues and organs [37]. Given this connection, by analyzing respiratory signals captured by the LIG sensors, we aimed to investigate their potential to serve as indicators for estimators of SpO2 levels.

Earlier work has shown that simple RR measurements do not correlate well with SpO2 in patients [38]. Nevertheless, here we have shown that the utilization of the chosen 10 respiratory features as input vectors in our ANN model overcomes this hurdle, with excellent matching to reference measurements of SpO2.

Although our study has used real-time respiratory data acquired with LIG sensors, we believe that the model could be applicable to different real-time respiratory data obtained from various types of sensors. Given that the number and type of available sensors are growing at a fast pace with the development of new technology [39], our model will certainly find numerous applications in wearable and bedside devices.

Direct comparisons with other studies that estimate SpO2 by ML/ANN models using LIG-obtained signals as input are not possible due to the lack of existing literature and databases with LIG-obtained respiratory signals and SpO2 references. Future work may encompass analysis of publicly available databases containing respiratory signals and SpO2 measured by other means, such as radio frequency signals [40]. For example, the National Institute of Health (USA) has compiled databases with various health parameters, including respiration signals and SpO2, measured during sleep [41].

One notable advantage of our approach is the flexibility in the length of the respiration signal as an input vector for the ANN. Also, the feature extraction process maintains the integrity of the respiratory information. Furthermore, our research revealed that preprocessing the LIG-obtained respiratory signals is unnecessary, in contrast to the preprocessing required for PPG signals used in SpO2 estimation. This finding holds great significance for practical applications with time-critical implications.

The compactness of the presented ANN model is a significant advantage since it provides quick and timely SpO2 level estimations. The model can be deployed in emergency situations as an integral component of a telemedicine system [18].

In terms of future work, there is potential to enhance the novel database by including more data, particularly focusing on SpO2 values below 90%, which may be obtained in a hospital setting from unhealthy patients. As indicated in Table 2, the largest disparity between the actual and predicted SpO2 values is observed for the value of 89%. By expanding the dataset to incorporate a wider range of SpO2 values, including those below 90%, the performance and accuracy of the ANN model can be further refined.

5. Conclusion

This paper explores the use of wearable mechanical deflection sensors made of LIG for estimating SpO2 as an alternative to PPG-based oximeters. These sensors measure mechanical deflection, such as chest or abdominal movements, to estimate respiration rate and volume. The LIG sensors are affixed to the subject’s chest to enable real-time monitoring of respiratory signals. An ANN model is built to estimate SpO2 from the respiratory measurements, leveraging the close relationship between oxygenation and respiration. To facilitate SpO2 estimation, we have created a unique database comprising respiratory signals and corresponding SpO2 values ranging from 86% to 100%. The neural network’s performance demonstrates promise, as evidenced by regression metrics such as MSE of 0.184, MAE of 0.301, RMSE of 0.429, and R2 of 0.804. The integration of mechanical respiration sensors with ANNs in biosensing presents exciting possibilities and the potential of using LIG sensors in SpO2 estimation with high precision and accuracy.

Appendix

The first four features from the 10 features described in Section 2.2 represent “Amplitude” features. In respiratory signals, the term “amplitude” refers to the magnitude or intensity of the variations observed in the signal waveform. Amplitude represents the difference between the baseline or resting state and the peak values of the respiratory signal. The following four amplitude features were identified:

Peak amplitude of the sinusoidal respiratory waveform is the maximum positive deviation from its zero reference level;

Trough amplitude is the lowest part of the respiratory signal. It represents the minimum value reached by the signal during the exhalation phase;

Mean amplitude refers to the average magnitude or intensity of the variations observed in the signal waveform over a certain period of time. It represents the average difference between the baseline or resting state and the peak values of the respiratory signal;

Standard deviation, in this case, measures how far the signal fluctuates from the mean.

The following two features belong to the time domain in signal processing:

Breath duration is the duration of one breath in the measured sequence;

Respiratory rate, RR, is calculated as the number of breaths per minute, usually obtained by frequency analysis or peak detection methods. We calculate RR by the latter.

The next two are variability features:

Respiratory rate variabilities (RRVs) are the mean variations within the respiratory rhythm. Normal breathing exhibits a relatively constant rate and tidal volume that together constitute normal respiratory rhythm. However, variations within this respiratory rhythm are labeled as RRV [42]. RRV is a useful predictor of the deterioration of a patient’s health. RRV Mean is calculated as the average value of respiratory rate variability across a specific time period;

RRV Std is the standard deviation or spread of the respiratory rate intervals around the mean value. A higher RRV Std indicates greater variability and a wider range of interval values [43].

The last two features are shape features. In respiratory signals, these values refer to characteristics of measurements that describe the overall shape or waveform pattern of the respiratory signal:

Rise time in respiratory signals is defined as the time during which airway pressure builds up from a baseline toward a preset maximum value. A short rise time will allow instantaneous delivery of flow at the start of a breath cycle, resulting in an immediate rise in pressure to a preset level [44];

Fall time is defined as the time that airway pressure takes to return from a maximum value toward a baseline or resting state.

Data Availability

The data used in this study are available upon request. Please contact Bojana Koteska at [email protected] for access to the dataset.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the NATO Science for Peace and Security Program under project SP4LIFE, number G5825. This work was supported in part by the Ministry of Science, Technological Development and Innovation of the Republic of Serbia, grant number 451-03-47/2023-01/200026.