Next Article in Journal
Multidisciplinary Forensic Approach in “Complex” Bodies: Systematic Review and Procedural Proposal
Previous Article in Journal
Advanced Time-Frequency Methods for ECG Waves Recognition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of the Differential Pathlength Factor for Human Skin Using Monte Carlo Simulations

Biomedical Engineering Department, College of Engineering, Imam Abdulrahman Bin Faisal University, Dammam 31441, Saudi Arabia
Diagnostics 2023, 13(2), 309; https://doi.org/10.3390/diagnostics13020309
Submission received: 10 November 2022 / Revised: 28 December 2022 / Accepted: 11 January 2023 / Published: 13 January 2023
(This article belongs to the Section Biomedical Optics)

Abstract

:
Near-infrared technology is an emerging non-invasive technique utilized for various medical applications. Recently, there have been many attempts to utilize NIR technology for the continues monitoring of blood glucose levels through the skin. Different approaches and designs have been proposed for non-invasive blood glucose measurements. Light photons penetrating the skin can undergo multiple scattering events, and the actual optical pathlength becomes larger than the source-to-detector separation (optode spacing) in the reflection-mode configuration. Thus, the differential pathlength factor (DPF) must be incorporated into the modified Beer–Lambert law. The accurate estimation of the DPF values will lead to an accurate quantification of the physiological variations within the tissue. In this work, the aim was to systematically estimate the DPF for human skin for a range of source-to-detector separations and wavelengths. The Monte Carlo (MC) method was utilized to mimic the different layers of human skin with different optical properties and blood and water volume fractions. This work could help improve the accuracy of the near-infrared technique in the measurement of physiological variations within skin tissue.

1. Introduction

Near-infrared (NIR) is an emerging non-invasive technique utilized for various medical applications. In the past two decades, NIR has been heavily investigated in the applications of breast and brain imaging. For breast imaging, NIR has been used to non-invasively detect breast tumors [1,2,3,4,5,6]. Moreover, NIR was successfully used to non-invasively differentiate between benign and malignant tumors [7,8]. Recently, comparisons have been made with optical systems and traditional breast imaging modalities [4,5,6]. To accurately localize breast lesions, NIR is normally combined with other medical imaging modalities such as ultrasound, X-ray, or magnetic resonance imaging (MRI) [9,10]. For brain imaging, NIR is used to assess brain activities [11,12]. It is now considered a potential alternative to existing brain-imaging modalities such as functional MRI and electroencephalography (EEG), since it is more compact and cheaper. Research groups have assessed multi-modality applications with functional NIR. Regular combinations with functional NIR include EEG or fMRI to concurrently measure brain activity [13]. NIR also has many potential applications in the brain–computer interface [14,15,16]. Recently, there have been many attempts to utilize NIR technology for the continuous monitoring of blood glucose levels [17,18]. Different approaches and designs have been proposed for non-invasive blood glucose measurement through the skin [19,20,21].
The transmission and reflection modes of NIR setups have been investigated [22]. The NIR reflection mode is preferred for thick and dense tissue, while the transmission mode is better for thin tissue. The reflection mode is most common, since the majority of non-invasive NIR applications involve thick and dense tissue. NIR systems for both brain and breast imaging are examples of the NIR reflection mode [7,23]. Blood glucose NIR sensors are also based on the reflection mode, since the glucose contents are found in the skin dermis layer. In a previous work [24], we proposed dual NIR sensors, wherein two channels were used to measure both the signal of interest and the interfering noise arising from the superficial skin layer. Furthermore, the factors affecting the NIR penetration depth and the suitable source-to-detector separation (SDS) for a range of NIR wavelengths were investigated. Recently, attempts were made by researchers [25] to combine both the NIR transmission and NIR reflection configurations for blood glucose measurements. The design used two reflection-based channels and one transmission-based channel. A deep neural network model was then used to predict the glucose blood contents.
NIR technology relies on illuminating the human tissue with a light source and detecting the escaping penetrating light at a distance. The physiological variations within the tissue are related to the light attenuation (A) according to the modified Beer–Lambert law [26]. For the tissue medium, the Beer–Lambert law can be stated as follows:
A λ = ε λ × C λ × d × D P F
where ε is the molar absorption coefficient, C is the concentration of the absorbing particles, and d is the source-to-detector separation. Depending on the optical properties of the different layers of the tissue and the operating wavelength λ of the light source, photons can experience multiple scattering events. Thus, the actual optical pathlength is much larger than source-to-detector separation (optode spacing), and so the differential pathlength factor (DPF) was incorporated into the above equation to relate the actual optical pathlength to the source-to-detector separation [27,28].
Many studies have estimated the optical pathlength and the DPF computationally or experimentally [27,28,29,30]. Studies have shown that the DPF can be estimated directly from direct time-of-flight measurements [27]. Phase-resolved optical spectroscopy has also been utilized to estimate the DPF [29,31]. Computationally, the Monte Carlo method was utilized to estimate the DPF in reflectance photoplethysmography [28] and blood volume and oxygen saturation [30] for two wavelengths. Moreover, researchers have investigated the DPF’s dependence on optode spacing, operating wavelength, and age [32,33]. These studies showed that the DPF depends on the source wavelength and the optical properties. Thus, different tissues have different DPFs at different wavelengths. The optical pathlength increases as the optode spacing (the source-to-detector distance) increases. The results show that if the optode spacing is greater than 2.5 cm, the DPF can be assumed to be constant [28]. Knowing the DPF is of high importance in the quantification of physiological variations within tissue, such as changes in the oxygen saturation in the blood. The significance of the accurate determination of the DPF has been reported for non-invasive near-infrared applications [34,35].
The aim of this study was to systematically investigate the DPF for human skin. More specifically, this work investigated the DPF’s relationship to different wavelengths, source-to-detectors distances, and skin layers. Additionally, we assessed the effect of different skin colors on the DPF values. The Monte Carlo (MC) method was utilized to mimic the different layers of human skin. The adopted MC model allowed us to consider the anatomical variations between skin layers, such as the blood and water volume fractions, melanin concentrations, and optical properties of the skin layers. This work could help improve the accuracy of near-infrared techniques in the measurement of physiological variations within tissue, such as blood glucose levels.

2. Methods

2.1. Human Skin Model

The MC skin model adopted in this study is illustrated in Figure 1. The model comprised a semi-infinite slab of seven skin layers. Each skin layer had different optical properties, including the absorption coefficient μa, scattering coefficient μs, and anisotropy factor g. In the three-dimensional cartesian coordinate system presented in Figure 1, the light source was placed at the center (0,0,0) of the model, while the detector was placed at a distance (d,0,0) away from the source. The simulation was run for a range of distances (d; 0.5 mm to 8 mm, with a step size of 0.5 mm). The maximum distance considered here (d = 8 mm) was enough for the NIR to reach the deep blood net dermis of the skin with an acceptable signal-to-noise ratio (SNR) [24]. In particular, the SNR for the shorter wavelengths (550 nm and 650 nm) dropped very quickly as the source-to-detector distance increased. Longer wavelengths (750 nm to 1050 nm) provided better SNRs at longer source-to-detector distances, which meant that deeper tissue layers could be measured [24]. For each distance (d), simulations were also conducted for a range of light-source wavelengths. The wavelength range in this model represented the so-called “diagnostic window” of the light spectrum, 450 nm to 1050 nm, with a step size of 100 nm. This window of the spectrum is characterized by its low water-absorption coefficient and, since water is the main absorber in the body, it allows the light to penetrate deeper into the tissue [12,36]. Moreover, the light in this spectrum range is usually utilized to measure oxygenated and deoxygenated blood, since their absorption coefficients are higher than those of water.
This skin model consisted of several layers. The optical properties and the thickness of each layer used to run the simulations are summarized in Table 1. The absorption coefficients of each layer were calculated as explained in [24]. In brief, the blood and water volume fractions Vblood and VH2O and the absorption coefficients for water, deoxyhemoglobin, and oxyhemoglobin for each layer were used to estimate the layer absorption coefficients. The epidermis layer, which is the skin’s outermost superficial layer, contains melanin, which is responsible for the color of the skin. The absorption coefficient of the melanin in the epidermis was estimated as follows [26,27]:
μ a m e l λ = 6.6 × 10 10 × λ 3.33
Thus, the absorption coefficient for the epidermis was evaluated as:
μ a e p i d e r m i s λ = V m e l   μ a m e l λ + V H 2 O   μ a H b O 2 λ + 1 V m e l + V H 2 O       μ a o t h e r λ
In this study, the melanin volume fraction V m e l was assumed to be 2% and 20%, representing people of light- and dark-colored skin, respectively.

2.2. Monte Carlo Simulation

The method of MC simulation is detailed in [31,32]. The light source is modeled as a pencil source directed in the z-direction (Figure 1). Light photons propagate in the tissue, and the reflected (escaped) photons are detected by a detector, which is placed on the same side of the source. In principle, the light photons diffuse in the tissue, following a “banana-shaped” path between the source and the detector [6,33]. Thus, the longer the source-to-detector separation, the deeper the photons can penetrate into the tissue and the higher the photons’ propagation pathlength. In an MC simulation, millions of photons are launched (here: 100 × 10 6 photons) at the start of each simulation run. Each photon packet at its initial launch is given a unity weight ( w = 1 ). The weight and the propagation of each photon is calculated based on the optical properties of the tissue. Between each successive scattering event, the photon free pathlength l probability distribution was calculated as:
p l = μ s e μ s . l  
During scattering, the direction of the photon is determined by two randomly generated angles, and φ , expressed as:
c o s   = 1 2 g 1 + g 2 1 g 2 1 g + 2 g ξ
φ = 2 π ξ
After scattering, part of the photon weight is absorbed by the tissue layer, Δ w = μ a μ a + μ s . If the photon packet has enough energy to reach the boundary of the next tissue layer, it will either be reflected or transmitted. Otherwise, the photon will be terminated. In the case that the photon is transmitted and exits the tissue layer at the detector site (d,0,0), the photon is detected. The optical path of the detected photon is also registered [36,37,38,39,40]. The DPF was calculated from the mean optical path l of the detected photon packets at wavelength λ for the source–detector separation d as follows:
D P F λ , d = l λ d
As explained above, this study considered some important parameters. First, we considered the effect of different NIR light-source wavelengths (450 to 1050 nm) on the calculation of the DPF. We also considered a range of source–detector separations (d = 0.5 to 8 mm). In addition, the effect of different skin colors (different melanin concentrations) on the DPF was studied. The melanin concentrations considered in this study were 2% and 20%, representing people of light- and dark-colored skin, respectively. In each scenario, the simulations were repeated (N = 20) by independently seeded MC simulations [40,41,42,43,44].

3. Results and Discussion

The optical pathlength of the photons for different skin layers with a melanin concentration of 2% are presented in Figure 2. This figure shows the optical pathlength for different wavelengths and skin layers. The figure also indicates the effect of the source-to-detector separation on the optical pathlength.
In addition, Figure 3 present the optical pathlength of the photons for skin with a 20% melanin concentration. The two figures are similar, showing that the impact of skin color on optical pathlength was not significant. The figures also show that as the wavelength increased from 450 nm to 1050 nm, the optical pathlength clearly decreased. In addition, as the source-to-detector separation increased, there was a clear increase in the optical pathlength. The combination of a short wavelength and a long source-to-detector separation provided the highest optical pathlength values. Similarly, the combination of a longer wavelength and a short source-to-detector separation presented the lowest values of the photon optical pathlength. Thus, in estimating the DPF, it is important to consider the effect of both the wavelength and the source-to-detector separation. It is also important to consider the impact of the wavelength and source-to-detector separation on the DPF when designing an NIR sensor. More specifically, to penetrate deeper tissue layers, one would opt for a longer wavelength.
In addition, the two figures show the details of the optical pathlength of photon propagation in each layer of the skin. In all skin layers, the shorter the wavelength, the higher the photon optical pathlength values. Clearly, the optical pathlength had higher values in the deeper layers of the skin. In comparing the different skin layers (Figure 2 and Figure 3), it is worth noting that the layers most affected were the stratum corneum and the epidermis of the skin. Thus, the effect of the skin color on the optical pathlength was not significant, since the epidermis is very thin in comparison to the other layers.
The optical pathlength of photons was larger than the source-to-detector separation. The larger the source-to-detector separation, the higher the penetration depth in the tissue. Figure 4 shows the spatial sensitivity profile of the photon distributions for three different source-to-detector separations: 2 mm, 4 mm, and 6 mm. This figure illustrates that the source-to-detector range used in this study was sufficient to penetrate into the deep blood net dermis, which is a necessary layer for blood content measurement, including glucose. The source-to-detector range verses the depth of penetration of different skin layers is discussed in [24].
Figure 5 shows how the optical pathlength changed for different source-to-detector separations. The figure also shows the number of trials run (N = 20) for each simulation. The number of trials refers to the number repetitions of independently seeded MC simulations. As expected, the standard deviation was very small for the short source-to-detector separations and became larger as the source-to-detector separation increased. Nevertheless, it is clear from the figure that the standard deviation was good enough to clearly differentiate between the different skin layers. Besides the source-to-detector separation, the standard deviation in Figure 5 was also affected by the number of photons launched in the MC simulations. For a lower number of launched photons, the standard deviation tended to be higher for all distances. A high number of launched photons provided more accurate simulations with a low standard deviation. This study included different numbers of launched photons: 10 × 10 5 ,   10 × 10 6 ,   and   100 × 10 6 (the results are not shown here). Launching 100 × 10 6 photons provided a low enough standard deviation for all wavelengths and source-to-detector separations (as shown in Figure 5).
Table 2 summarizes the estimated DPF values calculated based on Equation (7). The table shows the DPFs for the different investigated wavelengths. The wavelength was a significant factor when estimating the DPF values. Different wavelengths had different DPFs at different source-to-detector separations. As the source-to-detector separation increased, the DPF increased for short wavelengths (450 nm and 550 nm). For longer wavelengths, it was found that increasing the source-to-detector separation decreased the DPF values.

4. Conclusions

In conclusion, this study aimed to investigate the DPF for human skin for a range of source-to-detector separations and wavelengths. The study considered the NIR reflection configuration, since it is the most suited for the majority of non-invasive NIR applications, which are designed for thick and dense tissue. Light photons penetrating the skin can undergo multiple scattering events, and the actual optical pathlength becomes larger than the source-to-detector separation in the reflection-mode configuration. Thus, the DPF must be incorporated in the modified Beer–Lambert law. Several important factors were considered, such as the wavelength, the source-to-detector separation, and the skin color. The wavelength range in this model represented the so-called “diagnostic window” of the light spectrum. The MC simulation method was utilized to mimic the different layers of human skin. The adopted MC model allowed us to consider anatomical variations between the skin layers, such as the blood and water volume fractions, melanin concentration, and optical properties of the skin layers.
It was found that the wavelength and the source-to-detector separations were the most important factors affecting the optical pathlength. Shorter wavelengths had the strongest effect on DPFs when changing the source-to-detector separations. It was also found that the color of the skin was not a significant factor effecting the optical pathlength. Table 2 summarizes the estimated DPF values for the different source-to-detector separations and light-source wavelengths. These values could be used and incorporated into the modified Beer–Lambert law according to the configuration of the NIR sensor. The accurate estimation of the DPF values will lead to the accurate quantification of physiological variations within the tissue. This work could guide researchers in designing an NIR sensor for human skin by providing a better understanding of the actual photon pathlengths in the skin.

Funding

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education of Saudi Arabia, for funding this research under the project number IF-2022-046-Eng at Imam Abdulrahman bin Faisal University/College of Engineering.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in FigShare at https://doi.org/10.6084/m9.figshare.21895716.v1.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boas, D.; Brooks, D.; Miller, E.; DiMarzio, C.; Kilmer, M.; Gaudette, R.; Zhang, Q. Imaging the body with diffuse optical tomography. IEEE Signal Process. Mag. 2001, 18, 57–75. [Google Scholar] [CrossRef]
  2. Guven, M.; Yazici, B.; Intes, X.; Chance, B. Diffuse optical tomography with a priori anatomical information. Phys. Med. Biol. 2005, 50, 2837–2858. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Üncü, Y.A.; Sevim, G.; Canpolat, M. Approaches to preclinical studies with heterogeneous breast phantom using reconstruction and three-dimensional image processing algorithms for diffuse optical imaging. Int. J. Imaging Syst. Technol. 2022, 32, 343–353. [Google Scholar] [CrossRef]
  4. Üncü, Y.A.; Sevim, G.; Mercan, T.; Vural, V.; Durmaz, E.; Canpolat, M. Differentiation of tumoral and non-tumoral breast lesions using back reflection diffuse optical tomography: A pilot clinical study. Int. J. Imaging Syst. Technol. 2021, 31, 2023–2031. [Google Scholar] [CrossRef]
  5. Uddin, K.M.S.; Zhang, M.; Anastasio, M.; Zhu, Q. Optimal breast cancer diagnostic strategy using combined ultrasound and diffuse optical tomography. Biomed. Opt. Express 2020, 11, 2722–2737. [Google Scholar] [CrossRef]
  6. Chae, E.Y.; Kim, H.H.; Sabir, S.; Kim, Y.; Kim, H.; Yoon, S.; Ye, J.C.; Cho, S.; Heo, D.; Kim, K.H.; et al. Development of digital breast tomosynthesis and diffuse optical tomography fusion imaging for breast cancer detection. Sci. Rep. 2020, 10, 13127. [Google Scholar] [CrossRef]
  7. Vavadi, H.; Mostafa, A.; Zhou, F.; Uddin, K.M.S.; Althobaiti, M.; Xu, C.; Bansal, R.; Ademuyiwa, F.; Poplack, S.; Zhu, Q. Compact ultrasound-guided diffuse optical tomography system for breast cancer imaging. J. Biomed. Opt. 2019, 24, 021203. [Google Scholar] [CrossRef]
  8. Althobaiti, M.; Vavadi, H.; Zhu, Q. An Automated Preprocessing Method for Diffuse Optical Tomography to Improve Breast Cancer Diagnosis. Technol. Cancer Res. Treat. 2018, 17, 1533033818802791. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, Q.; Brukilacchio, T.J.; Li, A.; Stott, J.; Chaves, T.; Hillman, E.M.C.; Wu, T.; Chorlton, M.; Rafferty, E.; Moore, R.H.; et al. Coregistered tomographic X-ray and optical breast imaging: Initial results. J. Biomed. Opt. 2005, 10, 24033–240339. [Google Scholar] [CrossRef] [Green Version]
  10. Brooksby, B.; Jiang, S.; Dehghani, H.; Pogue, B.; Paulsen, K.D.; Weaver, J.; Kogel, C.; Poplack, S.P. Combining near-infrared tomography and magnetic resonance imaging to study in vivo breast tissue: Implementation of a Laplacian-type regularization to incorporate magnetic resonance structure. J. Biomed. Opt. 2005, 10, 51504–51510. [Google Scholar] [CrossRef] [Green Version]
  11. Gibson, A.; Hebden, J.C.; Arridge, S. Recent advances in diffuse optical imaging. Phys. Med. Biol. 2005, 50, R1–R43. [Google Scholar] [CrossRef]
  12. Althobaiti, M.; Al-Naib, I. Recent Developments in Instrumentation of Functional Near-Infrared Spectroscopy Systems. Appl. Sci. 2020, 10, 6522. [Google Scholar] [CrossRef]
  13. Ortega, P.; Faisal, A.A. Deep learning multimodal fNIRS and EEG signals for bimanual grip force decoding. J. Neural Eng. 2021, 18, 0460e6. [Google Scholar] [CrossRef]
  14. Khan, R.A.; Naseer, N.; Qureshi, N.K.; Noori, F.M.; Nazeer, H.; Khan, M.U. FNIRS-based Neurorobotic Interface for gait rehabilitation. J. Neuroeng. Rehabil. 2018, 15, 7. [Google Scholar] [CrossRef] [Green Version]
  15. Rea, M.; Rana, M.; Lugato, N.; Terekhin, P.; Gizzi, L.; Brötz, D.; Caria, A. Lower limb movement preparation in chronic stroke: A pilot study toward an fnirs-bci for gait rehabilitation. Neurorehabil. Neural Repair 2014, 28, 564–575. [Google Scholar] [CrossRef]
  16. Almulla, L.; Al-Naib, I.; Ateeq, I.S.; Althobaiti, M. Observation and motor imagery balance tasks evaluation: An fNIRS feasibility study. PLoS ONE 2022, 17, e0265898. [Google Scholar] [CrossRef]
  17. Delbeck, S.; Vahlsing, T.; Leonhardt, S.; Steiner, G.; Heise, H.M. Non-invasive monitoring of blood glucose using optical methods for skin spectroscopy—Opportunities and recent advances. Anal. Bioanal. Chem. 2019, 411, 63–77. [Google Scholar] [CrossRef]
  18. Alsunaidi, B.; Althobaiti, M.; Tamal, M.; Albaker, W.; Al-Naib, I. A Review of Non-Invasive Optical Systems for Continuous Blood Glucose Monitoring. Sensors 2021, 21, 6820. [Google Scholar] [CrossRef]
  19. Rachim, V.P.; Chung, W.-Y. Wearable-band type visible-near infrared optical biosensor for non-invasive blood glucose monitoring. Sensors Actuators B Chem. 2019, 286, 173–180. [Google Scholar] [CrossRef]
  20. Haxha, S.; Jhoja, J. Optical Based Noninvasive Glucose Monitoring Sensor Prototype. IEEE Photon. J. 2016, 8, 1–11. [Google Scholar] [CrossRef]
  21. Althobaiti, M.; Al-Naib, I. Optimization of Dual-Channel Near-Infrared Non-Invasive Glucose Level Measurement Sensors Based On Monte-Carlo Simulations. IEEE Photon. J. 2021, 13, 1–9. [Google Scholar] [CrossRef]
  22. Shokrekhodaei, M.; Quinones, S. Review of non-invasive glucose sensing techniques: Optical, electrical and breath acetone. Sensors 2020, 20, 1251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Fantini, S.; Sassaroli, A. Frequency-Domain Techniques for Cerebral and Functional Near-Infrared Spectroscopy. Front. Neurosci. 2020, 14, 300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Althobaiti, M. In Silico Investigation of SNR and Dermis Sensitivity for Optimum Dual-Channel Near-Infrared Glucose Sensor Designs for Different Skin Colors. Biosensors 2022, 12, 805. [Google Scholar] [CrossRef] [PubMed]
  25. Joshi, A.M.; Jain, P.; Mohanty, S.P.; Agrawal, N. iGLU 2.0: A New Wearable for Accurate Non-Invasive Continuous Serum Glucose Measurement in IoMT Framework. IEEE Trans. Consum. Electron. 2020, 66, 327–335. [Google Scholar] [CrossRef]
  26. Kocsis, L.; Herman, P.; Eke, A. The modified Beer-Lambert law revisited. Phys. Med. Biol. 2006, 51, N91. [Google Scholar] [CrossRef]
  27. Delpy, D.T.; Cope, M.; Van Der Zee, P.; Arridge, S.; Wray, S.; Wyatt, J. Estimation of optical pathlength through tissue from direct time of flight measurement. Phys. Med. Biol. 1988, 33, 1433–1442. [Google Scholar] [CrossRef] [Green Version]
  28. Chatterjee, S.; Abay, T.Y.; Phillips, J.P.; Kyriacou, P.A. Investigating optical path and differential pathlength factor in reflectance photoplethysmography for the assessment of perfusion. J. Biomed. Opt. 2018, 23, 075005. [Google Scholar] [CrossRef] [Green Version]
  29. Duncan, A.; Meek, J.H.; Clemence, M.; E Elwell, C.; Tyszczuk, L.; Cope, M.; Delpy, D. Optical pathlength measurements on adult head, calf and forearm and the head of the newborn infant using phase resolved optical spectroscopy. Phys. Med. Biol. 1995, 40, 295–304. [Google Scholar] [CrossRef] [Green Version]
  30. Chatterjee, S.; Kyriacou, P.A. Estimating the Dependence of Differential Pathlength Factor on Blood Volume and Oxygen Saturation using Monte Carlo method. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBS, Berlin, Germany, 23–27 July 2019. [Google Scholar]
  31. Delpy, D.T.; Cope, M. Quantification in tissue near–infrared spectroscopy. Philos. Trans. R. Soc. B Biol. Sci. 1997, 352, 649–659. [Google Scholar] [CrossRef]
  32. Scholkmann, F.; Wolf, M. General equation for the differential pathlength factor of the frontal human head depending on wavelength and age. J. Biomed. Opt. 2013, 18, 105004. [Google Scholar] [CrossRef] [Green Version]
  33. Kohl, M.; Nolte, C.; Heekeren, H.; Horst, S.; Scholz, U.; Obrig, H.; Villringer, A. Determination of the wavelength dependence of the differential pathlength factor from near-infrared pulse signals. Phys. Med. Biol. 1998, 43, 1771–1782. [Google Scholar] [CrossRef]
  34. Boas, D.A.; Gaudette, T.; Strangman, G.; Cheng, X.; Marota, J.J.A.; Mandeville, J.B. The accuracy of near infrared spectroscopy and imaging during focal changes in cerebral hemodynamics. Neuroimage 2001, 13, 76–90. [Google Scholar] [CrossRef] [Green Version]
  35. Uludağ, K.; Steinbrink, J.; Villringer, A.; Obrig, H. Separability and cross talk: Optimizing dual wavelength combinations for near-infrared spectroscopy of the adult head. Neuroimage 2004, 22, 583–589. [Google Scholar] [CrossRef]
  36. Jacques, S.L. Optical properties of biological tissues: A review. Phys. Med. Biol. 2013, 58, R37. [Google Scholar] [CrossRef]
  37. Meglinski, I.; Matcher, S. Computer simulation of the skin reflectance spectra. Comput. Methods Programs Biomed. 2003, 70, 179–186. [Google Scholar] [CrossRef] [Green Version]
  38. Nishidate, I.; Aizu, Y.; Mishina, H. Estimation of melanin and hemoglobin in skin tissue using multiple regression analysis aided by Monte Carlo simulation. J. Biomed. Opt. 2004, 9, 700–710. [Google Scholar] [CrossRef]
  39. Chatterjee, S.; Budidha, K.; A Kyriacou, P. Investigating the origin of photoplethysmography using a multiwavelength Monte Carlo model. Physiol. Meas. 2020, 41, 084001. [Google Scholar]
  40. Kessoku, S.; Maruo, K.; Okawa, S.; Masamoto, K.; Yamada, Y. Influence of blood glucose level on the scattering coefficient of the skin in near-infrared spectroscopy. In ASME/JSME 2011 8th Thermal Engineering Joint Conference, AJTEC 2011; ASME: Honolulu, HI, USA, 2011. [Google Scholar]
  41. Fang, Q.; Boas, D.A. Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units. Opt. Express 2009, 17, 20178–20190. [Google Scholar] [CrossRef] [Green Version]
  42. Yan, S.; Fang, Q. Hybrid mesh and voxel based Monte Carlo algorithm for accurate and efficient photon transport modeling in complex bio-tissues. Biomed. Opt. Express 2020, 11, 6262–6270. [Google Scholar]
  43. Brigadoi, S.; Cooper, R. How short is short? Optimum source–detector distance for short-separation channels in functional near-infrared spectroscopy. Neurophotonics 2015, 2, 025005. [Google Scholar] [CrossRef] [Green Version]
  44. Yuan, Y.; Yu, L.; Doğan, Z.; Fang, Q. Graphics processing units-accelerated adaptive nonlocal means filter for denoising three-dimensional Monte Carlo photon transport simulations. J. Biomed. Opt. 2018, 23, 121618. [Google Scholar] [CrossRef]
Figure 1. Illustration of the Monte Carlo model for human skin adopted in this study. The light source was placed at the center (0,0,0) of the model, and the detector was placed at a distance (d,0,0) away from the source.
Figure 1. Illustration of the Monte Carlo model for human skin adopted in this study. The light source was placed at the center (0,0,0) of the model, and the detector was placed at a distance (d,0,0) away from the source.
Diagnostics 13 00309 g001
Figure 2. Optical pathlength of light photons for different skin layers at 2% melanin concentration. The figure shows the effects of different light source wavelengths and source-to-detector separations on the optical pathlength of photons.
Figure 2. Optical pathlength of light photons for different skin layers at 2% melanin concentration. The figure shows the effects of different light source wavelengths and source-to-detector separations on the optical pathlength of photons.
Diagnostics 13 00309 g002
Figure 3. Optical pathlength of the light photons for different skin layers at 20% melanin concentration. The figure shows the effects of different light source wavelengths and source-to-detector separations on the optical pathlength of photons.
Figure 3. Optical pathlength of the light photons for different skin layers at 20% melanin concentration. The figure shows the effects of different light source wavelengths and source-to-detector separations on the optical pathlength of photons.
Diagnostics 13 00309 g003
Figure 4. Spatial sensitivity profile of the photon distributions for three different source-to-detector separations ((a) 2 mm, (b) 4 mm, and (c) 6 mm) for the wavelength 1050 nm. Photons follow a so-called “banana-shaped” path, which showed that the optical pathlength increased when the source-to-detector separation increased.
Figure 4. Spatial sensitivity profile of the photon distributions for three different source-to-detector separations ((a) 2 mm, (b) 4 mm, and (c) 6 mm) for the wavelength 1050 nm. Photons follow a so-called “banana-shaped” path, which showed that the optical pathlength increased when the source-to-detector separation increased.
Diagnostics 13 00309 g004
Figure 5. Optical pathlength for different layers at the wavelength of 550 nm. The MC simulations were run independently seeded for N = 20. The number of launched photons in the MC simulations was 100 × 10 6 .
Figure 5. Optical pathlength for different layers at the wavelength of 550 nm. The MC simulations were run independently seeded for N = 20. The number of launched photons in the MC simulations was 100 × 10 6 .
Diagnostics 13 00309 g005
Table 1. Values used in the Monte Carlo skin model based on [28,29,30] for different skin layers.
Table 1. Values used in the Monte Carlo skin model based on [28,29,30] for different skin layers.
Skin LayerThickness (mm)Scattering Coefficient (1/mm) VbloodVH2O
μ s , 450 n m μ s , 550 n m μ s , 650 n m μ s , 750 n m μ s , 850 n m μ s , 950 n m μ s , 1050 n m
Stratum Corneum 0.0260.6141.7422.3315.0513.329.57.5200.05
Epidermis0.25 00.2
Papillary dermis 0.10.040.5
Upper blood net dermis0.080.30.6
Reticular dermis0.20.040.7
Deep blood net dermis 0.30.10.7
Subcutaneous tissue20.050.7
Table 2. Estimated DPF based on Equation (7) for different light source wavelengths and source-to-detector separations.
Table 2. Estimated DPF based on Equation (7) for different light source wavelengths and source-to-detector separations.
SDS (mm) D P F 450 n m D P F 550 n m D P F 650 n m D P F 750 n m D P F 850 n m D P F 950 n m D P F 1050 n m
0.59.3272297.6136237.7345518.4458648.7316999.3213859.629079
114.266369.0170286.8491516.7429566.8275697.243327.628147
1.518.9221910.797736.9260076.1715036.0975956.1924366.473498
223.2021712.523197.2226625.9914585.7910095.6130655.772583
2.526.8585314.187097.5791915.9859235.6457955.2659595.288481
330.0962715.640197.9718146.0288395.612515.0414514.947292
3.532.8337116.987098.3314246.1026815.5961674.8969954.705665
435.2065518.157318.6670866.1934065.6231614.7935384.519836
4.537.3074719.131739.0212676.2975075.6622554.731434.384798
539.0632420.00429.2775626.3769035.6926624.655144.272771
5.540.2948820.772229.5411986.4395535.7316034.6255154.171025
641.3953821.380139.7488796.5252365.7732414.5819584.096541
6.542.3118322.131039.8851776.5589095.7777824.5359694.004704
743.0074822.3984910.019526.6235015.7975354.509063.953183
7.543.0623922.6481610.169626.628685.8281024.4747143.897391
843.2189622.8091110.24326.6363245.8132294.4451113.844256
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Althobaiti, M. Estimation of the Differential Pathlength Factor for Human Skin Using Monte Carlo Simulations. Diagnostics 2023, 13, 309. https://doi.org/10.3390/diagnostics13020309

AMA Style

Althobaiti M. Estimation of the Differential Pathlength Factor for Human Skin Using Monte Carlo Simulations. Diagnostics. 2023; 13(2):309. https://doi.org/10.3390/diagnostics13020309

Chicago/Turabian Style

Althobaiti, Murad. 2023. "Estimation of the Differential Pathlength Factor for Human Skin Using Monte Carlo Simulations" Diagnostics 13, no. 2: 309. https://doi.org/10.3390/diagnostics13020309

APA Style

Althobaiti, M. (2023). Estimation of the Differential Pathlength Factor for Human Skin Using Monte Carlo Simulations. Diagnostics, 13(2), 309. https://doi.org/10.3390/diagnostics13020309

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

Article Metrics

Back to TopTop