12 February 2014 Linear dependency of full scattering profile isobaric point on tissue diameter
Author Affiliations +
J. of Biomedical Optics, 19(2), 026007 (2014). doi:10.1117/1.JBO.19.2.026007
Abstract
Most methods for measuring light-tissue interaction focus on volume reflectance, while very few measure light transmission. In a previous work, we suggested investigating the influence of blood vessel diameter on photons exiting the tissue at all exit angles to receive the full scattering profile. By this method, we have shown that there is a central angle, i.e., the isobaric point, independent of blood vessel diameter. The vessel diameter changes the effective reduced scattering coefficient. However, both the scattering profile and the value of the isobaric point strongly depend on optical properties and the exact geometry of the tissue. In this study, we investigate the dependency of the isobaric point on tissue diameter and scattering coefficient in both two-dimensional and three-dimensional simulations. We show that the value of this point linearly depends on tissue diameter. The findings of this work solve the dilemma of whether to measure transmission or reflection since the isobaric point reduces by half the total amount of exiting photons. Furthermore, the full scattering profile is sensitive to changes in the scattering properties, but a single isobaric point to these changes is expected. If this point is not found, it is a diagnostic indication of an unexpected change in the tissue.
Duadi, Feder, and Fixler: Linear dependency of full scattering profile isobaric point on tissue diameter

1.

Introduction

Human tissue is one of the most complex optical media since it is turbid and nonhomogeneous. Hence, optical properties, such as absorption coefficient μa, scattering coefficient μs, and anisotropy factor g, are unknown and vary in different areas and physiological states. Furthermore, in respiration, blood vessels, which are the main cause for light absorption and also cause strong scattering, vary in size.

Measurement of light-tissue interaction uses either direct methods that need no theoretical models of light propagation or indirect methods. Direct methods include attenuation measurements in calibration solutions1 or integrating spheres for the measurement of total transmitted light through a blood sample and total reflected light from it.23.4 Indirect methods are based on a mathematical simulation model of photon migration. Photon migration simulations are commonly described by the diffusion theory.5 However, in the close region near the light source, diffusion theory does not accurately describe the light distribution.6 The Monte Carlo (MC) simulation method assumes that all photons begin as ballistic photons and that change in their direction is due to the scattering coefficient of the medium.78.9.10.11

Most optical-physiological diagnostic methods are based on the insertion of light with known parameters to a tested tissue, followed by the measurement of the re-emitted light. When investigating light–tissue interaction, human tissue is usually dealt with as a semi-infinite surface. Hence, the investigation is commonly based on the light reflected from the tissue (volume reflection).7,1213.14.15.16 Very few methods test the transmitted light23.4,17 or ballistic photons.8 It has been previously suggested18 looking for the first time at the full profile of light that is scattered (transmitted and reflected) at all possible exit angles from a circle of tissue, such as a fingertip joint, ear lobe, or pinched tissue. In this work, we found that the full scattering profile is affected by blood vessel diameter. This effect can be explained by the change in the effective scattering coefficient. This parameter is influenced by both tissue and blood vessels, where larger blood vessels are less effective in light absorption because of the shielding effect.17 However, we have also made a unique discovery that there is an angle (i.e., the isobaric point) that is indifferent to these changes. We have mentioned that the value of this isobaric point depends on optical properties and exact geometry; however, further investigation is necessary.

In this work, we establish the meaning and usefulness of this isobaric point and calculate the change in its value in both two-dimensional (2-D) and three-dimensional (3-D) simulations in a simplified tissue model. In this model, we replaced the change in blood vessel diameter with the change in the effective reduced scattering coefficient. Since absorption due to blood vessels does not change the shape of the full scattering profile, it can influence the total scattered intensity and is neglected in the current simulations. We calculate the full scattering profile via MC simulation of a circular cross-section of tissue13,18 in the near-infrared (NIR) wavelength of 850 nm since it is common in photoplethysmography (PPG) measurements due to its high penetration depth and minimal scattering.19,20 We investigate the dependency of the isobaric point in case of different tissue diameters and show that this phenomenon is independent of the value of the anisotropy factor. The new isobaric point can be implemented in NIR spectroscopy, PPG experiments, and analyzing exact oxygen saturation values.

2.

Materials and Methods

2.1.

2-D Model of Circular Tissue

To understand the influence of tissue diameter on the isobaric point, several simulations were done on a 2-D model of a cross-section of circular tissue (Fig. 1). A beam of photons with a waist of w0=1.5mm enters the circle of tissue parallel to the z direction. The propagation path of each photon is calculated from the scattering constant, assuming the absorption is negligible. An MC simulation of photon migration within irradiated tissues was built18 in order to calculate the full scattering profile at all possible exit angles. This simulation is based on the assumption that all photons reaching the tissue begin as ballistic photons. Given the current photon’s direction (θold), and the probability of a photon to scatter [1exp(μsdr)], if the photon scattered, its new direction (θnew) was calculated using

θnew=θold+s×cos(g),
where s is a random number from the group {1,1}.

Fig. 1

Scheme of simulated two-dimensional (2-D) model. A beam of photons (dashed arrows) with a waist of w0 was introduced at one end and collected when exiting the tissue at different central angles (θ).

JBO_19_2_026007_f001.png

This process is repeated until the photon exits the tissue, and then its location is saved.

Several repetitions of the simulation for each tissue diameter were held to determine the variations in transmission due to the change in the reduced scattering coefficient μs=μs(1g). The values of the reduced scattering coefficient (μs) were in the range of the human skin values of 2 to 26cm1.21

2.2.

3-D Model of Cylindrical Tissue

A 3-D model of a cylinder of tissue was created with the same method and properties as described in the 2-D model (Fig. 2). Ten repetitions of the 3-D MC simulation13 were built in order to calculate the full scattering profile at all possible exit angles. Simulations for each tissue diameter (4, 5, and 6 mm) were held to determine the variations in transmission due to the change in reduced scattering coefficient. The illumination was set to be a square of 1.5×1.5mm, and the transmitted photons were collected from a 1-mm-thick slice (gray area in Fig. 2) along the cylinder surface.

Fig. 2

The cross-section of simulated three-dimensional (3-D) tissue model. Photons are collected from the surface of a 1-mm-thick area of interest.

JBO_19_2_026007_f002.png

3.

Results

3.1.

Isobaric Point for Different Tissue Diameters in 2-D

The 2-D simulation described in Sec. 2.1 was held for several tissue diameters. The mean fraction of photons that exited the tissue at several central angles (θ in Fig. 1) from 0 to 180 deg is presented in Figs. 3(a), 3(b), and 3(c) for tissue diameters of 5, 10, and 15 mm, respectively. The mean of 10 different simulations was calculated. The value of the isobaric point, as extracted from these graphs for 5, 10, and 15 mm, is 115, 135, and 165 deg, respectively. As expected, as tissue diameter increases, so does the value of the isobaric point. In order to facilitate the presentation and to compare it to known results from the literature,9,18 reflected photons are defined as those that exit the tissue near the entrance point at a central angle of |θ|>135deg, and forward transmitted photons are defined as those that exit the tissue at a central angle of |θ|<45deg. Figures 4(a), 4(b), and 4(c) present the average fraction of the reflected (squares) and forward transmitted photons (diamonds), respectively, versus scattering coefficients, for several tissue diameters.

Fig. 3

The full scattering profile in 2-D simulation for various tissue diameters: (a) Dt=5mm; (b) Dt=10mm; and (c) Dt=15mm. Reduced scattering coefficients (in cm1) are noted in the figure legend on the right.

JBO_19_2_026007_f003.png

Fig. 4

The fixed angle transmission (diamonds) and reflection (squares) in 2-D simulation for different tissue diameters: (a) Dt=5mm; (b) Dt=10mm; and (c) Dt=15mm.

JBO_19_2_026007_f004.png

3.2.

Isobaric Point for Different Tissue Diameters in 3-D

The 3-D simulation described in Sec. 2.2 was held for several tissue diameters. The mean fraction of photons that exited the tissue at several central angles (θ in Fig. 2) from 0 to 180 deg is presented in Figs. 5(a), 5(b), and 5(c) for tissue diameters of 4, 5, and 6 mm, respectively. The value of the isobaric point, as extracted from these graphs for 4, 5, and 6 mm, is 155, 158, and 160 deg, respectively.

Fig. 5

The full scattering profile for various reduced scattering coefficients (square, triangle, x, and diamond in respect to 8, 10, 16, and 26cm1). Tissue diamater in 3-D simulation is (a) Dt=4mm; (b) Dt=5mm; and (c) Dt=6mm.

JBO_19_2_026007_f005.png

Since the photons are collected from a 1-mm-thick slice of tissue, there is no point at looking at the percent of forward transmitted photons since it is too low. Hence, a comparison of the reflection as defined in Sec. 3.1 is presented in Fig. 6. Note that although the change in tissue diameter hardly affects (<5%) the fixed angle reflection in Fig. 6, it significantly changes the maximal value in the full scattering profile (>30%) for the same scattering coefficient in Fig. 5.

Fig. 6

Reflection versus reduced scattering coefficient for different tissue diameters (triangle, diamond, and square in respect to 4, 5, and 6 mm).

JBO_19_2_026007_f006.png

In order to verify that the value of the isobaric point does not change via anisotropy factor (g), which is also a variable in photon migration, and changes between different tissue areas (for example, fate as opposed to blood), three different simulations were held. In these simulations for a 4-mm diameter of tissue, the reduced scattering coefficient was the same but g changed, and the scattering coefficient was set to compensate for this change using the relation μs=μs(1g). As presented in Fig. 7, this change significantly influences the fixed angle reflection (up to 10%), but does not change the value of the isobaric point (data not shown).

Fig. 7

Reflection versus reduced scattering coefficient for different anisotropy factors (square, diamond, and triangle in respect to 0.8, 0.9, and 0.95) for a 3-D simulation with a 4-mm tissue diameter.

JBO_19_2_026007_f007.png

Dependency of the isobaric point on the tissue diameter in both 2-D and 3-D simulations is presented in Fig. 8. As expected, the value of the isobaric point in the 3-D simulation is different from that obtained in 2-D due to the scattering in the additional dimension. According to these simulations, we have found a linear dependency (with a 98% fitting) to tissue diameter.

Fig. 8

Value of the isobaric point versus tissue diameter in 2-D (square) and 3-D (diamond) simulation.

JBO_19_2_026007_f008.png

4.

Discussion

In our current work, we have investigated the full scattering profile in different tissue diameters in 2-D, as well as 3-D simulations. From these profiles, we have extracted the value of the isobaric point to reduced scattering coefficient. Although the isobaric point value was different in the 2-D as well as 3-D simulations, a linear dependency was found to the tissue diameter.

When measuring light-tissue interactions, the first question that arises is whether to measure reflected or transmitted light. The common method for measuring blood saturation is from volume reflection7,1213.14.15.16 since there is a very high amount of transmitted light from thick tissues. Furthermore, diffusion reflection methods,67.8.9.10.11.12.13.14 which are used for the investigation of tissue’s optical parameters, are based on the assumption that the tissue is semi-infinite.

Very few methods test the transmitted light23.4,17 or ballistic photons8 because of low light intensity. By measuring the isobaric point, we solve this dilemma since this is the point that truly divides by half the total amount of exiting photons.

In our current work, we have chosen an NIR wavelength of 850 nm since it is commonly used in PPG measurements19,22 and because the range of scattering coefficients was according to known ranges in this wavelength.4,21,22 However, if one wants to investigate the full scattering coefficient in a lower reduced scattering coefficient, one may generate a completely opposite profile. For example, in our published results,18 the dependency between the reflectance and blood vessel diameter in a circular tissue illuminated at 850 nm was found to be opposite to that found by Jacques9 from a semi-infinite tissue illuminated at 585 nm since the tissue scattering is much lower at this wavelength.

In order to demonstrate the effect of lower scattering, the 2-D simulation where Dt=5mm (Fig. 9) was repeated for a wider range of reduced scattering coefficients (diamond, triangle, asterisk, cross, circle, x, and square in respect to 2, 6, 10, 14, 18, 22, and 26cm1). While in a lower reduced scattering coefficient (for example, 2cm1 marked by diamonds), the light is mainly transmitted; in a higher reduced scattering coefficient (for example, 26cm1 marked by squares), the light is mainly reflected. Furthermore, the isobaric point is only relevant in higher reduced scattering coefficients.

Fig. 9

The full scattering profile in 2-D simulation for a wide range of reduced scattering coefficients. For a tissue diamater of Dt=5mm reduced scattering coefficients (diamond, triangle, asterisk, cross, circle, x, and square in respect to 2, 6, 10, 14, 18, 22, and 26cm1) in a wider range may completely alter the profile.

JBO_19_2_026007_f009.png

In conclusion, this work allows predicting the isobaric point, which is optimal for different types of NIR spectroscopy measurements for a given tissue with known optical properties. Furthermore, we have presented a method that is sensitive to the reduced scattering coefficient. Therefore, if experimental measurements will show a different full scattering profile, which does not cross the isobaric point, this indicates an abnormal change in the optical tissue properties.

References

1. 

A. N. Yaroslavskyet al., Optics of Blood, SPIE Press, Bellingham, Washington (2002).Google Scholar

2. 

M. Friebelet al., “Influence of oxygen saturation on the optical scattering properties of human red blood cells in the spectral range 250 to 2000 nm,” J. Biomed. Opt. 14(3), 034001 (2009).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3127200Google Scholar

3. 

M. Friebelet al., “Determination of optical properties of human blood in the spectral range 250 to 1100 nm using Monte Carlo simulations with hematocrit-dependent effective scattering phase functions,” J. Biomed. Opt. 11(3), 034021 (2006).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2203659Google Scholar

4. 

A. Rogganet al., “Optical properties of circulatin human blood in the wavelength range 400–2500 nm,” J. Biomed. Opt. 4(1), 36–46 (1999).JBOPFO1083-3668http://dx.doi.org/10.1117/1.429919Google Scholar

5. 

G. ZaccantiS. Del BiancoF. Martelli, “Measurements of optical properties of high-density media,” Appl. Opt. 42(19), 4023–4030 (2003).APOPAI0003-6935http://dx.doi.org/10.1364/AO.42.004023Google Scholar

6. 

S. L. JacquesB. W. Pogue, “Tutorial on diffuse light transport,” J. Biomed. Opt. 13(4), 041302 (2008).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2967535Google Scholar

7. 

R. Ankriet al., “In-vivo tumor detection using diffusion reflection measurements of targeted gold nanorods—a quantitative study,” J. Biophotonics 5(3), 263–273 (2012).JBOIBX1864-063Xhttp://dx.doi.org/10.1002/jbio.201100120Google Scholar

8. 

R. VeredS. HavlinH. Taitelbaum, “Optical detection of hidden tumors,” Proc. SPIE 2389, 851–858 (1995).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.210035Google Scholar

9. 

S. L. Jacques, “Optical assessment of cutaneous blood volume depends on the vessel size distribution: a computer simulation study,” J. Biophotonics 3(1–2), 75–81 (2010).JBOIBX1864-063Xhttp://dx.doi.org/10.1002/jbio.200900085Google Scholar

10. 

L. WangS. L. Jacques, “Hybrid model of Monte Carlo simulation and diffusion theory for light reflectance by turbid media,” J. Opt. Soc. Am. A 10(8), 1746–1752 (1993).JOAOD60740-3232http://dx.doi.org/10.1364/JOSAA.10.001746Google Scholar

11. 

L. WangS. L. JacquesL. Zheng, “MCML—Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed. 47(2), 131–146 (1995).CMPBEK0169-2607http://dx.doi.org/10.1016/0169-2607(95)01640-FGoogle Scholar

12. 

R. AnkriH. TaitelbaumD. Fixler, “Reflected light intensity profile of two-layer tissues: phantom experiments,” J. Biomed. Opt. 16(8), 085001 (2011).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3605694Google Scholar

13. 

D. FixlerR. Ankri, “Subcutaneous gold nanorod detection with diffusion reflection measurement,” J. Biomed. Opt. 18(6), 061226 (2013).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.18.6.061226Google Scholar

14. 

R. A. J. GroenhuisA. H. FerwedaJ. J. Ten Bosch, “Scattering and absorption of turbid materials determined from reflection measurements. 1: Theory,” Appl. Opt. 22(16), 2456–2462 (1983).APOPAI0003-6935http://dx.doi.org/10.1364/AO.22.002456Google Scholar

15. 

D. Jakubowskiet al., “Quantitative absorption and scattering spectra in thick tissues using broadband diffuse optical spectroscopy,” Chapter 12 in Biomedical Optical Imaging, J. G. FujimotoD. L. Farkas, Eds., pp. 330–355, Oxford University Press, New York (2009).Google Scholar

16. 

T. H. Phamet al., “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).RSINAK0034-6748http://dx.doi.org/10.1063/1.1150665Google Scholar

17. 

L. ZhangA. ShiH. Lu, “Determination of optical coefficients of biological tissue from a single integrating-sphere,” J. Mod. Opt. 59(2), 121–125 (2012).JMOPEW0950-0340http://dx.doi.org/10.1080/09500340.2011.631052Google Scholar

18. 

H. DuadiD. FixlerandR. Popovtzer, “Dependence of light scattering profile in tissue on blood vessel diameter and distribution: a computer simulation study,” J. Biomed. Opt. 18(11), 111408 (2013).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.18.11.111408Google Scholar

19. 

M. NitzanS. Engelberg, “Three-wavelength technique for the measurement of oxygen saturation in arterial blood and in venous blood,” J. Biomed. Opt. 14(2), 024046 (2009).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3120496Google Scholar

20. 

S. Del BiancoF. MartelliG. Zaccanti, “Penetration depth of light re-emitted by a diffusive medium: theoretical and experimental investigation,” Phys. Med. Biol. 47(23), 4131–4144 (2002).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/47/23/301Google Scholar

21. 

T. ListerP. A. WrightP. H. Chappell, “Optical properties of human skin,” J. Biomed. Opt. 17(9), 090901 (2012).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.17.9.090901Google Scholar

22. 

A. N. Bashkatovet al., “Optical properties of human skin, subcutaneous and mucous tissues in the wavelength range from 400 to 2000 nm,” J. Phys. D: Appl. Phys. 38(15), 2543–2555 (2005).JPAPBE0022-3727http://dx.doi.org/10.1088/0022-3727/38/15/004Google Scholar

Biography

Hamootal Duadi received her PhD degree in 2011 from the Engineering Faculty, Bar-Ilan University, Israel. She specialized in optical information processing, super resolution, diffractive optical elements and beam shaping, 3D and range estimation. Her current research fields are pulse oximetry, tissue light interaction, and live cell imaging.

Idit Feder received her BSc degree in 2013 from the Faculty of Engineering of Bar-Ilan University, and currently she is working toward the Msc degree and researching the dependency of the scattering profile of tissues.

Dror Fixler received his PhD degree in 2003 from the Department of Physics, Bar-Ilan University, Israel. He is a member of the Faculty of Engineering and the Nano Center of Bar-Ilan University. He has published over 50 original research papers and holds over 9 issued patents. His research interests include fluorescence measurements (FLIM and anisotropy decay), optical super resolution, high-end electro-optical system engineering, and light-tissue interaction. He received several international awards and organized and presented at over 20 international conferences.

Hamootal Duadi, Idit Feder, Dror Fixler, "Linear dependency of full scattering profile isobaric point on tissue diameter," Journal of Biomedical Optics 19(2), 026007 (12 February 2014). http://dx.doi.org/10.1117/1.JBO.19.2.026007
Submission: Received ; Accepted
JOURNAL ARTICLE
6 PAGES


SHARE
KEYWORDS
Scattering

Tissue optics

Tissues

Photons

Light scattering

3D modeling

Blood vessels

Back to Top