High-density speckle contrast optical tomography of cerebral blood flow response to functional stimuli in the rodent brain

Abstract. Noninvasive, three-dimensional, and longitudinal imaging of cerebral blood flow (CBF) in small animal models and ultimately in humans has implications for fundamental research and clinical applications. It enables the study of phenomena such as brain development and learning and the effects of pathologies, with a clear vision for translation to humans. Speckle contrast optical tomography (SCOT) is an emerging optical method that aims to achieve this goal by directly measuring three-dimensional blood flow maps in deep tissue with a relatively inexpensive and simple system. High-density SCOT is developed to follow CBF changes in response to somatosensory cortex stimulation. Measurements are carried out through the intact skull on the rat brain. SCOT is able to follow individual trials in each brain hemisphere, where signal averaging resulted in comparable, cortical images to those of functional magnetic resonance images in spatial extent, location, and depth. Sham stimuli are utilized to demonstrate that the observed response is indeed due to local changes in the brain induced by forepaw stimulation. In developing and demonstrating the method, algorithms and analysis methods are developed. The results pave the way for longitudinal, nondestructive imaging in preclinical rodent models that can readily be translated to the human brain.

High-density speckle contrast optical tomography of cerebral blood flow response to functional stimuli in the rodent brain Tanja Dragojević, a, * , † Ernesto E. Vidal Rosas, a, † Joseph L. Hollmann, a Joseph P. Culver, b,c Carles Justicia, d,e and Turgut Durduran a,f 1 Introduction Microvascular cerebral blood flow (CBF) dynamics is an important biomarker of healthy brain development, pathology, and injury since it is a measure of the delivery of nutrients necessary for tissue metabolism and of the clearance of the by-products. The coupling between the neuronal activity and regional CBF is regarded as fundamental in the study of functional brain activity and is the basis of many functional neuroimaging techniques. 1,2 For example, it is responsible for the blood oxygen level-dependent (BOLD) signal that is commonly employed in functional magnetic resonance imaging (fMRI). Wide-field superficial 3 and deep tissue (>millimeters) diffuse optical 4-6 methods provide similar contrasts, and in comparison to large imaging modalities, their advantages are portability, fast and continuous monitoring, the use of nonionizing radiation, and being noninvasive. The main contrast to optical methods has traditionally been that of the oxy-and deoxy-hemoglobin concentration 7,8 but various methods utilizing laser speckle statistics have allowed for the direct measurement of blood flow. [9][10][11] Most common optical techniques for blood flow monitoring are based on the dynamic scattering of light and include laser Doppler flowmetry, 12,13 laser speckle contrast imaging, 9,11 and diffuse correlation spectroscopy (DCS). 14 Of these methods, DCS uses long-coherence laser light with point sources and displaced point detectors to measure the temporal fluctuations of a diffuse speckle which is translated into a blood flow index (BFI) from the deep tissues using model-based algorithms. 10 DCS is very attractive for CBF monitoring because its sensitivity can reach up to several centimeters. However, it also has several disadvantages, such as sophisticated hardware requirements and the associated lack of scalability thus reducing the spatial resolution and the ability to localize the signals. This is mainly because the high-density source-detector arrangements involving parallel acquisition from many source-detector pairs over a large range of separations are needed to achieve better resolution and location. This was, in fact, crucial to the success of diffuse optical tomography (DOT). [15][16][17] This is not, so far, accessible for diffuse correlation tomography (DCT). [18][19][20] To date, DCT has been limited to pilot studies on small animals 18,19 and coarse measurements on humans. 20 Recently, speckle contrast optical spectroscopy/tomography (SCOS/SCOT) was introduced to overcome some of these limitations. [21][22][23][24] It combines single-point illumination with a two-dimensional detector array to probe deep tissues. This method has successfully been used to monitor relative cerebral blood flow (rCBF) due to changes in a functional task in the adult human brain. 25 SCOT can achieve scalability using the speckle contrast with an appropriate correlation diffusion model to relax the requirement, which the data-type used by DCT has, i.e., relatively fast (tens of nanoseconds) output from the detectors, thus allowing for readily available detector arrays to be utilized in the measurements. [25][26][27][28] Here we highlight speckle contrast optical tomography (SCOT) which, by taking advantage of the scalability of SCOS detection, achieves tomographic measurements using scanning sources and array detectors. 27,28 Our recent developments in mice revealed comparable results to magnetic resonance imaging in following the infarct induced by CBF reduction. 27 A similar development 29 has been successfully used in the intraoperative room, where the three-dimensional (3-D) distribution of blood flow from a mastectomy skin flap was resolved. These applications highlight the potential of tomographic speckle contrast imaging in the clinic and in biomedical research. However, these examples are either too slow or use a too sparse array of detectors/sources to push the technique to its full potential.
On the other hand, traditional DOT has been pushed closer to its full potential by the introduction of high density (HD-DOT) and fast scanning. 15,17,30 HD-DOT allows increased spatial resolution and improves the accuracy of the detected changes by utilizing many, densely spaced overlapping samples of the probed tissue volume. By being fast, it also allows for the monitoring of the evolution of the functional changes and the regression of the confounding systemic signals. HD-DOT was shown to approach the resolution and accuracy of BOLD-fMRI in the adult brain. 31 However, HD-DOT measures the change in concentration of oxy-and deoxy-hemoglobin during functional activation in small animals and humans 15,30,32 without retrieving CBF information. In HD-DOT, CBF information can be inferred by the use of specific challenges or contrast agents, thus making the measurements complex.
Given the above developments, the stage is now set for the introduction of fast, high-density SCOT (HD-SCOT) to achieve high-density 3-D blood flow imaging. The first step toward that end is to apply HD-SCOT in the rat brain, transcranially, without thinning the skull. This is motivated by the considerable interest in small animal imaging, since they allow the investigation of neurovascular coupling in the normal brain, during the progression of diseases (e.g., Alzheimer's) or severe pathological events (e.g., stroke). Furthermore, HD-SCOT would allow to follow up on the effects of therapeutic drugs, the evolution of treatment and the recovery of brain function 33 over longer time periods (months to years). The measurement of the same parameter was considered by DCS in longitudinal monitoring of rodents, 34 demonstrating that the growth of the animal can impact the measurements in a detrimental manner unless a large field-ofview (FOV) measurement with tomographic reconstruction is achieved. In other words, both the magnitude and the location of the signal changes over time in a rapidly aging small animal model.
In this study, we have developed and demonstrated the capabilities of SCOT for the monitoring of the cerebral blood flow response due to functional forepaw stimulation-induced changes in the somatosensory cortex of the rat brain. Tomographic images with good spatial resolution were possible due to the availability of high-density data. Furthermore, in order to increase the temporal resolution, we have introduced a technique based on a moving window scheme, which effectively improves the sampling efficiency.

Animal Preparation
The experimental protocol was approved by the ethical committee of the University of Barcelona. Adult male Wistar rats (n ¼ 5) were anesthetized with 4% isoflurane and moved to a stereotaxic frame. The anesthesia was reduced to 2% isoflurane and the surgical procedure was performed, that consisted of the retraction of the skin and underlying muscle while keeping the skull intact, to avoid confounding effects. We note that this is not a necessary condition for the method and it could be applied fully noninvasively. 35 A well surrounding the imaging area was built with dental cement and was filled with mineral oil to prevent the skull from drying out. After surgery, the isoflurane infusion was stopped and anesthesia was maintained with a continuous intravenous administration of α-chloralose (40 mg∕kg∕h).

Functional Activation Paradigm
Electrical stimulation was achieved with a pair of subdermal needle electrodes that were inserted into right and left forepaws. The experiment consisted of electrically stimulating each forepaw individually with 200-μs pulses at 6-Hz repetition rate and an amplitude of 2 mA. The stimulation paradigm comprised of 18-s baseline, 30-s electrical stimulation, and 90-s recovery. This was repeated eight times for each forepaw. 16 Once the functional protocol was completed, four repetitions of sham measurements were taken for the same duration, during which no electrical stimulation was performed, but the subdermal needles were kept inside forepaws.

Magnetic Resonance Imaging
To compare SCOT, functional MRI was acquired for the same animals (n ¼ 5) under α-chloralose (40 mg∕kg∕h) anesthesia within 24 to 72 h after optical data were acquired. Experiments were conducted in a 7T MRI magnet (Bruker BioSpin, Germany). Coronal multislice spin-echo echo-planar images were acquired with a repetition time (TR) of 3 s and a bandwidth of 150 kHz. In total, five slices with 2-mm thickness were acquired with an FOV of 39.9 mm × 39.9 mm and a pixel size of 0.208 mm × 0.208 mm. Additionally, anatomical high-resolution T2-weighted images were obtained to align the fMRI and SCOT reconstructions (see Sec. 2.6). Anatomical data consisted of forty slices with 0.08-mm thickness, FOV of 30 mm × 30 mm, and pixel size 0.117 mm × 0.117 mm.

Speckle Contrast Optical Tomography Instrumentation
The instrumentation was described in a previous publication 27 and we review it in brief. The main components of the SCOT system are displayed in Fig. 1(a), where the scanning module consisted of a temperature-controlled continuous-wave laser at 785 nm (Thorlabs, Germany) that was focused by aspheric lenses onto a pair of galvo-mirrors (Nutfield), which were directing the light beam onto the rat head. 27 The reflected diffuse light was captured by a CCD camera (Basler, Germany, scA640-120 fm).

Optical Data Acquisition Protocol
Once the surgery was completed, the rat was moved under the imaging arm of the setup as shown in Fig. 1(a). During the alignment of the imaging domain, the central source of the array was adjusted to coincide with the location of a Bregma point [ Fig. 1(b)]. The number of source positions was fixed at nine for all animals organized in a 3 × 3 arrays [ Fig. 1(b)]. A calibration image was then taken where the physical imaged area per pixel was calibrated using a ruler placed on top of the rat's head. The f-number of the optical system was chosen to match the pixel to speckle size. Then the dental cement well was filled with mineral oil. Additionally, a green light LED (525 nm) was used to illuminate the rat's head and verify that the imaging area was free of contaminants (i.e., hair and bubbles). Otherwise, the mineral oil was withdrawn, the head was cleaned, and the oil was reapplied. Furthermore, to obtain the baseline blood flow, 20 exposure times that were logarithmically spaced between 50 μs and 10 ms were used to recover the baseline CBF.
Then the stimulation protocol (Sec. 2.2) was conducted while SCOT data were being acquired. Images were acquired with an exposure time of 1 ms and an FOV of 1.2 cm × 1 cm. One full tomographic data set consisted of the acquisition of a 100 images per source location, which was completed in 6 s.
Calculation of the speckle contrast was done in the temporal approach. 36 In the end, to account for the variance of the sensor dark noise, 900 images with an exposure time of 1 ms were collected with the laser light off.

Theory of speckle contrast
The speckle contrast squared κ 2 is defined as the ratio of the variance σ 2 to the squared mean value of the intensity I 2 as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 6 3 ; 1 5 5 κ 2 ðr; TÞ ¼ σ 2 ðr; TÞ Ihðr; TÞi 2 ; (1) where T is the exposure time and r is the source-detector separation. The mean Ihðr; TÞi 2 and the variance can be computed either temporally (over independent speckles in a single pixel over time) or spatially (over independent speckles in different, neighboring pixels). 36,37 Furthermore, the speckle contrast can be related to the normalized electric field autocorrelation function g 1 by E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 3 2 6 ; 5 0 7 where τ is the correlation delay time, and β is an optical system parameter proportional to the number modes/speckles in the detection system, which is ∼0.5 for unpolarized light used in this study. This equation is related to the optical parameters of the medium and BFI through modeling of the g 1 function using the correlation diffusion equation. 10 These methods were previously validated for SCOS/SCOT in phantoms, small animal models and adult, human subjects. 21,28

Speckle contrast correction and data analysis
Prior to the solution of the inverse problem for the tomographic reconstruction, the speckle contrast was corrected to account for the shot and dark noise variations associated with the camera sensor. The corrected speckle contrast is defined as 20 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 3 2 6 ; 2 8 7 where I c ¼ I − I d is the dark corrected intensity, I d and σ 2 d are the mean and variance of the dark noise over a number of frames, and σ 2 s denotes the variance of the estimated shot noise. In addition, to improve the signal-to-noise ratio (SNR) and to reduce the computational load, the speckle contrast of four adjacent pixels was binned into a single pixel [ Fig. 2(b)]. 27 The baseline BFI was determined using multi-exposure speckle imaging (MESI), 26,38 using the multi-exposure data that was acquired as described in Sec. 2.5 and Ref. 37. The MESI, together with Eq. (2), was used to fit for an averaged BFI 0 using the dependence of the speckle contrast on the exposure time as previously described for SCOS. 21,25,26 Valid source-detector pairs were selected based on two criteria, where first was that source-detector pairs with separations shorter than 0.15 cm were disregarded since the diffusion model fails to predict the transport of photon at short separations, 39 also, pixels near the source illumination point were discarded since they were saturated. In the future, a higher order approximation to the radiation transport equation can be used to retain all the nonsaturated data. Second, source-detector pairs with a separation larger than 0.65 cm were discarded to avoid contamination or biasing of the signals due to the limited FOV, where many detectors were reaching the dental cement surrounding well and the lower SNR.
where J is the Jacobian, κ 2 0;c is the baseline (Sec. 2.6.4), corrected speckle contrast squared, and κ 2 c is the corrected speckle contrast squared that was altered due to a perturbation. Under this approximation, the temporal changes in the measured speckle contrast are linearly related to spatial variations in the BFI (ΔBFI) in a 3-D space. The tomographic image of these values can be obtained by inverting the above relationship. Since the inversion process is an ill-posed problem, we have regularized the inversion using a spatially varying Tikhonov regularization scheme. 27,28,40 An illustration of the imaging domain is displayed in Fig. 2(a), which consisted of 5,280 voxels (22 × 20 × 12). A sample image of the input measurement data is shown in Fig. 2(b) together with the source locations (red circles) and the voxel positions. ΔrCBFðtÞ ¼ ΔBFIðtÞ BFI 0 × 100ð%Þ: To quantify the spatial extent of the activation at the depth with the largest change, the area of activation was defined as the number of pixels with rCBF >50% of the maximum response. 11 Additionally, the recovered ΔrCBF images were compared to the MRI data by co-registering the anatomical MRI and optical images using the Bregma point as a reference. Then optical data were transformed into the MRI coordinates taking the curvature of the brain into account for better visualization. 27 2.6.5 Sliding tomographic reconstruction scheme to improve the temporal resolution The standard reconstruction scheme, described in Fig. 3(a), involves solving the inverse problem once a full tomographic data set is available, i.e., after all sources have completed a scanning cycle. Then the process is repeated for subsequent scanning cycles. One disadvantage of this method is that rapidly evolving changes, such as those during somatosensory activation, of the signals are averaged over due to having a tomographic image every few seconds.
Here a reconstruction scheme based on a moving window was employed to improve the temporal resolution. It is similar to the jittering of the RT used in fMRI to improve the efficiency of the sampling. 41 The process is displayed in Fig. 3(b), where a tomographic image can be obtained as soon as data from a new source position is available. This is possible on the basis that the reference source can be arbitrarily assigned and the necessary conditions for the inversion are fulfilled as long as data from every source is available [ Fig. 3(b)]. Furthermore, since the time delay for switching between scanning positions is considerably smaller than the acquisition time, i.e., the scanning process can be considered continuous regardless of the selection of the source as a reference. The method does not produce independent samples, but rather a running average thus the reconstructions include data across different stages, Fig. 3(c) (dotted circles). In Fig. 3(c), the dense imaging achieved by the sliding window reconstruction scheme (void circles) is compared with the sparse temporal imaging obtained with the standard approach (filled circles). Since the "contamination" from samples, i.e., source positions, from different stages introduce artifacts due to the averaging of data at different conditions, we have also defined a period of influence to indicate the time points were the contribution of data from the stimulation stage is <50%. This period of influence is illustrated with the shaded regions before and after the stimulation stage in Fig. 3(c).

Data Conditioning Results
Before image reconstruction, the usable source-detector pairs needed to be identified as described in Sec. 2.6.2. As Fig. 4(a) shows, the speckle contrast data from source one, which displays several features including the location of the source, large vessels and the effect of the dental cement. In particular, the speckle contrast along the horizontal red line is plotted in Fig. 4(b), where the effect of the source is visible. Fig. 4 (a) Speckle contrast data from source one. (b) Speckle contrast squared as a function of the source-detector distance. The speckle contrast squared at the source location drops close to zero because these detectors were saturated where in this case, the variance is zero. As described in Sec. 2.6.2, only source-detector pairs with separations in the range of 0.15 to 0.65 cm were kept before proceeding to the inversion process. On average, considering the nine sources, a total of ∼50;000 valid source-detector pairs were used for image reconstruction for each animal. Figure 5(a) depicts a sample of the differential speckle contrast (Δκ 2 c ) during the stimulation stage from the source indicated with the red circle while the remaining source locations are shown in white. Particularly, there is a region indicated with the colored squares that shows a reduction of the speckle contrast during the stimulation stage. The relative temporal changes of Δκ 2 c ðtÞ with respect to baseline (κ 2 0;c ) at these locations are displayed in Fig. 5(b). A decrease of the speckle contrast during the stimulation can be observed followed by a return to the baseline. In other words, as expected, these changes reveal an increase of the CBF during the stimulation stage. Figure 6(a) shows representative volumetric images of rCBF during left forepaw stimulation. Each column corresponds to different time points, whereas rows indicate slices in depth. The first column during the baseline period shows ΔrCBF maps that are relatively homogeneous. However, during the stimulation stage, an active region at y ¼ −0.05 cm with higher ΔrCBF is distinguishably followed by maps similar to those at baseline. The locations of apparent changes of ΔrCBF at the skull level (y ¼ 0 cm) are due to the regularization operation which produces smooth responses rather than sharp changes. The time course of the activation region (0.1 cm × 0.1 cm) centered at the peak value, displayed in Fig. 6(b), demonstrates the increased blood flow during the stimulation stage. The onset of the increase of ΔrCBF before the start of the stimulation is an artifact produced by the sliding average effect, which covers data within the baseline and activation periods as described above.

Comparison to Functional Magnetic Resonance Imaging
The registration of the MRI and the optical data allowed a direct comparison of the location and the spatial extent of the activation. Figures 7(a) and 7

(b) show representative BOLD fMRI and
ΔrCBF images of an animal stimulated on the left forepaw. These images demonstrate that the activation is correctly located by speckle contrast. Furthermore, following the transformation of the recovered ΔrCBF into the MRI coordinates, Fig. 7(c) demonstrates that the maximum change of ΔrCBF occurs on the cortex region and not on the overlying skull. Similar images were obtained for the right forepaw stimulation [Figs. 7(d)-7(f)].  Fig. 8(d), where no activation but flat responses can be seen during the stimulation period. Table 1 summarizes the estimated activation area for each individual at the slice with the maximum change.

Discussion
In this work, we have demonstrated three-dimensional monitoring of the cerebral blood flow response due to functional stimuli using high-density SCOT through the intact skull of the rat's head. The results show that the location and the spatial extent of the activation maps are in agreement with BOLD fMRI studies under a similar stimulation paradigm. 42 Figure 6(a) shows a clear distinction between the different stages of the experiment and more importantly, the maximum change is correctly located at the cortex layer (y ¼ −0.05 cm). 42 The time trace of ΔrCBF [ Fig. 6(b)] shows a relatively rapid and sustained increase during the stimulation period. However, care must be taken to interpret the apparent increase before the start of the stimulation, which is introduced by the moving window reconstruction scheme. In the future, the temporal resolution can be improved by calculating the speckle contrast over space instead of time in sequential frames using a denser detection array structure while maintaining the high frame rate. If our method is further utilized, its effects can be minimized by considering an average of the forward and backward moving window, by explicitly accounting for the moving window using a weighing function in the postprocessing and by other methods  that are utilized in similar windowing schemes. This could be the topic of further investigation in the future. In this study, we have opted for the current detector array for simplicity. In fact, speckle data can be acquired at even higher frame rates in animals 43 and in adult humans 44 that may in the future allow the consideration of the estimation of other biomarkers in addition to CBF. For example, with the improved temporal resolution, it would be possible to map the pulsatility of the cerebral blood flow due to cardiac cycle as a marker for vascular stiffness 43 and intracranial pressure. 45 The spatial extent of the activation listed in Table 1 is in agreement with similar studies, 11 and also with the size of the activation reported by other diffuse optical modalities, for example, the spatial extent obtained with SCOT is 11.7 AE 1.7 mm 2 , which is very close to the activation area obtained with DOT  which is 12.3 AE 2 mm 2 . 16 However, to show the volumetric change instead of the activation area, one would require a statistical analysis that would take into account factors such as differential effects or regularization at different slides and the hemodynamic response function. 46 Figure 7 provides further confirmation of the accuracy of the SCOT technique in the location and the spatial extent of the activation. In particular, Figs. 7(c) and 7(f) display the maximum ΔrCBF at the cortex and a decrease as the depth increases. However, the smooth and continuous decrease of ΔrCBF at deeper layers (y > 0.15 cm) is partly a consequence of the regularization method. The purpose of regularization is to stabilize the inversion process by finding a trade-off between the accuracy of the solution and the influence of noise, where this is achieved by attenuating the effect of higher-order modes, which result in smoother solutions. Additionally, it is well-known that recovered perturbations in the reflectance geometry tend to be located closer to the surface. 40 Sharper changes can be obtained by adapting fully nonlinear iterative schemes which are readily available in open access software such as NIRFAST or TOAST for SCOT 47,48 which would be the topic for further research. SCOT provides a unique combination of deep tissue sensitivity together with a relatively high spatial resolution for the tomography of CBF changes. The possibility of high-density measurements contributed to an improvement in spatial resolution of 1 mm. Similar developments in DOT 49 provided substantial improvements in the spatial resolution and in the minimization of localization errors. High-density measurements allow the sampling of the brain tissue with a large number of overlapping measurements 50 as a result, each point in the imaging space is better constrained and also less sensitive to noise during the inversion process. On the contrary, high-density data do not scale well in DCS/DCT because the monetary cost is prohibitively large since it requires expensive fast detectors for every added detector position. 18,19,23 Figure 8 shows the performance of SCOT for the imaging of CBF in detail. As it is very typical of neuroimaging methods, intertrial, intra-animal, and inter-animal averaging were utilized. This was originally demonstrated for the measurement of CBF due to forepaw stimulation in Ref. 51, where they have characterized different sources of variability. In particular, Fig. 8(a) displays the individual responses for the left forepaw stimulation. The first point to note is that SCOT can characterize changes due to individual trials with an agreeable contrast and SNR. The variability is partly due to the variability introduced by the measurement method but also due to physiological variations in response to each trial including potential habituation effects. It is important to emphasize that the sharp increase in ΔrCBF at the beginning of the stimulation period occurs in all trials and contributes to obtaining a robust peak in the averaged response (black line). Figure 8(b) shows the ΔrCBF after averaging over these trials for each animal again due to left forepaw stimulation. As expected, there are intra-animal variations presumably due to differences in anatomy, the dynamics of its physiology and condition at a given moment. However, the consistent initial overshoot followed by the persistent activation during the stimulation period, observed across all animals, demonstrates that SCOT is capable of imaging CBF in a robust manner. This characteristic is further supported by Fig. 8(c), where the grand average is plotted against the individual responses for the five animals and from the stimulation of both hemispheres. Baseline, activation, and recovery stages are markedly distinguishable from each other. For comparison, Fig. 8(d) shows the sham responses of the same animals, where mainly flat responses characterized the trends and allow one to appreciate the repeatability of the method. We note that a Wilcoxon test has demonstrated that this response reaches statistical significance in all animals during stimuli but not during sham. However, a full-fledged statistical analysis using a general linear model is beyond the scope of this work. Taken together, Figs. 8(a)-8(d) demonstrate the feasibility of tomographic imaging of changes in ΔrCBF from a single trial level to a population level, with a considerable improvement of the response due to the effect of averaging. This is not always true for all neuroimaging methods.
Future works using multimodal approaches, e.g., combining DOT with SCOT and/or different fMRI methods with SCOT can be envisioned. These methods would enable a quantitative comparison of the CBF response to that of the oxygenation and derive the changes in cerebral metabolic rate of oxygen extraction. We have not carried this out in this manuscript since the fMRI data were not taken simultaneously and the protocol was not optimized for this purpose.
As partially demonstrated in this work, SCOT allows for a range of advantages over traditional DCT. SCOT can be implemented with arrays of detectors using existing technologies, which, as mentioned in Sec. 1 and Refs. 22-24, 27-29, allows for the improvement of the SNR, FOV, implementation of high-density source-detector pairs, and for potentially low cost, wearable devices. 25 The simplicity of the data-type allows for embedding the calculations in firmware, 25 which, in turn, allows for robust devices that could be utilized in the field without relying on fragile software and computers. The usability of cameras 24,27,29 and camera-like arrays 30,35 readily enable noncontact scanning which can bring substantial benefits to imaging. For example, it allows the collection of highly dense data at multiple locations with lesser mechanical limitations such as those due to the location and the physical diameter of the fiber cladding. 30 However, noncontact measurements have limitations and we note here, that contact measurements can be implemented with SCOT, 30,52 which would allow for measurements through the hair when translating to studies on humans. 30 Other implementations can also be envisioned combining spatial, spectral, and temporal information efficiently.
Our previous study demonstrated the possibility of tomographic brain imaging of stroke in a mouse model, 27 where the present work is an extension to a larger geometry, to deeper and more rapidly evolving signals and to a more challenging experimental paradigm. Combining both contributions, this technology provides new opportunities for the study of the rodent brain after an injury. For example, it would be possible to monitor not only pathological changes in the brain due to a stroke, but also to study the evolution of the functional brain activity after the injury 51,53,54 and in the long term, it would also be possible to study longitudinal effects 17,30,34 such as those due to rehabilitation. The unique innovation of HD-SCOT is that these studies can be carried out with the same imaging platform in a noninvasive manner. The unique innovation of HD-SCOT is that these analyses can be carried out with the same imaging platform.

Conclusion
The feasibility of functional brain monitoring using high-density SCOT has been demonstrated. Robust changes of ΔrCBF due to electrical stimulation of forepaw have been observed and the Joseph P. Culver is a professor of radiology at Washinton University. He obtained a PhD in physics at the University of Pennsylvania for work in ultrafast infrared laser spectroscopy. As a postdoc, he developed diffuse optical tomography systems for imaging breast cancer and imaging cerebal hemodynamics in animals. His lab explores ways of leveraging optical measurements for neuroimaging, with a focus on development of diffuse optical tomography techniques for mapping human brain function.