In this paper, we develop a random process theory to explain the laser speckle phenomena. The relation between the probability distribution of speckle’s integrated intensity random process and the relative velocity is derived. Based on the random process theory, traditional spatial or temporal laser speckle contrast analysis (i.e. spatial or temporal LASCA) can be derived as the spatial or temporal estimators respectively. Both spatial LASCA and temporal LASCA suffer from noise due to insufficient statistics and nonstationarity in either spatial or temporal domain. Furthermore, either LASCA results in a reduction of spatial or temporal resolution. A new random process estimator is proposed and able to overcome these drawbacks. In an in-vitro study, random process estimator outperforms either spatial LASCA or temporal LASCA by providing much higher SNR (random process estimator vs. spatial LASCA vs. temporal LASCA: 33.64±6.87 () vs. 9.08±2.85 vs. 3.83±1.05). In an in-vivo structural imaging study, random process estimator efficiently suppresses the noise in contrast image and thus improves the distinguishability of small vessels. In a functional imaging study of cerebral blood flow change in the somatosensory cortex induced by rat’s hind paw stimulation, random process estimator provides much lower estimation errors in single trial data (random process estimator vs. temporal LASCA: 0.31±0.03 vs. 1.36±0.09) and finally leads to higher resolution spatiotemporal patterns of cerebral blood flow.
© 2009 OSA
Structural and functional brain imaging has become a very important tool in neuroscience research. As a simple two-dimensional imaging method, laser speckle imaging (LSI)  has been widely used in imaging the cerebral vasculature and cerebral blood flow (CBF) , e.g. the vasculature in brain tumor  and retina , the functional changes of CBF under focal cerebral ischemia , hypotension  and/or peripheral electrical stimulation . Both spatial  and temporal  laser speckle contrast analysis methods, i.e. spatial LASCA and temporal LASCA respectively, have been used to estimate the contrast values () which is inversely proportional to the relative velocity (v) . However, spatial LASCA decreases the spatial resolution while temporal LASCA is under the cost of temporal resolution. Furthermore, the inevitable noise in both LASCA methods not only decreases the distinguishability of vessels, but also results in errors in estimating functional changes of blood flow.
In this paper, we develop a random process theory for explaining the speckle phenomena and derive the relation between the velocity of CBF and the probability distribution of speckle’s random process. With the help of the random process theory, both LASCA methods are directly derived as the spatial and temporal estimators, and the noise is also modeled and analyzed. Finally, based on the property of noise, a new random process estimator is proposed to improve the signal-noise ratio (SNR) of the estimation, which thus leads to a better structural and functional imaging of CBF.
2.1 Background and preparation
When coherent laser light illuminates a surface of rat’s cortex, the random interference patterns, i.e. laser speckles, are produced by the coherent addition of scattered laser lights with slightly difference in light-path lengths. The motion of scattering particles, i.e. the red cells in blood vessels, changes the temporal pattern of speckles. The principle of LASCA is to estimate the relative velocity of CBF by analyzing the changing speckle patterns.
Essentially, light is a kind of electromagnetic radiation. According to the electromagnetic theory, the area under coherent laser illumination can be modeled as a ‘complex electric field’ [10,11]. Following the dynamic light scattering (DLS) theory , the scatters’ velocityis related to’s autocorrelation function:13] proposed the first model: ( is the correlation time, ); based on diffusing wave spectroscopy (DWS) model , Bandyopadhyay et al.  suggested a more rigorous model: (γ is a constant parameter, ). In practice, it is very difficult to measure and directly, but it is common to study the intensity, defined as . Generally, the recorded laser speckle image contains the intensity information of a composition of speckles with same speckle size (λis the wavelength of laser light, is the f number of the imaging system) . For each speckle, the scattering particles’ velocity is related toand thus related to the statistical property of .
2.2 Random process theory for laser speckle phenomenon
Due to the randomness property of speckle phenomena, the intensityof one speckle is considered as a single realization of the continuous intensity random process. In this section, we are going to develop the random process theory to explain the laser speckle phenomenon and derive the relation betweenand velocity.
At any time, is a single realization of the corresponding intensity random variable whose probability distribution is related to the velocity. In practice, a CCD camera is used to acquire the intensity information of the imaging plane. By carefully adjusting the f number of the imaging system, we can make the pixel size of the recorded image equal to the speckle size so that there is only one speckle in one pixel. Therefore, the ‘intensity’ value of any pixel is actually equal to the integrated intensity over the exposure durationT:
Suppose the exposure time T is designed to be very short (), then at any time, the velocity can be approximately considered as a constant. Thus, the distributions of random variables are the same, and each time-limited continuous random process can be considered as a wide-sense stationarity (WSS) random process. The properties of WSS random process  include: (a) its mean value is a constant when; (b) its temporal average value of the product of any two samples, i.e. and, , depends only upon the time interval .
Based on the property (a), when the number of ensemble samples, we have:
According to the Siegert relation , we have:
Suppose we obtain infinite temporal realizations of, then we have:Eq. (5) is applied in the third step; (iii) Eq. (3) is used in the fifth step; (iv) [15,17] is used in the last step.
By simply manipulating Eq. (6), we conclude:
Usingand, we rewrite Eq. (7) and define the contrast of integrated intensity random variable, i.e. , as follows:18].
Based on different models describing the relation between and , the relation between and can also be derived directly [13,15,19,20]. For any single speckle, if we can obtain its ensemble realizations (samples) of the integrated intensity random variable at time , the ensemble samples can be used to estimate the and then the relative velocity v can be estimated.
2.3 The spatial and temporal estimators and the estimation noise
In practice, at any time, we can only get one realization of , i.e. , for each speckle (pixel). So we have to estimate in other ways. An alternative method is to estimate the and using the intensities of the pixel () and its surrounding pixels () based on the assumption that velocities at the surrounding pixels are very close. When this assumption is true, the distributions of integrated intensity variable at different pixels are the same. Thus, spatial samples can be considered as ensemble samples approximately: and . The estimations of , and can be obtained as follows:
It is noted that this Eq. (9) is actually the spatial LASCA method . Such estimation works well in most in-vitro simulation studies when the moving scattering particles in a large area have the same velocity. However, this assumption does not hold for a complex biological tissue, such as brain, where there is: 1) an abundance of scattering particles such as red cells, and 2) a large number of vessels with different diameters and velocities.
Next, we analyze the estimation noise in spatial LASCA. Generally, the estimation of and in a spatial window centered at the interesting pixel could be affected by (i) limited size of the window, which results in the estimation noises and; and (ii) the spatial inhomogeneity of velocities within the window, which could result in estimation noises and . Therefore, and can be modeled as follows:
Although a larger window (e.g. or larger) can result in smallerandbased on central limit theory, the andwill increase due to higher inhomogeneity of velocities. In contrast, a smaller window (e.g. ) will result in smallerand, but largerand.
Another problem of spatial LASCA is the loss of spatial resolution. When we need a high spatial resolution image of cerebral vessels, spatial LASCA is not appropriate. To preserve the spatial resolution, the temporal estimator was developed.
If we assume that the velocity in the interesting pixel is constant during the time interval, then the probability distributions of integrated intensity variables are the same. In practice, we use a CCD camera to record the temporal samples () of integrated intensity random process:
Again the temporal estimator is actually the temporal LASCA method . The temporal LASCA is also contaminated by the estimation noise, i.e. the noise components and in and as follows:
Temporal LASCA preserves the spatial resolution because and are calculated based on of the same pixel. Therefore, high resolution structural imaging of cerebral blood vessel network is achieved. When imaging rat’s CBF under stable condition, the velocity in each speckle is approximately stable during the recording, and the noise component can be ignored. However, if the blood velocity is changing, e.g. elicited in response to functional stimulation, the noise will be prominent.
To achieve a high performance of statistical calculation (low noise), usually more than 50 frames are used in the temporal LASCA. In practice, the frame rate of the CCD camera is limited, for example in this study. So, the temporal resolution is compromised. Fewer frames, i.e. small time window, will reduce the loss of temporal resolution, but leads to larger noise.
2.4 Random process estimator
Based on the properties of integrated intensity random process, we propose a new estimator, called random process estimator hereafter, which provides high SNR estimation and high spatiotemporal resolution for both structural and functional imaging.
The velocity of CBF can be considered as a constant within a short exposure duration (e.g. ). For an interesting pixel, we calculate the using a window from each frame. Because the window is small, there would be large noise component and small in . Then, Eq. (11) is simplified into: , where.
It is noted that each is related to the velocity and thus could change during the imaging procedure. The noise component is due to the limited number of samples, so it is affected by the number of samples. Based on the central limit theorem, can be hypothesized to follow a zero-mean Gaussian distribution, i.e. , which will be confirmed based on the simulation data in Discussion. Therefore, the discrete sequences can be considered as independent and identically-distributed (IID) random variables with the same Gaussian distribution.
Then we calculate the average of as the random process estimator for averaged:
is used to estimate the averaged in the random process estimator, so there is no loss of spatial resolution in the calculation of. Actually, each can be considered as a ‘spatial estimation’ of with a spatial window, which leads to a large noise component without noise. Ignoring the noise in Eq. (10), we have: . Then we get:
Similarly, we model as a sequence of IID random variables with a zero-mean Gaussian distribution based on the central limit theorem, i.e. (this assumption will also be confirmed based on the simulation data in Discussion). Therefore, when , and then, where corresponds to the averaged velocity within.
Finally, the new random process estimator is used to estimate for the averaged velocity during time interval as follows:
Compared with spatial LASCA, provides much higher spatial resolution because the calculation of uses only a window. Furthermore, because the frame number N only affects the averaging (denoising) performance of the random process estimator, we can use fewer frames (e.g. ) than the temporal LASCA (images) to achieve a higher temporal resolution. Even using only 10 frames, the denoising performance of random process estimator is still significant while temporal LASCA is contaminated by both and. Therefore, random process estimator is very appropriate for functional imaging where both spatial and temporal resolutions need to be preserved as much as possible.
Although temporal LASCA can provide high resolution image of blood vessel network based on a large number of raw images, the noise component is unavoidable. Furthermore, some physiologic changes, e.g. body temperature, heart rate, blood pressure and medicine, could also lead to a time-varying CBF and thus the noise, which will decrease the distinguishability of small vessels. On the other hand, larger N in random process estimator will have better denoising performance according to Eq. (16) and Eq. (17). Therefore, although random process estimator produces a slight loss in spatial resolution compared to temporal LASCA, it significantly improves the distinguishability of vessels.
3.1 Imaging setup
The imaging plane was illuminated with a He-Ne laser beam source (1508P-O, Uniphase, CA) which was reshaped by optical lens to expand the range of illumination. A 12-bit cooled CCD camera (Sensicam SVGA, Cooke, MI) with a macro (1:1 maximum reproduction ratio) lens was held and focused on the imaging plane. Exposure time of the CCD camera was set to . Frame rate was . The resolution of image recorded by the CCD camera was pixels corresponding to an imaging area about .
3.2 In-vitro simulation experiment
A phantom is created to mimic blood vessels in vitro by flowing blood through plastic tubing using a syringe pump (PHD2000, Harvard Apparatus, MA). Real blood was obtained by ex-sanguination of rats destined for sacrifice and mixed with heparin sodium (10 IU/ml, Sigma Chemical Co., St. Louis, MO) to prevent clotting. A syringe filled with such blood was fitted into the syringe pump and connected to a polyethylene tube (internal diameter, PE90). Blood was infused into the tube at . When the flow was steady after , the tube was imaged using the imaging setup described above. 600 frames were acquired for analysis.
3.3 In-vivo structural and functional LSI experiment
The experimental protocol used in this study has been approved by the Animal Care and Use Committee of the Johns Hopkins Medical Institutions. The female Sprague-Dawley rats () were used in this study. To avoid the influence of surgical preparation on the rats’ CBF, the thinned skull preparation was done on the first day. On the second day, the rats were taken back for the structural and functional LSI experiment.
On the first day, the rat was anesthetized with intraperitoneal (IP) injection of mixture of of Ketamine and of Xylazine. The animal was placed in a stereotactic frame (David Kopf Instruments, Tujunga, CA) throughout the experiment. The scalp was shaved and disinfected with ethanol and povidone-iodine solution. All procedures were performed using standard sterile precautions. After a midline scalp incision, the galea and periosteum overlying the parietal bone bilaterally were swept and retracted laterally. A area centered at lateral and posterior to the bregma over the right cortex was thinned using a high speed dental drill (Fine Science Tools Inc. North Vancouver, Canada), until the inner cortical layer of bone was encountered. Rectal temperature was maintained at during the experiment using a homeothermic blanket system (model TP-500 T/Pump, Gaymar Industries, Inc., USA).
On the second day, the rat was first held in a transparent chamber with isofluorane gas and room air flow until becoming drowsiness. Its mouth and nose were then placed within a customized anesthesia mask in a well-fitting rodent size, which was connected to a C-Pram circuit designed to deliver and evacuate the gas through one tube. A mixed flow of isofluorane, oxygen and room air was delivered to the mask at the flow of for anesthesia. 600 frames were recorded for the structural imaging experiment.
Afterwards, subcutaneous needle electrode pairs (Safelead F-E3-48, Grass-Telefactor, West Warwick, RI) were inserted just under the skin of left hindpaw, one on the top of the hindpaw between 2 and 3 digits, the other one on the back skin of the hind paw, without direct contact to the nerve bundle. An isolated constant current stimulator (DS3, Digitimer Ltd., Hertfordshire, England) was used for the electrical stimulation of the hindpaw. Positive current pulses of magnitude and duration at a frequency of were used for hind paw stimulation. Each stimulation trial consisted of pre-stimulus 200 frames (), 200 stimulation frames (), and 400 post-stimulus frames (), which was repeated ten times with several-minute in-between breaks.
4.1 The denoising performance of random process estimator
The denoising performance of random process estimator is evaluated in the in-vitro simulation experiment. During the experiment, the blood flow velocity in the polyethylene tube was strictly maintained at using the syringe pump. According to Eq. (16) and Eq. (17), the noise parts and in the random process estimator would be very closed to zeros with. Therefore, for any pixel, and based on all 600 raw images were considered as the true and. Figure 1(a) shows one frame of the 600 raw images. Figure 1(c) is the estimated contrast image using random process estimator. One pixel (the white point in the white box area in Fig. 1(c)) is selected arbitrarily to compare the estimation results based on spatial LASCA, temporal LASCA and random process estimator.
For the selected pixel in Fig. 1(c), the result of spatial LASCA is calculated using a spatial window at each second () while the results of temporal LASCA and random process estimator are calculated using a temporal window (10 frames) starting from each second. Figure 1(b) shows the estimations of based on spatial LASCA, temporal LASCA and random process estimator respectively. As expected, all estimations vary around the true value (white dashed line). Among all three methods, temporal LASCA produces the largest fluctuations because of large noise component (small temporal window) in Eq. (15). Although both LASCA methods are contaminated by noise, spatial
LASCA still provides a better estimation using a comparably bigger window size (spatial LASCA: ; temporal LASCA: ). Furthermore, according to the hydrodynamics, the blood velocity is still not distributed uniformly within the tube. Therefore, besides the noise, there is also noise in the result of both LASCA methods. Compared to the conventional LASCA methods, random process estimator provides the best estimation of in the sense of smallest fluctuations. Since the estimations of are the same in temporal LASCA and random process estimator, we did not show the results of.
Figure 1(d) shows the estimations of at each second () based on spatial LASCA, temporal LASCA and random process estimator where random process estimator outperformed the others again. To quantitatively compare the denoising performances of different methods, we define the signal-noise ratio (SNR) for the estimation of as follows:Fig. 1(c), the averaged SNR of random process estimator is () while the spatial LASCA and temporal LASCA are only with the averaged SNR and respectively. Clearly, random process estimator provides higher denoising performance than either spatial LASCA or temporal LASCA.
Although spatial LASCA provides a slightly better estimation than temporal LASCA, the loss of spatial resolution is unavoidable and usually unacceptable. Since small vessels play an important role in the in-vivo structural and functional studies, we’ll only compare temporal LASCA and random process estimator in the rest sections.
4.2 Random process estimator for structural imaging
In the in-vivo experiment, 600 raw images () are recorded for imaging the structure of blood vessel network. During the imaging procedure, the rats were anesthetized and kept in stable condition, so blood velocities in the imaging area are assumed to be stable. Figure 2 shows the contrast images estimated from the first 2 (Fig. 2(a,d)), 10 (Fig. 2(b,e)) and 80 (Fig. 2(c,f)) frames using random process estimator and temporal LASCA respectively. For either method, the distinguishability of small vessels is improved when more frames are used. For example, the small vessel designated with black arrow in Fig. 2 (b) and (c) is hard to see in Fig. 2(a).
According to Eq. (15), when more frames are used, the noise component decreases but still exists in the result of temporal LASCA. Furthermore, if there are small changes of blood velocity during the imaging procedure, the noise is also introduced in. However, using random process estimator, the noise part is rapidly averaged out according to Eq. (16), 17) provided that more frames are used. More precisely, based on the same number of frames, random process estimator always provides better estimation and thus better distinguishability than temporal LASCA. For example, the small vessels in the white circled area in Fig. 2(c) are more distinguishable than in Fig. 2(f).
The contrast values along the vertical white line in Fig. 2(c) estimated from different number of frames using the two methods is plotted in Fig. 3 . There are six vessels denoted as V1~V6 crossing the line. In Fig. 3, the results of random process estimator are always smoother than temporal LASCA. In particular, based on only 2 frames, the noise in the result of temporal LASCA was prominent (Fig. 3(a)), while the result of random process estimator still reveals the denoising aspect especially in the vessel areas V1 and V2. For 10 frames (Fig. 3(b)), the random process estimator already provides an efficient estimation of with high
SNR while the result of temporal LASCA is still very noisy. Furthermore, the denoising performance of random process estimator is better in the tissue areas compared to vessel areas.
Although more frames can be used in the calculation by random process estimator, no significant improvement is found compared to Fig. 2(c). Therefore, 80 images are enough for the structural imaging using random process estimator. In conclusion random process estimator provides better distinguishability of structural information of cerebral blood vessels.
4.3 Random process estimator for functional imaging
In the functional stimulation experiment, each trial includes recording 800 raw images () continuously, i.e. 200 pre-stimulus frames (), 200 stimulation frames () and 400 post-stimulus frames (). Images within each second (10 frames) are used to calculate one contrast image using random process estimator and temporal LASCA. Therefore, for each trial, we obtain 80 contrast images ().is used to normalize the change of contrast values as follows ():
Because the stimulation procedure is 20s, other cortex areas may also be activated besides the somatosensory cortex. In this study, the area covering the activated somatosensory cortex is selected as the region of interest in the following discussions.
To analyze the denoising performance for functional imaging, the of the interesting region in pre-stimulation stage of the first trial are shown in Fig. 4 using random process estimator (the first row) and temporal LASCA (the second row) respectively. Since the rat’s condition was kept stable before the stimulation, there should not be significant changes in the contrast images and thus the values in images are expected to be much closed to zero. However, as shown in Fig. 4, temporal LASCA leads more noises than random process estimator. To quantify this estimation error, we define the ‘averaged errors’ in each as follows:
Figure 5 shows the averaged errors of single trial estimation using temporal LASCA and random process estimator respectively. For random process estimator, the averaged errors in all is only (), in contrast with for temporal LASCA.
correspond to functional changes during the stimulation and post-stimulation procedure. By averaging the of all ten trials, the averaged spatiotemporal changing map () is obtained. Figure 6 and Fig. 7 show the averaged changing map, using random process estimator and temporal LASCA respectively. As a reference, the first frames in Fig. 6 and Fig. 7 are the structural image calculated by random process estimator.
In the functional imaging experiment, the blood flow in the somatosensory cortex changes in response to electrical stimulation. Therefore, besides the noise, the noise is also introduced into the result of temporal LASCA. Hence, the result of temporal LASCA (Fig. 7) is too noisy to precisely quantify the changes of CBF in small vessels. However, random process estimator (Fig. 6) provides high SNR and high spatiotemporal resolution for the functional changing pattern of CBF.
In Fig. 6, the response of CBF induced by functional stimulation starts from and reaches maximum at, then recovers to the baseline at. All significant changes can be located to small vessels (e.g. V1~V4 in Fig. 6). V1 is the vessel with the most significant changes, which is activated earlier than other vessels and shows longer recovery duration. With the help of high spatial resolution, we also find that not all vessels are activated in this area (e.g. V5).
One pixel in vessel V1 is chosen for quantitative analysis (the white circle in the first frame of Fig. 6). The averaged value and standard deviation of at this pixel across all 10 trials are plotted in Fig. 8 for both methods. The result of temporal LASCA exhibits a noisy fluctuations () even in the pre-stimulation stage, so that only changes greater than can be considered as significant change during stimulation. However, the result of random process estimator shows much smaller fluctuations and a more stable profile. Furthermore, the standard deviation of the over all 10 trials using random process estimator are always less than that using temporal LASCA (Fig. 8 (b)), implying that random process estimator also suppress the noisy variations over trials.
5.1 The Gaussian noise assumptions for random process estimator
In the development of random process estimator, the most important hypotheses are the Gaussian noise models in Eq. (16) and Eq. (17). Using the data acquired in the in-vitro experiment, we could test the hypotheses.
For the selected pixel in Fig. 1(c), 600 estimations of can be obtained using a spatial window. The noise component of each estimation is then calculated by subtracting the true value. The probability density function (pdf) is then calculated from all 600 by the kernel smoothing density estimation (Fig. 9(a) ). The pdf is very close to a Gaussian distribution with (). The P-P plot of all 600 in Fig. 9(b) is also very close to that of a Gaussian distribution. The hypothesis that the distribution of is a Gaussian distribution is further confirmed by the K-S test with a high z value ().
Similarly, the distribution of the noise component, i.e. , is also tested by K-S test with a high. The pdf of the distribution () and the P-P plot are shown in Fig. 9(c) and (d) respectively. These tests confirm that the distribution of is also a Gaussian distribution.
5.2 Comparison with hybrid temporal and spatial method
Recently, a hybrid temporal and spatial method was proposed to improve the SNR by spatially averaging (spatial window) the result of temporal LASCA which was successfully applied in imaging rat’s retinal blood flow . Similar to spatial and temporal LASCA, the hybrid method works well when blood flow is constant within both spatial and temporal window, i.e. noise does not exist. We further process the same data set in our study using the algorithm in  to compare the performances of random process estimator and hybrid method in SNR, distinguishability of small vessels and revealing the functional responses in small vessels.
In the in-vitro simulation experiment, for the selected pixel in Fig. 1(c), the hybrid method produces a higher SNR (16.48) than temporal LASCA (2.75) and spatial LASCA (10.34), which, however, is still significantly lower than the SNR by random process estimator (31.09), because the use of a spatial window increases thenoise while suppressing thenoise. In the structural imaging of the selected area in Fig. 10(a) , the details of small vessels (indicated with black arrows) would be hardly visible by the hybrid method with aspatial window (Fig. 10(b)). Although using a spatial window in hybrid method slightly improves the spatial resolution (Fig. 10(c)), but the distinguishability of the vessels is still not as good as that by random process estimator (Fig. 10(d)) due to the decrease of SNR. Furthermore, the hybrid method cannot suppress thenoise when there are significant blood flow changes within the temporal window. For example, noise is so large that the small functional changes are strongly contaminated by background noise in temporal LASCA (Fig. 11(b) ) under the functional stimulation. Hybrid method will filter out noise and the weak functional signals as well (Fig. 11(c)). However, random process estimator well reveals the functional changes in the small vessel while removing the noise (Fig. 11(d)).
In this paper, we develop the random process theory for the laser speckle phenomenon and thus propose a new random process estimator to estimate the true value of . Compared with traditional spatial LASCA and temporal LASCA, random process estimator provides higher SNR and higher spatiotemporal resolution contrast images for both structural and functional imaging of cerebral blood vessels and flow.
This work is supported by NIH/NIA1R01AG029681. S. Tong is also supported by the New Century Talent Program by the Ministry of Education of China, and Shanghai Shuguang Program (07SG13). P. Miao is also supported by the China Scholarship Council.
References and links
1. J. Briers, “Laser Doppler, speckle and related techniques for blood perfusion mapping and imaging,” Physiol. Meas. 22(4), 35–66 (2001). [CrossRef]
3. D. Zhu, W. Lu, Y. Weng, H. Cui, and Q. Luo, “Monitoring thermal-induced changes in tumor blood flow and microvessels with laser speckle contrast imaging,” Appl. Opt. 46(10), 1911–1917 (2007). [CrossRef] [PubMed]
4. H. Cheng, Y. Yan, and T. Duong, “Temporal statistical analysis of laser speckle images and its application to retinal blood-flow imaging,” Opt. Express 16(14), 214–219 (2008). [CrossRef]
5. A. Kharlamov, B. R. Brown, K. A. Easley, and S. C. Jones, “Heterogeneous response of cerebral blood flow to hypotension demonstrated by laser speckle imaging flowmetry in rats,” Neurosci. Lett. 368(2), 151–156 (2004). [CrossRef] [PubMed]
6. T. Durduran, M. G. Burnett, G. Yu, C. Zhou, D. Furuya, A. G. Yodh, J. A. Detre, and J. H. Greenberg, “Spatiotemporal Quantification of Cerebral Blood Flow During Functional Activation in Rat Somatosensory Cortex Using Laser-Speckle Flowmetry,” J. Cereb. Blood Flow Metab. 24(5), 518–525 (2004). [CrossRef] [PubMed]
7. J. Briers and S. Webster, “Laser speckle contrast analysis (LASCA): a nonscanning, full-field technique for monitoring capillary blood flow,” J. Biomed. Opt. 1(2), 174–179 (1996). [CrossRef]
9. P. Miao, M. Li, H. Fontenelle, A. Bezerianos, Y. Qiu, and S. Tong, “Imaging the Cerebral Blood Flow with Enhanced Laser Speckle Contrast Analysis (eLASCA) by Monotonic Point Transformation,” IEEE Trans. Biomed. Eng. 56(4), 1127–1133 (2009). [CrossRef] [PubMed]
10. P. Lemieux and D. Durian, “Investigating non-Gaussian scattering processes by using nth-order intensity correlation functions,” J. Opt. Soc. Am. A 16(7), 1651–1664 (1999). [CrossRef]
11. D. Boas and A. Yodh, “Spatially varying dynamical properties of turbid media probed with diffusing temporal light correlation,” J. Opt. Soc. Am. A 14(1), 192–215 (1997). [CrossRef]
12. B. Berne, and R. Pecora, Dynamic Light Scattering: with Applications to Chemistry, Biology and Physics (Dover Publications, 2000).
13. A. Fercher and J. Briers, “Flow visualization by means of single-exposure speckle photography,” Opt. Commun. 37(5), 326–330 (1981). [CrossRef]
15. R. Bandyopadhyay, A. Gittings, S. Suh, P. Dixon, and D. Durian, “Speckle-visibility spectroscopy: A tool to study time-varying dynamics,” Rev. Sci. Instrum. 76(9), 093110 (2005). [CrossRef]
16. G. Grimmett, and D. Stirzaker, Probability and random processes (Oxford University Press, 2001).
18. J. W. Goodman, Statistical Optics (Wiley 1985).
21. H. Cheng, Y. Yan, and T. Q. Duong, “Laser speckle imaging of rat retinal blood flow with hybrid temporal and spatial analysis method,” Proc. SPIE 7163, 716304 (2009). [CrossRef]