# Proposal for a skin layer-wise decomposition model of spatially-resolved diffuse reflectance spectra based on maximum depth photon distributions: A numerical study

Victor Colas, Walter Blondel, Grégoire Khairallah, Christian Daul, Marine Amouroux

► To cite this version:

Victor Colas, Walter Blondel, Grégoire Khairallah, Christian Daul, Marine Amouroux. Proposal for a skin layer-wise decomposition model of spatially-resolved diffuse reflectance spectra based on maximum depth photon distributions: A numerical study. Photonics, 2021, 8 (10), pp.444. 10.3390/photon-ics8100444. [https://arxiv.org/abs/2108.04447](https://arxiv.org/abs/2108.04447)

HAL Id: hal-03377401

https://hal.science/hal-03377401

Submitted on 25 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Proposal for a Skin Layer-Wise Decomposition Model of Spatially-Resolved Diffuse Reflectance Spectra Based on Maximum Depth Photon Distributions: A Numerical Study

Victor Colas <fnref n="1" />*, Walter Blondel <fnref n="1" />, Grégoire Khairallah <fnref n="1" />,<fnref n="2" />, Christian Daul <fnref n="1" /> and Marine Amouroux <fnref n="1" />

**Abstract:** In the context of cutaneous carcinoma diagnosis based on in vivo optical biopsy, Diffuse Reflectance (DR) spectra, acquired using a Spatially Resolved (SR) sensor configuration, can be analyzed to distinguish healthy from pathological tissues. The present contribution aims at studying the depth distribution of SR-DR-detected photons in skin from the perspective of analyzing how these photons contribute to acquired spectra carrying local physiological and morphological information. Simulations based on modified *Cuda Monte Carlo Modeling of Light* transport were performed on a five-layer human skin optical model with epidermal thickness, prototype and dermal blood content as variable parameters using (i) wavelength-resolved scattering and absorption properties and (ii) the geometrical configuration of a multi-optical fiber probe implemented on an SR-DR spectroscopic device currently used in clinics. Through histograms of the maximum probed depth and their exploitation, we provide numerical evidence linking the characteristic penetration depth of the detected photons to their wavelengths and four source-sensor distances, which made it possible to propose a decomposition of the DR signals related to skin layer contributions.

**Keywords:** Diffuse Reflectance Spectroscopy; probed depth; skin layer; Monte Carlo Simulation; skin optics

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

check for
updates

**Citation:** Colas, V.; Blondel, W.; Khairallah, G.; Daul, C.; Amouroux, M. Proposal for a Skin Layer-Wise Decomposition Model of Spatially-Resolved Diffuse Reflectance Spectra Based on Maximum Depth Photon Distributions: A Numerical Study. *Photonics* 2021, 8, 444. https://doi.org/10.3390/photonics8100444

Received: 27 August 2021
Accepted: 12 October 2021
Published: 14 October 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Diffuse Reflectance Spectroscopy (DRS) applied to biological tissues is a widely studied [1–4], non-destructive optical characterization technique. It consists of measuring back-reflected intensity spectra carrying information from light-tissue interactions, from which the optical and structural properties of the probed medium can then be extracted and analyzed. By illuminating the tissue with a broadband light, photons propagating in the medium undergo wavelength-dependent absorption and scattering events, before some of them may move back towards the surface and potentially be captured by a spectrometric detection system. An analysis of the amplitude and shape features of the latter intensity spectrum thus gives access to information of interest such as absorption bands and elastic scattering properties related to tissue-specific chromophores and structures. Its non-invasive nature, its sensitivity to metabolic and structural changes in tissues and its clinical applicability are very well suited to in vivo biomedical diagnostics.

Spatially Resolved DRS (SR-DRS) refers to fibered systems of DRS featuring several Detection Fibers (DF) located at different distances from the Source Fiber (SF) allowing for depth-resolved acquisitions. Indeed, photons detected at a short/large source-detector (SD) distance will be associated with shallower/deeper trajectories into the biological tissues under evaluation. The literature on this subject points to numerous applications

<footnote n="1">¹ CRAN, BioSiS Department, Université de Lorraine-CNRS, F-54000 Nancy, France; walter.blondel@univ-lorraine.fr (W.B.); gregoire.khairallah@univ-lorraine.fr (G.K.); Christian.Daul@univ-lorraine.fr (C.D.); marine.amouroux@univ-lorraine.fr (M.A.)</footnote>

<footnote n="2">² Department of Plastic, Aesthetic and Reconstructive Surgery, Metz-Thionville Regional Hospital, F-57530 Ars-Laquenexy, France</footnote>

<footnote>* Correspondence: victor.colas@univ-lorraine.fr</footnote>
of SR-DRS such as the measurement of blood glucose levels [5,6], hemoglobin oxygen saturation [7] and for pre-cancer and cancer diagnosis: lung [8], cervix [9], colon [10], bladder [11], brain [12], breast [13], oesophagus [14] or pancreas [15]. There are also a large number of applications of SR-DRS for the in vivo characterization of cutaneous pathological states, including melanoma [16–18] and carcinoma [19,20], and more generally for the characterization of skin layer optical properties [21,22].

In this context of skin optical biopsy, point spectroscopy devices exploiting light-tissue interactions have been designed and used in several clinical studies [11,23,24]. The estimation of diagnostically relevant optical parameters from in vivo SR-DR spectra, and especially their evolution around the lesion, can be linked to the tissue pathological state. As the metabolic modifications associated with carcinogenesis are local, the quantitative estimation of the penetration depth of the photons detected using point spectroscopy is at stake, as reported by Breslin et al. [25], who studied the sensitivity of DR spectra to elastosis (accumulation of elastotic material in the papillary dermis located underneath epidermis). Cancer growth is also related to the enzymatic degradation of collagen fibers in the dermis leading to collagen fluorescence decrease [26] and to the decrease of Nicotinamide Adenine Dinucleotide Hydrogen (NADH) fluorescence in the epidermis [27]. Finally, morphological changes may be significant for part of skin carcinomas such as ulceration and hyperkeratosis, which are respectively characterized by a thinning or skin thickening (~×10) of the stratum corneum superficial layer.

The need to explain DR signals acquired on the surface, resulting from light-tissue interactions in all skin layers, justifies the use of numerical simulation as a tool to interpret the acquired signals. The literature reports work by several teams using photon transport simulation in a multi-layer skin model to extract features that describe the propagation of the acquired photons [28,29]. First, in a very visual way, absorption density maps of detected photons in a three-layer skin model [30] provide information on a photon's trajectories. In a more quantitative manner, teams have presented their results on the penetration depth of photons. In a two-layer model, histograms of probed depths have been created at 450 and 650 nm to extract the average penetration depth [31], while the interrogation depth (mean depth of all detected photon scattering events) [32] or measurement depth (median of the cumulative absorption on fluence absorption maps) [30] can also be extracted from simulation to characterize the photon's penetration depth. Other quantities of interest are exploited to know more about the subcutaneous behavior of backscattered photons, such as the path length distribution with increasing SD distance [33,34]. This quantity was even declined in the mean partial path length [35,36], which consists of dividing the total path length into layer path contributions.

Although these approaches have made it possible to study the behavior of photons (path, depth) in different multi-layer skin models, none of them, to our knowledge, has led to the analytical formulation of an explicit decomposition of the “global” DR spectra, as a function of the SD distances and of the modeled skin layers. This is an essential issue to (i) directly interpret the measured spectra, especially in the case of clinical application and patients’ inter-individual variability in terms of epidermal thicknesses, skin type (healthy, carcinoma), phenotype and blood dermal content and to (ii) indirectly estimate the optical properties of the different layers in an efficient way. In order to address this gap, the present contribution deals with the acquisition of simulated DR spectra for a multi-layer medium considering (i) the absorption and scattering of the five main endogenous chromophores (ii) on the near UV-visible range and (iii) for a geometrical model of the probe adapted to the previous wavelengths, namely SD distances smaller than a millimeter. Histograms of the maximal probed depth for photons detected at several SD distances were created and analyzed to quantify light penetration in the media. The study was led for skin media with varying parameters: epidermal thickness, phenotype, dermal blood concentration and health status of the tissue (ulcerated and hyperkeratosed carcinomas). Finally, a DR signal decomposition according to the deepest layer probed is offered.
## 2. Materials and Methods

Section 2.1 describes the geometrical and spectral configuration of the probe, Section 2.2 provides an explanation of the numerical skin model developed, and Section 2.3 refers to the simulation description.

### 2.1. SR-DRS Probe—Geometrical and Spectral Features Implemented for Simulations

The geometrical and optical features of source and detection fibers implemented in our SR-DRS simulations were based on the technical characteristics of an existing device developed by our team [37], whose purpose is to assist a surgeon in establishing surgical margins during skin carcinoma resection [38]. Work with a similar goal has already been presented in the literature using DR [39] (elastic scattering) or Raman [40] (inelastic scattering) spectra. The multiple fiber optical probe of the clinical device consists of (i) one central 600 µm-core diameter Source Fiber (SF) and numerical aperture $NA_{SF} = 0.22$ delivering broadband 365–765 nm white light illumination to the tissue in contact and (ii) 24,200 µm core diameter detection fibers (DF) of $NA_{DF} = 0.22$ each, equally distributed on 4 concentric circles corresponding to the 4 SD center-to-center distances $D_n$ ($n \in \{1, ..., 4\}$) from 400 to 1000 µm by 200 µm steps. The photons reaching each of the six collection fibers of the $D_n$ distance are merged to constitute the diffuse reflectance signal $DR_{D_n}(\lambda)$. To maximize the signal to noise ratio in the simulation, cylindrical symmetry was exploited using a ring detection surface before multiplying the quantity by a factor that ensures the geometrical correspondence with the clinical device. Complementary simulations were conducted to quantify the error related to this geometric approximation (detailed results not shown). In the “worst” case—i.e., for small SD distances—the relative error was less than 0.5%. A schematic representation of the acquisition geometry is shown in Figure 1, while the numerical values of the different parameters used are summarized in Table 1. The optical probe geometrical configuration considered in our simulation does not take into account the epoxy resin and the polished stainless steel in which optical fibers (whose core and cladding are made of silica) are embedded in the clinical device. As mentioned by Zelinskyi et al., this may limit the accuracy of the results obtained [41].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Schematic representation of the SR-DRS fiber probe geometry and of the five-layer skin implemented for simulations: vertical section view (**left**) and top view (**right**). SF, DF and $DR_{D_n}(\lambda)$ stand for Source Fiber, Detection Fiber and Diffuse Reflectance spectrum at the distance $D_n$ ($n \in \{1, 2, 3, 4\}$), respectively.</figcaption></figure>
<table><caption>**Table 1.** Values of the geometrical and spectral source-detector features for the SR-DRS simulation. SF: source fiber; DF: detection fiber; Dn: SD distance; rSF/DF: fiber radius.</caption><thead><tr><th>[λ<sub>min</sub>, λ<sub>max</sub>] (nm)</th><th>D<sub>n</sub> (μm)</th><th>NA<sub>SF/DF</sub></th><th>r<sub>SF</sub> (μm)</th><th>r<sub>DF</sub> (μm)</th></tr></thead><tbody><tr><td>[365, 765]</td><td>400, 600, 800, 1000</td><td>0.22</td><td>300</td><td>100</td></tr></tbody></table>

## 2.2. Skin Five-Layered Model

### 2.2.1. Layer Thicknesses and Skin Types Modeled

The numerical model mimicking the skin was designed following a five planar multi-layered medium, representing (from superficial to deepest layers) the stratum corneum (1-SC), the living epidermis (2-LE), the papillary dermis (3-PD), the reticular dermis (4-RD) and the subcutaneous fat (5-SF). The total thickness was fixed to z<sup>tot</sup> = 6120 µm, while the total radius of the simulation medium was set to r<sup>tot</sup> = 3000 µm. The results presented hereafter are valid under the assumption that the optical properties of each of the skin layers are homogeneous in the cylinder thus described (cf. Figure 1). Considering the order of magnitude of skin optical properties and the position of the farthest collecting fiber (D<sub>4</sub> = 1000 µm), edge effects (i.e., photons that would be collected by DF after reflection on the side or bottom border of the medium) were negligible. This was verified numerically with absorption maps of detected photons: none of them reached the bottom (in the worst case—i.e., a wavelength in near IR, great SD distance and fair prototype—only 1% of detected photons were scattered in the SF), while less than 0.1% were collected after having reached the side border.

Consequently, five combinations of SC and LE thicknesses were defined, among which three featured a thin, medium and thick healthy epidermis and two represented some specific values encountered in carcinoma. SC and LE thicknesses may significantly vary according to age, gender, and body location [42–44] and also for different skin physiological states. For instance, the stratum corneum can either completely disappear or thicken 10-fold [45] in the case of ulceration or hyperkeratosis, respectively. Indeed, other epidermal layers are also affected by cell proliferation due to skin carcinoma. However, for the sake of comparison with healthy conditions, only SC thickness was varied to model these two sub-types of skin carcinomas; i.e., in the case of ulceration, the absence of SC was modeled through the removal this superficial layer. As shown in Figure 1, the thickness values of middle layers PD and RD were fixed (200 and 1800 µm, respectively), while those of the top layers (SC and LE) were varied to simulate different tissue conditions. The thickness of the underlying SF layer was consequently adjusted to keep the total thickness z<sup>tot</sup> of our simulation model constant. We considered the range of values of SC and LE thicknesses as well as their approximate ratio (~1/4) reported for sun-exposed skin sites only (head including face, scalp, ear, cheek and shoulders, back, arms, back of hand and legs), which are 10–30 µm and 60–120 µm for SC and LE layers, respectively. The numerical values associated with the elaboration of these three healthy and two other pathological epithelium (carcinoma) simulation media configurations are summarized in Table 2.

<table><caption>**Table 2.** Values of the SC (stratum corneum) and LE (living epidermis) layer thicknesses implemented in our skin model to simulate five different healthy skin and carcinoma conditions.</caption><thead><tr><th colspan="2">Skin Categories</th><th>SC Thickness z<sup>1</sup> (μm)</th><th>LE Thickness z<sup>2</sup> (μm)</th></tr></thead><tbody><tr><td rowspan="3">Healthy</td><td>Smallest thicknesses</td><td>10</td><td>60</td></tr><tr><td>Mean thicknesses</td><td>20</td><td>90</td></tr><tr><td>Highest thicknesses</td><td>30</td><td>120</td></tr><tr><td rowspan="2">Carcinoma</td><td>Ulcerated</td><td>0</td><td>90</td></tr><tr><td>Hyperkeratosed</td><td>300</td><td>90</td></tr></tbody></table>
### 2.2.2. Skin Optical Properties

The optical properties defined for each skin layer in terms of absorption and scattering are synthesized in the following paragraphs. Complementary details can also be found in [46].

#### Absorption Properties

Absorption coefficients $\mu_a^k(\lambda)$ were defined for each layer $k \in \{1, ..., 5\}$ of the numerical skin model taking into account the contribution of the five main endogenous absorbers i of human skin, which are blood (decomposed in oxygenated hemoglobin $HbO_2$, deoxygenated hemoglobin $Hb$, and bilirubin contributions), water, melanin (decomposed in eumelanin and pheomelanin), lipids and $\beta$-caroten, respectively, as well as their volume fraction $f_{v,i}^k$, in the mixture. A residual absorption baseline was also added according to [47] (cf. $\mu_{a,base}$ in Equation (1)). Thus, the total absorption coefficient of a skin layer was defined as the sum of the latter baseline and the contributions $\mu_{a,i}^k$ of each of the five chromophores, weighted by their respective volume fraction, in accordance with [48]:

$$ \begin{aligned} \mu_a^k(\lambda) = & \mu_{a,base}(\lambda) \\ & + f_{v,blood}^k[\gamma\mu_{a,HbO_2}^k(\lambda) + (1-\gamma)\mu_{a,Hb}^k(\lambda) + \mu_{a,bilirubin}^k(\lambda)] \\ & + f_{v,water}^k\mu_{a,water}^k(\lambda) \\ & + f_{v,mel}^k[\mu_{a,eum}^k(\lambda) + \mu_{a,phe}^k(\lambda)] \\ & + f_{v,lipids}^k\mu_{a,lipids}^k(\lambda) \\ & + f_{v,\beta-caroten}^k\mu_{a,\beta-caroten}^k(\lambda) \end{aligned} \quad (1) $$

where $\mu_{a,i}^k$ are all expressed in cm$^{-1}$, $\mu_{a,base}(\lambda) = 7.84 \times 10^8 \times \lambda^{-3225}$ ($\lambda$ in nm), and $\gamma$ is the ratio between $HbO_2$ and $Hb$, set to $\gamma = 0.9$ in our simulations. Lipids and water absorption spectra as well as extinction coefficients of the other compounds were directly implemented in our model from [49]. It is noticeable that, for $\beta$-caroten and blood (namely both hemoglobins and bilirubin), molar extinction coefficient $\epsilon_i^{mol}(\lambda)$ (in cm$^{-1}$/mol·L$^{-1}$) was given, while for melanin, the mass extinction coefficients $\epsilon_i^{mass}(\lambda)$ (in cm$^{-1}$/g·L$^{-1}$) of eu- and pheo-melanin were provided. These are linked to absorption coefficients according to the following formulas:

$$ \mu_{a,i}^k(\lambda) = \frac{\epsilon_i^{mol}(\lambda)}{M_i} c_i^k \qquad \mu_{a,i}^k(\lambda) = \epsilon_i^{mass}(\lambda) c_i^k \quad (2) $$

where $c_i^k$ is the mass concentration (in g/L) of chromophore i in the mixing model and $M_i$ the respective molar weight with the following values implemented in our simulations:

$$ M_{\beta-Caroten} = 587 \text{ g} \cdot \text{mol}^{-1} \qquad M_{\text{bilirubin}} = 585 \text{ g} \cdot \text{mol}^{-1} \qquad M_{HbO_2,Hb} = 66,500 \text{ g} \cdot \text{mol}^{-1} $$

In order to take into account skin pigmentation-related absorption variations, four different melanosome volume fraction values in the LE layer were used to simulate very light skin (phototype I) to dark skin (phototype IV) types. Considering our targeted clinical application of interest, it was chosen to firstly consider a melanosome volume fraction $f_{v,mel}^2 = 0.01$, corresponding to phototype I, which is much more prone to the development of skin carcinomas. In Section 3.2.3, darker skins are considered. The correspondence between phototypes and volume fractions of melanosome, adapted from [50], is summarized in Table 3.
<table><caption>**Table 3.** Values of the volume fraction of melanosome in the living epidermis used for simulating skin phototypes I to IV.</caption><thead><tr><th>Skin Type</th><th>Very Fair</th><th>Fair</th><th>Moderately Fair</th><th>Dark Skin</th></tr></thead><tbody><tr><td>Phototype</td><td>I</td><td>II</td><td>III</td><td>IV</td></tr><tr><td>Melanosome volume fraction f<sup>2</sup><sub>v,mel</sub></td><td>0.01</td><td>0.04</td><td>0.08</td><td>0.11</td></tr></tbody></table>

In our skin model, this parameter only influences the LE absorption coefficient, but its impact is significant. The absorption spectra resulting from the different calculations are finally shown in Figure 2. The top right subplot focuses on the absorption coefficient of LE.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Absorption ($\mu_a$ solid line) and scattering ($\mu_s$ dashed-dot line) coefficients of each layer (SC, LE, PD, RD and SF in blue, red, orange, violet and green, respectively) of the designed numerical medium set for skin phototype I and mean dermal blood volume fraction. The spectral curves of the LE absorption coefficient $\mu_a^2$ as a function of skin phototypes I, II, III and IV appear in the top-right subplot, while middle-right and bottom-right subplots focus on the dermal blood volume fraction impact over the PD and RD absorption coefficient $\mu_a^3$ and $\mu_a^4$, as well as their scattering coefficient $\mu_s^3$ and $\mu_s^4$. LE: living epidermis; SC: stratum corneum; PD: papillary dermis; RD: reticular dermis; SF: subcutaneous fat.</figcaption></figure>

Similarly, the impact of the amount of blood in the dermal layers, modified by the volume fraction of blood $f_{v,blood}^3$ and $f_{v,blood}^4$, was considered as an additional parameter of interest [51]. Simulations were then conducted for three pairs of blood volume fractions to represent three media; respectively, weakly, moderately and strongly vascularized. The numerical values of these configurations are shown in Table 4. Here again, the couple of mean values—i.e., $f_{v,blood}^3 = 0.024$ and $f_{v,blood}^4 = 0.02$—was used for all simulations, except in Section 3.2.4, which focuses on this blood amount effect.
<table><caption>**Table 4.** Values of the volume fraction of blood in the papillary dermis (PD) and reticular dermis (RD) used for simulating a dermis with more or less blood content.</caption><thead><tr><th>Blood Content</th><th>Weak</th><th>Mean</th><th>Strong</th></tr></thead><tbody><tr><td>PD blood volume fraction f<sup>3</sup><sub>v,blood</sub></td><td>0.0144</td><td>0.024</td><td>0.0336</td></tr><tr><td>RD blood volume fraction f<sup>4</sup><sub>v,blood</sub></td><td>0.012</td><td>0.02</td><td>0.028</td></tr></tbody></table>

This variation influences mainly the absorption coefficient of the two dermal layers (cf. Equation (1)) but slightly influences the scattering coefficient (cf. Equation (3)). The optical parameters of these three couples are then isolated in the middle-right ($\mu_a^3$ and $\mu_a^4$) and bottom-right ($\mu_s^3$ and $\mu_s^4$) subplot of Figure 2.

The values of the volume fraction and mass concentration of each chromophore applied in every layer of our numerical model are summarized in Table 5. They were chosen following (i) a literature review [47,49,52,53] and (ii) comparisons between simulated spectra and experimental data acquired in clinics [38] (see Figure 3) using our SR-DRS device [37].

<table><caption>**Table 5.** Values of the geometrical and optical parameters of our five layer skin model. SC, LE and SF thickness ranges, $f_{v, mel}^k$ melanosome fraction range in LE and $f_{v, blood}^k$ blood fraction in dermis layers refer to Tables 2–4, respectively.</caption><thead><tr><th>Name (Unit)</th><th>Symbol</th><th>SC</th><th>LE</th><th>PD</th><th>RD</th><th>SF</th></tr></thead><tbody><tr><td colspan="7">Absorption parameters</td></tr><tr><td rowspan="6">Volume fraction (dimensionless)</td><td>f<sup>k</sup><sub>v,water</sub></td><td>0.35</td><td>0.6</td><td>0.5</td><td>0.7</td><td>0.05</td></tr><tr><td>f<sup>k</sup><sub>v,lipids</sub></td><td>0.2</td><td>0.15</td><td>0.173</td><td>0.173</td><td>0.753</td></tr><tr><td>f<sup>k</sup><sub>v,β-caroten</sub></td><td>0.45</td><td>0.24</td><td>0.303</td><td>0.125</td><td>0.141</td></tr><tr><td>f<sup>k</sup><sub>v, mel</sub></td><td>0</td><td>See Table 3</td><td>0</td><td>0</td><td>0</td></tr><tr><td>f<sup>k</sup><sub>v, blood</sub></td><td>0</td><td>0</td><td>See Table 4</td><td>See Table 4</td><td>0.005</td></tr><tr><td rowspan="5">Concentration (g/L)</td><td>c<sup>k</sup><sub>β-caroten</sub></td><td>2.1 × 10<sup>-4</sup></td><td>2.1 × 10<sup>-4</sup></td><td>7.0 × 10<sup>-5</sup></td><td>7.0 × 10<sup>-5</sup></td><td>0</td></tr><tr><td>c<sup>k</sup><sub>eum</sub></td><td>0</td><td>80</td><td>0</td><td>0</td><td>0</td></tr><tr><td>c<sup>k</sup><sub>phe</sub></td><td>0</td><td>12</td><td>0</td><td>0</td><td>0</td></tr><tr><td>c<sup>k</sup><sub>bilirubin</sub></td><td>0</td><td>0</td><td>0.05</td><td>0.05</td><td>0.05</td></tr><tr><td>c<sup>k</sup><sub>HbO<sub>2</sub>/Hb</sub></td><td>0</td><td>0</td><td>150</td><td>150</td><td>150</td></tr><tr><td colspan="7">Scattering parameters</td></tr><tr><td>Volume fraction</td><td>f<sup>k</sup><sub>v,blood</sub></td><td>0</td><td>0</td><td>See Table 4</td><td>See Table 4</td><td>0.05</td></tr><tr><td>Vessel diam. (cm)</td><td>d<sup>k</sup><sub>vessels</sub></td><td>/</td><td>/</td><td>6.0 × 10<sup>-4</sup></td><td>1.5 × 10<sup>-3</sup></td><td>7.5 × 10<sup>-3</sup></td></tr><tr><td>Scat. coeff. (cm<sup>-1</sup>)</td><td>μ<sup>k</sup><sub>s,577</sub></td><td>400</td><td>300</td><td>120</td><td>120</td><td>130</td></tr><tr><td colspan="7">Other optical properties and medium geometry</td></tr><tr><td>Layer thick. (μm)</td><td>z<sup>k</sup></td><td>See Table 2</td><td>See Table 2</td><td>200</td><td>1800</td><td>See Table 2</td></tr><tr><td>Anisotropy factor</td><td>g<sup>k</sup></td><td>0.92</td><td>0.75</td><td>0.72</td><td>0.72</td><td>0.9</td></tr><tr><td>Refractive index</td><td>η<sup>k</sup></td><td>1.55</td><td>1.44</td><td>1.39</td><td>1.38</td><td>1.34</td></tr></tbody></table>

## Scattering Properties

The values of skin scattering coefficients $\mu_s^k(\lambda)$, expressed in cm$^{-1}$ in the various layers of our model, were obtained with Equation (3) given from [53] featuring two contributions: a blood scattering component and a bloodless “base” scattering component, such that

$$ \mu_s^k(\lambda) = (1 - f_{v,blood}^k) \mu_{s,577}^k \frac{577}{\lambda} + f_{v,blood}^k \frac{\mu_{s,blood}(\lambda)}{(1 + a(0.5\mu_{s,blood}(\lambda)d_{vessels}^k))^b} \quad (3) $$
In this first term, $\mu_{s,577}^k$ stands for the scattering coefficient value at 577 nm of the kth layer, and $\lambda$ in the denominator is in nm. For the second term, $\mu_{s,blood}(\lambda)$ is the scattering coefficient of pure blood in cm$^{-1}$ [49]. Moreover, we set $a = 1.007$ and $b = 1.228$ [53], while $d_{vessels}^k$ is the mean vessel diameter of the layer, expressed in cm. The numerical values of the parameters implemented in our simulation are summarized in Table 5, and the corresponding scattering spectra $\mu_s^k(\lambda)$ used in our model are presented for each layer in Figure 2.

## Other Optical Parameters

Both the anisotropy factor $g$ and refractive index $n$ were considered as wavelength-independent in our spectral range of study. For example, concerning near UV wavelengths, this assumption results in a maximum relative error of respectively 3% and 1.5% for the LE and two dermis layers [54]. The values of the former were set to $g^1 = 0.92$, $g^2 = 0.75$, $g^3 = 0.72$, $g^4 = 0.72$ and $g^5 = 0.9$ for each of the five layers [55], while the second values were set to $n^1 = 1.55$, $n^2 = 1.44$, $n^3 = 1.39$, $n^4 = 1.38$ and $n^5 = 1.34$ in each of the five layers, respectively [52].

## 2.3. Simulation Configuration, Method and Tools

### 2.3.1. Simulation Parameters Description

The simulation algorithm developed was based on the *CudaMCML* code from Alerstam et al. [56]. It was modified to consider the probe geometrical configuration mentioned in Section 2.1. Multi-thread calculations were performed on a DELL T630 device using the Intel Xeon E5-2683v4 2.1 GHz, AVX CPU and the GEFORCE GTX 1080 Ti GPU card cluster available within the High Performance Computing resources center “Explor” at the University of Lorraine. Technical aspects of the simulation are summarized in Table 6. The computation time $\Delta T_{sim}$ was different for the media simulated within an interval of 900–1500 s. This duration corresponds to the launch of the $10^8$ photons repeated for each of the 81 wavelengths used to build a DR spectrum between 365 and 765 nm with a spectral step of $\Delta \lambda = 5$ nm. This value was chosen to satisfy the best compromise between calculation time and the spectrum appearance. Ten simulations with $10^8$ photons were performed to quantify the trial by trial fluctuation for the healthy mean epidermal thickness medium of prototype I (see Section 3.1). This led to the conclusion that the relative error in DR spectra compared to the average of the 10 simulations was less than 0.5% in the worst case; i.e., when the least photons are detected, in the hemoglobin absorption peak (415 nm) at the $D_4$ SD distance. The spatial resolutions according to radial distance and depth were respectively set to $dz = 5$ µm and $dr = 10$ µm. Finally, concerning the photon packet stopping criteria, the weight threshold was set to $w_{th} = w_0 \times 10^{-4}$, with $w_0$ the initial photon packet weight, while the time of flight threshold was established at $t_{th} = 1$ s.

<table><caption>**Table 6.** Technical characteristics of the GPU-based parallel computing and MC simulation configurations.</caption><thead><tr><th>Photon Number</th><th>Block Number</th><th>Threads per Block</th><th>ΔT<sub>sim</sub> (s)</th><th>w<sub>th</sub></th><th>t<sub>th</sub> (s)</th></tr></thead><tbody><tr><td>81 × 10<sup>8</sup></td><td>28</td><td>320</td><td>[900–1500]</td><td>w<sub>0</sub> × 10<sup>-4</sup></td><td>1</td></tr></tbody></table>

Details of the photon propagation algorithm and detection conditions are presented in our previous work [46]. The changes introduced to the original code [56] allow for the acquisition of simulated spectra $DR_{D_n}(\lambda)$, the unit of which is the photon ratio, comprised of photons that have reached the DF located at distance $D_n$.

### 2.3.2. Probed Depth Histograms

Probed depth histograms were determined to enable the quantitative analysis of the detected photon penetration depth. The algorithm that allowed us to obtain the latter data set is presented below. It can be applied at a given $D_n$ distance and at a given wavelength $\lambda$.
The normalization line in Algorithm 1 was implemented to eliminate the signal level and thus allow for a proper comparison between histograms of different SD distances. The occurrence of histograms is then relative. A very similar approach was presented by Arifler et al. [31] and applied to resolve spectral information from epithelium and stroma layers in SR-DR spectra. The histogram created here has the particularity of considering the deepest coordinate of the trajectory of the detected photon. This can be associated with the deepest scattering event or to a reflection event between two consecutive layers.

<caption>**Algorithm 1:** Maximum depth histogram of weighted detected photons</caption>

**Result:** creation of the histogram *histoZmax*
**for every photon launched do**
| **while** photon is propagating in the medium **do**
| | store its maximal depth of scattering event $z_{max}$;
| **end**
| **if** photon is detected (i.e., collected by one of the DF) **then**
| | store its detected weight in variable *photonW*;
| | increment the number of *histoZmax* corresponding to $z_{max}$ by *photonW*<fnref n="1" />;
| **else**
| | photon is lost;
| | (definitively absorbed inside medium or exits the medium outside any DF)
| **end**
**end**
normalization of *histoZmax* by dividing by $DR_{D_n}(\lambda)$

### 3. Results and Discussion

Results of DR spectra and probed maximum depth distribution graphs are detailed in Section 3.1 for the healthy skin model with an intermediate thickness, mean dermal blood content and very fair prototype. Then, Section 3.2 focuses on the effect of the epidermis thickness, pathological state, prototype and dermal blood content on photon propagation. Based on these results and their analysis, a quantitative decomposition of the DR signal is then proposed in Section 3.3.

#### 3.1. Detailed Results for Healthy Skin Model with Intermediate Thicknesses

##### 3.1.1. SR-DR Simulated Spectra

Figure 3 (left) shows the DR spectra simulated for a model configuration corresponding to healthy skin of intermediate thickness (cf. Table 2), intermediate blood content (cf. Table 4) and prototype I (cf. Table 3) at each SD distance. DR intensity tends to increase toward the red spectral range due to the lower opacity of the near-infrared (near IR) wavelengths. Local minimums at 415, 540 and 575 nm correspond to hemoglobin absorption.

Figure 3 (right) shows the superposition of the normalized simulated spectra and the clinical spectra for each SD distance. Each of the subplots is composed of the the average clinical spectrum of three acquisitions $\overline{DR}^{clin}$ (dashed line) of this same quantity, to which we added the standard deviation (dash-dotted line) and finally the simulated spectrum (continuous line), for which model parameters (epidermal thickness, prototype, dermal blood content) were determined by inverse problem solving to fit clinical data. Those experimental DR spectra was acquired on the anterior skin surface of a patient's wrist. This skin site was clinically considered as healthy by the plastic surgeon in charge of the patient's inclusion in the *SpectroLive* clinical trial (ClinicalTrials.gov Identifier: NCT02956265) approved by the Institutional Board for Ethics and Patients Safety (CPP Est III authorization: 16.06.03) and by the French National Agency for Medicines and Health Products Safety (ANSM authorization: DMDPT-RIAL/MM/2016-A00608-43). Comparing spectra from

<footnote n="1">¹ Here, the maximal depth $z_{max}$ is discretized to allow the design of the histogram. The width of the length interval is the spatial resolution in depth i.e., $dz = 5$ µm.</footnote>
a MC simulation with clinical spectra is not an easy task [57]; thus, both sets of spectra were normalized (i.e., we set the minimum value to 0 and the maximum value to 1 in the 400–650 nm exploitable spectral band). The fitting between acquired and simulated spectra partly asserts the validity of the numerical skin model and more generally the fidelity of the simulation. A more complex skin model including skin intrinsic fluorophore absorption and emission would be interesting to increase the level of media accuracy, even if, for the considered spectral band, their influences seems to be negligible in comparison with other considered chromophores.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** DR spectra simulated for the four SD distances (**left**). The considered medium is the healthy skin model with intermediate thicknesses, mean dermal blood content and prototype I. The four other figures (**right**) show the superposition of the normalized $\overline{DR}^{\text{clin}}$ (dashed line) clinical spectra from *SpectroLive* acquisition on a patient's wrist, within a confidence interval calculated as $\overline{DR}^{\text{clin}} \pm \text{std}$ (dash-dotted line), as well as the estimated $DR^{\text{sim}}$ (continuous line) from our simulation.</figcaption></figure>

### 3.1.2. Depth Histogram and Exploitation

#### Depth Histogram

Quantitative information is presented in Figure 4 to describe and analyze the maximum penetration depth of the detected photons. The normalized y-axis of the histograms is the percentage of occurrence of detected photons for a given SD distance. The presence of abrupt transitions is noticeable in the histograms corresponding to a refractive index mismatch (i.e., reflection) between consecutive layers, represented by a black vertical line. Finally, a box-plot approach was chosen to extract a characteristic probed depth associated with the couple ($\lambda, D_n$). Thus, the three vertical colored lines for each of the histograms represent, from left to right, the first quartile (dotted) $z_{q1}$, the median (solid line) $z_{\text{med}}$ and the third quartile (dotted) $z_{q3}$. The colored area under the histogram, between the two quartiles, then represents 50% of the total area. In other words, half of the photons contributing to the $DR_{D_n}(\lambda)$ value have reached at most a depth between $z_{q1}$ and $z_{q3}$. For the considered wavelength—i.e., $\lambda = 565$ nm—a large part of the detected photons traveled through LE and RD for the shortest and the longest distance, respectively. For the second and third distance, the maximum penetration depths of detected photons are distributed between those layers; i.e., in the PD layer. A shorter distance would be of interest for this wavelength to focus on superficial SC and LE layers.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Normalized probed depth histograms of weighted detected DR photons at $\lambda = 565$ nm (central wavelength of the spectrum 365–765 nm) for the four SD distances. Vertical black lines represent the borders between layers. Colored vertical lines represent, from left to right, the first quartile (dashed-dot), the median (solid) and the third quartile (dashed-dot). The considered medium is the healthy skin model with intermediate thicknesses, prototype I and mean blood dermal content.</figcaption></figure>

Extraction of a Typical Probed Depth and Spectral Behavior

The histograms presented in Figure 4 are available at all wavelengths and for each of the four SD distances, from which the median and quartiles of the maximal penetration depth can be calculated. The choice to synthesize our depth histograms using a box-plot approach rather than the mean penetration depth [31] of these histograms is justified by its simplicity of interpretation from the point of view of the resulting DR signal. Indeed, by considering the signal captured for the pair ($\lambda, D_n$), we can decompose it into four equal contributions of different depths, for which the quartiles and medians represent the boundaries. Moreover, the interquartile range gives a good idea of the depth spread of the photons around the median, thus illustrating the enlargement of the probed volume in the three-dimensional medium.

Other alternatives were developed to extract a penetration depth for collected photons in SR-DR spectroscopy. In their study [32], Tseng et al. introduce a quantity called “interrogation depth”, defined as the z average of all detected photon scattering events in a two-layer skin model. The photon penetration was not the main focus of their study, but the paper revealed that their SD distances in a range of 1440–2880 µm allow researchers to probe between 293–949 µm within the 550–860 nm wavelength region. Other teams have adopted an original and experimentally transposable technique to extract a sampling depth $Z_S$ in a mono-layer medium defined by its thickness $Z_0$ and its optical coefficients $\mu_a$ and $\mu'_s = \mu_s(1 - g)$. This approach consists of varying $Z_0$ from 0 to a sufficiently large thickness so that it no longer has any impact on the DR signal. For a very thin thickness, photons escape through the bottom of the layer, which is responsible for a very weak signal. Thickening the medium allows more photons to reach the DF, until it becomes wide enough to converge to a maximum (semi-infinite medium). $Z_S$ is then defined as the thickness at which the DR signal reaches 50% of its final maximum value [58,59]. Unfortunately, this method can not be employed for multi-layered tissue. Finally, the depth measurement can
be defined as the median of the cumulative absorption of fluence absorption maps. This was applied in the study by Arimoto et al. [30] for a three-layer skin medium considering the water absorption and scattering in the near IR spectral range.

The results obtained from our box-plot approach from maximum depth histograms are presented in Figure 5, with the same graphic conventions (cf. Figure 4). The quasi-linearity of the probed depth can be seen as a function of the two varying parameters—i.e., the SD distance and the wavelength—except at the first hemoglobin absorption peak (415 nm). More importantly, it can be noticed that by choosing a couple ($\lambda, D_n$), it is possible to probe depths in a range of values covering the interval 50–900 µm. Moreover, one can read in that work that the photons that contributed to the spectrum $DR_{D_4}(\lambda)$ and verified $\lambda > 450$ nm mostly visited the RD layer. On the other hand, it is visible that the photons of $DR_{D_1}(\lambda)$ satisfying $\lambda < 450$ nm remained globally confined in SC, LE and PD layers. One can also observe that the inter-quartile range widens when $\lambda$ increases for all distances $D_n$. Indeed, the probed volume is very localized for the near UV wavelengths and spreads when it approaches the near IR.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Maximal probed depth values of detected photons in the healthy skin model with intermediate thicknesses of SC and LE layers and mean dermal blood content, as a function of wavelength and SD distances. Solid and dashed-dot lines represent median and quartiles values, respectively. The right column and the black horizontal lines illustrate the position of the border between the layers (SC: stratum corneum; LE: living epidermis; PD: papillary dermis; RD: reticular dermis; SF: subcutaneous fat (not represented here)).</figcaption></figure>

## 3.2. Analysis of Spectral and Probed Depth Data Related to Medium Configuration

### 3.2.1. Epidermis Thickness Modifications in Healthy Skin Model

Some studies have focused on the impact of the thickening of the most superficial skin layer on DR spectra acquisition in a skin multi-layered medium [60]. Among them, six levels of epidermis thickness in the range of 0–500 µm have been tested in the three-layer model (epidermis and two dermal layers) [36], leading to the conclusion that its impact on the simulated spectra is very significant. In [61], the authors quantified the accuracy of estimating the optical properties of the upper and lower layers in a bi-layer medium when the thickness of the most superficial layer varies. Thus, they showed that the accuracy of estimation of the upper (respectively lower) layer increases (decreases) when the thickness of the surface layer increases over the range 0–1000 µm. However, those studies do not show an evolution of the probed depth with the thickening of the superficial layer(s). The skin model configurations considered in this section are those described in Table 2. For the sake of conciseness and clarity, only the DR spectra and the graphs of probed depths are
presented (see Figure 6), where only the median of the maximum depth value is shown for each of the media (quartiles are not shown).

<figure><img src="image_9.png" /><figcaption>**Figure 6.** DR spectra (**left**) and maximum probed depth (**right**) for the thin (dotted line), mean (dashed) and thick epidermal thickness media (dashed-dot). The SD distance $D_n$ is indicated by the color code. The three right columns represent the position of the border between the layers of each skin model (thin, mean, thick).</figcaption></figure>

By looking at the DR spectra, it is noteworthy to state that (i) at $D_1$, the amplitude of the DR spectrum for the maximum thickness of the epidermis is (for 365–450 nm and 650–765 nm spectral ranges) greater than that for the minimum thickness of the epidermis, while (ii) along with longer SD distances, the amplitude of the DR spectra for the maximum thickness of the epidermis tends to be lower than the one for the minimum thickness of the epidermis. If we observe the probed depth at $D_1$ for the thickest epidermis model, we notice that most of the photons verify $\lambda < 550$ nm, wgucg contributes to $DR_{D_1}(\lambda)$ remaining confined in SC and LE layers. For the thinnest epidermis model, those photons largely visited the PD layer. The $DR_{D_1}(\lambda)$ spectra evolution for the different thicknesses may be related to the high scattering capacity of the SC and LE layers at the origin of a short path length return of many photons to the surface (and thus captured at $D_1$). Indeed, the figure of the probed depths shows the well-known fact [62] that, in general, the increase of the SD distance forces photons to dive deeper in the medium: the superficial back-scattering effect is no longer observed. Finally, we can notice that the wide spectral range of the excitation light and the geometry of the device (SD distances available) allow us to probe different depths without redundancy, and this is true for every epidermal thickness.

As indicated in Section 2.2, both dermal layers (PD and RD) were considered of fixed thickness. This choice results from a preliminary study in which the impact of the dermal thickness was considered secondary in comparison with the other parameters that we vary in this paper. As we have treated the impact of the epidermal thickness, we obtained results for three pairs ($z^3 = 150$; $z^4 = 1600$), ($z^3 = 200$; $z^4 = 1800$), ($z^3 = 250$; $z^4 = 2000$) µm of PD and RD thicknesses, respectively, representing minimum, mean and maximum dermal thicknesses. The DR and maximal probed depth appear in Figure 7. It can be noticed that the impact of this parameter on the signal (**left**) and on the detected photon penetration (**right**) is very small (variations smaller than 2% in DR for $D_1$ and for $\lambda > 600$ nm).
<figure><img src="image_10.png" /><figcaption>**Figure 7.** DR spectra (**left**) and maximum probed depth (**right**) for the thin (dotted line), mean (dashed) and thick dermal thickness media (dashed-dot). The SD distance $D_n$ is indicated by the color code. The three right columns represent the position of the border between the layers of each skin model (thin, mean, thick).</figcaption></figure>

### 3.2.2. Epidermis Modifications in Healthy and Pathological Skin Models

We present in Figure 8 the DR spectra obtained for skin models simulating two types of cutaneous carcinoma: ulcerated and hyperkeratosed (cf. Table 2). By way of comparison, the healthy medium of average thickness also appears in the various results.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** DR spectra (**left**) and maximum probed depth (**right**) for the simulated skin models: healthy (dotted line), hyperkeratosed (dashed) and ulcerated (dashed-dot). The SD distance $D_n$ is indicated by the color code. The three right columns represent the position of the border between the layers of each skin model (healthy, hyperkeratosed, ulcerated).</figcaption></figure>

DR spectra for healthy and ulcerated skin models are very similar; i.e., they are identical for $D_2$, $D_3$ and $D_4$. Slight distinct amplitudes are observed at $D_1$ but with a globally identical spectral shape. Indeed, the disappearance of the SC layer, which differs from both, can be seen as an overall decrease in the thickness of the epidermis $z^{Epi} = z^1 + z^2$. The back-scattering effect for $D_1$ is thus more important when $z^{Epi}$ is greater (healthy medium) than when it decreases (ulcerated medium). The associated signal increases and then fades when the SD distance increases.

Unlike healthy as well as ulcerated skin models, the results obtained for the hyperkeratosis model show that the majority of photons detected, whatever the SD distance, remained confined in SC (for $D_1$ and $D_2$) and LE (for $D_3$ and $D_4$). This information is confirmed by the shape of the DR spectra, since the dermal hemoglobin signature (pits at 415, 540 and 575 nm) appears only slightly. One can again notice the superficial back-scattering effect responsible for much larger amplitudes for every $D_n$, except for $D_1$ for the near IR
wavelengths. Such results emphasize the importance of taking into account the skin layer thickness of each pathological case to properly diagnose it.

### 3.2.3. Phototype Modifications in Healthy Skin Model

Figure 9 shows the DR spectra and probed depth graphs obtained for the healthy skin medium of intermediate thickness (cf. Table 2), mean dermal blood content (cf. Table 4) and for the four skin phototypes considered (cf. corresponding values of volume fraction of melanosome $f_{v,me1}^2$ in Table 3 and LE absorption coefficient $\mu_a^2$ in Figure 2, top right insert).

<figure><img src="image_12.png" /><figcaption>**Figure 9.** DR spectra (**left**) and maximum probed depth (**right**) for the four phototypes: I (line), II (dashed), III (dashed-dot) and IV (dotted). The SD distance $D_n$ is indicated by the color code. The right column represents the position of the border between the layers for the common geometrical configuration.</figcaption></figure>

The decay of the signal in the near UV and the beginning of the visible light spectrum, when the skin darkens, was expected and is the most obvious characteristic. Rather than simply diffusing, the LE becomes more absorbent with increasing pigmentation. The optical barrier of the epidermis (including SC+LE) becomes increasingly difficult to penetrate for photons verifying $\lambda < 450$ nm. For $D_1$ (shortest paths), the maximum depth probed by the detected photons decreases when $f_{v,me1}^2$ increases. The information carried by $DR_{D_1}$ ($\lambda < 450$) appears thus to be very superficial with $z_{max}$ confined in LE layer. The very weak signature of the hemoglobin peak in the spectra for the dark phototypes confirms the previous assertion, since photons that reach the vascularized layers are rare. However, for this same distance and for the rest of the spectrum ($\lambda > 500$ nm), it is interesting to note that the probed depth, in the PD, is virtually unchanged with the phenotype. The path (i.e., the trajectory) traveled by these photons seems then to be globally identical, and only the absorption during the forward passage—the way back in the LE—justifies the signal drop in the DR spectrum. This is partly explained by the fact that the change in phenotype in our model is simply considered by an increase in the melanosome volume fraction. A more complex model taking into account the volume distribution and shape of melanosomes varying with phenotype would probably have more impact on the reduced scattering coefficient and thus on the photon trajectories in the medium. This specific behavior observed at short distances ($D_1$), which is different for short ($\lambda < 450$ nm) and long ($\lambda > 500$ nm) wavelengths, might be of particular interest, namely to improve the robust and precise estimation of the absorption coefficient $\mu_a^2(\lambda)$ (and thus by extension a value for the volume fraction of melanosome) of the LE layer, specifically when solving the complex inverse problem of multiple optical properties in multilayer skin models.

For the other, longer SD distances, one can notice that all the detected photons reach at least the RD layer for $D_2, D_3, D_4$ and $\lambda > 500$ nm, and it is possible to observe the sensitivity of $z_{max}$ to the hemoglobin absorption peak (at 415 nm), especially for $D_2$ and $D_3$. The increasing probed depth with phenotype, for all longer SD distances $D_n$ ($n > 1$), can be explained by considering two complementary behaviors during the photon penetrating
trajectory from SC to RD layers. The first is to dive as little as possible into the medium and thus reduce the optical path traveled. The second is to quickly reach the PD in order to escape the LE with its high absorbing capacity. For fair phototypes, where the LE absorption is moderate, the first behavior seems to have more impact; that is, detected photons resulted from shallower paths. On the other hand, for dark phototypes, it is more probable to detect photons that cross the LE layer with the shortest path lengths; i.e., paths that are as vertical or straight as possible. These photons thus need more time to make their way back towards the detection fiber: the probed depth is then more important.

### 3.2.4. Blood Content in Dermal Layer Modifications in Healthy Skin Model

Figure 10 shows the DR spectra and probed depth graphs obtained for the healthy skin medium of intermediate thickness (cf. Table 2) for the fair phenotype I (cf. Table 3) and for the three media representing dermal layers with more or less blood content (cf. Table 4), for which absorption and scattering coefficients of the two concerned layers (PD and RD), respectively, appear on the middle-right and bottom-right subplots of Figure 2. The most remarkable and expected observation is the relative signal decrease in the three hemoglobin absorption peaks (415, 540 and 575 nm). This can be linked to the probed depth plot where, for distances $D_2$ to $D_4$, the spectral region where hemoglobin is absorbed (i.e., 365–630 nm) is affected by a decrease in the detected photon penetration with increasing blood concentration.

Concerning $D_1$, the interpretation is even simpler, because of the quasi-non-influence of the quantity of blood on the path of the detected photons. In $DR_{D_1}(\lambda)$, we notice that for the first part of the spectrum (365–600 nm), the blood concentration is responsible for an attenuation of the signal, while it is associated with a gain of intensity for the rest of the spectrum (600–765 nm). The analysis of the optical properties on these spectral bands leads us to the interpretation that the first phenomenon is caused by the absorption of hemoglobin in the dermal layers. For the second phenomenon, the absorption of blood becomes negligible, and only the scattering coefficients remain significantly different (about 15% between the two extremes maximum and minimum $f_{v,blood}$). This means then that the scattering power caused by the increase of the blood concentration of the dermal layers (particular PD for this shorter distance $D_1$) is at the origin of many short optical paths between SF and DF and therefore responsible for this non-obvious signal increase.

<figure><img src="image_13.png" /><figcaption>**Figure 10.** DR spectra (**left**) and maximum probed depth (**right**) for the three medium configurations of minimum (line), mean (dotted) and maximal (dashed) dermal blood contents. The SD distance $D_n$ is indicated by the color code. The right column represent the position of the border between the layers for the common geometrical configuration.</figcaption></figure>

### 3.3. A Proposal of a Model for DR Signal Decomposition by Skin Layer

In order to more quantitatively determine the content of the global DR signal in terms of layer contributions, we propose to exploit the results of analysis of the aforementioned developed maximum probed depth histograms and spectra to decompose the total signal
$DR_{D_n}(\lambda)$, as a sum of the DR contributions related to each of the five skin layers considered, such that for each SD distance $D_n$,

$$DR_{D_n}(\lambda) = DR_{D_n}^{SC}(\lambda) + DR_{D_n}^{LE}(\lambda) + DR_{D_n}^{PD}(\lambda) + DR_{D_n}^{RD}(\lambda) + DR_{D_n}^{SF}(\lambda) \quad (4)$$

where $DR_{D_n}^{layer}(\lambda)$ is composed of photons that reached at deepest the layer appearing in exponent verifying $d_{min}^{layer} < z_{max} < d_{max}^{layer}$ with $d_{min}^{layer}$ and $d_{max}^{layer}$, respectively, the minimum and maximum depth coordinates of the layer. In practice, this decomposition was obtained by integrating the histograms of probed depths (cf. Section 3.1.2 and Figure 4) between the upper and lower border of each skin layer. Figure 11 illustrates this decomposition for the four SD distances and the healthy skin model with intermediate thicknesses, mean dermal blood content for SC and LE layers and prototype I. This “deepest layer probed” separation criterion has already been proposed in the literature. In [36], Fredriksson and his team studied the path length distributions for detected photons in a three-layered skin model, although they did not propose any decomposition approach. However, a different “layer”-based approach of DR signal decomposition has been proposed in a few studies. In [31], the contribution of the detected photon associated with the k-layer was defined as the quotient of the number of scattering events undergone in this layer on the total number of scattering events between SF and DF. In [35], the authors proposed to explain the skin layer sensitivity by the mean partial path length, which consists of assigning to the layer k the relative weight $l_k/l_{tot}$ in the DR signal, with $l_k$ the path length in the layer k and $l_{tot}$ the total path length between SF and DF.

<figure><img src="image_14.png" /><figcaption>Figure 11. DR signal decomposition according to Equation (4) for the distances $D_1-D_4$ in the healthy skin model with intermediate thicknesses of SC and LE layers, mean dermal blood content and prototype I. The correspondence between symbols and layers considered is indicated in the legend on the right.</figcaption></figure>

The approach proposed in the present contribution differs from the latter by the following points. First, our proposed decomposition provides a straightforward and simple interpretation regarding the layer reached by the detected photons making the global DR signal, compared to less non-trivial interpretations linking the mean path length or the ratio of scattering events to the layer's weight in the DR signal. Second, as a consequence, the proposed decomposition should help in designing a method to estimate sequential optical properties by studying the layers from the shallowest to the deepest and by providing new,
a priori knowledge on the selection of the sensitive distances and wavelength range. This approach will bring a physical and anatomical sense to the estimation and should be less time consuming due to the successive isolation of the parameters.

The analysis to be drawn from Figure 11 is in accordance with that performed qualitatively in the Section 3.1 but allows us to better explain the DR signal in layer decomposition. To further compare the proportion of a signal from a layer, for a fixed couple (λ, Dn), we defined the weight fraction of the layer contribution in percentage, as follows:

$$w_{D_n}^{\text{layer}}(\lambda) = 100 \times \frac{DR_{D_n}^{\text{layer}}(\lambda)}{DR_{D_n}(\lambda)} \quad (\%) \qquad (5)$$

Finally, Figure 12 shows the associated evolution of $w_{D_n}^{\text{layer}}(\lambda)$ defined in Equation (5) as a function of the wavelength and the SD distance. This decomposition turns out to be very useful to extract the spectral band and SD distance of a DR signal that is susceptible to carry information associated with local depth information. For instance, the spectral “signature” of elastosis (accumulation of disordered elastic tissue in the PD due to carcinogenesis) should be looked for in the UV range ($365 < \lambda < 400$ nm) at the $D_2$ distance, where more than 40% of the DR signal is constituted of photons that have visited the concerned layer.

<figure><img src="image_15.png" /><figcaption>**Figure 12.** Weight of each layer (from SC at the top to SF at the bottom) in DR signal according to Equation (5) for distances $D_1-D_4$ in the healthy skin model with intermediate thicknesses of SC and LE layers, mean dermal blood content and prototype I. The correspondence between symbols and layers is still the same as indicated in Figure 11.</figcaption></figure>

It is worth noting that all the geometrical and optical parameters considered in the model described above were chosen within standard ranges of values, accounting for anatomical and physiological variability. The aim of the present work was not to perform a complete variation study taking into account those variability intervals but rather to investigate the variation of the dependence of epidermal thickness (see Section 3.2.1), pathological state (see Section 3.2.2), prototype (see Section 3.2.3) and dermal blood content (see Section 3.2.4) considering a unique set of basic but representative values for all other parameters of the model. A possible extended study could be led with other parameters, such as the individually variable scattering coefficient [63], the blood oxygen saturation [60] or the average blood vessel diameter [64], even if those parameter have a more limited impact on the layers' optical coefficients and thus on the penetration of photons in depth. It would also be interesting to consider 3D tumor models, since the assumption of semi-
infinite layers imposed by *CudaMCML* is a limitation. A study shows that the validity of this assumption depends on the type of tumor simulated (squamous cell carcinoma in epithelium or a basal cell carcinoma in upper dermis), the angular incidence of the exciting fiber, the height of the tumor and its width relative to the SD distance at which the signal is acquired [65].

All the results established so far allow us to extend our discussion towards the targeted clinical application of SR-DRS devices. Firstly, the in-depth elements characterizing the path of the detected photons within the various skin layers allow us to assert that the geometrical layout of the probe and its spectral configuration provide access to the different layers of the skin with a sensitivity depending on the chosen couple ($\lambda, D_n$) and on the physiological and morphological skin configuration. In particular, the results showed that the photons contributing to the DR spectra were likely to carry the discriminating local information that distinguishes healthy and pathological tissues. It should be noted, however, that the results presented here are valid for the case of fiber-based contact-type DRS configurations, and under the assumption that the optical properties of each of the skin layers do not significantly vary laterally within the ~2 mm diameter circular region covered by the multi-fiber optics probe tip. Moreover, for the estimation of diagnostically relevant optical parameters from in vivo SR-DRS measurements through an inverse problem solving method [66,67], such results should be useful to understand this iterative process better by defining a cost function adapted to the layer structure of the skin.

4. Conclusions and Prospects

In order to learn more about the depth behavior of detected photons and thus allow for a better interpretation of the simulated and measured diffuse reflectance signals, we decided to carry out the current numerical study. A photon propagation algorithm was adapted to match the measurement features of the *SpectroLive* device, which is a spatially-resolved diffuse reflectance spectroscopy device used in clinical applications. To do so, a numerical five-layered skin model, including realistic absorption and scattering properties from main intrinsic chromophores, was defined in order to reach a certain fidelity between the simulation and the experimental acquisition. The histograms of probed depths allowed us to extract the typical depths at which the photons took the return path to the collection fiber. The results dependent on the skin type placed under examination and the fate of photon behavior in terms of depth penetration for different physiological and pathological states were investigated. Indeed, diffuse reflectance spectra and probed depth graphs were analyzed for different sites on the body (by varying the thicknesses of the stratum corneum and the living epidermis), for different pathological states (ulcerated and hyperkeratosed carcinoma) and finally for different phototypes (by varying the volume fraction of melanosome in living epidermis).

Finally, to further the understanding of the diffuse reflectance signal, a layer decomposition model of the latter was proposed and applied to a healthy medium of intermediate thickness, mean dermal blood content and prototype I (most carcinomas occur in the fair-skinned population). By providing quantitative layer contributions in the resulting signal, these results refine the existing relationships between the wavelength, source-detector distance and probed skin layer. This relationship will be further used to interpret spectra recorded during a clinical trial that aims at evaluating the SR-DRS diagnostic accuracy for the differential diagnosis of skin carcinomas. More specifically, the exposed numerical evidence may facilitate the development of analytical photon propagation models that enable the inverse estimation of diagnostically relevant optical parameters or allow researchers to distinguish a healthy from a pathological tissue using spectral criteria, from in vivo clinical spectra. Such a numerical study may also be used to standardize measured spectra taking into account intra- and inter-individual variability. Currently, ratios of measured spectra are proposed to eliminate this variability, but such calculation methods can sometimes be impaired by artifacts of measurements.
**Author Contributions:** Conceptualization, V.C.; Data curation, V.C., G.K. and M.A.; Formal analysis, V.C., W.B. and M.A.; Funding acquisition, W.B. and M.A.; Investigation, V.C. and W.B.; Methodology, V.C. and W.B.; Project administration, W.B.; Resources, G.K.; Supervision, W.B., C.D. and M.A.; Validation, V.C.; Writing—original draft, V.C.; Writing—review and editing, V.C., W.B., C.D. and M.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** The authors would like to thank the French Région Grand-Est, the CNRS and the Metz-Thionville Regional Hospital (2016 clinical research award) in the frame of the SpectroLive project, the Contrat de Plan Etat—Région Grand Est 2015–2020 (CPER IT2MP, plateforme IMTI) including the European Regional Development Fund FEDER and the Ligue Contre le Cancer for their financial support; and the Lorraine Université d’Excellence (LUE) Doctorate 2019 under the PhD grant R01PJYRX-PHD19-COLA-CRAN (Victor Colas, 2019–2022).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Pure chromophore spectra mentioned in Section 2.2.2 were downloaded from OMLC database [49]. The basic version (CudaMCML) of the photon transport simulation that we have adapted is available on the Lund University website.

**Acknowledgments:** The authors want to express their gratitude for the high-performance computing resources provided by the EXPLOR centre hosted by the Université de Lorraine (project 2018AM2IX0877).

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Bydlon, T.M.; Nachabé, R.; Ramanujam, N.; Sterenborg, H.J.C.M.; Hendriks, B.H.W. Chromophore based analyses of steady-state diffuse reflectance spectroscopy: Current status and perspectives for clinical adoption. *J. Biophotonics* **2015**, *8*, 9–24. [CrossRef] [PubMed]

2. Bashkatov, A.N.; Genina, E.A.; Tuchin, V.V. Optical properties of skin, subcutaneous, and muscle tissues: A review. *J. Innov. Opt. Health Sci.* **2011**, *04*, 9–38. [CrossRef]

3. Kollias, N.; Seo, I.; Bargo, P.R. Interpreting diffuse reflectance for in vivo skin reactions in terms of chromophores. *J. Biophotonics* **2010**, *3*, 15–24. [CrossRef]

4. Wilson, R.H.; Nadeau, K.P.; Jaworski, F.B.; Tromberg, B.J.; Durkin, A.J. Review of short-wave infrared spectroscopy and imaging methods for biological tissue characterization. *J. Biomed. Opt.* **2015**, *20*, 030901. [CrossRef]

5. Bruulsema, J.T.; Essenpreis, M.; Heinemann, L.; Hayward, J.E.; Berger, M.; Gries, F.A.; Koschinsky, T.; Sandahl-Christiansen, J.; Orskov, H.; Farrell, T.J.; et al. Detection of Changes in Blood Glucose Concentration in vivo with Spatially Resolved Diffuse Reflectance. In *Biomedical Optical Spectroscopy and Diagnostics*; Optical Society of America: Orlando, FL, USA, 1996.

6. Hjalmarsson, P.; Thennadil, S.N. Determination of glucose concentration in tissue-like material using spatially resolved steady-state diffuse reflectance spectroscopy. In *Complex Dynamics and Fluctuations in Biomedical Photonics V*; Tuchin, V.V., Wang, L.V., Eds.; International Society for Optics and Photonics: San Jose, CA, USA, 2008; Volume 6855, pp. 48–56.

7. Finlay, J.; Foster, T. Hemoglobin oxygen saturations in phantoms and in vivo from measurements of steady-state diffuse reflectance at a single, short source-detector separation. *Med. Phys.* **2004**, *31*, 1949–1959. [CrossRef] [PubMed]

8. Spliethoff, J.W.; Evers, D.J.; Klomp, H.M.; van Sandick, J.W.; Wouters, M.W.; Nachabe, R.; Lucassen, G.W.; Hendriks, B.H.; Wesseling, J.; Ruers, T.J. Improved identification of peripheral lung tumors by using diffuse reflectance and fluorescence spectroscopy. *Lung Cancer* **2013**, *80*, 165–171. [CrossRef] [PubMed]

9. Arifler, D.; MacAulay, C.E.; Follen, M.; Richards-Kortum, R.R. Spatially resolved reflectance spectroscopy for diagnosis of cervical precancer: Monte Carlo modeling and comparison to clinical measurements. *J. Biomed. Opt.* **2006**, *11*, 064027. [CrossRef]

10. LaRiviere, B.; Garman, K.S.; Ferguson, N.L.; Fisher, D.A.; Jokerst, N.M. Spatially Resolved Diffuse Reflectance Spectroscopy of Tumor and Normal Excised Human Colon with Thin Film Si Sensors. In Proceedings of the Optical Sensors and Sensing Congress (ES, FTS, HISE, Sensors), San Jose, CA, USA, 25–27 June 2019; Optical Society of America: San Jose, CA, USA, 2019; p. STu2D.5. [CrossRef]

11. Pery, E.; Amouroux, M.; Diaz-Ayil, G.; Blondel, W.C.P.M.; Bourg-Heckly, G.; D'Hallewin, M.; Leroux, A.; Guillemin, F. Multimodality point spectroscopy: Interests, development and application to in vivo diagnosis of bladder and skin cancers. In Proceedings of the 2008 16th Mediterranean Conference on Control and Automation, Ajaccio, France, 25–27 June 2008; pp. 1353–1358. [CrossRef]

12. Nie, Z.; Du Le, V.N.; Cappon, D.; Provias, J.; Murty, N.; Hayward, J.E.; Farrell, T.J.; Patterson, M.S.; McMillan, W.; Fang, Q. Integrated Time-Resolved Fluorescence and Diffuse Reflectance Spectroscopy Instrument for Intraoperative Detection of Brain Tumor Margin. *IEEE J. Sel. Top. Quantum Electron.* **2016**, *22*, 49–57. [CrossRef]
13. Evers, D.J.; Nachabe, R.; Peeters, M.J.V.; van der Hage, J.A.; Oldenburg, H.S.; Rutgers, E.J.; Lucassen, G.W.; Hendriks, B.H.W.; Wesseling, J.; Ruers, T.J.M. Diffuse reflectance spectroscopy: towards clinical application in breast cancer. *Breast Cancer Res. Treat.* **2013**, *137*, 155–165. [CrossRef]

14. Georgakoudi, I.; Jacobson, B.C.; Van Dam, J.; Backman, V.; Wallace, M.B.; Müller, M.G.; Zhang, Q.; Badizadegan, K.; Sun, D.; Thomas, G.A.; et al. Fluorescence, reflectance, and light-scattering spectroscopy for evaluating dysplasia in patients with Barrett's esophagus. *Gastroenterology* **2001**, *120*, 1620–1629. [CrossRef]

15. Garcia-Uribe, A.; Chang, C.C.; Zou, J.; Banerjee, B.; Kuczynski, J.; Wang, L.V. Pancreatic tumor margin detection by oblique incidence diffuse reflectance spectroscopy. In *Proceedings of the Optical Biopsy IX*, San Francisco, CA, USA, 16 February 2011; Alfano, R.R., Ed.; International Society for Optics and Photonics: San Francisco, CA, USA, 2011; Volume 7895, pp. 110–114.

16. Garcia-Uribe, A.; Wang, L.V.; Zou, J.; Smith, E.B.; Duvic, M.D.; Prieto, V. In-vivo characterization of optical properties of pigmented skin lesions including melanoma using oblique incidence diffuse reflectance spectrometry. *J. Biomed. Opt.* **2011**, *16*, 020501. [CrossRef]

17. Pires, L.; Demidov, V.; Vitkin, I.A.; Bagnato, V.S.; Kurachi, C.; Wilson, B.C. Optical clearing of melanoma in vivo: Characterization by diffuse reflectance spectroscopy and optical coherence tomography. *J. Biomed. Opt.* **2016**, *21*, 081210. [CrossRef] [PubMed]

18. Amouroux, M.; Blondel, W.; Granel-Brocard, F.; Marchal, F.; Guillemin, F. A preliminary study on skin phantoms to test spatially resolved-Diffuse Reflectance Spectroscopy as a tool to help diagnose cutaneous melanoma: A non-invasive measurement of Breslow index. *Bio-Med. Mater. Eng.* **2008**, *18*, 339–343. [CrossRef]

19. Khairallah, G.; Blondel, W.; Rakotomanga, P.; Soussen, C.; Delconte, A.; Plenat, F.; Marchal, F.; Amouroux, M. Impact of standardization of autofluorescence and diffuse reflectance spectra on diagnosis accuracy of optical spectroscopy used for skin carcinomas diagnosis. In *Proceedings of the Saratov Fall Meeting SFM'18 International Symposium"Optics and Biophotonics-V*, Saratov, Russia, 24–29 September 2018.

20. Amouroux, M.; Díaz-Ayil, G.; Blondel, W.C.P.M. Classification of ultraviolet irradiated mouse skin histological stages by bimodal spectroscopy: Multiple excitation autofluorescence and diffuse reflectance. *J. Biomed. Opt.* **2009**, *14*, 14011. [CrossRef] [PubMed]

21. Koenig, A.; Roig, B.; Digabel, J.L.; Josse, G.; Dinten, J.M. Accessing deep optical properties of skin using diffuse reflectance spectroscopy. In *Clinical and Biomedical Spectroscopy and Imaging IV*; Brown, J.Q., Deckert, V., Eds.; International Society for Optics and Photonics: Munich, Germany, 2015; Volume 9537, pp. 42–49.

22. Naglič, P.; Pernuš, F.; Likar, B.; Bürmen, M. Estimation of optical properties by spatially resolved reflectance spectroscopy in the subdiffusive regime. *J. Biomed. Opt.* **2016**, *21*, 095003. [CrossRef]

23. Sharma, M.; Marple, E.; Reichenberg, J.; Tunnell, J.W. Design and characterization of a novel multimodal fiber-optic probe and spectroscopy system for skin cancer applications. *Rev. Sci. Instruments* **2014**, *85*, 083101. [CrossRef] [PubMed]

24. Koenig, A.; Grande, S.; Dahel, K.; Planat-Chrétien, A.; Poher, V.; Goujon, C.; Dinten, J.M. Diffuse reflectance spectroscopy: A clinical study of tuberculin skin tests reading. In *Biomedical Applications of Light Scattering VII*; Wax, A.P., Backman, V., Eds.; International Society for Optics and Photonics: San Francisco, CA, USA, 2013; Volume 8592, pp. 98–105. [CrossRef]

25. Breslin, T.M.; Ramanujam, N. Autofluorescence and Diffuse Reflectance Properties Malignant and Benign Breast Tissues. *Ann. Surg. Oncol.* **2003**, *11*, 65–70. [CrossRef]

26. Ye Yuan, P.R. Enzymatic degradation of human skin dermis revealed by fluorescence and reflectance spectroscopy. *Opt. Exp.* **2008**, *16*, 9857–9868.

27. Borisova, E.G.; Angelova, L.P.; Pavlova, E.P. Endogenous and Exogenous Fluorescence Skin Cancer Diagnostics for Clinical Applications. *IEEE J. Sel. Top. Quantum Electron.* **2014**, *20*, 211–222. [CrossRef]

28. Zhang, H.; Salo, D.C.; Kim, D.M.; Komarov, S.; Tai, Y.C.; Berezin, M.Y. Penetration depth of photons in biological tissues from hyperspectral imaging in shortwave infrared in transmission and reflection geometries. *J. Biomed. Opt.* **2016**, *21*, 126006. [CrossRef]

29. Martelli, F.; Binzoni, T.; Pifferi, A.; Spinelli, L.; Farina, A.; Torricelli, A. There's plenty of light at the bottom: Statistics of photon penetration depth in random media. *Sci. Rep.* **2016**, *6*, 27057. [CrossRef]

30. Hidenobu Arimoto, M.E.; Yamada, Y. Depth profile of diffuse reflectance near-infrared spectroscopy for measurement of water content in skin. *Skin Res. Technol.* **2005**, *11*, 27–35. [CrossRef]

31. Arifler, D.; Schwarz, R.A.; Chang, S.K.; Richards-Kortum, R. Reflectance spectroscopy for diagnosis of epithelial precancer: Model-based analysis of fiber-optic probe designs to resolve spectral information from epithelium and stroma. *Appl. Opt.* **2005**, *44*, 4291–4305. [CrossRef]

32. Tseng, S.H.; Tzeng, S.Y.; Liaw, Y.K.; Hsu, C.K.; Lee, J.; Chen, W.R. Noninvasive evaluation of collagen and hemoglobin contents and scattering property of in vivo keloid scars and normal skin using diffuse reflectance spectroscopy: Pilot study. *J. Biomed. Opt.* **2012**, *17*, 077005. [CrossRef] [PubMed]

33. Nilsson, H.; Larsson, M.; Nilsson, G.; Strömberg, T. Photon pathlength determination based on spatially resolved diffuse reflectance. *J. Biomed. Opt.* **2002**, *7*, 475–485. [CrossRef]

34. Larsson, M.; Nilsson, H.; Strömberg, T. In vivo determination of local skin optical properties and photon path length by use of spatially resolved diffuse reflectance with applications in laser Doppler flowmetry. *Appl. Opt.* **2003**, *42*, 124–134. [CrossRef]

35. Liebert, A.; Wabnitz, H.; Steinbrink, J.; Obrig, H.; Möller, M.; Macdonald, R.; Villringer, A.; Rinneberg, H. Time-resolved multidistance near-infrared spectroscopy of the adult head: Intracerebral and extracerebral absorption changes from moments of distribution of times of flight of photons. *Appl. Opt.* **2004**, *43*, 3037–3047. [CrossRef]
36. Fredriksson, I.; Larsson, M.; Strömberg, T. Inverse Monte Carlo method in a multilayered tissue model for diffuse reflectance spectroscopy. *J. Biomed. Opt.* **2012**, *17*, 047004. [CrossRef] [PubMed]

37. Blondel, W.; Delconte, A.; Khairallah, G.; Marchal, F.; Gavoille, A.; Amouroux, M. Spatially-Resolved Multiply-Excited Autofluorescence and Diffuse Reflectance Spectroscopy: SpectroLive Medical Device for Skin In Vivo Optical Biopsy. *Electronics* **2021**, *10*, 243. [CrossRef]

38. Khairallah, G.; Amouroux, M.; Plénat, F.; Rakotomanga, P.; Soussen, C.; Marchal, F.; Delconte, A.; Chen, H.; Blondel, W. Spatially resolved spectroscopy for guiding margin delineation during human skin carcinomas resection: first clinical results on diffuse reflectance and autofluorescence spectra and in vivo skin optical properties. In *Biophotonics: Photonic Solutions for Better Health Care VI*; Popp, J., Tuchin, V.V., Pavone, F.S., Eds.; International Society for Optics and Photonic: Strasbourg, France, 2018; Volume 10685, pp. 200–206.

39. Dahlstrand, U.; Sheikh, R.; Nguyen, C.D.; Hult, J.; Reistad, N.; Malmsjö, M. Identification of tumor margins using diffuse reflectance spectroscopy with an extended-wavelength spectrum in a porcine mode. *Skin Res. Technol.* **2018**, *24*, 667–671. [CrossRef]

40. Feng, X.; Fox, M.C.; Reichenberg, J.S.; Lopes, F.C.P.S.; Sebastian, K.R.; Markey, M.K.; Tunnell, J.W. Biophysical basis of skin cancer margin assessment using Raman spectroscopy. *Biomed. Opt. Express* **2019**, *10*, 104–118. [CrossRef]

41. Zelinskyi, Y.; Naglič, P.; Pernuš, F.; Likar, B.; Bürmen, M. Fast and accurate Monte Carlo simulations of subdiffusive spatially resolved reflectance for a realistic optical fiber probe tip model aided by a deep neural network. *Biomed. Opt. Express* **2020**, *11*, 3875–3889. [CrossRef] [PubMed]

42. Maiti, R.; Duan, M.; Danby, S.G.; Lewis, R.; Matcher, S.J.; Carré, M.J. Morphological parametric mapping of 21 skin sites throughout the body using optical coherence tomography. *J. Mech. Behav. Biomed. Mater.* **2020**, *102*, 103501. [CrossRef]

43. Chopra, K.; Calva, D.; Sosin, M.; Tadisina, K.; Banda, A.; Cruz, C.D.L.; Chaudhry, M.; Legesse, T.; Drachenberg, C.; Manson, P.; et al. A comprehensive examination of topographic thickness of skin in the human face. *Aesthet. Surg. J.* **2015**, *35*, 1007–1013. [CrossRef] [PubMed]

44. Robertson, K.; Re, J.L. Variation in Epidermal Morphology in Human Skin at Different Body Sites as Measured by Reflectance Confocal Microscopy. *Med. J. Ltd.* **2010**, *90*, 368–373.

45. Amouroux, M. Caractérisation de la Transformation Néoplasique de la peau par Spectroscopies Optiques sur Fantôme de Mélanome et Carcinome Epidermoïde murin Photo-Induit. Ph.D. Thesis, Université Henri Poincaré, Nancy, France, 2008.

46. Colas, V.; Daul, C.; Khairallah, G.; Amouroux, M.; Blondel, W. Spatially resolved diffuse reflectance and autofluorescence photon depth distribution in human skin spectroscopy: A modeling study. In *Optics in Health Care and Biomedical Optics X*; Luo, Q., Li, X., Gu, Y., Zhu, D., Eds.; International Society for Optics and Photonics: Bellingham, Washington, USA, 2020; Volume 11553, pp. 82–96. [CrossRef]

47. Meglinski, I.V.; Matcher, S.J. Quantitative assessment of skin layers absorption and skin reflectance spectra simulation in the visible and near-infrared spectral regions. *Physiol. Meas* **2002**, *23*, 741. [CrossRef]

48. Jacques, S.L. Optical properties of biological tissues: A review. *Phys. Med. Biol.* **2013**, *58*, R37–R61. [CrossRef]

49. Prahl, S. A Compendium of Tissue Optical Properties. Available online: http://omlc.ogi.edu/spectra/ (accessed on 1 December 2020).

50. Karsten, A.E.; Smit, J.E. Modeling and Verification of Melanin Concentration on Human Skin Type. *Photochem. Photobiol.* **2012**, *88*, 469–474. [CrossRef]

51. Nishidate, I.; Maeda, T.; Aizu, Y.; Niizeki, K. Visualizing depth and thickness of a local blood region in skin tissue using diffuse reflectance images. *J. Biomed. Opt.* **2007**, *12*, 054006. [CrossRef]

52. Krishnaswamy, A.; Baranoski, G. *Light and Skin Interactions*, 1st ed.; Elsevier: Burlington, MA, USA, 2010.

53. Altshuler, G.; Smirnov, M.; Yaroslavsky, I. Lattice of optical islets: A novel treatment modality in photomedicine. *J. Phys. D Appl. Phys.* **2005**, *38*, 2732. [CrossRef]

54. Ding, H.; Lu, J.; Wooden, W.; Kragel, P.; Hu, X. Refractive indices of human skin tissues at eight wavelengths and estimated dispersion relations between 300 and 1600 nm. *Phys. Med. Biol.* **2006**, *51*, 1479–1489. [CrossRef]

55. Tsui, S.Y.; Wang, C.Y.; Huang, T.H.; Sung, K.B. Modelling spatially-resolved diffuse reflectance spectra of a multi-layered skin model by artificial neural networks trained with Monte Carlo simulations. *Biomed. Opt. Express* **2018**, *9*, 1531–1544. [CrossRef]

56. Alerstam, E.; Svensson, T.; Andersson-Engels, S. CUDAMCML: User Manual and Implementation Notes. 2009. Available online: http://www.lth.se/fileadmin/atomfysik/Biophotonics/Software/CUDAMCML.pdf (accessed on 1 October 2019).

57. Sun, C.; Aernouts, B.; Saeys, W. Bridging the gap between measurement-based and simulation-based metamodels for deriving bulk optical properties from spatially-resolved reflectance profiles: effect of illumination and detection geometry. *Opt. Express* **2021**, *29*, 15882–15905. [CrossRef]

58. Hennessy, R.; Goth, W.; Sharma, M.; Markey, M.K.; Tunnell, J.W. Effect of probe geometry and optical properties on the sampling depth for diffuse reflectance spectroscopy. *J. Biomed. Opt.* **2014**, *19*, 107002. [CrossRef] [PubMed]

59. Greening, G.; Mundo, A.; Rajaram, N.; Muldoon, T.J. Sampling depth of a diffuse reflectance spectroscopy probe for in-vivo physiological quantification of murine subcutaneous tumor allografts. *J. Biomed. Opt.* **2018**, *23*, 085006. [CrossRef] [PubMed]

60. Yudovsky, D.; Pilon, L. Rapid and accurate estimation of blood saturation, melanin content, and epidermis thickness from spectral diffuse reflectance. *Appl. Opt.* **2010**, *49*, 1707–1719. [CrossRef] [PubMed]
61. Sharma, M.; Hennessy, R.; Markey, M.K.; Tunnell, J.W. Verification of a two-layer inverse Monte Carlo absorption model using multiple source-detector separation diffuse reflectance spectroscopy. *Biomed. Opt. Express* **2014**, *5*, 40–53. [CrossRef] [PubMed]

62. Feng, S.; Zeng, F.A.; Chance, B. Photon migration in the presence of a single defect: A perturbation analysis. *Appl. Opt.* **1995**, *34*, 3826–3837. [CrossRef]

63. Jonasson, H.; Fredriksson, I.; Bergstrand, S.; Östgren, C.J.; Larsson, M.; Strömberg, T. In vivo characterization of light scattering properties of human skin in the 475- to 850-nm wavelength range in a Swedish cohort. *J. Biomed. Opt.* **2018**, *23*, 121608. [CrossRef]

64. Jacques, S.L. Optical assessment of cutaneous blood volume depends on the vessel size distribution: A computer simulation study. *J. Biophotonics* **2010**, *3*, 75–81.

65. Zhu, C.; Liu, Q. Validity of the semi-infinite tumor model in diffuse reflectance spectroscopy for epithelial cancer diagnosis: a Monte Carlo study. *Opt. Express* **2011**, *19*, 17799–17812. [CrossRef] [PubMed]

66. Palmer, G.M.; Ramanujam, N. Monte Carlo-based inverse model for calculating tissue optical properties. Part I: Theory and validation on synthetic phantoms. *Appl. Opt.* **2006**, *45*, 1062–1071. [CrossRef] [PubMed]

67. Cen, H.; Lu, R.; Dolan, K. Optimization of inverse algorithm for estimating the optical properties of biological materials using spatially-resolved diffuse reflectance. *Inverse Probl. Sci. Eng.* **2010**, *18*, 853–872. [CrossRef]