Laser Speckle Contrast Imaging for Measuring Cerebral Blood Flow Changes Caused by Electrical Sensory Stimulation
Laser Speckle Contrast Imaging for Measuring Cerebral Blood Flow Changes Caused by Electrical Sensory Stimulation
Journal of the Optical Society of Korea. 2016. Feb, 20(1): 88-93
Copyright © 2016, Optical Society of Korea
  • Received : December 10, 2015
  • Accepted : January 27, 2016
  • Published : February 25, 2016
Export by style
Cited by
About the Authors
Ahra, Cho
School of Mechatronics, Gwangju Institute of Science and Technology, Gwangju 61005, Korea
Chanmi, Yeon
Department of Medical System Engineering, Gwangju Institute of Science and Technology, Gwangju 61005, Korea
Donghyeon, Kim
School of Information and Communications, Gwangju Institute of Science and Technology, Gwangju 61005, Korea
Euiheon, Chung
Department of Medical System Engineering, Gwangju Institute of Science and Technology, Gwangju 61005, Korea

Recently laser speckle contrast (LSC) imaging has become a widely used optical method for in vivo assessment of blood flow in the animal brain. LSC imaging is useful for monitoring brain hemodynamics with relatively high spatio-temporal resolution. A speckle contrast imaging system has been implemented with electrical sensory stimulation apparatus. LSC imaging is combined with optical intrinsic signal imaging in order to measure changes in cerebral blood flow as well as neural activity in response to electrical sensory stimulation applied to the hindlimb region of the mouse brain. We found that blood flow and oxygen consumption are correlated and both sides of hindlimb activation regions are symmetrically located. This apparatus could be used to monitor spatial or temporal responses of cerebral blood flow in animal disease models such as ischemic stroke or cortical spreading depression.
Optical approaches can bring structural and functional information together in preclinical brain research. High-resolution and real-time imaging of cerebral neocortex has been accelerated with optical tools to monitor vascular function and microvascular structure of the brain in normal and pathophysiological conditions [1] . Currently several research groups have demonstrated functional brain imaging of the cerebral cortical network with high spatio-temporal resolution [2 - 4] .
There are several optical methods for vascular imaging of the brain including white light imaging, fluorescence imaging, and laser speckle contrast imaging. White light imaging displays intrinsic reddish vasculature as hemoglobin molecules within red blood cells absorb a specific spectrum of illumination light. Fluorescence imaging provides high vascular contrast by blocking excitation light and passing only the emission light from functional vasculature using fluorescent vascular tracer that is intravenously injected into the circulation prior to imaging [5] . However, both white light and fluorescence imaging have limitations: the former being of limited use in discerning functional and non-functional blood vessels and the latter requiring exogenous fluorescent tracers.
To overcome these limitations, laser speckle contrast imaging can be used in the field of biomedical imaging. Unlike other techniques, laser speckle contrast imaging provides label free high contrast and real time functional imaging without the need of exogenous contrast agent, although it is sensitive to vibrations and requires stable fixation of the target region to minimize any perturbation caused by breathing and heartbeat. Therefore, laser speckle contrast imaging could be used for superficial vascular imaging on skin, retina or cerebral cortex [6] .
Laser speckle is a random granular interference pattern due to the roughness - on the scale of an optical wavelength - of surface or to multiple scattering events within tissue through which the photons travel. When coherent light illuminates a sample area, backscattered light fields are superposed together resulting in constructive or destructive interference. The resultant fluctuations of intensity, or laser speckles, are imaged onto a camera as grainy patterns [7] . When the scattered particles are in motion in some instance such as particle suspensions or blood cells in circulation, the speckle pattern can change during the integrated camera exposure time. If the particles are moving fast, the speckle image at the corresponding region appears blurred because the changing speckle patterns are averaged out during the camera integration time. Hence, the intensity fluctuations of speckles in terms of normalized standard deviation either in space or in time can provide dynamic information of scatterers within the illuminated area. Moreover, qualitative blood flow can be measured by analyzing the blurring of speckles or speckle contrast [8 - 11] . The speckle contrast can be calculated in a cortical area to provide 2-dimensional cerebral blood flow map with relatively high spatiotemporal resolution [1] . Multiple modality imaging (i.e. intrinsic signal) of functional responses could provide a better picture of the hemodynamic responses to sensory stimulation both in spatial and temporal domains.
In this paper, we delineated optimal stimulation parameters to maximizing laser speckle contrast with a home-built laser speckle contrast (LSC) imaging setup, and applied it for measuring cerebral blood flow changes in response to electrical sensory stimulation with an experimental mouse model. The comparison with optical intrinsic signal (OIS) imaging demonstrates an exemplary correlation of blood flow and neural activity-dependent changes.
Animal experiments were performed according to the institutional guidelines of the Gwangju Institute of Science and Technology, which are compliant with the ARRIVE guidelines ( ). C57BL/6 mice with cranial windows were used for experiments. Mice with 25-30 g body weight and ages of 8-10 weeks were anesthetized with Zoletil/Xylazine mixture in saline solution (60/10 mg/kg body weight). Body temperature was kept at 37℃ throughout the experiments. The cranial window with a customized angel ring was stably mounted on a stereotaxic animal holding and heating apparatus (Customized ring type in vivo heating system, Live Cell Instrument, Seoul, Korea).
The laser speckle imaging instrumentation consists of a helium-neon laser ( λ = 632.8 nm, 17 mW, Research Electro-Optics, Inc., Colorado USA), a scientific CMOS camera (Neo 5.5 sCMOS, Andor Technology Ltd., Belfast, UK), and a microscope body (BX-51, Olympus Co.) having an objective lens (RMS4X, NA 0.10, Olympus) as shown in Fig. 1 . The laser beam was coupled into a single mode fiber (core diameter = 4.3 μm, Thorlabs) to illuminate onto mouse cortex at an angle of approximately 30 - 45 degree from vertical direction. To manipulate the angle and position of the fiber tip, a micromanipulator (MN-153, NARISHEGE, Tokyo, Japan) was used.
PPT Slide
Lager Image
Schematic of laser speckle imaging instrumentation. A He:Ne laser (λ = 632.8 nm) is coupled into a single mode fiber to illuminate the cortex area using a micro-manipulator. The illuminated area is imaged onto a sCMOS camera through the microscope. (L1: objective lens, L2: tube lens, f=180 mm)
To analyze the speckle data, speckle contrast K is calculated as,
PPT Slide
Lager Image
where σ and < I > represent the standard deviation and mean value, respectively. The value of speckle contrast lies between 0 and 1: The value of 1 means static condition while the value close to 0 signifies fast motion of scatterers such as moving blood cells [12] . For the spatial domain analysis, the speckle contrast was determined with 7 by 7 window pixels which contain enough pixels for sensitive detection of perfused blood vessels [13 , 14] .
The schematic of electrical sensory stimulation apparatus synchronized with camera acquisition is illustrated in Fig. 2 . To monitor functional activation of somatosensory cortex, a couple of needle electrodes made of stainless steel were inserted just underneath the skin of the hindlimb [6] . During electrical stimulation experiment, the animals were under gas anesthesia (isofluorane at 1.2-1.5%). A current output module (NI 9265, 20 mA, National Instruments, USA) was used to supply electrical pulses (0.5 mA, 1 msec duration at 3 Hz) of current. The current output was controlled with LabView program and was synchronized with camera acquisition using a bidirectional digital I/O module (NI 9403, National Instruments, USA). Twenty repeated trials were averaged in order to increase the signal to noise ratio (SNR) [15] . Each experimental trial lasted for 60 sec including 3 sec of baseline and 5 sec of stimulation ( Fig. 2 ) during which the image acquisition took 20 sec at 5 Hz with 5 msec exposure time.
PPT Slide
Lager Image
Schematic of electrical sensory stimulation apparatus synchronized with camera acquisition.
To represent the cerebral blood flow, the relative changes of speckle contrast can be computed as
PPT Slide
Lager Image
where Kt and
PPT Slide
Lager Image
mean the speckle contrast value at time t and averaged contrast value during the baseline measurement, respectively [15] . Although the speckle contrast, the spatial blurring in a small region of the image, is not directly proportional to the flow velocity, the change of speckle contrast could be considered as the changes in cerebral blood flow [16 , 17] . On the other hand, the exact relationship between the speckle contrast and the underlying blood flow is nonlinear and is currently an active area of research as speckle contrast-based mapping only provides qualitative estimates of blood flow. There has been better understanding to link the speckle contrast value to the characteristic correlation time of the speckles via improved quantitative modeling as the correlation time is known to be inversely proportional to the blood flow [10 , 17 , 18] . A threshold below 3% was regarded as random noises to be eliminated in the image of changes in contrast as the noise levels of ex-vivo and in-vivo samples are 2-3% in our experimental system.
In order to monitor the response of tissue perfusion that is related to the blood flow as well as activity-dependent oxygenation due to the functional activation, LSC imaging was combined with optical intrinsic signal (OIS) imaging as this reflects local blood volume and oxygen saturation signals through the back-reflection of light when the surface of the brain is illuminated with red LED (625 nm, M625L3, Thorlabs Inc.) [19] . Therefore, the OIS imaging can guide the activation region by analyzing oxygen consumption of hemoglobin [20] .
For OIS imaging, twenty trials are averaged to generate 200 sequential frames at 10 Hz for 20 sec. The frames for the initial 3 seconds of the sequence, taken before stimulation, are averaged together to give a single baseline image. The frames for 5 seconds taken during stimulation, are averaged together to give a single stimulation image. The ratio can be calculated to generate a functional map. Each frame in the temporal sequence was normalized to the mean baseline image, yielding percent changes in light reflectance [21 , 22] .
PPT Slide
Lager Image
OIS is an appropriate tool to obtain information about the functional response to somatosensory stimulation although a single wavelength illumination allows only qualitative mapping and does not give separate information on oxygenation and blood volume. However, in our study, single wavelength OIS illumination was sufficient to provide a localized region of interest caused by changes in deoxyhemoglobin. The red wavelength we used for OIS imaging provides maps that are highly spatially correlated with neuronal activity as the absorption of oxyhemoglobin is much lower than that of deoxyhemoglobin in 600~630 nm range to produce high signal-to-noise ratio compared with near infrared wavelength [21] . By analyzing reflectance change between baseline and stimulation period using OIS imaging, the left hindlimb activation region was exposed to electrical sensory stimulation for 5 seconds. The spatial changes in reflected light intensity from LED illumination due to hindlimb stimulation is shown in Fig. 3 . Based on the reflectance change, a rectangular region of interest (ROI, size of 0.5 mm × 0.5 mm) was chosen at the right side of the bregma.
PPT Slide
Lager Image
Functional activation of left hindlimb monitored by optical intrinsic signal imaging. The blue circle indicates the location of the bregma and color scale bar shows the OIS ratio (%).
Because speckle contrast is related to the velocity of scatterers in blood flow, the changes of speckle contrast are considered as the changes in cerebral blood flow [23] . The images of the changes in speckle contrast including the 5 sec of hindlimb electrical sensory stimulation are shown in Fig. 4 . The series of images are shown at 0.4 sec intervals. The changes in speckle contrast are displayed for all pixels and the localized region identified in the OIS imaging ( Fig. 3 .) was also activated in LSC imaging which is indicated by the yellow box in the Fig. 4 . Especially, a peak of relative increase in speckle contrast due to cerebral blood flow was over 5%. The time course of the contrast changes averaged in the ROI that is centered over the activated area is shown in Fig. 5 . The peak relative changes in speckle contrast (~3%) during the 5 sec stimulation occurs at 6.7 sec after the initiation of stimulation. After reaching the peak point, changes in cerebral blood flow or speckle contrast were steadily decreased and returned to the baseline level within 6 sec.
PPT Slide
Lager Image
Changes in speckle contrast (cerebral blood flow) in response to the electrical sensory stimulation of the hindlimb from 0 to 5 sec based on green bar laser speckle contrast imaging. The speckle contrast images over 3% threshold are overlaid with a white light image at 0.4 sec interval.
PPT Slide
Lager Image
Time courses of percent change in speckle contrast averaged in left hindlimb sensory region (ROI) after low pass filtering (LPF). The stimulation intervals are shown as a green bar.
The use of 633 nm HeNe laser and not conventional near-infrared light sources for our LSCI was to better match with the intrinsic signal information. Not only does the longer coherence length have advantage for laser speckle, this wavelength can further emphasize the changes in optical properties during our hemodynamic measurement. The molar extinction coefficient of deoxyhemoglobin is about an order of magnitude larger than that of oxyhemoglobin and thus HeNe illumination provides sensitive detection of hemoglobin changes [9] .
Time courses of percent changes in speckle contrast and intrinsic signal averaged in the left sensory region is illustrated in Fig. 6 . The starting times of speckle contrast and intrinsic signals right after electrical stimulation are comparable. Subsequent to sensory stimulation, the cerebral blood flow (speckle contrast) starts to increase at approximately 1.4 sec, while the highest oxygen consumption or the early reflectance decay known as initial dip, occurred at a similar time. This phenomenon can be explained as correlation between oxygen consumption and blood flow. The initial dip is caused by deoxyhemoglobin increase and oxyhemoglobin decrease due to oxygen consumption in response to neural activity [6 , 24 , 25] . Upon functional activation, neurons require more oxygen than their resting state and subsequently the blood flow is increased via neurovascular coupling [2] .
PPT Slide
Lager Image
Time courses of percent change in speckle contrast (Red line) and reflectance (Blue line) averaged in left hindlimb (LH) sensory region after low pass filtering. The stimulation intervals are shown as a green bar.
An equivalent experimental set was performed on the right hindlimb of the same animal to find sensory region of the opposite side using both OIS imaging reflectance and LSC imaging contrast as shown Fig. 7 . We found that both sides of the hindlimb region were symmetrically located in 1.2-1.5 mm apart from the center of the bregma, in a manner consistent with a previous craniotomy study [26] .
PPT Slide
Lager Image
Changes in contrast and OIS ratio images of left and right hindlimb sensory region. Both sides of hindlimb region are symmetrically located in 1.2-1.5 mm apart from the center of the bregma (at 5.2 sec.)
To eliminate the random noise of laser speckle due to the vibration or breathing, the threshold technique (3%) was used [27] . However, the weak signal of changes in blood flow at 0.3 mA was not easily noticeable as it was not distinguishable from random noise as in Fig. 8 . It is important to find optimal amplitude of electrical current in order to minimize any detrimental and unexpected effect on the experimental animal such as motor reaction by stimulus or edema from invasive electrode. Each experimental block includes 40 repeated trials. Each trial consists of baseline (3 sec) and stimulation time (5 sec). To use electrical current of 0.3 mA, it was necessary to minimize the noise level. As the major noise sources were environmental vibrations, mainly generated by the cooling fan of the camera or neighboring equipment, further suppression of environmental or mechanical noise was required.
PPT Slide
Lager Image
Optimization of stimulation parameter using change in contrast by laser speckle contrast imaging with respect to electrical current (a) control, (b) 0.3 mA, (c) 0.5 mA, (d) 0.7 mA.
As spatial speckle analysis was used in our laser speckle contrast imaging, spatial resolution was sacrificed. To improve spatial resolution, a sliding window technique can be used [3] . However this requires N 2 times longer time than simple spatial analysis (N=7 in our case) and the trade-off between the computational burden and finer spatial analysis should be taken into account [18] .
LSC imaging provides relatively high spatio-temporal resolution of cerebral blood flow using simple optical instrumentation by analyzing the laser speckle pattern caused by the motion of scatterers that are mostly blood cells. This study demonstrates that the combined imaging of laser speckle contrast and intrinsic optical signal can provide correlative information with electrical sensory stimulation of a limb towards preclinical stroke research. To induce hemodynamic response by functional activation, an electrical sensory stimulation apparatus was built and synchronized with camera acquisition. The LSC imaging was combined with optical intrinsic signal (OIS) imaging in order to monitor response of blood flow and oxygenation by functional activation. The LSC and OIS imaging was alternatively performed to obtain combined data. We found that changes in speckle contrast and optical intrinsic signal have notable correlation: the temporal trajectories were similar while the latter shows initial dip at the beginning of stimulation. In addition, the hindlimb region in the cortex was symmetrically located in approximately 1.2-1.5 mm from the center of the bregma in mouse cortex.
As a future direction, a fully combined multi-modal simultaneous imaging system of both LSC and OIS would provide simultaneous hemodynamic information on the changes in blood flow and oxygenation in response to electrical sensory stimulation.
Color versions of one or more of the figures in this paper are available online.
The authors thank Jihee Ryu for technical help and A. Mohaimen Safi of GIST for his careful reading of the manuscript. The work is supported by the Bio & Medical Technology Development Program and Basic Science Research Program through the National Research Foundation (NRF) funded by the Ministry of Science, ICT & Future Planning (No. 2011-0019633), and by a grant from the Institute of Medical System Engineering (iMSE) at GIST, the “GIST-Caltech Research Collaboration” Project through a grant provided by GIST.
Dunn A. K. , Bolay H. , Moskowitz M. A. , Boas D. A. (2001) “Dynamic imaging of cerebral blood flow using laser speckle,” Journal of Cerebral Blood Flow & Metabolism 21 195 - 201
Durduran T. , Burnett M. G. , Yu G. , Zhou C. , Furuya D. , Yodh A. G. , Detre J. A. , Greenberg J. H. (2004) “Spatiotemporal quantification of cerebral blood flow during functional activation in rat somatosensory cortex using laser-speckle flowmetry,” Journal of Cerebral Blood Flow & Metabolism 24 518 - 525
Li N. , Jia X. , Murari K. , Parlapalli R. , Rege A. , Thakor N. V. (2009) “High spatiotemporal resolution imaging of the neurovascular response to electrical stimulation of rat peripheral trigeminal nerve as revealed by in vivo temporal laser speckle contrast,” Journal of Neuroscience Methods 176 230 - 236    DOI : 10.1016/j.jneumeth.2008.07.013
Weber B. , Burger C. , Wyss M. T. , von Schulthess G. K. , Scheffold F. , Buck A. (2004) “Optical imaging of the spatiotemporal dynamics of cerebral blood flow and oxidative metabolism in the rat barrel cortex,” Eur. J. Neurosci. 20 2664 - 2670    DOI : 10.1111/j.1460-9568.2004.03735.x
Weber B. , Helmchen F. 2014 Optical Imaging of Neocortical Dynamics Springer
Ances B. M. , Buerk D. G. , Greenberg J. H. , Detre J. A. (2001) “Temporal dynamics of the partial pressure of brain tissue oxygen during functional forepaw stimulation in rats,” Neurosci. Lett. 306 106 - 110    DOI : 10.1016/S0304-3940(01)01868-7
Goodman J. W. 2007 Speckle Phenomena in Optics: Theory and Applications Roberts and Company Publishers
Kazmi S. M. S. , Parthasarthy A. B. , Song N. E. , Jones T. A. , Dunn A. K. (2013) “Chronic imaging of cortical blood flow using multi-exposure speckle imaging,” Journal of Cerebral Blood Flow & Metabolism 33 798 - 808    DOI : 10.1038/jcbfm.2013.57
Mazhar A. , Cuccia D. J. , Rice T. B. , Carp S. A. , Durkin A. J. , Boas D. A. , Choi B. , Tromberg B. J. (2011) “Laser speckle imaging in the spatial frequency domain,” Biomed. Opt. Express 2 1553 - 1563    DOI : 10.1364/BOE.2.001553
Kazmi S. M. , Richards L. M. , Schrandt C. J. , Davis M. A. , Dunn A. K. (2015) “Expanding applications, accuracy, and interpretation of laser speckle contrast imaging of cerebral blood flow,” J. Cereb. Blood Flow Metab. 35 1076 - 1084    DOI : 10.1038/jcbfm.2015.84
Nadort A. , Woolthuis R. G. , van Leeuwen T. G. , Faber D. J. (2013) “Quantitative laser speckle flowmetry of the in vivo microcirculation using sidestream dark field microscopy,” Biomed. Opt. Express 4 2347 - 2361    DOI : 10.1364/BOE.4.002347
Briers J. D. , Webster S. (1996) “Laser speckle contrast analysis(LASCA): a nonscanning, full-field technique for monitoring capillary blood flow,” J. Biomed. Opt. 1 174 - 179    DOI : 10.1117/12.231359
Briers J. D. , He X.-W. (1998) “Laser speckle contrast analysis(LASCA) for blood flow visualization: improved image processing,” Proc. BiOS'98 International Biomedical Optics Symposium 26 - 33
Choi B. , Kang N. M. , Nelson J. S. (2004) “Laser speckle imaging for monitoring blood flow dynamics in the in vivo rodent dorsal skin fold model,” Microvasc. Res. 68 143 - 146    DOI : 10.1016/j.mvr.2004.04.003
Dunn A. K. , Devor A. , Dale A. M. , Boas D. A. (2005) “Spatial extent of oxygen metabolism and hemodynamic changes during functional activation of the rat somatosensory cortex,” Neuroimage 27 279 - 290    DOI : 10.1016/j.neuroimage.2005.04.024
Dunn A. K. , Devor A. , Bolay H. , Andermann M. L. , Moskowitz M. A. , Dale A. M. , Boas D. A. (2003) “Simultaneous imaging of total cerebral hemoglobin concentration, oxygenation, and blood flow during functional activation,” Opt. Lett. 28 28 - 30    DOI : 10.1364/OL.28.000028
Briers J. D. (2007) “Laser speckle contrast imaging for measuring blood flow,” Optica Applicata 37 139 -
Li H. , Liu Q. , Lu H. , Li Y. , Zhang H. F. , Tong S. (2014) “Directly measuring absolute flow speed by frequency-domain laser speckle imaging,” Opt. Express 22 21079 - 21087    DOI : 10.1364/OE.22.021079
Frostig R. D. , Lieke E. E. , Ts’o D. Y. , Grinvald A. (1990) “Cortical functional architecture and local coupling between neuronal activity and the microcirculation revealed by in vivo high-resolution optical imaging of intrinsic signals,” Proc. National Academy of Sciences 87 6082 - 6086    DOI : 10.1073/pnas.87.16.6082
Scott N. A. , Murphy T. H. (2012) “Hemodynamic responses evoked by neuronal stimulation via channelrhodopsin-2 can be independent of intracortical glutamatergic synaptic transmission,” PloS One 7 e29859 -    DOI : 10.1371/journal.pone.0029859
Harrison T. C. , Sigler A. , Murphy T. H. (2009) “Simple and cost-effective hardware and software for functional brain mapping using intrinsic optical signal imaging,” Journal of Neuroscience Methods 182 211 - 218    DOI : 10.1016/j.jneumeth.2009.06.021
Bouchard M. B. , Chen B. R. , Burgess S. A. , Hillman E. M. (2009) “Ultra-fast multispectral optical imaging of cortical oxygenation, blood flow, and intracellular calcium dynamics,” Opt. Express 17 15670 - 15678    DOI : 10.1364/OE.17.015670
Cheng H. , Luo Q. , Zeng S. , Chen S. , Cen J. , Gong H. (2003) “Modified laser speckle imaging method with improved spatial resolution,” J. Biomed. Opt. 8 559 - 564    DOI : 10.1117/1.1578089
Chen-Bee C. H. , Agoncillo T. , Xiong Y. , Frostig R. D. (2007) “The triphasic intrinsic signal: implications for functional imaging,” J. Neurosci. 27 4572 - 4586    DOI : 10.1523/JNEUROSCI.0326-07.2007
Malonek D. , Dirnagl U. , Lindauer U. , Yamada K. , Kanno I. , Grinvald A. (1997) “Vascular imprints of neuronal activity: relationships between the dynamics of cortical blood flow, oxygenation, and volume changes following sensory stimulation,” Proc. National Academy of Sciences 94 14826 - 14831    DOI : 10.1073/pnas.94.26.14826
Lim D. H. , Mohajerani M. H. , LeDue J. , Boyd J. , Chen S. , Murphy T. H. (2012) “In vivo large-scale cortical mapping using channelrhodopsin-2 stimulation in transgenic mice reveals asymmetric and reciprocal relationships between cortical areas,” Frontiers in Neural Circuits 6
Li M. , Miao P. , Zhu Y. , Tong S. (2011) “Functional laser speckle imaging of cerebral blood flow under hypothermia,” J. Biomed. Opt. 16 086011 -    DOI : 10.1117/1.3610995