Fast and precise image generation of blood vessels embedded in skin

Christian Zoller; Alwin Kienle

doi:10.1117/1.JBO.24.1.015002

29 January 2019 Fast and precise image generation of blood vessels embedded in skin

Christian Zoller, Alwin Kienle

Author Affiliations +

Journal of Biomedical Optics, Vol. 24, Issue 1, 015002 (January 2019). https://doi.org/10.1117/1.JBO.24.1.015002

Abstract

A software for fast rendering the visual appearance of a blood vessel located in human skin was developed based on a numerical solution of the radiative transfer equation. The user can specify geometrical properties, such as the depth and the diameter of the vessel, and physiological properties, such as the oxygen saturation of the vessel or the blood concentration in the skin. From these data, the spatially and spectrally resolved reflectance from the skin containing the blood vessel is calculated via Monte Carlo simulations, by which a two-dimensional image is generated. The short calculation time of about a second is achieved by precalculating and storing the spatially resolved reflectance for a variety of combinations of the optical and geometrical properties. This concept gives the user the opportunity to rapidly explore the influence of the physiological and geometrical properties of the investigated blood vessel on its visual appearance. The correctness of the lookup table was validated by comparison with independent Monte Carlo simulations. Rendering examples of different blood vessels in human skins are given. The current version of the software can be downloaded at https://www.ilm-ulm.de/software.

1. Introduction

For a long time, color phenomena in nature, such as the changing colors of the sky from sunset to twilight, have been fascinating and inspiring people. Knowing the basic physics of the involved effects, such as the single or multiple scattering of light by particles or inhomogeneities, allows us to understand these phenomena.¹ About 20 years ago, the physical explanation for the bluish appearance of veins, besides the fact that blood itself is red, was quantitatively given.² It was found that the spectral characteristics of light propagation in combination with Land’s retinex theory can be used to explain the color perception of blood vessels. Interestingly, the reason for the bluish color of veins is not a larger reflection of blue light compared to red light from the skin above the vein. It rather is the larger decrease of the reflection for the red light than for the blue light from the skin above the vessel compared to the surrounding skin. Red light has a larger penetration depth compared to blue light in the involved tissues and is, thus, affected to a larger extent by a sufficiently deep vessel. Furthermore, it was calculated that not only veins but also arteries can have a bluish appearance for certain conditions.² In general, however, arteries look less bluish compared to veins having the same radius and depth due to the decreased absorption of oxygenated hemoglobin compared to deoxygenated hemoglobin in the red wavelength range. Recently, it was postulated that Rayleigh scattering in the skin has a major impact on the bluish appearance of blood vessels.³ Studies of other authors focused on related topics, for example, the multiple influences on the color appearance of skin, such as melanin concentration as well as different melanin types, the oiliness of skin or the blood concentration in skin.⁴^,⁵ Besides the interests to render skin correctly, diagnostics methods, such as the evaluation of erythema for different skin types, efficient port wine stain treatment, or forensic applications, where the colors of bruises are studied, also benefit from the understanding of light propagation in human skin.⁶^–⁹

Researchers have applied different techniques to calculate the light propagation in biological tissue. One of the most commonly used methods, especially in the field of computer graphics, is the so called diffusion approximation.¹⁰^–¹² The diffusion approximation has its benefits in the easy implementation and high computational speed, but it lacks in accuracy for certain circumstances, for example, if low scattering and high absorption is present. If a higher precision is needed, the radiative transfer equation (RTE) has to be solved exactly. Analytical solutions have been found for certain geometries, such as layered media with and without fluorescence.¹³ Unfortunately, the RTE cannot be solved analytically yet for more complex geometries, such as a vessel embedded in a layered skin model. The Monte Carlo (MC) simulation as an, in principle, exact numerical solution of the RTE has been established as the gold standard to calculate the light propagation for such problems.¹⁴^,¹⁵ A big disadvantage of the MC method is the potentially enormous computational time, due to its statistical nature. To reduce this handicap, massive parallelization using graphic cards (GPU), as well as multiple variance reduction methods, have been developed.¹⁶^–¹⁹

In this study, a software is presented to render the visual appearance of a blood vessel embedded in skin. The skin is modeled by two layers representing the epidermis and the dermis. A MATLAB (Mathworks Inc., Natick, Massachusetts) interface is provided, where multiple physiological and geometrical parameters, such as the oxygen saturation of the vessel or the radius of the vessel, can be specified by the user. The image generation is based on the spectrally and spatially resolved reflectance above the vessel, within the visible range of light, which is converted to an RGB-image using the CIE-XYZ color space.²⁰ The presented software can calculate the spatially and spectrally resolved reflectance, and hence, RGB-image in roughly one second using a lookup table (LUT). To accelerate the generation of the LUT, a variance reduction method was implemented. In the second part of this study, exemplary images are generated by the proposed method. Hereby, the results found by Kienle et al.² could be confirmed, such as the fact that a small vessel close to the surface will appear red but can look bluish by increasing its diameter, or the finding that veins are more likely to appear bluish compared to arteries. Beyond that, further interesting cases are investigated. First, the influence of the radius and depth for veins and arteries is systematically investigated. Second, the appearance of superficial vessels depending on the radius and the oxygen saturation of the blood is discussed. Finally, the influence of different spectral dependencies of the reduced scattering coefficient are presented.

2. Method

The visual appearance of a blood vessel depends on multiple properties, such as its depth, its radius, or the absorption and scattering coefficients of the tissue close to the vessel and of the vessel itself. To predict the appearance of a blood vessel and the surrounding skin, the spatially and spectrally resolved reflectance for the skin surface near the vessel is needed in the visible wavelength range between 400 and 800 nm. In this section, the generation of an efficient and accurate LUT, delivering the spectrally and spatially resolved reflectance $I_{srr} (λ, i)$ detected at the skin surface along a vessel, is described. A variance reduction method was applied to accelerate the generation of the LUT. Curve fitting was performed to reduce the amount of data stored in the LUT. A MATLAB user interface is provided, where the operator can set multiple geometrical and physiological properties. These properties are converted to an image in sRGB color space using the generated LUT. Thus, the appearance of the blood vessel can be rendered in a two-dimensional image.

2.1.

LUT Generation

A graphic card-based Monte Carlo (MC) simulation, as presented previously,²¹ was used to generate an LUT, storing the spatially resolved reflectance $I_{srr} (i)$ for different scenarios of a three-dimensional model as shown in Fig. 1. Note that the model is symmetric along its $y$ axis. The epidermal layer was assumed to be a slab, the dermis to be semi-infinitely thick, and the blood vessel to be extended to infinity along the $y$ axis. The blood vessel’s $x$ -position was set to be $x = 0$ . Photons were started homogeneously on a line from $x = - 0.5 l_{s}$ to $x = 0.5 l_{s}$ at $y = 0$ in $z$ direction where $l_{s}$ is the light source size. Reflected photons were recorded with respect to their $x$ -position within a certain range. The light source size $l_{s}$ was set three times as long as the detection range to approach an infinitely extended light source, such that problems at the detector-boundaries were avoided. For the following simulations, $l_{s}$ was fixed at $l_{s} = 45 mm$ and the detection bin size was kept constant at $| x_{i + 1} - x_{i} | = 0.1 mm$ . The spatially resolved reflectance

Eq. (1)

I_{srr} (i) = \frac{P (i)}{P} 3 N

can be derived, where

P (i)

is the number of photons received between

x_{i}

and

x_{i + 1}

,

P

is the total number of photons started, and

N

is the number of spatial bins. The factor 3 is due to the three times larger light source size compared with the detection range.

Fig. 1

Simulation setup of a blood vessel embedded in a two-layered model, representing the epidermis (brown) and the dermis (light red). The blood vessel was always placed in the dermis and extends to infinity in $\pm y$ direction. The epidermis and dermis extend to infinity in $\pm x$ and $\pm y$ directions. Photons were started uniformly on the surface along the dotted line, perpendicular to the blood vessel. Reflected photons were detected with respect to their $x$ -position.

Even though a variance reduction method was used, as described in Sec. 2.2, only a limited number of dimensions and corresponding grid points could be chosen for the LUT, to keep the computational time as well as the memory demand of the final LUT in an acceptable range. The LUT has six-dimensions representing properties that have a major influence on the color impression:

1. Absorption coefficient of the epidermis $μ_{a, 1}$ in ${mm}^{- 1}$ .
2. Absorption coefficient of the dermis $μ_{a, 2}$ in ${mm}^{- 1}$ .
3. Absorption coefficient of the blood vessel $μ_{a, 3}$ in ${mm}^{- 1}$ .
4. Reduced scattering coefficient of the dermis $μ_{s}^{'}$ in ${mm}^{- 1}$ . The same reduced scattering coefficient was used for the blood vessel. A two times higher reduced scattering coefficient was applied for the epidermis.⁴^,²²
5. Depth of the blood vessel $z$ in mm. The depth is defined as the shortest distance of the blood vessel (top of the vessel) to the surface. The minimal depth is limited by the thickness of the epidermis since the vessel is always placed in the dermis.
6. Radius of the blood vessel $r$ in mm.

The used ranges of the dimensions as well as the corresponding number of grid points are listed in Table 1. The grid points for $z$ and $r$ were spaced linearly across the corresponding ranges. The grid points for $μ_{s}^{'}$ were also spaced linearly, but a finer linear spacing for $μ_{s}^{'} < 1 {mm}^{- 1}$ than $μ_{s}^{'} > 1 {mm}^{- 1}$ was chosen. Values for the absorption coefficients $μ_{a, m}$ , $m = 1,2, 3$ were spaced logarithmically such that the difference between two grid points is smaller for small $μ_{a, m}$ and is increasing for larger $μ_{a, m}$ .

Table 1

Value ranges and the corresponding number of grid points for each dimension of the LUT.

Variable	Range	# grid points
$μ_{a, 1}$	$10^{- 4} - 20 {mm}^{- 1}$	16
$μ_{a, 2}$	$10^{- 4} - 20 {mm}^{- 1}$	20
$μ_{a, 3}$	$10^{- 2} - 525 {mm}^{- 1}$	18
$μ_{s}^{'}$	$0.1 - 16 {mm}^{- 1}$	19
$z$	0.1 to 2 mm	15
$r$	0.02 to 1 mm	10

The values for the fixed parameters can be found in Table 2.

Table 2

Value ranges and the corresponding number of grid points for each dimension of the LUT.

Variable	Value
Refractive index $n$ (all tissues)	1.4
Anisotropy factor $g$ (all tissues)	0.9
Thickness of epidermis²³	0.1 mm
Number of photons $P$	$10^{7}$

2.2.

Scaled Monte Carlo

A variance reduction method, the so called scaled MC simulation, was used to reduce the computational time when the LUT is created.¹⁶^,¹⁹ Implementing this method allows us to generate $I_{srr} (i)$ for arbitrary combinations of $μ_{a, m}$ with a single MC simulation for given $z$ , $r$ , and $μ_{s}^{'}$ . At first, an MC simulation with rather small absorptions for each tissue type, namely ${\tilde{μ}}_{a, 1} = 10^{- 4} {mm}^{- 1}$ , ${\tilde{μ}}_{a, 2} = 10^{- 4} {mm}^{- 1}$ , and ${\tilde{μ}}_{a, 3} = 10^{- 2} {mm}^{- 1}$ , was performed. During the simulation the path lengths $s_{m}$ a photon spent in the corresponding tissues $m = 1,2, 3$ were saved for each photon. Then

Eq. (2)

I_{srr} (i) = [\sum \prod_{m = 1}^{3} \exp (s_{m} ({\tilde{μ}}_{a, m} - μ_{a, m}))] \frac{1}{P} 3 N

was calculated consecutively for all possible

16 \times 20 \times 18

combinations of

μ_{a, m}

.¹⁶^,¹⁹

2.3.

Curve Fitting

Saving the whole spatially resolved reflectance with the reported resolution would result in roughly 10 gigabytes of memory for the LUT in single precision. To reduce the memory demand as well as to increase the interpolation speed, the function

Eq. (3)

f (i) = a + \sum_{k = 1}^{K} - b (k) \exp (- \frac{x_{i}^{2}}{2 c^{2} (k)})

Eq. (4)

s.t. 0.0 \leq a \leq 1.0

Eq. (5)

b (k) \geq 0.0 \forall k if μ_{a, 3} > μ_{a, 2}

Eq. (6)

b (k) \leq 1.0 \forall k

Eq. (7)

c (k) \geq 0.0 \forall k

was fitted on the spatially resolved reflectance

I_{srr} (i)

, generated by the scaled MC simulation [see Eq. (2)]. The constraint 5 can be applied since a dip of

I_{srr} (i)

at

x = 0 mm

is expected, if the absorption of the blood vessel

μ_{a, 3}

is larger than the absorption of the dermis

μ_{a, 2}

. In total

2 K + 1

individual parameters were stored within the LUT for one spatially resolved reflectance. Curve fitting was performed in C++ using the library ceres.²⁴ It was found that

K = 3

results in a good trade-off between accuracy and memory demand. In Fig. 2, the fitting results for different

K

values are exemplary shown. To guarantee a good interpolation result, the ordering of

\vec{c}

and the corresponding

\vec{b}

within the LUT were chosen such that

c [k] \leq c [k + 1]

\forall k

. Nonunique solutions are possible, if

Eq. (8)

c (l) \approx c (k), l \neq k

since all corresponding

b (l)

and

b (k)

fulfilling

Eq. (9)

b (l) + b (k) \approx const .

lead to the same, correct, fit result. To improve the robustness of the interpolation, a postprocessing step was performed if

c (l) \approx c (k)

where

Eq. (10)

c (l) = c (k) = 0.5 (c (l) + c (k))

Eq. (11)

b (l) = b (k) = 0.5 (b (l) + b (k)) .

Fig. 2

Influence of $K$ on the fit result of the function $f$ on $I_{srr} (i)$ calculated by an MC simulation, for a certain set of $z$ , $r$ , $μ_{s}^{'}$ , and $μ_{a, m}$ . The insets show close-up views at the indicated positions.

2.4.

RGB-Image Generation

Based on the wavelength-dependent optical properties $μ_{s}^{'} (λ)$ and $μ_{a, m} (λ)$ , the visual appearance, in terms of RGB-values, of a blood vessel and the surrounding tissue can be displayed on a screen. At first a spectrally and spatially resolved reflectance $I_{srr} (λ, i)$ is delivered by the LUT. Hereby linear, $n$ -dimensional interpolation is applied separately for all $2 K + 1$ parameters for each wavelength by the MATLAB function interpn (Mathworks Inc., Natick, Massachusetts). Based on $I_{srr} (λ, i)$ , the CIE XYZ-tristimulus is calculated²⁰^,²⁵ with

Eq. (12)

X (i) = \frac{1}{W} \sum_{λ = 400 nm}^{800 nm} S (λ) I_{srr} (λ, i) \bar{x} (λ)

Eq. (13)

Y (i) = \frac{1}{W} \sum_{λ = 400 nm}^{800 nm} S (λ) I_{srr} (λ, i) \bar{y} (λ)

Eq. (14)

Z (i) = \frac{1}{W} \sum_{λ = 400 nm}^{800 nm} S (λ) I_{srr} (λ, i) \bar{z} (λ)

Eq. (15)

W = \sum_{λ = 400 nm}^{800 nm} S (λ) \bar{y} (λ),

where

S (λ)

is the spectrum of the incident light and

\bar{x} (λ)

,

\bar{y} (λ)

, and

\bar{z} (λ)

are the CIE standard observer color matching functions.²⁰^,²⁵ Throughout this study 5-nm steps were used for the calculations. RGB-values

Eq. (16)

RGB (i) = M \cdot X Y Z (i)

with

M = (\begin{matrix} 3.24 & - 1.54 & - 0.50 \\ - 0.97 & 1.88 & 0.04 \\ 0.06 & - 0.20 & 1.06 \end{matrix})

are obtained by a matrix multiplication, where the matrix M depends on the spectrum of the incident light

S (λ)

and the desired RGB-color space. Here, a CIE-D65 light source as well as the sRGB color space were chosen.²⁶ A gamma correction was applied before displaying the image on the screen.²⁷^,²⁶

2.5.

User Interface

Typically, the six-dimensions of the LUT are rather unintuitive to discuss the visual appearance of a blood vessel. It is more convenient to use properties that are directly linked to the physiology of skin, such as melanin concentration of the epidermis or oxygen saturation of the blood vessel. Therefore, a user interface was designed, where such physiological properties can be specified.

Similar to the existing literature, some assumptions about the chromophores of the skin model were made.⁴^,²²^,²⁸^–³⁰ It was assumed that the main absorber present in the epidermis is melanin, the main absorber inside the dermis is hemoglobin, which also is the only absorber for the blood vessel. The absorptions of the dermis and epidermis were extended by a baseline-absorption.

The following eight parameters $p_{1} - p_{8}$ can be entered by the user in the current version of the MATLAB interface.

• $p_{1}$ specifies the magnitude of the reduced scattering coefficient of the epidermis:⁴^,³⁰
Eq. (17)
$μ_{s}^{'} (λ) = p_{1} (14.74 λ^{- 0.22} + 2.2 \times 10^{11} λ^{- 4}) {mm}^{- 1}$

The reduced scattering coefficient of the dermis is set to be always half of $μ_{s}^{'}$ of the epidermis.

• $p_{2}$ specifies the magnitude of the baseline absorption:⁴^,²²
Eq. (18)
$μ_{a}^{base} (λ) = p_{2} (0.0244 + 8.53 \exp (- \frac{λ - 154}{66.2})) {mm}^{- 1}$
• $p_{3}$ specifies the melanin concentration inside the epidermis and thus $μ_{a, 1}$ :⁴^,²²^,²⁸^–³⁰
Eq. (19)
$μ_{a}^{mel} (λ) = 6.6 \times 10^{11} λ^{- 3.33} {mm}^{- 1}$
Eq. (20)
$μ_{a, 1} (λ) = p_{3} μ_{a}^{mel} + (1 - p_{3}) μ_{a}^{base} {mm}^{- 1}$
• $p_{4}$ specifies the hemoglobin concentration inside the dermis.
• $p_{5}$ specifies the oxygen saturation inside the dermis. $p_{4}$ and $p_{5}$ define $μ_{a, 2}$ :⁴^,²²^,²⁸^,³⁰^,³¹
Eq. (21)
$μ_{a}^{blood} (λ) = p_{5} μ_{a}^{oxy} + (1 - p_{5}) μ_{a}^{deoxy} {mm}^{- 1}$
Eq. (22)
$μ_{a, 2} (λ) = p_{4} μ_{a}^{blood} + (1 - p_{4}) μ_{a}^{base} {mm}^{- 1}$
• $p_{6}$ specifies the oxygen saturation inside the blood vessel and thus $μ_{a, 3}$ :³¹
Eq. (23)
$μ_{a, 3} (λ) = p_{6} μ_{a}^{oxy} + (1 - p_{6}) μ_{a}^{deoxy} {mm}^{- 1}$
• $p_{7}$ specifies the depth of the blood vessel: $z = p_{7} mm$ .
• $p_{8}$ specifies the radius of the blood vessel: $r = p_{8} mm$ .

Note that the proposed equations can easily be changed within the provided MATLAB script as desired since the LUT is generated completely independent of the parameters $p_{1} - p_{6}$ . An example of the user interface can be seen in Fig. 3 with ${p_{1}, \dots, p_{8}} = {0.66, 0.9, 0.005, 0.006, 0.80, 0.60, 0.64, 0.55}$ .

Fig. 3

User interface with sliders to control $p_{1} - p_{8}$ on the right. The final image is displayed on the left. The three plots in the center column show (a) $μ_{s}^{'} (λ)$ , (b) $μ_{a, m} (λ)$ , and (c) $RGB (i)$ .

3. LUT Quality Test

To verify the LUT and especially to check the quality of the $n$ -dimensional interpolation, 450 sample sets ${z, r, μ_{s}^{'}, μ_{a, m}}$ were randomly drawn among the allowed six-dimensional parameter space as listed in Table 1, such that all combinations in between the grid points could be chosen. For each set a regular, independent MC simulation was performed with the exact six parameters and the spatially resolved reflectance $I_{srr}^{MC}$ was recorded using Eq. (1). Next, the LUT was used to generate the 450-corresponding data $I_{srr}^{LUT}$ . The mean relative difference

Eq. (24)

δ_{I} = \frac{1}{N} \sum_{i = 0}^{N} \frac{| I_{srr}^{MC} (i) - I_{srr}^{LUT} (i) |}{I_{srr}^{MC} (i)}

was computed for all sample sets. A histogram of the

δ_{I}

distribution, considering all 450 sample sets, can be found in Fig. 4(a). In Fig. 4(b), the

δ_{I}

distribution for all sample sets with

μ_{s}^{'} > 0.5 {mm}^{- 1}

is shown. For small scattering coefficients, the interpolation may occasionally lead to relative errors up to 9%, but only 4 out of 450 sample sets had a relative error of

δ_{I} > 2 %

. The four outliers were the results for combinations of high absorption coefficients with small reduced scattering coefficients of the dermis, causing

I_{srr}

to be smaller than 0.05, which is an irrelevant range of the reflected intensity for this study. Overall, a very good agreement of the spatially resolved reflectance calculated with a separate MC simulation and the corresponding reflectance obtained by the LUT was found for the vast majority of the parameter range.

Fig. 4

Histogram of the relative errors $δ_{I}$ between $I_{srr}^{MC}$ computed by a conventional MC simulation and $I_{srr}^{LUT}$ obtained by the LUT. The $δ_{I}$ distribution for all 450 randomly drawn samples is shown in (a) and the $δ_{I}$ distribution for the set, where $μ_{s}^{'} > 0.5 {mm}^{- 1}$ , is shown in (b). The insets show close-up views at the indicated positions.

4. Sample Images

In this section, a variety of influences on the visual appearance of a blood vessel are presented using the method introduced before. Further investigations of interest can easily be performed by each user individually after downloading the LUT and the user interface.

4.1.

Depth and Radius of Blood Vessels

A major influence on the appearance of a blood vessel have its geometrical properties, namely its depth $z$ and radius $r$ , especially to elaborate on the question, why veins appear blue although blood itself is red. First, nine images were generated, where all combinations of three distinct depths and three radii were varied. The results for veins with an oxygen saturation of 0.6 for the vessel are shown in Fig. 5 and the images for arteries with an oxygen saturation of 1.0 for the vessel can be found in Fig. 6. Figure 5 shows that for the bluish appearance of a vein a sufficient depth and radius are necessary, compare Figs. 5(e), 5(f), 5(h), and 5(i). This behavior can be explained by the fact that the reduced scattering coefficient and, more pronounced, the absorption coefficient is larger for the blue part of the visible light compared to the red part. First of all, the increasing fraction $\frac{μ_{s}^{'} (λ)}{μ_{a} (λ)}$ for larger wavelengths causes the typical light reddish color of the skin, for areas without a vessel, where more red light is reflected compared to blue light. To understand the bluish appearance of veins, it is important to notice that the spectral course of the optical properties lead to a larger penetration depth of red light and thereby result in a higher probability to reach the vessel compared with blue light. The vein has a high absorption, even for red light, compared to skin and thus absorbs a large part of the light incident onto the vein. Combining these observations induces that the reflection of red light is affected to a larger extend by the vessel compared to blue light in a way that a dip of the fraction $\frac{R (i)}{B (i)}$ is formed for positions $i$ above the vessel. Although the reflection in the red $R (i)$ above the vessel might be larger then the reflection in the blue $B (i)$ above the blood vessel, a vein typically appears bluish according to Land’s retinex theory due to the dip of the fraction $\frac{R (i)}{B (i)}$ , see for details Refs. 2 and 32.

Fig. 5

Color impression of veins depending on the depth of the vessel [left to right: (a)–(c) 0.2 mm, (d)–(f) 0.5 mm, and (g)–(i) 0.8 mm] and the radius of the vessel [top to bottom: (a), (d), and (g) 0.1 mm, (b), (e), and (h) 0.5 mm, and (c), (f), and (i) 0.9 mm]. The remaining parameters are ${p_{1}, \dots, p_{6}} = {0.66, 0.9, 0.005, 0.006, 0.8, 0.6}$ .

Fig. 6

Color impression of an artery depending on the depth of the vessel [left to right: (a)–(c) 0.2 mm, (d)–(f) 0.5 mm, and (g)–(i) 0.8 mm] and the radius of the vessel [top to bottom: (a), (d), and (g) 0.1 mm, (b), (e), and (h) 0.5 mm, and (c), (f), and (i) 0.9 mm]. The remaining parameters are ${p_{1}, \dots, p_{6}} = {0.66, 0.9, 0.005, 0.006, 0.8, 1.0}$ .

Blood vessels that are very close to the surface and have a small radius, such as Figs. 5(a) and 6(a), appear reddish. An explanation is given in Sec. 4.2, where superficial vessels are investigated in detail. For arteries, it can be seen that the bluish appearance is not that pronounced but still present in Figs. 6(h) and 6(i). In general, the arteries appear less bluish compared to veins because the absorption in the vessel is decreasing in the red wavelength range when the oxygen saturation is increased. On top of the reddish appearance, the decreasing absorption for red light is also a reason for the overall worse visibility of arteries compared to veins. For arteries as well as veins, vessels with a rather small radius become almost invisible for larger depths, see Figs. 5(b) and 5(c) as well as Figs. 6(b) and 6(c).

4.2.

Oxygen Saturation and Radius of Superficial Vessels

Next, superficial blood vessels are investigated with a fixed depth at $z = 0.1 mm$ . The oxygen saturation as well as the radius of the vessels were varied in three steps. The resulting images of the nine combinations can be found in Fig. 7. As described previously, the color of a vessel can turn from a bluish to a reddish appearance, when the oxygen saturation is increased. Superficial vessels appear darker by increasing the radius because more light is absorbed in total. In addition, the amount of red light that is absorbed increases at a faster rate compared to blue light, when the radius is increased, due to the larger penetration depth of red light and the color is eventually becoming more and more purple for arteries and bluish for veins. For superficial vessels with a small diameter, as shown in Figs. 7(a)–7(c), the vessels always appear reddish, independent of the oxygen saturation because major parts of the blue light penetrate to the vessel, where they are strongly absorbed. Contrarily, red light is not as strongly absorbed and has a higher chance to transmit through the vessel and ultimately reflect from the tissue.

Fig. 7

Color impression of a superficial blood vessel ( $z = 0.1 mm$ ) depending on the oxygen saturation of the vessel [left to right: (a)–(c) 0.5, (d)–(f) 0.75, and (g)–(i) 1.0] and the radius of the vessel [top to bottom: (a), (d), (g) 0.1 mm, (b), (e), (h) 0.5 mm, and (c), (f), (i) 0.9 mm]. The remaining parameters are ${p_{1}, \dots, p_{5}, p_{7}} = {0.66, 0.9, 0.005, 0.006, 0.8, 0.1}$ .

4.3.

Spectral Dependency of the Reduced Scattering Coefficient

To evaluate the influence of the spectral dependency of the reduced scattering coefficient on the appearance of the blood vessel, different equations found in literature were applied for $μ_{s}^{'}$ :

• Donner06:⁴ $μ_{s}^{'} = 14.74 λ^{- 0.22} + 2.2 \times 10^{11} λ^{- 4} {mm}^{- 1}$
• Jacques98:³³ $μ_{s}^{'} = 2 \times 10^{4} λ^{- 1.5} + 2.2 \times 10^{11} λ^{- 4} {mm}^{- 1}$
• Verdel18:²² $μ_{s}^{'} = 5.5 {(\frac{λ}{500 nm})}^{- 1.5} {mm}^{- 1}$
• Leeuwen18:³ $μ_{s}^{'} = 5.5 {(\frac{λ}{500 nm})}^{- 4} {mm}^{- 1}$
which are plotted in Fig. 8. “Donner06” was used throughout this study. The equation “Leeuwen18” is not explicitly given in their publication. Here, the spectral dependency was kept at $λ^{- 4}$ and the magnitude was adjusted. The corresponding images can be found in Fig. 9. The spectral dependency of $μ_{s}^{'} (λ)$ found by Verdel et al. has only a $λ^{- 1.5}$ dependency, contrarily to the first two, which are given by a combination of a $λ^{- 4}$ dependency with a $λ^{- 0.22}$ dependency and a $λ^{- 1.5}$ dependency, respectively. Therefore, the reduced scattering coefficient “Verdel18” has the smallest difference between the red and blue wavelength range, resulting in an overall lighter reddish appearance of the skin compared to the other images. In addition, the blood vessel appears more blurred since the reduced scattering coefficient “Verdel18” is the largest for the red part of the spectrum. The reduced scattering coefficients “Jacques98” and “Leeuwen18” show the biggest difference between the red and blue wavelength ranges, resulting in an overall bluish touch of the skin. Leeuwen et al. postulated that the Rayleigh scattering of skin has a major influence on the typical bluish color of veins.³ By comparing Figs. 9(c) and 9(d), one can clearly see that even a $λ^{- 1.5}$ dependency of $μ_{s}^{'} (λ)$ leads to a bluish appearance of the vessel. Although, it was stated that the $λ^{- 4}$ dependency of $μ_{s}^{'} (λ)$ is the decisive reason that leads to a significantly higher reflection for the blue channel B over the red channel R above the vessel. However, as Figs. 9(a)–9(d) show, this is not important, because the vessel appears bluish in all images.

5. Discussion

5.1.

Model Selection

In principle, skin consists of multiple layers, such as the stratum corneum, the stratum granulosum, dermal papillae, followed by different subcutaneous tissue layers.²³ In addition, a blood vessel is not only a tube but also is a multilayer system considering the vein wall. On the one hand, this study focuses on the major parameters that influence the visual appearance, such as oxygen saturation or depth and radius of the vessel. On the other hand, an efficient way to obtain a spatially resolved spectrum within a second should be achieved. The model was chosen to optimize these two opposing goals.

Fig. 8

Different spectral dependencies of the reduced scattering coefficient from literature.

Fig. 9

Visual appearance of a blood vessel depending on different spectral dependencies of $μ_{s}^{'} (λ)$ [(a) Donner06, (b) Jacques98, (c) Verdel18, and (d) Leeuwen18]. The chosen parameters are ${p_{1}, \dots, p_{8}} = {1, 0.9, 0.005, 0.006, 0.8, 0.6, 0.4, 0.55}$ .

One of the most debatable restrictions of the used model is the assumption of an infinitely thick dermis layer. In order to investigate the influence of a subcutaneous layer, two MC simulations were performed. The first simulation calculated an image applying the three different tissue types, epidermis, dermis, and blood vessel, as it was used throughout this study. For the second simulation, a subcutaneous layer was added such that the dermis becomes a slab with a thickness of 1.0 mm and the subcutaneous layer was assumed to be semi-infinitely thick. The optical properties for the subcutaneous layer were set to $μ_{a} (λ) = 0.01 {mm}^{- 1}$ and $μ_{s}^{'} (λ) = 1 {mm}^{- 1}$ for all wavelengths.³⁴ For both simulations, the same parameters ${p_{1}, \dots, p_{8}} = {1, 0.9, 0.005, 0.006, 0.8, 0.6, 0.1, 0.35}$ were used and the results are shown in Fig. 10. No major difference of the visual appearance can be seen in Fig. 10(a) when a subcutaneous layer is included (bottom) compared to the original model (top). Looking at the spatially resolved RGB-values in Fig. 10(b) and the corresponding relative difference $\frac{RGB - {RGB}_{sub}}{RGB}$ in Fig. 10(c), only small alterations were found. Therefore, it can be concluded that a subcutaneous layer has only a minor influence on the visual appearance of the model.

Fig. 10

Influence of a subcutaneous layer on the color impression. The color appearance of a blood vessel (a) without subcutaneous layer and (b) with subcutaneous layer. (c) The spatially resolved RGB values for both cases (dotted lines include a subcutaneous layer). (d) The relative differences between the two models.

5.2.

Computational Time

The time needed in MATLAB to generate one image is roughly one second on a single core of a standard computer. A future goal is the usage of the presented LUT to fit a measurement hyperspectrally along a blood vessel. Therefore, an implementation, e.g., in c++ of the multidimensional interpolation might be needed to reduce the fitting time. Another possibility would be the use of neural networks. Previous studies successfully showed the determination of optical properties of tissue using neural networks.³⁵^,³⁶

6. Conclusion

In this study, a software is presented to render the visual appearance of a blood vessel in skin for a variety of relevant parameters of the vessel and the surrounding tissue. To achieve fast image generation, an LUT was generated that stores the spatially resolved reflectance $I_{srr} (i)$ from a two-layered skin model containing a blood vessel. Six parameters can be chosen arbitrarily, namely the depth and radius of the vessel, the scattering coefficient of the epidermis/dermis, and the absorption coefficients for each of the three materials. To significantly reduce the computational time when generating the LUT, a variance reduction method was applied, where the results for all absorptions are generated from a single MC simulation. To reduce the memory demand as well as to increase the interpolation speed, a function, which is the sum of three Gaussian, was fitted to each $I_{srr} (i)$ resulting in an almost exact representation with seven parameters. By obtaining $I_{srr} (λ, i)$ for $λ = 400 nm$ to 800 nm the visual appearance was calculated using the CIE-XYZ color space. A MATLAB user interface is provided, which allows the user to enter physiological and geometrical properties of the model, such as oxygen saturation of the vessel or its radius, which are converted and displayed as an sRGB-image. The correctness of the LUT, especially of the $n$ -dimensional interpolation, was proven. The LUT was used to demonstrate a variety of influences on the visual appearance of veins and arteries.

A main goal of this study was to introduce the new software and demonstrate some applications to give an idea of the opportunities. The combination of the fast evaluation as well as the physical correctness granted by more than 16 millions MC simulations makes this framework applicable for future applications, i.e., for diagnostics based on hyperspectral imaging or educational purposes. Since the LUT depends only on basic physical properties, four optical properties, and two geometrical properties, the users can adjust and implement their own equations linking physiological properties to the physical ones. The short calculation time makes the software usable in real-time applications.

Disclosures

The authors declare that there are no conflicts of interest related to this article.

Acknowledgments

We acknowledge the support of Deutsche Forschungsgemeinschaft (DFG).

References

1.

J. Strutt, “XV. On the light from the sky, its polarization and colour,” London Edinburgh Dublin Philos. Mag. J. Sci., 41 (271), 107 –120 (1871). https://doi.org/10.1080/14786447108640452 Google Scholar

2.

A. Kienle et al., “Why do veins appear blue? A new look at an old question,” Appl. Opt., 35 (7), 1151 –1160 (1996). https://doi.org/10.1364/AO.35.001151 APOPAI 0003-6935 Google Scholar

3.

S. Leeuwen and G. Baranoski, “Elucidating the contribution of Rayleigh scattering to the bluish appearance of veins,” J. Biomed. Opt., 23 (2), 025001 (2018). https://doi.org/10.1117/1.JBO.23.2.025001 JBOPFO 1083-3668 Google Scholar

4.

C. Donner and H. W. Jensen, “A spectral BSSRDF for shading human skin,” in Proc. 17th Eurographics Conf. Rendering Tech., 409 –418 (2006). Google Scholar

5.

E. Angelopoulou, “Understanding the color of human skin,” Proc. SPIE, 4299 243 –251 (2001). https://doi.org/10.1117/12.429495 PSISDG 0277-786X Google Scholar

6.

A. Kienle and R. Hibst, “A new optimal wavelength for treatment of port wine stains?,” Phys. Med. Biol., 40 (10), 1559 –1576 (1995). https://doi.org/10.1088/0031-9155/40/10/001 PHMBA7 0031-9155 Google Scholar

7.

L. Svaasand et al., “Tissue parameters determining the visual appearance of normal skin and port-wine stains,” Lasers Med. Sci., 10 (1), 55 –65 (1995). https://doi.org/10.1007/BF02133165 Google Scholar

8.

R. Olson, J. Gaylor and M. Everett, “Skin color, melanin, and erythema,” Arch. Dermatol., 108 (4), 541 –544 (1973). https://doi.org/10.1001/archderm.1973.01620250029008 Google Scholar

9.

L. Randeberg et al., “The optics of bruising,” Optical-Thermal Response of Laser-Irradiated Tissue, 825 –858 Springer, Dordrecht (2010). Google Scholar

10.

E. d’Eon, D. Luebke and E. Enderton, “Efficient rendering of human skin,” in Proc. 18th Eurographics Conf. Rendering Tech., (2007). Google Scholar

11.

R. Habel, P. Christensen and W. Jarosz, “Photon beam diffusion: a hybrid Monte Carlo method for subsurface scattering,” Comput. Graphics Forum, 32 (4), 27 –37 (2013). https://doi.org/10.1111/cgf.12148 CGFODY 0167-7055 Google Scholar

12.

A. Liemert and A. Kienle, “Light diffusion in N-layered turbid media: frequency and time domains,” J. Biomed. Opt., 15 (2), 025002 (2010). https://doi.org/10.1117/1.3368682 JBOPFO 1083-3668 Google Scholar

13.

A. Liemert, D. Reitzle and A. Kienle, “Analytical solutions of the radiative transport equation for turbid and fluorescent layered media,” Sci. Rep., 7 3819 (2017). https://doi.org/10.1038/s41598-017-02979-4 SRCEC3 2045-2322 Google Scholar

14.

S. Jacques, “MCxyz,” https://omlc.org/software/mc/mcxyz/index.html Google Scholar

15.

A. Doronin and I. Meglinski, “Online object oriented Monte Carlo computational tool for the needs of biomedical optics,” Biomed. Opt. Express, 2 (9), 2461 –2469 (2011). https://doi.org/10.1364/BOE.2.002461 BOEICL 2156-7085 Google Scholar

16.

E. Alerstam, S. Andersson-Engels and T. Svensson, “White Monte Carlo for time-resolved photon migration,” J. Biomed. Opt., 13 (4), 041304 (2008). https://doi.org/10.1117/1.2950319 JBOPFO 1083-3668 Google Scholar

17.

E. Alerstam, T. Svensson and S. Andersson-Engels, “Parallel computing with graphics processing units for high-speed Monte Carlo simulation of photon migration,” J. Biomed. Opt., 13 (6), 060504 (2008). https://doi.org/10.1117/1.3041496 Google Scholar

18.

Q. Fang and D. Boas, “Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units,” Opt. Express, 17 (22), 20178 –20190 (2009). https://doi.org/10.1364/OE.17.020178 OPEXFF 1094-4087 Google Scholar

19.

A. Kienle and M. Patterson, “Determination of the optical properties of turbid media from a single Monte Carlo simulation,” Phys. Med. Biol., 41 (10), 2221 –2227 (1996). https://doi.org/10.1088/0031-9155/41/10/026 PHMBA7 0031-9155 Google Scholar

20.

E. Hecht, Optik, Walter de Gruyter GmbH & Co KG, Berlin (2018). Google Scholar

21.

C. Zoller et al., “Parallelized Monte Carlo software to efficiently simulate the light propagation in arbitrarily shaped objects and aligned scattering media,” J. Biomed. Opt., 23 (6), 065004 (2018). https://doi.org/10.1117/1.JBO.23.6.065004 JBOPFO 1083-3668 Google Scholar

22.

N. Verdel et al., “Noninvasive assessment of skin structure by combined photothermal radiometry and optical spectroscopy: coregistration with multiphoton microscopy,” Appl. Opt., 57 (18), D117 –D122 (2018). https://doi.org/10.1364/AO.57.00D117 APOPAI 0003-6935 Google Scholar

23.

P. Kolarsick, M. Kolarsick and C. Goodwin, “Anatomy and physiology of the skin,” J. Dermatol. Nurses’ Assoc., 3 (4), 203 –213 (2011). https://doi.org/10.1097/JDN.0b013e3182274a98 Google Scholar

24.

S. Agarwal et al., “Ceres solver,” (2018) http://ceres-solver.org Google Scholar

25.

G. Wyszecki and W. S. Stiles, Color Science, 8 Wiley, New York (1982). Google Scholar

26.

M. Anderson et al., “Proposal for a standard default color space for the internet SRGB,” in Color and Imaging Conf., (1996). Google Scholar

27.

B. J. Lindbloom, “Gamma correction,” (2018) http://www.brucelindbloom.com Google Scholar

28.

G. I. Petrov et al., “Human tissue color as viewed in high dynamic range optical spectral transmission measurements,” Biomed. Opt. Express, 3 (9), 2154 –2161 (2012). https://doi.org/10.1364/BOE.3.002154 BOEICL 2156-7085 Google Scholar

29.

T. F. Chen et al., “Hyperspectral modeling of skin appearance,” ACM Trans. Graphics, 34 (3), 1 –14 (2015). https://doi.org/10.1145/2774971 ATGRDF 0730-0301 Google Scholar

30.

I. V. Meglinski and S. J. Matcher, “Quantitative assessment of skin layers absorption and skin reflectance spectra simulation in the visible and near-infrared spectral regions,” Physiol. Meas., 23 (4), 741 –753 (2002). https://doi.org/10.1088/0967-3334/23/4/312 PMEAE3 0967-3334 Google Scholar

31.

S. Jacques and S. Prahl, “Hemoglobin absorption,” (2018) http://omlc.org/ Google Scholar

32.

E. H. Land and J. J. McCann, “Lightness and retinex theory,” J. Opt. Soc. Am. A, 61 (1), 1 –11 (1971). https://doi.org/10.1364/JOSA.61.000001 JOAOD6 0740-3232 Google Scholar

33.

S. Jacques, “Skin optics,” (2018) http://omlc.ogi.edu/news/jan98/skinoptics.html Google Scholar

34.

R. Simpson et al., “Near-infrared optical properties of ex vivo human skin and subcutaneous tissues measured using the Monte Carlo inversion technique,” Phys. Med. Biol., 43 (9), 2465 –2478 (1998). https://doi.org/10.1088/0031-9155/43/9/003 PHMBA7 0031-9155 Google Scholar

35.

I. Matic et al., “Efficient estimation of subdiffusive optical parameters in real time from spatially resolved reflectance by artificial neural networks,” Opt. Lett., 43 (12), 2901 –2904 (2018). https://doi.org/10.1364/OL.43.002901 OPLEDP 0146-9592 Google Scholar

36.

A. Kienle et al., “Spatially resolved absolute diffuse reflectance measurements for noninvasive determination of the optical scattering and absorption coefficients of biological tissue,” Appl. Opt., 35 (13), 2304 –2314 (1996). https://doi.org/10.1364/AO.35.002304 APOPAI 0003-6935 Google Scholar

Biographies of the authors are not available.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Christian Zoller and Alwin Kienle "Fast and precise image generation of blood vessels embedded in skin," Journal of Biomedical Optics 24(1), 015002 (29 January 2019). https://doi.org/10.1117/1.JBO.24.1.015002

Received: 19 November 2018; Accepted: 9 January 2019; Published: 29 January 2019

Access the abstract

JOURNAL ARTICLE
9 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 8 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Blood vessels

Skin

Veins

Absorption

Monte Carlo methods

Reflectivity

Scattering

1.

Introduction

2.

Method