|
1.IntroductionLight that has traveled through tissue carries information about tissue properties, such as structure and biochemical composition. The distance that light travels determines the scale on which the tissue is investigated, as optical properties are averaged over the sampling volume. To detect localized, small-scale changes, techniques with small sampling volumes are, therefore, appropriate. Consequently, short source–detector separations are required to ensure most detected photons will have scattered only once or a few times. In this letter, we will investigate the reflectance measured in this so-called subdiffusive regime, specifically for overlapping source–detector geometries, relevant to, e.g., single fiber spectroscopy measurements. In the diffusive regime, photon direction is randomized and reflectance is independent of the exact shape of the phase function [, the probability distribution of scattering angles]. In contrast, in the subdiffusive regime, the photon direction is not fully randomized; therefore, measurements are sensitive to the phase function.1,2 Advances have been made to derive analytic models for subdiffusive reflectance.3–5 However, for overlapping source–detector geometries, challenges remain regarding incorporation of the tissue phase function.3,4 To model light transport, solutions to the radiative transport equation (RTE) involve expanding the radiance into a series of spherical harmonics and the phase function into Legendre polynomials. The latter are weighted with their moments , where is commonly referred to as the scattering anisotropy. For , the diffusion approximation to the RTE is obtained, with the similarity relation , which expresses that tissues with a different scattering coefficient and a different scattering anisotropy , but equal yield the same reflectance. However, reflectance at short source–detector separations is inadequately described for . Improvement is possible by increasing , giving additional similarity relations. For , Bevilacqua and Depeursinge6 suggested an alternative form of these similarity relations using the parameter Various studies have investigated describing subdiffusive reflectance using in combination with , , and the detector diameter .7–10 However, recent work showed that a range of values (for the same , , and ) can result in the same reflectance.11,12 Therefore, Bodenschatz et al.13 incorporated more similarity relations into the parameter , employing all phase function moments where depends on both the measurement geometry and the scattering regime (empirically, was found to give the best results). The authors concluded that predicts reflectance better than ; however, becomes less predictive for lower numerical apertures (NAs).We consider subdiffusive reflectance as the sum of a diffusive component and a semiballistic component: the diffusive component depends on , and, for the semiballistic component, only photons that have experienced a single backscattering event in combination with an arbitrary number of small-angle forward scattering events are considered. To model this contribution, not all details of the phase function are needed. Rather, we consider the phase function within the NA of the detector. The NA characterizes the range of angles over which a detector can accept incoming photons. We propose a theoretically derived parameter , which is related to the NA of the detector and the integral of the phase function over the corresponding angular interval in the backward () and forward () directions, as We show the theoretical derivation of and the improved prediction of subdiffusive reflectance based on Monte Carlo simulations by compared to and . , therefore, has the potential to improve models used in subdiffusive measurements to extract optical properties.2.Theoretical Derivation ofWe consider semiballistic photons that are backscattered once in combination with an arbitrary number of forward directed scattering events and assume that these photons will only be detected if subsequent scattering events occur at angles smaller or equal to the acceptance angle (Fig. 1), where and is the refractive index of the sample. We introduce the probabilities and as The contribution of semiballistic photons to subdiffusive reflectance scales with the probability of a single backscattering event in the pathlength interval ; this probability equals . To determine the effect of forward semiballistic scattering, we use the result obtained by Wang and Wilson,14 which gives the probability of a photon traveling a pathlength while experiencing forward scattering eventsSumming Eq. (6) over all possible numbers of forward scattering events, from 0 to , we obtain Thus, the combined probability of one backscatter event and an arbitrary number of forward scattering events, each occurring within , is given by integrating over all possible pathlengths, 0 to 3.MethodsFor our Monte Carlo simulations, we modified the software of Prahl et al.15 (which was the core programming later used in Monte Carlo model of steady-state light transport in multilayered tissues (MCML)16 software), to allow the use of arbitrary phase functions using the method of Zijp and ten Bosch.17 The modified code was benchmarked against standard MCML using the Henyey–Greenstein phase function. We modeled geometries where the source and detector geometry were equal and overlapping, where photons were launched from a location based on a uniform distribution across the source with an angle from a uniform angular distribution within the NA. We performed simulations using 15 modified Henyey–Greenstein (mHG),6 144 double Henyey–Greenstein (dbHG),18 8 modified power of cosines (MPC),6 and 46 Reynolds McCormick19 phase functions (RMC, which is equivalent to the Gegenbauer kernel phase function) employing the parameters specified in Table 1 and applying the restrictions and to exclude biologically unreasonable phase functions. Table 1Parameters employed in the selection of phase functions.
For each set of phase functions, simulations were performed for three values of (, ), 1 (, ), and 9 (, ) in combination with an NA of 0.22 or 0.5. An absorption coefficient of and refractive indices inside and outside the sample of 1.35 and 1, respectively, were used. 4.ResultsFigure 2 shows the simulated reflectance versus , , and for an NA of 0.22 and 0.5, and three different values of (0.1, 1, and 9). Reflectance correlates with all three parameters for of 0.1 and 1. This correlation is far less pronounced for . To compare , , and , we determined the spread in each of these three parameters for a chosen reflectance () relative to the total range of the parameter (Table 2). For example, for , NA = 0.22 and a reflectance of 0.001 (), ranges from 1.43 to 1.66. The total range of is 0.68 to 2.43, so the variability is calculated as 0.23/1.75 = 0.13. For all combinations of NA (0.22 and 0.5) and values (0.1 and 1), we determined these variability values for three reflectance values and found the variability to be lowest for compared to and . Table 2Variability of RpNA, σ, and γ for μs′ddet=0.1 and μs′ddet=1, defined as the spread in RpNA, σ, and γ values for a chosen reflectance (±10%) relative to the total range of each parameter.
5.DiscussionWe have theoretically derived the parameter to model subdiffusive light scattering for overlapping source–detector geometries, and we have shown with Monte Carlo simulations that the reflectance depends on . In comparison to and , improves prediction of the reflectance. Our findings indicate that the reflectance does not depend on the details of the phase function—knowledge of all Legendre moments is not required—but only on the magnitude of the phase function within the acceptance angle of the detector (in the backward and forward directions). Currently, the parameter is widely used to model subdiffusive reflectance. Since improves prediction of the measured reflectance, incorporating into subdiffusive models is expected to improve their reliability and thereby the estimation of other optical properties, such as . This improvement is significant because subdiffusive measurements have the potential to detect small-scale tissue changes. For higher values of , the reflectance becomes more diffusive and, therefore, depends less on , , and . The contribution of diffuse photons is a remaining challenge for describing subdiffusive scattering. In our simulations, we kept the absorption coefficient constant. For subdiffusive reflectance, path lengths are short; therefore, the effect of will be minor. While the absorption affects the diffuse contribution to the reflectance, it has a minor effect on semiballistic photons. To implement for subdiffusive measurements, models have to be developed incorporating , the diffusive reflectance, and absorption for specific measurement geometries, such as single fiber reflectance spectroscopy. 6.ConclusionIn conclusion, the theoretically derived parameter predicts subdiffusive light scattering for overlapping source–detector geometries. Consequently, the reflectance does not depend on the details of the entire phase function, but on the phase function within the acceptance angles of the detector. Since improves prediction of the measured reflectance compared to and , the use of is expected to improve derivation of optical properties from subdiffusive measurements. AcknowledgmentsThis research was supported by the Netherlands Organization for Scientific Research (Technology Foundation STW, iMIT-PROSPECT Grant No. 12707) and the Dutch Cancer Society/Alpe d’HuZes (Project No. 2014-7009). We would like to thank the Organization Committee of the Spinoza Chair from the Academic Medical Center of the University of Amsterdam for funding a “Spinoza” lectureship for Steven L. Jacques. ReferencesJ. R. Mourant et al.,
“Influence of the scattering phase function on light transport measurements in turbid media performed with small source-detector separations,”
Opt. Lett., 21 546
–548
(1996). https://doi.org/10.1364/OL.21.000546 OPLEDP 0146-9592 Google Scholar
A. Kienle, F. K. Forster and R. Hibst,
“Influence of the phase function on determination of the optical properties of biological tissue by spatially resolved reflectance,”
Opt. Lett., 26 1571
–1573
(2001). https://doi.org/10.1364/OL.26.001571 OPLEDP 0146-9592 Google Scholar
E. Vitkin et al.,
“Photon diffusion near the point-of-entry in anisotropically scattering turbid media,”
Nat. Commun., 2 587
–604
(2012). https://doi.org/10.1038/ncomms1599 Google Scholar
M. Jia et al.,
“Virtual-source diffusion approximation for enhanced near-field modeling of photon-migration in low-albedo medium,”
Opt. Express, 23 1337
–1352
(2015). https://doi.org/10.1364/OE.23.001337 OPEXFF 1094-4087 Google Scholar
D. Piao and S. Patel,
“Simple empirical master-slave dual-source configuration within the diffusion approximation enhances modeling of spatially resolved diffuse reflectance at short-path and with low scattering from a semi-infinite homogeneous medium,”
Appl. Opt., 56 1447
(2017). https://doi.org/10.1364/AO.56.001447 APOPAI 0003-6935 Google Scholar
F. Bevilacqua and C. Depeursinge,
“Monte Carlo study of diffuse reflectance at source-detector separations close to one transport mean free path,”
J. Opt. Soc. Am. A, 16 2935
(1999). https://doi.org/10.1364/JOSAA.16.002935 JOAOD6 0740-3232 Google Scholar
S. C. Kanick, H. J. C. M. Sterenborg and A. Amelink,
“Empirical model of the photon path length for a single fiber reflectance spectroscopy device,”
Opt. Express, 17 860
–871
(2009). https://doi.org/10.1364/OE.17.000860 OPEXFF 1094-4087 Google Scholar
U. A. Gamm et al.,
“Quantification of the reduced scattering coefficient and phase-function-dependent parameter of turbid media using multidiameter single fiber reflectance spectroscopy: experimental validation,”
Opt. Lett., 37 1838
–1840
(2012). https://doi.org/10.1364/OL.37.001838 OPLEDP 0146-9592 Google Scholar
S. C. Kanick et al.,
“Sub-diffusive scattering parameter maps recovered using wide-field high-frequency structured light imaging,”
Biomed. Opt. Express, 5 3376
(2014). https://doi.org/10.1364/BOE.5.003376 BOEICL 2156-7085 Google Scholar
F. van Leeuwen-van Zaane et al.,
“In vivo quantification of the scattering properties of tissue using multi-diameter single fiber reflectance spectroscopy,”
Biomed. Opt. Express, 4 696
–708
(2013). https://doi.org/10.1364/BOE.4.000696 BOEICL 2156-7085 Google Scholar
N. Bodenschatz et al.,
“Model-based analysis on the influence of spatial frequency selection in spatial frequency domain imaging,”
Appl. Opt., 54 6725
(2015). https://doi.org/10.1364/AO.54.006725 APOPAI 0003-6935 Google Scholar
P. Naglic et al.,
“Estimation of optical properties by spatially resolved reflectance spectroscopy in the subdiffusive regime,”
J. Biomed. Opt., 21 095003
(2016). https://doi.org/10.1117/1.JBO.21.9.095003 JBOPFO 1083-3668 Google Scholar
N. Bodenschatz et al.,
“Quantifying phase function influence in subdiffusively backscattered light,”
J. Biomed. Opt., 21 035002
(2016). https://doi.org/10.1117/1.JBO.21.3.035002 JBOPFO 1083-3668 Google Scholar
R. K. Wang and M. J. Wilson,
“Vertex/propagator model for least-scattered photons traversing a turbid medium,”
J. Opt. Soc. Am. A, 18 224
–231
(2001). https://doi.org/10.1364/JOSAA.18.000224 Google Scholar
S. A. Prahl et al.,
“A Monte Carlo model of light propagation in tissue,”
Dosimetry of Laser Radiation in Medicine and Biology, IS 102
–111 SPIE Press, Bellingham, Washington
(1989). Google Scholar
L. V. Wang, S. L. Jacques and L. Zheng,
“MCML–Monte Carlo modeling of light transport in multi-layered tissues,”
Comput. Methods Programs Biomed., 47 131
–146
(1995). https://doi.org/10.1016/0169-2607(95)01640-F CMPBEK 0169-2607 Google Scholar
J. R. Zijp and J. J. ten Bosch,
“Use of tabulated cumulative density functions to generate pseudorandom numbers obeying specific distributions for Monte Carlo simulations,”
Appl. Opt., 33 533
(1994). https://doi.org/10.1364/AO.33.000533 APOPAI 0003-6935 Google Scholar
P. Saccomandi et al.,
“Estimation of anisotropy coefficient of swine pancreas, liver and muscle at 1064 nm based on goniometric technique,”
J. Biophotonics, 8
(5), 422
–428
(2015). https://doi.org/10.1002/jbio.v8.5 Google Scholar
L. O. Reynolds and N. J. McCormick,
“Approximate two-parameter phase function for light scattering,”
J. Opt. Soc. Am., 70 1206
–1212
(1980). https://doi.org/10.1364/JOSA.70.001206 JOSAAH 0030-3941 Google Scholar
|