Velocity gradients in spatially resolved laser Doppler flowmetry and dynamic light scattering with confocal and coherence gating

Dual beam optical coherence tomography angiography for decoupling axial velocity gradient

Axial velocity gradient (AVG) in the optical coherence tomography angiography (OCTA) signal affects measurement accuracy when the flow is not perpendicular to the scanning beam. We developed a dual beam OCTA method to decouple the contribution of AVG from the decorrelation signal. Decoupling is first verified by phantom experiments which reduces measurement uncertainty from 1.5 to 0.7% (standard deviation). We also tested the method in human skin in vivo and the results indicate that the contribution of AVG to decorrelation signal is reduced.

Experimental and theoretical model of microvascular network remodeling and blood flow redistribution following minimally invasive microvessel laser ablation

Overly dense microvascular networks are treated by selective reduction of vascular elements. Inappropriate manipulation of microvessels could result in loss of host tissue function or a worsening of the clinical problem. Here, experimental, and computational models were developed to induce blood flow changes via selective artery and vein laser ablation and study the compensatory collateral flow redistribution and vessel diameter remodeling. The microvasculature was imaged non-invasively by bright-field and multi-photon laser microscopy, and optical coherence tomography pre-ablation and up to 30 days post-ablation. A theoretical model of network remodeling was developed to compute blood flow and intravascular pressure and identify vessels most susceptible to changes in flow direction. The skin microvascular remodeling patterns were consistent among the five specimens studied. Significant remodeling occurred at various time points, beginning as early as days 1-3 and continuing beyond day 20. The remodeling patterns included collateral development, venous and arterial reopening, and both outward and inward remodeling, with variations in the time frames for each mouse. In a representative specimen, immediately post-ablation, the average artery and vein diameters increased by 14% and 23%, respectively. At day 20 post-ablation, the maximum increases in arterial and venous diameters were 2.5× and 3.3×, respectively. By day 30, the average artery diameter remained 11% increased whereas the vein diameters returned to near pre-ablation values. Some arteries regenerated across the ablation sites via endothelial cell migration, while veins either reconnected or rerouted flow around the ablation site, likely depending on local pressure driving forces. In the intact network, the theoretical model predicts that the vessels that act as collaterals after flow disruption are those most sensitive to distant changes in pressure. The model results correlate with the post-ablation microvascular remodeling patterns.

Experimental and Theoretical Model of Single Vessel Minimally Invasive Micro-Laser Ablation: Inducing Microvascular Network Remodeling and Blood Flow Redistribution Without Compromising Host Tissue Function

Overly dense microvascular networks are treated by selective reduction of vascular elements. Inappropriate manipulation of microvessels could result in loss of host tissue function or a worsening of the clinical problem. Here, experimental, and computational models were developed to induce blood flow changes via selective artery and vein laser ablation and study the compensatory collateral flow redistribution and vessel diameter remodeling. The microvasculature was imaged non-invasively by bright-field and multi-photon laser microscopy, and Optical Coherence Tomography pre-ablation and up to 30 days post-ablation. A theoretical model of network remodeling was developed to compute blood flow and intravascular pressure and identify vessels most susceptible to changes in flow direction. The skin microvascular remodeling patterns were consistent among the five specimens studied. Significant remodeling occurred at various time points, beginning as early as days 1-3 and continuing beyond day 20. The remodeling patterns included collateral development, venous and arterial reopening, and both outward and inward remodeling, with variations in the time frames for each mouse. In a representative specimen, immediately post-ablation, the average artery and vein diameters increased by 14% and 23%, respectively. At day 20 post-ablation, the maximum increases in arterial and venous diameters were 2.5x and 3.3x, respectively. By day 30, the average artery diameter remained 11% increased whereas the vein diameters returned to near pre-ablation values. Some arteries regenerated across the ablation sites via endothelial cell migration, while veins either reconnected or rerouted flow around the ablation site, likely depending on local pressure driving forces. In the intact network, the theoretical model predicts that the vessels that act as collaterals after flow disruption are those most sensitive to distant changes in pressure. The model results match the post-ablation microvascular remodeling patterns.

Localized transverse flow measurement with dynamic light scattering line-scan OCT

A novel decorrelation-based approach for measuring localized transverse flow velocity using line-scan (LS) optical coherence tomography (OCT) is proposed. The new approach allows for separation of the flow velocity component along the line-illumination direction of the imaging beam from other orthogonal velocity components, from particle diffusion motion, and from noise-induced distortion in the OCT signal's temporal autocorrelation. The new method was verified by imaging flow in a glass capillary and a microfluidic device and mapping the spatial distribution of the flow velocity within the beam's illumination plane. This method can be extended in the future to map the three-dimensional flow velocity fields for both ex-vivo and in-vivo applications.

Sub-diffusion flow velocimetry with number fluctuation optical coherence tomography

We have implemented number fluctuation dynamic light scattering optical coherence tomography (OCT) for measuring extremely slow, sub-diffusion flows of dilute particle suspensions using the second-order autocovariance function. Our method has a lower minimum measurable velocity than conventional correlation-based OCT or phase-resolved Doppler OCT, as the velocity estimation is not affected by the particle diffusion. Similar to non-dilute correlation-based OCT, our technique works for any Doppler angle. With our analysis we can quantitatively determine the concentration of particles under flow. Finally, we demonstrate 2D sub-diffusion flow imaging with a scanning OCT system at high rate by performing number fluctuation correlation analysis on subsequent B-scans.

Using the dynamic forward scattering signal for optical coherence tomography based blood flow quantification

To our knowledge, all existing optical coherence tomography approaches for quantifying blood flow, whether Doppler-based or decorrelation-based, analyze light that is back-scattered by moving red blood cells (RBCs). This work investigates the potential advantages of basing these measurements on light that is forward-scattered by RBCs, i.e., by looking at the signals back-scattered from below the vessel. We show experimentally that flowmetry based on forward-scattering is insensitive to vessel orientation for vessels that are approximately orthogonal to the imaging beam. We further provide proof-of-principle demonstrations of dynamic forward-scattering (DFS) flowmetry in human retinal and choroidal vessels.

Functional Ultrasound Speckle Decorrelation-Based Velocimetry of the Brain

A high-speed, contrast-free, quantitative ultrasound velocimetry (vUS) for blood flow velocity imaging throughout the rodent brain is developed based on the normalized first-order temporal autocorrelation function of the ultrasound field signal. vUS is able to quantify blood flow velocity in both transverse and axial directions, and is validated with numerical simulation, phantom experiments, and in vivo measurements. The functional imaging ability of vUS is demonstrated by monitoring the blood flow velocity changes during whisker stimulation in awake mice. Compared to existing Power-Doppler- and Color-Doppler-based functional ultrasound imaging techniques, vUS shows quantitative accuracy in estimating both axial and transverse flow speeds and resistance to acoustic attenuation and high-frequency noise.

A Neural Network Approach to Quantify Blood Flow from Retinal OCT Intensity Time-Series Measurements

Many diseases of the eye are associated with alterations in the retinal vasculature that are possibly preceded by undetected changes in blood flow. In this work, a robust blood flow quantification framework is presented based on optical coherence tomography (OCT) angiography imaging and deep learning. The analysis used a forward signal model to simulate OCT blood flow data for training of a neural network (NN). The NN was combined with pre- and post-processing steps to create an analysis framework for measuring flow rates from individual blood vessels. The framework’s accuracy was validated using both blood flow phantoms and human subject imaging, and across flow speed, vessel angle, hematocrit levels, and signal-to-noise ratio. The reported flow rate of the calibrated NN framework was measured to be largely independent of vessel angle, hematocrit levels, and measurement signal-to-noise ratio. In vivo retinal flow rate measurements were self-consistent across vascular branch points, and approximately followed a predicted power-law dependence on the vessel diameter. The presented OCT-based NN flow rate estimation framework addresses the need for a robust, deployable, and label-free quantitative retinal blood flow mapping technique.

Real-time co-localized OCT surveillance of laser therapy using motion corrected speckle decorrelation

We present a system capable of real-time delivery and monitoring of laser therapy by imaging with optical coherence tomography (OCT) through a double-clad fiber (DCF). A double-clad fiber coupler is used to inject and collect OCT light into the core of a DCF and inject the therapy light into its larger inner cladding, allowing for both imaging and therapy to be perfectly coregistered. Monitoring of treatment depth is achieved by calculating the speckle intensity decorrelation occurring during tissue coagulation. Furthermore, an analytical noise correction was used on the correlation to extend the maximum monitoring depth. We also present a method for correcting motion-induced decorrelation using a lookup table. Using the value of the noise- and motion-corrected correlation coefficient in a novel approach, our system is capable of identifying the depth of thermal coagulation in real time and automatically shut the therapy laser off when the targeted depth is reached. The process is demonstrated ex vivo in rat tongue and abdominal muscles for depths ranging from 500 µm to 1000 µm with induced motion in real time.

Noise and bias in optical coherence tomography intensity signal decorrelation

Functional optical coherence tomography (OCT) imaging based on the decorrelation of the intensity signal has been used extensively in angiography and is finding use in flowmetry and therapy monitoring. In this work, we present a rigorous analysis of the autocorrelation function, introduce the concepts of contrast bias, statistical bias and variability, and identify the optimal definition of the second-order autocorrelation function (ACF) g ⁽²⁾ to improve its estimation from limited data. We benchmark different averaging strategies in reducing statistical bias and variability. We also developed an analytical correction for the noise contributions to the decorrelation of the ACF in OCT that extends the signal-to-noise ratio range in which ACF analysis can be used. We demonstrate the use of all the tools developed in the experimental determination of the lateral speckle size depth dependence in a rotational endoscopic probe with low NA, and we show the ability to more accurately determine the rotational speed of an endoscopic probe to implement NURD detection. We finally present g ⁽²⁾-based angiography of the finger nailbed, demonstrating the improved results from noise correction and the optimal bias mitigation strategies.

Forward multiple scattering dominates speckle decorrelation in whole-blood flowmetry using optical coherence tomography

Quantitative blood flow measurements using optical coherence tomography (OCT) have a wide potential range of medical research and clinical applications. Flowmetry based on the temporal dynamics of the OCT signal may have the ability to measure three-dimensional flow profiles regardless of the flow direction. State-of-the-art models describing the OCT signal temporal statistics are based on dynamic light scattering (DLS), a model which is inherently limited to single scattering regimes. DLS methods continue to be applied to OCT despite the knowledge that red blood cells produce strong forward multiple scattering. Here, we postulate that forward multiple scattering is the primary mechanism causing the rate of speckle-decorrelation derived from data acquired in vivo to deviate from the rate of decorrelation determined in phantom experiments. We also postulate that multiple scattering contributions to decorrelation are only present when the sample exhibits velocity field inhomogeneities larger than the scale of a resolution volume and are thus absent in rigid bulk motion. To test these hypotheses, we performed a systematic study of the effects of forward multiple scattering on OCT signal decorrelation with phantom experiments under physiologically relevant flow conditions and relative bulk motion. Our experimental results confirm that the amount of forward multiple scattering affects the proportionality between lateral flow and decorrelation. We propose that multiply scattered light carries information from different locations in the sample and each location imprints scattering dynamics on the scattered light causing increased decorrelation rates. Our analysis confirms that the detection of forward scattered light inside the vessel lumen causes an increase in the rate of decorrelation which results in an overestimation of blood flow velocities at depths as shallow as 40 µm into whole blood for OCT systems with typical numerical apertures used in retinal imaging.

Optical coherence tomography velocimetry based on decorrelation estimation of phasor pair ratios (DEPPAIR)

Quantitative velocity estimations in optical coherence tomography requires the estimation of the axial and lateral flow components. Optical coherence tomography measures the depth resolved complex field reflected from a sample. While the axial velocity component can be determined from the Doppler shift or phase shift between a pair of consecutive measurements at the same location, the estimation of the lateral component for in vivo applications is still challenging. One approach to determine lateral velocity is multiple simultaneous measurements at different angles. In another approach the lateral component can be retrieved through repeated measurements at (nearly) the same location by an analysis of the decorrelation over time. In this paper we follow the latter approach. We describe a model for the complex field changes between consecutive measurements and use it to predict the uncertainties for amplitude-based, phase-based and complex algorithms. The uncertainty of the flow estimations follows from a statistical analysis and is determined by the number of available measurements and the applied analysis method. The model is verified in phantom measurements and the dynamic range of velocity estimations is investigated. We demonstrate that phase-based and complex (phasor) based lateral flow estimation methods are superior to amplitude-based algorithms.

Impact of velocity gradient in Poiseuille flow on the statistics of coherent radiation scattered by flowing Brownian particles in optical coherence tomography

A closed-form expression is obtained for the temporal correlation function of the scattered radiation detected in optical coherence tomography (OCT), taking into account the flow velocity gradient across the OCT detection volume in the suspension of flowing Brownian particles. The analytical approach we use includes both the laser beam and wavefront curvature radii changing over the depth. Also, we compare our results with a previously obtained theoretical model, partially an empirical approach. Our findings suggest the importance of the flow velocity gradient for accurate measurements of flow velocity vector, particle diffusivity, shear-induced diffusion, and potentially other OCT applications.

Normalized field autocorrelation function-based optical coherence tomography three-dimensional angiography

Optical coherence tomography angiography (OCTA) has been widely used for en face visualization of the microvasculature, but is challenged for real three-dimensional (3-D) topologic imaging due to the "tail" artifacts that appear below large vessels. Further, OCTA is generally incapable of differentiating descending arterioles from ascending venules. We introduce a normalized field autocorrelation function-based OCTA (g1-OCTA), which minimizes the tail artifacts and is capable of distinguishing penetrating arterioles from venules in the 3-D image. g1   (  τ  )   is calculated from repeated optical coherence tomography (OCT) acquisitions for each spatial location. The decay amplitude of g1   (  τ  )   is retrieved to represent the dynamics for each voxel. To account for the small g1   (  τ  )   decay in capillaries where red blood cells are flowing slowly and discontinuously, Intralipid is injected to enhance the OCT signal. We demonstrate that the proposed technique realizes 3-D OCTA with negligible tail projections and the penetrating arteries are readily identified. In addition, compared to regular OCTA, the proposed g1-OCTA largely increased the depth-of-field. This technique provides a more accurate rendering of the vascular 3-D anatomy and has the potential for more quantitative characterization of vascular networks.

Computationally effective 2D and 3D fast phase unwrapping algorithms and their applications to Doppler optical coherence tomography

We propose a simplification for a robust and easy to implement fast phase unwrapping (FPU) algorithm that is used to solve the phase wrapping problem encountered in various fields of optical imaging and metrology. We show that the number of necessary computations using the algorithm can be reduced compared to its original version. FPU can be easily extended from two to three spatial dimensions. We demonstrate the applicability of the two- and three-dimensional FPU algorithm for Doppler optical coherence tomography (DOCT) in numerical simulations, and in the imaging of a flow phantom and blood flow in the human retina in vivo. We introduce an FPU applicability plot for use as a guide in the selection of the most suitable version of the algorithm depending on the phase noise in the acquired data. This plot uses the circular standard deviation of the wrapped phase distribution as a measure of noise and relates it to the root-mean-square error of the recovered, unwrapped phase.

Shear-induced diffusion of red blood cells measured with dynamic light scattering-optical coherence tomography

Quantitative measurements of intravascular microscopic dynamics, such as absolute blood flow velocity, shear stress and the diffusion coefficient of red blood cells (RBCs), are fundamental in understanding the blood flow behavior within the microcirculation, and for understanding why diffuse correlation spectroscopy (DCS) measurements of blood flow are dominantly sensitive to the diffusive motion of RBCs. Dynamic light scattering-optical coherence tomography (DLS-OCT) takes the advantages of using DLS to measure particle flow and diffusion within an OCT resolution-constrained three-dimensional volume, enabling the simultaneous measurements of absolute RBC velocity and diffusion coefficient with high spatial resolution. In this work, we applied DLS-OCT to measure both RBC velocity and the shear-induced diffusion coefficient within penetrating venules of the somatosensory cortex of anesthetized mice. Blood flow laminar profile measurements indicate a blunted laminar flow profile and the degree of blunting decreases with increasing vessel diameter. The measured shear-induced diffusion coefficient was proportional to the flow shear rate with a magnitude of ~0.1 to 0.5 × 10^-6  mm² . These results provide important experimental support for the recent theoretical explanation for why DCS is dominantly sensitive to RBC diffusive motion.

Intravascular optical coherence tomography [Invited]

Shortly after the first demonstration of optical coherence tomography for imaging the microstructure of the human eye, work began on developing systems and catheters suitable for intravascular imaging in order to diagnose and investigate atherosclerosis and potentially to monitor therapy. This review covers the driving considerations of the clinical application and its constraints, the major engineering milestones that enabled the current, high-performance commercial imaging systems, the key studies that laid the groundwork for image interpretation, and the clinical research that traces intravascular optical coherence tomography (OCT) from early human pilot studies to current clinical trials.