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

Dual-exposure temporal laser speckle imaging for simultaneously accessing microvascular blood perfusion and angiography

Laser speckle contrast imaging (LSCI) has gained significant attention in the biomedical field for its ability to map the spatio-temporal dynamics of blood perfusion in vivo. However, LSCI faces difficulties in accurately resolving blood perfusion in microvessels. Although the transmissive detecting geometry can improve the spatial resolution of tissue imaging, ballistic photons directly transmitting forward through tissue without scattering will cause misestimating in the flow speed by LSCI because of the lack of a quantitative theoretical model of transmissvie LSCI. Here, we develop a model of temporal LSCI which accounts for the effect of nonscattered light on estimating decorrelation time. Based on this model, we further propose a dual-exposure temporal laser speckle imaging method (dEtLSCI) to correct the overestimation of background speed when performing traditional transmissive LSCI, and reconstruct microvascular angiography using the scattered component extracted from total transmitted light. Experimental results demonstrated that our new method opens an opportunity for LSCI to simultaneously resolve the blood vessels morphology and blood flow speed at microvascular level in various contexts, ranging from the drug-induced vascular response to angiogenesis and the blood perfusion monitoring during tumor growth.

Tracking the dynamics of thrombus formation in a blood vessel-on-chip with visible-light optical coherence tomography

Thrombus formation is a physiological response to damage in a blood vessel that relies on a complex interplay of platelets, coagulation factors, immune cells, and the vessel wall. The dynamics of thrombus formation are essential for a deeper understanding of many disease processes, like bleeding, wound healing, and thrombosis. However, monitoring thrombus formation is challenging due to the limited imaging options available to analyze flowing blood. In this work, we use a visible-light optical coherence tomography (vis-OCT) system to monitor the dynamic process of the formation of thrombi in a microfluidic blood vessel-on-chip (VoC) device. Inside the VoC, thrombi form in a channel lined with a monolayer of endothelial cells and perfused by human whole blood. We show that the correlation of the vis-OCT signal can be utilized as a marker for thrombus formation. By thresholding the correlation during thrombus formation, we track and quantify the growth of the thrombi over time. We validate our results with fluorescence microscopic imaging of fibrin and platelet markers at the end of the blood perfusion assay. In conclusion, we demonstrate that the correlation of the vis-OCT signal can be used to visualize both the spatial and temporal behavior of the thrombus formation in flowing human whole blood.

Mapping optical scattering properties to physical particle information in singly and multiply scattering samples

Optical coherence tomography (OCT) leverages light scattering by biological tissues as endogenous contrast to form structural images. Light scattering behavior is dictated by the optical properties of the tissue, which depend on microstructural details at the cellular or sub-cellular level. Methods to measure these properties from OCT intensity data have been explored in the context of a number of biomedical applications seeking to access this sub-resolution tissue microstructure and thereby increase the diagnostic impact of OCT. Most commonly, the optical attenuation coefficient, an analogue of the scattering coefficient, has been used as a surrogate metric linking OCT intensity to subcellular particle characteristics. To record attenuation coefficient data that is accurately representative of the underlying physical properties of a given sample, it is necessary to account for the impact of the OCT imaging system itself on the distribution of light intensity in the sample, including the numerical aperture (NA) of the system and the location of the focal plane with respect to the sample surface, as well as the potential contribution of multiple scattering to the reconstructed intensity signal. Although these considerations complicate attenuation coefficient measurement and interpretation, a suitably calibrated system may potentiate a powerful strategy for gaining additional information about the scattering behavior and microstructure of samples. In this work, we experimentally show that altering the OCT system geometry minimally impacts measured attenuation coefficients in samples presumed to be singly scattering, but changes these measurements in more highly scattering samples. Using both depth-resolved attenuation coefficient data and layer-resolved backscattering coefficients, we demonstrate the retrieval of scattering particle diameter and concentration in tissue-mimicking phantoms, and the impact of presumed multiple scattering on these calculations. We further extend our approach to characterize a murine brain tissue sample and highlight a tumor-bearing region based on increased scattering particle density. Through these methods, we not only enhance conventional OCT attenuation coefficient analysis by decoupling the independent effects of particle size and concentration, but also discriminate areas of strong multiple scattering through minor changes to system topology to provide a framework for assessing the accuracy of these measurements.

Adaptive dynamic analysis-based optical coherence tomography angiography for blood vessel projection artifact suppression

Optical coherence tomography angiography (OCTA) for blood vessel 3-D structure imaging suffers from blood vessel projection artifacts/tail artifacts when using a long decorrelation time (e.g., repeat B-scan acquisition in regular OCTA) or loss of micro vessel signal when using a short decorrelation time. In this work, we developed an adaptive first-order field autocorrelation function (g₁) analysis-based technique to suppress the projection artifacts under macro vessels while enhancing the dynamic signal of micro vessels. The proposed method is based on the differences of the decorrelation rate and the phase variations of g₁ between the vessel voxels and the artifacts regions. A short or long decorrelation time was applied to obtain the dynamic index of the projection artifacts region or the blood vessel region, respectively. Compared to the slab subtraction-based post-image processing-based techniques, the proposed approach addresses this problem on a physical basis and shows the ability to suppress the projection artifacts while enhancing the detection of the micro vessels.

Using beam-offset optical coherence tomography to reconstruct backscattered photon profiles in scattering media

Raster scanning imaging technologies capture least scattered photons (LSPs) and reject multiple scattered photons (MSPs) in backscattered photons to image the underlying structures of a scattering medium. However, MSPs can still squeeze into the images, resulting in limited imaging depth, degraded contrast, and significantly reduced lateral resolution. Great efforts have been made to understand how MSPs affect imaging performance through modeling, but the techniques for visualizing the backscattered photon profile (BSPP) in scattering media during imaging are unavailable. Here, a method of reconstructing BSPP is demonstrated using beam-offset optical coherence tomography (OCT), in which OCT images are acquired at offset positions from the illumination beam. The separation of LSPs and MSPs based on the BSPP enables quantification of imaging depth, contrast, and lateral resolution, as well as access to the depth-resolved modulated transfer function (MTF). This approach presents great opportunities for better retrieving tissue optical properties, correctly interpreting images, or directly using MTF as the feedback for adaptive optical imaging.

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.

Signal-carrying speckle in optical coherence tomography: a methodological review on biomedical applications

SIGNIFICANCE

Speckle has historically been considered a source of noise in coherent light imaging. However, a number of works in optical coherence tomography (OCT) imaging have shown that speckle patterns may contain relevant information regarding subresolution and structural properties of the tissues from which it is originated.

AIM

The objective of this work is to provide a comprehensive overview of the methods developed for retrieving speckle information in biomedical OCT applications.

APPROACH

PubMed and Scopus databases were used to perform a systematic review on studies published until December 9, 2021. From 146 screened studies, 40 were eligible for this review.

RESULTS

The studies were clustered according to the nature of their analysis, namely static or dynamic, and all features were described and analyzed. The results show that features retrieved from speckle can be used successfully in different applications, such as classification and segmentation. However, the results also show that speckle analysis is highly application-dependant, and the best approach varies between applications.

CONCLUSIONS

Several of the reviewed analyses were only performed in a theoretical context or using phantoms, showing that signal-carrying speckle analysis in OCT imaging is still in its early stage, and further work is needed to validate its applicability and reproducibility in a clinical context.

Pulsatile tissue deformation dynamics of the murine retina and choroid mapped by 4D optical coherence tomography

Irregular ocular pulsatility and altered mechanical tissue properties are associated with some of the most sight-threatening eye diseases. Here we present 4D optical coherence tomography (OCT) for the quantitative assessment and depth-resolved mapping of pulsatile dynamics in the murine retina and choroid. Through a pixel-wise analysis of phase changes of the complex OCT signal, we reveal spatiotemporal displacement characteristics across repeated frame acquisitions. We demonstrate in vivo fundus elastography (FUEL) imaging in wildtype mouse retinas and in a mouse model of retinal neovascularization and uncover subtle structural deformations related to ocular pulsation. Our data in mouse eyes hold promise for a powerful retinal elastography technique that may enable a new paradigm of OCT-based measurements and image contrast.

Layer-based, depth-resolved computation of attenuation coefficients and backscattering fractions in tissue using optical coherence tomography

Structural optical coherence tomography (OCT) images of tissue stand to benefit from greater functionalization and quantitative interpretation. The OCT attenuation coefficient µ, an analogue of the imaged sample's scattering coefficient, offers potential functional contrast based on the relationship of µ to sub-resolution physical properties of the sample. Attenuation coefficients are computed either by fitting a representative µ over several depth-wise pixels of a sample's intensity decay, or by using previously-developed depth-resolved attenuation algorithms by Girard et al. [Invest. Ophthalmol. Vis. Sci.52, 7738 (2011). 10.1167/iovs.10-6925] and Vermeer et al. [Biomed. Opt. Express5, 322 (2014). 10.1364/BOE.5.000322]. However, the former method sacrifices axial information in the tomogram, while the latter relies on the stringent assumption that the sample's backscattering fraction, another optical property, does not vary along depth. This assumption may be violated by layered tissues commonly observed in clinical imaging applications. Our approach preserves the full depth resolution of the attenuation map but removes its dependence on backscattering fraction by performing signal analysis inside individual discrete layers over which the scattering properties (e.g., attenuation and backscattering fraction) vary minimally. Although this approach necessitates the detection of these layers, it removes the constant-backscattering-fraction assumption that has constrained quantitative attenuation coefficient analysis in the past, and additionally yields a layer-resolved backscattering fraction, providing complementary scattering information to the attenuation coefficient. We validate our approach using automated layer detection in layered phantoms, for which the measured optical properties were in good agreement with theoretical values calculated with Mie theory, and show preliminary results in tissue alongside corresponding histological analysis. Together, accurate backscattering fraction and attenuation coefficient measurements enable the estimation of both particle density and size, which is not possible from attenuation measurements alone. We hope that this improvement to depth-resolved attenuation coefficient measurement, augmented by a layer-resolved backscattering fraction, will increase the diagnostic power of quantitative OCT imaging.

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.