Effect of dependent scattering on the optical properties of Intralipid tissue phantoms

doi:10.1364/BOE.2.002265

« Show journal navigation

Effect of dependent scattering on the optical properties of Intralipid tissue phantoms

Paola Di Ninni, Fabrizio Martelli, and Giovanni Zaccanti »View Author Affiliations

Biomedical Optics Express, Vol. 2, Issue 8, pp. 2265-2278 (2011)
http://dx.doi.org/10.1364/BOE.2.002265

View Full Text Article

Acrobat PDF (1210 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Materials and Methods
3. Experimental Results
4. Discussion and Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (28)
Cited By
Figures (3)
Metrics

Abstract

The calibration of optical tissue-simulating phantoms remains an open question in spite of the many techniques proposed for accurate measurements of optical properties. As a consequence, a reference phantom with well known optical properties is still missing. As a first step towards a reference phantom we have recently proposed to use dilutions of Intralipid 20%. In this paper we discuss a matter that is commonly ignored when dilutions are prepared, i.e., the possibility of deviations from the simple linear relationships between the optical properties of the dilution and the Intralipid concentration due to the effects of dependent scattering. The results of an experimental investigation showed that dependent scattering does not affect absorption. As for the reduced scattering coefficient the effect can be described adding a term proportional to the square of the concentration. However, for concentrations of interest for tissue optics deviations from linearity remain within about 2%. The experimental investigation also showed that the microphysical properties of Intralipid are not affected by dilution. These results show the possibility to easily obtain a liquid diffusive phantom whose optical properties are known with error smaller than about 1%. Due to the intrinsic limitations of the different techniques proposed for measuring the optical properties it seems difficult to obtain a similar accuracy for solid phantoms.

1. Introduction

Many different methods have been presented over the past two decades for measuring the optical properties of diffusive media based on both measurements in the time domain [1

S. Madsen, B. Wilson, M. Patterson, Y. Park, S. Jacques, and Y. Hefetz, “Experimental tests of a simple diffusion model for the estimation of scattering and absorption coefficients of turbid media from time-resolved diffuse reflectance measurements,” Appl. Opt. 31, 3509–3517 (1992). [CrossRef] [PubMed]

–3

E. Alerstam, S. Andersson-Engels, and T. Svensson, “Improved accuracy in time-resolved diffuse reflectance spectroscopy,” Opt. Express 16, 10440–10454 (2008). [CrossRef]

], in the frequency domain [4

S. Fantini, M. Franceschini, and E. Gratton, “Semi-infinite-geometry boundary problem for light migration in highly scattering media: a frequency-domain study in the diffusion approximation,” J. Opt. Soc. Am. B 11, 2128–2138 (1994). [CrossRef]

, 5

H. Xu and M. Patterson, “Determination of the optical properties of tissue-simulating phantoms from interstitial frequency domain measurements of relative fluence and phase difference,” Opt. Express 14, 6485–6501 (2006). [CrossRef]

], and in the CW domain [6

P. Marquet, F. P. Bevilacqua, C. D. Depeursinge, and E. B. de Haller, “Determination of reduced scattering and absorption coefficients by a single charge-coupled-device array measurement, part I: comparison between experiments and simulations,” Opt. Eng. 34, 2055–2063 (1995). [CrossRef]

–9

N. Rajaram, T. H. Nguyen, and J. W. Tunnell, “Lookup table-based inverse model for determining optical properties of turbid media,” J. Biomed. Opt. 13, 050501 (2008). [CrossRef]

]. However, as also shown by recently published results [10

A. Pifferi, A. Torricelli, A. Bassi, P. Taroni, R. Cubeddu, H. Wabnitz, D. Grosenick, M. Mller, R. Macdonald, J. Swartling, T. Svensson, S. Andersson-Engels, R. L. P. van Veen, H. J. C. M. Sterenborg, J.-M. Tualle, H. L. Nghiem, S. Avrillier, M. Whelan, and H. Stamm, “Performance assessment of photon migration instruments: the MEDPHOT protocol,” Appl. Opt. 44, 2104–2114 (2005). [CrossRef]

, 11

J. Bouchard, I. Veilleux, R. Jedidi, I. Noiseux, M. Fortin, and O. Mermut, “Reference optical phantoms for diffuse optical spectroscopy. Part 1 – Error analysis of a time resolved transmittance characterization method,” Opt. Express 18, 11495–11507 (2010). [CrossRef]

], in spite of the many methodologies and instrumentation proposed, the accurate determination of the optical properties of diffusive media remains a difficult task. Ref. [10

] reports the results of a multi-laboratory experiment: The application of the MED-PHOT protocol in measuring the optical properties of the same phantom with eight different instrumentation revealed differences of about 40% for the reduced scattering coefficient (μ′_s ) and 30% for the absorption coefficient (μ_a ). In Ref. [11

], devoted to the development of a solid reference phantom for diffuse optical spectroscopy, refined measurements of time resolved transmittance have been inverted using a lookup table based on Monte Carlo results and the different sources of error have been carefully investigated. From the results presented it seems difficult to understand the ultimate reason for the relatively large error (about 6% for the reduced scattering coefficient and 11% for the absorption coefficient) that remains on the results.

Difficulties in obtaining accurate measurements of optical properties of diffusive media are intrinsic to the relationships between μ_a and μ′_s and the measured quantities. Usually μ_a and μ′_s are obtained from measurements of diffuse reflectance or of diffuse transmittance, and for the difficulties in calibrating the instrumentation measurements are usually available only in arbitrary units. Without the information on the absolute values for both time domain, frequency domain, and multidistance CW measurements the strong correlation between μ_a and μ′_s makes the non-linear inversion procedure very sensitive to experimental errors. The intrinsic difficulties related to the inversion procedures can be overcome if relative measurements of diffused light can be repeated after changing in a well controlled way the absorption or the reduced scattering coefficient of the medium, adding for instance a calibrated absorber, or changing the concentration of scatterers. In this case, between the optical properties to be determined and the relative variations of the measured quantities a linear relationship can often be obtained with a low correlation between μ_a and μ′_s leading to robust and accurate results [12

F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]

, 13

L. Spinelli, F. Martelli, A. Farina, A. Pifferi, A. Torricelli, R. Cubeddu, and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. Time-resolved method,” Opt. Express 15, 6589–6604 (2007). [CrossRef]

As a consequence of the difficulties in measuring the optical properties a reference phantom with well known optical properties is still missing, in spite of the many different materials proposed especially for tissue-simulating phantoms [14

B. W. Pogue and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt. 11, 041102 (2006). [CrossRef]

]. We point out that reference phantoms are necessary tools for studying photon migration, for the development, validation, and calibration of biomedical optical instrumentation [14

B. W. Pogue and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt. 11, 041102 (2006). [CrossRef]

] and, more in general, of devices for quality control of agricultural [15

M. Zude (Ed.), Optical Monitoring of Fresh and Processed Agricultural Crops , (Contemporary Food Engineering Series), (CRC Press, Boca Raton, Florida, 2009).

], pharmaceutical products [16

J. Johansson, S. Folestad, M. Josefson, A. Sparn, C. Abrahamsson, S. Anderson-Engels, and S. Svanberg, “Time-resolved NIR/Vis spectroscopy for analysis of solids: Pharmaceutical tablets,” Appl. Spectrosc. 56, 725–731 (2002). [CrossRef]

], etc. In a recently published paper [17

P. Di Ninni, F. Martelli, and G. Zaccanti, “Intralipid: towards a diffusive reference standard for optical tissue phantoms,” Phys. Med. Biol. 56, N21–N28 (2011). [CrossRef]

] we proposed the use of Intralipid 20% as a first step towards a reference diffusive medium for tissue-simulating phantoms. The proposal originated from measurements of optical properties on many samples from nine different batches produced over a period of ten years that showed a high temporal stability and surprisingly small batch-to-batch variations: At NIR wavelengths, for the reduced scattering coefficient, deviations from the value averaged over the nine batches were within 1.5%, and the results for the absorption coefficient were slightly smaller with respect to absorption of pure water [17

P. Di Ninni, F. Martelli, and G. Zaccanti, “Intralipid: towards a diffusive reference standard for optical tissue phantoms,” Phys. Med. Biol. 56, N21–N28 (2011). [CrossRef]

]. It might be therefore useful to measure as precisely as possible the optical properties of Intralipid. In Ref. [12

F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]

] a method has been presented based on measurements of effective attenuation coefficient (μ_eff ) for different concentrations of Intralipid in water (method of water absorption). Making the assumption that the optical properties of the dilution are related to the concentration by linear relationships, Intralipid can be calibrated exploiting the knowledge of the absorption coefficient of water. The accuracy of the method is ultimately limited by the error on the absorption coefficient of water and by the assumption of the linear relationships between the optical properties and the concentration.

The possibility of deviations from the linear relationships is commonly ignored in evaluating the optical properties of dilutions, but it should be reminded that the linear relationships are expected to be applicable only if 1) the microphysical properties of suspended particles (size distribution, shape, and refractive index) are not affected by dilution in water, and 2) the independent scattering approximation is fulfilled [18

M. I. Mishchenko, L. D. Travis, and A. A. Lacis, Scattering, Absorption, and Emission of Light by Small Particles , (Cambridge University Press, Cambridge, UK, 2002).

, 19

A. Ishimaru and Y. Kuga , “Attenuation constant of a coherent field in a dense distribution of particles,” J. Opt. Soc. Am. 72 1317–20 (1982). [CrossRef]

]. The first assumption ensures that the optical properties of suspended single particles (scattering function, scattering and absorption cross section) are independent of the concentration of diffusive medium. The second assumption ensures that suspended particles act as independent and uncorrelated scatterers, so that the power scattered (absorbed) by the unit volume is the sum of the power scattered (absorbed) by each particle. If both the assumptions are fulfilled the scattering and absorption coefficients due to suspended particles are proportional to the number of particles in the unit volume and thus to the volume concentration of particles, whereas the scattering function is independent of the concentration.

In Ref. [12

F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]

] the effect of dependent scattering on measurements carried out at λ = 751 nm has been discussed making use of experimental results obtained at concentrations significantly higher than concentrations of interest for tissue simulating phantoms, and at a shorter wavelength (λ = 632.8 nm) [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

], whereas the possibility that the microphysical properties of scattering particles depend on the concentration has not been discussed at all. This paper is devoted to a deeper investigation of these effects. Experimental results showed that dilution of Intralipid in water does not appreciably change the microphysical properties of suspended particles, and that dependent scattering does not affect the absorption coefficient. For the reduced scattering coefficient small deviations from linearity have been observed that can be accounted for with a term proportional to the square of the concentration. Anyway, for concentrations necessary to obtain the values of μ′_s of practical interest for tissue phantoms (μ′_s < 2 mm⁻¹) deviations remain within about 2%.

The procedure used to investigate the effect of dependent scattering needs to use an absorber: We used India ink, after having checked that it can be added without affecting the microphysical properties of Intralipid.

The methodology we used for the experimental investigation is described in section 2; the experimental results are presented in section 3 and conclusions are in section 4.

2. Materials and Methods

2.1. Monitoring the Microphysical Properties of Diluted Intralipid

Intralipid consists of fat droplets suspended in water that make it a highly scattering medium [21

R. Michels, F. Foschum, and A. Kienle, “Optical properties of fat emulsions,” Opt. Express 16, 5907–5925 (2008). [CrossRef]

]. To investigate if the microphysical properties of fat droplets change when Intralipid 20% is diluted in water or when an absorber is added, we carried out measurements of extinction coefficient. With reference to spherical droplets, and making the assumption of independent scattering [18

M. I. Mishchenko, L. D. Travis, and A. A. Lacis, Scattering, Absorption, and Emission of Light by Small Particles , (Cambridge University Press, Cambridge, UK, 2002).

], the extinction coefficient μ_e due to Intralipid is related to the scattering and absorption properties of suspended droplets by [22

C. F. Bohren and D. R. Huffman, Absorption and scattering of light by small particles , (John Wiley and Sons, New York, 1983).

]

μ_{e} (ρ_{i l}) = N \int_{0}^{\infty} C_{e i l} (\frac{r}{λ}, n) f (r) d r

(1)

where f (r) is the particle size distribution, N f (r)dr is the number of droplets in the unit volume with radius between r and r + dr, and C_eil is the extinction cross section. C_eil depends on the relative refractive index n of fat droplets with respect to water, and on the ratio between the size of the droplet and the wavelength. Since N is proportional to the volume concentration ρ_il of Intralipid 20% diluted in water, μ_eil can be written as

μ_{e} (ρ_{i l}) = ɛ_{e i l} ρ_{i l}

(2)

where the coefficient ɛ_eil is the specific extinction coefficient of Intralipid 20%, determined by the concentration of fat droplets, by their size distribution, and their refractive index.

Measurements of ɛ_eil can therefore be used to monitor any change in the microphysical properties of fat droplets due to the dilution of Intralipid 20% in water or to the addition of absorbers: If dilution or addition of absorbers affects the size distribution or the refractive index of fat droplets we expect to obtain different values of ɛ_eil from measurements carried out on dilutions prepared in different ways.

The specific extinction coefficient has been obtained from measurements of collimated transmittance on samples with moderate turbidity (μ_e < 0.05 mm⁻¹). Intralipid was diluted inside a scattering cell and the transmitted power P(ρ_il ) has been measured for different concentrations. Inverting the Lambert-Beer law the value of the extinction coefficient has been obtained for each concentration as

μ_{e} (ρ_{i l}) = \frac{1}{L} ln \frac{P (ρ_{i l} = 0)}{P (ρ_{i l})},

(3)

where L is the thickness of the cell, and ɛ_eil has been determined from the slope of the straight line that best fits μ_e (ρ_il ) as a function of ρ_il .

2.2. Dependent Scattering

The independent scattering approximation is applicable provided the average distance between scatterers is sufficiently large with respect to the wavelength [18

M. I. Mishchenko, L. D. Travis, and A. A. Lacis, Scattering, Absorption, and Emission of Light by Small Particles , (Cambridge University Press, Cambridge, UK, 2002).

]. Experimental investigations [19

A. Ishimaru and Y. Kuga , “Attenuation constant of a coherent field in a dense distribution of particles,” J. Opt. Soc. Am. 72 1317–20 (1982). [CrossRef]

, 20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

] showed that for particles with radius smaller than the wavelength the effect of dependent scattering becomes significant when the volume concentration of scatterers is larger than 0.01. Since the volume concentration of fat droplets of Intralipid 20% necessary to obtain the values of μ′_s of interest for tissue phantoms can be larger than 0.01 (as an example, the volume concentration of Intralipid 20% necessary to obtain μ′_s = 2 mm⁻¹ at λ = 751 nm is ρ_il ≅ 0.1, that corresponds to a volume concentration of fat droplets of ≅ 0.022), to obtain phantoms with well calibrated optical properties it is therefore necessary to check if the dependent scattering leads to significant deviations from the linearity.

The procedure we used to investigate the effect of dependent scattering at NIR wavelengths is based on measurements of effective attenuation coefficient

μ_{e f f} = \sqrt{3 μ_{a} {μ'}_{s}}

of dilutions with high turbidity of Intralipid and India ink in water. μ_eff has been obtained from multidistance measurements of fluence rate ϕ(r) carried out in an infinite medium illuminated by an isotropic CW source [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

]. From the solution of the diffusion equation for the fluence at distance r from a point source emitting a unitary power:

ϕ (r) = \frac{3 {μ'}_{s}}{4 π r} exp (- μ_{e f f} r),

(4)

we obtain

ln [r ϕ (r)] = ln \frac{3 {μ'}_{s}}{4 π} - μ_{e f f} r,

(5)

and μ_eff is obtained from the slope of the straight line that best fits ln[rϕ(r)] as a function of the source-receiver distance r. The overall procedure can be summarized in 6 steps.

Step 1: The optical properties of Intralipid 20% are first calibrated using the method of water absorption [12

F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]

], which consists in measuring μ_eff as a function of the concentration of Intralipid diluted in water. With the assumption that the independent scattering approximation is applicable and that dilution does not affect the microphysical properties of fat droplets, the reduced scattering and the absorption coefficients of the dilution are related to ρ_il by

{μ'}_{s} (ρ_{i l}) = {ɛ'}_{s i l} ρ_{i l}

(6)

μ_{a} (ρ_{i l}) = ɛ_{a i l} ρ_{i l} + ɛ_{a H_{2} O} (1 - ρ_{i l})

(7)

and

μ_{e f f}^{2} (ρ_{i l}) = 3 {ɛ'}_{s i l} ɛ_{a H_{2} O ρ_{i l}} + 3 {ɛ'}_{s i l} (ɛ_{a i l} - ɛ_{a H_{2} O}) ρ_{i l}^{2}

(8)

where ɛ′_sil and ɛ_ail are the specific reduced scattering coefficient and the absorption coefficient of Intralipid 20% and ɛ _{aH
₂
O} is the absorption coefficient of pure water. Exploiting the knowledge of the absorption coefficient of water, ɛ′_sil and ɛ_ail are obtained as

{ɛ'}_{s i l} = \frac{I_{i l}}{3 ɛ_{a H_{2} O}}

(9)

ɛ_{a i l} = \frac{S_{i l}}{3 {ɛ'}_{s i l}} + ɛ_{a H_{2} O} = ɛ_{a H_{2} O} (1 + \frac{S_{i l}}{I_{i l}})

(10)

where I_il and S_il are respectively the intercept and the slope of the straight line that best fits

μ_{e f f}^{2} (ρ_{i l}) / ρ_{i l}

as a function of ρ_il . These relations show that to get accurate measurements of ɛ′_sil and ɛ_ail the value of ɛ _{aH
₂
O} must be known with high accuracy.

Step 2: A sample of the Intralipid calibrated at step 1 is used to calibrate the absorption coefficient of India ink. The absorption coefficient is obtained from measurements of effective attenuation coefficient on a dilution of the calibrated Intralipid, as a function of the ink concentration ρ_ink . Denoting with μ′_sil (ρ_il ) and μ_ail (ρ_il ) the reduced scattering and the absorption coefficient of the diluted Intralipid and with ɛ_aink the absorption coefficient of non diluted ink,

μ_{e f f}^{2}

can be written as

μ_{e f f}^{2} (ρ_{i l}, ρ_{i n k}) = 3 {μ'}_{s} (ρ_{i l}) (μ_{a} (ρ_{i l}) + ɛ_{a i n k} ρ_{i n k})

(11)

where we assumed that the added ink does not appreciably change the Intralipid concentration and we neglected the contribution due to the scattering of ink. ɛ_aink is obtained as

ɛ_{a i n k} = \frac{S_{i n k}}{3 {μ'}_{s} (ρ_{i l})}

(12)

where S_ink is the slope of the straight line that best fits

μ_{e f f}^{2} (ρ_{i l}, ρ_{i n k})

as a function of ρ_ink . We point out that if the effect of dependent scattering is significant Eq. (8) is approximated and the error on μ′_s (ρ_il ) affects the value of ɛ_aink .

Step 3: Making use of the India ink calibrated at step 2, the method of adding absorption [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

] is repeatedly applied to measure the optical properties of dilutions with different concentrations of Intralipid. For a given concentration ρ_il the method consists in measuring μ_eff after the addition of small quantities of the calibrated ink that changes the absorption coefficient of the dilution of known quantities Δμ_a (ρ_ink ) = ɛ_ainkρ_ink . Using Eq. (11), μ′_s (ρ_il ) and μ_a (ρ_il ) are obtained as

{μ'}_{s} (ρ_{i l}) = \frac{S_{1 i n k}}{3 ɛ_{a i n k}}

(13)

μ_{a} (ρ_{i l}) = ɛ_{a i n k} \frac{I_{1 i n k}}{S_{1 i n k}}

(14)

where S _1ink and I _1ink are respectively the slope and the intercept of the straight line that best fits

μ_{e f f}^{2} (ρ_{i l}, ρ_{i n k})

as a function of ρ_ink .

Step 4: To study deviations from linearity in Eqs. (6) and (7) due to dependent scattering, the values of μ′_s (ρ_il ) and μ_a (ρ_il ) measured at step 3 are plotted as a function of ρ_il . Since, as it will be shown in Sect. 3.1, the microphysical properties of fat droplets are not affected by the dilution of Intralipid in water and by the addition of India ink, deviations from linearity can be ascribed to dependent scattering. As it will be shown by experimental results of Sect. 3.2, for concentrations of Intralipid 20% of practical interest for tissue phantoms deviations from the linearity of μ′_s as a function of ρ_il are well described by

{μ'}_{s} (ρ_{i l}) = {ɛ'}_{s 1 i l} ρ_{i l} + {ɛ'}_{s 2 i l} ρ_{i l}^{2},

(15)

with ɛ′ _s2il significantly smaller than ɛ′ _s1il, while for μ_a (ρ_il ) no appreciable deviation from the linear behavior (Eq. (7)) has been observed. The coefficients ɛ′ _s1il and ɛ′ _s2il are obtained again with a linear fit from the intercept I _1il and the slope S _1il of the straight line that best fits μ′_s (ρ_il )/ρ_il as a function of ρ_il . In particular the ratio ɛ′ _s2il/ɛ′ _s1il can be obtained as:

\frac{{ɛ'}_{s 2 i l}}{{ɛ'}_{s 1 i l}} = \frac{S_{1 i l}}{I_{1 i l}} .

(16)

We point out that the ratio ɛ′ _s2il/ɛ′ _s1il is not affected by the error on ɛ_aink . In fact, the error on ɛ_aink causes a systematic overestimation or underestimation of μ′_s (ρ_il ) by a constant factor that does not affect the ratio ɛ′ _s2il/ɛ′ _s1il.

From the intercept I ₂ _il and the slope S ₂ _il of the straight line that best fits μ_a (ρ_il ) as a function of ρ_il , according to Eq. (7) we obtain: ɛ _{aH
₂
O} = I _2il and ɛ_ail = S _2il + I _2il. We notice that these values are affected by the error on the absorption coefficient of India ink calibrated at step 2.

Step 5: The method of water absorption (step 1) is applied again, but using for μ′_s (ρ_il ) the expression that includes the effect of dependent scattering, i.e., using Eq. (15) instead of Eq. (6). Equation (8) changes to:

μ_{e f f}^{2} (ρ_{i l}) = 3 {ɛ'}_{s 1 i l} ɛ_{a H_{2} O} ρ_{i l} + 3 [{ɛ'}_{s 1 i l} (ɛ_{a i l} - ɛ_{a H_{2} O}) + {ɛ'}_{s 2 i l} ɛ_{a H_{2} O}] ρ_{i l}^{2} + 3 {ɛ'}_{s 2 i l} (ɛ_{a i l} - ɛ_{a H_{2} O}) ρ_{i l}^{3} .

(17)

With the values obtained for ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail (see Sect. 3.2) it is possible to show that the contribution of the term proportional to

ρ_{i l}^{3}

remains smaller than about 1% for all values of ρ_il at which measurements with the method of absorption of water have been carried out. This term can be therefore neglected in Eq. (17) and a new linear relationship between

μ_{e f f}^{2} (ρ_{i l}) / ρ_{i l}

and ρ_il is obtained, but with different coefficients with respect to Eq. (8). Using this new relationship, from the slope and the intercept of the line that best fits

μ_{e f f}^{2} (ρ_{i l}) / ρ_{i l}

as a function of ρ_il we obtain

{ɛ'}_{s 1 i l} = \frac{I_{i l}}{3 ɛ_{a H_{2} O}}

(18)

ɛ_{a i l} = ɛ_{a H_{2} O} (1 + \frac{S_{i l}}{I_{i l}} - \frac{{ɛ'}_{s 2 i l}}{{ɛ'}_{s 1 i l}}) .

(19)

Comparison of Eqs. (18) and (19) with Eqs. (9) and (10) shows that the value of ɛ′_sil obtained using the method of water absorption with the assumption of independent scattering, actually represents the term ɛ′ _s1il of Eq. (18), and that the correct value of ɛ_ail can be obtained using for the ratio ɛ′ _s2il/ɛ′ _s1il in Eq. (19) the value obtained at step 4.

In conclusion, to include the effect of dependent scattering three parameters are necessary: ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail . ɛ′ _s1il is obtained from Eq. (18), ɛ_ail from Eq. (19) using for ɛ′ _s2il/ɛ′ _s1il the value obtained with the method of adding absorption, and ɛ′ _s2il from the values of ɛ′ _s1il and of ɛ′ _s2il/ɛ′ _s1il. Of the results obtained with the method of adding absorption we therefore use only the ratio ɛ′ _s2il/ɛ′ _s1il that, as previously pointed out, is independent of the error on the calibration of India ink.

Step 6: Finally, a more accurate calibration for the specific absorption coefficient of the India ink can be obtained using in Eq. (12) (Step 2) the value of μ′_sil (ρ_il ) obtained from Eq. (15) that includes the effect of dependent scattering.

2.3. Experimental Setup

Measurements have been carried out using two CW laser diodes at λ = 751 nm and 833 nm. The emitted power was 3 and 30 mW respectively. To measure the extinction coefficient we have carried out measurements of collimated transmittance at different concentrations. We have used a scattering cell 34.5 mm thick, and transmittance was measured with a photodiode and a lock-in amplifier. The experimental setup was similar to that of Ref. [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

] with an acceptance angle of the detection system of 7 mrad. With this small acceptance angle the error on the extinction coefficient due to the unavoidable fraction of scattered received power was negligible, and the specific extinction coefficient ɛ_eil has been obtained with an error smaller than 0.5%.

The setup for multidistance measurements of fluence in the infinite medium was similar to that of Ref. [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

]. Two thin fibers with a small diffusive tip (outer diameter 0.5 mm) having a substantially uniform radiation pattern (US Patent Number 6,071,302 ‘Phototherapeutic apparatus for wide-angle diffusion’) were used to illuminate the medium and to measure the fluence. The interfibre distance was varied between 10 and 35 mm (1 mm step) with a computer controlled translation stage and received photons were measured with a photomultiplier and a lock-in amplifier. The small error on the measured fluence (random error of about 0.1%) enabled us to obtain the effective attenuation coefficient with high accuracy (the standard deviation due to random errors was of less than 0.05%). Measurements have been carried out inside a cylindrical vessel with volume of about 3L. For the interfibre distances and for the values of μ′_s and μ_a at which measurements have been carried out this volume was sufficiently large to act as an infinite medium. The volume concentration ρ_il has been obtained from the weight of Intralipid and water using the value 0.988 for the relative density of Intralipid 20% with respect to water [21

R. Michels, F. Foschum, and A. Kienle, “Optical properties of fat emulsions,” Opt. Express 16, 5907–5925 (2008). [CrossRef]

]. The absorber we have used was a sample of prediluted Higgins India ink. Predilution was necessary since, due to its high value of absorption coefficient (hundreds of mm⁻¹), it would be difficult to weigh with good accuracy the very small quantities of ink necessary to obtain the concentrations of interest for our measurements. For more detailed information on the use of India ink we refer to Ref. [23

P. Di Ninni, F. Martelli, and G. Zaccanti, “The use of India ink in tissue-simulating phantoms,” Opt. Express 18, 26854–26865 (2010). [CrossRef]

]. The concentrations both for Intralipid and for ink have been determined with error smaller than 0.1%.

3. Experimental Results

All the measurements have been carried out using different bags of Intralipid 20% from the same batch. As also shown in Ref. [17

P. Di Ninni, F. Martelli, and G. Zaccanti, “Intralipid: towards a diffusive reference standard for optical tissue phantoms,” Phys. Med. Biol. 56, N21–N28 (2011). [CrossRef]

] we did not observe any significant difference in the optical properties of different bags.

3.1. Monitoring the Microphysical Properties of Diluted Intralipid

We have measured the specific extinction coefficient at λ = 751 nm and 833 nm carrying out measurements on dilutions with similar concentrations but prepared in different ways. Examples of results are reported in Fig. 1. Panels a)–c) pertain to λ = 751 nm, panels d)–f) to 833 nm. The error bars are not shown since smaller than the marks. Panels a) and d) show the results for a dilution prepared taking a sample directly from the bag of Intralipid 20%. To obtain accurate values for the concentrations the sample has been prediluted (1.185 g of Intralipid in 98.348 g of water).

Fig. 1 Examples of measurements of collimated transmittance on dilutions of Intralipid 20% prepared in different ways. Panels a) and d): predilution with 1.185 g of Intralipid in 98.348 g of water; panels b) and e): predilution with 267.0 g of Intralipid in 2478 g of water; panels c) and f): predilution with 80.11 g of Intralipid and 2.49 g of prediluted ink in 2458 g of water. Panels a)–c) pertain to λ = 751 nm, panels d)–f) to 833 nm.

Panels b) and e) show the results of extinction measurements on a sample of the dilution used for measuring the optical properties of Intralipid 20% with the method of water absorption (step 1). We have used a sample of the dilution in the cylindrical vessel at the end of measurements at λ = 751 nm (267.0 g of Intralipid 20% in 2478 g of water).

Panels c) and f) show the results of measurements on a sample of the dilution of Intralipid and India ink taken from the cylindrical vessel at the end of measurements at λ = 751 nm for step 2 (80.11 g of Intralipid 20% and 2.49 g of prediluted ink in 2458 g of water). The contribution to the extinction coefficient due to the small quantity of prediluted ink was negligible (smaller than 0.1%). The value expected for ɛ_eil from these measurements is therefore the same obtained from measurements on other samples provided the added ink does not change the microphysical properties of Intralipid. The results we have obtained for ɛ_eil from measurements at 751 nm (panels a)–c)) were 66.3±0.3, 66.3±0.3, and 66.2±0.3 mm⁻¹ respectively, and at 833 nm (panels d)–f)) 50.9±0.3, 51.1±0.3, and 50.8±0.3 mm⁻¹ respectively. Within the standard deviation the results are indistinguishable both at 751 nm and at 833 nm, showing that the microphysical properties of diluted Intralipid do not depend on the way the dilution has been prepared. In particular, they are not affected by the addition of the India ink used for measurements at steps 2 and 3.

Measurements of collimated transmittance have been also carried out on dilutions of Intralipid 20% prepared using buffered water. To adjust the pH of water to the value (pH ≅ 8) measured for Intralipid [21

R. Michels, F. Foschum, and A. Kienle, “Optical properties of fat emulsions,” Opt. Express 16, 5907–5925 (2008). [CrossRef]

] a small quantity of NaOH has been added. The results we have obtained were indistinguishable with respect to those obtained using purified water.

3.2. Dependent Scattering

The experimental results at λ = 751 nm and 833 nm are summarized in Figs. 2 and 3 respectively. Different panels show the results pertaining to the different steps of the measuring procedure together with the straight lines that best fit the results.

Fig. 2 Dependent scattering: Experimental results for λ = 751 nm. Each panel reports the experimental results (marks) together with the straight line that best fits the results. The error bars are shown only when larger than the marks. Panel a): results for

μ_{e f f}^{2}

as a function of ρ_il used to obtain ɛ′_sil and ɛ_ail with the method of absorption of water (step 1) and ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail (step 5). Panel b) results for

μ_{e f f}^{2}

as a function of ρ_ink used to obtain ɛ_aink (steps 2 and 6). Panels c) and d): examples of results for

μ_{e f f}^{2}

as a function of ρ_ink used to obtain μ′_s (ρ_il ) and μ_a (ρ_il ) with the method of adding absorption (step 3) for ρ_il = 0.0204 and 0.1408 respectively. Panels e) and f): the results for μ′_s (ρ_il ) and μ_a (ρ_il ) of step 3 are plotted as a function of ρ_il . These results are used to evaluate the effect of dependent scattering (step 4).

Fig. 3 Same as Fig. 2, but for λ = 833 nm. The results displayed in panels c) and d) refer to ρ_il = 0.0214 and 0.1082 respectively.

The results of panel a) are used to obtain ɛ′_sil and ɛ_ail with the method of water absorption (Eqs. (9) and (10), step 1). For the absorption coefficient of water, since published data show an appreciable spread of values and often the information on the standard error is not provided (see Ref. [12

F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]

] for a deeper discussion), we have used the values obtained from measurements of attenuation in pure water. Using the method described in [12

F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]

] we obtained ɛ _{aH
₂
O} = (2.77±0.02)×10⁻³ mm⁻¹ and (3.55±0.02)×10⁻³ mm⁻¹ at 751 and 833 nm respectively.

The results used to measure ɛ_aink (step 2) are plotted in panel b). The concentration of Intralipid 20% was ρ_il = 0.0289 at 751 nm and 0.0306 at 833 nm, and was not appreciably changed by the addition of the prediluted ink (ρ_il changed less than 0.1%).

Panels c) and d) show two examples of using the method of adding absorption to measure μ′_s (ρ_il ) and μ_a (ρ_il ) (step 3). Panel c) reports measurements for the smallest value of ρ_il we considered (ρ_il = 0.0204 at 751 nm and 0.0214 at 833 nm) and panel d) for the largest value (ρ_il = 0.1408 at 751 nm and 0.1082 at 833 nm). The values of μ′_s (ρ_il ) and μ_a (ρ_il ) measured at step 3 are displayed in panels e) and f) respectively. These data are used to evaluate the effect of dependent scattering (step 4). To highlight the small deviations from the linearity observed for μ′_s in panel e) we displayed the ratio μ′_s (ρ_il )/ρ_il . The ratio would be independent of ρ_il if the microphysical properties of Intralipid are not affected by the dilution or by the addition of India ink and the assumption of independent scattering (Eq. (6)) is fulfilled. Since it has been shown in Sect. 3.1 that dilution does not affect the microphysical properties, variations of the ratio μ′_s (ρ_il )/ρ_il are therefore ascribed to the effect of dependent scattering. The experimental results show that μ′_s (ρ_il )/ρ_il slightly decreases as ρ_il increases, i.e., the dependent scattering slightly decreases the reduced scattering efficiency. As anticipated in Sect. 2.2 the results are fitted reasonably well by a straight line, and the ratio S _1il/I _1il is used to obtain ɛ′ _s2il/ɛ′ _s1il (Eq. (15)). Panel f) shows μ_a (ρ_il ) as a function of ρ_il . Also these results are fitted reasonably well by a straight line. This indicates that the dependent scattering does not appreciably affect absorption and also that absorption of Intralipid is smaller than absorption of water. From the slope and the intercept of the straight line the values of ɛ_ail and ɛ _{aH
₂
O} are obtained using Eq. (7).

The results on the effect of dependent scattering are then used to reanalyze the data of panel a) and the final results for ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail are obtained (step 5). Finally, using for μ′_s (ρ_il ) the value evaluated with the results obtained at step 5), a more accurate value for ɛ_aink is obtained from the data of panel b) (step 6).

The results we have obtained are summarized in Table 1. The table shows the values obtained at different steps together with the corresponding standard deviation. The standard deviation has been evaluated considering only the contribution of random errors. We point out that Table 1 shows for ɛ_ail both the results obtained at steps 4 and 5. The two values are consistent, however since the value obtained at step 4 is affected by the error on the calibration of India ink, the value obtained at step 5 is more reliable. We also note that the value of ɛ _{aH
₂
O} obtained at step 4 is consistent with that obtained from direct attenuation measurements in pure water (section 3.2 and Ref. [12

F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]

]).

Table 1 Summary of the Experimental Results Obtained at Different Steps

	λ (nm)	751	833
STEP 1	ɛ′_sil (mm⁻¹)	21.16±0.04	18.57±0.03
STEP 1	ɛ_ail (mm⁻¹)	(1.45±0.09)×10⁻³	(1.67±0.09)×10⁻³
STEP 2	ɛ_aink (mm⁻¹)	315.0±0.4	286.1±0.8
STEP 4	ɛ′ _s2il/ɛ′ _s1il	−0.27±0.05	−0.29±0.05
	ɛ_ail (mm⁻¹)	(2.14±0.06)×10⁻³	(2.46±0.36)×10⁻³
	ɛ _{aH ₂ O} (mm⁻¹)	(2.77±0.01)×10⁻³	(3.56±0.02)×10⁻³
STEP 5	ɛ′ _s1il (mm⁻¹)	21.16±0.04	18.57±0.03
	ɛ′ _s2il (mm⁻¹)	−5.8±0.9	−5.3±0.9
	ɛ_ail (mm⁻¹)	(2.2±0.2)×10⁻³	(2.7±0.3)×10⁻³
STEP 6	ɛ_aink (mm⁻¹)	319±1	289±2

4. Discussion and Conclusions

Intralipid 20% has been recently proposed as a first step towards a reference diffusive medium for tissue-simulating phantoms [17

P. Di Ninni, F. Martelli, and G. Zaccanti, “Intralipid: towards a diffusive reference standard for optical tissue phantoms,” Phys. Med. Biol. 56, N21–N28 (2011). [CrossRef]

]. In this paper we have investigated the limits of applicability of the linear relationships between the optical properties of diluted Intralipid and the volume concentration. The validity of these linear relationships is commonly assumed when liquid tissue phantoms are prepared, but they are applicable only if the microphysical properties of scattering particles are not affected by dilution, and if the independent scattering approximation is fulfilled. The results presented in section 3 showed that the microphysical properties of suspended particles are not appreciably affected by the dilution in water or by the addition of India ink. As for the effect of dependent scattering experimental results showed that the absorption coefficient is not affected, while the reduced scattering efficiency slightly decreases when concentration increases. Deviations from linearity can be accounted for with a corrective term proportional to the square of the concentration. However, for concentrations of practical interest for tissue phantoms (ρ_il < 0.1 to obtain μ′_s < 2 mm⁻¹) deviations remain within about 2%.

The effect of dependent scattering on the extinction coefficient has been investigated in many papers [19

A. Ishimaru and Y. Kuga , “Attenuation constant of a coherent field in a dense distribution of particles,” J. Opt. Soc. Am. 72 1317–20 (1982). [CrossRef]

,20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

,24

B. L. Drolen and C. L. Tien, “Independent and dependent scattering in packed-sphere systems,” J. Thermophys. 1, 63–68 (1987). [CrossRef]

–26

J. Kalkman, A. Bykov, D. Faber, and T. van Leeuwen, “Multiple and dependent scattering effects in Doppler optical coherence tomography,” Opt. Express 18, 3883–3892 (2010). [CrossRef]

] but to our knowledge results on μ_a and μ′_s have been reported only in [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

]. Measurements reported in Ref. [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

] have been carried out at λ = 632.8 nm and were focussed on high concentrations of Intralipid 20% (up to ρ_il = 1). The results reported in this paper are in agreement with those of Ref. [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

] both for μ_a and μ′_s : absorption is not affected by dependent scattering while the efficiency of reduced scattering decreases as ρ_il increases. From measurements at 751 and 833 nm for moderate concentrations (ρ_il < 0.15) we obtained for the ratio ɛ′ _s2il/ɛ′ _s1il the values −0.27±0.04 and −0.29±0.05 respectively. The value obtained from Ref. [20

G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]

] for λ = 632.8 nm and high concentrations was −0.37.

The results we obtained for ɛ_eil , ɛ_ail , and ɛ′_sil are in excellent agreement with the results we reported in Ref. [23

P. Di Ninni, F. Martelli, and G. Zaccanti, “The use of India ink in tissue-simulating phantoms,” Opt. Express 18, 26854–26865 (2010). [CrossRef]

] for different batchs of Intralipid 20%. For comparisons with other published data we refer to the discussion in Ref. [17

P. Di Ninni, F. Martelli, and G. Zaccanti, “Intralipid: towards a diffusive reference standard for optical tissue phantoms,” Phys. Med. Biol. 56, N21–N28 (2011). [CrossRef]

The results presented in this paper show that the optical properties of dilutions of Intralipid can be predicted with high accuracy for all concentrations of practical interest for tissue phantoms once the parameters ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail have been determined. The accuracy is ultimately limited by random and systematic errors on the parameters ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail and on the concentration ρ_il . With our setup the parameters ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail have been determined with a random error of 0.2%, 15%, and 10%, respectively. As for systematic errors, from the discussion in Ref. [12

F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]

] comes that the main contribution is due to the uncertainty in positioning the fibers and to the error on the absorption coefficient of water. The absorption coefficient of water has been measured with error of 0.5%. The corresponding systematic errors on ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail are of about 0.5%. In positioning the source and the receiving fibers we estimated an uncertainty of 0.1 mm. The corresponding systematic errors are of 1%, 1%, and 5% for ɛ′ _s1il, ɛ′ _s2il, and ɛ_ail respectively, when measurements of fluence are carried at interfibre distances between 10 and 35 mm, as is the case for the results reported in section 3. However, this error can be easily reduced if measurements are carried out for larger interfibre distances. As an example, with an uncertainty of 0.1 mm the errors are reduced to 0.2%, 0.2%, and 1%, i.e., smaller than random errors, if measurements are carried out between 20 and 45 mm. The disadvantage in carrying out measurements at larger interfibre distances is the larger volume of diffusive medium necessary to mimic the infinite medium.

We point out that the intrinsic optical properties of Intralipid have been presented referring to the volume concentration ρ_il . Since the density of Intralipid 20% is close to the density of water, it is possible with negligible error to obtain the optical properties referred to the weight concentration ρ_wil taking into account that for concentrations of practical interest is ρ_il ≅ ρ_wil /0.988.

In conclusion we have shown that by using the method of water absorption, and including the corrective term to take into account the effect of dependent scattering on the scattering properties, it is possible to characterize the reduced scattering coefficient of diluted Intralipid with an overall error of less than 1% on μ′_sil , and that at the NIR wavelengths we investigated for dilutions of practical interest the absorption coefficient is practically equal to the absorption of pure water. Furthermore, using the calibrated Intralipid it is possible to measure the absorption coefficient of India ink with the same accuracy we have on the reduced scattering coefficient of Intralipid. Mixing calibrated Intralipid as a scattering medium and calibrated India ink as an absorbing medium it is therefore possible to obtain liquid phantoms with the desired scattering and absorption properties with errors smaller than 1%.

Also in the light of a recently published paper [11

], it seems difficult to measure the optical properties of a solid phantom with accuracy similar to that obtained for the Intralipid and India ink phantoms. Liquid phantoms with well known optical properties can be therefore useful to investigate and to understand the ultimate accuracy obtainable with different techniques proposed for measuring the optical properties of diffusive media, and more in general for studying photon migration. Using black containers with small transparent windows for input and output of emitted and of received light it is possible to obtain homogeneous diffusive media with well defined and known boundary conditions whose optical properties can be easily adjusted changing the concentration of scatterers and of absorbers. Layered media can be also made by using thin foils of slightly scattering material to separate different layers [27

S. Del Bianco, F. Martelli, F. Cignini, G. Zaccanti, A. Pifferi, A. Torricelli, A. Bassi, P. Taroni, and R. Cubeddu, “Liquid phantom for investigating light propagation through layered diffusive media,” Opt. Express 12, 2102–2111 (2004). [CrossRef]

]. Furthermore, inhomogeneities can be easily put inside the liquid to mimic the effect of abnormalities on photon migration. With liquid phantoms, being possible to repeat measurements before and after putting the inhomogeneity inside without changing anything else, also small perturbations on light propagation can be measured with high accuracy [28

G. Zaccanti, L. Alianelli, C. Blumetti, and S. Carraresi, “Method for measuring the mean time of flight spent by photons inside a volume element of a highly diffusing medium,” Opt. Lett. 24, 1290–1292 (1999). [CrossRef]

Acknowledgments

The research leading to these results has received funding from the European Community’s Seventh Framework Programme [FP7/2007–2013] under grant agreement n FP7-HEALTH-F5-2008-201076. The authors thank Danilo Marcucci for his help during the experiments.

References and links

1.	S. Madsen, B. Wilson, M. Patterson, Y. Park, S. Jacques, and Y. Hefetz, “Experimental tests of a simple diffusion model for the estimation of scattering and absorption coefficients of turbid media from time-resolved diffuse reflectance measurements,” Appl. Opt. 31, 3509–3517 (1992). [CrossRef] [PubMed]
2.	R. Cubeddu, M. Musolino, A. Pifferi, P. Taroni, and G. Valentini, “Time-resolved reflectance: a systematic study for application to the optical characterization of tissues,” IEEE J. Quantum Electron. 30, 2421–2430 (1994). [CrossRef]
3.	E. Alerstam, S. Andersson-Engels, and T. Svensson, “Improved accuracy in time-resolved diffuse reflectance spectroscopy,” Opt. Express 16, 10440–10454 (2008). [CrossRef]
4.	S. Fantini, M. Franceschini, and E. Gratton, “Semi-infinite-geometry boundary problem for light migration in highly scattering media: a frequency-domain study in the diffusion approximation,” J. Opt. Soc. Am. B 11, 2128–2138 (1994). [CrossRef]
5.	H. Xu and M. Patterson, “Determination of the optical properties of tissue-simulating phantoms from interstitial frequency domain measurements of relative fluence and phase difference,” Opt. Express 14, 6485–6501 (2006). [CrossRef]
6.	P. Marquet, F. P. Bevilacqua, C. D. Depeursinge, and E. B. de Haller, “Determination of reduced scattering and absorption coefficients by a single charge-coupled-device array measurement, part I: comparison between experiments and simulations,” Opt. Eng. 34, 2055–2063 (1995). [CrossRef]
7.	R. L. P. van Veen, W. Verkruysse, and H. J. C. M. Sterenborg, “Diffuse-reflectance spectroscopy from 500 to 1060 nm by correction for inhomogeneously distributed absorbers,” Opt. Lett. 27, 246–248 (2002). [CrossRef]
8.	C. Chen, J. Q. Lu, H. Ding, K. Jacobs, Y. Du, and X. H. Hu, “A primary method for determination of optical parameters of turbid samples and application to intralipid between 550 and 1630nm,” Opt. Express 14, 7420–7435 (2006). [CrossRef]
9.	N. Rajaram, T. H. Nguyen, and J. W. Tunnell, “Lookup table-based inverse model for determining optical properties of turbid media,” J. Biomed. Opt. 13, 050501 (2008). [CrossRef]
10.	A. Pifferi, A. Torricelli, A. Bassi, P. Taroni, R. Cubeddu, H. Wabnitz, D. Grosenick, M. Mller, R. Macdonald, J. Swartling, T. Svensson, S. Andersson-Engels, R. L. P. van Veen, H. J. C. M. Sterenborg, J.-M. Tualle, H. L. Nghiem, S. Avrillier, M. Whelan, and H. Stamm, “Performance assessment of photon migration instruments: the MEDPHOT protocol,” Appl. Opt. 44, 2104–2114 (2005). [CrossRef]
11.	J. Bouchard, I. Veilleux, R. Jedidi, I. Noiseux, M. Fortin, and O. Mermut, “Reference optical phantoms for diffuse optical spectroscopy. Part 1 – Error analysis of a time resolved transmittance characterization method,” Opt. Express 18, 11495–11507 (2010). [CrossRef]
12.	F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]
13.	L. Spinelli, F. Martelli, A. Farina, A. Pifferi, A. Torricelli, R. Cubeddu, and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. Time-resolved method,” Opt. Express 15, 6589–6604 (2007). [CrossRef]
14.	B. W. Pogue and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt. 11, 041102 (2006). [CrossRef]
15.	M. Zude (Ed.), Optical Monitoring of Fresh and Processed Agricultural Crops , (Contemporary Food Engineering Series), (CRC Press, Boca Raton, Florida, 2009).
16.	J. Johansson, S. Folestad, M. Josefson, A. Sparn, C. Abrahamsson, S. Anderson-Engels, and S. Svanberg, “Time-resolved NIR/Vis spectroscopy for analysis of solids: Pharmaceutical tablets,” Appl. Spectrosc. 56, 725–731 (2002). [CrossRef]
17.	P. Di Ninni, F. Martelli, and G. Zaccanti, “Intralipid: towards a diffusive reference standard for optical tissue phantoms,” Phys. Med. Biol. 56, N21–N28 (2011). [CrossRef]
18.	M. I. Mishchenko, L. D. Travis, and A. A. Lacis, Scattering, Absorption, and Emission of Light by Small Particles , (Cambridge University Press, Cambridge, UK, 2002).
19.	A. Ishimaru and Y. Kuga , “Attenuation constant of a coherent field in a dense distribution of particles,” J. Opt. Soc. Am. 72 1317–20 (1982). [CrossRef]
20.	G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]
21.	R. Michels, F. Foschum, and A. Kienle, “Optical properties of fat emulsions,” Opt. Express 16, 5907–5925 (2008). [CrossRef]
22.	C. F. Bohren and D. R. Huffman, Absorption and scattering of light by small particles , (John Wiley and Sons, New York, 1983).
23.	P. Di Ninni, F. Martelli, and G. Zaccanti, “The use of India ink in tissue-simulating phantoms,” Opt. Express 18, 26854–26865 (2010). [CrossRef]
24.	B. L. Drolen and C. L. Tien, “Independent and dependent scattering in packed-sphere systems,” J. Thermophys. 1, 63–68 (1987). [CrossRef]
25.	G. Göbel, J. Kuhn, and J. Fricke, “Dependent scattering effects in latex-sphere suspensions and scattering powders,” Waves in Random Media 5, 413–426 (1995). [CrossRef]
26.	J. Kalkman, A. Bykov, D. Faber, and T. van Leeuwen, “Multiple and dependent scattering effects in Doppler optical coherence tomography,” Opt. Express 18, 3883–3892 (2010). [CrossRef]
27.	S. Del Bianco, F. Martelli, F. Cignini, G. Zaccanti, A. Pifferi, A. Torricelli, A. Bassi, P. Taroni, and R. Cubeddu, “Liquid phantom for investigating light propagation through layered diffusive media,” Opt. Express 12, 2102–2111 (2004). [CrossRef]
28.	G. Zaccanti, L. Alianelli, C. Blumetti, and S. Carraresi, “Method for measuring the mean time of flight spent by photons inside a volume element of a highly diffusing medium,” Opt. Lett. 24, 1290–1292 (1999). [CrossRef]

OCIS Codes
(160.4760) Materials : Optical properties
(170.6510) Medical optics and biotechnology : Spectroscopy, tissue diagnostics
(170.7050) Medical optics and biotechnology : Turbid media

ToC Category:
Calibration, Validation and Phantom Studies

History
Original Manuscript: June 8, 2011
Revised Manuscript: June 29, 2011
Manuscript Accepted: July 6, 2011
Published: July 14, 2011

Citation
Paola Di Ninni, Fabrizio Martelli, and Giovanni Zaccanti, "Effect of dependent scattering on the optical properties of Intralipid tissue phantoms," Biomed. Opt. Express 2, 2265-2278 (2011)
http://www.opticsinfobase.org/boe/abstract.cfm?URI=boe-2-8-2265

Sort: Author | Year | Journal | Reset

References

S. Madsen, B. Wilson, M. Patterson, Y. Park, S. Jacques, and Y. Hefetz, “Experimental tests of a simple diffusion model for the estimation of scattering and absorption coefficients of turbid media from time-resolved diffuse reflectance measurements,” Appl. Opt. 31, 3509–3517 (1992). [CrossRef] [PubMed]
R. Cubeddu, M. Musolino, A. Pifferi, P. Taroni, and G. Valentini, “Time-resolved reflectance: a systematic study for application to the optical characterization of tissues,” IEEE J. Quantum Electron. 30, 2421–2430 (1994). [CrossRef]
E. Alerstam, S. Andersson-Engels, and T. Svensson, “Improved accuracy in time-resolved diffuse reflectance spectroscopy,” Opt. Express 16, 10440–10454 (2008). [CrossRef]
S. Fantini, M. Franceschini, and E. Gratton, “Semi-infinite-geometry boundary problem for light migration in highly scattering media: a frequency-domain study in the diffusion approximation,” J. Opt. Soc. Am. B 11, 2128–2138 (1994). [CrossRef]
H. Xu and M. Patterson, “Determination of the optical properties of tissue-simulating phantoms from interstitial frequency domain measurements of relative fluence and phase difference,” Opt. Express 14, 6485–6501 (2006). [CrossRef]
P. Marquet, F. P. Bevilacqua, C. D. Depeursinge, and E. B. de Haller, “Determination of reduced scattering and absorption coefficients by a single charge-coupled-device array measurement, part I: comparison between experiments and simulations,” Opt. Eng. 34, 2055–2063 (1995). [CrossRef]
R. L. P. van Veen, W. Verkruysse, and H. J. C. M. Sterenborg, “Diffuse-reflectance spectroscopy from 500 to 1060 nm by correction for inhomogeneously distributed absorbers,” Opt. Lett. 27, 246–248 (2002). [CrossRef]
C. Chen, J. Q. Lu, H. Ding, K. Jacobs, Y. Du, and X. H. Hu, “A primary method for determination of optical parameters of turbid samples and application to intralipid between 550 and 1630nm,” Opt. Express 14, 7420–7435 (2006). [CrossRef]
N. Rajaram, T. H. Nguyen, and J. W. Tunnell, “Lookup table-based inverse model for determining optical properties of turbid media,” J. Biomed. Opt. 13, 050501 (2008). [CrossRef]
A. Pifferi, A. Torricelli, A. Bassi, P. Taroni, R. Cubeddu, H. Wabnitz, D. Grosenick, M. Mller, R. Macdonald, J. Swartling, T. Svensson, S. Andersson-Engels, R. L. P. van Veen, H. J. C. M. Sterenborg, J.-M. Tualle, H. L. Nghiem, S. Avrillier, M. Whelan, and H. Stamm, “Performance assessment of photon migration instruments: the MEDPHOT protocol,” Appl. Opt. 44, 2104–2114 (2005). [CrossRef]
J. Bouchard, I. Veilleux, R. Jedidi, I. Noiseux, M. Fortin, and O. Mermut, “Reference optical phantoms for diffuse optical spectroscopy. Part 1 – Error analysis of a time resolved transmittance characterization method,” Opt. Express 18, 11495–11507 (2010). [CrossRef]
F. Martelli and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method,” Opt. Express 15, 486–500 (2007). [CrossRef]
L. Spinelli, F. Martelli, A. Farina, A. Pifferi, A. Torricelli, R. Cubeddu, and G. Zaccanti, “Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. Time-resolved method,” Opt. Express 15, 6589–6604 (2007). [CrossRef]
B. W. Pogue and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt. 11, 041102 (2006). [CrossRef]
M. Zude (Ed.), Optical Monitoring of Fresh and Processed Agricultural Crops , (Contemporary Food Engineering Series), (CRC Press, Boca Raton, Florida, 2009).
J. Johansson, S. Folestad, M. Josefson, A. Sparn, C. Abrahamsson, S. Anderson-Engels, and S. Svanberg, “Time-resolved NIR/Vis spectroscopy for analysis of solids: Pharmaceutical tablets,” Appl. Spectrosc. 56, 725–731 (2002). [CrossRef]
P. Di Ninni, F. Martelli, and G. Zaccanti, “Intralipid: towards a diffusive reference standard for optical tissue phantoms,” Phys. Med. Biol. 56, N21–N28 (2011). [CrossRef]
M. I. Mishchenko, L. D. Travis, and A. A. Lacis, Scattering, Absorption, and Emission of Light by Small Particles , (Cambridge University Press, Cambridge, UK, 2002).
A. Ishimaru and Y. Kuga , “Attenuation constant of a coherent field in a dense distribution of particles,” J. Opt. Soc. Am. 721317–20 (1982). [CrossRef]
G. Zaccanti, S. Del Bianco, and F. Martelli, “Measurements of optical properties of high density media,” Appl. Opt. 42, 4023–4030 (2003). [CrossRef]
R. Michels, F. Foschum, and A. Kienle, “Optical properties of fat emulsions,” Opt. Express 16, 5907–5925 (2008). [CrossRef]
C. F. Bohren and D. R. Huffman, Absorption and scattering of light by small particles , (John Wiley and Sons, New York, 1983).
P. Di Ninni, F. Martelli, and G. Zaccanti, “The use of India ink in tissue-simulating phantoms,” Opt. Express 18, 26854–26865 (2010). [CrossRef]
B. L. Drolen and C. L. Tien, “Independent and dependent scattering in packed-sphere systems,” J. Thermophys. 1, 63–68 (1987). [CrossRef]
G. Göbel, J. Kuhn, and J. Fricke, “Dependent scattering effects in latex-sphere suspensions and scattering powders,” Waves in Random Media 5, 413–426 (1995). [CrossRef]
J. Kalkman, A. Bykov, D. Faber, and T. van Leeuwen, “Multiple and dependent scattering effects in Doppler optical coherence tomography,” Opt. Express 18, 3883–3892 (2010). [CrossRef]
S. Del Bianco, F. Martelli, F. Cignini, G. Zaccanti, A. Pifferi, A. Torricelli, A. Bassi, P. Taroni, and R. Cubeddu, “Liquid phantom for investigating light propagation through layered diffusive media,” Opt. Express 12, 2102–2111 (2004). [CrossRef]
G. Zaccanti, L. Alianelli, C. Blumetti, and S. Carraresi, “Method for measuring the mean time of flight spent by photons inside a volume element of a highly diffusing medium,” Opt. Lett. 24, 1290–1292 (1999). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper

Figures


Fig. 1	Fig. 2	Fig. 3

« Previous Article | Next Article »

OSA is a member of CrossRef.

OpticsInfoBase

Biomedical Optics Express

Effect of dependent scattering on the optical properties of Intralipid tissue phantoms

Article Tools

Article Navigation

Abstract

1. Introduction

2. Materials and Methods

2.1. Monitoring the Microphysical Properties of Diluted Intralipid

2.2. Dependent Scattering

2.3. Experimental Setup

3. Experimental Results

3.1. Monitoring the Microphysical Properties of Diluted Intralipid

3.2. Dependent Scattering

4. Discussion and Conclusions

Acknowledgments

References and links

References

Appl. Opt.

Appl. Spectrosc.

IEEE J. Quantum Electron.

J. Biomed. Opt.

J. Opt. Soc. Am.

J. Opt. Soc. Am. B

J. Thermophys.

Opt. Eng.

Opt. Express

Opt. Lett.

Phys. Med. Biol.

Waves in Random Media

Other

Cited By

Figures

Select a Conference
Session or Paper #
Year