# Contact, high-resolution spatial diffuse reflectance imaging system for skin condition diagnosis: a first-in-human clinical trial

Anne Koenig, Nils Petitdidier, Henri Grateau, S. Characoun, Abdallah Ghaith, Samuel Verges, Stéphane Doutreleau, Sadok Gharbi, Rémi Gerbelot-Barillon, Sylvain Gioux, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Anne Koenig, Nils Petitdidier, Henri Grateau, S. Characoun, Abdallah Ghaith, et al. Contact, high-resolution spatial diffuse reflectance imaging system for skin condition diagnosis: a first-in-human clinical trial. Journal of Biomedical Optics, 2021, 26 (1), 10.1117/1.JBO.26.1.012706. http://dx.doi.org/10.1117/1.JBO.26.1.012706

HAL Id: hal-03797882

https://hal.science/hal-03797882v1

Submitted on 6 Jan 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Contact, high-resolution spatial diffuse reflectance imaging system for skin condition diagnosis: a first-in-human clinical trial

Anne Koenig<sup>a,*</sup> Nils Petitdidier<sup>o, a, b, *</sup> Henri Grateau,<sup>a</sup>

Samarmar Characoun,<sup>c</sup> Abdallah Ghaith,<sup>c</sup> Samuel Verges,<sup>c</sup>

Stéphane Doutreleau<sup>c</sup>,<sup>c</sup> Sadok Gharbi,<sup>a</sup> Rémi Gerbelot,<sup>a</sup>

Sylvain Gioux,<sup>b,d</sup> and Jean-Marc Dinten<sup>a</sup>

<sup>a</sup>Université Grenoble Alpes, France—CEA, LETI, MINATEC Campus, Grenoble, France

<sup>b</sup>Université de Strasbourg, Strasbourg, France

<sup>c</sup>Université Grenoble Alpes, CHU Grenoble Alpes, INSERM, Hôpital Sud, Laboratoire HP2 (U 1042), Avenue Kimberley, Echirolles, France

<sup>d</sup>Télécom Physique Strasbourg, Laboratoire ICube, Illkirch, France

## Abstract

**Significance:** Oxygenation is one of the skin tissue physiological properties to follow for patient care management. Furthermore, long-term monitoring of such parameters is needed at the patient bed as well as outside the hospital. Diffuse reflectance spectroscopy has been widely used for this purpose.

**Aim:** The aim of the study is to propose a low-cost system for the long-term measurement of skin physiological parameters in contact.

**Approach:** We have developed a low-cost, wearable, CMOS-based device. We propose an original method for processing diffuse reflectance data to calculate the tissue oxygen saturation (StO<sub>2</sub>).

**Results:** We tested the device for the assessment of tissue oxygenation during a first-in-human clinical trial that took place at the Grenoble University Hospital France.

**Conclusions:** The results of this clinical trial show a good accordance between our sensor and commercial devices used a reference.

© The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI. [DOI: 10.1117/1.JBO.26.1.012706]

**Keywords:** skin characterization; optical properties; diffuse reflectance; multipixel sensor; contact imaging; tissue oxygen saturation; wearable device.

Paper 200286R received Sep. 1, 2020; accepted for publication Dec. 15, 2020; published online Jan. 29, 2021.

## 1 Introduction

Knowledge of tissue properties such as oxygen saturation are important in patient care management in a great number of pathologies. Long-term monitoring of such parameters is currently needed at the patient bed as well as outside the hospital. The development of tools capable of monitoring physiological parameters could be of great interest also for applications *in situ* at the hospital such as tissue oxygenation monitoring toward postoperative assessment of superficial flaps viability.¹ In addition, the use of these instruments outside the hospital will become more and more necessary toward the development of remote medicine of tomorrow. Unfortunately, as of today, few instruments are able to assess the tissue viability continuously, in the long term and at low cost. Near-infrared spectroscopy is one of the available techniques providing relative

<footnote>*Address all correspondence to Anne Koenig, anne.koenig@cea.fr; Nils Petitdidier, nils.petitdidier@etu.unistra.fr</footnote>
concentrations of deoxygenated and oxygenated hemoglobin and calculation of an index of tissue oxygenation.² This technique is used in a research setting and in some specialized medical units.³

Diffuse reflectance spectroscopy (DRS) has been widely used in the field of biological tissue characterization with various modalities. One of these modalities consists of measuring the spatially resolved diffuse reflectance (srDR), collecting light at multiple distances from the excitation point. The obtained reflectance decay curve is used to determine scattering and absorption properties of the tissue,⁴ which are directly related to tissue content and structure. Existing spatially resolved diffuse reflectance spectroscopy (srDRS) systems usually use fiber-optic probes to collect light reflected from the tissue and transfer it to an optical sensor.⁵ Recently, efforts have been made to miniaturize the technology through the use of custom photodiodes to detect photons at the tissue surface.⁶ As part of the previous work, we have developed a wearable device based on an innovative architecture for srDRS using a commercially available complementary metal oxide semiconductor (CMOS) sensor placed in contact with the tissue.⁷ In this paper, we present results obtained using this CMOS-DRS device to monitor the oxygenation of tissue *in vivo* during a first-in-human clinical trial where desaturation was controlled and followed by an ischemia. We first describe the instrument, the data processing methodology, and the protocol of the clinical trial. Then, we detail and compare results obtained *in vivo* using this sensor and commercial reference systems.

# 2 Material and Method

## 2.1 Sensor Design

srDRS is a powerful tool for quantitative tissue analysis using cost-effective and compact systems. As a way to overcome the shortcomings of current systems, usually using fibered probes, we proposed an innovative srDRS architecture to achieve contact, high-resolution imaging of the diffuse reflectance from tissue.⁷ While previous studies necessitated the development of custom photodiode arrays,⁶ we proposed to implement a system using low-cost, commercially available components including a CMOS image sensor and light-emitting diodes. The limitations of this approach concerning the access to short source detector separation (SDS) were identified, and an instrumental solution to circumvent these limitations was proposed. This srDRS sensor architecture herein referred to as CMOS-DRS is based on the insertion of a fiber-optic plate (FOP) between the sensor and tissue.⁸

The CMOS-DRS system is shown in Fig. 1. The custom FOP with 6-µm unit fiber core diameter (OS-ST, SZPhoton) is inserted between the sensor [CMOS sensor (UI-1492 IDS imaging Inc.) with 6.413 mm × 4.589 mm area and 1.67 µm pixel pitch] and the analyzed medium. The FOP is positioned to ensure alignment of the pixel area borders with the FOP edges. Four fiber-coupled luminescent electronic diodes (LED) (Kingbright Inc.) emitting at λ = 515, 611, and 660 nm are placed next to the FOP side and are used as excitation sources. They were chosen because of their small lateral dimension (0.8 mm), which was required to provide access to sufficiently small SDS. The sources and FOP are hosted in a specifically designed packaging

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Layout and photograph of the CMOS-based srDRS prototype.</figcaption></figure>
module. The packaging was 3-D printed in a biocompatible polyamide (PA2200). A central hole corresponding to the imaging area was left free for insertion of the FOP. The module is mounted on the camera using plastic screws at the corners of the printed circuit board. For *in vivo* measurements, a dedicated fixation ring is first attached to the tissue using standard medical adhesive tapes. The device is then inserted through the ring and set with a screw to maintain a constant pressure on the tissue. The outer surface of the packaging was painted in black to limit unwanted reflections at the tissue-exterior interface.

In summary, the CMOS-DRS system lighting module has the characteristics listed in Table 1.

The sensor was connected via USB connection to a computer (Dell Latitude E6530, 64 bits, 4 Go RAM, Intel® CoreTM i3). The acquisition was controlled using a dedicated LabVIEW interface (National Instruments Inc., Austin, Texas).

## 2.2 Image Acquisition Sequence

Diffuse reflectance images were stored and processed on the computer. All imaging processing steps were performed using Matlab (Mathworks Inc., Natick, Massachusetts). Raw reflectance images and dark noise images (background subtraction) were first recorded at multiple exposure times. Typically, for the clinical trial experiments, we used two exposure times. The obtained images were recombined to expand the sensor dynamic (dynamic expansion). Then, the geometry of our collecting light system using the CMOS camera was exploited by ring averaging. The pixel values were averaged over concentric annular areas surrounding the illumination center to extract the radial profile of raw diffuse reflectance.⁷

The LEDs were lit in the following order: orange, red, green, with the first exposure time for the first image acquisition, then with the second exposure time (equal to the first exposure time × 8) for the second image acquisition. The exposure times ranges are presented in Table 2.

## 2.3 Method for Data Processing

In the srDR technique, the measured diffuse reflectance decay curves are exploited to determine the scattering and absorption properties of the analyzed tissue.⁷ Data processing was conducted according to the flowchart shown in Fig. 2.

A Monte Carlo (MC)⁹ algorithm was used to simulate the propagation of photons in tissue. The source diameter and numerical aperture (NA) were accounted for in the computation. The theoretical reflectance profiles calculated for various $\mu_a$ and $\mu_s'$ were stored in a look-up table. Following an approach described in previous studies,¹⁰ measured and theoretical reflectance profiles were matched through merit function minimization to extract the optical properties

<table><caption>**Table 1** Lighting module characteristics.</caption><thead><tr><td>Color</td><td>Wavelength (nm)</td><td>SDS min (μm)</td><td>SDS max (mm)</td></tr></thead><tbody><tr><td>Green</td><td>515</td><td>760</td><td>7</td></tr><tr><td>Orange</td><td>611</td><td>820</td><td>7</td></tr><tr><td>Red</td><td>660</td><td>690</td><td>7</td></tr></tbody></table>

<table><caption>**Table 2** Exposure times.</caption><thead><tr><td>Color</td><td>Time 1 (ms)</td><td>Time 2 (ms)</td></tr></thead><tbody><tr><td>Green</td><td>1.3 &lt; t1 &lt; 7.28</td><td>10.52 &lt; t2 &lt; 26.08</td></tr><tr><td>Orange</td><td>1.15 &lt; t1 &lt; 3.26</td><td>9.5 &lt; t2 &lt; 76.1</td></tr><tr><td>Red</td><td>7.5 &lt; t1 &lt; 18.04</td><td>60.3 &lt; t2 &lt; 144.3</td></tr></tbody></table>
of the tissue. Performances of the designed prototype and of the method to extract optical properties (i.e., $\mu_a$ and $\mu_s'$ coefficients) were evaluated with tissue-simulating phantoms.⁷

## 2.4 Calibration Phantoms

To calibrate the system, diffuse reflectance measurements were performed on a phantom having known optical properties. We used a range of homemade solid tissue-simulating phantoms manufactured in polydimethylsiloxane (PDMS). Ethanol-dissolved nigrosin (Ng-EtOH) and titanium dioxide ($TiO_2$) particles were added in the phantom mixture to generate absorption and scattering.⁸,¹¹

Phantoms optical properties at 665 and 860 nm were characterized through spatial frequency domain imaging (SFDI). The obtained $\mu_a$ and $\mu_s'$ optical coefficients were then extrapolated at 511, 615, and 660 nm. Several authors have measured scattering spectra of $TiO_2$-based PDMS phantoms, showing a nearly linear behavior in the 500 to 900 nm range.¹²⁻¹⁴ Therefore, linear extrapolation using the values of $\mu_s'$ available at 665 and 860 nm was used to calculate the reduced scattering coefficient of the calibration phantom at our wavelengths of interest, i.e., 511, 615, and 660 nm. Extrapolation of $\mu_a$ coefficients was based on spectra measured with an srDRS fiber-optic probe previously developed at CEA laboratory.⁵

The calibration phantom measurement was used to yield a scaling factor through comparison with the theoretical data (box “calibration” of the flowchart Fig. 2). The scaling factor was applied to unknown object (phantom or tissue) measurements to convert raw diffuse reflectance profiles into absolute unit.

During this study, we used two calibration phantoms with properties presented in Table 3. Depending on subjects’ prototype, we used one or the other of these phantoms. Several authors pointed out that errors in optical properties determination should be minimal when the scattering coefficients of the calibration phantom and the analyzed medium are the same.¹⁵ In the case of using multiple calibration phantoms, the work done by Sorgato et al.¹⁶ could help with choosing the best adapted.

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Flowchart of the data processing steps for optical properties quantification.</figcaption></figure>

<table><caption>**Table 3** Calibration phantoms optical properties at 665 nm.</caption><thead><tr><th>Name</th><th>Nigrosin concentration (mg/g)</th><th>TiO<sub>2</sub> concentration (mg/g)</th><th>μ<sub>a</sub> @ 665 nm (cm<sup>-1</sup>)</th><th>μ<sub>s</sub>' @ 665 nm (cm<sup>-1</sup>)</th></tr></thead><tbody><tr><td>B1</td><td>0.08</td><td>3</td><td>0.52</td><td>19.4</td></tr><tr><td>B2</td><td>0.239</td><td>3</td><td>0.99</td><td>18.2</td></tr></tbody></table>
## 2.5 Derivation of the Tissue Oxygen Saturation

We have implemented two original methods to extract tissue oxygen saturation from optical coefficients estimates. The first method hereafter referenced as “method A” uses data measured at two wavelengths only and requires a secondary reference system for calibration. To avoid the need of a reference system, a second approach using a three-wavelengths acquisition and processing scheme was investigated, it is referenced as “method B.”

### 2.5.1 Method A

The measured wavelength-dependent absorption coefficient was exploited to determine the relative concentration of oxy- and deoxygenated blood in tissue, which is represented by the tissue oxygen saturation parameter (StO₂). In method A, images recorded at wavelengths λ₁ = 660 nm and λ₂ = 611 nm were used.

Using modeling based on volume fractions proposed by Jacques,¹⁷ absorption in skin was modeled as

$$ \mu_a(\lambda, t) = B(t) * \{\text{StO}_2(t) * \mu_a^{\text{HbO}_2}(\lambda) + [1 - \text{StO}_2(t)] * \mu_a^{\text{Hb}}(\lambda)\} + [1 - B(t)] * \mu_a^{\text{other}}(\lambda, t), \quad (1) $$

where μₐ(λ, t) is the absorption coefficient measured by the CMOS-DRS device over time t (cm⁻¹), B(t) stands for the blood volume fraction in tissue as defined by Jacques¹⁷ or Meglinski and Matcher,¹⁸ StO₂(t) represents the tissue saturation, μₐ<sup>HbO₂</sup>(λ) and μₐ<sup>Hb</sup>(λ) are the tabulated spectra of oxygenated and deoxygenated blood absorption¹⁹ (cm⁻¹), μₐ<sup>other</sup>(λ, t) gathers the contributions of “other” chemical species to the skin absorption (melanin, bilirubin, ...) (cm⁻¹).

Meglinski and Matcher proposed an analytic formula to express a “mean” μₐ<sup>other</sup> as a function of wavelength:¹⁸

$$ \mu_a^{\text{other}}(\lambda) = 7.84 * 10^8 * \lambda^{-3.255}, \quad (2) $$

where λ is the wavelength expressed in nanometers.

To account for intersubject variations in background absorption, especially related to varying melanin content, we introduced an additional calibration factor x(t), the exponential shape of absorption decay with wavelength of Eq. (2) being conserved, μₐ<sup>other</sup>(λ, t) thus writes:

$$ \mu_a^{\text{other}}(\lambda, t) = x(t) * 7.84 * 10^8 * \lambda^{-3.255}, \quad (3) $$

where x(t) is characteristic of each subject.

In order to calculate the StO₂, we solve a system of two equations with two wavelengths:

$$ \left\{ \begin{array}{l} \mu_a(\lambda_1, t) = B(t) * \{\text{StO}_2(t) * \mu_a^{\text{HbO}_2}(\lambda_1) + [1 - \text{StO}_2(t)] * \mu_a^{\text{Hb}}(\lambda_1)\} + [1 - B(t)] * \mu_a^{\text{other}}(\lambda_1, t) \\ \mu_a(\lambda_2, t) = B(t) * \{\text{StO}_2(t) * \mu_a^{\text{HbO}_2}(\lambda_2) + [1 - \text{StO}_2(t)] * \mu_a^{\text{Hb}}(\lambda_2)\} + [1 - B(t)] * \mu_a^{\text{other}}(\lambda_2, t) \end{array} \right. \quad (4) $$

As described in Sec. 2.2, diffuse reflectance images for wavelengths λ₁ and λ₂ were recorded sequentially during *in vivo* experiments. Prior to StO₂ calculation, μₐ measurements for λ₂ were therefore readjusted to the acquisition times of λ₁. The resolution of this system gives the following relations:

With:

$$ \left\{ \begin{array}{l} C_{1i} = \mu_a(t, \lambda_i) - \mu_a^{\text{other}}(\lambda_i, t) \\ C_{2i} = \mu_a^{\text{Hb}}(\lambda_i) - \mu_a^{\text{other}}(\lambda_i, t), \\ C_{3i} = \mu_a^{\text{HbO}_2}(\lambda_i) - \mu_a^{\text{Hb}}(\lambda_i) \end{array} \right. $$
$$
\left\{
\begin{aligned}
\text{StO}_2(t) &= \frac{C_{11}C_{22} - C_{21}C_{12}}{C_{31}C_{12} - C_{32}C_{11}} \\
B(t) &= \frac{C_{11}(C_{31}C_{12} - C_{32}C_{11})}{C_{31}(C_{11}C_{22} - C_{21}C_{12}) + C_{21}(C_{31}C_{12} - C_{32}C_{11})}
\end{aligned}
\right.
\qquad (5)
$$

Calibration of the system was based on the comparison of readings recorded during the initial stable period of the protocol (ambient air step as mentioned in Sec. 2.6) by the CMOS-DRS device and a reference StO₂ sensor (PortaLite, Artinis Medical Systems, The Netherlands, referred hereafter as “Artinis”²⁰). The average StO₂ value provided by the Artinis noted Sbsl, and μa value extracted by the CMOS-DRS system noted μabsl yield the calibration factor x specific to each subject. Note that with this method A, we do not take the variation of x over time into account but consider an average value of x over a stable time interval:

$$
x = \frac{\{S_{bsl} * [\mu_{absl}(\lambda_2) * C_{31} - \mu_{absl}(\lambda_1) * C_{32}] + \mu_{absl}(\lambda_2) * \mu_a^{Hb}(\lambda_1) - \mu_{absl}(\lambda_1) * \mu_a^{Hb}(\lambda_2)\}}{7.84 * 10^8 * \{\lambda_1^{-3.255} * [\mu_{absl}(\lambda_2) - \mu_a^{Hb}(\lambda_2) + S_{bsl} * C_{32}] - \lambda_2^{-3.255} * [\mu_{absl}(\lambda_1) - \mu_a^{Hb}(\lambda_1) + S_{bsl} * C_{31}]\}} \quad (6)
$$

2.5.2 Method B

As a second approach, we used the absorption coefficient measured at a third wavelength (511 nm) to determine the calibration factor x(t), thereby avoiding the use of a reference sensor.
Using Eqs. (1) and (3), and with the following notation, Eq. (1) leads to Eq. (7):

$$
L(\lambda) = 7.84 * 10^8 * \lambda^{-3.255} (\lambda \text{ in nm})
$$

$$
\mu_a(\lambda, t) = B(t) * \text{StO}_2(t)[\mu_a^{\text{HbO}_2}(\lambda) - \mu_a^{\text{Hb}}(\lambda)] + B(t) * \mu_a^{\text{Hb}}(\lambda) + [1 - B(t)] * x(t) * L(\lambda). \quad (7)
$$

With:

$$
\begin{align*}
C_{11} &= \mu_a(\lambda_1, t)L(\lambda_2) \\
       &\quad - \mu_a(\lambda_2, t)L(\lambda_1) \\
C_{12} &= \mu_a(\lambda_2, t)L(\lambda_3) \\
       &\quad - \mu_a(\lambda_3, t)L(\lambda_2) \\
C_{13} &= \mu_a(\lambda_3, t)L(\lambda_1) \\
       &\quad - \mu_a(\lambda_1, t)L(\lambda_3)
\end{align*}
$$

$$
\begin{align*}
C_{21} &= \mu_a^{\text{Hb}}(\lambda_1)L(\lambda_2) \\
&\quad - \mu_a^{\text{Hb}}(\lambda_2)L(\lambda_1) \\
C_{21} &= \mu_a^{\text{Hb}}(\lambda_2)L(\lambda_3) \\
&\quad - \mu_a^{\text{Hb}}(\lambda_3)L(\lambda_2) \\
C_{21} &= \mu_a^{\text{Hb}}(\lambda_3)L(\lambda_1) \\
&\quad - \mu_a^{\text{Hb}}(\lambda_1)L(\lambda_3)
\end{align*}
$$

$$
C_{31} = \mu_a^{\text{HbO}_2}(\lambda_1) - \mu_a^{\text{Hb}}(\lambda_1) \qquad C_{32} = \mu_a^{\text{HbO}_2}(\lambda_2) - \mu_a^{\text{Hb}}(\lambda_2) \qquad C_{33} = \mu_a^{\text{HbO}_2}(\lambda_3) - \mu_a^{\text{Hb}}(\lambda_3)
$$

$$
\begin{align*}
C_{41} &= \mu_a^{\text{Hb}}(\lambda_1)\mu_a(\lambda_2, t) \\
&\quad - \mu_a^{\text{Hb}}(\lambda_2)\mu_a(\lambda_1, t) \\
C_{41} &= \mu_a^{\text{Hb}}(\lambda_2)\mu_a(\lambda_3, t) \\
&\quad - \mu_a^{\text{Hb}}(\lambda_3)\mu_a(\lambda_2, t) \\
C_{41} &= \mu_a^{\text{Hb}}(\lambda_3)\mu_a(\lambda_1, t) \\
&\quad - \mu_a^{\text{Hb}}(\lambda_1)\mu_a(\lambda_3, t)
\end{align*}
$$

The resolution of this new system gives the following relations for StO₂ calculation:

$$
\left\{
\begin{aligned}
B(t) &= \frac{C_{31}C_{12} + C_{32}C_{13} + C_{33}C_{11}}{C_{31}C_{22} + C_{32}C_{23} + C_{33}C_{21}} \\
\text{StO}_2(t) &= \frac{\mu_1 C_{22} + \mu_2 C_{23} + \mu_3 C_{21}}{C_{31}C_{12} + C_{32}C_{13} + C_{33}C_{11}} \\
x(t) &= \frac{C_{31}C_{42} + C_{32}C_{43} + C_{33}C_{41}}{C_{31}(C_{22} - C_{12}) + C_{32}(C_{23} - C_{13}) + C_{33}(C_{21} - C_{11})}
\end{aligned}
\right.
\qquad (8)
$$

2.6 First-in-Human Clinical Trial

We planned a clinical trial that began at the beginning of June 2019. This trial aimed at testing
two sensors that will eventually be on the same platform to measure two physiological param-
eters: partial pressure of carbon dioxide (PtCO₂) and StO₂.²¹ The clinical trial has been
performed in collaboration with INSERM & CHU Grenoble University Hospital. It was planned
to acquire data on a panel of 20 healthy volunteers [male and female, nonsmoker (attested by
exhaled CO < 8 ppb), free from any disease and medication except contraceptives, IMC (body
mass index) < 27] aged from 20 to 50 years old. The subjects' prototype varied from light to
dark skin. Power assessment was based on a minimum expected difference in tissue saturation
index of 2% between the two methods. With an α (type I error rate) level of 5% and power of
80%, 15 subjects were required. To take account of potential technical issues, 20 subjects were
included in the trial.

Our institutional review board (CPP Grenoble Sud Est III) approved this study (CTCO,
38RC18.176) and all subjects included signed a written informed consent. Experiments were
performed according to the Declaration of Helsinki (clinical trial registration: NCT03992651).

The protocol is shown in Fig. 3. For each volunteer, the experiment was divided in several
steps (from A to K) corresponding to different levels of hypoxia, hypocapnia, hypercapnia, or
ischemia. The subject breathed through a nasobucal mask allowing the inhalation of different gas
mixtures controlled by the Altitrainer (ALTITRAINER200®, S.M.TEC, Geneva, Switzerland)²²
system. Hypoxic, hypercapnic, gas mixture are obtained by adding nitrogen, carbon dioxide, to
ambient air.

Before and after each volunteer measurements, acquisitions on reference phantoms were
performed (not mentioned in Fig. 3).

- Three levels of hypoxia were targeted (steps A, B, C). The first target was ambient air at an altitude of 200 m, the subject breathed normally through the nasobucal mask. Then, the gas mixture was changed to simulate breathing at altitudes of 3500 and 4500 m, this last target corresponding to a pulse oxygen saturation level (SpO₂) of ~80%. Once reached, the target hypoxic levels were kept during 3 min at least. To reach distinct levels of SpO₂ for each subject, various levels of inspiratory oxygen fractions (FiO₂) were controlled by the Altitrainer system. It is worth noting that this approach led to some interindividual heterogeneity in SpO₂ levels, however, allowing us to compare the different devices (CMOS-DRS, SpO₂, PortaLite) throughout a relatively large spectrum of SpO₂ levels.

- Ambient air step (D) was done to retrieve normal parameters.

- Three consecutive hypercapnia or hypocapnia were targeted (E, F, G), for the test of the capnometry sensor.²³

- Ambient air step (H) was done to retrieve normal parameters.

- Ischemia of the right arm (step I) was induced using a cuff inflated at 250 mmHg pressure. This step aimed at evaluating the sensitivity of the device to tissue oxygenation changes, following an approach used in previous studies.² During a last step starting from the pressure cuff release, recording of skin parameters was maintained to monitor the reoxygenation of tissue. Two acquisitions were made, one right after cuff release (J) and the other at the end of the experiment (K).

By restricting the supply of blood to the vessels in the arm, upstream of the sensor, ischemia
aims to simulate a quick shortage of oxygen in the forearm. This step of ischemia is the true test
of our system to detect flap failure.

The photographs presented in Figs. 4 and 5 illustrate this experimental protocol. The left arm
was dedicated to capnometry sensors while StO₂ sensors were placed on the right arm.

<figure><img src="image_5.png" /><figcaption>**Fig. 3** Protocol steps summary.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Schematic and photograph of the sensors placement on the volunteer.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 5** Installation steps of sensors for in vivo StO₂ measurement. (a) Attachment of the CMOS-DRS system was performed using a fixation ring and medical adhesive tapes. (b) The wearable CMOS-DRS system was placed on the inner forearm and away from the Artinis illumination field to avoid parasite light signals. (c) An ischemia of the upper right arm was performed using a cuff inflated at 250-mmHg pressure.</figcaption></figure>

Continuous monitoring of tissue optical properties was carried out on the volar right forearm using the CMOS-DRS device and an NIRS system simultaneously. The NIRS PortaLite sensor, which is in routine use in our hospital to measure tissue oxygenation in muscles,²⁴ was chosen as reference system for this study. Particular attention was paid to the installation of the devices to ensure that optical paths of the two systems would not interfere. Diffuse reflectance acquisitions were carried out for all the protocol duration, including hypo- and hypercapnia. The total duration of the test was about 60 min.

The so-named CAPNO device introduces an innovative noninvasive transcutaneous PtCO₂ measurement method (or TcpCO₂). The measurement principle relies on the absorption of infrared radiation by the carbon dioxide released through skin.²⁵ The skin surface is heated to 42°C. CAPNO device drives the carbon dioxide emanating from the skin through an optical measurement cell (thermopile and near-infrared source). The amount of radiation received by the thermopile is related to the concentration of transcutaneous carbon dioxide. During the experiment, measurements were compared with a commercially available electrochemical sensor:(SenTec, V-Sign, Switzerland).²⁶ This system is composed of a transcutaneous electrochemical tcPCO₂
sensor and of a reflectance optical SpO₂ sensor. It allows continuous measurement of tcPCO₂, pulse oximetry (SpO₂, pulse rate). In addition to tcPCO₂ reference measurements, SpO₂ readings were recorded on the left arm using the SenTec sensor. Results concerning CAPNO test are the subject of a separate communication.²³

## 3 Results

We have completed experiments with 20 volunteers, among which 19 were exploitable. We present the results with an emphasis on one subject: subject 14 (Figs. 6 and 7). We present results for the entire study in Figs. 8 and 9. Tables of results for each subject at each step are given in appendix.

### 3.1 Optical Properties

As described in Sec. 2.3, the first step of data processing consists of extracting absorption and scattering properties from reflectance measurements. For most subjects (15), we used as reference PDMS phantom, phantom B1. For darker skin (4), we had to use phantom B2, which is more absorbent to avoid saturation in resulting tissue absorption spectra. With our calibration method, the properties are better estimated if the properties of the reference phantom are close to those of the subject.¹⁵

As an illustration, the optical coefficients $\mu_a$ and $\mu_s'$ measured by the CMOS-DRS device on volunteer 14 at the three acquisition wavelengths are shown in Figs. 6(a) and 6(b), respectively. Letters on top of figure correspond to protocol steps as shown in Fig. 3.

### 3.2 Oxygen Saturation, Method A

The StO₂ levels measured by the Artinis and the CMOS-DRS systems for volunteer 14, using method A are shown in Fig. 7(a). SpO₂ readings recorded by the SenTec system on the left arm are also displayed. We observed that the SpO₂ measurements are much more sensitive to desaturation (A, B, C) than StO₂ readings. Changes in arterial blood versus tissue oxygenation in response to changes in inspiratory oxygen pressure can indeed show distinct patterns, as previously shown by our group.²⁷,²⁸ Tissue vascular (increase in regional blood flow or capillary recruitment) and metabolic (e.g., reduction in tissue metabolism in hypoxia) responses can lead to smaller reductions in local oxygen saturation compared with arterial oxygen saturation. Therefore, in the context of this study inducing hypoxic conditions by decreasing inspiratory

<figure><img src="image_8.png" /><figcaption>**Fig. 6** Optical properties: (a) absorption coefficient, (b) reduced scattering coefficient, measured for subject 14 using phantom B2 as reference, at 515 nm (green), 611 nm (orange), and 660 nm (red). Letters on top of figures correspond to protocol steps, and black dashed lines correspond to ischemia steps I and J.</figcaption></figure>
oxygen pressure (which can be different from tissue hypoxia due to local ischemia for instance), this distinct pattern of $SpO_2$ versus $StO_2$ seems to be expected.

Conversely, the ischemia (steps I, J) has no effect on the $SpO_2$ level as oxygenation changes are induced only locally by the cuff inflation. This confirms that the $StO_2$ decrease detected by the Artinis and the CMOS-DRS sensor is exclusively induced by the cuff inflation. $StO_2$ levels yielded by the CMOS-DRS device match the values obtained from reference measurements within 5%.

The $StO_2$ levels measured by the Artinis and the CMOS-DRS systems were compared for all subjects, as for subject 14. To resume the comparison between the two sensors, we calculate representative figures over the whole subject's panel. The $StO_2$ evolution obtained from the CMOS-DRS and Artinis system averaged for each step, and averaged over all 19 subjects is shown in Fig. 8(a). In Fig. 8(aa), we plotted the $StO_2$ value distribution for each step and for the two systems. On each box (red for reference sensor and blue for CMOS-DRS sensor), the central mark indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles of the measurements, respectively. The whiskers extend to the most extreme data points considered as outliers (Matlab® boxplot function). Results show that mean values are in accordance [Fig. 8(a)] and that the distribution of values is typically more extensive for the CMOS-DRS sensor [Fig. 8(aa)].

During desaturation steps, a slow decrease (up to 10% with the Artinis sensor for the subject of Fig. 7) was observed. With the CMOS-DRS sensor, the first hypoxia step (B) is not visible; a slight decrease is observed for the second hypoxia step (C).

The desaturation steps are slightly better detected in Fig. 8(a) with the Artinis reference sensor than with our CMOS-DRS sensor. However, the recovery to normal saturation (ambient and capnia steps) follows the same slope.

The ability of the CMOS-DRS sensor to monitor the decrease of tissue oxygenation during ischemia is shown in Figs. 7(a) and 8(a). Moreover, the slight rise in $StO_2$ value, typically observed after cuff release,² is also detected both for the CMOS-DRS device and reference sensor. The hypoxia decrease (I) is of the same order as for the Artinis sensor, the increase after recovery (J) is detected.

The deviation of CMOS-DRS mean values with respect to reference sensor step values over all subjects at each protocol step is shown in Fig. 9(a). The deviation to the Artinis reference sensor is under 5% for all protocol steps excluding the ischemia step, where the $StO_2$ decrease

<figure><img src="image_9.png" /><figcaption>**Fig. 7** (a) Oxygen saturation readings for subject 14 using method A, (b) using method B. Measurements of the tissue oxygen saturation obtained using the CMOS-DRS (blue line) and Artinis (red line) systems are plotted along with the pulsed oxygen saturation measured by the SenTec sensor (red dotted line). At the bottom of figure are indicated in blue mean CMOS-DRS $StO_2$ values for last 30 s of each protocol step, and near the curves min and max of $StO_2$ value at ischemia for both sensors (in blue for CMOS-DRS and in red for Artinis). Vertical black dashed lines and letters on top correspond to protocol steps.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 8** (a), (b) StO₂ mean value for all subjects measured with the reference (red) and CMOS-DRS (blue) StO₂ sensors. (aa), (bb) Value distribution for each step as measured with the reference (red) and CMOS-DRS sensor (blue). On each box, the median value is represented by the central mark. Bottom and top edges of the box indicate the 25th and 75th percentiles of the measurements, respectively. (a, aa) using method A. (b), (bb) using method B.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 9** Mean deviation between StO₂ values measured with the CMOS-DRS sensor and the Artinis. For each step, deviations are averaged over all subjects and expressed in percent of the expected value (Artinis measurement). Error bars correspond to the standard deviation of this mean deviation over all subjects. (a) Using method A. (b) Using method B.</figcaption></figure>
is typically overestimated by 10%. The largest uncertainty lies on the first steps, we noted that our sensor requires a time of stabilization after installation on the skin. Furthermore, for the ischemia step, for two subjects, we observed an overestimation of the decrease of the StO₂ reaching almost 0. This explains the huge error bar of step I.

## 3.3 Oxygen Saturation, Method B

We present in Fig. 7(b) the results of the StO₂ calculations for volunteer 14, using method B, which is independent of an anterior calibration step using a reference sensor. For subject 14 shown in Fig. 7(b), the two hypoxia steps (B and C) are now visible in comparison with Fig. 7(a) where they were not. The Artinis values go from 75% (step A) to 69% (step C), whereas CMOS-DRS values go from 70% (step A) to 62% (step C).

For ischemia, the same observations as for the method A are made. The ischemia decrease (J) is of the same order as for the Artinis sensor from 75% (step I) to 49% (step J) for the Artinis and from 70% (step I) to 50% (step J) for the CMOS-DRS sensor. The increase after recovery (K) is well detected in both cases.

As for method A, we plotted in Figs. 8(b) and 8(bb) the mean StO₂ value for all 19 subjects and value distribution after normalization for each step using method B for calculation. We observed a difference between Artinis and CMOS-DRS curves for each step. CMOS-DRS values are for all steps below Artinis ones, by an offset of almost 8%. However, desaturation and ischemia decreases follow the same slope.

We plotted in Fig. 9(b), the deviation to the reference sensor calculated for all subjects at each protocol step. Using method B, the deviation to the Artinis reference sensor is now around 10% for all protocol steps because of the negative offset. The greatest uncertainty lies again on the first steps and ischemia (for two subjects the decrease is overestimated).

# 4 Discussion

Despite the great accordance between our CMOS-DRS sensor and the Artinis, we observe discrepancies for certain subjects. The methods are quite dependent on the optical properties computation, the first step of data processing. This result is itself dependent on the calibration process, using a phantom of known optical properties as reference. For example, we have observed, for certain subjects, similar values for $\mu_a$ parameter but a cross-talk in $\mu_s'$ parameter, values at 611 and 660 nm being inversed. This phenomenon suggests that the reference phantom was not adapted. In those cases, we adjusted the reference phantom, using B1 for lighter skin and B2 for darker ones.

One limitation identified and discussed in the first publication⁷ was the choice of wavelengths. We recommend for a second future version of the CMOS-DRS sensor to use an LED in the near-infrared range (around 940 nm as in most oxymeters) instead of 611 nm. Actually in the range 590 to 640 nm, the absorption of oxyhemoglobin varies very rapidly and the measurement is therefore highly sensitive to the chosen wavelength. The distribution of values, typically more extensive for the CMOS-DRS sensor, may be related to uncertainties in the optical properties determination prior to StO₂ computations.

A first analysis after the acquisition on the first six subjects showed that our sensor requires a time of stabilization after installation on the skin, estimated to 5 min. This is clearly visible on subject 14 shown in Fig. 7(a). For the first hypoxia level (step B), we can see the StO₂ value goes on increasing and starting to decrease only upon second hypoxia level (step C), whereas for the reference Artinis sensor the hypoxia is clearly visible for the two steps B and C. This may be due to a slightly local heating of the skin by the LEDs. Note that the Artinis system was turned on immediately upon installation on the arm and prior to the CMOS-DRS device, which was turned on only at the beginning of step A. To confirm this hypothesis, we started CMOS-DRS measurements prior to the first step A for the two last subjects of the trial, for which the StO₂ decrease was detected at step B. Note that method B, which uses the green acquisition to learn the “other” absorption coefficient, allows compensation for this stabilization time, by correcting the signal through a time dependent $\mu_a^{\text{other}}$ parameter. For subject 14, as shown
in Fig. 7(b), the two hypoxia steps (B and C) are then visible. In this trial, we generate hypoxia and ischemia over a short enough time to test first the sensitivity of the system and second the capacity of a short acquisition time. For continuous monitoring (of the order of a few seconds), we are effectively limited by this stabilization time. In the monitoring of oxygenation of the flaps, thrombosis can occur quite “suddenly” but the order of magnitude is more like 1 h.¹,³ In this case, an acquisition every 5 min might be sufficient and stabilization time is less of a problem.

Another difference with the Artinis reference sensor is the overestimation of ischemia. For 9/19 subjects, the CMOS-DRS sensor overestimates the StO₂ decrease as it can be seen on mean values of Figs. 8(a) and 8(b). However, it is generally in accordance with the Artinis for the increase after recovery. We can remark that method B allows the signals to be readjusted in a given oxygenation state, yet a negative offset of almost 8% is observed. This difference in absolute values may be explained by the fact that the two sensors do not probe the tissue at the same depths. The NIRS has a measurement depth of about half the distance between the transmitter and the receiver (which is in our case ~2 cm), meaning that both the dermis and muscles can be probed. Conversely, the CMOS DRS was designed to probe the dermis only, i.e., limited in depth to a few millimeters.

This work shows that the CMOS-DRS is a system allowing qualitative continuous monitoring of StO₂. In addition, the use of the third wavelength (method B) makes our system autonomous and independent of a reference system. However, the discrepancies observed between quantitative measurements of StO₂ between CMOS-DRS system and PortaLite Artinis reference system may not entirely explained by the different probing depths of devices. Further studies need to be engaged to fully assess the ability of CMOS-DRS approach to provide quantitative measurements. In particular, a complementary study involving oxygenated and deoxygenated hemoglobin phantoms with data on a larger spectral range should be conducted. The experiments presented in this paper, however, demonstrate that our approach is suitable for qualitative detection of ischemia, showing potential for the detection of free flaps failure in postoperative monitoring.

## 5 Conclusion

This clinical trial that aimed at testing two sensors for the monitoring of PtCO₂ and StO₂ allowed us to investigate their implementation on the same platform.²¹ Regarding the CMOS-DRS sensor, this clinical trial and the comparison with reference sensors helped us to understand the μₐ estimation discrepancy previously seen.⁷ Results obtained assess the potential of the system for in vivo tissue oxygenation monitoring in clinical settings.

In the future, further work has to be done toward the development of a fully integrated sensor. This will notably involve the development of an embedded software for the online computation of the StO₂ and eventually the use of Bluetooth technology to suppress the USB connection.

Besides, this work has been done with the underlying objective of layered media analysis. A theoretical validation of the approach has already been done through simulations of multi-layer propagation with both vertical and angled sources.⁸ In particular, an angled source or detector incidence enhances the sensitivity to shallow depths. For this purpose, a fourth source having 35-deg incidence angle has been implemented in the device for experimental verification. Further studies will be carried out to validate the approach, which may lead to improvements in measurements of optical parameters in the dermis independently of the epidermis.

In summary, a CMOS-based contact imaging system for srDR measurements has been developed. This first proof-of-concept has led to the development of a low-cost, wearable prototype designed for the monitoring of tissue oxygenation in vivo. The instrument was evaluated on tissue-simulating phantoms in comparison with SFDI in a preliminary work⁷ and to NIRS systems in this study. The preliminary study assessed the capacity of our approach for the qualitative in vivo monitoring of tissue oxygenation, whereas the clinical trial demonstrates a real prospect for quantitative results. In particular, our approach based on a three wavelengths
acquisition shows great potential for quantitative results independently of any calibration reference device.

Others advantages of the CMOS-DRS technique over the “noncontact imaging” modality are low cost, small footprint, a reduced sensitivity to movement, an easier positioning of the system on the patient that make it easier to use for nursing staff.

Future work will address further validation for quantitative results, the integration aspect of the system, through software embedment and implementation of wireless data transfer, making eventually a low cost, easy-of-use system for continuous StO₂ monitoring. The characterization work will be continued to make the measurement quantitative and allow in the future the use of the CMOS-DRS system for the monitoring of flaps and the detection of flap failure. Moreover, the coupling with the measurement of ptCO₂ on the same platform may be a real advantage in the monitoring of respiratory pathologies such as sleep apnea for example.

# 6 Appendix

The StO₂ results for each subject at each protocol step, are given in Tables 4–6 for respectively, the CMOS-DRS method A, CMOS-DRS method B, and the reference device Artinis PortaLite. In these tables, StO₂ values are calculated by taking the average over the last 30 s of each step; MD stands for missing data.

<table><caption>**Table 4** StO₂ (%) for each subject at each step. CMOS-DRS method A.</caption><thead><tr><th>Subject</th><th>A</th><th>B</th><th>C</th><th>D</th><th>E</th><th>F</th><th>G</th><th>H</th><th>I</th><th>J</th><th>K</th></tr></thead><tbody><tr><td>01</td><td>63.7</td><td>51.5</td><td>69.9</td><td>69.2</td><td>70.1</td><td>72.9</td><td>72.8</td><td>75.7</td><td>52.0</td><td>82.8</td><td>76.4</td></tr><tr><td>02</td><td>68.6</td><td>65.7</td><td>64.6</td><td>66.5</td><td>66.9</td><td>66.8</td><td>66.6</td><td>68.0</td><td>63.5</td><td>69.3</td><td>66.6</td></tr><tr><td>03</td><td>72.1</td><td>84.3</td><td>71.1</td><td>71.2</td><td>81.2</td><td>66.8</td><td>65.6</td><td>61.0</td><td>42.7</td><td>91.6</td><td>90.5</td></tr><tr><td>04</td><td>78.7</td><td>87.4</td><td>77.9</td><td>72.5</td><td>73.6</td><td>71.4</td><td>71.5</td><td>72.2</td><td>58.6</td><td>99.8</td><td>72.0</td></tr><tr><td>05</td><td>74.9</td><td>72.3</td><td>70.8</td><td>72.3</td><td>70.6</td><td>70.2</td><td>72.6</td><td>69.9</td><td>57.9</td><td>72.8</td><td>71.4</td></tr><tr><td>06</td><td>68.4</td><td>66.7</td><td>65.2</td><td>67.7</td><td>68.4</td><td>68.6</td><td>68.6</td><td>68.3</td><td>55.0</td><td>68.5</td><td>67.9</td></tr><tr><td>07</td><td>66.1</td><td>67.9</td><td>73.9</td><td>74.8</td><td>72.9</td><td>74.3</td><td>77.2</td><td>73.8</td><td>71.5</td><td>88.1</td><td>78.4</td></tr><tr><td>08</td><td>77.7</td><td>77.0</td><td>75.8</td><td>77.4</td><td>77.2</td><td>77.3</td><td>77.8</td><td>77.3</td><td>67.3</td><td>77.9</td><td>77.5</td></tr><tr><td>09</td><td>75.1</td><td>74.4</td><td>79.6</td><td>76.3</td><td>80.5</td><td>77.4</td><td>75.8</td><td>75.5</td><td><strong>0.0</strong></td><td>81.1</td><td>72.7</td></tr><tr><td>10</td><td>76.0</td><td>71.9</td><td>64.9</td><td>72.0</td><td>72.0</td><td>73.1</td><td>63.2</td><td>72.4</td><td>40.7</td><td>74.5</td><td>72.3</td></tr><tr><td>11</td><td>61.3</td><td>58.5</td><td>57.0</td><td>66.0</td><td>66.4</td><td>66.6</td><td>66.4</td><td>66.7</td><td>47.6</td><td>66.8</td><td>59.6</td></tr><tr><td>12</td><td>74.3</td><td>73.1</td><td>72.2</td><td>72.6</td><td>72.8</td><td>72.7</td><td>73.2</td><td>72.0</td><td>66.3</td><td>72.8</td><td>72.5</td></tr><tr><td>13</td><td>92.1</td><td>91.8</td><td>89.6</td><td>91.4</td><td>91.1</td><td>91.3</td><td>92.0</td><td>90.8</td><td>44.7</td><td>91.5</td><td>90.9</td></tr><tr><td>14</td><td>85.2</td><td>88.4</td><td>78.7</td><td>76.7</td><td>82.7</td><td>87.8</td><td>85.6</td><td>83.1</td><td>41.5</td><td>90.6</td><td>81.0</td></tr><tr><td>15</td><td>69.9</td><td>66.1</td><td>52.1</td><td>75.2</td><td>80.2</td><td>82.1</td><td>76.4</td><td>57.8</td><td>34.0</td><td>72.0</td><td>70.7</td></tr><tr><td>16</td><td><strong>MD</strong></td><td>72.2</td><td>76.7</td><td>72.9</td><td>81.4</td><td>79.5</td><td>79.3</td><td>82.3</td><td>71.9</td><td>87.2</td><td>73.0</td></tr><tr><td>17</td><td>73.2</td><td>78.2</td><td>79.7</td><td>77.0</td><td>74.6</td><td>75.6</td><td>76.0</td><td>72.0</td><td>53.8</td><td>81.3</td><td>65.3</td></tr><tr><td>18</td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td><td><strong>MD</strong></td></tr><tr><td>19</td><td>91.6</td><td>69.6</td><td>73.8</td><td>70.6</td><td>80.8</td><td>84.2</td><td>76.9</td><td>76.0</td><td>51.6</td><td>88.3</td><td>86.7</td></tr><tr><td>20</td><td>67.7</td><td>76.1</td><td>74.5</td><td>74.0</td><td>83.0</td><td>87.1</td><td>71.7</td><td><strong>MD</strong></td><td>62.1</td><td>89.4</td><td>78.1</td></tr></tbody></table>
<table><caption>**Table 5** StO₂ (%) for each subject at each step. CMOS-DRS method B.</caption><thead><tr><th>Subject</th><th>A</th><th>B</th><th>C</th><th>D</th><th>E</th><th>F</th><th>G</th><th>H</th><th>I</th><th>J</th><th>K</th></tr></thead><tbody><tr><td>01</td><td>77.3</td><td>76.0</td><td>66.3</td><td>74.4</td><td>74.3</td><td>78.2</td><td>78.3</td><td>74.0</td><td>73.8</td><td>89.3</td><td>83.0</td></tr><tr><td>02</td><td>77.6</td><td>71.7</td><td>68.4</td><td>72.6</td><td>73.6</td><td>73.6</td><td>72.1</td><td>74.6</td><td>63.5</td><td>81.0</td><td>70.2</td></tr><tr><td>03</td><td>79.5</td><td>74.1</td><td>77.7</td><td>77.6</td><td>77.4</td><td>75.4</td><td>80.0</td><td>79.8</td><td>42.4</td><td>81.9</td><td>62.1</td></tr><tr><td>04</td><td>30.4</td><td>32.7</td><td>42.1</td><td>43.0</td><td>45.3</td><td>43.7</td><td>44.1</td><td>40.5</td><td>6.1</td><td>59.3</td><td>40.9</td></tr><tr><td>05</td><td>61.8</td><td>55.5</td><td>54.2</td><td>65.6</td><td>65.0</td><td>65.7</td><td>65.7</td><td>63.1</td><td>51.4</td><td>67.1</td><td>60.1</td></tr><tr><td>06</td><td>50.7</td><td>49.2</td><td>45.8</td><td>52.5</td><td>57.2</td><td>56.9</td><td>54.5</td><td>56.7</td><td>36.6</td><td>63.4</td><td>58.9</td></tr><tr><td>07</td><td>57.4</td><td>58.0</td><td>57.1</td><td>58.1</td><td>65.5</td><td>60.5</td><td>55.7</td><td>56.3</td><td>46.1</td><td>61.2</td><td>55.6</td></tr><tr><td>08</td><td>86.2</td><td>76.3</td><td>71.7</td><td>76.8</td><td>79.3</td><td>80.1</td><td>76.7</td><td>78.3</td><td>48.7</td><td>84.0</td><td>76.4</td></tr><tr><td>09</td><td>67.5</td><td>66.9</td><td>55.7</td><td>64.5</td><td>61.7</td><td>66.6</td><td>67.8</td><td>62.3</td><td>50.4</td><td>86.3</td><td>67.3</td></tr><tr><td>10</td><td>76.9</td><td>76.4</td><td>76.0</td><td>77.2</td><td>74.1</td><td>73.4</td><td>82.1</td><td>72.9</td><td>58.3</td><td>84.2</td><td>73.4</td></tr><tr><td>11</td><td>58.3</td><td>56.9</td><td>61.4</td><td>63.0</td><td>64.3</td><td>64.5</td><td>64.5</td><td>66.3</td><td>60.9</td><td>77.3</td><td>66.1</td></tr><tr><td>12</td><td>62.7</td><td>61.1</td><td>61.3</td><td>65.0</td><td>66.0</td><td>66.7</td><td>66.3</td><td>65.8</td><td>44.4</td><td>69.4</td><td>66.9</td></tr><tr><td>13</td><td>70.5</td><td>65.5</td><td>62.6</td><td>64.0</td><td>68.2</td><td>68.1</td><td>68.9</td><td>66.3</td><td>58.6</td><td>79.3</td><td>67.0</td></tr><tr><td>14</td><td>66.8</td><td>65.0</td><td>65.3</td><td>69.4</td><td>68.8</td><td>68.7</td><td>70.5</td><td>68.8</td><td>53.6</td><td>75.5</td><td>68.3</td></tr><tr><td>15</td><td>74.6</td><td>74.8</td><td>68.3</td><td>74.1</td><td>77.1</td><td>78.1</td><td>78.5</td><td>75.2</td><td>35.5</td><td>78.5</td><td>77.2</td></tr><tr><td>16</td><td>MD</td><td>67.1</td><td>67.1</td><td>67.4</td><td>62.3</td><td>63.1</td><td>66.3</td><td>59.9</td><td>55.2</td><td>72.6</td><td>68.4</td></tr><tr><td>17</td><td>65.3</td><td>63.0</td><td>61.8</td><td>73.7</td><td>78.0</td><td>78.4</td><td>74.1</td><td>75.9</td><td>57.2</td><td>83.3</td><td>76.1</td></tr><tr><td>18</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td></tr><tr><td>19</td><td>58.5</td><td>64.2</td><td>62.7</td><td>64.4</td><td>62.9</td><td>60.9</td><td>62.5</td><td>65.6</td><td>59.2</td><td>71.1</td><td>60.1</td></tr><tr><td>20</td><td>62.8</td><td>57.4</td><td>58.6</td><td>58.9</td><td>61.6</td><td>56.4</td><td>63.6</td><td>59.2</td><td>49.6</td><td>67.1</td><td>63.7</td></tr></tbody></table>

<table><caption>**Table 6** StO₂ (%) for each subject at each step. Artinis PortaLite.</caption><thead><tr><th>Subject</th><th>A</th><th>B</th><th>C</th><th>D</th><th>E</th><th>F</th><th>G</th><th>H</th><th>I</th><th>J</th><th>K</th></tr></thead><tbody><tr><td>01</td><td>67.6</td><td>65.4</td><td>63.2</td><td>63.2</td><td>68.4</td><td>66.7</td><td>66.9</td><td>66.9</td><td>50.7</td><td>72.3</td><td>68</td></tr><tr><td>02</td><td>66.8</td><td>63.7</td><td>62.2</td><td>65.4</td><td>65.9</td><td>66.7</td><td>69.5</td><td>65.3</td><td>51.4</td><td>77.0</td><td>67.2</td></tr><tr><td>03</td><td>74.1</td><td>73.6</td><td>70.3</td><td>74</td><td>76</td><td>75.5</td><td>76.4</td><td>74</td><td>62.7</td><td>80.8</td><td>74.3</td></tr><tr><td>04</td><td>71.3</td><td>70.1</td><td>69.7</td><td>71.5</td><td>71.7</td><td>73.2</td><td>72.7</td><td>70.7</td><td>67.9</td><td>75</td><td>71.3</td></tr><tr><td>05</td><td>72.8</td><td>71.3</td><td>70</td><td>71.4</td><td>71.9</td><td>71.8</td><td>71.2</td><td>71</td><td>62.8</td><td>76.4</td><td>70.5</td></tr><tr><td>06</td><td>68.6</td><td>66.3</td><td>64.8</td><td>66.1</td><td>69.2</td><td>68.9</td><td>69.9</td><td>66.2</td><td>51.8</td><td>77.6</td><td>68.3</td></tr><tr><td>07</td><td>74.7</td><td>72</td><td>68.2</td><td>72.3</td><td>73.7</td><td>74.1</td><td>74.1</td><td>75.1</td><td>61.8</td><td>80.8</td><td>74.8</td></tr><tr><td>08</td><td>75</td><td>75</td><td>74</td><td>76</td><td>79</td><td>79</td><td>76</td><td>78</td><td>70</td><td>79.6</td><td>74</td></tr><tr><td>09</td><td>82</td><td>75</td><td>75</td><td>78</td><td>80</td><td>79</td><td>78</td><td>80</td><td>45</td><td>85.7</td><td>75</td></tr><tr><td>10</td><td>75</td><td>75</td><td>71</td><td>75</td><td>76</td><td>80</td><td>80</td><td>74</td><td>58</td><td>85.9</td><td>76</td></tr></tbody></table>
<table><caption>**Table 6 (Continued).**</caption><thead><tr><th>Subject</th><th>A</th><th>B</th><th>C</th><th>D</th><th>E</th><th>F</th><th>G</th><th>H</th><th>I</th><th>J</th><th>K</th></tr></thead><tbody><tr><td>11</td><td>63</td><td>65</td><td>61</td><td>66</td><td>65</td><td>64</td><td>65</td><td>67</td><td>43</td><td>85.1</td><td>68</td></tr><tr><td>12</td><td>66</td><td>64</td><td>62</td><td>64</td><td>72</td><td>70</td><td>69</td><td>71</td><td>49</td><td>94.5</td><td>75</td></tr><tr><td>13</td><td>94</td><td>86</td><td>85</td><td>88</td><td>92</td><td>89</td><td>95</td><td>92</td><td>56</td><td>96.7</td><td>92</td></tr><tr><td>14</td><td>77</td><td>74</td><td>72</td><td>74</td><td>78</td><td>78</td><td>76</td><td>77</td><td>53</td><td>86.3</td><td>75</td></tr><tr><td>15</td><td>75</td><td>70</td><td>65</td><td>60</td><td>72</td><td>70</td><td>70</td><td>67</td><td>27</td><td>94.47</td><td>70</td></tr><tr><td>16</td><td>MD</td><td>78</td><td>77</td><td>77</td><td>78</td><td>79</td><td>77</td><td>80</td><td>65</td><td>86.41</td><td>79</td></tr><tr><td>17</td><td>77</td><td>77</td><td>75</td><td>77</td><td>75</td><td>76</td><td>78</td><td>76</td><td>61</td><td>84.78</td><td>77</td></tr><tr><td>18</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td><td>MD</td></tr><tr><td>19</td><td>74</td><td>73</td><td>74</td><td>74</td><td>77</td><td>79</td><td>77</td><td>77</td><td>65</td><td>80.27</td><td>76</td></tr><tr><td>20</td><td>74</td><td>73</td><td>72</td><td>72</td><td>71</td><td>75</td><td>73</td><td>MD</td><td>64</td><td>79.56</td><td>72</td></tr></tbody></table>

## Disclosures

The authors have no relevant financial interests in this article and no potential conflicts of interest to disclose. Nils Petitdidier and Anne Koenig are coinventors of a patent devoted to the srDRS system architecture presented in this manuscript. All rights on this patent are owned by the Commissariat à l'Energie Atomique et aux Energies Alternatives (CEA).

## Acknowledgments

This work was supported by the Commissariat à l'Energie Atomique et aux Energies Alternatives (CEA). The Capnometry PtCO₂ sensor research is funded by L3 Medical Society. Part of the clinical trial has been funded by L3 Medical Society.

## References

1. J. M. Smit et al., "Advancements in free flap monitoring in the last decade: a critical review," *Plast. Reconstr. Surg.* **125**(1), 177-185 (2010).

2. D. T. Delpy and M. Cope, "Quantification in tissue near-infrared spectroscopy," *Philos. Trans. R. Soc. Lond. B* **352**(1354), 649-659 (1997).

3. N. F. Jones, "Intraoperative and postoperative monitoring of microsurgical free tissue transfers," *Clin. Plast. Surg.* **19**(4), 783-797 (1992).

4. T. J. Farrell, M. S. Patterson, and B. Wilson, "A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo," *Med. Phys.* **19**(4), 879-888 (1992).

5. A. Koenig et al., "Diffuse reflectance spectroscopy: a clinical study of tuberculin skin tests reading," *Proc. SPIE* **8592**, 85920S (2013).

6. O. Senlik and N. M. Jokerst, "Concentric multipixel silicon photodiode array probes for spatially resolved diffuse reflectance spectroscopy," *IEEE J. Sel. Top. Quantum Electron.* **22**(3), 7-12 (2016).

7. N. Petitdidier et al., "Contact, high-resolution spatial diffuse reflectance imaging system for skin condition diagnosis," *J. Biomed. Opt.* **23**(11), 115003 (2018).

8. N. Petitdidier, "Imaging system for the characterization of skin layers using diffuse reflectance (Système d'imagerie pour la caractérisation en couches de la peau par réflectance diffuse)," PhD Thesis, Université de Strasbourg (2018).
9. Q. Liu and N. Ramanujam, "Scaling method for fast Monte Carlo simulation of diffuse reflectance spectra from multilayered turbid media," *J. Opt. Soc. Am. A* **24**(4), 1011 (2007).

10. N. Rajaram, T. H. Nguyen, and J. W. Tunnell, "Lookup table-based inverse model for determining optical properties of turbid media," *J. Biomed. Opt.* **13**(5), 050501 (2008).

11. B. Roig et al., "Multilayered phantoms with tunable optical properties for a better understanding of light/tissue interactions," *Proc. SPIE* **9325**, 93250B (2015).

12. R. B. Saager et al., "Multilayer silicone phantoms for the evaluation of quantitative optical techniques in skin imaging," *Proc. SPIE* **7567**, 756706 (2010).

13. G. J. Greening et al., "Characterization of thin poly(dimethylsiloxane)-based tissue-simulating phantoms with tunable reduced scattering and absorption coefficients at visible and near-infrared wavelengths," *J. Biomed. Opt.* **19**(11), 115002 (2014).

14. D. M. de Bruin et al., "Optical phantoms of varying geometry based on thin building blocks with controlled optical properties," *J. Biomed. Opt.* **15**(2), 025001 (2010).

15. J. E. Bender et al., "Noninvasive monitoring of tissue hemoglobin using UV-VIS diffuse reflectance spectroscopy: a pilot study," *Opt. Express* **17**(26), 23396-23409 (2009).

16. V. Sorgato et al., "ACA-pro: calibration protocol for quantitative diffuse reflectance spectroscopy. Validation on contact and noncontact probe- and CCD-based systems," *J. Biomed. Opt.* **21**(6), 065003 (2016).

17. S. L. Jacques, "Quick analysis of optical spectra to quantify epidermal melanin and papillary dermal blood content of skin," *J. Biophotonics* **8**(4), 309-316 (2015).

18. I. V. Meglinski and S. J. Matcher, "Computer simulation of the skin reflectance spectra," *Comput. Methods Programs Biomed.* **70**(2), 179-186 (2003).

19. S. L. Jacques, "Optical properties of biological tissues: a review," *Phys. Med. Biol.* **58**(11), R37 (2013).

20. Artinis Medical Systems, "Artinis Medical Systems | fNIRS devices | NIRS devices-PortaLite," Artinis Medical Systems | fNIRS devices | NIRS devices, 2019, https://www.artinis.com/portalite (accessed 26 November 2019).

21. W. van Weteringen et al., "Novel transcutaneous sensor combining optical $tcPO_2$ and electrochemical $tcPCO_2$ monitoring with reflectance pulse oximetry," *Med. Biol. Eng. Comput.* **58**(2), 239-247 (2020).

22. Sports & Medical Technologies SA, "AltiTrainer Hypoxia, Hyperoxia, Hypercapnia," Sports & Medical Technologies S.A., 2019, https://www.smtec.net/en/products/altitrainer/ (accessed 26 November 2019).

23. P. Grangeat et al., "Evaluation in healthy subjects of a transcutaneous carbon dioxide monitoring wristband during hypo and hypercapnia conditions," in *42nd Annu. Int. Conf. IEEE Eng. Med. and Biol. Soc.*, Montreal, Canada (2020).

24. M. C. van der Sluijs et al., "New and highly sensitive continuous-wave near-infrared spectrophotometer with multiple detectors," *Proc. SPIE* **3194**, 63-72 (1998).

25. P. Grangeat et al., "First evaluation of a transcutaneous carbon dioxide monitoring wristband device during a cardiopulmonary exercise test," in *41st Annu. Int. Conf. IEEE Eng. Med. and Biol. Soc.*, Berlin, Germany, pp. 3352-3355 (2019).

26. SenTec Inc., "SenTec V-SignTM System—SenTec," SenTec products, 2019, https://www.sentec.com/products/sentec-v-signtm-system/ (accessed 26 November 2019).

27. T. Rupp et al., "Cerebral and muscle oxygenation during intermittent hypoxia exposure in healthy humans," *Sleep* **39**(6), 1197-1199 (2016).

28. T. Rupp et al., "Tissue deoxygenation kinetics induced by prolonged hypoxic exposure in healthy humans at rest," *J. Biomed. Opt.* **18**(9), 095002 (2013).

**Anne Koenig** graduated as an automation and computer engineer in 1986 then joined CEA-LETI where she received her PhD in artificial intelligence for robotics. She got the accreditation to supervise research degrees (HDR) in 2015. Until 2005, her research has been focused on radiographic image processing for nondestructive testing and 3D dynamic CT, then in the field of fluorescence diffuse optical tomography. She is now involved in skin characterization by diffuse reflectance spectroscopy and wearable sensors.
**Nils Petitdidier** graduated as an optical engineer from IOGS, Institut d'Optique Graduate School, Palaiseau, France, and received his MS degree in engineering physics from KTH, Kungliga Tekniska Högskolan, Stockholm, Sweden, in 2015. He received his PhD in imaging physics from the Université de Strasbourg, Strasbourg, France, in 2018. His research interests are focused on the development of wearable, low-cost semiconductor spectral imaging systems for skin condition diagnosis.

Biographies of the other authors are not available.