Home | About | Journals | Submit | Contact Us | Français |

**|**Opt Express**|**PMC3601731

Formats

Article sections

- Abstract
- 1. Introduction
- 2. Materials and methods
- 3. Results
- 4. Discussion
- 5. Conclusion
- References and links

Authors

Related links

Optics Express

Opt Express. 2012 September 24; 20(20): 22262–22277.

Published online 2012 September 13. doi: 10.1364/OE.20.022262

PMCID: PMC3601731

Jonghwan Lee,^{
1
,}^{
*
} Weicheng Wu,^{
1
} James Y. Jiang,^{
2
} Bo Zhu,^{
3
} and David A. Boas^{
1
}

Received 2012 June 15; Revised 2012 August 18; Accepted 2012 August 24.

Copyright ©2012 Optical Society of America

This article has been cited by other articles in PMC.

We introduce an integration of dynamic light scattering (DLS) and optical coherence tomography (OCT) for high-resolution 3D imaging of heterogeneous diffusion and flow. DLS analyzes fluctuations in light scattered by particles to measure diffusion or flow of the particles, and OCT uses coherence gating to collect light only scattered from a small volume for high-resolution structural imaging. Therefore, the integration of DLS and OCT enables high-resolution 3D imaging of diffusion and flow. We derived a theory under the assumption that static and moving particles are mixed within the OCT resolution volume and the moving particles can exhibit either diffusive or translational motion. Based on this theory, we developed a fitting algorithm to estimate dynamic parameters including the axial and transverse velocities and the diffusion coefficient. We validated DLS-OCT measurements of diffusion and flow through numerical simulations and phantom experiments. As an example application, we performed DLS-OCT imaging of the living animal brain, resulting in 3D maps of the absolute and axial velocities, the diffusion coefficient, and the coefficient of determination.

Dynamic light scattering (DLS) is widely used to quantify the dynamics of scattering particles by analyzing the autocorrelation function of light scattered from the particles [1–3]. Optical coherence tomography (OCT) uses coherence gating to enable the measurement and analysis of light back-scattered from small volumes (typically on the order of a few wavelengths), facilitating high-resolution 3D imaging of the sample structure [4–6]. Therefore, an integration of these two techniques would provide high-resolution 3D imaging of particle dynamics (e.g., a 3D map of the diffusion coefficient quantifying the diffusive motion of particles at every voxel).

This paper describes the integration of DLS and OCT leading to a novel technique for high-resolution 3D imaging of particle dynamics, named DLS-OCT. For this integration, we propose and validate a theory that considers three concerns simultaneously. First, as OCT provides a very small probing volume, the theory considers the effect of such a wavelength-scale probing volume on the autocorrelation function of scattered light. Second, the theory considers the non-ergodic effect, since static and moving particles can be mixed within the OCT resolution volume in a highly heterogeneous sample. Finally, as both diffusion and flow can appear in a sample though spatially separated (e.g., the brain where both translational blood flow and diffusive intracellular organelle motions appear), the theory should be able to distinguish diffusive motions from translational flow and measure either the diffusion coefficient or the velocity according to the distinction. Although each of these concerns on DLS have been studied separately in the literature: the effect of the finite sample volume [7–9], the non-ergodic effect [10,11], and the mixture of diffusion and flow [12]; no theory has considered all of these simultaneously and has related these to OCT. Therefore, we derived the field autocorrelation function directly from the phase-resolved OCT signal under the assumption that static and moving particles are mixed in the resolution volume of OCT and the moving particles can exhibit either diffusive or translational motion.

This paper describes the derivation of the DLS-OCT theory, from the OCT signal when a single particle moves, to the field autocorrelation function when many particles exhibit heterogeneous dynamics. We propose a fitting algorithm to accurately estimate the diffusion coefficient and velocity from the autocorrelation function. The DLS-OCT theory has been validated through numerical simulations starting with the numerical position data of particles. Further, measurements of the axial and transverse velocities and the diffusion coefficient have been validated with phantom experiments. Based on these validations, we present an example application where DLS-OCT imaging was performed on the living rodent brain cortex, resulting in a velocity map showing blood flow and a diffusion map showing non-translational motions at the boundary of vessels.

We used a spectral-domain optical coherence tomography (SD-OCT) system (Thorlabs, Inc) and optimized it for DLS-OCT imaging. The details of the OCT system are described in our previous publication [13]. SD-OCT generally has several advantages over time-domain OCT, one of which is to provide phase information that is very sensitive to the motion of particles [14,15]. The system employs a large-bandwidth near-infrared light source (1310 ± 170 nm in wavelength) for a large imaging depth and high spatial resolution (an axial resolution of 3.5 μm in biological tissue). The transverse resolution depends on which objective lens is used. Our 10X objective lens was used for the phantom and animal experiments described in this paper, leading to the transverse resolution of 3.5 μm in biological tissue. The scanning speed is 47,000 A-scans/s.

A detailed description of the animal preparation is given in our previous publication [13]. All experimental procedures were approved by the Massachusetts General Hospital Subcommitee on Research Animal Care. Sprague Dawley rats (250-300 g) were anesthetized with isofluorane and ventilated with a mixture of air and oxygen during surgical procedures. A craniotomy was performed, and a 3 mm × 3 mm area over the somatosensory cortex was exposed. The dura was carefully removed, and then the brain surface was covered with agarose gel and a glass cover slip. Physiological signs such as heart rate, body temperature and blood pressure were continuously monitored during surgery and during the experiment.

Our derivation starts from the SD-OCT signal when only one particle exists in the sample. SD-OCT generally measures the spectrum of the interference between the reference beam and light reflected from the sample. With reasonable approximations and the assumption of Gaussian-shaped source light spectrum, the Fourier transform of the interference spectrum results in the signal arising from a voxel (i.e., the complex-valued field reflectivity) where a single particle moves as *z*_{1}(*t*) along the z axis with respect to the center of the voxel:

$${R}_{1}(t)={R}_{01}{e}^{-{\scriptscriptstyle \frac{1}{2}}{(\Delta q)}^{2}{z}_{1}^{2}(t)}{e}^{iq{z}_{1}(t)}$$

(1)

where the amplitude *R*_{01} is proportional to the field reflectivity of the particle. The other quantities are defined in Table 1
. The spectral width of the light source Δ*q* determines the axial length of the coherence gating of SD-OCT (i.e. the Gaussian function in the magnitude of Eq. (1)). The transverse resolution is generally determined by the objective lens used in the focusing optics. Therefore, the signal from a particle moving in the three-dimensional space is

$${R}_{1}(t)={R}_{01}{e}^{-2{h}_{t}^{2}[{x}_{1}^{2}(t)+{y}_{1}^{2}(t)]-2{h}^{2}{z}_{1}^{2}(t)}{e}^{iq{z}_{1}(t)}.$$

(2)

Now, we derive the field autocorrelation function of the OCT signal (Eq. (2)) when a single particle moves in the voxel. The movement of a particle can be defined by the self-part of the Van Hove space-time correlation function [16], that is, the probability that the particle which was initially at position **r _{i}** at time

$${P}_{1}({r}_{f},{t}_{f}|{r}_{i},{t}_{i})=\frac{1}{2\sqrt{\pi D({t}_{f}-{t}_{i})}}{e}^{-\frac{{\left|{r}_{f}-{r}_{i}-v({t}_{f}-{t}_{i})\right|}^{2}}{4D({t}_{f}-{t}_{i})}}.$$

(3)

Then, the field autocorrelation function of the signal is

$$\begin{array}{c}g(\tau )=E\left[\frac{{\u3008{R}_{1}^{*}(t){R}_{1}(t+\tau )\u3009}_{t}}{{\u3008{R}_{1}^{*}(t){R}_{1}(t)\u3009}_{t}}\right]\\ =\frac{1}{1+4{h}_{t}^{2}D\tau}\frac{1}{\sqrt{1+4{h}^{2}D\tau}}{e}^{-{h}_{t}^{2}\frac{{({v}_{t}\tau )}^{2}}{1+4{h}_{t}^{2}D\tau}-{h}^{2}\frac{{({v}_{z}\tau )}^{2}}{1+4{h}^{2}D\tau}}{e}^{-{q}^{2}\frac{D\tau}{1+4{h}^{2}D\tau}}{e}^{iq\frac{{v}_{z}\tau}{1+4{h}^{2}D\tau}}\end{array}$$

(4)

where *E*[ ] means the average over random initial positions of the particle and
${v}_{t}^{2}={v}_{x}^{2}+{v}_{y}^{2}$ is the squared transverse speed. Here, we assumed that the initial position of the particle is evenly distributed through the whole space. In the regime of this study, where the maximum time lag *τ*_{max} = 25 Δ*t* limits the maximum measureable velocity to 1/*h*/*τ*_{max} ~5 mm/s and the corresponding maximum measurable diffusion coefficient to (*h*/*q*)^{2}*v*_{max}^{2}*τ*_{max} ~10 μm^{2}/s, the factors that couple the diffusive movement and the small voxel size,
$1/(1+4{h}_{t}^{2}D\tau )$ and
$1/(1+4{h}^{2}D\tau )$, are very close to 1 (the minimum was 0.9969). Therefore, Eq. (4) becomes

$$g(\tau )={e}^{-{h}_{t}^{2}{v}_{t}^{2}{\tau}^{2}-{h}^{2}{v}_{z}^{2}{\tau}^{2}}{e}^{-{q}^{2}D\tau}{e}^{iq{v}_{z}\tau}.$$

(5)

When there are many particles exhibiting an identical mixture of translational and diffusive motions in the voxel, the OCT signal can be expressed by the superposition of signals from the individual particles:

$$R(t)={\displaystyle \sum _{j=1}^{N}{R}_{j}(t)}={\displaystyle \sum _{j=1}^{N}{R}_{0j}{e}^{-2{h}_{t}^{2}[{x}_{j}^{2}(t)+{y}_{j}^{2}(t)]-2{h}^{2}{z}_{j}^{2}(t)}{e}^{iq{z}_{j}(t)}}$$

(6)

where *N* is the number of particles. The field autocorrelation function of this signal is the same as Eq. (5),

$$g(\tau )={e}^{-{h}_{t}^{2}{v}_{t}^{2}{\tau}^{2}-{h}^{2}{v}_{z}^{2}{\tau}^{2}}{e}^{-{q}^{2}D\tau}{e}^{iq{v}_{z}\tau}$$

(7)

where we used the assumption that (i) the movement of the j-th particle is uncorrelated with that of the k-th particle so that $P({r}_{k}^{\text{'}},t+\tau |{r}_{j},t)=0$; and (ii) the magnitude of reflectivity does not markedly vary across particles so that ${\left|{R}_{0j}\right|}^{2}\cong {\overline{R}}_{0}^{2}$. It is not surprising that the field autocorrelation function of many homogenous particles is close to that of a single particle because the latter has already been averaged over random initial particle positions.

When particles exhibit heterogeneous dynamics in the voxel, those dynamics can be categorized into three groups (static (S), diffusing or flowing (F), and entering or exiting (E); see Fig. 1(a)
). We now assume an isotropic voxel (*h _{t}* =

$$R(t)={\displaystyle \sum _{j=1}^{{N}_{S}}{R}_{S0j}{e}^{-2{h}^{2}{r}_{Sj}^{2}}{e}^{iq{z}_{Sj}}}+{\displaystyle \sum _{j=1}^{{N}_{F}}{R}_{F0j}{e}^{-2{h}^{2}{r}_{Fj}^{2}(t)}{e}^{iq{z}_{Fj}(t)}}+{\displaystyle \sum _{j=1}^{{N}_{E}}{R}_{E0j}{e}^{-2{h}^{2}{r}_{Ej}^{2}(t)}{e}^{iq{z}_{Ej}(t)}}$$

(8)

where *N*_{Ω} is the number of Ω-type particles in the voxel (Ω = S, F, or E). Here, the *r _{Ej}*(

$$R(t)={\displaystyle \sum _{j=1}^{{N}_{S}}{e}^{-2{h}^{2}{r}_{Sj}^{2}}{e}^{iq{z}_{Sj}}}+{\displaystyle \sum _{j=1}^{{N}_{F}}{e}^{-2{h}^{2}{r}_{Fj}^{2}(t)}{e}^{iq{z}_{Fj}(t)}}+{\displaystyle \sum _{j=1}^{{N}_{E}}{e}^{-2{h}^{2}{r}_{Ej}^{2}(t)}{e}^{iq{z}_{Ej}(t)}}\equiv {R}_{S}+{R}_{F}(t)+{R}_{E}w(t)$$

(9)

where *w*(*t*) is complex-valued white noise. In practice, the third term may include stochastic fluctuations in the OCT signal due to imperfections of the OCT system.

Conceptual illustration of the DLS-OCT theory. (A) Particles within the OCT resolution volume can be categorized into three groups: static, flowing or diffusing, and entering or exiting particles. For clarification, entering/exiting particles enter into **...**

The field autocorrelation function of this OCT signal has several groups of coupled terms. Note that static and moving particles are mixed in this case (i.e., non-ergodic) and thus the system is not stationary; that is, $\u3008{R}^{*}(t)R(t+\tau )\u3009\ne {R}^{*}(0)R(\tau )$.

$$\begin{array}{l}\u3008{R}^{*}(t)R(t+\tau )\u3009\\ =\u3008[{R}_{S}^{*}+{R}_{F}^{*}(t)+{R}_{E}^{*}{w}^{*}(t)][{R}_{S}+{R}_{F}(t+\tau )+{R}_{E}w(t+\tau )]\u3009\\ ={\left|{R}_{S}\right|}^{2}+\u3008{R}_{F}^{*}(t){R}_{F}(t+\tau )\u3009+{\left|{R}_{E}\right|}^{2}\u3008{w}^{*}(t)w(t+\tau )\u3009(\text{self}\text{coupled}\text{terms})\\ +{R}_{S}^{*}\u3008{R}_{F}(t+\tau )\u3009+{R}_{S}\u3008{R}_{F}^{*}(t)\u3009(\mathrm{S}-\mathrm{F}\text{coupled}\text{terms})\\ +{R}_{S}^{*}{R}_{E}\u3008w(t+\tau )\u3009+{R}_{S}{R}_{E}^{*}\u3008{w}^{*}(t)\u3009(\mathrm{S}-\mathrm{E}\text{coupled}\text{terms})\\ +{R}_{E}\u3008{R}_{F}^{*}(t)w(t+\tau )\u3009+{R}_{E}^{*}\u3008{R}_{F}(t+\tau ){w}^{*}(t)\u3009(\mathrm{F}-\mathrm{E}\text{coupled}\text{terms})\end{array}$$

(10)

Each coupled term leads to

$$\begin{array}{l}E[{\left|{R}_{S}\right|}^{2}]={N}_{S}\frac{\pi \sqrt{\pi}}{8{h}^{3}}\\ E[\u3008{R}_{F}^{*}(t){R}_{F}(t+\tau )\u3009]={N}_{F}\frac{\pi \sqrt{\pi}}{8{h}^{3}}{e}^{-{h}^{2}{v}^{2}{\tau}^{2}}{e}^{-{q}^{2}D\tau}{e}^{iq{v}_{z}\tau}\\ E[{\left|{R}_{E}\right|}^{2}\u3008{w}^{*}(t)w(t+\tau )\u3009]={N}_{E}\delta (\tau )\\ E[{R}_{S}^{*}\u3008{R}_{F}(t+\tau )\u3009]={N}_{S}\frac{\pi \sqrt{\pi}}{2\sqrt{2}{h}^{3}}{e}^{-\frac{3}{8}{(\frac{q}{h})}^{2}}\left\{{\displaystyle \sum _{k=1}^{{N}_{F}}\frac{1}{T}{\displaystyle \underset{\tau}{\overset{\tau +T}{\int}}dt\text{'}}}{\displaystyle \prod _{j=x,y,z}\frac{{e}^{-\frac{{[qDt\text{'}-{\scriptscriptstyle \frac{i}{2}}({r}_{jFk0}+{v}_{j}t\text{'})]}^{2}}{8{h}^{2}{D}^{2}t{\text{'}}^{2}+Dt\text{'}}}}{\sqrt{1+8{h}^{2}Dt\text{'}}}}\right\}\underset{{(\frac{q}{h})}^{2}>>1}{\to}0\\ E[{R}_{S}\u3008{R}_{F}^{*}(t)\u3009]\underset{{(\frac{q}{h})}^{2}>>1}{\to}0\\ E[{R}_{S}^{*}{R}_{E}\u3008w(t+\tau )\u3009]=E[{R}_{S}{R}_{E}^{*}\u3008{w}^{*}(t)\u3009]=0\\ E[{R}_{E}\u3008{R}_{F}^{*}(t)w(t+\tau )\u3009]={N}_{S}{N}_{E}\frac{{\pi}^{3}}{8{h}^{6}}{e}^{-\frac{3}{4}{(\frac{q}{h})}^{2}}\underset{{(\frac{q}{h})}^{2}>>1}{\to}0\\ E[{R}_{E}^{*}\u3008{R}_{F}(t+\tau ){w}^{*}(t)\u3009]=\underset{{(\frac{q}{h})}^{2}>>1}{\to}0\end{array}$$

(11)

where *δ*(*τ*) is the delta function, *r _{jFk0}* is the

$$P({r}_{Fk},t+\tau |{r}_{Sm},t)=\frac{P[({r}_{Fk},t+\tau )\cap ({r}_{Sm},t)]}{P({r}_{Sm},t)}=P({r}_{Fk},t+\tau ).$$

(12)

Note that the S-F coupled terms do not necessarily vanish when the particle movements are uncorrelated. For the terms to vanish, the deviation in the probing light spectrum should be relatively small, (*q*/Δ*q*)^{2} >> 1, which is satisfied in this study as (*q*/Δ*q*)^{2} > 300. Consequently, the field autocorrelation function of mixed dynamics is given by

$$g(\tau )={M}_{S}+{M}_{F}{e}^{-{h}_{t}^{2}{v}_{t}^{2}{\tau}^{2}-{h}^{2}{v}_{z}^{2}{\tau}^{2}}{e}^{-{q}^{2}D\tau}{e}^{iq{v}_{z}\tau}+(1-{M}_{S}-{M}_{F})\delta (\tau ).$$

(13)

In summary, under the assumption that static and moving particles are mixed in the OCT resolution volume and the moving particles can exhibit either translational or diffusive motion, the 4D (space and timelag) field autocorrelation function of the 4D (space and time) SD-OCT signal can be described as

$$\begin{array}{c}g(r,\tau )=\frac{{\u3008{R}^{*}(r,t)R(r,t+\tau )\u3009}_{t}}{{\u3008{R}^{*}(r,t)R(r,t)\u3009}_{t}}\\ ={M}_{S}(r)+{M}_{F}(r){e}^{-{h}_{t}^{2}{v}_{t}^{2}(r){\tau}^{2}-{h}^{2}{v}_{z}^{2}(r){\tau}^{2}}{e}^{-{q}^{2}D(r)\tau}{e}^{iq{v}_{z}(r)\tau}\\ +[1-{M}_{S}(r)-{M}_{F}(r)]\delta (\tau ).\end{array}$$

(14)

*M _{S}*(

The general behavior of this autocorrelation function is illustrated in Fig. 1(b). For a given voxel, the constant term (*M _{S}*) is the center of rotation, and

We developed and optimized an algorithm for estimating the dynamic parameters through fitting the theory to the autocorrelation data. The fitting has to determine five independent coefficients (*M _{S}*,

$${R}^{2}=1-\frac{\u3008{\left|g(\tau )-{M}_{S}-{M}_{F}{e}^{-{h}_{t}^{2}{v}_{t}^{2}{\tau}^{2}-{h}^{2}{v}_{z}^{2}{\tau}^{2}}{e}^{-{q}^{2}D\tau}{e}^{iq{v}_{z}\tau}-(1-{M}_{S}-{M}_{F})\delta (\tau )\right|}^{2}\u3009}{\u3008{\left|g(\tau )-\u3008g(\tau )\u3009\right|}^{2}\u3009}.$$

(15)

As shown in Fig. 2
, we found it inaccurate to use a simple fitting algorithm containing a single simplex search method [17] to find a global minimum in the five-dimensional parameter space. For this reason, we developed an advanced fitting algorithm. First, we reduced the number of parameters by determining *D*, *v _{z}* and

$$\begin{array}{l}{v}_{z}=\frac{1}{q}\frac{{\u3008\mathrm{arg}[g(\tau )-{M}_{S}]\tau \left|g(\tau )-{M}_{S}\right|\u3009}_{\tau 0}}{{\u3008{\tau}^{2}\left|g(\tau )-{M}_{S}\right|\u3009}_{\tau 0}}\\ {g}_{1}(\tau )=\frac{1}{{M}_{F}}\left|g(\tau )-{M}_{S}\right|{e}^{{h}^{2}{v}_{z}^{2}{\tau}^{2}}\\ \left[\begin{array}{cc}{\u3008{\tau}^{4}{g}_{1}(\tau )\u3009}_{\tau 0}& {\u3008{\tau}^{3}{g}_{1}(\tau )\u3009}_{\tau 0}\\ {\u3008{\tau}^{3}{g}_{1}(\tau )\u3009}_{\tau 0}& {\u3008{\tau}^{2}{g}_{1}(\tau )\u3009}_{\tau 0}\end{array}\right]\left[\begin{array}{c}-{h}_{t}^{2}{v}_{t}^{2}\\ -{q}^{2}D\end{array}\right]=\left[\begin{array}{c}{\u3008{\tau}^{2}{g}_{1}(\tau )\mathrm{ln}{g}_{1}(\tau )\u3009}_{\tau 0}\\ {\u3008\tau {g}_{1}(\tau )\mathrm{ln}{g}_{1}(\tau )\u3009}_{\tau 0}\end{array}\right].\end{array}$$

(16)

During this determination, only the nonzero timelag data were used to exclude the term of the entering/exiting dynamics, (1-*M _{S}*-

The performance of the simple fitting algorithm for an isotropic voxel (A) and an anisotropic voxel (B). *h*_{t} = 4*h* was used for the anisotropic voxel. The performance was tested for a total of 25,200 combinations (36 number densities, 10 diffusion coefficients, **...**

Even when fitting for the two remaining parameters (*M _{S}* and

For the initial guesses of *M _{S}* and

$$\begin{array}{l}\left[\begin{array}{ccc}{\tau}_{1}^{2}& {\tau}_{1}& 1\\ {\tau}_{2}^{2}& {\tau}_{2}& 1\\ {\tau}_{3}^{2}& {\tau}_{3}& 1\end{array}\right]\left[\begin{array}{c}{c}_{1}\\ {c}_{2}\\ {c}_{3}\end{array}\right]=\left[\begin{array}{c}\left|g({\tau}_{1})\right|\\ \left|g({\tau}_{2})\right|\\ \left|g({\tau}_{3})\right|\end{array}\right]\\ {M}_{E}=1-{c}_{3}\end{array}$$

(17)

where *τ _{k}* =

For the first initial guess of *M _{S}*, the center of rotation was used. We derived an equation to find the center of rotation (

$$\begin{array}{l}\left[\begin{array}{cc}\u3008{a}_{k}^{2}\u3009-{\overline{a}}^{2}& \u3008({a}_{k}-\overline{a}){b}_{k}\u3009\\ \u3008({b}_{k}-\overline{b}){a}_{k}\u3009& \u3008{b}_{k}^{2}\u3009-{\overline{b}}^{2}\end{array}\right]\left[\begin{array}{c}A\\ B\end{array}\right]=\left[\begin{array}{c}\u3008{a}_{k}^{3}\u3009-\overline{a}\u3008{a}_{k}^{2}\u3009+\u3008({a}_{k}-\overline{a}){b}_{k}^{2}\u3009\\ \u3008{b}_{k}^{3}\u3009-\overline{b}\u3008{b}_{k}^{2}\u3009+\u3008({b}_{k}-\overline{b}){a}_{k}^{2}\u3009\end{array}\right]\\ {M}_{S}=A\end{array}$$

(18)

where < > is the average over *k* > *N _{τ}*/2 and

For the second initial guess of *M _{S}* and

For each initial guess, the simplex search method was used to find *M _{S}* and

We tested if the algorithm results in better estimation of the parameters than the simple algorithm shown in Fig. 2. Autocorrelation data were generated for various values of *M _{S}*,

Our DLS-OCT theory has been validated using numerical simulations. As summarized in Fig. 5
, we numerically generated two-dimensional (transverse and axial axes) position data of 100 particles for 100 time steps (Δ*t* = 20 μs) for various parameter values (Table 2
). Position data were used to generate SD-OCT signals, and autocorrelation data were obtained from the signals. The autocorrelation data were fit using the algorithm described in the above section, leading to the estimation of the dynamic parameters.

Process of numerical simulation for validating the DLS-OCT theory. True values of the diffusion coefficient and velocity were determined by fitting the mean square displacement (MSD) of the numerical position data.

First, we validated the estimation of the flow velocity for various parameters while fixing the diffusion coefficient to zero. As can be seen in Fig. 6(a)
, the autocorrelation function behaved as our theory predicted across several combinations of the parameters. Note that the center of rotation moved when the number of static particles changed. In this study, we mainly investigated the diffusion coefficient and velocities (*D*, *v*, and *v _{z}*) multiplied by

Numerical validation of the velocity estimation. (A) Examples of the autocorrelation functions obtained from the numerical position data and their fitting results. *D* = 0 μm^{2}/s. (B) Results of the estimation of the absolute velocity and the axial **...**

Next, we validated the estimation of the diffusion coefficient. The position data of diffusing particles was generated with the fractional Brownian motion synthesis algorithm [18]. Our model was fit to the numerical autocorrelation function data (Fig. 7(a)
), and the *M _{F}*-weighted diffusion coefficient was estimated close to the true values (Fig. 7(b)). In this simulation, the absolute velocity was estimated negligible (0.24 ± 0.47 mm/s).

Numerical validation of the diffusion estimation. (A) Examples of decays of the autocorrelation functions obtained from the numerical position data and their fitting results. *v* = 0 mm/s and *v*_{z} = 0 mm/s. (B) Results of the estimation of the diffusion coefficient. **...**

Results from these numerical simulations confirm that the DLS-OCT theory enables us to effectively distinguish diffusion from flow for a given autocorrelation function and accurately estimate the diffusion coefficient or the absolute and axial velocities according to the distinction. The distinction is mainly attributed to the difference in the waveforms of decay of the *M _{F}*-term (Fig. 7(c)) as predicted in Eq. (14).

DLS-OCT measurement of the velocity and diffusion coefficient was further validated through phantom experiments. We used the SD-OCT system described in the Materials and methods. In order to verify measurement of the flow velocity, a piezoelectrically actuated static sample was used to simulate axial movements of particles while transverse movements were implemented by galvanometric lateral scanning of OCT beam. For each combination of the axial and transverse velocities, we acquired the OCT signal and obtained 3D autocorrelation function data, *g*(*x*,*z*,*τ*). The autocorrelation function was averaged over neighboring 3 × 3 voxels and then was used to estimate the dynamic parameters, leading to 2D maps of *M _{S}*(

Based on the validated DLS-OCT measurement of the velocity and diffusion, we performed DLS-OCT imaging of the living rodent cortex. A-scans were repeated 100 times at a fixed position, and the position was moved in a raster manner to scan the volume of 600 μm × 600 μm × 300 μm (*x* × *y* × *z*) of the cortical surface. This FOV consisted of 400 × 400 positions, leading to the scanning step size of 1.5 μm and a total scanning time of ~10 min. We chose 100 A-scans per position because our fitting algorithm worked well when the autocorrelation function data had ≥25 time points, and the measurement time of ~2 ms corresponding to 100 A-scans was sufficiently shorter than the characteristic time constants of the primary sources of motion artifacts (i.e., cardiac and respiratory motions) as revealed in our previous study [13]. The 4D (space and time) complex-valued field reflectivity of the sample was obtained, and then was used to produce the 4D (space and timelag) autocorrelation function data. The autocorrelation function was averaged over 3 × 3 × 3 neighboring voxels, since the theory considers the autocorrelation function being the average over not only the measurement time but also various initial positions of particles. The spatial averaging of the autocorrelation function is based on the assumption that the initial position of the particles will vary across the neighboring voxels whereas their dynamics will be similar. For this assumption to be as valid as possible, we chose the scanning step size slightly smaller than the conventional Nyquist sampling frequency. Analysis of the 4D autocorrelation function data led to 3D maps of the transverse and axial velocities, the diffusion coefficient, and the coefficient of determination (*R*^{2}). This processing took ~3 hours with our 200-core parallel processing cluster system.

As shown in Fig. 9(a) , the absolute velocity map clearly revealed the structure of the vascular anatomy and cerebral blood flow. The axial velocity map showed the axial component of the flow velocity, which looked very similar to conventional Doppler OCT images. The flow direction determined by the axial and transverse velocities agreed with the structural direction of vessels. This agreement between the structural and flow direction can be more clearly seen in the cross-sectional maps (Fig. 9(a), right).

DLS-OCT imaging of the living rodent cortex. (A) The first image presents the maximum projection (MP) of the 3D map of the absolute velocity along the depth (i.e., *en face*), and the second image presents the *en face* signed maximum projection (SMP) of **...**

When the diffusion map is overlaid with the absolute velocity map (Fig. 9(b)), we observed high-diffusion and low-velocity voxels at the boundary of vessels. Interestingly, the vessel boundaries also exhibited a characteristic low coefficient of determination (Fig. 9(c)). Single-plane images clearly showed the result that the vessel boundaries exhibit high diffusion, low velocity, and low *R*^{2}. At the circular cross-sections of the vessels (Fig. 9(c), blue arrows), high-velocity and high-*R*^{2} blood flow are surrounded by the high-diffusion, low-velocity, and low-*R*^{2} dynamics of the vessel boundaries. In particular, the low *R*^{2} means that the motion was neither translational nor diffusive; we hypothesize that it might be oscillatory due to the interaction between blood flow and the tension of vessel walls. As can be seen in the examples of the autocorrelation functions (Fig. 9(d)), the composition ratio of static particles (*M _{S}*; the center of rotation in the complex plane) increased as the voxel is located far from the center of the vessel. When the voxel is located at the vessel boundary, it exhibited the characteristic dynamics with a low

DLS-OCT has several advantages over conventional techniques. As for the application to blood flow imaging, laser Doppler flowmetry is used to measure blood flow at a fixed point [19,20], and its imaging corollary provides the 2D map of blood flow [21]. Doppler OCT enables 3D imaging of the axial flow velocity with microscopic resolution [22]. Compared to these techniques, DLS-OCT can simultaneously and independently measure the axial and transverse components of the flow velocity. Further, DLS-OCT can distinguish diffusion from flow, leading to a separate 3D map of the diffusion coefficient. In addition, DLS-OCT may offer additional contrast mechanisms (e.g., those based on *R*^{2} and *M _{F}*). The

Each term in the DLS-OCT model we derived is similar to those predicted in several DLS studies. The velocity-dependent decay term is similar to that predicted in the study of the effect of the finite sample volume on DLS analysis [7–9], although the autocorrelation function in our theory was directly derived from the OCT signal whereas the finite size of the sample volume in the literature resulted from the illuminating and collecting optics. The linear decomposition of the static component and the flowing/diffusing component is similar to the decomposition shown in the study of DLS in non-ergodic media [10,11]. The combined diffusion-oriented and flow-oriented decay terms were similarly introduced in the study of DLS where diffusion and flow are mixed [12]. Therefore, our theory can be understood as a mathematical combination of the phase-resolved OCT signal, the finite sample volume DLS model, the non-ergodic DLS model, and the model for the mixture of diffusion and flow.

Our DLS-OCT theory assumed that the composition ratios of the static and diffusing/flowing particles (*M _{S}* and

This study used the DLS-OCT theory for distinguishing whether the motion is translational or diffusive, not for measuring the mixture of the two motions. Nevertheless, the model seems to be able to measure the mixture of translational and diffusive motions. It was reported, however, that the diffusion is estimated inaccurately when mixed with large flow [12]. We also observed this tendency in a phantom experiment, where the diffusion coefficient was overestimated by ≈1.5 times when it was mixed with flow whose velocity is larger than >1 mm/s. In this study, the relative magnitudes of flow and diffusion can be estimated based on how much they contribute to the decay in the autocorrelation function. With the parameters used in this study, *v* = 0.1 and 1 mm/s approximately correspond to *D* = 0.004 and 0.4 μm^{2}/s, respectively. The coupling between diffusion and large flow will not be a critical problem in the studies investigating a sample where diffusion and flow are spatially separated. Meanwhile, considerable caution would be required if one wants to apply the present technique to the measurement of diffusion that is mixed with large flow within the resolution volume. On the other hand, the voxels with large blood flow in vessels often exhibited high diffusion, which can be either a real diffusive motion or decorrelation of the OCT signal that can be quantified by the exponential decay with the diffusion coefficient. This high diffusion in vessels might be attributed to the mixture of non-translational motion of red blood cells including rotation, turbulence in blood flow, and the variance in the velocity distribution within the resolution volume. Since this interpretation has not been yet validated, and since the present study used the DLS-OCT theory for determining whether the motion is translational or not, we overlaid the diffusion map with the velocity map as in Figs. 9(B) and 9(C). This overlay means that our analysis gave priority to flow over diffusion so that diffusion is of interest only at the voxels with low flow. Therefore, the coupling between diffusion and large flow was not an important concern at those vessel boundaries.

OCT uses coherence gating to collect light only scattered from the resolution volume and is known to effectively exclude multiply scattered photons [23]. However, strong multiple scattering can give rise to a distortion in the OCT signal, which often causes an undesired shadow of large vessels [13]. This multiple scattering also affected DLS-OCT estimation of dynamics at the voxels located beneath the surface vessels. This effect can be seen clearly by comparing the merged images in Fig. 9, where Fig. 9(b) presents the maximum projection over the whole depth while Fig. 9(c) only presents a single plane at the depth of 120 μm. High velocity and diffusion appeared in the shadow of the large surface vessels as shown in the single-plane image. This multiple scattering effect also can be found in the cross-section image of the velocity map of Fig. 9(a), where the absolute velocity was estimated large in the vessel shadow but the axial velocity was not. Although this multiple scattering would not cause a serious problem as one may generally build an *en face* image by including surface vessels as in Fig. 9(b), in the future it will be interesting to derive a DLS-OCT model that considers the effect of multiple scattering.

Future theoretical efforts should consider the possibility of measuring a mixture of translational and diffusive motions, and the effect of multiple scattering. Although there can be various modified versions of the DLS-OCT theory, the one described here will work well for 3D imaging of dynamics in a highly heterogeneous sample where static and moving particles can be mixed within the micrometer-scale resolution volume and the moving particles can exhibit either diffusion or flow, as long as it is used in the single-scattering regime.

We demonstrated a technique based on the integration of DLS and OCT for high-resolution 3D imaging of heterogeneous diffusion and flow. A theory for this purpose was proposed and validated with numerical simulations and phantom experiments. The DLS-OCT theory enabled us to simultaneously measure the axial and transverse velocities and the diffusion coefficient with the micrometer-scale resolution. We tested the utility of this technique by applying it to 3D *in vivo* imaging of translational blood flow and non-translational motion of the vessel boundary in the living animal brain cortex.

This study was supported by the NIH ( R01-NS057476, R01-EB000790, P01NS055104) and the AFOSR ( MFEL FA9550-07-1-0101).

1. Clark N. A., Lunacek J. H., Benedek G. B., “A study of Brownian motion using light scattering,” Am. J. Phys. 38(5), 575–585 (1970).10.1119/1.1976408 [Cross Ref]

2. Durian D. J., Weitz D. A., Pine D. J., “Multiple light-scattering probes of foam structure and dynamics,” Science 252(5006), 686–688 (1991).10.1126/science.252.5006.686 [PubMed] [Cross Ref]

3. Boas D. A., Yodh A. G., “Spatially varying dynamical properties of turbid media probed with diffusing temporal light correlation,” J. Opt. Soc. Am. A 14(1), 192–215 (1997).10.1364/JOSAA.14.000192 [Cross Ref]

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

5. Tearney G. J., Brezinski M. E., Bouma B. E., Boppart S. A., Pitris C., Southern J. F., Fujimoto J. G., “In vivo endoscopic optical biopsy with optical coherence tomography,” Science 276(5321), 2037–2039 (1997).10.1126/science.276.5321.2037 [PubMed] [Cross Ref]

6. Chen Y., Vuong L. N., Liu J., Ho J., Srinivasan V. J., Gorczynska I., Witkin A. J., Duker J. S., Schuman J., Fujimoto J. G., “Three-dimensional ultrahigh resolution optical coherence tomography imaging of age-related macular degeneration,” Opt. Express 17(5), 4046–4060 (2009).10.1364/OE.17.004046 [PMC free article] [PubMed] [Cross Ref]

7. Edwards R. V., Angus J. C., French M. J., “Spectral analysis of the signal from the laser Doppler flowmeter: time-independent systems,” J. Appl. Phys. 42(2), 837–850 (1971).10.1063/1.1660102 [Cross Ref]

8. Chowdhury D. P., Sorensen C. M., Taylor T. W., Merklin J. F., Lester T. W., “Application of photon correlation spectroscopy to flowing Brownian motion systems,” Appl. Opt. 23(22), 4149–4154 (1984).10.1364/AO.23.004149 [PubMed] [Cross Ref]

9. Taylor T. W., Sorensen C. M., “Gaussian beam effects on the photon correlation spectrum from a flowing Brownian motion system,” Appl. Opt. 25(14), 2421–2426 (1986).10.1364/AO.25.002421 [PubMed] [Cross Ref]

10. Pusey P. N., Van Megen W., “Dynamic light scattering by non-ergodic media,” Physica A 157(2), 705–741 (1989).10.1016/0378-4371(89)90063-0 [Cross Ref]

11. Joosten J. G. H., Geladé E. T. F., Pusey P. N., “Dynamic light scattering by nonergodic media: Brownian particles trapped in polyacrylamide gels,” Phys. Rev. A 42(4), 2161–2175 (1990).10.1103/PhysRevA.42.2161 [PubMed] [Cross Ref]

12. Leung A. B., Suh K. I., Ansari R. R., “Particle-size and velocity measurements in flowing conditions using dynamic light scattering,” Appl. Opt. 45(10), 2186–2190 (2006).10.1364/AO.45.002186 [PubMed] [Cross Ref]

13. Lee J., Srinivasan V., Radhakrishnan H., Boas D. A., “Motion correction for phase-resolved dynamic optical coherence tomography imaging of rodent cerebral cortex,” Opt. Express 19(22), 21258–21270 (2011).10.1364/OE.19.021258 [PMC free article] [PubMed] [Cross Ref]

14. de Boer J. F., Cense B., Park B. H., Pierce M. C., Tearney G. J., Bouma B. E., “Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography,” Opt. Lett. 28(21), 2067–2069 (2003).10.1364/OL.28.002067 [PubMed] [Cross Ref]

15. Leitgeb R., Hitzenberger C., Fercher A., “Performance of fourier domain vs. time domain optical coherence tomography,” Opt. Express 11(8), 889–894 (2003).10.1364/OE.11.000889 [PubMed] [Cross Ref]

16. Van Hove L., “Correlations in space and time and Born approximation scattering in systems of interacting particles,” Phys. Rev. 95(1), 249–262 (1954).10.1103/PhysRev.95.249 [Cross Ref]

17. Lagarias J. C., Reeds J. A., Wright M. H., Wright P. E., “Convergence properties of the Nelder-Mead simplex method in low dimensions,” Siam J Optimiz 9(1), 112–147 (1998).10.1137/S1052623496303470 [Cross Ref]

18. Abry P., Sellan F., “The wavelet-based synthesis for fractional Brownian motion proposed by F. Sellan and Y. Meyer: remarks and fast implementation,” Appl. Comput. Harmon. Anal. 3(4), 377–383 (1996).10.1006/acha.1996.0030 [Cross Ref]

19. Stern M. D., “In vivo evaluation of microcirculation by coherent light scattering,” Nature 254(5495), 56–58 (1975).10.1038/254056a0 [PubMed] [Cross Ref]

20. Dirnagl U., Kaplan B., Jacewicz M., Pulsinelli W., “Continuous measurement of cerebral cortical blood flow by laser-Doppler flowmetry in a rat stroke model,” J. Cereb. Blood Flow Metab. 9(5), 589–596 (1989).10.1038/jcbfm.1989.84 [PubMed] [Cross Ref]

21. Ances B. M., Greenberg J. H., Detre J. A., “Laser Doppler imaging of activation-flow coupling in the rat somatosensory cortex,” Neuroimage 10(6), 716–723 (1999).10.1006/nimg.1999.0510 [PubMed] [Cross Ref]

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

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |