OSA's Digital Library

Virtual Journal for Biomedical Optics

Virtual Journal for Biomedical Optics

| EXPLORING THE INTERFACE OF LIGHT AND BIOMEDICINE

  • Editor: Gregory W. Faris
  • Vol. 2, Iss. 11 — Nov. 26, 2007
« Show journal navigation

Discrimination between Doppler-shifted and non-shifted light in coherence domain path length resolved measurements of multiply scattered light

B. Varghese, V. Rajan, T.G. van Leeuwen, and W. Steenbergen  »View Author Affiliations


Optics Express, Vol. 15, Issue 20, pp. 13340-13350 (2007)
http://dx.doi.org/10.1364/OE.15.013340


View Full Text Article

Acrobat PDF (157 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We show a novel technique to distinguish between Doppler shifted and unshifted light in multiple scattering experiments on mixed static and dynamic media. With a phase modulated low coherence Mach-Zehnder interferometer, optical path lengths of shifted and unshifted light and path length dependent Doppler broadening are measured in a two-layer tissue phantom, with a superficial static layer of different thickness covering a semi-infinite dynamic medium having identical optical properties. No Doppler broadening is observed until a certain optical path length depending on the thickness of the superficial static layer. From the minimum optical path length corresponding to the Doppler-shifted light the thickness of the static layer that overlies the dynamic layer can be estimated. Validation of the experimentally determined thickness of the static layer is done with the Doppler Monte Carlo technique. This approach has potential applications in discriminating between statically and dynamically scattered light in the perfusion signal and in determining superficial burn depths.

© 2007 Optical Society of America

1. Introduction

Laser Doppler flowmetry (LDF), a non-invasive technique for monitoring blood microcirculation, characterizes the time-varying signal arising from the temporal variations in the speckle pattern to estimate the perfusion in biological tissues [1

1. A. P. Shepherd and P. Å. öberg, Laser-Doppler Blood Flowmetry (Kluwer Academic, Boston, 1990).

]. In this technique, the coherent light delivered into the tissue through an optical fiber follows different trajectories and reaches different depths and is guided to the detector through a spatially separated light-collecting fiber. Light which is Doppler shifted by moving red blood cells will be mixed with the non-frequency shifted light scattered from surrounding static tissue matrices, resulting in photodetector current fluctuations. Perfusion is measured based on the spectral analysis of the frequency components of this fluctuating photocurrent. However, skin perfusion measurements are complex not only due to the heterogeneous vascular structure and but also due to the variance in the penetration depth of the detected photons. The depth of light penetration in the tissue depends on the wavelength of the light source, the distance between the transmitting and receiving fibers and the optical properties of the tissue. Thus for a constant perfusion, the LDF output signal is affected by probe induced variations and by the changes in the tissue optical properties in terms of absorption and scattering [2

2. M. Larsson, W. Steenbergen, and T. Strömberg, “Influence of optical properties and fibre separation on laser Doppler flowmetry,” J. Biomed. Opt. 7, 236–243 (2001). [CrossRef]

]. This dependence leads to uncertainties in interpreting the Doppler shifted and non-shifted fraction of photons and also in discriminating the fraction of light scattered from superficial (nutrient) and deeper (thermoregulatory) layers of skin. As such, depth resolved perfusion information could be obtained by varying the distance between transmitting and receiving fibers [2

2. M. Larsson, W. Steenbergen, and T. Strömberg, “Influence of optical properties and fibre separation on laser Doppler flowmetry,” J. Biomed. Opt. 7, 236–243 (2001). [CrossRef]

] or by using a multichannel laser Doppler probe [3

3. A. Liebert, M. Leahy, and R. Maniewski, “Multichannel laser-Doppler probe for blood perfusion measurements with depth discrimination,” Med. Bio. Eng. Comp. 36, 740–747 (1998). [CrossRef]

]. However, these methods still give no control over the optical path length traveled by the detected light. Coherence domain path length resolved optical Doppler perfusion monitoring, of which the basic technique is presented in this work, may overcome this limitation of conventional LDF techniques and will enable to correctly interpret the inter-and intra-individual variations in the LDF readings introduced by the variance in individual photon path lengths (e.g., length and depth) due to changes in tissue optical properties.

The technique involves a phase modulated low coherence Mach-Zehnder interferometer with two spatially separated light-delivering and light-collecting fibers as used in conventional laser Doppler perfusion monitors [4

4. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Path length resolved measurements of multiple scattered photons in static and dynamic turbid media using phase modulated low coherence interferometry,” J. Biomed. Opt. 12, 024020 (2007). [CrossRef] [PubMed]

]. Compared to both the Michelson [6

6. A. Wax, C. Yang, R.R. Dasari, and M.S. Feld, “Path-length-resolved dynamic light scattering: modeling the transition from single to diffusive scattering,” Appl. Opt. 40, 4222–4227 (2001). [CrossRef]

,7

7. K.K. Bizheva, A.M. Siegel, and D.A. Boas, “Path-length-resolved dynamic light scattering in highly scattering random media: The transition to diffusing wave spectroscopy,” Phys. Rev. E , 58 (1998). [CrossRef]

] and the Mach-Zehnder [8

8. A.L. Petoukhova, W. Steenbergen, and F.F.M. de Mul, “Path-length distribution and path-length resolved Doppler measurements of multiply scattered photons by use of low-coherence interferometer,” Opt. Lett. 26, 1492–1494 (2001). [CrossRef]

,9

9. A.L. Petoukhova, W. Steenbergen, T.G. van Leeuwen, and F. F. M. de Mul, “Effects of absorption on coherence domain path length resolved dynamic light scattering in the diffuse regime,” Appl. Phys. Lett. 81, 595–597 (2002). [CrossRef]

] based interferometric measurements that depended only on the Doppler shifted fraction of photons, our method is able to measure path length resolved information of non-shifted and Doppler shifted fractions of photons. In a preliminary study, we showed that path length resolved Doppler information from static and dynamic turbid media can be obtained with increased signal-to-noise ratio by sinusoidally modulating the phase of the light in the reference arm [4

4. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Path length resolved measurements of multiple scattered photons in static and dynamic turbid media using phase modulated low coherence interferometry,” J. Biomed. Opt. 12, 024020 (2007). [CrossRef] [PubMed]

]. An example of the practical relevance of this approach is in measuring path length resolved information from mixed media, such as flow of moving red blood cells within static tissue matrices. Furthermore, we validated the optical path length distributions and path length dependent diffusion broadening of multiple scattered light with Monte Carlo simulations and Diffusive wave spectroscopy respectively [5

5. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Quantification of optical Doppler broadening and optical path lengths of multiply scattered light by phase modulated low coherence interferometry,” Opt. Express. 15(15), 9157–9165 (2007). [CrossRef] [PubMed]

].

In this study, we will demonstrate that using phase modulated low coherence interferometry, Doppler-shifted and non-shifted multiply scattered light from a mixed static and dynamic scattering medium can be distinguished and the optical path length of each can be measured. The technique was used on a two-layer tissue phantom, with a superficial static layer of different thickness covering a semi-infinite dynamic medium with identical optical properties. From the shape of the Doppler broadened phase modulation peak in the power spectrum, the Doppler shifted fraction of photons and their Doppler distribution can be estimated. Since the photons with large optical path lengths have greater probability to penetrate deep into the dynamic layer and thus to get Doppler shifted, we expect that the shortest optical path length at which the Doppler shifted light can be detected should correlate with the thickness of the static layer. We have performed Monte Carlo simulations to predict this transition optical path length, which is dependent on the thickness of the superficial static layer.

2. Materials and methods

The scattering medium (Fig. 1) that is used to study the depth sensitivity of our method consists of a semi-infinite dynamic layer of thickness 20 mm under a static layer of varying thickness (0.1–0.9 mm). The static and dynamic scattering phantoms are prepared with polystyrene microspheres of ∅ 0.77 μm to get the same optical properties (μ’s=2.0 mm-1, g=0.85), based on scattering cross sections following from Mie theory calculations. The static layer is prepared by mixing 4% (by weight) aqueous solution of poly (vinyl alcohol) (PVA) with the water suspension of polystyrene microspheres and 4% (by weight) aqueous sodium borate (borax) solution (cross-linking agent). PVA with a degree of hydrolysis greater than 99%, and an average molecular weight (MW) of 85000–140000 (Sigma-Aldrich, catalog nr 36 314-6), is used to prepare the aqueous solutions. The amount of borax solution added to the whole volume of the PVA-water-PS suspension (12 ml) is about 1.0 ml to produce a static gel medium. The static layer is placed in between two thin glass slides of thickness 150 μm and two spacers are placed to get the desired thickness for the static layer. Optical phantoms based on PVA gel are used for tissue mimicking phantoms in diffuse optical tomography [10

10. J. C. Hebden, B.D. Price, A. P. Gibson, and G. Royle, “A soft deformable tissue-equivalent phantom for diffuse optical tomography,” Phys. Med. Bio. 51, 5581–5590 (2006). [CrossRef]

]. Tissue perfusion within deeper layers was represented by particles in Brownian motion.

The path length distribution and average Doppler shift is measured from the noise corrected zero order moment M0 (area under the peak) and FWHM of the Doppler broadened interference peak appearing at the modulation frequency in the photodetector signal power spectrum [4

4. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Path length resolved measurements of multiple scattered photons in static and dynamic turbid media using phase modulated low coherence interferometry,” J. Biomed. Opt. 12, 024020 (2007). [CrossRef] [PubMed]

,5

5. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Quantification of optical Doppler broadening and optical path lengths of multiply scattered light by phase modulated low coherence interferometry,” Opt. Express. 15(15), 9157–9165 (2007). [CrossRef] [PubMed]

]. The power spectrum of photocurrent fluctuations of light that is scattered by a monodisperse suspension of particles undergoing Brownian motion and that is heterodyne detected is a Lorentzian distribution. The average Doppler shift, which depends on the scattering properties and photon path length within the strongly multiply scattering medium can be obtained from the FWHM [4

4. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Path length resolved measurements of multiple scattered photons in static and dynamic turbid media using phase modulated low coherence interferometry,” J. Biomed. Opt. 12, 024020 (2007). [CrossRef] [PubMed]

,5

5. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Quantification of optical Doppler broadening and optical path lengths of multiply scattered light by phase modulated low coherence interferometry,” Opt. Express. 15(15), 9157–9165 (2007). [CrossRef] [PubMed]

,11

11. A. G. Yodh, P. D. Kaplan, and D. J. Pine, “Pulsed diffusing-wave spectroscopy - high resolution through nonlinear optical gating,” Phys. Rev. B. 42, 4744–4747 (1990). [CrossRef]

]. The zeroth moment M0 of the heterodyne spectrum calculated around the modulation frequency of 22 kHz within a bandwidth of 2 kHz is the total power of the photocurrent fluctuations caused by the interference of reference light with the scattered light from the sample and it corresponds to the total intensity of photons with a certain path length [4

4. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Path length resolved measurements of multiple scattered photons in static and dynamic turbid media using phase modulated low coherence interferometry,” J. Biomed. Opt. 12, 024020 (2007). [CrossRef] [PubMed]

]. The FWHM of the interference signal is about 50 Hz when a mirror or a statically scattering medium replaces the sample. Thus, the noise corrected M0 in a bandwidth of 50 Hz at the modulation frequency corresponds to the intensity of non-shifted photons. The fraction of Doppler shifted photons (f) can be obtained by filtering out this contribution of non-shifted photons at the modulation frequency from the total intensity of photons. For a given thickness of the superficial static layer, the fraction of Doppler shifted and non-shifted photons for all optical path lengths can be estimated from the area of the corresponding optical path length distributions.

Fig. 1. Typical modulation peak appearing in the power spectra measured for short and large optical path lengths (Top). Schematic diagram of the light propagation through a two-layer scattering phantom, consisting of a top static layer and a bottom dynamic layer with same optical properties (Bottom). The model used for the estimation of the thickness (Inset).

In the Monte Carlo simulation [12

12. F. F. M. de Mul, M. H. Koelink, M. L. Kok, P. J. Harmsma, J. Greve, R. Graaff, and J. G. Aarnoudse, “Laser Doppler velocimetry and Monte Carlo simulations on models for blood perfusion in tissue,” Appl. Opt. 34, 6595–6611 (1995). [CrossRef] [PubMed]

], a five-layer system is defined in which light (λ =832nm) from a fiber with a diameter of 100 μm is delivered into a scattering medium consisting of a glass layer, a static scattering layer, a glass layer, and a dynamic layer with the same optical properties as the static scattering layer. Since the scattering properties and the particle number density of calibrated polystyrene sphere suspensions as used in this study are well known, it is possible to exactly mimic these properties in simulations. The parameters used for the static and dynamic scattering samples are refractive index, 1.33; absorption coefficient, 0.001 mm-1; Henyey-Greenstein scattering function, g=0.85 (∅0.77), reduced scattering coefficients (μs') of 2.0 mm-1 (g=0.85). The thickness of the two thin glass slides and the semi-infinite dynamic media were 150 μm and 20 mm respectively. The thickness of the superficial static layer was varied from 0.1 to 0.9 mm. Brownian motion of the particles within the dynamic layer was modeled by assigning a Gaussian velocity distribution with an average velocity of 0.1mm/s with random direction. The fifth layer is defined with a high absorption (μa =10 mm-1) to increase the simulation speed by removing from the simulation a minority of long path length photons. But we estimated from the statistics of the maximum scattering depth of Monte Carlo simulated detected photons that there are no photons reaching a depth beyond 10 mm in a scattering sample with lowest scattering level. So in this case there is no influence of highly absorbing second layer on the detected photons. Both the fibers are defined inside the scattering sample. In all simulations, 1*107 photons are injected into the sample, and each photon returning to the detection fiber (fiber diameter= 100 μm, fiber separation=300 μm, NA=0.29) is assumed to be detected, and its optical path length is recorded.

3. Experiments and results

Figure 1 shows the typical modulation peak appearing in the power spectra measured with a fiber distance of 300 μm for optical path lengths of 0.8 and 2.4 mm in the sample. The scattering medium consists of a static layer of thickness 0.1 mm overlying a dynamic layer of thickness 20 mm. For short optical path lengths, the reference light interferes with the photons scattered from the superficial static layer and are not Doppler shifted (circles and solid line). For large optical path length, the width of the peak is broadened due to the Doppler shifts imparted to the deeply penetrated multiply scattered photons by the Brownian motion of particles (squares and dashed line).

The optical path length distributions measured for the scattering media consisting of a dynamic layer (20 mm) under a static layer of thickness 0.1 and 0.7 mm with identical optical properties are shown in fig. 2 and fig. 3, respectively. By separating the 50 Hz bandwidth signal from the 2 kHz bandwidth signal, the optical path length distributions of the Doppler and non-Doppler shifted photons could be determined. For a static layer thickness of 0.1 mm, for optical path lengths up to 0.7 mm, the contributions of Doppler-shifted and non-shifted fractions of photons are equally strong, whereas the majority of the photons will be Doppler shifted for larger optical path lengths (Fig. 2). This is expected, since the photon mean free scattering path of 78 μm indicates that 28% of the photons will be scattered for the first time in the dynamic layer and thus will get Doppler-shifted. Many photons scattered for the first time in the static layer will be forwardly scattered into the dynamic layer. For a sample with a static layer thickness of 0.7 mm, the optical path length distributions of the Doppler and non-Doppler shifted photons are changed (Fig. 3). For this sample, the non-shifted fraction of photons exceeds the Doppler-shifted fraction up to 2.1 mm. For larger optical path lengths the majority of the light is Doppler shifted. When both fractions of photons averaged over all optical path lengths are estimated, the Doppler-shifted photons is decreased from 0.78 to 0.40 when the thickness of the static layer is increased from 0.1 to 0.7 mm. The fractions of Doppler shifted and non-shifted photons estimated for different thicknesses of the static layer are shown in fig. 4. As expected, the non-shifted fraction of photons exceeds the Doppler shifted fraction as the thickness of the superficial static layer increases.

Fig. 2. Optical path length distributions of non-shifted and Doppler shifted photons, measured for a two-layer tissue phantom consisting of a top static layer of thickness 0.1 mm and a bottom dynamic layer of thickness 20mm.

The average Doppler shift, measured from the FWHM of Doppler broadened phase modulation interference peak is represented in fig. 5 as a function of the optical path length. As depicted in fig. 5, in dynamic media, the average Doppler shift increases with the optical path length. But, when a thin layer of statically scattering medium is introduced on top of the dynamic media, the width measured for shorter optical path lengths corresponds to statically scattered light. For longer optical path lengths, the width of the peak increases with path length as the photons scattered from the dynamically scattering layer contribute to the interference signal. The optical path length at which the transition from the statically scattered to dynamically scattered light is detected depends on the thickness of the static layer.

Fig. 3. Optical path length distributions of non-shifted and Doppler shifted photons, measured for a two-layer tissue phantom consisting of a top static layer of thickness 0.7 mm and a bottom dynamic layer of thickness 20mm.
Fig. 4. The fraction of Doppler shifted and non-shifted photons (summed over all optical path lengths) estimated for different thickness of the superficial static layer.
Fig. 5. The average Doppler shift extracted from the FWHM of the phase modulation peak, as a function of the optical path length through the scattering medium, for various thickness of the static layer (0–.0.9 mm).
Fig. 6. The path length dependent Doppler broadening of multiply scattered light obtained from Monte Carlo simulations for various static layer thickness (0–0.9 mm).
Fig. 7. Thickness of the superficial static layer (including two glass layers) estimated from the path length dependent Doppler broadening. Experimental results (squares) and Monte Carlo simulations (circles).

To verify the dependence of minimum optical path length of Doppler-shifted photons on the thickness of the superficial static layer, we have performed Doppler Monte Carlo simulations. The results of simulations are shown in fig. 6. For short optical path lengths, no Doppler shifted light is measured and the minimum optical path lengths corresponding to the dynamically scattered light increases with increasing thickness of the static layer. The total thickness of the static layer and two glass slides estimated from the minimum optical path length corresponding to the Doppler shifted light is in good agreement with the results of Monte Carlo simulations (fig. 7). The error bars in fig. 7 indicate the variations in the shortest optical path length that may arise due to the path length resolution of ± 100 μm used in the experiment and Monte Carlo simulations. The straight line represents the real total layer thickness of the static layer and two glass slides.

4. Discussion

In this paper we have presented path length resolved measurements of multiply scattered light from mixed static and dynamic turbid media. Optical path length distributions and path length resolved Doppler shifts of multiply scattered light are measured from a two-layer scattering phantom, with a superficial static layer of various thicknesses covering a dynamic medium with identical optical properties. In the case of a static medium like a visco-elastic gel the width of the interference peak does not broaden with optical path length, whereas in the case of dynamic media more and more power is set to frequencies around the phase modulation frequency, resulting in Doppler broadening [4

4. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Path length resolved measurements of multiple scattered photons in static and dynamic turbid media using phase modulated low coherence interferometry,” J. Biomed. Opt. 12, 024020 (2007). [CrossRef] [PubMed]

]. Thus, the area of the Doppler broadened peak, excluding the statically scattered light contribution at the modulation frequency, forms an estimation of the fraction of Doppler shifted light. As expected, the Doppler fraction exceeds the complimentary non-shifted part for large optical path lengths. For increasing thickness of the superficial static layer, the estimated non-shifted fraction of photons increases.

The observed path length dependent line width broadening shown in fig. 5 results from the detection of multiply scattered photons and as the number of scattering events increases, the width of the peak increases with optical path length. However, the average Doppler shift per unit optical path length decreases with increase in the thickness of the static layer. For optical path lengths greater than 2.5 mm, the average Doppler shift measured with a static layer of thickness 0.1 mm is nearly equal to that measured in a completely dynamic medium, indicating that the contribution from the statically scattered light to the interference signal is significantly lower for large optical path lengths. A remarkable observation is that an average Doppler shift equivalent to the static medium is measured until a certain optical path length, which in turn depends on the thickness of the static layer. This corresponds to the photons that are scattered from the superficial static layer. After this transition optical path length, Doppler broadening of the peak is observed, since the photons with large optical path lengths penetrate into the dynamic layer. From the minimum optical path length (Lmin) of the Doppler-shifted light, we have estimated the thickness of the superficial static layer.

We validated the estimated thickness of the superficial static layer with Monte Carlo simulations. The experimentally determined transition optical path length and the thickness of the static layer estimated from this optical path length are in good agreement with the simulation results. The estimated thickness is marginally lower than the real one, except for the smallest and largest thickness where a larger difference is found. The marginal difference can be explained from the methodological accuracy. The assumed relation between the position of the retroreflector and the actual optical path length in the sample arm is based on the zero optical length determined with a reflecting mirror with accuracy ±50 μm. Another source of error might be in the simplification of assuming linear photon paths without taking refraction into account. However, when the refraction effects though glass layers are taken into account, the estimated shift is about 2 μm in a glass layer of 150 μm, which is not significant. The larger differences for the smallest and largest layer thickness cannot be explained from the applied method.

In spite of its simplicity, the accuracy of our static layer estimation is good. When this approach is extended to in vivo measurements, it seems to be justified to regard the tissue as a single static layer of thickness d with homogeneous average refractive index (ntissue) and from geometrical considerations, the thickness of the superficial static layer can be calculated. Sadhwani et al. reported an estimation of superficial static layer thickness that overlies a dynamic layer, by monitoring the speckle pattern at the surface as light propagates radially away from the point where the tissue was illuminated with a focused beam [13

13. A. Sadhwani, K. T. Schomacker, G. J. Tearney, and N. S. Nishioka, “Determination of Teflon thickness with laser speckle. I. Potential for burn depth diagnosis,” Appl. Opt. 35, 5727–5735 (1996). [CrossRef] [PubMed]

]. They showed that the diameter of the static speckle pattern increases linearly with the static layer thickness and the thickness of the static layer was empirically calculated from the speckle lobe size obtained by the convolution of a step function of the detector and a Gaussian function of the speckle lobe. As an application Sadhwani et al. mentioned the assessment of the depth of burn wounds. They modeled a burn wound as a superficial static layer with little or no blood flow overlying a perfused tissue and the quantitative estimation of the depth of the burn can be made from the radial extent of a speckle pattern that arises from burned tissue. Watts et al. showed the correlation between the laser Doppler measurements of dermal flow and the burn depth measured based on the histological assessment [14

14. A.M.I. Watts, M.P.H. Tyler, M.E. Perry, A.H.N. Roberts, and D.A. McGrouther, “Burn depth and its histological measurement,” Burns 27, 154–160 (2001). [CrossRef] [PubMed]

]. In the histological measurements, they studied the interface of the deepest blocked vessels and the most superficial patent as a demarcation of damaged and healthy tissue. However, burn wounds are much more complex than the two-layer model used in our study, and the perfusion in the tissue layer underlying the burn will cause less dynamic light scattering than the particle suspension used in this study. Hence, practical use of this approach in determining the depth of burn requires further in vivo measurements.

5. Conclusion

We aim to apply this method to perform path length resolved perfusion measurements in skin, which will overcome the inherent limitation of conventional LDPM that restrict its clinical usefulness, where perfusion signal depends on photon path length. In contrast to Doppler OCT systems measuring singly scattered photons [15

15. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science 254, 1178–1181 (1991). [CrossRef] [PubMed]

], we aim for multiply scattered light having traveled few millimeters through the tissue, and we quantitatively analyze the path length dependent Doppler broadening in a non-confocal two-fiber configuration as used in conventional LDPMs. Single-mode fiber optic Michelson interferometric systems used in Doppler OCT [15

15. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science 254, 1178–1181 (1991). [CrossRef] [PubMed]

] adopt on axis back reflection, whereas we use a multimode fiber-optic Mach-Zehnder interferometer, with positions for illumination and detection separated by a distance of ten to a hundred times the scattering mean free path length. This scheme offers greater flexibility, since the distance between the illumination and the detection fiber can be varied, giving one some control over the relative depth sensitivity of the system. Thus our measurements explore the regime of multiply scattered photons. In a normal skin, less than 10% of the photons are Doppler shifted due to the relative lower amount of red blood cells in the total skin volume. However, in richly perfused tissues this fraction of Doppler shifted photons increases due to the higher amount of red blood cells present in the skin capillary network. The results presented here show that this approach has potential applications in discriminating between statically and dynamically scattered light in the perfusion signal. Also, the path length resolved Doppler approach presented here will enable to correctly interpret the inter-and intra-individual variations in the LDF readings introduced by the variance in individual photon path lengths. Another important feature of this approach is the tunable depth resolved perfusion information that can be achieved. By translating the optical path length in the reference arm, the photons migrated deeper into the tissue can be made to interfere with the reference light and thus enable to discriminate between perfusion signal from superficial and deeper layers of tissue. Determination of superficial burn depth may be an important application of our technique.

Acknowledgements

References and links

1.

A. P. Shepherd and P. Å. öberg, Laser-Doppler Blood Flowmetry (Kluwer Academic, Boston, 1990).

2.

M. Larsson, W. Steenbergen, and T. Strömberg, “Influence of optical properties and fibre separation on laser Doppler flowmetry,” J. Biomed. Opt. 7, 236–243 (2001). [CrossRef]

3.

A. Liebert, M. Leahy, and R. Maniewski, “Multichannel laser-Doppler probe for blood perfusion measurements with depth discrimination,” Med. Bio. Eng. Comp. 36, 740–747 (1998). [CrossRef]

4.

B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Path length resolved measurements of multiple scattered photons in static and dynamic turbid media using phase modulated low coherence interferometry,” J. Biomed. Opt. 12, 024020 (2007). [CrossRef] [PubMed]

5.

B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, “Quantification of optical Doppler broadening and optical path lengths of multiply scattered light by phase modulated low coherence interferometry,” Opt. Express. 15(15), 9157–9165 (2007). [CrossRef] [PubMed]

6.

A. Wax, C. Yang, R.R. Dasari, and M.S. Feld, “Path-length-resolved dynamic light scattering: modeling the transition from single to diffusive scattering,” Appl. Opt. 40, 4222–4227 (2001). [CrossRef]

7.

K.K. Bizheva, A.M. Siegel, and D.A. Boas, “Path-length-resolved dynamic light scattering in highly scattering random media: The transition to diffusing wave spectroscopy,” Phys. Rev. E , 58 (1998). [CrossRef]

8.

A.L. Petoukhova, W. Steenbergen, and F.F.M. de Mul, “Path-length distribution and path-length resolved Doppler measurements of multiply scattered photons by use of low-coherence interferometer,” Opt. Lett. 26, 1492–1494 (2001). [CrossRef]

9.

A.L. Petoukhova, W. Steenbergen, T.G. van Leeuwen, and F. F. M. de Mul, “Effects of absorption on coherence domain path length resolved dynamic light scattering in the diffuse regime,” Appl. Phys. Lett. 81, 595–597 (2002). [CrossRef]

10.

J. C. Hebden, B.D. Price, A. P. Gibson, and G. Royle, “A soft deformable tissue-equivalent phantom for diffuse optical tomography,” Phys. Med. Bio. 51, 5581–5590 (2006). [CrossRef]

11.

A. G. Yodh, P. D. Kaplan, and D. J. Pine, “Pulsed diffusing-wave spectroscopy - high resolution through nonlinear optical gating,” Phys. Rev. B. 42, 4744–4747 (1990). [CrossRef]

12.

F. F. M. de Mul, M. H. Koelink, M. L. Kok, P. J. Harmsma, J. Greve, R. Graaff, and J. G. Aarnoudse, “Laser Doppler velocimetry and Monte Carlo simulations on models for blood perfusion in tissue,” Appl. Opt. 34, 6595–6611 (1995). [CrossRef] [PubMed]

13.

A. Sadhwani, K. T. Schomacker, G. J. Tearney, and N. S. Nishioka, “Determination of Teflon thickness with laser speckle. I. Potential for burn depth diagnosis,” Appl. Opt. 35, 5727–5735 (1996). [CrossRef] [PubMed]

14.

A.M.I. Watts, M.P.H. Tyler, M.E. Perry, A.H.N. Roberts, and D.A. McGrouther, “Burn depth and its histological measurement,” Burns 27, 154–160 (2001). [CrossRef] [PubMed]

15.

D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science 254, 1178–1181 (1991). [CrossRef] [PubMed]

OCIS Codes
(120.3180) Instrumentation, measurement, and metrology : Interferometry
(170.3340) Medical optics and biotechnology : Laser Doppler velocimetry
(170.3890) Medical optics and biotechnology : Medical optics instrumentation
(170.4580) Medical optics and biotechnology : Optical diagnostics for medicine

ToC Category:
Instrumentation, Measurement, and Metrology

History
Original Manuscript: July 18, 2007
Revised Manuscript: September 26, 2007
Manuscript Accepted: September 26, 2007
Published: September 28, 2007

Virtual Issues
Vol. 2, Iss. 11 Virtual Journal for Biomedical Optics

Citation
B. Varghese, V. Rajan, T. G. van Leeuwen, and W. Steenbergen, "Discrimination between Doppler-shifted and non-shifted light in coherence domain path length resolved measurements of multiply scattered light," Opt. Express 15, 13340-13350 (2007)
http://www.opticsinfobase.org/vjbo/abstract.cfm?URI=oe-15-20-13340


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. A. P. Shepherd and P. Å. Öberg, Laser-Doppler Blood Flowmetry (Kluwer Academic, Boston, 1990).
  2. M. Larsson, W. Steenbergen, and T. Strömberg, "Influence of optical properties and fibre separation on laser Doppler flowmetry," J. Biomed. Opt. 7, 236-243 (2001). [CrossRef]
  3. A. Liebert, M. Leahy and R. Maniewski, "Multichannel laser-Doppler probe for blood perfusion measurements with depth discrimination," Med. Bio. Eng. Comp. 36, 740-747 (1998).Q1 [CrossRef]
  4. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, "Path length resolved measurements of multiple scattered photons in static and dynamic turbid media using phase modulated low coherence interferometry," J. Biomed. Opt. 12, 024020 (2007). [CrossRef] [PubMed]
  5. B. Varghese, V. Rajan, T.G. Van Leeuwen, and W. Steenbergen, "Quantification of optical Doppler broadening and optical path lengths of multiply scattered light by phase modulated low coherence interferometry," Opt. Express. 15(15), 9157-9165 (2007). [CrossRef] [PubMed]
  6. A. Wax, C. Yang, R.R. Dasari, and M.S. Feld, "Path-length-resolved dynamic light scattering: modeling the transition from single to diffusive scattering," Appl. Opt. 40, 4222-4227 (2001). [CrossRef]
  7. K.K. Bizheva, A.M. Siegel, and D.A. Boas, "Path-length-resolved dynamic light scattering in highly scattering random media: The transition to diffusing wave spectroscopy," Phys. Rev. E,  58 (1998).Q2 [CrossRef]
  8. A.L. Petoukhova, W. Steenbergen, and F.F.M. de Mul, "Path-length distribution and path-length resolved Doppler measurements of multiply scattered photons by use of low-coherence interferometer," Opt. Lett. 26, 1492-1494 (2001). [CrossRef]
  9. A.L. Petoukhova, W. Steenbergen, T.G. van Leeuwen and F. F. M. de Mul, "Effects of absorption on coherence domain path length resolved dynamic light scattering in the diffuse regime," Appl. Phys. Lett. 81, 595-597 (2002). [CrossRef]
  10. J. C. Hebden, B.D. Price, A. P. Gibson and G. Royle, "A soft deformable tissue-equivalent phantom for diffuse optical tomography," Phys. Med. Bio. 51, 5581-5590 (2006).Q3 [CrossRef]
  11. A. G. Yodh, P. D. Kaplan, and D. J. Pine, 'Pulsed diffusing-wave spectroscopy - high resolution through nonlinear optical gating,' Phys. Rev. B. 42, 4744-4747 (1990). [CrossRef]
  12. F. F. M. de Mul, M. H. Koelink, M. L. Kok, P. J. Harmsma, J. Greve, R. Graaff, and J. G. Aarnoudse, "Laser Doppler velocimetry and Monte Carlo simulations on models for blood perfusion in tissue," Appl. Opt. 34, 6595-6611 (1995). [CrossRef] [PubMed]
  13. A. Sadhwani, K. T. Schomacker, G. J. Tearney, and N. S. Nishioka, "Determination of Teflon thickness with laser speckle. I. Potential for burn depth diagnosis," Appl. Opt. 35, 5727-5735 (1996). [CrossRef] [PubMed]
  14. A.M.I. Watts, M.P.H. Tyler, M.E. Perry, A.H.N. Roberts, D.A. McGrouther, "Burn depth and its histological measurement," Burns 27, 154-160 (2001). [CrossRef] [PubMed]
  15. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, "Optical coherence tomography," Science 254, 1178-1181 (1991). [CrossRef] [PubMed]

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.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited