Search tips
Search criteria 


Logo of boeaboutauthor infoeditorial boardsearchOptics InfoBaseboeThis article
Biomed Opt Express. 2012 May 1; 3(5): 1086–1100.
Published online 2012 April 24. doi:  10.1364/BOE.3.001086
PMCID: PMC3342184

Scattering phase function spectrum makes reflectance spectrum measured from Intralipid phantoms and tissue sensitive to the device detection geometry


Reflectance spectra measured in Intralipid (IL) close to the source are sensitive to wavelength-dependent changes in reduced scattering coefficient equation m1 and scattering phase function (PF). Experiments and simulations were performed using device designs with either single or separate optical fibers for delivery and collection of light in varying concentrations of IL. Spectral reflectance is not consistently linear with varying IL concentration, with PF-dependent effects observed for single fiber devices with diameters smaller than ten transport lengths and for separate source-detector devices that collected light at less than half of a transport length from the source. Similar effects are thought to be seen in tissue, limiting the ability to quantitatively compare spectra from different devices without compensation.

OCIS codes: (170.3660) Light propagation in tissues, (170.7050) Turbid media, (170.6510) Spectroscopy, tissue diagnostics

1. Introduction

Reflectance spectroscopy provides estimates of optical properties within a turbid medium such as tissue. This modality has an important role in the field of biomedical optics, as reflectance spectra can be used to extract information about aspects of the tissue physiology (e.g. vascular oxygen saturation and blood volume fraction) [1,2] and morphology (e.g. cell, organelle size and density) [36] and can be used to correct raw fluorescence measurements for the influences of background scattering and absorption [79]. Spectra measured in reference media such as the fat emulsion commercialized as Intralipid (IL) are often used to calibrate or understand instrument performance [10], yet the spectra recorded by different devices can be dramatically affected by device-specific sensitivities to the scattering properties. In this study, the influence of detection geometry with respect to IL concentration is studied in detail using measurements and Monte Carlo studies, to illustrate and interpret how the spectra can change.

Fiber optic reflectance devices are designed to sample light from a location of interest where typically the size and shape is customized for the biological application. Devices that collect reflectance remitted at large distances from the source can achieve (relatively) deep sampling and therefore are capable of monitoring volume averaged quantities in bulk tissues such as blood oxygenation in the brain [11] or breast [12]. Conversely, devices that collect reflectance close to the light source sample highly localized volumes, and therefore can provide selective sampling of superficial tissue layers [1319]. Validation of these devices with tissue-simulating optical phantoms in vitro is critical to allow characterization of the observed reflectance in response to controlled optical properties within the phantom. Phantom measurements may also be a valuable method to calibrate the wavelength-specific device response for use in analyzing measurements of media with unknown optical properties, such as tissue. Perhaps the most commonly adopted scattering phantom is IL [10,20], because it is stable across batches, readily purchased and has been extensively investigated to characterize the wavelength-dependent optical properties [2124]. Dilutions of IL can yield wavelength-dependent reduced scattering coefficients (equation m2) that are similar to values reported in tissue. These characteristics have led to wide acceptance of IL as a tissue-simulating scattering standard for the investigation of many types of reflectance devices [13,1519]. However, as this study shows, spectra measured in the same IL phantom by different reflectance geometries may show substantially different spectral features.

To understand how scattering properties influence reflectance spectra measured by different device geometries differently, it is important to consider the influence of scattering on the transport of light being collected during measurements. Light that travels relatively large distances in scattering media is termed diffuse, with photon transport occurring in the direction of fluence (concentration) gradients within the medium [25,26]. Collected diffuse reflectance intensities can be described analytically by the diffusion approximation to the radiative transport equation, which allows estimation of both equation m3 and the absorption coefficient equation m4 [18,25]. The heuristic guideline for diffuse transport is light that has travelled 1 to 2 mean free transport lengths from the source [26]; here, the term transport length is defined as equation m5. Analysis of diffuse light is insensitive to the form of the scattering phase function (PF), as changes to the scatter frequency (given by equation m6) and directionality (given by PF) interchangeably influence diffuse light transport in the form of the reduced scattering coefficient equation m7, where g is the anisotropy given as equation m8. As the distance between source and collection location decreases, the sampled transport regime shifts from the collection of multiply scattered light to the collection of photons that may have undergone few (or single) scattering events [27,28]. Light transport in this semi- to non-diffuse regime is not exclusively characterized by equation m9, but is also influenced by the probability of large angle backscattering events, which is defined by the scattering PF [28,29].

The influence of PF on light transport collected at various source-detector separations (ρ) is illustrated in Fig. 1 . Light collected at large ρ is multiply scattered, making it insensitive to any individual directional scattering event. Light collected at small ρ is sensitive to large angle scattering events. The frequency of these events is proportional to the magnitude of the backscatter-lobe of the PF. In IL phantoms, it is important to consider that both equation m10and PF change considerably over the ultraviolet-visible (UV-vis) wavelength range; these changes are displayed in Fig. 2 (A) and (B) , respectively. While it is well known that equation m11exhibits a power law relationship with wavelength [30,31], changes to PF vs. wavelength are less commonly understood. In Fig. 2 (B), PFs from selected wavelengths are plotted as the probability density distributions of scattering angles; these plots show marked differences between wavelengths [24]. Previous studies have developed metrics to describe the average directional scattering behavior associated with a PF [28,32]. Specifically, these studies utilized Legendre moments of the PF distributions to characterize the average direction of scatter, given by the anisotropy (equation m12), which is shown in Fig. 2 (C), and the likelihood of large angle backscatter events, given by a parameter equation m13 which includes information about the first two moments of the PF, equation m14and equation m15, respectively and is shown in Fig. 2 (D). For forward-directed scattering media, such as IL, γ can be interpreted as follows: decreases in values of γ are associated with an increased likelihood of large angle backscatter events. A simple approach to characterizing the differences in the likelihood of backscatter between PFs is to consider the probability of a large angle backscatter event, which we define as an scattering angle that would fall within the cone of acceptance for the fiber probes included in this study. Consideration of these large angle events is especially relevant for the single scattering regime. Figure 2 (E) shows that for IL, the likelihood of large angle backscattering events increases substantially across the UV-vis wavelength range; correspondingly, γ decreases over the same range. Previous studies have utilized the moments of the PF to describe collected light intensity near the source [28,3234]. But to the best of the author’s knowledge, these findings have not been translated to the characterization of scattering-induced spectral features for devices that measure over different length scales close to the source, even for a well-characterized scattering standard such as IL.

Fig. 1
Fiber optic geometries of reflectance devices that measure light close to the source fiber. Single fiber measurements obtained using same fiber optic to deliver and collect light with diameter given by df. Separate source detector devices have core-to-core ...
Fig. 2
Optical properties of Intralipid optical phantoms. A) shows reduced scattering coefficient vs. wavelength for selected IL dilutions. B) scattering phase function of IL at selected wavelengths. Right column shows wavelength-dependent metrics of PF: C) ...

This study characterizes scattering features observed in reflectance spectra measured in IL optical phantoms. Experimental and simulation data are used to identify and investigate the spectral response of reflectance measured in a wide range of IL dilutions as measured by fiber optic device geometries with 1) overlapping source and detector with a range of fiber diameters, and 2) separate source and detector fibers with a wide range of source-detector separations. The underlying mechanisms that link the observed spectral response to the wavelength-dependent scattering properties are presented, and the potential errors associated with a failure to consider PF-related effects for measurements close to the source are discussed.

2. Methods

2.1. Single fiber reflectance experimental setup

The experimental apparatus for the reflectance device that uses a single fiber optic to deliver and collect light during measurements has been described in detail previously [35,36]. The device consists of a halogen light source (HL-2000-FHSA, Ocean Optics, Duiven, Netherlands), a spectrometer (SD2000, Ocean Optics, Duiven, Netherlands) and a solid core single fiber probe that is connected to the light source and the spectrometer via a bifurcated fiber. The probe tip is polished at an angle of 15 degrees to remove specular reflections within the collected reflectance intensity due to the index of refraction mismatch at the probe/medium interface. This study used a calibration procedure to account for variations in the output of the lamp, transmission characteristics of the fiber, and sensitivity of the spectrometer. Calibration measurements were performed in a calibration phantom containing 20% IL diluted to a volume particle concentration of 1.5%, which yielded equation m16 (800 nm) = 1.2 mm−1 (equation m17), and in a dark container filled with water (equation m18):

equation m19

Here, equation m20is the raw signal of a sampled optical phantom, and equation m21 is the absolute single fiber reflectance (SFR) signal, expressed in % of incident photons, obtained by Monte Carlo (MC) simulations of the IL calibration phantom. Importantly, for these simulations the wavelength-specific PFs and equation m22 values reported by [24] were used. The calculation results in the calibrated single fiber reflectance intensity, equation m23.

2.2. Intralipid tissue simulating phantom preparation

The liquid phantoms were prepared by mixing 0.9% NaCl (Baxter, Utrecht, Netherlands) with different amounts of Intralipid 20% (Fresenius Kabi, s-Hertogenbosch, Netherlands); this stock solution has a volume particle concentration percentage (ILp) of 22.7% [22]. Phantoms were prepared in volumes of 20 ml with dilutions resulting in ILp values of [0.455, 0.568, 0.908, 1.135, 1.476, 2.951, 4.540] %. The reduced scattering coefficient of the Intralipid stock solution was determined to be equation m24 (800 nm) = 18.6 mm−1 by using the formula given by Michels et al. [24]. The dependence of equation m25 on the Intralipid concentration was assumed to be linear (for the range of concentrations investigated) and the resulting wavelength-dependent equation m26 is given for selected phantoms in Fig. 2 (A). The dilutions selected for this study resulted in a range of equation m27 = [0.4-8] mm−1 that is representative of the values reported in tissue [37].

Reflectance measurements were performed by submerging the probe tip a few millimeters into the phantom and suspending it multiple centimeters from the container bottom or sides, such that the probe would not collect reflections from container surfaces. The measurement was repeated 3 times for each phantom, with these spectra averaged for each phantom. All measurements were performed with 5 fiber probes of diameter: [0.2, 0.4, 0.6, 0.8, 1.0] mm. The phantom was gently agitated before collecting each spectrum to ensure homogeneity. Measurements were reproducible, with variation between repeated measurements of < 3%.

2.3. Monte Carlo simulations of reflectance measurements

The Monte Carlo code utilized in this study is based on the steady state light transport in multi-layered structures (MCML) approach to stochastically simulate photon propagation within a turbid medium [35,38,39]. The code was adapted to simulate the measurement of reflectance intensity from a turbid medium with homogeneously distributed optical properties; this model of photon propagation has previously been described in detail and will be reviewed here briefly. Photons were initialized by selecting a location on the fiber face in contact with the medium and were launched into a direction within the fiber cone of acceptance, given as equation m28, with the fiber numerical aperture (NA) and the index of refraction of the medium given by n. Both location and direction were sampled from uniform distributions. The fiber probes were modeled in contact with the turbid medium at an air interface, with n of the medium set to 1.33 and that of the fiber to 1.45. The fiber NA was specified as 0.22. Reflection and refraction at the boundary between fiber face and medium were accounted for by using Fresnel’s equation and Snell’s law. The photons were propagated through the medium in discrete steps selected from an exponential distribution that was weighted by the scattering coefficient. Each scattering event resulted in a change in propagation direction that was sampled from a discrete user-specified PF distribution. Photons that contacted the detector fiber within the cone of acceptance were collected. For all simulations, photons that propagated far from the fiber face did not contribute to the collected reflectance intensity and were terminated at a hemispherical limit from the source fiber at a distance great enough to not influence the model outputs for the range of optical properties investigated in this study. The collection area for the single fiber simulations overlapped the source fiber location, with equation m29 specified in the range [0.2, 0.4, 0.6, 0.8, 1.0] mm. For separate source detector devices, source and detector fiber diameters were fixed to 0.2 mm, and photons were collected using annular areas at appropriate radial distance such that the core-to-core separation ρ was specified in the range [0.2, 0.4, 0.6, 1.0, 2.0] mm. Simulations returned collected reflectance intensity in the absence of absorption, and also returned the weighted average number of scattering events experienced by collected photons from each measurement. Reflectance intensity returned for separate source detector simulations was calculated by the ratio of the total number of photons collected (TPC) by the total number of photons launched (TPL) and corrected for differences in collection area (A), as

equation m30

The optical properties of the simulated medium were selected to mimic IL phantoms. The reduced scattering coefficient was specified by linearly scaling the equation m31 of the IL 20% stock solution, given in [24] for the selected IL dilutions. Simulations of single fiber reflectance probes were performed to match the measured optical phantoms, and simulations for the separate source detector device were performed for IL dilutions resulting in ILp values of [0.454, 0.568, 1.135, 2.270, 3.405, 4.540, 5.675]. Wavelength-dependent PFs were specified using empirical formulas also reported in [24], with simulations performed at discrete wavelengths between 400 and 800 nm at 100 nm increments. Each simulation initialized between 10 and 100 million photons.

3. Results

3.1. Reflectance measurement of IL phantoms with a single fiber device

Figure 3 shows a comparison of reflectance spectra measured in IL optical phantoms by reflectance devices that use a single fiber optic to deliver and collect light. Panels (A), (B) and (C) contain spectra measured in a dilution series of IL phantoms, with ILp in the range [0.454–4.54]%, with each panel showing spectra obtained by a different fiber diameter from the series [0.2, 0.4, 1.0] mm. Initial inspection of the data show expected trends of increased equation m32 in response to increases in ILp, which is observed by comparison of spectra within each panel. equation m33 also increases with increases in equation m34, which is observed by comparison of spectra between different panels. A closer inspection of the spectra shows that the wavelength-dependent spectral shape of each curve changes with respect to changes in ILp and equation m35. Figure 3 (D) and (E) show normalized equation m36spectra measured by 3 different fiber diameters in selected ILp phantoms of 4.54% and 0.454%, respectively. Spectra are normalized to the intensity at 400 nm to highlight the substantial differences in spectral shapes returned by measurements with different fiber diameters. First we consider the spectra measured by equation m37 = 1 mm in the ILp = 4.54% phantom, whereequation m38 monotonically decreases vs. wavelength. Reflectance spectra measured by smaller fibers in the same phantom show a transition to a more complicated biphasic behavior. Specifically for equation m39 = 0.2 mm, equation m40 monotonically decreases over 400-500 nm, but exhibits a plateauing effect over 500-800 nm. For the lower scattering phantom, with ILp = 0.454%, reflectance spectra exhibit the biphasic behavior for equation m41 = [0.2, 0.4, 1.0] mm, with equation m42 decreasing over 400-500 nm followed by equation m43 increasing from 500 to 800 nm. The data in panels (D) and (E) show that the wavelength-dependent changes in equation m44 measured in the same phantom are substantially different between different fiber diameters. For example, reflectance spectra measured in the 4.54% phantom show a decrease in equation m45from 400 to 800 nm of 25% for equation m46 = 0.2 mm, and 45% for equation m47 = 1.0 mm. Moreover, the equation m48data show substantially different spectral responses for measurements of phantoms with different IL dilutions with the same fiber diameter. For example, reflectance spectra measured with a fiber with equation m49 = 0.4 mm in the 4.54 and 0.454% ILp phantoms show of decrease of 40% and an increase of 5%, respectively.

Fig. 3
Analysis of spectral scattering features observed in experimentally measured reflectance collected by single fiber reflectance device. Left hand column presents calibrated reflectance spectra as measured with single fiber diameters of A) 0.2 (x marks), ...

Understanding the mechanisms underlying these spectral differences requires consideration of two contributing factors: 1) how scattering influences light propagation during measurement and, 2) how these influences change over the wavelength range sampled in the IL phantoms. It is understood that reflectance devices sample photons that have experienced a range of scattering events [27]. For reflectance measured near the source, decreases in equation m50or equation m51are likely to decrease the number of scattering events experienced by collected photons [40], which in turn increases the sensitivity of the collected reflectance intensity to PF-effects. In IL the scattering properties that influence light transport (equation m52and PF) change substantially across the UV-vis wavelength range. Figure 2 (A) shows equation m53 over 400-800 nm. The PF also changes over the same wavelength range, with equation m54 showing that the likelihood of backscatter events increases vs. wavelength. To characterize the link between these wavelength-dependent optical properties and the observed reflectance, we considered the scattering histories of photons collected across the sampled wavelength range. For this purpose, a Monte Carlo model was used to return the average number of scattering events for photons collected at discrete wavelengths along each spectrum reported Fig. 3 (D) and (E); the corresponding scattering event data are shown in panels (F) and (G), respectively. These scattering data are plotted on a log-scale and include annotations to highlight two important thresholds: 1) the single scattering baseline, and 2) an approximate limit to the diffuse regime of 10 scattering events. We first consider the reflectance spectrum that exhibited single-phase behavior, (equation m55 = 1 mm measurement of the 4.54% ILp phantom). Light collected across this spectrum was, on average, multiply scattered (with >10). The resulting diffuse signal is expected to be insensitive to PF-effects, a claim that is supported by agreement in wavelength-dependent trends of equation m56and equation m57 over the sampled wavelength range. Inspection of scattering data for reflectance curves that exhibit the convex biphasic behavior shows that the collected photons often experienced few (<10) scattering events. These measurements contained contributions from photons sampled in the semi-diffuse regime and are expected to be sensitive to the wavelength-dependent changes in PF. This claim is supported by considering that the observed increase in equation m58vs. wavelength in the 500–800 nm region for these spectra is correlated with the wavelength-dependent changes in backscatter probability described by equation m59 and PF(equation m60) as observed in Fig. 2. These data indicate that the spectral features observed in equation m61spectra measured in IL phantoms are attributable to the sensitivity of the reflectance device to the combined wavelength dependent changes in equation m62and PF in IL. While the combined influence of equation m63 and PF on equation m64at a single wavelength has been modeled in detail previously [34,41]; these data present a qualitative representation of how wavelength-dependent features in reflectance spectra are affected by the wavelength dependent scattering properties and device geometry.

The spectral features observed in single fiber reflectance spectra measured in IL phantoms highlight the changing influence of PF across the wavelengths of a single spectrum; these changes can have profound effects on the relationship between equation m65and equation m66. The left column of Fig. 4 shows equation m67vs. equation m68 measured in all IL phantoms with panels (A), (B), and (C), showing measurements by equation m69 = [0.2, 0.4, 1.0] mm, respectively. Here, colors are used to indicate wavelength. These data show that equation m70 cannot be used to uniquely identify equation m71, because PF-effects influence the proportionality between equation m72 and equation m73 differently at each wavelength. This statement is evidenced by the data in Figs. 3 (D) and 3(E), where equation m74 does not explicitly follow equation m75, but is also influenced by equation m76. Comparison of measurements from different fiber diameters showed an underlying relationship between equation m77 and the dimensionless reduced scattering coefficient, given as equation m78. Figure 4 panel (D) contains measurements from all investigated fiber diameters, equation m79 [0.2, 0.4, 0.6, 0.8, 1.0] mm, in investigated phantoms with ILp in the range of [0.454–4.54]%. Inspection of the equation m80 data show that for increasing equation m81 values, less stratification is observed between different wavelengths, meaning that the collected intensities are less influenced by PF. This is consistent with previous work that identified a threshold limit of PF-independent behavior for values of equation m82>10 [35]. This plot shows that reflectance intensity collected by a single fiber optic is interchangeably influenced by equation m83 andequation m84 and highlights the influence of the wavelength-dependent PF on equation m85

Fig. 4
Relationship between experimentally measured single fiber reflectance intensity and scattering properties in IL phantoms. Left hand column shows reflectance intensities vs. reduced scattering coefficient for fiber diameters of A) 0.2, B) 0.4, and C) 1.0 ...

3.2. Reflectance measurement of IL phantoms with separate source-detector fibers

Figure 5 shows simulated reflectance measurements in IL optical phantoms performed with separate source and detector fibers. Here, panels (A), (B), and (C) contain spectra measured in a dilution series of IL phantoms, in the range [0.454–5.675]%, with each panel showing data obtained by a source-detector separation from the series ρ = [0.2, 0.4, 2.0] mm. Initial inspection of these data show a complicated dependence of equation m86 on the phantom ILp. For adjacent source-detector fibers (ρ = 0.2 mm), an increase in ILp results in an increased reflectance intensity across the entire spectrum. This is not the case for source detectors with large separation (ρ = 2.0 mm), as increases in ILp from 0.4 to 1.0% result in increased equation m87, but increases in ILp from 2 to 5% result in decreased equation m88; this nonlinear dependence of equation m89 on equation m90 is described by diffusion theory and has been characterized in detail previously [18]. In order to understand the influence of scattering properties on the spectral response, we first compare spectra measured in the same phantom with different ρ. Figure 5 panels (D) and (E) show reflectance spectra measured at 3 different source detector separations in phantoms with ILp of 4.54% and 0.454%. The panels (F) and (G) show the average number of scattering events corresponding to discrete wavelengths in the spectra plotted on panels (D) and (E), respectively. The normalized spectra from the 4.54% ILp phantom show decreases in equation m91 vs. wavelength for the small fiber separations (ρ = [0.2, 0.4] mm), but an increase of equation m92 vs. wavelength for the large ρ (2.0 mm). The latter spectrum can be described by increases in equation m93 vs. wavelength increasingly attenuating the collected reflectance intensity; therefore equation m94 follows the inverse trend of equation m95 vs. wavelength over the sampled wavelength range. Inspection of the scattering data for these spectra in panel (F) show that these measurements sample diffuse light, and are independent of PF. Figure 5 panel (E) shows normalized reflectance spectra measured in the 0.454% ILp phantom at different ρ. These data show a transition from equation m96 monotonically decreasing vs. wavelength (for ρ = 2.0 mm) to biphasic behavior (for ρ = 0.2). Inspection of the number of scattering events experienced by collected photons in panel (G), shows a transition from sampling multiply scattered light in the diffuse regime (for ρ = 2.0 mm) to the sampling of fewer scattering events in the semi-diffuse regime (for ρ = 0.2 mm). These observations are consistent with the concept that measurements close to the source are more likely to be sensitive to the frequency of large angle backscatter events, and therefore the collected spectrum is influenced by PF-effects.

Fig. 5
Analysis of spectral scattering features observed in simulated reflectance spectra collected by devices with separate source and detector. Left hand column shows dimensionless reflectance spectra returned using selected source detector separations of ...

Figure 6 shows analysis of the relationship between equation m97 and equation m98 for separate source detector devices. Panels (A), (B), and (C) show equation m99vs. equation m100 for measurements of phantoms with a wide range of ILp, with each panel showing measurements from a device with ρ = [0.2, 0.4, 2.0] mm. equation m101measurements close to the source (ρ = 0.2 and 0.4 mm) show a sensitivity to PF effects for small equation m102 values; this is characterized by the non-unique relationship between equation m103and equation m104 in the low scattering region. For ρ = 0.2 and 0.4 mm devices, as equation m105 increases the equation m106shows less stratification due to PF-effects and the data transition to PF-insensitivity for large equation m107 values; this transition follows the shift from sampling predominantly in the semi- diffuse regime to the diffuse regime as scattering increases. The data for the large source detector separation of ρ = 2.0 mm shows a distinct nonlinear dependence on equation m108, with the equation m109vs. equation m110 showing a monotonically increasing region that transitions through a peak into a monotonically decreasing region; this transition occurs for measurements sampling diffuse light that are independent of PF effects. Inspection of the reflectance data measured at different ρ shows an underlying relationship that is identified by conversion of the x-axis to the dimensionless reduced scattering coefficient, here given as equation m111, and conversion of the y-axis to equation m112. Figure 6 panel (D) contains measurements from all investigated source-detector separationsρ = [0.2, 0.4, 0.6, 1.0, 2.0] mm in all simulated IL phantoms ILp = [0.454-5.675] %. These data show that equation m113 is influenced similarly by equation m114 for different ρ.

Fig. 6
Relationship between reflectance intensity returned by simulations of separate source-detector geometries and scattering properties in IL phantoms. Left hand column shows reflectance intensities vs. reduced scattering coefficient for source-detector separations ...

Inspection of these data showed that for equation m115, equation m116showed less than 10% variation about the mean value, and are therefore considered independent of PF influences at different wavelengths. This threshold was used to identify the smallest equation m117 that could be sampled independent of PF for each ρ. For separate source detectors of ρ = [0.2, 0.4, 2.0] mm, the lower bound of equation m118 was estimated as [2.5, 1.25, 0.25] mm−1, respectively. These thresholds are included as vertical dashed lines on plots in Fig. 6 to allow visualization of scattering regions influenced by PF and independent of PF.

4. Discussion

IL optical phantoms have substantial wavelength-dependent changes in equation m119 and PF in the UV-vis region. This study characterizes how these scattering properties influence spectral features of reflectance measured close to the source. Experiments and simulations were used to investigate reflectance spectra returned in a set of IL phantoms that contained a wide range of equation m120 values. Reflectance spectra were measured or simulated for device geometries that delivered and collected with either a single fiber or separate fibers. Data show that changes to the observed spectral remission characteristics occur in response to slight changes in either measurement geometry or IL dilution. These changes are related to the light transport regime sampled, with spectral features showing a (non-linear) dependence on equation m121 for measurements that collect diffusely scattered light, and on top of that a dependence on PF for measurements that collect semi-diffuse light.

Tissue simulating optical phantoms are commonly used for two purposes: 1) characterization of observed reflectance response to controlled optical properties within the phantom, and 2) calibration of the wavelength-specific reflectance device response to allow system independent interpretation of measurements made in samples with uncontrolled optical properties. The results presented in this study have implications for both of these uses. Device characterization studies may seek to establish or validate a relationship between observed reflectance and equation m122 [1618]. However, the data presented in this study show that scattering PF can profoundly influence observed relationships between reflectance intensity and equation m123. The data in Figs. 4 and and66 show that for reflectance measurements that sample dimensionless reduced scattering coefficients below empirically identified threshold values (equation m124 for single fiber devices andequation m125 for separate source detector devices) the collected reflectance intensities are sensitive to PF-effects. Reflectance measurements performed below these thresholds collect a substantial component of light that has undergone few, or single, scattering events, making the probability of collecting light in that regime (at least somewhat) dependent on the backscattering features of the scattering PF. These data indicate that a proportionality established between reflectance intensity and equation m126 in IL at a single wavelength will not hold at other wavelengths. Moreover, proportionalities are not directly extensible to another scattering medium that may contain an unknown wavelength-dependent PF, such as tissue.

IL optical phantoms are also commonly used for system calibration purposes [10]. While this approach can be valid, care must be taken to consider influences of the PF. As shown in Figs. 3 and and5,5, spectral remission features can depend strongly on the IL dilution of the phantom. It is possible to account for these features by characterizing the response of the investigated device to the moments of the PF. This can be investigated either by experimental measurements in another scattering medium with a controlled wavelength-dependent PF, such as a phantom constructed with polystyrene beads [42], or simulation, where variation of PF can be completely controlled. The experimental data for single fiber reflectance spectra in this study were calibrated with measurements in IL, by using a Monte Carlo analysis to convert the measured reflectance spectra to an absolute scale (in units of % of incident photons) [35], a conversion that accounts for the wavelength-dependent PF influence. Such a calibration becomes independent of ILp. This approach to spectral analysis requires sufficient information to estimate PF effects; such information could be provided by spectral information measured at an additional length scale such as with multiple fiber diameters [34,35,41] or at multiple distances from the source [14,33], or spectral information measured with an angular resolution [43,44]. The use of IL as a calibration phantom without correction for PF-effects may introduce a substantial bias into quantitative estimates of scattering parameters, especially if applied to measurements in a sample of unknown optical properties, such as tissue.

This study identified an interesting difference in the relationship between reflectance intensity and equation m127 for device geometries that utilize a single fiber as a source and detector compared with a device that utilizes separate source and detector fibers. Reflectance intensity collected by single fiber devices shows a monotonically increasing response to increases in equation m128, as displayed in Fig. 4. In contrast, reflectance intensity collected by separate source detectors shows a nonlinear response that increases or decreases based on equation m129, as shown in Fig. 6. These limiting behaviors can be intuitively explained. The single fiber geometry allows collection of light that has scattered an infinitesimally small distance into the medium prior to collection. Therefore, in the limit of highly scattering media, the single fiber reflectance device is no less likely to collect light in response to increased equation m130. In contrast, the collection geometry with separate source and detector fibers has a geometrical minimum path that light must travel within the medium prior to collection. Therefore, in the limit of highly scattering media, the collected signal will be attenuated in response to increases in equation m131. These differences in the required light path travelled by light collected by either device leads to different threshold sensitivities to PF, which empirically are equation m132 andequation m133 for each device-type, respectively. This means that for single fiber measurements, the transport length, defined by the inverse of the reduced scattering coefficient, must be an order of magnitude greater than the collection length to be independent of PF effects. A device with separate source detectors is (in general) less sensitive to PF, with collection at locations up to half of the transport length showing insensitivity to PF. This threshold is expected to be dependent on the source and detector fiber diameters, which were fixed at 0.2 mm in the current study; the exact value of the threshold may vary for different fibers sizes. It should be noted that in IL phantoms constructed to mimic the range of equation m134 expected in tissue, in the range of [0.5-2.0] mm−1 [37], devices that collect with ρ as large as 1 mm can be influenced by PF effects. Moreover, inspection of the data in Figs. 5 and and66 show that measurements made in such phantoms with adjacent source and detector fibers (equation m135 mm) may contain a spectral signature that shows little to no wavelength-dependence on equation m136.

While this study specifically investigates the spectral response of specific reflectance geometries, the findings are important for other measurement devices as well. This study characterized the response over a set range of geometrical perturbations to single fiber size and source-detector distance. This study does not perform a comprehensive evaluation of all possible geometrical configurations of fiber probes (e.g. all possible source and detector fiber sizes, inter-fiber distances, or numerical apertures). The results of this study are applicable to reflectance spectra measured by any probe configuration: the sampled transport regime will determine sensitivity to wavelength-dependent changes in PF. The results are also important for devices that sample reflectance at a single wavelength. In this case, care must be taken to correctly interpret the collected signal to avoid mischaracterizing the observed reflectance intensity as influenced exclusively by equation m137. For example, the use of raw reflectance intensity at a single wavelength to correct fluorescence for the influence of optical properties will be problematic considering that equation m138 and PF will influence reflectance and fluorescence intensities differently [9]. Ongoing work will characterize the PF-specific influences of fluorescence correction algorithms.

5 Conclusions

This study investigates the influence of the wavelength-dependent scattering properties equation m139 and PF on spectral features of reflectance collected near the source location in IL optical phantoms. Experiments and simulations are used to characterize features in the 400-800 nm wavelength range from IL optical phantoms with ILp in the range 0.454-5.675%, which relates to a wide range of equation m140 = [0.3–8] mm−1. Data were measured or simulated for reflectance geometries that 1) utilized a single fiber optic for both source and detector with the equation m141 = [0.2–1.0] mm, and 2) separate source and detector fibers with core to core separations in the range ρ = [0.2–2.0] mm. Results for both collection configurations showed a dependence of spectral scattering features on the transport regime sampled during measurement. PF-dependent effects were observed in reflectance spectra measured below a dimensionless threshold of equation m142 (i.e. ten transport lengths) for single fiber devices andequation m143 (i.e. half a transport length) for separate source detector devices. The results show the need to consider the influence of PF when using measurements of IL to either characterize or calibrate a reflectance device that collects light close to the source.


This work has been funded by NIH grants PO1 CA084203 and PO1 CA080139.

References and Links

1. Amelink A., Kaspers O. P., Sterenborg H. J. C. M., van der Wal J. E., Roodenburg J. L. N., Witjes M. J. H., “Non-invasive measurement of the morphology and physiology of oral mucosa by use of optical spectroscopy,” Oral Oncol. 44(1), 65–71 (2008). doi: 10.1016/j.oraloncology.2006.12.011. [PubMed] [Cross Ref]
2. Brown J. Q., Vishwanath K., Palmer G. M., Ramanujam N., “Advances in quantitative UV-visible spectroscopy for clinical and pre-clinical application in cancer,” Curr. Opin. Biotechnol. 20(1), 119–131 (2009). doi: 10.1016/j.copbio.2009.02.004. [PMC free article] [PubMed] [Cross Ref]
3. Mourant J. R., Hielscher A. H., Eick A. A., Johnson T. M., Freyer J. P., “Evidence of intrinsic differences in the light scattering properties of tumorigenic and nontumorigenic cells,” Cancer 84(6), 366–374 (1998). doi: 10.1002/(SICI)1097-0142(19981225)84:6<366::AID-CNCR9>3.0.CO;2-7. [PubMed] [Cross Ref]
4. Sokolov K., Follen M., Richards-Kortum R., “Optical spectroscopy for detection of neoplasia,” Curr. Opin. Chem. Biol. 6(5), 651–658 (2002). doi: 10.1016/S1367-5931(02)00381-2. [PubMed] [Cross Ref]
5. Bigio I. J., Bown S. G., “Spectroscopic sensing of cancer and cancer therapy: current status of translational research,” Cancer Biol. Ther. 3(3), 259–267 (2004). doi: 10.4161/cbt.3.3.694. [PubMed] [Cross Ref]
6. Garcia-Allende P. B., Krishnaswamy V., Hoopes P. J., Samkoe K. S., Conde O. M., Pogue B. W., “Automated identification of tumor microscopic morphology based on macroscopically measured scatter signatures,” J. Biomed. Opt. 14(3), 034034 (2009). doi: 10.1117/1.3155512. [PMC free article] [PubMed] [Cross Ref]
7. Weersink R., Patterson M. S., Diamond K., Silver S., Padgett N., “Noninvasive measurement of fluorophore concentration in turbid media with a simple fluorescence /reflectance ratio technique,” Appl. Opt. 40(34), 6389–6395 (2001). doi: 10.1364/AO.40.006389. [PubMed] [Cross Ref]
8. Müller M. G., Georgakoudi I., Zhang Q., Wu J., Feld M. S., “Intrinsic fluorescence spectroscopy in turbid media: disentangling effects of scattering and absorption,” Appl. Opt. 40(25), 4633–4646 (2001). doi: 10.1364/AO.40.004633. [PubMed] [Cross Ref]
9. Kanick S. C., Robinson D. J., Sterenborg H. J. C. M., Amelink A., “Semi-empirical model of the effect of scattering on single fiber fluorescence intensity measured on a turbid medium,” Biomed. Opt. Express 3(1), 137–152 (2012). doi: 10.1364/BOE.3.000137. [PMC free article] [PubMed] [Cross Ref]
10. Pogue B. W., Patterson M. S., “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt. 11(4), 041102 (2006). doi: 10.1117/1.2335429. [PubMed] [Cross Ref]
11. Boas D. A., Dale A. M., Franceschini M. A., “Diffuse optical imaging of brain activation: approaches to optimizing image sensitivity, resolution, and accuracy,” Neuroimage 23(Suppl 1), S275–S288 (2004). doi: 10.1016/j.neuroimage.2004.07.011. [PubMed] [Cross Ref]
12. McBride T. O., Pogue B. W., Gerety E. D., Poplack S. B., Osterberg U. L., Paulsen K. D., “Spectroscopic diffuse optical tomography for the quantitative assessment of hemoglobin concentration and oxygen saturation in breast tissue,” Appl. Opt. 38(25), 5480–5490 (1999). doi: 10.1364/AO.38.005480. [PubMed] [Cross Ref]
13. Mourant J. R., Fuselier T., Boyer J., Johnson T. M., Bigio I. J., “Predictions and measurements of scattering and absorption over broad wavelength ranges in tissue phantoms,” Appl. Opt. 36(4), 949–957 (1997). doi: 10.1364/AO.36.000949. [PubMed] [Cross Ref]
14. Thueler P., Charvet I., Bevilacqua F., St. Ghislain M., Ory G., Marquet P., Meda P., Vermeulen B., Depeursinge C., “In vivo endoscopic tissue diagnostics based on spectroscopic absorption, scattering, and phase function properties,” J. Biomed. Opt. 8(3), 495–503 (2003). doi: 10.1117/1.1578494. [PubMed] [Cross Ref]
15. Amelink A., Sterenborg H. J. C. M., “Measurement of the local optical properties of turbid media by differential path-length spectroscopy,” Appl. Opt. 43(15), 3048–3054 (2004). doi: 10.1364/AO.43.003048. [PubMed] [Cross Ref]
16. Johns M., Giller C. A., German D. C., Liu H., “Determination of reduced scattering coefficient of biological tissue from a needle-like probe,” Opt. Express 13(13), 4828–4842 (2005). doi: 10.1364/OPEX.13.004828. [PubMed] [Cross Ref]
17. Reif R., A’Amar O., Bigio I. J., “Analytical model of light reflectance for extraction of the optical properties in small volumes of turbid media,” Appl. Opt. 46(29), 7317–7328 (2007). doi: 10.1364/AO.46.007317. [PubMed] [Cross Ref]
18. Kim A., Roy M., Dadani F., Wilson B. C., “A fiberoptic reflectance probe with multiple source-collector separations to increase the dynamic range of derived tissue optical absorption and scattering coefficients,” Opt. Express 18(6), 5580–5594 (2010). doi: 10.1364/OE.18.005580. [PubMed] [Cross Ref]
19. Krishnaswamy V., Laughney A. M., Paulsen K. D., Pogue B. W., “Dark-field scanning in situ spectroscopy platform for broadband imaging of resected tissue,” Opt. Lett. 36(10), 1911–1913 (2011). doi: 10.1364/OL.36.001911. [PubMed] [Cross Ref]
20. Ninni P. D., Martelli F., Zaccanti G., “Intralipid: towards a diffusive reference standard for optical tissue phantoms,” Phys. Med. Biol. 56(2), N21–N28 (2011). doi: 10.1088/0031-9155/56/2/N01. [PubMed] [Cross Ref]
21. van Staveren H. J., Moes C. J., van Marie J., Prahl S. A., van Gemert M. J., “Light scattering in Intralipid-10% in the wavelength range of 400-1100 nm,” Appl. Opt. 30(31), 4507–4514 (1991). doi: 10.1364/AO.30.004507. [PubMed] [Cross Ref]
22. Zaccanti G., Del Bianco S., Martelli F., “Measurements of optical properties of high-density media,” Appl. Opt. 42(19), 4023–4030 (2003). doi: 10.1364/AO.42.004023. [PubMed] [Cross Ref]
23. Martelli F., Zaccanti G., “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15(2), 486–500 (2007). doi: 10.1364/OE.15.000486. [PubMed] [Cross Ref]
24. Michels R., Foschum F., Kienle A., “Optical properties of fat emulsions,” Opt. Express 16(8), 5907–5925 (2008). doi: 10.1364/OE.16.005907. [PubMed] [Cross Ref]
25. Farrell T. J., Patterson M. S., Wilson B., “A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo,” Med. Phys. 19(4), 879–888 (1992). doi: 10.1118/1.596777. [PubMed] [Cross Ref]
26. Jacques S. L., Pogue B. W., “Tutorial on diffuse light transport,” J. Biomed. Opt. 13(4), 041302 (2008). doi: 10.1117/1.2967535. [PubMed] [Cross Ref]
27. Mourant J. R., Boyer J., Hielscher A. H., Bigio I. J., “Influence of the scattering phase function on light transport measurements in turbid media performed with small source-detector separations,” Opt. Lett. 21(7), 546–548 (1996). doi: 10.1364/OL.21.000546. [PubMed] [Cross Ref]
28. Bevilacqua F., Depeursinge C., “Monte Carlo study of diffuse reflectance at source–detector separations close to one transport mean free path,” J. Opt. Soc. Am. A 16(12), 2935–2945 (1999). doi: 10.1364/JOSAA.16.002935. [Cross Ref]
29. Kienle A., Forster F. K., Hibst R., “Influence of the phase function on determination of the optical properties of biological tissue by spatially resolved reflectance,” Opt. Lett. 26(20), 1571–1573 (2001). doi: 10.1364/OL.26.001571. [PubMed] [Cross Ref]
30. Saidi I. S., Jacques S. L., Tittel F. K., “Mie and Rayleigh modeling of visible-light scattering in neonatal skin,” Appl. Opt. 34(31), 7410–7418 (1995). doi: 10.1364/AO.34.007410. [PubMed] [Cross Ref]
31. Bargo P. R., Prahl S. A., Goodell T. T., Sleven R. A., Koval G., Blair G., Jacques S. L., “In vivo determination of optical properties of normal and tumor tissue with white light reflectance and an empirical light transport model during endoscopy,” J. Biomed. Opt. 10(3), 034018 (2005). doi: 10.1117/1.1921907. [PubMed] [Cross Ref]
32. Hull E. L., Foster T. H., “Steady-state reflectance spectroscopy in the P3 approximation,” J. Opt. Soc. Am. A 18(3), 584–599 (2001). doi: 10.1364/JOSAA.18.000584. [Cross Ref]
33. Bevilacqua F., Piguet D., Marquet P., Gross J. D., Tromberg B. J., Depeursinge C., “In vivo local determination of tissue optical properties: applications to human brain,” Appl. Opt. 38(22), 4939–4950 (1999). doi: 10.1364/AO.38.004939. [PubMed] [Cross Ref]
34. Kanick S. C., Gamm U. A., Sterenborg H. J. C. M., Robinson D. J., Amelink A., “Method to quantitatively estimate wavelength-dependent scattering properties from multidiameter single fiber reflectance spectra measured in a turbid medium,” Opt. Lett. 36(15), 2997–2999 (2011). doi: 10.1364/OL.36.002997. [PubMed] [Cross Ref]
35. Kanick S. C., Gamm U. A., Schouten M., Sterenborg H. J. C. M., Robinson D. J., Amelink A., “Measurement of the reduced scattering coefficient of turbid media using single fiber reflectance spectroscopy: fiber diameter and phase function dependence,” Biomed. Opt. Express 2(6), 1687–1702 (2011). doi: 10.1364/BOE.2.001687. [PMC free article] [PubMed] [Cross Ref]
36. Kanick S. C., Sterenborg H. J. C. M., Amelink A., “Empirical model of the photon path length for a single fiber reflectance spectroscopy device,” Opt. Express 17(2), 860–871 (2009). doi: 10.1364/OE.17.000860. [PubMed] [Cross Ref]
37. Cheong W. F., Prahl S. A., Welch A. J., “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26(12), 2166–2185 (1990). doi: 10.1109/3.64354. [Cross Ref]
38. Wang L., Jacques S. L., Zheng L., “MCML--Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed. 47(2), 131–146 (1995). doi: 10.1016/0169-2607(95)01640-F. [PubMed] [Cross Ref]
39. Kanick S. C., Robinson D. J., Sterenborg H. J. C. M., Amelink A., “Monte Carlo analysis of single fiber reflectance spectroscopy: photon path length and sampling depth,” Phys. Med. Biol. 54(22), 6991–7008 (2009). doi: 10.1088/0031-9155/54/22/016. [PubMed] [Cross Ref]
40. Pogue B. W., Burke G., “Fiber-optic bundle design for quantitative fluorescence measurement from tissue,” Appl. Opt. 37(31), 7429–7436 (1998). doi: 10.1364/AO.37.007429. [PubMed] [Cross Ref]
41. Gamm U. A., Kanick S. C., Sterenborg H. J. C. M., Robinson D. J., Amelink A., “Measurement of tissue scattering properties using multi-diameter single fiber reflectance spectroscopy: in silico sensitivity analysis,” Biomed. Opt. Express 2(11), 3150–3166 (2011). doi: 10.1364/BOE.2.003150. [PMC free article] [PubMed] [Cross Ref]
42. Gélébart B., Tinet E., Tualle J., Avrillier S., “Phase function simulation in tissue phantoms: a fractal approach,” Pure Appl. Opt. 5(4), 377–388 (1996). doi: 10.1088/0963-9659/5/4/005. [Cross Ref]
43. Pyhtila J., Graf R., Wax A., “Determining nuclear morphology using an improved angle-resolved low coherence interferometry system,” Opt. Express 11(25), 3473–3484 (2003). doi: 10.1364/OE.11.003473. [PubMed] [Cross Ref]
44. Kim Y. L., Liu Y., Wali R. K., Roy H. K., Goldberg M. J., Kromin A. K., Chen K., Backman V., “Simultaneous measurement of angular and spectral properties of light scattering for characterization of tissue microarchitecture and its alteration in early precancer,” IEEE J. Sel. Top. Quantum Electron. 9(2), 243–256 (2003). doi: 10.1109/JSTQE.2003.814183. [Cross Ref]

Articles from Biomedical Optics Express are provided here courtesy of Optical Society of America