The Cousa objective: a long-working distance air objective for multiphoton imaging in vivo.

Multiphoton microscopy can resolve fluorescent structures and dynamics deep in scattering tissue and has transformed neural imaging, but applying this technique in vivo can be limited by the mechanical and optical constraints of conventional objectives. Short working distance objectives can collide with compact surgical windows or other instrumentation and preclude imaging. Here we present an ultra-long working distance (20 mm) air objective called the Cousa objective. It is optimized for performance across multiphoton imaging wavelengths, offers a more than 4 mm2 field of view with submicrometer lateral resolution and is compatible with commonly used multiphoton imaging systems. A novel mechanical design, wider than typical microscope objectives, enabled this combination of specifications. We share the full optical prescription, and report performance including in vivo two-photon and three-photon imaging in an array of species and preparations, including nonhuman primates. The Cousa objective can enable a range of experiments in neuroscience and beyond.

The development of bi-directionally coupled self-organizing neurovascular networks captures orientation-selective neural and hemodynamic cortical responses.

Networks of neurons are the primary substrate of information processing. Conversely, blood vessels in the brain are generally viewed to have physiological functions unrelated to information processing, such as the timely supply of oxygen, and other nutrients to the neural tissue. However, recent studies have shown that cerebral microvessels, like neurons, exhibit tuned responses to sensory stimuli. Tuned neural responses to sensory stimuli may be enhanced with experience-dependent Hebbian plasticity and other forms of learning. Hence, it is possible that the microvascular network might also be subject to some form of competitive learning rules during early postnatal development such that its fine-scale structure becomes optimized for metabolic delivery to a given neural micro-architecture. To explore the possibility of adaptive lateral interactions and tuned responses in cerebral microvessels, we modelled the cortical neurovascular network by interconnecting two laterally connected self-organizing networks. The afferent and lateral connections of the neural and vascular networks were defined by trainable weights. By varying the topology of lateral connectivity in the vascular network layer, we observed that the partial correspondence of feature selectivity between neural and hemodynamic responses could be explained by lateral coupling across local blood vessels such that the central domain receives an excitatory drive of more blood flow and a distal surrounding region where blood flow is reduced. Critically, our simulations suggest a new role for feedback from the vascular to the neural network because the radius of vascular perfusion determines whether the cortical neural map develops into a clustered vs. salt-and-pepper organization.

Invited Session IV: Studies of the visual cortex with sub-millimeter resolution: Three-photon imaging reveals the neural basis of fMRI across cortical layers.

V1 is ideally suited to study the spatial organization of neurovascular coupling at the level of synapses, neurons, individual blood vessels and laminar-resolution fMRI. This is because at least in layers 1 and 2/3 of V1, the functional micro-architecture for neurons, synapses and blood vessels has been determined using 2-photon imaging (Ohki et al 2005 Nature; Kara and Boyd 2009 Nature; O'Herron et al 2016 Nature). Hence, feature selectivity, e.g., orientation and direction selectivity of spiking, synaptic and hemodynamic activity in layer 2/3 is known. However, the micro-architecture of layer 4 neural activity (spiking and synaptic) along with individual blood vessel responses is unknown because conventional 2-photon imaging cannot access deeper cortical layers. The organizing principles of neural maps and the selectivity of hemodynamic responses is of paramount importance for laminar processing because the thalamic inputs arriving into layer 4 are untuned. 3-photon imaging triples the imaging depth compared to 2-photon imaging. Using this optical technique and high-resolution fMRI, we have determined the extent to which different types of neural (spiking, synaptic) and vascular signals (blood flow from individual vessels and fMRI voxels) are coupled across cortical layers in the cat V1. Our data show systematic changes in selectivity of hemodynamic signals across cortical layers that have clear underpinnings in neural circuitry and the propagation of hemodynamic signals.

Cortical layer-specific differences in stimulus selectivity revealed with high-field fMRI and single-vessel resolution optical imaging of the primary visual cortex.

The mammalian neocortex exhibits a stereotypical laminar organization, with feedforward inputs arriving primarily into layer 4, local computations shaping response selectivity in layers 2/3, and outputs to other brain areas emanating via layers 2/3, 5 and 6. It cannot be assumed a priori that these signatures of laminar differences in neuronal circuitry are reflected in hemodynamic signals that form the basis of functional magnetic resonance imaging (fMRI). Indeed, optical imaging of single-vessel functional responses has highlighted the potential limits of using vascular signals as surrogates for mapping the selectivity of neural responses. Therefore, before fMRI can be employed as an effective tool for studying critical aspects of laminar processing, validation with single-vessel resolution is needed. The primary visual cortex (V1) in cats, with its precise neuronal functional micro-architecture, offers an ideal model system to examine laminar differences in stimulus selectivity across imaging modalities. Here we used cerebral blood volume weighted (wCBV) fMRI to examine if layer-specific orientation-selective responses could be detected in cat V1. We found orientation preference maps organized tangential to the cortical surface that typically extended across depth in a columnar fashion. We then examined arterial dilation and blood velocity responses to identical visual stimuli by using two- and three- photon optical imaging at single-vessel resolution-which provides a measure of the hemodynamic signals with the highest spatial resolution. Both fMRI and optical imaging revealed a consistent laminar response pattern in which orientation selectivity in cortical layer 4 was significantly lower compared to layer 2/3. This systematic change in selectivity across cortical layers has a clear underpinning in neural circuitry, particularly when comparing layer 4 to other cortical layers.

Neural correlates of single-vessel haemodynamic responses in vivo.

Neural activation increases blood flow locally. This vascular signal is used by functional imaging techniques to infer the location and strength of neural activity. However, the precise spatial scale over which neural and vascular signals are correlated is unknown. Furthermore, the relative role of synaptic and spiking activity in driving haemodynamic signals is controversial. Previous studies recorded local field potentials as a measure of synaptic activity together with spiking activity and low-resolution haemodynamic imaging. Here we used two-photon microscopy to measure sensory-evoked responses of individual blood vessels (dilation, blood velocity) while imaging synaptic and spiking activity in the surrounding tissue using fluorescent glutamate and calcium sensors. In cat primary visual cortex, where neurons are clustered by their preference for stimulus orientation, we discovered new maps for excitatory synaptic activity, which were organized similarly to those for spiking activity but were less selective for stimulus orientation and direction. We generated tuning curves for individual vessel responses for the first time and found that parenchymal vessels in cortical layer 2/3 were orientation selective. Neighbouring penetrating arterioles had different orientation preferences. Pial surface arteries in cats, as well as surface arteries and penetrating arterioles in rat visual cortex (where orientation maps do not exist), responded to visual stimuli but had no orientation selectivity. We integrated synaptic or spiking responses around individual parenchymal vessels in cats and established that the vascular and neural responses had the same orientation preference. However, synaptic and spiking responses were more selective than vascular responses--vessels frequently responded robustly to stimuli that evoked little to no neural activity in the surrounding tissue. Thus, local neural and haemodynamic signals were partly decoupled. Together, these results indicate that intrinsic cortical properties, such as propagation of vascular dilation between neighbouring columns, need to be accounted for when decoding haemodynamic signals.

The shape of dendritic arbors in different functional domains of the cortical orientation map.

The neocortex is organized into macroscopic functional maps. However, at the microscopic scale, the functional preference and degree of feature selectivity between neighboring neurons can vary considerably. In the primary visual cortex, adjacent neurons in iso-orientation domains share the same orientation preference, whereas neighboring neurons near pinwheel centers are tuned to different stimulus orientations. Moreover, several studies have found greater orientation selectivity in iso-orientation domains than in pinwheel centers. These differences suggest that neurons sample local inputs in a spatially homogenous fashion and independently of the location of their soma on the orientation map. Here we determine whether dendritic geometry is affected by neuronal position on the orientation map. We labeled individual layer 2/3 pyramidal neurons with fluorescent dyes in specific domains of the orientation map in cat primary visual cortex and imaged their dendritic trees in vivo by two-photon microscopy. We found that the circularity and uniformity of dendritic trees is independent of somatic position on the orientation map. Moreover, the dendrites of neurons located close to pinwheel centers extend across all orientation domains in an unbiased fashion. Thus, unbiased dendritic trees appear to provide an anatomical substrate for the systematic variations in feature selectivity across the orientation map.

An artery-specific fluorescent dye for studying neurovascular coupling.

We demonstrate that Alexa Fluor 633 hydrazide (Alexa Fluor 633) selectively labels neocortical arteries and arterioles by binding to elastin fibers. We measured sensory stimulus-evoked arteriole dilation dynamics in mouse, rat and cat visual cortex using Alexa Fluor 633 together with neuronal activity using calcium indicators or blood flow using fluorescein dextran. Arteriole dilation decreased fluorescence recorded from immediately underlying neurons, representing a potential artifact during neuronal functional imaging experiments.

A micro-architecture for binocular disparity and ocular dominance in visual cortex.

In invertebrate predators such as the praying mantis and vertebrate predators such as wild cats the ability to detect small differences in inter-ocular retinal disparities is a critical means for accurately determining the depth of moving objects such as prey. In mammals, the first neurons along the visual pathway that encode binocular disparities are found in the visual cortex. However, a precise functional architecture for binocular disparity has never been demonstrated in any species, and coarse maps for disparity have been found in only one primate species. Moreover, the dominant approach for assaying the developmental plasticity of binocular cortical neurons used monocular tests of ocular dominance to infer binocular function. The few studies that examined the relationship between ocular dominance and binocular disparity of individual cells used single-unit recordings and have provided conflicting results regarding whether ocular dominance can predict the selectivity or sensitivity to binocular disparity. We used two-photon calcium imaging to sample the response to monocular and binocular visual stimuli from nearly every adjacent neuron in a small region of the cat visual cortex, area 18. Here we show that local circuits for ocular dominance always have smooth and graded transitions from one apparently monocular functional domain to an adjacent binocular region. Most unexpectedly, we discovered a new map in the cat visual cortex that had a precise functional micro-architecture for binocular disparity selectivity. At the level of single cells, ocular dominance was unrelated to binocular disparity selectivity or sensitivity. When the local maps for ocular dominance and binocular disparity both had measurable gradients at a given cortical site, the two gradient directions were orthogonal to each other. Together, these results indicate that, from the perspective of the spiking activity of individual neurons, ocular dominance cannot predict binocular disparity tuning. However, the precise local arrangement of ocular dominance and binocular disparity maps provide new clues regarding how monocular and binocular depth cues may be combined and decoded.

High fidelity sensory-evoked responses in neocortex after intravenous injection of genetically encoded calcium sensors.

Two-photon imaging of genetically-encoded calcium indicators (GECIs) has traditionally relied on intracranial injections of adeno-associated virus (AAV) or transgenic animals to achieve expression. Intracranial injections require an invasive surgery and result in a relatively small volume of tissue labeling. Transgenic animals, although they can have brain-wide GECI expression, often express GECIs in only a small subset of neurons, may have abnormal behavioral phenotypes, and are currently limited to older generations of GECIs. Inspired by recent developments in the synthesis of AAVs that readily cross the blood brain barrier, we tested whether an alternative strategy of intravenously injecting AAV-PhP.eB is suitable for two-photon calcium imaging of neurons over many months after injection. We injected young (postnatal day 23 to 31) C57BL/6J mice with AAV-PhP.eB-Synapsin-jGCaMP7s via the retro-orbital sinus. After allowing 5 to 34 weeks for expression, we performed conventional and widefield two-photon imaging of layers 2/3, 4 and 5 of the primary visual cortex. We found reproducible trial-by-trial neural responses and tuning properties consistent with known feature selectivity in the visual cortex. Thus, intravenous injection of AAV-PhP.eB does not interfere with the normal processing in neural circuits. In vivo and histological images show no nuclear expression of jGCaMP7s for at least 34 weeks post-injection.

High fidelity sensory-evoked responses in neocortex after intravenous injection of genetically encoded calcium sensors.

Two-photon imaging of genetically-encoded calcium indicators (GECIs) has traditionally relied on intracranial injections of adeno-associated virus (AAV) or transgenic animals to achieve expression. Intracranial injections require an invasive surgery and result in a relatively small volume of tissue labeling. Transgenic animals, although they can have brain-wide GECI expression, often express GECIs in only a small subset of neurons, may have abnormal behavioral phenotypes, and are currently limited to older generations of GECIs. Inspired by recent developments in the synthesis of AAVs that readily cross the blood brain barrier, we tested whether an alternative strategy of intravenously injecting AAV-PHP.eB is suitable for two-photon calcium imaging of neurons over many months after injection. We injected C57BL/6 J mice with AAV-PHP.eB-Synapsin-jGCaMP7s via the retro-orbital sinus. After allowing 5 to 34 weeks for expression, we performed conventional and widefield two-photon imaging of layers 2/3, 4 and 5 of the primary visual cortex. We found reproducible trial-by-trial neural responses and tuning properties consistent with known feature selectivity in the visual cortex. Thus, intravenous injection of AAV-PHP.eB does not interfere with the normal processing in neural circuits. In vivo and histological images show no nuclear expression of jGCaMP7s for at least 34 weeks post-injection.

Highly ordered arrangement of single neurons in orientation pinwheels.

In the visual cortex of higher mammals, neurons are arranged across the cortical surface in an orderly map of preferred stimulus orientations. This map contains 'orientation pinwheels', structures that are arranged like the spokes of a wheel such that orientation changes continuously around a centre. Conventional optical imaging first demonstrated these pinwheels, but the technique lacked the spatial resolution to determine the response properties and arrangement of cells near pinwheel centres. Electrophysiological recordings later demonstrated sharply selective neurons near pinwheel centres, but it remained unclear whether they were arranged randomly or in an orderly fashion. Here we use two-photon calcium imaging in vivo to determine the microstructure of pinwheel centres in cat visual cortex with single-cell resolution. We find that pinwheel centres are highly ordered: neurons selective to different orientations are clearly segregated even in the very centre. Thus, pinwheel centres truly represent singularities in the cortical map. This highly ordered arrangement at the level of single cells suggests great precision in the development of cortical circuits underlying orientation selectivity.

From connectomes to projectomes: an interview with Clay Reid.

Neurophotonics Associate Editor Prakash Kara (Univ. of Minnesota) interviewed his colleague Clay Reid (Allen Institute for Brain Science) about his pioneering work in neuroscience.

3D optogenetic control of arteriole diameter in vivo.

Modulation of brain arteriole diameter is critical for maintaining cerebral blood pressure and controlling regional hyperemia during neural activity. However, studies of hemodynamic function in health and disease have lacked a method to control arteriole diameter independently with high spatiotemporal resolution. Here, we describe an all-optical approach to manipulate and monitor brain arteriole contractility in mice in three dimensions using combined in vivo two-photon optogenetics and imaging. The expression of the red-shifted excitatory opsin, ReaChR, in vascular smooth muscle cells enabled rapid and repeated vasoconstriction controlled by brief light pulses. Two-photon activation of ReaChR using a spatial light modulator produced highly localized constrictions when targeted to individual arterioles within the neocortex. We demonstrate the utility of this method for examining arteriole contractile dynamics and creating transient focal blood flow reductions. Additionally, we show that optogenetic constriction can be used to reshape vasodilatory responses to sensory stimulation, providing a valuable tool to dissociate blood flow changes from neural activity.

Functional imaging with cellular resolution reveals precise micro-architecture in visual cortex.

Neurons in the cerebral cortex are organized into anatomical columns, with ensembles of cells arranged from the surface to the white matter. Within a column, neurons often share functional properties, such as selectivity for stimulus orientation; columns with distinct properties, such as different preferred orientations, tile the cortical surface in orderly patterns. This functional architecture was discovered with the relatively sparse sampling of microelectrode recordings. Optical imaging of membrane voltage or metabolic activity elucidated the overall geometry of functional maps, but is averaged over many cells (resolution >100 microm). Consequently, the purity of functional domains and the precision of the borders between them could not be resolved. Here, we labelled thousands of neurons of the visual cortex with a calcium-sensitive indicator in vivo. We then imaged the activity of neuronal populations at single-cell resolution with two-photon microscopy up to a depth of 400 microm. In rat primary visual cortex, neurons had robust orientation selectivity but there was no discernible local structure; neighbouring neurons often responded to different orientations. In area 18 of cat visual cortex, functional maps were organized at a fine scale. Neurons with opposite preferences for stimulus direction were segregated with extraordinary spatial precision in three dimensions, with columnar borders one to two cells wide. These results indicate that cortical maps can be built with single-cell precision.

Improving laser standards for three-photon microscopy.

Significance: Three-photon excitation microscopy has double-to-triple the penetration depth in biological tissue over two-photon imaging and thus has the potential to revolutionize the visualization of biological processes in vivo. However, unlike the plug-and-play operation and performance of lasers used in two-photon imaging, three-photon microscopy presents new technological challenges that require a closer look at the fidelity of laser pulses. Aim: We implemented state-of-the-art pulse measurements and developed innovative techniques for examining the performance of lasers used in three-photon microscopy. We then demonstrated how these techniques can be used to provide precise measurements of pulse shape, pulse energy, and pulse-to-pulse intensity variability, all of which ultimately impact imaging. Approach: We built inexpensive tools, e.g., a second harmonic generation frequency-resolved optical gating (SHG-FROG) device and a deep-memory diode imaging (DMDI) apparatus to examine laser pulse fidelity. Results: First, SHG-FROG revealed very large third-order dispersion (TOD). This extent of phase distortion prevents the efficient temporal compression of laser pulses to their theoretical limit. Furthermore, TOD cannot be quantified when using a conventional method of obtaining the laser pulse duration, e.g., when using an autocorrelator. Finally, DMDI showed the effectiveness of detecting pulse-to-pulse intensity fluctuations on timescales relevant to three-photon imaging, which were otherwise not captured using conventional instruments and statistics. Conclusions: The distortion of individual laser pulses caused by TOD poses significant challenges to three-photon imaging by preventing effective compression of laser pulses and decreasing the efficiency of nonlinear excitation. Moreover, an acceptably low pulse-to-pulse amplitude variability should not be assumed. Particularly for low repetition rate laser sources used in three-photon microscopy, pulse-to-pulse variability also degrades image quality. If three-photon imaging is to become mainstream, our diagnostics may be used by laser manufacturers to improve system design and by end-users to validate the performance of their current and future imaging systems.

An Unexpected Dependence of Cortical Depth in Shaping Neural Responsiveness and Selectivity in Mouse Visual Cortex.

Two-photon imaging studies in mouse primary visual cortex (V1) consistently report that around half of the neurons respond to oriented grating stimuli. However, in cats and primates, nearly all neurons respond to such stimuli. Here we show that mouse V1 responsiveness and selectivity strongly depends on neuronal depth. Moving from superficial layer 2 down to layer 4, the percentage of visually responsive neurons nearly doubled, ultimately reaching levels similar to what is seen in other species. Over this span, the amplitude of neuronal responses also doubled. Moreover, stimulus selectivity was also modulated, not only with depth but also with response amplitude. Specifically, we found that orientation and direction selectivity were greater in stronger responding neurons, but orientation selectivity decreased with depth whereas direction selectivity increased. Importantly, these depth-dependent trends were found not just between layer 2/3 and layer 4 but at different depths within layer 2/3 itself. Thus, neuronal depth is an important factor to consider when pooling neurons for population analyses. Furthermore, the inability to drive the majority of cells in superficial layer 2/3 of mouse V1 with grating stimuli indicates that there may be fundamental differences in the micro-circuitry and role of V1 between rodents and other mammals.

In vivo two-photon imaging of neuronal and brain vascular responses in mice chronically exposed to ethanol.

The effects of ethanol on brain function have been extensively studied using a variety of in vitro and in vivo techniques. For example, electrophysiological studies using brain slices from rodents and non-human primates have demonstrated that acute and chronic exposure to ethanol alters the intrinsic excitability and synaptic signaling of neurons within cortical and sub-cortical areas of the brain. In humans, neuroimaging studies reveal alterations in measures of brain activation and connectivity in subjects with alcohol use disorder. While complementary, these methods are inherently limited due to issues related to either disruption of normal sensory input (in vitro slice studies) or resolution (whole brain imaging). In the present study, we used 2-photon laser scanning microscopy in intact animals to assess the impact of chronic ethanol exposure on sensory-evoked neuronal and vascular responses. Adult male C57BL/6J mice were exposed to four weekly cycles of chronic intermittent ethanol (CIE) exposure, while control mice were exposed to air. After withdrawal (≥72 h), a cranial window was placed over the primary visual cortex (V1), and sensory-evoked responses were monitored using the calcium indicator OGB-1. CIE exposure produced small but significant changes in response amplitude (decrease) and orientation selectivity of V1 neurons (increase). While arteriole diameter did not differ between control and CIE mice under baseline conditions, sensory-evoked dilation was enhanced in vessels from CIE-exposed mice as compared to controls. This was accompanied by a reduced latency in response to stimulation. In separate experiments, pial arteriole diameter was measured in the barrel cortex of control and CIE-exposed mice. Baseline diameter of barrel cortex arterioles was similar between control and CIE-exposed mice, but unlike vessels in V1, sensory-evoked dilation of barrel cortex arterioles was similar between the two groups. Together, the results of these studies suggest that chronic exposure to alcohol induces changes in neurovascular coupling that are region-dependent.

Three-photon imaging of synthetic dyes in deep layers of the neocortex.

Multiphoton microscopy has emerged as the primary imaging tool for studying the structural and functional dynamics of neural circuits in brain tissue, which is highly scattering to light. Recently, three-photon microscopy has enabled high-resolution fluorescence imaging of neurons in deeper brain areas that lie beyond the reach of conventional two-photon microscopy, which is typically limited to ~ 450 µm. Three-photon imaging of neuronal calcium signals, through the genetically-encoded calcium indicator GCaMP6, has been used to successfully record neuronal activity in deeper neocortical layers and parts of the hippocampus in rodents. Bulk-loading cells in deeper cortical layers with synthetic calcium indicators could provide an alternative strategy for labelling that obviates dependence on viral tropism and promoter penetration, particularly in non-rodent species. Here we report a strategy for visualized injection of a calcium dye, Oregon Green BAPTA-1 AM (OGB-1 AM), at 500-600 µm below the surface of the mouse visual cortex in vivo. We demonstrate successful OGB-1 AM loading of cells in cortical layers 5-6 and subsequent three-photon imaging of orientation- and direction- selective visual responses from these cells.

Efficacy of retinal spikes in driving cortical responses.

How does a single retinal ganglion cell (RGC) affect the firing of simple cells in the visual cortex? Although much is known of the functional connections between the retina and the lateral geniculate nucleus (LGN) and between LGN and visual cortex, it is hard to infer the effect of disynaptic connections from retina to visual cortex. Most importantly, there is considerable divergence from retina to LGN, so cortical neurons might be influenced by ganglion cells through multiple feedforward pathways. We recorded simultaneously from ganglion cells in the retina and cortical simple cells in the striate cortex with overlapping receptive fields and evaluated disynaptic connections with cross-correlation analysis. In all disynaptically connected pairs, the retinal receptive field center and overlapping cortical subregion always shared the same sign (either both ON or both OFF). Connected pairs were similar in other respects, such as relative position and timing of their receptive fields, and thus obeyed the same rules of connectivity found previously for retinothalamic and thalamocortical connections. We found that a single RGC directly contributed on average to approximately 3% of the activity of its cortical target. The relative timing of pairs of spikes from the retinal cell affected their efficacy in driving the cortical cell. When two retinal spikes were closely spaced (<10 msec), the second spike was several times more likely to drive the cortical target. The relative magnitude of this disynaptic paired spike enhancement was considerably larger than has been found previously for retinogeniculate and geniculocortical connections. The amplified paired spike enhancement from retina to cortex ensures that signal transmission from retina to cortex is particularly effective when the retina fires a series of closely spaced action potentials.

Improved blood velocity measurements with a hybrid image filtering and iterative Radon transform algorithm.

Neural activity leads to hemodynamic changes which can be detected by functional magnetic resonance imaging (fMRI). The determination of blood flow changes in individual vessels is an important aspect of understanding these hemodynamic signals. Blood flow can be calculated from the measurements of vessel diameter and blood velocity. When using line-scan imaging, the movement of blood in the vessel leads to streaks in space-time images, where streak angle is a function of the blood velocity. A variety of methods have been proposed to determine blood velocity from such space-time image sequences. Of these, the Radon transform is relatively easy to implement and has fast data processing. However, the precision of the velocity measurements is dependent on the number of Radon transforms performed, which creates a trade-off between the processing speed and measurement precision. In addition, factors like image contrast, imaging depth, image acquisition speed, and movement artifacts especially in large mammals, can potentially lead to data acquisition that results in erroneous velocity measurements. Here we show that pre-processing the data with a Sobel filter and iterative application of Radon transforms address these issues and provide more accurate blood velocity measurements. Improved signal quality of the image as a result of Sobel filtering increases the accuracy and the iterative Radon transform offers both increased precision and an order of magnitude faster implementation of velocity measurements. This algorithm does not use a priori knowledge of angle information and therefore is sensitive to sudden changes in blood flow. It can be applied on any set of space-time images with red blood cell (RBC) streaks, commonly acquired through line-scan imaging or reconstructed from full-frame, time-lapse images of the vasculature.