Absorption and reduced scattering coefficients in epidermis and dermis from a Swedish cohort study.

Significance: Knowledge of optical properties is important to accurately model light propagation in tissue, but in vivo reference data are sparse. Aim: The aim of our study was to present in vivo skin optical properties from a large Swedish cohort including 3809 subjects using a three-layered skin model and spatially resolved diffuse reflectance spectroscopy (Periflux PF6000 EPOS). Approach: Diffuse reflectance spectra (475 to 850 nm) at 0.4 and 1.2 mm source-detector separations were analyzed using an inverse Monte Carlo method. The model had one epidermis layer with variable thicknesses and melanin-related absorptions and two dermis layers with varying hemoglobin concentrations and equal oxygen saturations. The reduced scattering coefficient was equal across all layers. Results: Median absorption coefficients (mm-1) in the upper dermis ranged from 0.094 at 475 nm to 0.0048 at 850 nm and similarly in the lower dermis from 0.059 to 0.0035. The reduced scattering coefficient (mm-1) ranged from 3.22 to 1.20, and the sampling depth (mm) ranged from 0.23 to 0.38 (0.4 mm separation) and from 0.49 to 0.68 (1.2 mm separation). There were differences in optical properties across sex, age groups, and BMI categories. Conclusions: Reference material for skin optical properties is presented.

Workload and sex effects in comprehensive assessment of cutaneous microcirculation.

INTRODUCTION: Workload and sex-related differences have been proposed as factors of importance when evaluating the microcirculation. Simultaneous assessments with diffuse reflectance spectroscopy (DRS) and laser Doppler flowmetry (LDF) enable a comprehensive evaluation of the microcirculation. The aim of the study was to compare the response between sexes in the microcirculatory parameters red blood cell (RBC) tissue fraction, RBC oxygen saturation, average vessel diameter, and speed-resolved perfusion during baseline, cycling, and recovery, respectively. METHODS: In 24 healthy participants (aged 20 to 30 years, 12 females), cutaneous microcirculation was assessed by LDF and DRS at baseline, during a workload generated by cycling at 75 to 80 % of maximal age-predicted heart rate, and recovery, respectively. RESULTS: Females had significantly lower RBC tissue fraction and total perfusion in forearm skin microcirculation at all phases (baseline, workload, and recovery). All microvascular parameters increased significantly during cycling, most evident in RBC oxygen saturation (34 % increase on average) and perfusion (9-fold increase in total perfusion). For perfusion, the highest speeds (>10 mm/s) increased by a factor of 31, whereas the lowest speeds (<1 mm/s) increased by a factor of 2. CONCLUSION: Compared to a resting state, all studied microcirculation measures increased during cycling. For perfusion, this was mainly due to increased speed, and only to a minor extent due to increased RBC tissue fraction. Skin microcirculatory differences between sexes were seen in RBC concentration and total perfusion.

In vivo dose-response analysis to acetylcholine: pharmacodynamic assessment by polarized reflectance spectroscopy.

Transdermal iontophoresis offers an in vivo alternative to the strain-gauge model for measurement of vascular function but is limited due to lack of technical solutions for outcome assessment. The aims of this study were to, after measurement by polarized reflectance spectroscopy (PRS), use pharmacodynamic dose-response analysis on responses to different concentrations of acetylcholine (ACh); and to examine the effect of three consecutively administered iontophoretic current pulses. The vascular responses in 15 healthy volunteers to iontophorised ACh (5 concentrations, range 0.0001% to 1%, three consecutive pulses of 0.02 mA for 10 min each) were recorded using PRS. Data were fitted to a four-parameter logistic dose response model and compared. Vascular responses were quantifiable by PRS. Similar pharmacodynamic dose response curves could be generated irrespectively of the ACh concentration. Linearly increasing maximum vasodilatory responses were registered with increasing concentration of ACh. A limited linear dose effect of the concentration of ACh was seen between pulses. Polarized reflectance spectroscopy is well suited for measuring vascular responses to iontophoretically administrated ACh. The results of this study support further development of iontophoresis as a method to study vascular function and pharmacological responses in vivo.

Machine learning for direct oxygen saturation and hemoglobin concentration assessment using diffuse reflectance spectroscopy.

SIGNIFICANCE: Diffuse reflectance spectroscopy (DRS) is frequently used to assess oxygen saturation and hemoglobin concentration in living tissue. Methods solving the inverse problem may include time-consuming nonlinear optimization or artificial neural networks (ANN) determining the absorption coefficient one wavelength at a time. AIM: To present an ANN-based method that directly outputs the oxygen saturation and the hemoglobin concentration using the shape of the measured spectra as input. APPROACH: A probe-based DRS setup with dual source-detector separations in the visible wavelength range was used. ANNs were trained on spectra generated from a three-layer tissue model with oxygen saturation and hemoglobin concentration as target. RESULTS: Modeled evaluation data with realistic measurement noise showed an absolute root-mean-square (RMS) deviation of 5.1% units for oxygen saturation estimation. The relative RMS deviation for hemoglobin concentration was 13%. This accuracy is at least twice as good as our previous nonlinear optimization method. On blood-intralipid phantoms, the RMS deviation from the oxygen saturation derived from partial oxygen pressure measurements was 5.3% and 1.6% in two separate measurement series. Results during brachial occlusion showed expected patterns. CONCLUSIONS: The presented method, directly assessing oxygen saturation and hemoglobin concentration, is fast, accurate, and robust to noise.

Handheld multispectral imager for quantitative skin assessment in low-resource settings.

SIGNIFICANCE: Spatial frequency domain imaging (SFDI) is a quantitative imaging method to measure absorption and scattering of tissue, from which several chromophore concentrations (e.g., oxy-/deoxy-/meth-hemoglobin, melanin, and carotenoids) can be calculated. Employing a method to extract additional spectral bands from RGB components (that we named cross-channels), we designed a handheld SFDI device to account for these pigments, using low-cost, consumer-grade components for its implementation and characterization. AIM: With only three broad spectral bands (red, green, blue, or RGB), consumer-grade devices are often too limited. We present a methodology to increase the number of spectral bands in SFDI devices that use RGB components without hardware modification. APPROACH: We developed a compact low-cost RGB spectral imager using a color CMOS camera and LED-based mini projector. The components' spectral properties were characterized and additional cross-channel bands were calculated. An alternative characterization procedure was also developed that makes use of low-cost equipment, and its results were compared. The device performance was evaluated by measurements on tissue-simulating optical phantoms and in-vivo tissue. The measurements were compared with another quantitative spectroscopy method: spatial frequency domain spectroscopy (SFDS). RESULTS: Out of six possible cross-channel bands, two were evaluated to be suitable for our application and were fully characterized (520 ± 20 nm; 556 ± 18 nm). The other four cross-channels presented a too low signal-to-noise ratio for this implementation. In estimating the optical properties of optical phantoms, the SFDI data have a strong linear correlation with the SFDS data (R2 = 0.987, RMSE = 0.006 for µa, R2 = 0.994, RMSE = 0.078 for µs'). CONCLUSIONS: We extracted two additional spectral bands from a commercial RGB system at no cost. There was good agreement between our device and the research-grade SFDS system. The alternative characterization procedure we have presented allowed us to measure the spectral features of the system with an accuracy comparable to standard laboratory equipment.

Assessment of oxygenation with polarized light spectroscopy enables new means for detecting vascular events in the skin.

INTRODUCTION: Impaired oxygenation in the skin may occur in disease states and after reconstructive surgery. We used tissue viability imaging (TiVi) to measure changes in oxygenation and deoxygenation of haemoglobin in an in vitro model and in the dermal microcirculation of healthy individuals. MATERIALS AND METHODS: Oxygenation was measured in human whole blood with different levels of oxygenation. In healthy subjects, changes in red blood cell concentration (CRBC,TiVi), oxygenation (ΔCOH,TiVi) and deoxygenation (ΔCDOH,TiVi) of haemoglobin were measured during and after arterial and venous occlusion using TiVi and were compared with measurements from the enhanced perfusion and oxygen saturation system (EPOS). RESULTS: During arterial occlusion, CRBC,TiVi remained unchanged while ΔCOH,TiVi decreased to -44.2 (10.4) AU (p = 0.04), as compared to baseline. After release, CRBC,TiVi increased to 39.2 (18.8) AU (p < 0.001), ΔCOH,TiVi increased to 38.5. During venous occlusion, CRBC,TiVi increased to 28.9 (11.2) AU (p < 0.001), ΔCOH,TiVi decreased to -52.2 (46.1) AU (p < 0.001) compared to baseline after 5 min of venous occlusion. There was a significant correlation between the TiVi Oxygen Mapper and EPOS, for arterial (r = 0.92, p < 0.001) and venous occlusion (r = 0.87, p < 0.001), respectively. CONCLUSION: This study shows that TiVi can measure trends in oxygenation and deoxygenation of haemoglobin during arterial and venous stasis in healthy individuals.

Validation of speed-resolved laser Doppler perfusion in a multimodal optical system using a blood-flow phantom.

The PeriFlux 6000 EPOS system combines diffuse reflectance spectroscopy (DRS) and laser Doppler flowmetry (LDF) for the assessment of oxygen saturation (expressed in percentage), red blood cell (RBC) tissue fraction (expressed as volume fraction, %RBC), and perfusion (%RBC × mm / s) in the microcirculation. It also allows the possibility of separating the perfusion into three speed regions (0 to 1, 1 to 10, and >10 mm / s). We evaluate the speed-resolved perfusion components, i.e., the relative amount of perfusion within each speed region, using a blood-flow phantom. Human blood was pumped through microtubes with an inner diameter of 0.15 mm. Measured DRS and LDF spectra were compared to Monte Carlo-simulated spectra in an optimization routine, giving the best-fit parameters describing the measured spectra. The root-mean-square error for each of the three speed components (0 to 1, 1 to 10, and >10 mm / s, respectively) when describing the blood-flow speed in the microtubes was 2.9%, 8.1%, and 7.7%. The presented results show that the system can accurately discriminate blood perfusion originating from different blood-flow speeds, which may enable improved measurement of healthy and dysfunctional microcirculatory flow.

Machine learning in multiexposure laser speckle contrast imaging can replace conventional laser Doppler flowmetry.

Laser speckle contrast imaging (LSCI) enables video rate imaging of blood flow. However, its relation to tissue blood perfusion is nonlinear and depends strongly on exposure time. By contrast, the perfusion estimate from the slower laser Doppler flowmetry (LDF) technique has a relationship to blood perfusion that is better understood. Multiexposure LSCI (MELSCI) enables a perfusion estimate closer to the actual perfusion than that using a single exposure time. We present and evaluate a method that utilizes contrasts from seven exposure times between 1 and 64 ms to calculate a perfusion estimate that resembles the perfusion estimate from LDF. The method is based on artificial neural networks (ANN) for fast and accurate processing of MELSCI contrasts to perfusion. The networks are trained using modeling of Doppler histograms and speckle contrasts from tissue models. The importance of accounting for noise is demonstrated. Results show that by using ANN, MELSCI data can be processed to LDF perfusion with high accuracy, with a correlation coefficient R = 1.000 for noise-free data, R = 0.993 when a moderate degree of noise is present, and R = 0.995 for in vivo data from an occlusion-release experiment.

In vivo characterization of light scattering properties of human skin in the 475- to 850-nm wavelength range in a Swedish cohort.

We have determined in vivo optical scattering properties of normal human skin in 1734 subjects, mostly with fair skin type, within the Swedish CArdioPulmonary bioImage Study. The measurements were performed with a noninvasive system, integrating spatially resolved diffuse reflectance spectroscopy and laser Doppler flowmetry. Data were analyzed with an inverse Monte Carlo algorithm, accounting for both scattering, geometrical, and absorbing properties of the tissue. The reduced scattering coefficient was found to decrease from 3.16 ± 0.72 to 1.13 ± 0.27 mm-1 (mean ± SD) in the 475- to 850-nm wavelength range. There was a negative correlation between the reduced scattering coefficient and age, and a significant difference between men and women in the reduced scattering coefficient as well as in the fraction of small scattering particles. This large study on tissue scattering with mean values and normal variation can serve as a reference when designing diagnostic techniques or when evaluating the effect of therapeutic optical systems.

Evaluation of a pointwise microcirculation assessment method using liquid and multilayered tissue simulating phantoms.

A fiber-optic probe-based instrument, designed for assessment of parameters related to microcirculation, red blood cell tissue fraction (fRBC), oxygen saturation (SO2), and speed resolved perfusion, has been evaluated using state-of-the-art tissue phantoms. The probe integrates diffuse reflectance spectroscopy (DRS) at two source-detector separations and laser Doppler flowmetry, using an inverse Monte Carlo method for identifying the parameters of a multilayered tissue model. Here, we characterize the accuracy of the DRS aspect of the instrument using (1) liquid blood phantoms containing yeast and (2) epidermis-dermis mimicking solid-layered phantoms fabricated from polydimethylsiloxane, titanium oxide, hemoglobin, and coffee. The root-mean-square (RMS) deviations for fRBC for the two liquid phantoms were 11% and 5.3%, respectively, and 11% for the solid phantoms with highest hemoglobin signatures. The RMS deviation for SO2 was 5.2% and 2.9%, respectively, for the liquid phantoms, and 2.9% for the solid phantoms. RMS deviation for the reduced scattering coefficient (µs'), for the solid phantoms was 15% (475 to 850 nm). For the liquid phantoms, the RMS deviation in average vessel diameter (D) was 1 µm. In conclusion, the skin microcirculation parameters fRBC and SO2, as well as, µs' and D are estimated with reasonable accuracy.

Temporal and spatiotemporal variability in comprehensive forearm skin microcirculation assessment during occlusion protocols.

Forearm skin hyperemia during release after brachial occlusion has been proposed for evaluating peripheral arterial disease and endothelial dysfunction. We used a novel fiberoptic system integrating Laser Doppler Flowmetry and Diffuse Reflectance Spectroscopy for a comprehensive pointwise model based microcirculation characterization. The aim was to evaluate and compare the temporal and the spatiotemporal variabilities in forearm skin microcirculation parameters (speed resolved perfusion; low speed <1mm/s, PerfSR, <1; mid-speed 1-10mm/s, high speed >10mm/s, and total perfusion (PerfSR, tot); the concentration and oxygenation of red blood cells, CRBC and SO2). Ten healthy subjects underwent arterial and venous forearm occlusions (AO, VO), repeated within one week. The repeatability was calculated as the coefficient of variation (CV) and the agreement as the intra-class correlation coefficient (ICC). The temporal CVs for conventional perfusion, Perfconv, PerfSR, tot, CRBC and SO2 were 14%, 12%, 9% and 9%, respectively, while the ICC were >0.75 (excellent). The perfusion measures generally had a higher spatiotemporal than temporal variability, which was not the case for SO2 and CRBC. The corresponding spatiotemporal CVs were 33%, 32%, 18% and 15%, respectively. During VO, CRBC had a CV<35% and ICC>0.40 (fair-good), and after release this was the case for CRBC (AO and VO), SO2 (VO) and PerfSR, <1 (VO). In conclusion, the skin microcirculation parameters showed excellent temporal repeatability, while the spatiotemporal repeatability especially for perfusion was poorer. The parameters with acceptable repeatability and fair-good agreement were: CRBC during and after release of VO, the PerfSR, <1 after release of VO, the SO2 and the CRBC after release of AO. However, the value of these parameters in discriminating endothelial function remains to be studied.

Assessment of advanced glycated end product accumulation in skin using auto fluorescence multispectral imaging.

Several studies have shown that advanced glycation end products (AGE) play a role in both the microvascular and macrovascular complications of diabetes and are closely linked to inflammation and atherosclerosis. AGEs accumulate in skin and can be detected using their auto fluorescence (AF). A significant correlation exists between AGE AF and the levels of AGEs as obtained from skin biopsies. A commercial device, the AGE Reader, has become available to assess skin AF for clinical purposes but, while displaying promising results, it is limited to single-point measurements performed in contact to skin tissue. Furthermore, in vivo imaging of AGE accumulation is virtually unexplored. We proposed a non-invasive, contact-less novel technique for quantifying fluorescent AGE deposits in skin tissue using a multispectral imaging camera setup (MSI) during ultraviolet (UV) exposure. Imaging involved applying a region-of-interest mask, avoiding specular reflections and a simple calibration. Results of a study conducted on 16 subjects with skin types ranging from fair to deeply pigmented skin, showed that AGE measured with MSI in forearm skin was significantly correlated with the AGE reference method (AGE Reader on forearm skin, R=0.68, p=0.005). AGE measured in facial skin was borderline significantly related to AGE Reader on forearm skin (R=0.47, p=0.078). These results support the use of the technique in devices for non-touch measurement of AGE content in either facial or forearm skin tissue over time.

Oxygen saturation, red blood cell tissue fraction and speed resolved perfusion - A new optical method for microcirculatory assessment.

We have developed a new fiber-optic system that combines diffuse reflectance spectroscopy (DRS) and laser Doppler Flowmetry (LDF) for a multi-modal assessment of the microcirculation. Quantitative data is achieved with an inverse Monte Carlo algorithm based on an individually adaptive skin model. The output parameters are calculated from the model and given in absolute units: hemoglobin oxygen saturation (%), red blood cell (RBC) tissue fraction (%), and the speed resolved RBC perfusion separated into three speed regions; 0-1mm/s, 1-10mm/s and above 10mm/s (% mm/s). The aim was to explore microcirculatory parameters using the new optical method, integrating DRS and LDF in a joint skin model, during local heating of the dorsal foot and venous and arterial occlusion of the forearm in 23 healthy subjects (age 20-28years). There were differences in the three speed regions in regard to blood flow changes due to local heating, where perfusion for high speeds increased the most. There was also a high correlation between changes in oxygenation and changes in perfusion for higher speeds. Oxygen saturation at baseline was 44% on foot, increasing to 83% at plateau after heating. The larger increase in perfusion for higher speeds than for lower speeds together with the oxygenation increase during thermal provocation, shows a local thermoregulatory blood flow in presumably arteriolar dermal vessels. In conclusion, there are improved possibilities to assess microcirculation using integrated DRS and LDF in a joint skin model by enabling both oxygenation and speed resolved blood flow assessment simultaneously and in the same skin site. Output parameters in absolute units may also yield new insights about the microcirculatory system.

Microcirculation assessment using an individualized model for diffuse reflectance spectroscopy and conventional laser Doppler flowmetry.

Microvascular assessment would benefit from co-registration of blood flow and hemoglobin oxygenation dynamics during stimulus response tests. We used a fiber-optic probe for simultaneous recording of white light diffuse reflectance (DRS; 475-850 nm) and laser Doppler flowmetry (LDF; 780 nm) spectra at two source-detector distances (0.4 and 1.2 mm). An inverse Monte Carlo algorithm, based on a multiparameter three-layer adaptive skin model, was used for analyzing DRS data. LDF spectra were conventionally processed for perfusion. The system was evaluated on volar forearm recordings of 33 healthy subjects during a 5-min systolic occlusion protocol. The calibration scheme and the optimal adaptive skin model fitted DRS spectra at both distances within 10%. During occlusion, perfusion decreased within 5 s while oxygenation decreased slowly (mean time constant 61 s; dissociation of oxygen from hemoglobin). After occlusion release, perfusion and oxygenation increased within 3 s (inflow of oxygenized blood). The increased perfusion was due to increased blood tissue fraction and speed. The supranormal hemoglobin oxygenation indicates a blood flow in excess of metabolic demands. In conclusion, by integrating DRS and LDF in a fiber-optic probe, a powerful tool for assessment of blood flow and oxygenation in the same microvascular bed has been presented.

Inverse Monte Carlo in a multilayered tissue model: merging diffuse reflectance spectroscopy and laser Doppler flowmetry.

The tissue fraction of red blood cells (RBCs) and their oxygenation and speed-resolved perfusion are estimated in absolute units by combining diffuse reflectance spectroscopy (DRS) and laser Doppler flowmetry (LDF). The DRS spectra (450 to 850 nm) are assessed at two source-detector separations (0.4 and 1.2 mm), allowing for a relative calibration routine, whereas LDF spectra are assessed at 1.2 mm in the same fiber-optic probe. Data are analyzed using nonlinear optimization in an inverse Monte Carlo technique by applying an adaptive multilayered tissue model based on geometrical, scattering, and absorbing properties, as well as RBC flow-speed information. Simulations of 250 tissue-like models including up to 2000 individual blood vessels were used to evaluate the method. The absolute root mean square (RMS) deviation between estimated and true oxygenation was 4.1 percentage units, whereas the relative RMS deviations for the RBC tissue fraction and perfusion were 19% and 23%, respectively. Examples of in vivo measurements on forearm and foot during common provocations are presented. The method offers several advantages such as simultaneous quantification of RBC tissue fraction and oxygenation and perfusion from the same, predictable, sampling volume. The perfusion estimate is speed resolved, absolute (% RBC×mm/s), and more accurate due to the combination with DRS.

Inverse Monte Carlo for estimation of scattering and absorption in liquid optical phantoms.

A spectroscopic probe with multiple detecting fibers was used for quantifying absorption and scattering in liquid optical phantoms. The phantoms were mixtures of Intralipid and red and blue food dyes. Intensity calibration for the detecting fibers was undertaken using either a microsphere suspension (absolute calibration) or a uniform detector illumination (relative calibration between detectors). Two different scattering phase functions were used in an inverse Monte Carlo algorithm. Data were evaluated for residual spectra (systematic deviations and magnitude) and accuracy in estimation of scattering and absorption. Spectral fitting was improved by allowing for a 10% intensity relaxation in the optimization algorithm. For a multi-detector setup, non-systematic residual spectrum was only found using the more complex Gegenbauer-kernel phase function. However, the choice of phase function did not influence the accuracy in the estimation of absorption and scattering. Similar estimation accuracy as in the multi-detector setup was also obtained using either two relative calibrated detectors or one absolute calibrated detector at a fiber separation of 0.46 mm.

Inverse Monte Carlo method in a multilayered tissue model for diffuse reflectance spectroscopy.

Model based data analysis of diffuse reflectance spectroscopy data enables the estimation of optical and structural tissue parameters. The aim of this study was to present an inverse Monte Carlo method based on spectra from two source-detector distances (0.4 and 1.2 mm), using a multilayered tissue model. The tissue model variables include geometrical properties, light scattering properties, tissue chromophores such as melanin and hemoglobin, oxygen saturation and average vessel diameter. The method utilizes a small set of presimulated Monte Carlo data for combinations of different levels of epidermal thickness and tissue scattering. The path length distributions in the different layers are stored and the effect of the other parameters is added in the post-processing. The accuracy of the method was evaluated using Monte Carlo simulations of tissue-like models containing discrete blood vessels, evaluating blood tissue fraction and oxygenation. It was also compared to a homogeneous model. The multilayer model performed better than the homogeneous model and all tissue parameters significantly improved spectral fitting. Recorded in vivo spectra were fitted well at both distances, which we previously found was not possible with a homogeneous model. No absolute intensity calibration is needed and the algorithm is fast enough for real-time processing.

Forced detection Monte Carlo algorithms for accelerated blood vessel image simulations.

Two forced detection (FD) variance reduction Monte Carlo algorithms for image simulations of tissue-embedded objects with matched refractive index are presented. The principle of the algorithms is to force a fraction of the photon weight to the detector at each and every scattering event. The fractional weight is given by the probability for the photon to reach the detector without further interactions. Two imaging setups are applied to a tissue model including blood vessels, where the FD algorithms produce identical results as traditional brute force simulations, while being accelerated with two orders of magnitude. Extending the methods to include refraction mismatches is discussed.

Measurement depth and volume in laser Doppler flowmetry.

A new method for estimating the measurement depth and volume in laser Doppler flowmetry (LDF) is presented. The method is based on Monte Carlo simulations of light propagation in tissue. The contribution from each individual Doppler shift is calculated and thereby multiple Doppler shifts are handled correctly. Different LDF setups for both probe based (0.0, 0.25, 0.5, and 1.2 mm source-detector separation) and imaging systems (0.5 and 2.0 mm beam diameter) are considered, at the wavelengths 543 nm, 633 nm, and 780 nm. Non-linear speckle pattern effects are accounted for in the imaging system setups. The effects of tissue optical properties, blood concentration, and blood oxygen saturation are evaluated using both homogeneous tissue models and a layered skin model. The results show that the effect on the measurement depth of changing tissue properties is comparable to the effect of changing the system setup, e.g. source-detector separation and wavelength. Skin pigmentation was found to have a negligible effect on the measurement depth. Examples of measurement depths are (values are given for a probe based system with 0.25 mm source-detector separation and an imaging system with a 0.5 mm beam diameter, respectively, both operating at 780 nm): muscle - 0.55/0.79 mm; liver - 0.40/0.53 mm; gray matter - 0.48/0.68 mm; white matter - 0.20/0.20 mm; index finger pulp - 0.41/0.53 mm; forearm skin - 0.53/0.56 mm; heat provoked forearm skin - 0.66/0.67 mm.

Optical microcirculatory skin model: assessed by Monte Carlo simulations paired with in vivo laser Doppler flowmetry.

An optical microvascular skin model, valid at 780 nm, was developed. The model consisted of six layers with individual optical properties and variable thicknesses and blood concentrations at three different blood flow velocities. Monte Carlo simulations were used to evaluate the impact of various model parameters on the traditional laser Doppler flowmetry (LDF) measures. A set of reference Doppler power spectra was generated by simulating 7000 configurations, varying the thickness and blood concentrations. Simulated spectra, at two different source detector separations, were compared with in vivo recorded spectra, using a nonlinear search algorithm for minimizing the deviation between simulated and measured spectra. The model was validated by inspecting the thickness and blood concentrations that generated the best fit. These four parameters followed a priori expectations for the measurement situations, and the simulated spectra agreed well with the measured spectra for both detector separations. Average estimated dermal blood concentration was 0.08% at rest and 0.63% during heat provocation (44 degrees C) on the volar side of the forearm and 1.2% at rest on the finger pulp. The model is crucial for developing a technique for velocity-resolved absolute LDF measurements with known sampling volume and can also be useful for other bio-optical modalities.