Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Opt Lett. Author manuscript; available in PMC Jun 29, 2010.
Published in final edited form as:
PMCID: PMC2894034
Simplified laser-speckle-imaging analysis method and its application to retinal blood flow imaging
Haiying Cheng1 and Timothy Q. Duong2
Yerkes Imaging Center, Department of Neurology, Emory University, Atlanta, Georgia 30329, USA
Laser speckle imaging (LSI) is widely used to study blood flow at high spatiotemporal resolution. Several papers recently pointed out that the commonly used LSI equation involves an approximation that could result in incorrect data analysis. We investigated the impact of such an approximation and introduced a simplified analysis method to improve computation time. Flow phantom studies were performed for validation. Moreover, we demonstrated a novel LSI application by imaging blood flow of rat retinas under normal and physiologic-challenge conditions. Because blood-flow abnormality is implicated in many retinal diseases, LSI could provide valuable physiologic, and potentially diagnostic, information.
Laser speckle imaging (LSI) [1] can be used to image instantaneous velocity distribution in vivo at very high spatiotemporal resolution, an improvement over the single point techniques, such as laser Doppler flowmetry [2]. Laser speckle is an optical interference effect. The speckle pattern fluctuates if the illuminated area contains moving particles such as moving red blood cells. By integrating the intensity fluctuations of the speckle pattern over a finite time, information about the motion of the scattering particles could be derived. LSI has recently been applied to image changes in cerebral blood flow associated with focal brain ischemia and cortical spreading depression in rats [3]. LSI has since proven to be a cost-effective technique for measuring dynamic blood flow changes at very high spatiotemporal resolution [4,5]. However, several papers [6-8] recently pointed out that the commonly used LSI equation involves an approximation that could result in incorrect data analysis.
In this study, we investigated the contribution of such approximation and its impact on LSI data analysis and proposed a simplified LSI analysis method to speed up computation time. For validation, we performed flow phantoms experiments at different physiological flow rates and different camera exposure times. Moreover, we demonstrated a novel in vivo application by imaging blood flow of the rat retinas in which the animals breathed air or oxygen. The latter was used to modulate blood flow for testing sensitivity. LSI blood-flow index maps and physiologically induced percent-change maps were analyzed.
Speckle contrast (K) is defined in terms of the standard deviation (σs) and the mean local spatial speckle intensity fluctuations (left angle bracketIright angle bracket). Per Fercher and Briers' formalism, the relation between speckle contrast and spatial variance of the time-averaged speckle pattern is [1]
equation M1
where g(τ) is the normalized autocorrelation function, τ is the lag time, and T is the CCD camera exposure time. Assuming a Lorentzian spectrum, g(τ) =exp(−|τ|τc), where τc is the correlation time, Eq. (1a) can be rewritten as
equation M2
where Tc is proportional to the mean velocity [9] V and is the blood-flow index.
However, in Goodman's original master LSI equation [10], the speckle contrast is related to the spatial variance in the time-averaged speckle pattern by
equation M3
equation M4
equation M5
Note that the factor (1−τ/T) was approximated to unity in Eq. (1a) to simplify the integration. Assuming a Lorentzian spectrum, Eq. (2b) can be rewritten as [11]
equation M6
It is unknown how such an approximation impacts LSI data analysis [6-8].
A typical laser diode in LSI experiments [3-5,7,12] has a wavelength spread Δλ>0.1 nm and a mean wavelength [lambda with macron]~780 nm, which yields a coherence time τc=[lambda with macron]2/(cΔλ) of ~10−12 s. For T=1–40 ms, Tc is >109. Substituting x for Tc and y for 1/K2, Eq. (1b) becomes
equation M7
and Eq. (2c) becomes
equation M8
Because 1/K2 is measured, the blood-flow index Tc calculated from the asymptotic approximation (x → ∞) of Method I is half of that of Method II. However, the percent changes between a stimulation and basal condition calculated using Eqs. (3) and (4) are identical. This analysis shows that the approximation of (1−τ /T) to unity is valid for calculating blood-flow change when Tc is large. Figure 1 graphically demonstrates this point and the asymptotic approximation is valid for T/τ >100. This approximation should also hold for non-Lorentzian (such as Gaussian) spectra [11].
Fig. 1
Fig. 1
(Color online) 1/K2 versus Tc. 1/K2 versus flow velocity index of Method I [Eq. (1b)] is depicted by the thinner solid curve and that of Method II [Eq. (2c)] by the thicker solid curve. Their corresponding asymptotes [Eqs. (3) and (4)] are shown (more ...)
The Newtonian iterative method is often utilized to calculate the blood-flow index Tc via Eq. (1b) or (2c). This calculation is computationally expensive for large time-series data sets. In contrast, the asymptotic approximation shows that the blood-flow index Tc is linearly related to 1/K2 [Eqs. (3) and (4)]. We therefore propose that this asymptotic approximation could be used to simplify the LSI analysis method when calculating blood-flow change.
To validate this asymptotic approximation and the simplified LSI analysis method, we performed flow phantom experiments at different physiologic flow rates and different exposure times. The flow phantom consisted of 1% intralipid solution flowing in the polyethylene tubing (PE-50, 0.58 mm inner diameter) with mean velocities of 0 to 5 mm/s, generated via a digital syringe pump. LSI was performed on a vibration-isolation table using a 780 nm laser diode with Δλ ~1 nm and a modified commercial imager (Imager 3001, Optical Imaging, Mountainside, N.J.) at 25 Hz and 7 μm resolution. Figure 2 shows the normalized equation M9 versus known flow velocity V of the flow phantoms obtained at different exposure times T. equation M10 versus V were linearly correlated (R2=0.97 to 0.99) for all exposure times and velocity ranges studied, validating the asymptotic approximation and the simplified analysis method. The slope was largest for T=10 ms, indicating the optimal exposure time.
Fig. 2
Fig. 2
(Color online) Normalized value of equation M11 versus mean velocity V at different camera exposure times T obtained from flow phantoms. Normalization was performed with respect to the speckle contrast at the lowest velocity index (K0) for each exposure time such (more ...)
To demonstrate a novel application, we applied this approach to image blood flow in the retinas of three rats in which the animals breathed air followed by oxygen. Rats were anesthetized under ~1% isoflurane, mechanically ventilated, and paralyzed with pancuronium bromide (3 mg/kg first dose followed by 1 mg/kg/h, intraperitoneal administration). End-tidal carbon dioxide, heart rate, arterial oxygen saturation, and rectal temperature were monitored and maintained within normal physiologic ranges. LSI experiments were performed on a vibration-isolation table. Time-series LSI data using T=10 ms were collected over 3 min of air breathing and 3 min oxygen breathing. Two to four trials were acquired for each animal. Figure 3 shows the steady-state blood flow index maps of an animal breathing air and oxygen. The blood-flow index obtained by Newtonian iterative analysis of Method I [Eq. (1b), Fig. 3(A)] was half of that by Newtonian iterative analysis of Method II [Eq. (2c), Fig. 3(B)], as indicated by differences in scale bars. Both analysis Methods I and II, however, yielded identical percent-change maps between oxygen and air breathing. Importantly, the blood-flow index maps and the percent-change maps calculated using the simplified analysis of Method II [Eq. (4), Fig. 3(C)] were indistinguishable from Newtonian iterative analysis of Method II [Fig. 3(B)], confirming the validity of the asymptotic approximation and the simplified analysis method.
Fig. 3
Fig. 3
(Color online) LSI blood-flow index and percent-change maps of a rat retina. LSI was performed while the animal was breathing air or oxygen. Percent-change maps are the differences between oxygen and air breathing. LSI maps were obtained using A, Newtonian (more ...)
Moreover, blood flow in the retina was found to markedly decrease under oxygen breathing relative to air. This is because oxygen potently vasoconstricts retinal vessels, consistent with previous reports using laser Doppler flowmetry [13,14]. These results demonstrated the feasibility of imaging physiologically evoked blood flow changes in the retina at very high spatiotemporal resolution in a totally noninvasive manner.
It is interersting to note that Tc ranges from 100 to 400 in typical in vivo LSI experiments [3-5,7,12] including our study. This value is much smaller than the theoretical prediction because speckle contrast is influenced by many factors, including scattering properties of the tissue [2], degree of polarization [11], illumination angle [10], and the ratio of pixel size to speckle size [6]. These factors are reflected in the beta term in [6]. This beta term and thus its effects on speckle contrast are expected to cancel out when calculating percentage change between the stimulation and basal conditions. Nonetheless, the blood-flow index Tc estimated using the asymptotic approximation [Eq. (4)] are within 99.5% of that using the full equation (2c), consistent with Fig. 3. Moreover, the percent changes are essentially identical.
In summary, we demonstrated that the commonly used approximation in the LSI equation is valid when calculating blood-flow changes and introduced a time-efficient LSI analysis method. Moreover, we demonstrated a novel application of LSI by imaging blood flow of the rat retinas. Because blood flow is a critical physiological parameter and its perturbation has been implicated in the early stages of many retinal diseases, LSI has the potential to provide valuable physiologic information and to serve as an early surrogate marker for many retinal diseases.
This work is supported by the Department of Veterans Affairs (VISN7 VA Career Development Award) and the NIH (NEI R01-EY014211d).
OCIS codes: 170.6480, 170.3880, 170.4470.
1. Fercher AF, Briers JD. Opt. Commun. 1981;37:326.
2. Dirnagl U, Kaplan B, Jacewicz M, Pulsinelli W. J. Cereb. Blood Flow Metab. 1989;9:589. [PubMed]
3. Dunn AK, Bolay H, Moskowitz MA, Boas DA. J. Cereb. Blood Flow Metab. 2001;21:195. [PubMed]
4. Bolay H, Reuter U, Dunn AK, Huang Z, Boas DA, Moskowitz MA. Nat. Med. 2002;8:136. [PubMed]
5. Cheng H, Luo Q, Liu Q, Lu Q, Gong H, Zeng S. Phys. Med. Biol. 2004;49:1347. [PubMed]
6. Bandyopadhyay R, Gittings AS, Suh SS, Dixon PK, Durian DJ. Rev. Sci. Instrum. 2005;76:093110.
7. Wang Z, Hughes S, Dayasundara S, Menon RS. Invest. Ophthalmol. Vis. Sci. 2007;27:258.
8. Zakharov P, Volker A, Buck A, Weber B, Scheffold F. Opt. Lett. 2006;31:3465. [PubMed]
9. Bonner R, Nossal R. Appl. Opt. 1981;20:2097. [PubMed]
10. Goodman JW. Proc. IEEE. 1965;53:1688.
11. Goodman JW. Statistical Optics. Wiley; 1985.
12. Cheng H, Luo Q, Wang Z, Gong H, Chen S, Liang W, Zeng S. Appl. Opt. 2003;42:5759. [PubMed]
13. Clermont AC, Brittis M, Shiba T, Mcgovern T, King GL, Bursell SE. Invest. Ophthalmol. Visual Sci. 1994;35:981. [PubMed]
14. Riva CE, Cranstoun SD, Mann RM, Barnes GE. Invest. Ophthalmol. Visual Sci. 1994;35:608. [PubMed]