Next Article in Journal
Trajectory Clustering-Based Anomaly Detection in Indoor Human Movement
Previous Article in Journal
High-Precision 3D-DIC Measurement Method Based on Improved Forward Newton Iteration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An In-Ear PPG-Based Blood Glucose Monitor: A Proof-of-Concept Study

Department of Electrical and Electronic Engineering, Imperial College London, London SW7 2AZ, UK
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(6), 3319; https://doi.org/10.3390/s23063319
Submission received: 2 February 2023 / Revised: 12 March 2023 / Accepted: 17 March 2023 / Published: 21 March 2023
(This article belongs to the Section Wearables)

Abstract

:
Monitoring diabetes saves lives. To this end, we introduce a novel, unobtrusive, and readily deployable in-ear device for the continuous and non-invasive measurement of blood glucose levels (BGLs). The device is equipped with a low-cost commercially available pulse oximeter whose infrared wavelength (880 nm) is used for the acquisition of photoplethysmography (PPG). For rigor, we considered a full range of diabetic conditions (non-diabetic, pre-diabetic, type I diabetic, and type II diabetic). Recordings spanned nine different days, starting in the morning while fasting, up to a minimum of a two-hour period after eating a carbohydrate-rich breakfast. The BGLs from PPG were estimated using a suite of regression-based machine learning models, which were trained on characteristic features of PPG cycles pertaining to high and low BGLs. The analysis shows that, as desired, an average of 82% of the BGLs estimated from PPG lie in region A of the Clarke error grid (CEG) plot, with 100% of the estimated BGLs in the clinically acceptable CEG regions A and B. These results demonstrate the potential of the ear canal as a site for non-invasive blood glucose monitoring.

1. Introduction

Diabetes is a chronic illness that affects 537 million adults around the world and this number is expected to rise to 783 million people by 2045 [1]. People suffering from diabetes experience abnormal blood glucose levels (BGL), mainly owing to the absence or unregulated production of insulin. For example, in the case of type II diabetes, which accounts for 90% of the cases, the body does not produce enough insulin or is resistant to it, hence causing prolonged high levels of blood glucose [2]. When left unmanaged and untreated, diabetes causes serious health problems and complications that can affect the eyes, heart, kidneys, and blood vessels. Diabetes may even have life-threatening consequences [3], with 6.7 million people worldwide losing their lives to diabetes in 2021 [4]. However, this chronic illness can be cured or delayed when it is in its pre-diabetic stage, and its complications prevented under the correct diet, improved lifestyle, and constant BGL monitoring, a subject of this study.
To minimize the complications and adverse health effects of diabetes, it is crucial for people with diabetes (PWD) to monitor their BGLs several times per day (4–10 times), typically before and after meals, and especially if they are under insulin treatment [5]. Currently available out-of-clinic blood glucose monitoring (BGM) devices include invasive and minimally invasive ones. Invasive BGM involves pricking the finger to extract a small amount of capillary blood using a lancet (small needle). The blood is then collected on a strip which is inserted into a portable machine that measures and displays the BG level [6]. This method is not only painful, but also carries the risk of infection, and does not provide a continuous means for BGM. In addition, it is unsuitable for a wide range of populations, the elderly for example, who might lack fine control of their hands and fingers or have poor eyesight which prevents them from performing all the required steps of the measurement [7].
Minimally invasive devices, on the other hand, do provide a more continuous form of BGM. However, these devices typically involve the temporary implantation of a micro-needle that senses BGL from the interstitial fluid. The main problems with such devices include drifts and lags between the BGL from the interstitial fluid versus true BGL, in addition to the high running cost associated with an implanted needle that needs to be disposed of every few days [8].
With the ever-growing number of PWD and the shortcomings of current BGM systems, there is an urgent need for readily wearable, non-invasive devices for continuous BGM. In recent years, several such non-invasive methods have been investigated; these either rest upon the sampling of bodily fluids, such as urine, tears, and sweat, or are based on shining a light of a certain wavelength through the skin [9]. While each method has its advantages and disadvantages [10], optical methods, especially those pertaining to near-infrared (NIR) spectroscopy, has attracted much interest. The NIR region of the electromagnetic spectrum lies between 700–2500 nm. It penetrates the skin in the mm range, depending on the wavelength [11], and has been adopted for the development of continuous non-invasive BGM devices at various body locations [12,13,14,15,16,17,18]. These devices mainly derive the glucose level from features extracted from photoplethysmography (PPG) waveforms, which can be acquired from any part of the skin that has blood vessels [19]—most commonly from the fingers, wrist, and earlobes. The PPG sensor used typically comprises light-emitting diodes (LED) of a certain wavelength(s) and a photodiode that receives the reflected light and is used to measure changes in blood volume in the tissue. It typically involves shining a light onto the skin and measuring the light intensity that is reflected back. This measurement is then used to calculate the amount of blood in the tissue and can provide information about heart rate, blood pressure, and other physiological parameters.
Recent advancements in wearable technology for vital signs monitoring have given rise to in-ear wearable devices, the so-called hearables [20,21]. Given the privileged position of the head on the human body, the ear provides a convenient and stable site for physiological measurement, while the ear canal also acts as an insulator of external electrical noise [21]. Several studies have used the ear canal for general physiological sensing which includes electroencephalography (EEG) [22,23,24], electrocardiography (ECG) [25], and PPG [26,27,28], and have reported promising results when compared to clinical devices.
This work introduces the proof of concept for a continuously wearable in-ear non-invasive glucose monitoring system that uses a single NIR LED (880 nm) to acquire the PPG signal. For rigor, we considered a whole range of diabetic conditions—a non-diabetic control, pre-diabetic, type I, and type II diabetic subject with the considered BGL range varying from 63–345 mg/dL. The blood glucose levels were estimated based on 18 discriminating features, which were employed by machine learning algorithms to predict the BGL values from the corresponding PPG waveforms. The input variable selection (IVS) across all subjects indicates that the most discriminative features are related to the energy of the PPG cycle. As desired, all of the predicted BGLs fell in regions A and B of the Clarke error (CEG) plots—the regions with no adverse effects on clinical outcomes. To the best of our knowledge, the proposed system is the first to utilize the ear canal as a site of non-invasive blood glucose measurement through NIR spectroscopy and opens up avenues for unobtrusive 24/7 blood glucose monitoring in the community.

2. Materials and Experimental Design

2.1. Working Principle

When the NIR light is emitted through the blood vessels in the ear canal, it interacts with the medium and is reflected, as shown in Figure 1. Its intensity is reduced due to both absorption and scattering [17] so that pulsed-periodic signals are recorded, known as a photoplethysmograph (PPG). The principle of PPG is that blood with a higher concentration of glucose will have different absorbance than blood with lower glucose levels. This will affect the shape of the PPG signal which can then be used to indirectly estimate the blood glucose level through a mapping onto a set of features that are fed into the machine learning algorithms [12,13,14,15,17,18,29].

2.2. Hardware Design

In this study, we used the infrared (880 nm) LED of the MAX30101 digital PPG chip by Maxim Integrated (San Jose, CA, USA), which has been used in a number of studies based on in-ear PPG acquisition [26,27,28]. The PPG chip was positioned on a thin rectangular printed circuit board with decoupling capacitors and was then embedded onto a viscoelastic foam earpiece, and attached to a silicon-based ear hook as shown in Figure 2. For user comfort and sensor stability, the ear hook was available in three different sizes; small, medium, and large. The output of the PPG chip was then connected to a data acquisition board with a Bluetooth-enabled microcontroller (Nordic nRF52840), Nordic Semiconductor, Trondheim, Norway which transferred the data in real time to a PC where it was stored and processed. The collection of reference values of BGLs was performed by a commercially available glucometer (Accu-Check Performa Nano), whose accuracy meets the DIN EN ISO 15197:2013 requirements of a measuring interval; that is, 10–600 mg/dL.

2.3. Experimental Design

Four subjects, aged between 30 and 59 years, were recruited for this study, with one non-diabetic female, one pre-diabetic male, one type I diabetic male, and one type II diabetic male (Table 1). The participants were given detailed instructions and explanations about the recording protocol and were asked to attend the recording sessions in the morning while still fasting. Each participant attended the recordings over the course of 9 days. The BGLs were measured using the gold standard glucometer, while simultaneously the in-ear PPG signal was recorded for two minutes (starting within one minute of the glucometer measurement), and each recording was repeated twice. The BGL and PPG were first measured when a subject was fasting, then a varied and carbohydrate-rich (where permissible) breakfast was offered. Afterward, the BGL and PPG signals were measured at various intervals and at least up to 2 h after breakfast, as depicted in Figure 3. The total number of two-minute recordings acquired was 244, and the glucose levels varied between 65 and 345 mg/dL. It is important to note that during each recording day, the subject was asked to keep the sensor in their ear to minimize the variability between recordings and examine the comfort and ergonomic aspects of the proposed device. In addition, the lancets and test strips were safely disposed of after every finger prick and the subjects were asked to wash their hands with lukewarm water before and after each recording. All subjects gave a written informed consent for inclusion before they participated in the study. The recordings were performed under Imperial College’s ethics committee approval JRCO 20IC6414.

3. Signal Processing/Methods

3.1. Pre-Processing of PPG Waveforms

The recorded raw PPG signals were filtered through a 10th-order Butterworth low-pass filter with a cut-off frequency of 0.01 Hz, to acquire the DC component of PPG. For the AC component of the PPG signal, the raw data were filtered through an 8th-order Butterworth band-pass filter between 0.3 and 10 Hz. The signals were upsampled using linear interpolation to better locate the peaks and troughs. The lower envelope of the filtered signal was then detected using spline interpolation between the local minima and was subtracted to remove the baseline wander, so as to enable the identification of the peaks and troughs. Figure 4 shows a sample recording with the pre-processing steps.

3.2. PPG Cycle Analysis and Feature Extraction

In order to understand the variation of the PPG signal with the glucose level in order to extract features that best capture this relationship, we performed an analysis of the PPG cycles across the glucose levels. The PPG signal was divided into windows centered around the peak of each cycle, with each window containing a single cycle. The PPG cycles were selected for further analysis using a matched filter—a template cycle created by taking the mean of the middle ten cycles of each recording. After each PPG cycle was correlated with the template cycle, the cycles with a correlation value of r < 0.90 were deemed to be too noisy and removed from the analysis, as shown in Figure 5. In this way, 10–20% of the cycles were removed from recordings for each subject. The remaining valid cycles were then divided into those pertaining to low- and high-glucose cycles. For each subject, the average glucose level measured during fasting was used as a threshold to determine periods of low and high glucose levels. Low glucose cycles then correspond to the glucose levels below the threshold and high glucose cycles represent the glucose levels above a subject’s threshold.
The ensemble average PPG cycle was computed from high and low glucose cycles and compared between the two groups. Significant differences in the PPG cycle between the low and high glucose groups were detected by applying a two-sample t-test at each point in the cycle. Figure 6 shows that for all the subjects, the PPG cycles were found to have higher amplitudes and narrower peaks during low glucose periods. Based on this observation and the available literature [12,13,14,15,16,17,18,29,30,31], we created features that capture this change in the dynamics at different glucose levels, as shown in Figure 7. In addition to these features, the energy profile of each cycle was calculated using the Teager–Kaiser energy operator (TKEO), a feature commonly used to quantify the instantaneous energy profile in periodic components of signals [15,16,17,18]. Table 2 summarises all features extracted from every valid PPG cycle.

3.3. Glucose Prediction

Based on the set of features in Table 2, a suite of machine learning models was trained on each subject individually. Subject-specific training was necessary for this pilot study because each subject belongs to a different BG condition category. This was achieved using the first seven days as the training data (with a five-fold cross-validation) and the remaining two days as unseen test data. The models employed include linear regression, support vector machines (SVM), Gaussian process regression, kernel approximation models, ensemble models, and neural networks (NNs). The models were trained in the MATLAB programming environment, using machine learning and deep learning toolboxes [32]. The ensemble regression trees model resulted in the lowest root mean squared error (RMSE) value and was chosen for the analysis of the unseen data. To test the clinical significance of the test data, the Clarke error grid (CEG) was used [33], along with the mean absolute relative difference (MARD), the correlation coefficient r, and the RMSE.

4. Results

The gold standard glucose levels of the recorded subjects, shown in Figure 8, were consistent with their health conditions. For the non-diabetic subject (S1) the glucose ranges were within the normal limits, typically 90 mg/dL while fasting, and with a maximum of 142 mg/dL typically half an hour after the meal, after which the BGL dropped back to around 90 mg/dL. Subject 2 (S2), a pre-diabetic, had a higher fasting glucose level of around 100 mg/dL and a maximum glucose value of 194 mg/dL. High BGL remained constant over a longer period of time, typically returning to the fasting level within about two hours from breakfast. The type I diabetic subject (S3) had hypoglycemic BGLs. This subject was using insulin to regulate the BGL, typically right before the meal, after an estimate of the number of carbohydrates in food. For this type of diabetes, as evident from the data, the glucose levels tend to be quite unpredictable; for example, they can be low both before and after the mean, or equally, they can be high both before and after the meal. The highest BGL can be seen in the type II diabetic subject (S4) where the fasting glucose was around 200 mg/dL and the maximum BGL after the meal was around 354 mg/dL. It is important to note that this subject does not inject any insulin, hence the constantly elevated BGL.
The PPG signals obtained from each subject using the in-ear PPG sensor showed clear systolic peaks, with an average of 80% of the cycles (no. of valid cycles in Table 3) having a high correlation (r > 0.9) with the template cycle of the recording in hand, as shown in Figure 5. The highest number of discarded cycles was for subject S3; this is likely due to the subject moving their mouth, grinding their teeth, yawning, coughing, or due to sensor displacement. The clean cycles were split between the training and test sets, with the first seven days used to train the model and the last two days for testing. The ten most significant features which yielded statistically significant differences in high versus low blood glucose PPG cycles are presented in Figure 9. Observe that the highest percentage differences are found for S1 and S2, while the lowest percentage difference can be observed in S4. Features calculated from the energy feature TKEO were constantly highly representative of BGLs in all subjects. The features extracted from training data were used to train different machine learning regression models. The best-performing models were ensemble regression trees; these models were optimized for each subject with hyperparameter ranges and validation results, as shown in Table 4.
Figure 10 shows the CEG plot for each subject. Observe that, as desired, all the predicted BGL values reside within regions A and B of the CEG plots thus indicating no adverse effect on clinical decisions based on the proposed device. Moreover, on average 80% of the data for all subjects lie within region A of the CEG. The best prediction results, with around 90% of the predicted points in region A, were obtained for S1, followed by 83.3% for S2, while the lowest predicted points of 77% in region A were for S3 and S4, as shown in Table 5. This is consistent with the fact that the lowest significant difference in high versus low BGL PPG cycles were observed for the untreated type II diabetes subject S4. Least-squares regression analysis indicated a significant positive correlation between the predicted and reference glucose concentrations in all subjects (Table 6).

5. Discussion

The potential of the proposed method for continuous BGL monitoring is evidenced by 100% of the estimated BGL values falling within the clinically acceptable regions A and B of the CEG plots for all four subjects, as shown in Figure 10. The estimated CEG values in Region A fall within 20% of the reference (correct) BGL value, which has no negative effect on the clinical decisions, while the values in Region B of CEG are more than 20% apart from the reference BGL value, but this has low to no impact on clinical outcome [33] and does not lead to inappropriate treatment. This level of accuracy is particularly important for diagnosing hypoglycemia and hyperglycemia, where immediate intervention is required, as clinically wrong diagnosis may lead to the potentially fatal decision to inject insulin in patients with extremely low blood glucose levels. According to the tolerance range of the latest ISO standard (ISO: 15197:2015), the estimated blood glucose levels should be within ±15 mg/dL when the true value is <100 mg/dL, whereas values ≥100 mg/dL should be within ±20 mg/dL [34]. Observe from Figure 10 that the results for S1 and S2 meet these criteria with an average accuracy of 15 mg/dL (±9), while the estimated values for the type I diabetic S3 (accuracy: 19 ± 17 mg/dL) and type II diabetic S4 (35 ± 22 mg/dL) marginally fall outside the recommended range. This can be explained by the range of BGLs in this study being 2–3 times larger than those considered in the existing studies in the literature, with BGLs as high as 345.7 mg/dL for subject S4. For rigor, we also performed the analysis based on the median of PPG cycles, with the results very similar to those reported here which were obtained based on the means of PPG cycles.
Overall, the results obtained using the in-ear sensor are comparable to those obtained in previous studies for other PPG sites, as shown in Table 7. In [17] four different wavelengths were tested, with the best results obtained when all four wavelength measurements were used for estimating BGLs. The results in this work are consistent with the results in [17] when considering data within their tested glucose range (70–152 mg/dL) and the case where they use only two LEDs (880 nm and 940 nm) rather than four. While the choice of an optimal NIR wavelength for BGL monitoring is a subject of our continuing research, it is important to note that other glucose absorption peaks do exist at higher wavelengths [11,30,31], but some constituents of blood (e.g., water) also have absorption peaks at this range and would therefore mask the absorption by the glucose molecules.
In this work, the ear canal was used as a measurement site, owing to its unique and distinguishing advantages over body-worn sensing locations. These include:
  • Stable body/skin temperature, which in turn means that in-ear PPG measurements do not suffer from the adverse effects of vasoconstriction, the phenomenon of the constriction of blood vessels due to local body temperature changes, which in turn affects the amplitude of the PPG signal and thus may bias the estimated BG values [35];
  • Constant pressure/tension between the skin and the sensor yields stable and robust measurements (minimal day-to-day variation in the recordings);
  • The relative position between the head and the body remains mainly constant in most daily activities, which allows for truly continuous and stable long-term PPG measurements;
  • The in-ear sensor is generic and its viscoelastic nature means that it fits all ears;
  • Affordability and scalability, with only one off-the-shelf PPG sensor and only one wavelength needed for clinically acceptable accuracy.

6. Conclusions

This proof-of-concept study introduced an in-ear device for continuous blood glucose monitoring using an off-the-shelf pulse oximeter (PPG) sensor. An LED with a wavelength of 880 nm was used to acquire the PPG signal from which the blood glucose levels were estimated through machine learning. The subjects recruited for this study cover the whole spectrum of non-diabetic, pre-diabetic, and diabetic (type I and type II) conditions. The data were recorded over nine days while fasting and after the first daily meal (breakfast), over at least two hours. Statistically significant features were extracted from the PPG signals in order to capture the relationship between the PPG waveform and the reference glucose measurement from a commercially available glucometer. These features, when used for training ensemble regressors, resulted in the mean accuracy of 82% of the predicted BGLs residing within region A of the Clarke error grid (CEG). The highest clinical accuracy was achievable at around 90% of predicted values residing within region A of the CEG, while for all the subjects 100% of the predicted values reside within the CEG regions A and B—indicating no clinical misdiagnosis. Overall, this study demonstrated the potential of the ear canal as a site for non-invasive blood glucose monitoring. Our future work will focus on a large-scale study over a bigger cohort of subjects and longer monitoring hours (including between meals), which will facilitate the use of more advanced, data-hungry, machine learning techniques for BGL prediction from PPG.

Author Contributions

Conceptualization, G.H.; methodology, G.H.; software, G.H.; validation, G.H. and D.P.M.; formal analysis, G.H.; investigation, G.H.; resources, G.H. and D.P.M.; data curation, G.H.; writing—original draft preparation, G.H.; writing—review and editing, G.H. and D.P.M.; visualization, G.H.; supervision, D.P.M.; project administration, D.P.M.; funding acquisition, D.P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the UK Dementia Research Institute at Imperial College London.

Institutional Review Board Statement

The recordings were performed under Imperial College ethics committee approval JRCO 20IC6414.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data are not publicly available due to ethical restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lin, X.; Xu, Y.; Pan, X.; Xu, J.; Ding, Y.; Sun, X.; Song, X.; Ren, Y.; Shan, P.F. Global, regional, and national burden and trend of diabetes in 195 countries and territories: An analysis from 1990 to 2025. Sci. Rep. 2020, 10, 14790. [Google Scholar] [CrossRef]
  2. Roglic, G. WHO Global report on diabetes: A summary. Int. J. Noncommun. Dis. 2016, 1, 3–8. [Google Scholar] [CrossRef]
  3. Ravaut, M.; Sadeghi, H.; Leung, K.K.; Volkovs, M.; Kornas, K.; Harish, V.; Watson, T.; Lewis, G.F.; Weisman, A.; Poutanen, T.; et al. Predicting adverse outcomes due to diabetes complications with machine learning using administrative health data. NPJ Digit. Med. 2021, 4, 24. [Google Scholar] [CrossRef]
  4. Sun, H.; Saeedi, P.; Karuranga, S.; Pinkepank, M.; Ogurtsova, K.; Duncan, B.B.; Stein, C.; Basit, A.; Chan, J.C.N.; Mbanya, J.C.; et al. IDF Diabetes Atlas: Global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res. Clin. Pract. 2021, 183, 109119. [Google Scholar] [CrossRef]
  5. Renard, E. Monitoring glycemic control: The importance of self-monitoring of blood glucose. Am. J. Med. 2005, 118, 12–19. [Google Scholar] [CrossRef]
  6. Vashist, S.K.; Zheng, D.; Al-Rubeaan, K.; Luong, J.H.; Sheu, F.S. Technology behind commercial devices for blood glucose monitoring in diabetes management: A review. Anal. Chim. Acta 2011, 703, 124–136. [Google Scholar] [CrossRef]
  7. Heinemann, L. Finger pricking and pain: A never ending story. J. Diabetes Sci. Technol. 2008, 2, 919–921. [Google Scholar] [CrossRef] [Green Version]
  8. Garg, S.K.; Voelmle, M.; Gottlieb, P.A. Time lag characterization of two continuous glucose monitoring systems. Diabetes Res. Clin. Pract. 2009, 87, 348–353. [Google Scholar] [CrossRef]
  9. Shang, T.; Zhang, J.Y.; Thomas, A.; Arnold, M.A.; Vetter, B.N.; Heinemann, L.; Klonoff, D.C. Products for monitoring glucose levels in the human body with noninvasive optical, noninvasive fluid sampling, or minimally invasive technologies. J. Diabetes Sci. Technol. 2022, 16, 168–214. [Google Scholar] [CrossRef]
  10. Dziergowska, K.; Labowska, M.B.; Gasior-Glogowska, M.; Kmiecik, B.; Detyna, J. Modern noninvasive methods for monitoring glucose levels in patients: A review. Bio-Algorithms Med-Syst. 2019, 15, 20190052–20190074. [Google Scholar] [CrossRef]
  11. Yadav, J.; Rani, A.; Singh, V.; Murari, B.M. Prospects and limitations of non-invasive blood glucose monitoring using near-infrared spectroscopy. Biomed. Signal Process. Control 2015, 18, 214–227. [Google Scholar] [CrossRef]
  12. Islam, T.T.; Ahmed, M.S.; Hassanuzzaman, M.; Bin Amir, S.A.; Rahman, T. Blood Glucose Level Regression for Smartphone PPG Signals Using Machine Learning. Appl. Sci. 2021, 11, 618. [Google Scholar] [CrossRef]
  13. Lee, E.; Lee, C.Y. PPG-Based Smart Wearable Device With Energy-Efficient Computing for Mobile Health-Care Applications. IEEE Sens. J. 2021, 21, 13564–13573. [Google Scholar] [CrossRef]
  14. Susana, E.; Ramli, K.; Murfi, H.; Apriantoro, N.H. Non-Invasive Classification of Blood Glucose Level for Early Detection Diabetes Based on Photoplethysmography Signal. Information 2022, 13, 59. [Google Scholar] [CrossRef]
  15. Hina, A.; Saadeh, W. A Noninvasive Glucose Monitoring SoC Based on Single Wavelength Photoplethysmography. IEEE Trans. Biomed. Circuits Syst. 2020, 14, 504–515. [Google Scholar] [CrossRef]
  16. Hina, A.; Nadeem, H.; Saadeh, W. A Single LED Photoplethysmography-Based Noninvasive Glucose Monitoring Prototype System. In Proceedings of the 2019 IEEE International Symposium on Circuits and Systems (ISCAS), Sapporo, Japan, 26–29 May 2019; pp. 1–5. [Google Scholar]
  17. Rachim, V.P.; Chung, W.Y. Wearable-band type visible-near infrared optical biosensor for non-invasive blood glucose monitoring. Sens. Actuators B Chem. 2019, 286, 173–180. [Google Scholar] [CrossRef]
  18. Yadav, J.; Rani, A.; Singh, V.; Mohan Murari, B. Investigations on Multisensor-Based Noninvasive Blood Glucose Measurement System. J. Med. Devices 2017, 11, 1932–6181. [Google Scholar] [CrossRef]
  19. Nilsson, L.; Goscinski, T.; Kalman, S.; Lindberg, L.G.; Johansson, A. Combined photoplethysmographic monitoring of respiration rate and pulse: A comparison between different measurement sites in spontaneously breathing subjects. Acta Anaesthesiol. Scand. 2007, 51, 1250–1257. [Google Scholar] [CrossRef]
  20. Goverdovsky, V.; Looney, D.; Kidmose, P.; Mandic, D.P. In-Ear EEG From Viscoelastic Generic Earpieces: Robust and Unobtrusive 24/7 Monitoring. IEEE Sens. J. 2016, 16, 271–277. [Google Scholar] [CrossRef]
  21. Goverdovsky, V.; von Rosenberg, W.; Nakamura, T.; Looney, D.; Sharp, D.J.; Papavassiliou, C.; Morrell, M.J.; Mandic, D.P. Hearables: Multimodal physiological in-ear sensing. Sci. Rep. 2017, 7, 6948. [Google Scholar] [CrossRef] [Green Version]
  22. Nakamura, T.; Alqurashi, Y.D.; Morrell, M.J.; Mandic, D.P. Hearables: Automatic Overnight Sleep Monitoring with Standardized In-Ear EEG Sensor. IEEE Trans. Biomed. Eng. 2020, 67, 203–212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Looney, D.; Kidmose, P.; Park, C.; Ungstrup, M.; Rank, M.; Rosenkranz, K.; Mandic, D. The In-the-Ear Recording Concept: User-Centered and Wearable Brain Monitoring. IEEE Pulse 2012, 3, 32–42. [Google Scholar] [CrossRef] [PubMed]
  24. Hammour, G.M.; Mandic, D.P. Hearables: Making Sense from Motion Artefacts in Ear-EEG for Real-Life Human Activity Classification. In Proceedings of the 43rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Virtual, 1–5 November 2021; pp. 6889–6893. [Google Scholar] [CrossRef]
  25. Hammour, G.; Yarici, M.; Rosenberg, W.V.; Mandic, D.P. Hearables: Feasibility and Validation of In-Ear Electrocardiogram. In Proceedings of the 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Berlin, Germany, 23–27 July 2019; pp. 5777–5780. [Google Scholar] [CrossRef]
  26. Davies, H.J.; Williams, I.; Peters, N.S.; Mandic, D.P. In-Ear SpO2: A Tool for Wearable, Unobtrusive Monitoring of Core Blood Oxygen Saturation. Sensors 2020, 20, 4879. [Google Scholar] [CrossRef]
  27. Davies, H.J.; Williams, I.; Hammour, G.; Yarici, M.; Seemungal, B.M.; Mandic, D.P. In-Ear SpO2 for Classification of Cognitive Workload. IEEE Trans. Cogn. Dev. Syst. 2022. accepted. [Google Scholar] [CrossRef]
  28. Davies, H.J.; Bachtiger, P.; Williams, I.; Molyneaux, P.L.; Peters, N.S.; Mandic, D. Wearable In-Ear PPG: Detailed Respiratory Variations Enable Classification of COPD. IEEE Trans. Biomed. Eng. 2022. in print. [Google Scholar] [CrossRef] [PubMed]
  29. Ramasahayam, S.; Arora, L.; Roy Chowdhury, S. FPGA Based Smart System for Non Invasive Blood Glucose Sensing Using Photoplethysmography and Online Correction of Motion Artifact. Smart Sens. Meas. Instrum. 2017, 22, 1–21. [Google Scholar] [CrossRef]
  30. Yadav, J.; Rani, A.; Singh, V.; Murari, B.M. Design of Low Cost Blood Glucose Sensing System Using Diffused Reflectance Near-Infrared Light. In Networking Communication and Data Knowledge Engineering; Perez, G.M., Mishra, K.K., Tiwari, S., Trivedi, M.C., Eds.; Springer: Singapore, 2018; pp. 197–216. [Google Scholar]
  31. Al-Dhaheri, M.A.; Mekkakia-Maaza, N.E.; Mouhadjer, H.; Lakhdari, A. Noninvasive blood glucose monitoring system based on near-infrared method. Int. J. Electr. Comput. Eng. 2020, 10, 1736–1746. [Google Scholar] [CrossRef]
  32. MATLAB Machine Learning and Deep Learning Toolbox, MATLAB R2022a; The MathWorks: Natick, MA, USA, 2022.
  33. Sengupta, S.; Handoo, A.; Haq, I.; Dahiya, K.; Mehta, S.; Kaushik, M. Clarke Error Grid Analysis for Performance Evaluation of Glucometers in a Tertiary Care Referral Hospital. Indian J. Clin. Biochem. 2022, 37, 199–205. [Google Scholar] [CrossRef]
  34. BS ENISO 15197:2015; In Vitro Diagnostic Test Systems. Requirements for Blood-Glucose Monitoring Systems for Self-Testing in Managing Diabetes Mellitus. British Standards Institute: London, UK, 2015.
  35. Fine, J.; Branan, K.L.; Rodriguez, A.J.; Boonya-Ananta, T.; Ajmal; Ramella-Roman, J.C.; McShane, M.J.; Coté, G.L. Sources of Inaccuracy in Photoplethysmography for Continuous Cardiovascular Monitoring. Biosensors 2021, 11, 126. [Google Scholar] [CrossRef]
Figure 1. The principle of light reflectance and absorption through the different layers of the skin.
Figure 1. The principle of light reflectance and absorption through the different layers of the skin.
Sensors 23 03319 g001
Figure 2. The in-ear PPG sensor used in our study. (a) The different components of the in-ear sensor. (b) The PPG sensor chip embedded on a viscoelastic earbud attached to an ear hook. (c) Sensor placement within the ear canal.
Figure 2. The in-ear PPG sensor used in our study. (a) The different components of the in-ear sensor. (b) The PPG sensor chip embedded on a viscoelastic earbud attached to an ear hook. (c) Sensor placement within the ear canal.
Sensors 23 03319 g002
Figure 3. The recording protocol in this study.
Figure 3. The recording protocol in this study.
Sensors 23 03319 g003
Figure 4. Pre–processing steps for the raw PPG. Top to bottom: raw signal; band–pass signal with the lower envelope detected; upsampled signal with the lower envelope removed; pre–processed signal with the peaks and troughs identified.
Figure 4. Pre–processing steps for the raw PPG. Top to bottom: raw signal; band–pass signal with the lower envelope detected; upsampled signal with the lower envelope removed; pre–processed signal with the peaks and troughs identified.
Sensors 23 03319 g004
Figure 5. Matched filtering of PPG cycles. (Top left) Exemplar template used for matched filtering of PPG cycles. (Top right) Overlay plot of all PPG cycles acquired during one recording. (Bottom left) Exemplar low–quality PPG cycles, with a correlation < 0.9 with the template. (Bottom right) Overlay plot of the remaining PPG cycles after matched filtering.
Figure 5. Matched filtering of PPG cycles. (Top left) Exemplar template used for matched filtering of PPG cycles. (Top right) Overlay plot of all PPG cycles acquired during one recording. (Bottom left) Exemplar low–quality PPG cycles, with a correlation < 0.9 with the template. (Bottom right) Overlay plot of the remaining PPG cycles after matched filtering.
Sensors 23 03319 g005
Figure 6. Averaged low–glucose (blue) and high–glucose (red) PPG cycles with the corresponding shaded areas indicating the confidence intervals. The green dots designate time instances where the null hypothesis of no statistical difference between the low–glucose and high–glucose PPG cycles is rejected.
Figure 6. Averaged low–glucose (blue) and high–glucose (red) PPG cycles with the corresponding shaded areas indicating the confidence intervals. The green dots designate time instances where the null hypothesis of no statistical difference between the low–glucose and high–glucose PPG cycles is rejected.
Sensors 23 03319 g006
Figure 7. Features extracted from every PPG cycle.
Figure 7. Features extracted from every PPG cycle.
Sensors 23 03319 g007
Figure 8. The distribution of the true glucose levels of the calibration and prediction data for each subject (S1–S4), recorded using the gold standard glucometer. Histograms (horizontal bars) show the relative number of measurements at different glucose levels (y-axis) and the box plots indicate the summary statistics of all recordings for each subject.
Figure 8. The distribution of the true glucose levels of the calibration and prediction data for each subject (S1–S4), recorded using the gold standard glucometer. Histograms (horizontal bars) show the relative number of measurements at different glucose levels (y-axis) and the box plots indicate the summary statistics of all recordings for each subject.
Sensors 23 03319 g008
Figure 9. Percentage difference between the high and low glucose levels, based on the values of the ten most significant features across the four subjects. The % difference indicates the importance of a feature.
Figure 9. Percentage difference between the high and low glucose levels, based on the values of the ten most significant features across the four subjects. The % difference indicates the importance of a feature.
Sensors 23 03319 g009
Figure 10. Clarke error grid analysis (EGA) of the glucose predictions for a non-diabetic subject S1, pre-diabetic subject S2, type I diabetic subject S3, and type II diabetic subject S4 (black dots: Day 8 samples, red dots: Day 9 samples, orange line: linear least-squares regression of the reference versus the predicted glucose concentrations).
Figure 10. Clarke error grid analysis (EGA) of the glucose predictions for a non-diabetic subject S1, pre-diabetic subject S2, type I diabetic subject S3, and type II diabetic subject S4 (black dots: Day 8 samples, red dots: Day 9 samples, orange line: linear least-squares regression of the reference versus the predicted glucose concentrations).
Sensors 23 03319 g010
Table 1. Dataset description.
Table 1. Dataset description.
Subject1234
GenderFemaleMaleMaleMale
Age (years)30473859
Glucose range (mg/dL)75.6–142.286.4–194.463.0–174.6183.6–345.7
No. of days between first and last recording28453211
No. of finger pricks56604484
No. of PPG recordings56604484
No. of PPG cycles65086435664812476
Health statusNon-diabeticPre diabeticType I diabeticType II diabetic
Table 2. Summary of features extracted from each PPG cycle.
Table 2. Summary of features extracted from each PPG cycle.
FeatureDescription
DC valueAverage value of the low-passed PPG signal
Peak valueAmplitude of the systolic peak [12]
Delta ACDifference between the peak value and the trough on the left
Rise slopeSlope of the line between the trough on the left of the peak to the peak value
Rise timeTime difference between the trough on the left of the peak and the peak value
Fall levelDifference between the values of the peak and the trough on the right
Fall slopeSlope of the line between the trough on the right of the peak and the peak value
Fall timeTime between the trough on the right of the peak and the peak value
Width at 1/4 heightValue at 1/4 of the height of Delta AC
Width at 3/4 heightValue at 3/4 of the height of Delta AC
Base widthDifference between the trough on the right of the peak and the trough on its left
AreaArea enclosed between the two troughs
Optical density O D = log ( 1 + d e l t a A C D C ) [17]
TKEO meanMean, variance, standard
TKEO variancedeviation, skewness and
TKEO stdkurtosis of the TKEO calculated
TKEO skewnessfrom each cycle using:
TKEO kurtosis x [ n ] = x 2 [ n ] x [ n + 1 ] x [ n 1 ]
Table 3. Description of the training and test data.
Table 3. Description of the training and test data.
Subject1234
Total No. of cycles65086435664812,476
No. of valid cycles5473541944339935
No. of training cycles3524385933557607
No. of testing cycles1949156010782328
Training glucose range (mg/dL)75.6–129.686.4–194.463.0–160.2183.6–342
Testing glucose range (mg/dL)88.2–142.295.4–15375.6–174.6199.8–345.7
Table 4. Optimal parameters of the ensemble regression models and the RMSE of the 5-fold cross-validation results.
Table 4. Optimal parameters of the ensemble regression models and the RMSE of the 5-fold cross-validation results.
Subject1234
Ensemble method (Bag/LSBoost)BagBagLSBoostLSBoost
No. of learners (10–500)26304411
Minimum leaf size (1–18)4811
No. of predictors to sample (1–19)111919
RMSE (mg/dL)11.9827.9320.5939.03
Table 5. Predictions of the blood glucose levels.
Table 5. Predictions of the blood glucose levels.
Subject% of Predictions in Region A of CEG% of Predictions in Region A & B of CEGRMSE (mg/dL)Correlation (r) between the Predicted and Reference ValueMean Absolute Relative Difference (MARD) (%)
S189.4810017.540.6412.8
S283.3310017.780.5413.3
S377.7810025.250.6818.44
S477.7710042.780.5613.5
Table 6. The regressed coefficients (slope and y-intercept), standard error (SE), and coefficient of determination ( R 2 ).
Table 6. The regressed coefficients (slope and y-intercept), standard error (SE), and coefficient of determination ( R 2 ).
SubjectSlope (95% C.I.)y-Intercept (95% C.I.)Standard Error (mg/dL) R 2
S10.078 (0.029, 0.127)94.05 (88.48, 99.62)6.610.40
S20.351 (0.060, 0.641)84.7 (47.99, 121.4)47.540.29
S30.419 (0.010, 0.829)70.59 (22.45, 118.7)45.950.46
S40.255 (0.057, 0.454)189.8 (135.1, 244.4)78.280.32
Table 7. Related work.
Table 7. Related work.
ReferenceWavelength (nm)CEG A (%)Glucose Range (mg/dL)Sensor Placement
Rachim et al. [17]535, 660, 850, 95010070–152Wrist
Ramasahayam et al. [29]935, 950, 107095.3870–180Finger
Lee et al. [13]94088.5775–230Wrist
Hina et al. [16]94097.580–180Finger
This work88082.0563.0–345.7Ear canal
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hammour, G.; Mandic, D.P. An In-Ear PPG-Based Blood Glucose Monitor: A Proof-of-Concept Study. Sensors 2023, 23, 3319. https://doi.org/10.3390/s23063319

AMA Style

Hammour G, Mandic DP. An In-Ear PPG-Based Blood Glucose Monitor: A Proof-of-Concept Study. Sensors. 2023; 23(6):3319. https://doi.org/10.3390/s23063319

Chicago/Turabian Style

Hammour, Ghena, and Danilo P. Mandic. 2023. "An In-Ear PPG-Based Blood Glucose Monitor: A Proof-of-Concept Study" Sensors 23, no. 6: 3319. https://doi.org/10.3390/s23063319

APA Style

Hammour, G., & Mandic, D. P. (2023). An In-Ear PPG-Based Blood Glucose Monitor: A Proof-of-Concept Study. Sensors, 23(6), 3319. https://doi.org/10.3390/s23063319

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop