Estimation of melanin and hemoglobin in skin tissue using multiple regression analysis aided by Monte Carlo simulation

Engineering vascularized skin-mimetic phantom for non-invasive Raman spectroscopy

Recent advances in Raman spectroscopy have shown great potential for non-invasive analyte sensing, but the lack of a standardized optical phantom for these measurements has hindered further progress. While many research groups have developed optical phantoms that mimic bulk optical absorption and scattering, these materials typically have strong Raman scattering, making it difficult to distinguish metabolite signals. As a result, solid tissue phantoms for spectroscopy have been limited to highly scattering tissues such as bones and calcifications, and metabolite sensing has been primarily performed using liquid tissue phantoms. To address this issue, we have developed a layered skin-mimetic phantom that can support metabolite sensing through Raman spectroscopy. Our approach incorporates millifluidic vasculature that mimics blood vessels to allow for diffusion akin to metabolite diffusion in the skin. Furthermore, our skin phantoms are mechanically mimetic, providing an ideal model for development of minimally invasive optical techniques. By providing a standardized platform for measuring metabolites, our approach has the potential to facilitate critical developments in spectroscopic techniques and improve our understanding of metabolite dynamics in vivo.

Measuring and imaging of transcutaneous bilirubin, hemoglobin, and melanin based on diffuse reflectance spectroscopy

Significance

Evaluation of biological chromophore levels is useful for detection of various skin diseases, including cancer, monitoring of health status and tissue metabolism, and assessment of clinical and physiological vascular functions. Clinically, it is useful to assess multiple different chromophores in vivo with a single technique or instrument.

Aim

To investigate the possibility of estimating the concentration of four chromophores, bilirubin, oxygenated hemoglobin, deoxygenated hemoglobin, and melanin from diffuse reflectance spectra in the visible region.

Approach

A new diffuse reflectance spectroscopic method based on the multiple regression analysis aided by Monte Carlo simulations for light transport was developed to quantify bilirubin, oxygenated hemoglobin, deoxygenated hemoglobin, and melanin. Three different experimental animal models were used to induce hyperbilirubinemia, hypoxemia, and melanogenesis in rats.

Results

The estimated bilirubin concentration increased after ligation of the bile duct and reached around 18    mg / dl at 50 h after the onset of ligation, which corresponds to the reference value of bilirubin measured by a commercially available transcutaneous bilirubin meter. The concentration of oxygenated hemoglobin and that of deoxygenated hemoglobin decreased and increased, respectively, as the fraction of inspired oxygen decreased. Consequently, the tissue oxygen saturation dramatically decreased. The time course of melanin concentration after depilation of skin on the back of rats was indicative of the supply of melanosomes produced by melanocytes of hair follicles to the growing hair shaft.

Conclusions

The results of our study showed that the proposed method is capable of the in vivo evaluation of percutaneous bilirubin level, skin hemodynamics, and melanogenesis in rats, and that it has potential as a tool for the diagnosis and management of hyperbilirubinemia, hypoxemia, and pigmented skin lesions.

Deep learning-based optical approach for skin analysis of melanin and hemoglobin distribution

SIGNIFICANCE

Melanin and hemoglobin have been measured as important diagnostic indicators of facial skin conditions for aesthetic and diagnostic purposes. Commercial clinical equipment provides reliable analysis results, but it has several drawbacks: exclusive to the acquisition system, expensive, and computationally intensive.

AIM

We propose an approach to alleviate those drawbacks using a deep learning model trained to solve the forward problem of light-tissue interactions. The model is structurally extensible for various light sources and cameras and maintains the input image resolution for medical applications.

APPROACH

A facial image is divided into multiple patches and decomposed into melanin, hemoglobin, shading, and specular maps. The outputs are reconstructed into a facial image by solving the forward problem over skin areas. As learning progresses, the difference between the reconstructed image and input image is reduced, resulting in the melanin and hemoglobin maps becoming closer to their distribution of the input image.

RESULTS

The proposed approach was evaluated on 30 subjects using the professional clinical system, VISIA VAESTRO. The correlation coefficients for melanin and hemoglobin were found to be 0.932 and 0.857, respectively. Additionally, this approach was applied to simulated images with varying amounts of melanin and hemoglobin.

CONCLUSION

The proposed approach showed high correlation with the clinical system for analyzing melanin and hemoglobin distribution, indicating its potential for accurate diagnosis. Further calibration studies using clinical equipment can enhance its diagnostic ability. The structurally extensible model makes it a promising tool for various image acquisition conditions.

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

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.

Discrete emotions discovered by contactless measurement of facial blood flows

ABSTRACTExperiential and behavioural aspects of emotions can be measured readily but developing a contactless measure of emotions' physiological aspects has been a major challenge. We hypothesised that different emotion-evoking films can produce distinctive facial blood flow patterns that can serve as physiological signatures of discrete emotions. To test this hypothesis, we created a new Transdermal Optical Imaging system that uses a conventional video camera to capture facial blood flows in a contactless manner. Using this and deep machine learning, we analysed videos of the faces of people as they viewed film clips that elicited joy, sadness, disgust, fear or a neutral state. We found that each of these elicited a distinct blood flow pattern in the facial epidermis, and that Transdermal Optical Imaging is an effective contactless and inexpensive tool to the reveal physiological correlates of discrete emotions.

In Silico Investigation of SNR and Dermis Sensitivity for Optimum Dual-Channel Near-Infrared Glucose Sensor Designs for Different Skin Colors

Diabetes is a serious health condition that requires patients to regularly monitor their blood glucose level, making the development of practical, compact, and non-invasive techniques essential. Optical glucose sensors-and, specifically, NIR sensors-have the advantages of being non-invasive, compact, inexpensive, and user-friendly devices. However, these sensors have low accuracy and are yet to be adopted by healthcare providers. In our previous work, we introduced a non-invasive dual-channel technique for NIR sensors, in which a long channel is utilized to measure the glucose level in the inner skin (dermis) layer, while a short channel is used to measure the noise signal of the superficial skin (epidermis) layer. In this work, we investigated the use of dual-NIR channels for patients with different skin colors (i.e., having different melanin concentrations). We also adopted a Monte Carlo simulation model that takes into consideration the differences between different skin layers, in terms of blood content, water content, melanin concentration in the epidermis layer, and skin optical proprieties. On the basis of the signal-to-noise ratio, as well as the sensitivities of both the epidermis and dermis layers, we suggest the selection of wavelengths and source-to-detector separation for optimal NIR channels under different skin melanin concentrations. This work facilitates the improved design of a compact and non-invasive NIR glucose sensor that can be utilized by patients with different skin colors.

Repeatability and reproducibility of a hyperspectral imaging system for in vivo color evaluation

BACKGROUND

Color imaging is a tried and true method for the evaluation of cosmetic and dermatological effects, but it fails to capture all the information in a scene's spectral reflectance. For this reason, there has been in recent years increasing interest in the use of imaging spectrometers for clinical studies and product evaluation.

MATERIAL AND METHODS

We developed a novel HyperSpectral Imager (HSI) able to take in vivo full-face format images as a next generation instrument for skin color measurement and beyond. Here, we report part of the results of our first full-scale validation test of the HSI. We replicated a make-up foundation screening test by applying three products to a panel of 9 models and evaluated the product L^∗ , a^∗ , b^∗ , and ∆E effect immediately after application relative to the bare skin condition. We repeated this test twice in order to study the repeatability of the HSI as an evaluation instrument and during each test two different operators duplicated the data acquisition so we can assess the reproducibility of the measurements.

RESULTS

We find that the measurements from the HSI provide repeatability and reproducibility as good or better than those of our previous benchmark devices.

CONCLUSION

From these results, we conclude that not only is the HSI suitable for use in color evaluation studies, but also that it gives operational advantages over the previous generation of evaluation instruments, as it provides a spectral measurement combined with good spatial resolution. This allows for analysis of color over an area and post hoc selection of study regions and so opens new possibilities for studies of complex in vivo phenomena which neither non-imaging spectrometers nor conventional cameras can pursue. This study also raises points for future work concerning proper inclusion of instrument uncertainty in comparisons of results between instruments and handling of systematic uncertainties from analyses based on a single area.

An ultrasonography-based approach for tissue modelling to inform photo-therapy treatment strategies

Currently, diagnostic medicine uses a multitude of tools ranging from ionising radiation to histology analysis. With advances in piezoelectric crystal technology, high-frequency ultrasound imaging has developed to achieve comparatively high resolution without the drawbacks of ionising radiation. This research proposes a low-cost, non-invasive and real-time protocol for informing photo-therapy procedures using ultrasound imaging. We combine currently available ultrasound procedures with Monte Carlo methods for assessing light transport and photo-energy deposition in the tissue. The measurements from high-resolution ultrasound scans are used as input for optical simulations. Consequently, this provides a pipeline that will inform medical practitioners for better therapy strategy planning. While validating known inferences of light transport through biological tissue, our results highlight the range of information such as temporal monitoring and energy deposition at varying depths. This process also retains the flexibility of testing various wavelengths for individual-specific geometries and anatomy.

Visualizing veins from color images under varying illuminations for medical applications

SIGNIFICANCE

Effective vein visualization is critically important for several clinical procedures, such as venous blood sampling and intravenous injection. Existing technologies using infrared device or ultrasound rely on professional equipment and are not suitable for daily medical care. A regression-based vein visualization method is proposed.

AIM

We visualize veins from conventional RGB images to provide assistance in venipuncture procedures as well as clinical diagnosis of some venous insufficiency.

APPROACH

The RGB images taken by digital cameras are first transformed to spectral reflectance images using Wiener estimation. Multiple regression analysis is then applied to derive the relationship between spectral reflectance and the concentrations of pigments. Monte Carlo simulation is adopted to get prior information. Finally, vein patterns are visualized from the spatial distribution of pigments. To minimize the effect of illumination on skin color, light correction and shading removal operations are performed in advance.

RESULTS

Experimental results from inner forearms of 60 subjects show the effectiveness of the regression-based method. Subjective and objective evaluations demonstrate that the clarity and completeness of vein patterns can be improved by light correction and shading removal.

CONCLUSIONS

Vein patterns can be successfully visualized from RGB images without any professional equipment. The proposed method can assist in venipuncture procedures. It also shows promising potential to be used in clinical diagnosis and treatment of some venous insufficiency.

Design and performance of a hyperspectral camera for full-face in vivo imaging

Red, green, blue color photography is a mature technology and a powerful tool for the evaluation and understanding of the way an object reflects light and its related optical properties, but color photography fails to give a complete picture of these effects due to its inherent lack of spectral resolution. In this work, we update the L'OREAL reference device for skin color measurement, the Chromasphere, by replacing its current color camera system with an imaging spectrometer. This imaging spectrometer must provide a spatial resolution on par with the previous color cameras and a spectral resolution commensurate with a spectroradiometer while also achieving a time resolution suitable for in vivo studies of the human face. Due to these requirements, common spatial scanning techniques are not suitable for this application, and so we utilized a spectral-scanning approach based on a tunable liquid-crystal birefringent filter. We present the design and performance tests of a working prototype that is capable of measuring the spectrum in each of 4 MP with a nominal spectral resolution of 10 nm across the wavelength range from 420 to 730 nm in a total imaging time of less than 10 s. We cross-compared the spectral and color measurements obtained with this prototype, an industry-standard spectroradiometer, and a charge-coupled device color camera in order to assess the prototype's performance, and the results of this comparison show that our prototype is capable of taking spectral measurements near enough in quality to those of a spectroradiometer to successfully bridge the divide between such devices and conventional color cameras. Doing so, this instrument opens new possibilities for studies of complex in vivo phenomena that neither non-imaging spectrometers nor conventional cameras can pursue.

In Vivo Transcutaneous Monitoring of Hemoglobin Derivatives Using a Red-Green-Blue Camera-Based Spectral Imaging Technique

Cyanosis is a pathological condition that is characterized by a bluish discoloration of the skin or mucous membranes. It may result from a number of medical conditions, including disorders of the respiratory system and central nervous system, cardiovascular diseases, peripheral vascular diseases, deep vein thrombosis, and regional ischemia. Cyanosis can also be elicited from methemoglobin. Therefore, a simple, rapid, and simultaneous monitoring of changes in oxygenated hemoglobin and deoxygenated hemoglobin is useful for protective strategies against organ ischemic injury. We previously developed a red-green-blue camera-based spectral imaging method for the measurements of melanin concentration, oxygenated hemoglobin concentration (C_HbO), deoxygenated hemoglobin concentration (C_HbR), total hemoglobin concentration (C_HbT) and tissue oxygen saturation (StO₂) in skin tissues. We leveraged this approach in this study and extended it to the simultaneous quantifications of methemoglobin concentration (C_metHb), C_HbO, C_HbR, and StO₂. The aim of the study was to confirm the feasibility of the method to monitor C_metHb, C_HbO, C_HbR, C_HbT, and StO₂. We performed in vivo experiments using rat dorsal skin during methemoglobinemia induced by the administration of sodium nitrite (NaNO₂) and changing the fraction of inspired oxygen (FiO₂), including normoxia, hypoxia, and anoxia. Spectral diffuse reflectance images were estimated from an RGB image by the Wiener estimation method. Multiple regression analysis based on Monte Carlo simulations of light transport was used to estimate C_HbO, C_HbR, C_metHb, C_HbT, and StO₂. C_metHb rapidly increased with a half-maximum time of less than 30 min and reached maximal values nearly 60 min after the administration of NaNO₂, whereas StO₂ dramatically dropped after the administration of NaNO₂, indicating the temporary production of methemoglobin and severe hypoxemia during methemoglobinemia. Time courses of C_HbT and StO₂, while changing the FiO₂, coincided with well-known physiological responses to hyperoxia, normoxia, and hypoxia. The results indicated the potential of this method to evaluate changes in skin hemodynamics due to loss of tissue viability and vitality.

Transcutaneous monitoring of hemoglobin derivatives during methemoglobinemia in rats using spectral diffuse reflectance

SIGNIFICANCE

Untreated methemoglobinemia may cause severe hypoxemia and even death when methemoglobin levels in the blood stream exceed 70%. Although CO-oximetry can be used to monitor the response to treatment for methemoglobinemia, it is costly and requires an invasive procedure for collecting blood samples from patients. A pulse CO-oximeter with a contact probe can be used to continuously and non-invasively measure the percentage of methemoglobin, as well as the percutaneous oxygen saturation. In terms of the prevention of infectious diseases, however, it is desirable to monitor methemoglobin and oxygen saturation levels in a non-contact manner. Diffuse reflectance spectral imaging is promising as a non-contact, non-invasive, and cost-effective clinical diagnostic tool for methemoglobinemia.

AIM

To demonstrate the feasibility of visible spectral diffuse reflectance for in vivo monitoring of hemoglobin derivatives and evaluating methemoglobin production and reduction as well as hypoxemia during methemoglobinemia in rats.

APPROACH

A new imaging approach based on the multiple regression analysis aided by Monte Carlo simulations for light transport was developed to quantify methemoglobin, oxygenated hemoglobin, and deoxygenated hemoglobin using a hyperspectral imaging system. An in vivo experiment with rats exposed to sodium nitrite (NaNO2) at different doses was performed to confirm the feasibility of the method for evaluating the dynamics of methemoglobin, oxygenated hemoglobin, and deoxygenated hemoglobin during methemoglobinemia. Systemic physiological parameters, including the percutaneous arterial oxygen saturation, heart rate (HR), and pulse distention, were measured by a commercially available pulse oximeter, and the results were compared to those obtained by the proposed method.

RESULTS

Both the methemoglobin concentration and methemoglobin saturation rapidly increased with a half-maximum time of <20  min. They reached their maximal values nearly 60 min after the administration of NaNO2. Tissue oxygen saturation dramatically dropped to a minimum of 33.7  %    ±  0.4  %  , 23.1  %    ±  5.6  %  , 8.8  %    ±  1.7  %  , and 9.7  %    ±  5.1  %   on average for NaNO2 doses of 25, 37.5, 50, and 75 mg/kg, respectively. Changes in methemoglobin concentration and tissue oxygen saturation are indicative of the temporary production of methemoglobin and severe hypoxemia during methemoglobinemia. Profound increases in the HR and pulse distention implied an elevated cardiac output caused by tachycardia and the resultant increase in peripheral blood volume to compensate for the hypoxia and hypoxemia during methemoglobinemia. This was in agreement with the time course of the peripheral hemoglobin volume concentration obtained by the proposed method.

CONCLUSIONS

The proposed method is capable of the in vivo non-contact simultaneous evaluation of methemoglobin levels and hypoxemia during methemoglobinemia, and that it has potential as a tool for the diagnosis and management of methemoglobinemia.

Phase Velocity of Facial Blood Volume Oscillation at a Frequency of 0.1 Hz

Facial blood flow, which typically exhibits distinctive oscillation at a frequency of around 0.1 Hz, has been extensively studied. Although this oscillation may include important information about blood flow regulation, its origin remains unknown. The spatial phase distribution of the oscillation is thus desirable. Therefore, we visualized facial blood volume oscillation at a frequency of around 0.1 Hz using a digital camera imaging method with an improved approximation equation, which enabled precise analysis over a large area. We observed a slow spatial movement of the 0.1-Hz oscillation. The oscillation phase was not synchronized, but instead moved slowly. The phase velocity varies with person, measurement location, and time. An average phase velocity of 3.8 mm/s was obtained for several subjects. The results are consistent with previous studies; however, the conventional explanation that the blood flow at a certain point oscillates independently of adjacent areas should be corrected. If the primary origin of the movement is myogenic activity, the movement may ascend along a blood vessel toward the upstream. Otherwise, the oscillation and its propagation can be considered to be related to Mayer waves. By determining the mechanism, some questions regarding Mayer waves can be answered. The direction of the wave (upstream or downstream) provides important information.

Simple and affordable imaging of multiple physiological parameters with RGB camera-based diffuse reflectance spectroscopy

We propose a simple and affordable imaging technique to evaluate transcutaneously multiple physiological parameters by using a digital red-green-blue camera. In this method, the RGB-values were converted into tristimulus values in the CIE (Commission Internationale de l'Eclairage) XYZ color space, which is compatible with the common color spaces. Monte Carlo simulation for light transport in biological tissue was then performed to specify the relationship among the XYZ-values and the concentrations of oxygenated hemoglobin, deoxygenated hemoglobin, bilirubin, and melanin. The concentration of total hemoglobin and tissue oxygen saturation were also calculated from the estimated concentrations of oxygenated and deoxygenated hemoglobin. In vivo experiments with bile duct ligation in rats demonstrated that the estimated bilirubin concentration increased after ligation of the bile duct and reached around 22 mg/dl at 116 h after the onset of ligation, which corresponds to the ground truth value of bilirubin measured by a commercially available transcutaneous bilirubinometer. Experiments with rats while varying the fraction of inspired oxygen demonstrated that oxygenated hemoglobin and deoxygenated hemoglobin decreased and increased, respectively, as the fraction of inspired oxygen decreased. Consequently, tissue oxygen saturation dramatically decreased. We further extended the method to a non-contact imaging photo-plethysmograph and estimation of the percutaneous oxygen saturation. An empirical formula to estimate percutaneous oxygen saturation was derived from the pulse wave amplitudes of oxygenated and deoxygenated hemoglobin. The estimated percutaneous oxygen saturation dropped remarkably when a faction of inspired oxygen was below 19%, indicating the onset of hypoxemia due to hypoxia, whereas the tissue oxygen saturation decreased gradually according to the reduction of the faction of inspired oxygen. The results in this study indicate the potential of this method for imaging of multiple physiological parameters in skin tissue and evaluating an optical biomedical imaging technique that enables cost-effective, easy-to-use, portable, remotely administered, and/or point-of-care solutions.

Three-dimensional maps of human skin properties on full face with shadows using 3-D hyperspectral imaging

Hyperspectral imaging has shown great potential for optical skin analysis by providing noninvasive, pixel-by-pixel surface measurements from which, applying an optical model, information such as melanin concentration and total blood volume fraction can be mapped. Such applications have been successfully performed on small flat skin areas, but existing methods are not suited to large areas such as an organ or a face, due to the difficulty of ensuring homogeneous illumination on complex three-dimensional (3-D) objects, which leads to errors in the maps. We investigate two methods to account for these irradiance variations on a face. The first one relies on a radiometric correction of the irradiance, using 3-D information on the face's shape acquired by combining the hyperspectral camera with a 3-D scanner; the second relies on an optimization metric used in the map computation, which is invariant to irradiance. We discuss the advantages and drawbacks of the two methods, after having presented in detail the whole acquisition setup, which has been designed to provide high-resolution images with a short acquisition time, as required for live surface measurements of complex 3-D objects such as the face.

Validation of an Inverse Fitting Method of Diffuse Reflectance Spectroscopy to Quantify Multi-Layered Skin Optical Properties

Skin consists of epidermis and dermis layers that have distinct optical properties. The quantification of skin optical properties is commonly achieved by modeling photon propagation in tissue using Monte Carlo (MC) simulations and iteratively fitting experimentally measured diffuse reflectance spectra. In order to speed up the inverse fitting process, time-consuming MC simulations have been replaced by artificial neural networks to quickly calculate reflectance spectra given tissue geometric and optical parameters. In this study the skin was modeled to consist of three layers and different scattering properties of the layers were considered. A new inverse fitting procedure was proposed to improve the extraction of chromophore-related information in the skin, including the hemoglobin concentration, oxygen saturation and melanin absorption. The performance of the new inverse fitting procedure was evaluated on 40 sets of simulated spectra. The results showed that the fitting procedure without knowing the epidermis thickness extracted chromophore information with accuracy similar to or better than fitting with known epidermis thickness, which is advantageous for practical applications due to simpler and more cost-effective instruments. In addition, the melanin volume fraction multiplied by the thickness of the melanin-containing epidermis layer was estimated more accurately than the melanin volume fraction itself. This product has the potential to provide a quantitative indicator of melanin absorption in the skin. In-vivo cuff occlusion experiments were conducted and skin optical properties extracted from the experiments were comparable to the results of previously reported in vivo studies. The results of the current study demonstrated the applicability of the proposed method to quantify the optical properties related to major chromophores in the skin, as well as scattering coefficients of the dermis. Therefore, it has the potential to be a useful tool for quantifying skin optical properties in vivo.

Fundus-simulating phantom for calibration of retinal vessel oximetry devices

Retinal vessel oxygen supply is important for retinal tissue metabolism. Commonly used retinal vessel oximetry devices are based on dual-wavelength spectral measurement of oxyhemoglobin and deoxyhemoglobin. However, there is no traceable standard for reliable calibration of these devices. In this study, we developed a fundus-simulating phantom that closely mimicked the optical properties of human fundus tissues. Microchannels of precisely controlled topological structures were produced by soft lithography to simulate the retinal vasculature. Optical properties of the phantom were adjusted by adding scattering and absorption agents to simulate different concentrations of fundus pigments. The developed phantom was used to calibrate the linear correlation between oxygen saturation (SO₂) level and optical density ratio in a dual-wavelength oximetry device. The obtained calibration factors were used to calculate the retinal vessel SO₂ in both eyes of five volunteers aged between 24 and 27 years old. The test results showed that the mean arterial and venous SO₂ levels after phantom calibration were coincident with those after empirical value calibration, indicating the potential clinical utility of the produced phantom as a calibration standard.

Noninvasive evaluation of hemodynamics and light scattering property during two-stage mouse cutaneous carcinogenesis based on multispectral diffuse reflectance images at isosbestic wavelengths of hemoglobin

We investigate a multispectral imaging method to evaluate spatiotemporal changes in both cutaneous hemoglobin concentration and light scattering parameter in mouse skin through diffuse reflectance spectroscopy using the reflectance images acquired at isosbestic wavelengths of hemoglobin (420, 450, 500, and 585 nm). In the proposed approach, Monte Carlo simulation-based empirical formulas are introduced to extract the scattering power b representing the wavelength dependence of light scattering spectrum of skin tissue, as well as the total hemoglobin concentration Cth in dermal vasculatures. The use of isosbestic wavelengths of hemoglobin enables the values of Cth and b to be estimated independently of the oxygenation of hemoglobin. Experiments using in vivo mice two-stage chemical carcinogenesis model are performed to confirm the feasibility of the proposed method for evaluating the changes in cutaneous vasculatures and tissue morphology during tumor initiation, promotion, and progression processes. The experimental results reveal that the changes in scattering power b of back skin are significantly reduced and followed by the increase in total hemoglobin concentration Cth in the carcinogenesis mice group, which indicates morphological changes in skin tissue such as edema and cell swelling caused by tumor promotion and successive angiogenesis along with tumor progression. The results suggest that the potential of the present method to detect cutaneous carcinogenesis in an early stage and monitor physiological changes during promotion and progression process of nonmelanoma tumors.

Agarose-based Tissue Mimicking Optical Phantoms for Diffuse Reflectance Spectroscopy

This protocol describes how to make agarose-based tissue-mimicking phantoms and demonstrates how to determine their optical properties using a conventional optical system with an integrating sphere. Measuring systems for the acquisition of the diffuse reflectance and total transmittance spectra are constructed with a broadband white light source, a light guide, an achromatic lens, an integrating sphere, a sample holder, an optical fiber probe, and a multi-channel spectrometer. An acrylic mold consisting of two rectangular acrylic pieces and a U-shaped acrylic piece is constructed to create an epidermal phantom and a dermal phantom with whole blood. The application of a sodium dithionite (Na2S2O4) solution to the dermal phantom enables the researcher to deoxygenate hemoglobin in red blood cells distributed in the dermal phantom. The inverse Monte Carlo simulation with the diffuse reflectance and total transmittance spectra measured by a spectrometer with an integrating sphere is performed to determine the absorption coefficient spectrum µa(λ) and the reduced scattering coefficient spectrum µs'(λ) of each layer phantom. A two-layered phantom mimicking the diffuse reflectance of human skin tissue is also demonstrated by piling up the epidermal phantom on the dermal phantom.

Transdermal Optical Imaging Reveal Basal Stress via Heart Rate Variability Analysis: A Novel Methodology Comparable to Electrocardiography

The present study examined the validity of a novel physiological measurement technology called transdermal optical imaging (TOI) technology at assessing basal stress. This technology conveniently, contactlessly, and remotely measures facial blood flow changes using a conventional digital video camera. We compared data from TOI against the pulse data collected from the FDA approved BIOPAC system. One hundred thirty-six healthy adults participated in the study. We found that TOI measurements of heart rate and heart rate variability (HRV), which reflects basal stress, corresponded strongly to those obtained from BIOPAC. These findings indicate that TOI technology is a viable method to monitor heart rate and HRV not only accurately but also conveniently, contactlessly, and remotely. Further, measures of HRV obtained via TOI serves as a valid index of basal stress. Potential applications of this technology in psychological research and other fields are discussed.

Subcutaneous veins depth measurement using diffuse reflectance images

Intravenous (IV) procedures are often difficult due to the poor visualization of subcutaneous veins. Because existing vein locators lack the ability to assess depth, and also because mis-punctures and poor vascular access remain problematic, we propose an imaging system that employs diffuse reflectance images at three isosbestic wavelengths to measure both the depth and thickness of subcutaneous veins. This paper describes the proposed system as well as proof-of-principle experimental demonstrations. We initially introduce the working principle and structure of the system. All measurements were based on the Monte Carlo (MC) method and accomplished by referring an optical density (OD) ratio to a multi-layer diffuse reflectance model. Results were all validated by comparative ultrasound measurements. Experimental trials included 11 volunteers who were subjected to both ultrasound measurements and the proposed optical process to validate the system's applicability. However, the unreliability of the "thickness" measurement of the vein may be due to the fact that the veins have collapsible walls - so excess pressure by the transducer will give a false thickness.

Skin collagen can be accurately quantified through noninvasive optical method: Validation on a swine study

BACKGROUND/PURPOSE

Diffuse reflectance spectroscopy (DRS) is a noninvasive optical technology characterized by relatively low system cost and high efficiency. In our previous study, we quantified the relative concentration of collagen for the individual keloid patient. However, no actual value of collagen concentration can prove the reliability of collagen detection by our DRS system.

METHODS

Skin-mimicking phantoms were prepared using different collagen and coffee concentrations, and their chromophore concentrations were quantified using the DRS system to analyze the influence of collagen and other chromophores. Moreover, we used the animal study to compare the DRS system with the collagen evaluation of biopsy section by second-harmonic generation (SHG) microscopy at four different skin parts.

RESULTS

In the phantom study, the result showed that coffee chromophore did not severely interfere with collagen concentration recovery. In the animal study, a positive correlation (r=.902) between the DRS system and collagen evaluation with SHG microscopy was found.

CONCLUSIONS

We have demonstrated that the DRS system can quantify the actual values of collagen concentration and excluded the interference of other chromophores in skin-mimicking phantoms. Furthermore, a high positive correlation was found in the animal study with SHG microscopy. We consider that the DRS is a potential technique and can evaluate skin condition objectively.

Evaluation of Cerebral Hemodynamics and Tissue Morphology of In Vivo Rat Brain Using Spectral Diffuse Reflectance Imaging

We investigated a quantitative imaging of reduced scattering coefficients μ_s'( λ) and the absorption coefficients μ_a( λ) of in vivo cortical tissues in the range from visible to near-infrared (NIR) wavelengths based on diffuse reflectance spectral imaging technique. In this method, diffuse reflectance images of in vivo cortical tissue are acquired at nine wavelengths (500, 520, 540, 560, 570, 580, 600, 730, and 760 nm). A multiple regression analysis aided by the Monte Carlo simulation for the absorbance spectra is then utilized to estimate the optical coefficients of cortical tissue. This analysis calculates the concentration of oxygenated hemoglobin and that of deoxygenated hemoglobin, the scattering amplitude a and the scattering power b. The spectrum of absorption coefficient is deduced from the estimated concentrations of oxygenated hemoglobin and deoxygenated hemoglobin. The spectrum of reduced scattering coefficient is determined by the estimated scattering amplitude and scattering power. The particle size distribution of microstructure is calculated from the estimated scattering power b for evaluating the morphological change in brain tissue quantitatively. Animal experiments with in vivo exposed brain of rats demonstrated that the responses of the absorption properties to hyperoxic and anoxic conditions are in agreement with the expected well-known cortical hemodynamics. The average particle size was significantly reduced immediately after the onset of anoxia and then it was changed into an increase, which implied the swelling and shrinkage of the cellular and subcellular structures induced by loss of tissue viability in brain tissue.

Spectrophotometric analysis at the single-cell level: elucidating dispersity within melanic immortalized cell populations

The widely held notion that melanin-containing cells are uniform in both size and optical characteristics is demonstrably false.

Robust near real-time estimation of physiological parameters from megapixel multispectral images with inverse Monte Carlo and random forest regression

PURPOSE

Multispectral imaging can provide reflectance measurements at multiple spectral bands for each image pixel. These measurements can be used for estimation of important physiological parameters, such as oxygenation, which can provide indicators for the success of surgical treatment or the presence of abnormal tissue. The goal of this work was to develop a method to estimate physiological parameters in an accurate and rapid manner suited for modern high-resolution laparoscopic images.

METHODS

While previous methods for oxygenation estimation are based on either simple linear methods or complex model-based approaches exclusively suited for off-line processing, we propose a new approach that combines the high accuracy of model-based approaches with the speed and robustness of modern machine learning methods. Our concept is based on training random forest regressors using reflectance spectra generated with Monte Carlo simulations.

RESULTS

According to extensive in silico and in vivo experiments, the method features higher accuracy and robustness than state-of-the-art online methods and is orders of magnitude faster than other nonlinear regression based methods.

CONCLUSION

Our current implementation allows for near real-time oxygenation estimation from megapixel multispectral images and is thus well suited for online tissue analysis.

Fabrication of various optical tissue phantoms by the spin-coating method

Although numerous studies have been performed to fabricate various optical tissue phantom (OTP) models, the fabrication of OTPs that simulate skin layers is laborious and time-consuming owing to the intricate characteristics of skin tissue. This study presents various OTP models that optically and structurally simulate the epidermis–dermis skin layer. The spin-coating method was employed to reproduce a uniform thin layer that mimics the epidermis layer, and the fabrication parameters were optimized for epoxy and silicone reference materials. Various OTP models simulating blood vessels and hyperpigmentation lesions were fabricated using the two reference materials to determine their feasibility. The suitability of each of the two reference materials for OTP fabrication was qualitatively evaluated by comparing the quality of the OTP models.

Rapid calculation of diffuse reflectance from a multilayered model by combination of the white Monte Carlo and adding-doubling methods

To rapidly derive a result for diffuse reflectance from a multilayered model that is equivalent to that of a Monte-Carlo simulation (MCS), we propose a combination of a layered white MCS and the adding-doubling method. For slabs with various scattering coefficients assuming a certain anisotropy factor and without absorption, we calculate the transition matrices for light flow with respect to the incident and exit angles. From this series of precalculated transition matrices, we can calculate the transition matrices for the multilayered model with the specific anisotropy factor. The relative errors of the results of this method compared to a conventional MCS were less than 1%. We successfully used this method to estimate the chromophore concentration from the reflectance spectrum of a numerical model of skin and in vivo human skin tissue.

Estimation of melanin and hemoglobin using spectral reflectance images reconstructed from a digital RGB image by the Wiener estimation method

A multi-spectral diffuse reflectance imaging method based on a single snap shot of Red-Green-Blue images acquired with the exposure time of 65 ms (15 fps) was investigated for estimating melanin concentration, blood concentration, and oxygen saturation in human skin tissue. The technique utilizes the Wiener estimation method to deduce spectral reflectance images instantaneously from an RGB image. Using the resultant absorbance spectrum as a response variable and the extinction coefficients of melanin, oxygenated hemoglobin and deoxygenated hemoglobin as predictor variables, multiple regression analysis provides regression coefficients. Concentrations of melanin and total blood are then determined from the regression coefficients using conversion vectors that are numerically deduced in advance by the Monte Carlo simulations for light transport in skin. Oxygen saturation is obtained directly from the regression coefficients. Experiments with a tissue-like agar gel phantom validated the method. In vivo experiments on fingers during upper limb occlusion demonstrated the ability of the method to evaluate physiological reactions of human skin.

Update on techniques for the quantitation of facial skin characteristics

The purpose of this article is to review the strategies and methods for quantifying treatment outcomes, perhaps defined by the consumer/patient as a "decrease in perceived age." The demand for the rejuvenation of facial skin is expected to increase as the population ages and seeks optimal outcomes from the array of available treatment modalities. This information will be of value to the plastic surgeon in collaborating with patients on evaluation and treatment strategies.

In vivo spatial frequency domain spectroscopy of two layer media

Monitoring of tissue blood volume and local oxygen saturation can inform the assessment of tissue health, healing, and dysfunction. These quantities can be estimated from the contribution of oxyhemoglobin and deoxyhemoglobin to the absorption spectrum of the dermis. However, estimation of blood related absorption in skin can be confounded by the strong absorption of melanin in the epidermis and epidermal thickness and pigmentation varies with anatomic location, race, gender, and degree of disease progression. Therefore, a method is desired that decouples the effect of melanin absorption in the epidermis from blood absorption in the dermis for a large range of skin types and thicknesses. A previously developed inverse method based on a neural network forward model was applied to simulated spatial frequency domain reflectance of skin for multiple wavelengths in the near infrared. It is demonstrated that the optical thickness of the epidermis and absorption and reduced scattering coefficients of the dermis can be determined independently and with minimal coupling. Then, the same inverse method was applied to reflectance measurements from a tissue simulating phantom and in vivo human skin. Oxygen saturation and total hemoglobin concentrations were estimated from the volar forearms of weakly and strongly pigmented subjects using a standard homogeneous model and the present two layer model.

A review of in-vivo optical properties of human tissues and its impact on PDT

A thorough understanding of optical properties of biological tissues is critical to effective treatment planning for therapies such as photodynamic therapy (PDT). In the last two decades, new technologies, such as broadband diffuse spectroscopy, have been developed to obtain in vivo data in humans that was not possible before. We found that the in vivo optical properties generally vary in the ranges μ(a) = 0.03-1.6 cm⁻¹ and μ'(s) = 1.2-40 cm⁻¹, although the actual range is tissue-type dependent. We have also examined the overall trend of the absorption spectra (for μ(a) and μ'(s)) as a function of wavelength within a 95% confidence interval for various tissues in vivo. The impact of optical properties on light fluence rate is also discussed for various light application geometries including superficial, interstitial, and within a cavity.

Noninvasive spectral imaging of skin chromophores based on multiple regression analysis aided by Monte Carlo simulation

In order to visualize melanin and blood concentrations and oxygen saturation in human skin tissue, a simple imaging technique based on multispectral diffuse reflectance images acquired at six wavelengths (500, 520, 540, 560, 580 and 600 nm) was developed. The technique utilizes multiple regression analysis aided by Monte Carlo simulation for diffuse reflectance spectra. Using the absorbance spectrum as a response variable and the extinction coefficients of melanin, oxygenated hemoglobin, and deoxygenated hemoglobin as predictor variables, multiple regression analysis provides regression coefficients. Concentrations of melanin and total blood are then determined from the regression coefficients using conversion vectors that are deduced numerically in advance, while oxygen saturation is obtained directly from the regression coefficients. Experiments with a tissue-like agar gel phantom validated the method. In vivo experiments with human skin of the human hand during upper limb occlusion and of the inner forearm exposed to UV irradiation demonstrated the ability of the method to evaluate physiological reactions of human skin tissue.

Noninvasive imaging of human skin hemodynamics using a digital red-green-blue camera

In order to visualize human skin hemodynamics, we investigated a method that is specifically developed for the visualization of concentrations of oxygenated blood, deoxygenated blood, and melanin in skin tissue from digital RGB color images. Images of total blood concentration and oxygen saturation can also be reconstructed from the results of oxygenated and deoxygenated blood. Experiments using tissue-like agar gel phantoms demonstrated the ability of the developed method to quantitatively visualize the transition from an oxygenated blood to a deoxygenated blood in dermis. In vivo imaging of the chromophore concentrations and tissue oxygen saturation in the skin of the human hand are performed for 14 subjects during upper limb occlusion at 50 and 250 mm Hg. The response of the total blood concentration in the skin acquired by this method and forearm volume changes obtained from the conventional strain-gauge plethysmograph were comparable during the upper arm occlusion at pressures of both 50 and 250 mm Hg. The results presented in the present paper indicate the possibility of visualizing the hemodynamics of subsurface skin tissue.

Skin color independent assessment of aging using skin autofluorescence

Skin autofluorescence (AF) for the non-invasive assessment of the amount of accumulated tissue Advanced Glycation Endproducts (AGEs) increases with aging. In subjects with darker skin colors, measurements typically result in lower AF values than in subjects with fair skin colors, e.g. due to selective absorption by skin compounds. Our aim was to provide a new method for calculating skin AF, yielding values that are independent of skin color. The deviation of skin AF of healthy subjects with various darker skin types (N = 99) compared to reference values from Caucasians showed to be a function of various parameters that were derived from reflectance and emission spectra in the UV and visible range (adjusted R(2) = 80%). Validation of the new algorithm, based on these findings, in a separate dataset (N = 141) showed that results of skin AF can now be obtained to assess skin AGEs independently of skin color.

Determination of optical properties of turbid media spanning visible and near-infrared regimes via spatially modulated quantitative spectroscopy

We present a novel, noncontact method for the determination of quantitative optical properties of turbid media from 430 to 1050 nm. Through measuring the broadband reflectance from an unknown sample as a function of the spatial frequency of the projected illumination patterns, the absolute absorption and reduced scattering coefficients can be calculated without a priori assumptions of the chromophores present. This technique, which is called spatially modulated quantitative spectroscopy (SMoQS), was validated through the quantification of optical properties of homogenous liquid phantoms with known concentrations of absorbers and scatterers. The properties of the phantoms were recovered across the range of values prepared with R(2) values of 0.985 and 0.996 for absorption and reduced scattering, respectively. A measurement was also performed on skin tissue as a demonstration of the method's performance in vivo. The resultant absorption spectrum was well described by a multichromophore fit, and the quantitative values for oxy- and deoxyhemoglobin, water, and melanin were within published ranges for skin.

Chromophore concentrations, absorption and scattering properties of human skin in-vivo

Absorption and reduced scattering coefficients of in-vivo human skin provide critical information on non-invasive skin diagnoses for aesthetic and clinical purposes. To date, very few in-vivo skin optical properties have been reported. Previously, we reported absorption and scattering properties of in-vivo skin in the wavelength range from 650 to 1000 nm using the diffusing probe in the "modified two-layer geometry". In this study, we determine the spectra of skin optical properties continuously in the range from 500 to 1000 nm. It was found that the concentration of chromophores, such as oxy-hemoglobin, deoxy-hemoglobin, and melanin, calculated based on the absorption spectra of eighteen subjects at wavelengths above and below 600 nm were distinct because of the inherent difference in the interrogation region. The scattering power, which is related to the average scatterer's size, demonstrates a clear contrast between skin phototypes, skin sites, and wavelengths. We also applied venous occlusion on forearms and found that the concentrations of oxy- and deoxy-hemoglobin as assessed at wavelengths above and below 600 nm were different. Our results suggest that diffuse reflectance techniques with the visible and near infrared light sources can be employed to investigate the hemodynamics and optical properties of upper dermis and lower dermis.

Visualizing of skin chromophore concentrations by use of RGB images

A method is proposed for visualizing simply the concentrations of melanin, oxygenated blood, and deoxygenated blood in skin tissue using digital RGB images. The total blood concentration and oxygen saturation can also be reconstructed. Monte Carlo simulation of light transport specifies a relation between the chromophore concentrations and Commission Internationale de l'Eclairage XYZ, which are compatible with the common RGB working space. Experiments with a tissuelike agar gel phantom demonstrated the possibility of the method. In vivo imaging of a human hand during forearm occlusion demonstrated the ability of the method to evaluate hemodynamics of skin tissue.

In vivo determination of skin near-infrared optical properties using diffuse optical spectroscopy

We develop a superficial diffusing probe with a 3 mm source-detector separation that can be used in combination with diffuse optical spectroscopic (DOS) methods to noninvasively determine full-spectrum optical properties of superficial in vivo skin in the wavelength range from 650 to 1000 nm. This new probe uses a highly scattering layer to diffuse photons emitted from a collimated light source and relies on a two-layer diffusion model to determine tissue absorption coefficient mu a and reduced scattering coefficient mu's. By employing the probe to measure two-layer phantoms that mimic the optical properties of skin, we demonstrate that the probe has an interrogation depth of 1 to 2 mm. We carry out SSFDPM (steady state frequency-domain photon migration) measurements using this new probe on the volar forearm and palm of 15 subjects, including five subjects of African descent, five Asians, and five Caucasians. The optical properties of in vivo skin determined using the superficial diffusing probe show considerable similarity to published optical properties of carefully prepared ex vivo epidermis+dermis.

Visualizing depth and thickness of a local blood region in skin tissue using diffuse reflectance images

A method is proposed for visualizing the depth and thickness distribution of a local blood region in skin tissue using diffuse reflectance images at three isosbestic wavelengths of hemoglobin: 420, 585, and 800 nm. Monte Carlo simulation of light transport specifies a relation among optical densities, depth, and thickness of the region under given concentrations of melanin in epidermis and blood in dermis. Experiments with tissue-like agar gel phantoms indicate that a simple circular blood region embedded in scattering media can be visualized with errors of 6% for the depth and 22% for the thickness to the given values. In-vivo measurements on human veins demonstrate that results from the proposed method agree within errors of 30 and 19% for the depth and thickness, respectively, with values obtained from the same veins by the conventional ultrasound technique. Numerical investigation with the Monte Carlo simulation of light transport in the skin tissue is also performed to discuss effects of deviation in scattering coefficients of skin tissue and absorption coefficients of the local blood region from the typical values of the results. The depth of the local blood region is over- or underestimated as the scattering coefficients of epidermis and dermis decrease or increase, respectively, while the thickness of the region agrees well with the given values below 1.2 mm. Decreases or increases of hematocrit value give over- or underestimation of the thickness, but they have almost no influence on the depth.

Depth visualization of a local blood region in skin tissue by use of diffuse reflectance images

A simple method is proposed for visualizing the depth distribution of a local blood region in skin tissue by using diffuse reflectance images at two isosbestic wavelengths of hemoglobin, 420 and 585 nm. Monte Carlo simulation of light transport specifies a relation between optical densities and the depth of the region under given concentrations of melanin in the epidermis and blood in the dermis. Phantom and in vivo experiments were performed to show the usefulness of the method.