Noise Equalization for Ultrafast Plane Wave Microvessel Imaging

Capillary-scale Microvessel Imaging with High-frequency Ultrasound Localization Microscopy in Mouse Brain

Ultrasound localization microscopy is a super-resolution vascular imaging technique which has garnered substantial interest as a tool for small animal neuroimaging, neuroscience research, and the characterization of vascular pathologies. In the pursuit of increasingly high-fidelity reconstructions of microvasculature, there remains several outstanding questions concerning this sub-diffraction imaging technology, including the accurate reconstruction of microvessels approaching the capillary scale and the pragmatic challenges associated with long data acquisition times. In the context of small animal neurovascular imaging, we posit that increasing the ultrasound imaging frequency is a straightforward approach to enable higher concentrations of microbubble contrast agents, thus increasing the likelihood of microvascular/capillary mapping and decreasing the imaging duration. We demonstrate that higher frequency imaging results in improved ULM fidelity and more efficient microbubble localization due to a smaller microbubble point-spread function that is easier to localize, and which can achieve a higher localizable concentration within the same unit volume of tissue. A select example of in vivo capillary-level vascular reconstruction is demonstrated for the highest frequency imaging probe, which has substantial implications for neuroscientists investigating microvascular function in disease states, regulation, and brain development. High frequency ULM yielding a spatial resolution of 7.1μm, as measured by Fourier ring correlation, throughout the entire depth of the brain, highlighting this technology as a highly relevant tool for neuroimaging research.

High-Resolution Power Doppler Using Null Subtraction Imaging

To improve the spatial resolution of power Doppler (PD) imaging, we explored null subtraction imaging (NSI) as an alternative beamforming technique to delay-and-sum (DAS). NSI is a nonlinear beamforming approach that uses three different apodizations on receive and incoherently sums the beamformed envelopes. NSI uses a null in the beam pattern to improve the lateral resolution, which we apply here for improving PD spatial resolution both with and without contrast microbubbles. In this study, we used NSI with three types of singular value decomposition (SVD)-based clutter filters and noise equalization to generate high-resolution PD images. An element sensitivity correction scheme was also proposed as a crucial component of NSI-based PD imaging. First, a microbubble trace experiment was performed to evaluate the resolution improvement of NSI-based PD over traditional DAS-based PD. Then, both contrast-enhanced and contrast free ultrasound PD images were generated from the scan of a rat brain. The cross-sectional profile of the microbubble traces and microvessels were plotted. FWHM was also estimated to provide a quantitative metric. Furthermore, iso-frequency curves were calculated to provide a resolution evaluation metric over the global field of view. Up to six-fold resolution improvement was demonstrated by the FWHM estimate and four-fold resolution improvement was demonstrated by the iso-frequency curve from the NSI-based PD microvessel images compared to microvessel images generated by traditional DAS-based beamforming. A resolvability of [Formula: see text] was measured from the NSI-based PD microvessel image. The computational cost of NSI-based PD was only increased by 40 percent over the DAS-based PD.

Low-rank prior-based Fast-RPCA for clutter filtering and noise suppression in non-contrast ultrasound microvascular imaging

Accurate and real-time separation of blood signal from clutter and noise signals is a critical step in clinical non-contrast ultrasound microvascular imaging. Despite the widespread adoption of singular value decomposition (SVD) and robust principal component analysis (RPCA) for clutter filtering and noise suppression, the SVD's sensitivity to threshold selection, along with the RPCA's limitations in undersampling conditions and heavy computational burden often result in suboptimal performance in complex clinical applications. To address those challenges, this study presents a novel low-rank prior-based fast RPCA (LP-fRPCA) approach to enhance the adaptability and robustness of clutter filtering and noise suppression with reduced computational cost. A low-rank prior constraint is integrated into the non-convex RPCA model to achieve a robust and efficient approximation of clutter subspace, while an accelerated alternating projection iterative algorithm is developed to improve convergence speed and computational efficiency. The performance of the LP-fRPCA method was evaluated against SVD with a tissue/blood threshold (SVD1), SVD with both tissue/blood and blood/noise thresholds (SVD2), and the classical RPCA based on the alternating direction method of multipliers algorithm through phantom and in vivo non-contrast experiments on rabbit kidneys. In the slow flow phantom experiment of 0.2 mm/s, LP-fRPCA achieved an average increase in contrast ratio (CR) of 10.68 dB, 9.37 dB, and 8.66 dB compared to SVD1, SVD2, and RPCA, respectively. In the in vivo rabbit kidney experiment, the power Doppler results demonstrate that the LP-fRPCA method achieved a superior balance in the trade-off between insufficient clutter filtering and excessive suppression of blood flow. Additionally, LP-fRPCA significantly reduced the runtime of RPCA by up to 94-fold. Consequently, the LP-fRPCA method promises to be a potential tool for clinical non-contrast ultrasound microvascular imaging.

Adaptive higher-order singular value decomposition clutter filter for ultrafast Doppler imaging of coronary flow under non-negligible tissue motion

BACKGROUND AND OBJECTIVE

With the development of advanced clutter-filtering techniques by singular value decomposition (SVD) and leveraging favorable acquisition settings such as open-chest imaging by a linear high-frequency probe and plane waves, several studies have shown the feasibility of cardiac flow measurements during the entire cardiac cycle, ranging from coronary flow to myocardial perfusion. When applying these techniques in a routine clinical setting, using transthoracic ultrasound imaging, new challenges emerge. Firstly, a smaller aperture is needed that can fit between ribs. Consequently, diverging waves are employed instead of plane waves to achieve an adequate field of view. Secondly, to ensure imaging at a larger depth, the maximum pulse repetition frequency has to be reduced. Lastly, in comparison to the open-chest scenario, tissue motion induced by the heartbeat is significantly stronger. The latter complicates substantially the distinction between clutter and blood signals.

METHODS

This study investigates a strategy to overcome these challenges by diverging wave imaging with an optimal number of tilt angles, in combination with dedicated clutter-filtering techniques. In particular, a novel, adaptive, higher-order SVD (HOSVD) clutter filter, which utilizes spatial, temporal, and angular information of the received ultrasound signals, is proposed to enhance clutter and blood separation.

RESULTS

When non-negligible tissue motion is present, using fewer tilt angles not only reduces the decorrelation between the received waveforms but also allows for collecting more temporal samples at a given ensemble duration, contributing to improved Doppler performance. The addition of a third angular dimension enables the application of HOSVD, providing greater flexibility in selecting blood separation thresholds from a 3-D tensor. This differs from the conventional threshold selection method in a 2-D spatiotemporal space using SVD. Exhaustive threshold search has shown a significant improvement in Contrast and Contrast-to-Noise ratio for Power Doppler images filtered with HOSVD compared to the SVD-based clutter filter.

CONCLUSION

With the improved settings, the obtained Power Doppler images show the feasibility of measuring coronary flow under the influence of non-negligible tissue motion in both in vitro and ex vivo.

Super-Resolution Ultrasound Reveals Cerebrovascular Impairment in a Mouse Model of Alzheimer's Disease

Increasing evidence has suggested a link between cerebrovascular disease and the cognitive impairment associated with Alzheimer's disease. However, detailed descriptions of microvascular changes across brain regions and how they relate to other more traditional pathology have been lacking. Additionally, the efforts to elucidate the interplay between cerebral microvascular function and Alzheimer's disease progression are complicated by the necessity of probing deep-brain structures since early-stage Alzheimer's disease typically involves hippocampal pathology. The purpose of this study was to examine changes in microvascular dynamics in a mouse model of Alzheimer's disease using cohorts that were age-matched to wild-type controls. Data from both sexes were included in this study. Super-resolution ultrasound localization microscopy revealed microvascular functional and structural features throughout the whole brain depth to visualize and quantify. We found that functional decreases in hippocampal and entorhinal flow velocity preceded structural derangements in regional vascular density. Co-registered histological sectioning confirmed the regionalized perfusion deficits seen on ultrasound imaging, which were co-localized with amyloid beta plaque deposition. In addition to providing global vascular quantifications of deep brain structures with a high local resolution, this technology also permitted velocity-profile analysis of individual vessels and, in some cases, allowed for decoupling of arterial and venous flow contributions. These data suggest that microvascular pathology is an early and pervasive feature of Alzheimer's disease and may represent a novel therapeutic target for this disease.

High-Spatiotemporal-Resolution Ultrasound Flow Imaging to Determine Cerebrovascular Hemodynamics in Alzheimer's Disease Mice Model

Although the relationships of cerebrovascular hemodynamic dysfunction with neurodegenerative diseases remain unclear, many studies have indicated that poor cerebral perfusion accelerates the progression of neurodegenerative diseases, such as Alzheimer's disease (AD). Small animal models are widely used in AD research. However, providing an imaging modality with a high spatiotemporal resolution and sufficiently large field of view to assess cerebrovascular hemodynamics in vivo remains a challenge. The present study proposes a novel technique for high-spatiotemporal-resolution vector micro-Doppler imaging (HVμDI) based on contrast-free ultrafast high frequency ultrasound imaging to visualize the cerebrovascular hemodynamics of the mouse, with a data acquisition time of 0.4 s, a minimal detectable vessel size of 38 µm, and a temporal resolution of 500 Hz. In vivo experiments are conducted on wild-type and AD mice. Cerebrovascular hemodynamics are quantified using the cerebral vascular density, diameter, velocity, tortuosity, cortical flow pulsatility, and instant flow direction variations. Results reveal that AD significantly change the cerebrovascular hemodynamics. HVμDI offers new opportunities for in vivo analysis of cerebrovascular hemodynamics in neurodegenerative pathologies in preclinical animal research.

Longitudinal Awake Imaging of Deep Mouse Brain Microvasculature with Super-resolution Ultrasound Localization Microscopy

Super-resolution ultrasound localization microscopy (ULM) is an emerging imaging modality that resolves capillary-scale microvasculature in deep tissues. However, existing preclinical ULM applications are largely constrained to anesthetized animals, introducing confounding vascular effects such as vasodilation and altered hemodynamics. As such, ULM quantifications (e.g., vessel diameter, density, and flow velocity) may be confounded by the use of anesthesia, undermining the usefulness of ULM in practice. Here we introduce a method to address this limitation and achieve ULM imaging in awake mouse brain. Pupillary monitoring was used to confirm the awake state during ULM imaging. ULM revealed that veins showed a greater degree of vascularity reduction from anesthesia to awake states than did arteries. The reduction was most significant in the midbrain and least significant in the cortex. ULM also revealed a significant reduction in venous blood flow velocity across different brain regions under awake conditions. Serial in vivo imaging of the same animal brain at weekly intervals demonstrated the highly robust longitudinal imaging capability of the proposed technique. This is the first study demonstrating longitudinal ULM imaging in the awake mouse brain, which is essential for many ULM brain applications that require awake and behaving animals.

Multi-layered adaptive neoangiogenesis Intra-Operative quantification (MANIOQ)

Quantification of vascularization volume can provide valuable information for diagnosis and prognosis in vascular pathologies. It can be adapted to inform the surgical management of gliomas, aggressive brain tumors characterized by exuberant sprouting of new blood vessels (neoangiogenesis). Filtered ultrafast Doppler data can provide two main parameters: vascularization index (VI) and fractional moving blood volume (FMBV) that clinically reflect tumor micro vascularization. Current protocols lack robust, automatic, and repeatable filtering methods. We present a filtrating method called Multi-layered Adaptive Neoangiogenesis Intra-Operative Quantification (MANIOQ). First, an adaptive clutter filtering is implemented, based on singular value decomposition (SVD) and hierarchical clustering. Second a method for noise equalization is applied, based on the subtraction of a weighted noise profile. Finally, an in vivo analysis of the periphery of the B-mode hyper signal area allows to measure the vascular infiltration extent of the brain tumors. Ninety ultrasound acquisitions were processed from 23 patients. Compared to reference methods in the literature, MANIOQ provides a more robust tissue filtering, and noise equalization allows for the first time to keep axial and lateral gain compensation (TGC and LGC). MANIOQ opens the way to an intra-operative clinical analysis of gliomas micro vascularization.

Improved Ultrafast Power Doppler Imaging Using United Spatial-Angular Adaptive Scaling Wiener Postfilter

Ultrafast power Doppler imaging (uPDI) using high-frame-rate plane-wave transmission is a new microvascular imaging modality that offers high Doppler sensitivity. However, due to the unfocused transmission of plane waves, the echo signal is subject to interference from noise and clutter, resulting in a low signal-to-noise ratio (SNR) and poor image quality. Adaptive beamforming techniques are effective in suppressing noise and clutter for improved image quality. In this study, an adaptive beamformer based on a united spatial-angular adaptive scaling Wiener (uSA-ASW) postfilter is proposed to improve the resolution and contrast of uPDI. In the proposed method, the signal power and noise power of the Wiener postfilter are estimated by uniting spatial and angular signals, and a united generalized coherence factor (uGCF) is introduced to dynamically adjust the noise power estimation and enhance the robustness of the method. Simulation and in vivo data were used to verify the effectiveness of the proposed method. The results show that the uSA-ASW can achieve higher resolution and significant improvements in image contrast and background noise suppression compared with conventional delay-and-sum (DAS), coherence factor (CF), spatial-angular CF (SACF), and adaptive scaling Wiener (ASW) postfilter methods. In the simulations, uSA-ASW improves contrast-to-noise ratio (CNR) by 34.7 dB (117.3%) compared with DAS, while reducing background noise power (BNP) by 52 dB (221.4%). The uSA-ASW method provides full-width at half-maximum (FWHM) reductions of [Formula: see text] (59.5%) and [Formula: see text] (56.9%), CNR improvements of 25.6 dB (199.9%) and 42 dB (253%), and BNP reductions of 46.1 dB (319.3%) and 12.9 dB (289.1%) over DAS in the experiments of contrast-free human neonatal brain and contrast-free human liver, respectively. In the contrast-free experiments, uSA-ASW effectively balances the performance of noise and clutter suppression and enhanced microvascular visualization. Overall, the proposed method has the potential to become a reliable microvascular imaging technique for aiding in more accurate diagnosis and detection of vascular-related diseases in clinical contexts.

Localization Free Super-Resolution Microbubble Velocimetry Using a Long Short-Term Memory Neural Network

Ultrasound localization microscopy is a super-resolution imaging technique that exploits the unique characteristics of contrast microbubbles to side-step the fundamental trade-off between imaging resolution and penetration depth. However, the conventional reconstruction technique is confined to low microbubble concentrations to avoid localization and tracking errors. Several research groups have introduced sparsity- and deep learning-based approaches to overcome this constraint to extract useful vascular structural information from overlapping microbubble signals, but these solutions have not been demonstrated to produce blood flow velocity maps of the microcirculation. Here, we introduce Deep-SMV, a localization free super-resolution microbubble velocimetry technique, based on a long short-term memory neural network, that provides high imaging speed and robustness to high microbubble concentrations, and directly outputs blood velocity measurements at a super-resolution. Deep-SMV is trained efficiently using microbubble flow simulation on real in vivo vascular data and demonstrates real-time velocity map reconstruction suitable for functional vascular imaging and pulsatility mapping at super-resolution. The technique is successfully applied to a wide variety of imaging scenarios, include flow channel phantoms, chicken embryo chorioallantoic membranes, and mouse brain imaging. An implementation of Deep-SMV is openly available at https://github.com/chenxiptz/SR_microvessel_velocimetry, with two pre-trained models available at https://doi.org/10.7910/DVN/SECUFD.

Contrast-free ultrasound imaging for blood flow assessment of the lower limb in patients with peripheral arterial disease: a feasibility study

While being a relatively prevalent condition particularly among aging patients, peripheral arterial disease (PAD) of lower extremities commonly goes undetected or misdiagnosed due to its symptoms being nonspecific. Additionally, progression of PAD in the absence of timely intervention can lead to dire consequences. Therefore, development of non-invasive and affordable diagnostic approaches can be highly beneficial in detection and treatment planning for PAD patients. In this study, we present a contrast-free ultrasound-based quantitative blood flow imaging technique for PAD diagnosis. The method involves monitoring the variations of blood flow in the calf muscle in response to thigh-pressure-cuff-induced occlusion. Four quantitative metrics are introduced for analysis of these variations. These metrics include post-occlusion to baseline flow intensity variation (PBFIV), total response region (TRR), Lag0 response region (L0RR), and Lag4 (and more) response region (L4 + RR). We examine the feasibility of this method through an in vivo study consisting of 14 PAD patients with abnormal ankle-brachial index (ABI) and 8 healthy volunteers. Ultrasound data acquired from 13 legs in the patient group and 13 legs in the healthy group are analyzed. Out of the four utilized metrics, three exhibited significantly different distributions between the two groups (p-value < 0.05). More specifically, p-values of 0.0015 for PBFIV, 0.0183 for TRR, and 0.0048 for L0RR were obtained. The results of this feasibility study indicate the diagnostic potential of the proposed method for the detection of PAD.

High-Quality Ultrafast Power Doppler Imaging Based on Spatial Angular Coherence Factor

The morphological and hemodynamic changes of microvessels are demonstrated to be related to the diseased conditions in tissues. Ultrafast power Doppler imaging (uPDI) is a novel modality with a significantly increased Doppler sensitivity, benefiting from the ultrahigh frame rate plane-wave imaging (PWI) and advanced clutter filtering. However, unfocused plane-wave transmission often leads to a low imaging quality, which degrades the subsequent microvascular visualization in power Doppler imaging. Coherence factor (CF)-based adaptive beamformers have been widely studied in conventional B-mode imaging. In this study, we propose a spatial and angular coherence factor (SACF) beamformer for improved uPDI (SACF-uPDI) by calculating the spatial CF across apertures and the angular CF across transmit angles, respectively. To identify the superiority of SACF-uPDI, simulations, in vivo contrast-enhanced rat kidney, and in vivo contrast-free human neonatal brain studies were conducted. Results demonstrate that SACF-uPDI can effectively enhance contrast and resolution and suppress background noise simultaneously, compared with conventional uPDI methods based on delay-and-sum (DAS) (DAS-uPDI) and CF (CF-uPDI). In the simulations, SACF-uPDI can improve the lateral and axial resolutions compared with those of DAS-uPDI, from 176 to [Formula: see text] of lateral resolution, and from 111 to [Formula: see text] of axial resolution. In the in vivo contrast-enhanced experiments, SACF achieves 15.14- and 5.6-dB higher contrast-to-noise ratio (CNR), 15.25- and 3.68-dB lower noise power, and 240- and 15- [Formula: see text] narrower full-width at half-maximum (FWHM) than DAS-uPDI and CF-uPDI, respectively. In the in vivo contrast-free experiments, SACF achieves 6.11- and 1.09-dB higher CNR, 11.93- and 4.01-dB lower noise power, and 528- and 160- [Formula: see text] narrower FWHM than DAS-uPDI and CF-uPDI, respectively. In conclusion, the proposed SACF-uPDI method can efficiently improve the microvascular imaging quality and has the potential to facilitate clinical applications.

Region-based SVD processing of high-frequency ultrafast ultrasound to visualize cutaneous vascular networks

Observing alterations in cutaneous vasculature in response to any disease or pathology is considered as a potential diagnostic marker in the progression and cure of a disease. To observe skin morphologies and tissue conditions, high-frequency ultrasound (HFUS) has been used in dermatology, although its ability to selectively visualize micro-vessels is limited due to insufficient Doppler sensitivity to peripheral slow-speed blood flow. In recent studies, this issue has been improved by increasing the sensitivity of Doppler imaging to slow flow, leveraging advanced cutter filtering approaches based on singular value decomposition (SVD) techniques that aid to effectively extract flow signals overlapped with tissue echo signals. Nevertheless, in skin imaging, variations in flow speed, diameter, and depth of the blood vessels at different skin layers can make clutter filtering challenging because these variations are problematic in selecting the optimal cut-off value for the SVD filtering. In this study, we aimed to devise a novel region-based SVD filtering approach for ultrafast HFUS data to visualize cutaneous vascular networks. The proposed method divides the acquired high-framerate data into two regions based on B-mode cutaneous morphological identification (dermis layer and subcutaneous tissue). Singular value decomposition processing was performed on each region to effectively extract the desired flow signal, and the processed regions were merged to generate a single power Doppler image, thereby highlighting the appearance of a complete cutaneous vascular network. Finally, top-hat transform was applied to the power Doppler image to further suppress the background noises and enhances the visibility of the micro-vessels. Experimental observations of the human cutaneous circulation showed that the image quality (contrast-to-noise ratio) through the region-based SVD filtering was measured to be 4.1 dB (before top-hat filtering) and 5.2 dB (after top-hat filtering), which were improved from 3.4 dB and 4.0 dB obtained using the global SVD approach with and without top-hat filtering, respectively. We envisioned that this approach would provide diverse applications in the diagnosis of cutaneous disorders.

Resistivity index mapping in Kidney based on ultrasensitive Pulsed-Wave Doppler and automatic spectrogram envelope detection

In recent years, ultrasensitive Pulsed-Wave Doppler (uPWD) ultrasound (US) has emerged as an alternative imaging approach for microcirculation imaging and as a complementary tool to other imaging modalities, such as positron emission tomography (PET). uPWD is based on the acquisition of a large set of highly spatiotemporally coherent frames, which allows high-quality images of a wide field of view to be obtained. In addition, these acquired frames allow calculation of the resistivity index (RI) of the pulsatile flow detected over the entire field of view, which is of great interest to clinicians, for example, in monitoring the transplanted kidney course. This work aims to develop and evaluate a method to automatically obtain an RI map of the kidney based on the uPWD approach. The effect of time gain compensation (TGC) on the visualization of vascularization and aliasing on the blood flow frequency response, was also assessed. A pilot study conducted in patients referred for renal transplant Doppler examination showed that the proposed method provided relative errors of about 15% for RI measurements with respect to conventional pulsed-wave (PW) Doppler.

Adaptive noise reduction for power Doppler imaging using SVD filtering in the channel domain and coherence weighting of pixels

Objective. Ultrafast power Doppler (UPD) is an ultrasound method that can image blood flow at several thousands of frames per second. In particular, the high number of data provided by UPD enables the use of singular value decomposition (SVD) as a clutter filter for suppressing tissue signal. Notably, is has been demonstrated in various applications that SVD filtering increases significantly the sensitivity of UPD to microvascular flows. However, UPD is subjected to significant depth-dependent electronic noise and an optimal denoising approach is still being sought.Approach. In this study, we propose a new denoising method for UPD imaging: the Coherence Factor Mask (CFM). This filter is first based on filtering the ultrasound time-delayed data using SVD in the channel domain to remove clutter signal. Then, a spatiotemporal coherence mask that exploits coherence information between channels for identifying noisy pixels is computed. The mask is finally applied to beamformed images to decrease electronic noise before forming the power Doppler image. We describe theoretically how to filter channel data using a single SVD. Then, we evaluate the efficiency of the CFM filter for denoisingin vitroandin vivoimages and compare its performances with standard UPD and with three existing denoising approaches.Main results. The CFM filter gives gains in signal-to-noise ratio and contrast-to-noise ratio of up to 22 dB and 20 dB, respectively, compared to standard UPD and globally outperforms existing methods for reducing electronic noise. Furthermore, the CFM filter has the advantage over existing approaches of being adaptive and highly efficient while not requiring a cut-off for discriminating noise and blood signals nor for determining an optimal coherence lag.Significance. The CFM filter has the potential to help establish UPD as a powerful modality for imaging microvascular flows.

Quantitative tissue perfusion imaging using nonlinear ultrasound localization microscopy

Ultrasound localization microscopy (ULM) is a recent advancement in ultrasound imaging that uses microbubble contrast agents to yield vascular images that break the classical diffraction limit on spatial resolution. Current approaches cannot image blood flow at the tissue perfusion level since they rely solely on differences in velocity to separate tissue and microbubble signals; lower velocity microbubble echoes are removed during high pass wall filtering. To visualize blood flow in the entire vascular tree, we have developed nonlinear ULM, which combines nonlinear pulsing sequences with plane-wave imaging to segment microbubble signals independent of their velocity. Bubble localization and inter-frame tracking produces super-resolved images and, with parameters derived from the bubble tracks, a rich quantitative feature set that can describe the relative quality of microcirculatory flow. Using the rat spinal cord as a model system, we showed that nonlinear ULM better resolves some smaller branching vasculature compared to conventional ULM. Following contusion injury, both gold-standard histological techniques and nonlinear ULM depicted reduced in-plane vessel length between the penumbra and contralateral gray matter (-16.7% vs. -20.5%, respectively). Here, we demonstrate that nonlinear ULM uniquely enables investigation and potential quantification of tissue perfusion, arguably the most important component of blood flow.

Three-Dimensional Ultrasound Localization Microscopy with Bipartite Graph-Based Microbubble Pairing and Kalman-Filtering-Based Tracking on a 256-Channel Verasonics Ultrasound System with a 32 × 32 Matrix Array

Three-dimensional (3D) ultrasound localization microscopy (ULM) using a 2-D matrix probe and microbubbles (MBs) has been recently proposed to visualize microvasculature beyond the ultrasound diffraction limit in three spatial dimensions. However, 3D ULM suffers from several limitations: (1) high system complexity due to numerous channel counts, (2) complex MB flow dynamics in 3D, and (3) extremely long acquisition time. To reduce the system complexity while maintaining high image quality, we used a sub-aperture process to reduce received channel counts. To address the second issue, a 3D bipartite graph-based method with Kalman filtering-based tracking was used in this study for MB tracking. An MB separation approach was incorporated to separate high concentration MB data into multiple, sparser MB datasets, allowing better MB localization and tracking for a limited acquisition time. The proposed method was first validated in a flow channel phantom, showing improved spatial resolutions compared with the contrasted enhanced power Doppler image. Then the proposed method was evaluated with an in vivo chicken embryo brain dataset. Results showed that the reconstructed 3D super-resolution image achieved a spatial resolution of around 52 μm (smaller than the wavelength of around 200 μm). Microvessels that cannot be resolved clearly using localization only, can be well identified with the tailored 3D pairing and tracking algorithms. To sum up, the feasibility of the 3D ULM is shown, indicating the great possibility in clinical applications.

Adaptive and Robust Vessel Quantification in Contrast-Free Ultrafast Ultrasound Microvessel Imaging

The morphological features of vasculature in diseased tissue differ significantly from those in normal tissue. Therefore, vasculature quantification is crucial for disease diagnosis and staging. Ultrasound microvessel imaging (UMI) with ultrafast ultrasound acquisitions has been determined to have potential in clinical applications given its superior sensitivity in blood flow detection. However, the presence of spatial-dependent noise caused by a low imaging signal-to-noise ratio and incoherent clutter artifacts caused by moving hyperechoic scatterers degrades the performance of UMI and the reliability of vascular quantification. To tackle these issues, we proposed an improved UMI technique along with an adaptive vessel segmentation workflow for robust vessel identification and vascular feature quantification. A previously proposed sub-aperture cross-correlation technique and a normalized cross-correlation technique were applied to equalize the spatially dependent noise level and suppress the incoherent clutter artifact. A square operator and non-local means filter were then used to better separate the blood flow signal from residual background noise. On the de-noised ultrasound microvessel image, an automatic and adaptive vessel segmentation method was developed based on the different spatial patterns of blood flow signal and background noise. The proposed workflow was applied to a CIRS phantom, to a Doppler flow phantom and to an inflammatory bowel, kidney and liver, to validate its feasibility. Results revealed that automatic adaptive, and robust vessel identification performance can be achieved using the proposed method without the subjectivity caused by radiologists/operators.

Development of an Imaging Framework for Visualization of Cutaneous Micro-Vasculature by using High Frequency Ultrafast Ultrasound Imaging

Visualization of cutaneous micro-vasculatures is a determined approach in the diagnosis of skin vascular disorders. Clinically, high frequency ultrasound (HFUS) modalities have been used for cutaneous morphological and structural imaging, but visualization of micro-vessels has always been remained a daunting task. These tiny structures might be visualized by devising a highly sensitive Doppler technique for HFUS systems. In this study, we proposed an imaging framework using HFUS (30 MHz) ultrafast Doppler imaging along with SVD clutter filtering that is proficient in detection of micro-scale circulation. The performance of the devised framework was examined on a 200-micron flow phantom made of poly-vinyl alcohol under four different flow rates (56 - 18 ul/min) and visualized the micro-structure with averaged detected diameter of 93 - 170 µm. The results indicated that the devised framework has sufficient sensitivity and resolvability to visualize the micro-vasculatures in dermis layer of skin (depth ≤ 4 mm). Clinical Relevance - This study brings an insight to visualize in-vivo cutaneous micro-vasculatures with ultrafast Doppler imaging in clinical applications for better assessment of cutaneous disorders.

Improved Ultrafast Power Doppler Imaging by Using Spatiotemporal Non-Local Means Filtering

The change of microvasculature is associated with the occurrence and development of many diseases. Ultrafast power Doppler imaging (uPDI) is an emerging technology for the visualization of microvessels due to the development of ultrafast plane wave (PW) imaging and advanced clutter filters. However, the low signal-to-noise ratio (SNR) caused by unfocused transmit of PW imaging deteriorates the subsequent imaging of microvasculature. Nonlocal means (NLM) filtering has been demonstrated to be effective in the denoising of both natural and medical images, including ultrasound power Doppler images. However, the feasibility and performance of applying an NLM filter on the ultrasound radio frequency (RF) data have not been investigated so far. In this study, we propose to apply an NLM filter on the spatiotemporal domain of clutter filtered blood flow RF data (St-NLM) to improve the quality of uPDI. Experiments were conducted to compare the proposed method with three different methods (under various similarity window sizes), including conventional uPDI without NLM filtering (Non-NLM), NLM filtering on the obtained power Doppler images (PD-NLM), and NLM filtering on the spatial domain of clutter filtered blood flow RF data (S-NLM). Phantom experiments, in vivo contrast-enhanced human spinal cord tumor experiments, and in vivo contrast-free human liver experiments were performed to demonstrate the superiority of the proposed St-NLM method over the other three methods. Qualitative and quantitative results show that the proposed St-NLM method can effectively suppress the background noise, improve the contrast between vessels and background, and preserve the details of small vessels at the same time. In the human liver study, the proposed St-NLM method achieves 31.05-, 24.49-, and 11.15-dB higher contrast-to-noise ratios (CNRs) and 36.86-, 36.86-, and 15.22-dB lower noise powers than Non-NLM, PD-NLM, and S-NLM, respectively. In the human spinal cord tumor, the full-width at half-maximums (FWHMs) of vessel cross Section are 76, 201, and [Formula: see text] for St-NLM, Non-NLM, and S-NLM, respectively. The proposed St-NLM method can enhance the microvascular visualization in uPDI and has the potential for the diagnosis of many microvessel-change-related diseases.

Super-Resolution Ultrasound Localization Microscopy for Visualization of the Ocular Blood Flow

OBJECTIVE

The ocular vascular system plays an important role in preserving the visual function. Alterations in either anatomy or hemodynamics of the eye may have adverse effects on vision. Thus, an imaging approach that can monitor alterations of ocular blood flow of the deep eye vasculature ranging from capillary-level vessels to large supporting vessels would be advantageous for detection of early stage retinal and optic nerve diseases.

METHODS

We propose a super-resolution ultrasound localization microscopy (ULM) technique that can assess both the microvessel and flow velocity of the deep eye with high resolution. Ultrafast plane wave imaging was acquired using an L22-14v linear array on a high frequency Verasonics Vantage system. A robust microbubble localization and tracking technique was applied to reconstruct ULM images. The experiment was first performed on pre-designed flow phantoms in vitro and then tested on a New Zealand white rabbit eye in vivo calibrated to various intraocular pressures (IOP) - 10 mmHg, 30 mmHg and 50 mmHg.

RESULTS

We demonstrated that retinal/choroidal vessels, central retinal artery, posterior ciliary artery, and vortex vein were all visible at high resolution. In addition, reduction of vascular density and flow velocity were observed with elevated IOPs.

CONCLUSION

These results indicate that super-resolution ULM is able to image the deep ocular tissue while maintaining high resolution that is comparable with optical coherence tomography angiography.

SIGNIFICANCE

Capability to detect subtle changes of blood flow may be clinically important in detecting and monitoring eye diseases such as glaucoma.

Improved Ultrasound Localization Microscopy Based on Microbubble Uncoupling via Transmit Excitation

Ultrasound localization microscopy (ULM) demonstrates great potential for visualization of tissue microvasculature at depth with high spatial resolution. The success of ULM heavily depends on robust localization of isolated microbubbles (MBs), which can be challenging in vivo especially within larger vessels where MBs can overlap and cluster close together. While MB dilution alleviates the issue of MB overlap to a certain extent, it drastically increases the data acquisition time needed for MBs to populate the microvasculature, which is already on the order of several minutes using recommended MB concentrations. Inspired by optical super-resolution imaging based on stimulated emission depletion (STED), here we propose a novel ULM imaging sequence based on MB uncoupling via transmit excitation (MUTE). MUTE "silences" MB signals by creating acoustic nulls to facilitate MB separation, which leads to robust localization of MBs especially under high concentrations. The efficiency of localization accomplished via the proposed technique was first evaluated in simulation studies with conventional ULM as a benchmark. Then, an in-vivo study based on the chorioallantoic membrane (CAM) of chicken embryos showed that MUTE could reduce the data acquisition time by half, thanks to the enhanced MB separation and localization. Finally, the performance of MUTE was validated in an in vivo mouse brain study. These results demonstrate the high MB localization efficacy of MUTE-ULM, which contributes to a reduced data acquisition time and improved temporal resolution for ULM.

Estimating the neovascularity of human finger tendon through high frequency ultrasound micro-Doppler imaging

OBJECTIVE

Neovascularization of injured tendons prolongs the proliferative phase of healing, but prolonged neovascularization may cause improper healing and pain. Currently, ultrasound Doppler imaging is used for measuring the neovascularization of injured tendons (e.g., Achilles tendon). However, the resolution of state-of-the-art clinical ultrasound machines is insufficient for visualizing the neovascularization in finger tendons. In this study, a high-resolution micro-Doppler imaging (HFDI) based on 40-MHz ultrafast ultrasound imaging was proposed for visualizing the neovascularization in injured finger tendons during multiple rehabilitation phases.

METHOD

The vessel visibility was enhanced through a block-wise singular value decomposition filter and several curvilinear structure enhancement strategies, including the bowler-hat transform and Hessian-based vessel enhancement filtering. HFDI was verified through small animal kidney and spleen imaging because the related vessel structure patterns of mice are well studied. Five patients with finger tendon injuries underwent HFDI examination at various rehabilitation phases after surgery (weeks 1156), and finger function evaluations were performed for comparisons.

RESULTS

The results of small animal experiments revealed that the proposed HFDI provides excellent microvasculature imaging performance; the contrast-to-noise ratio of HFDI was approximately 15 dB higher than that of the conventional singular value decomposition filter, and the minimum detectable vessel size for mouse kidney was 35 m without the use of contrast agent. In the human study, neovascularization was clearly observed in injured finger tendons during the early phase of healing (weeks 1121), but it regressed from week 52 to 56. Finger rehabilitation appears to help reduce neovascularization; neovascular density decreased by approximately 1.8%8.0% in participants after 4 weeks of rehabilitation.

CONCLUSION

The experimental results verified the performance of HFDI for microvasculature imaging and its potential for injured finger tendon evaluations.

Aging-related cerebral microvascular changes visualized using ultrasound localization microscopy in the living mouse

Aging-related cognitive decline is an emerging health crisis; however, no established unifying mechanism has been identified for the cognitive impairments seen in an aging population. A vascular hypothesis of cognitive decline has been proposed but is difficult to test given the requirement of high-fidelity microvascular imaging resolution with a broad and deep brain imaging field of view, which is restricted by the fundamental trade-off of imaging penetration depth and resolution. Super-resolution ultrasound localization microscopy (ULM) offers a potential solution by exploiting circulating microbubbles to achieve a vascular resolution approaching the capillary scale without sacrificing imaging depth. In this report, we apply ULM imaging to a mouse model of aging and quantify differences in cerebral vascularity, blood velocity, and vessel tortuosity across several brain regions. We found significant decreases in blood velocity, and significant increases in vascular tortuosity, across all brain regions in the aged cohort, and significant decreases in blood volume in the cerebral cortex. These data provide the first-ever ULM measurements of subcortical microvascular dynamics in vivo within the context of the aging brain and reveal that aging has a major impact on these measurements.

Compressed Sensing-Based Super-Resolution Ultrasound Imaging for Faster Acquisition and High Quality Images

GOAL

Typical SRUS images are reconstructed by localizing ultrasound microbubbles (MBs) injected in a vessel using normalized 2-dimensional cross-correlation (2DCC) between MBs signals and the point spread function of the system. However, current techniques require isolated MBs in a confined area due to inaccurate localization of densely populated MBs. To overcome this limitation, we developed the ℓ1-homotopy based compressed sensing (L1H-CS) based SRUS imaging technique which localizes densely populated MBs to visualize microvasculature in vivo.

METHODS

To evaluate the performance of L1H-CS, we compared the performance of 2DCC, interior-point method based compressed sensing (CVX-CS), and L1H-CS algorithms. Localization efficiency was compared using axially and laterally aligned point targets (PTs) with known distances and randomly distributed PTs generated by simulation. We developed post-processing techniques including clutter reduction, noise equalization, motion compensation, and spatiotemporal noise filtering for in vivo imaging. We then validated the capabilities of L1H-CS based SRUS imaging technique with high-density MBs in a mouse tumor model, kidney, and zebrafish dorsal trunk, and brain.

RESULTS

Compared to 2DCC and CVX-CS algorithms, L1H-CS achieved faster data acquisition time and considerable improvement in SRUS image quality.

CONCLUSIONS AND SIGNIFICANCE

These results demonstrate that the L1H-CS based SRUS imaging technique has the potential to examine microvasculature with reduced acquisition and reconstruction time to acquire enhanced SRUS image quality, which may be necessary to translate it into clinics.

Characterization of Anti-angiogenic Chemo-sensitization via Longitudinal Ultrasound Localization Microscopy in Colorectal Carcinoma Tumor Xenografts

OBJECTIVE

Super-resolution ultrasound localization microscopy (ULM) has unprecedented vascular resolution at clinically relevant imaging penetration depths. This technology can potentially screen for the transient microvascular changes that are thought to be critical to the synergistic effect(s) of combined chemotherapy-antiangiogenic agent regimens for cancer.

METHODS

In this paper, we apply this technology to a high-throughput colorectal carcinoma xenograft model treated with either the antiangiogenic agent sorafenib, FOLFOX-6 chemotherapy, a combination of the two treatments, or vehicle control.

RESULTS

Longitudinal ULM demonstrated morphological changes in the antiangiogenic treated cohorts, and evidence of vascular disruption caused by chemotherapy. Gold-standard histological measurements revealed reduced levels of hypoxia in the sorafenib treated cohort for both of the human cell lines tested (HCT-116 and HT-29). Therapy resistance was associated with an increase in tumor vascular fractal dimension as measured by a box-counting technique on ULM images.

CONCLUSION

These results imply that the morphological changes evident on ULM signify a functional change in the tumor microvasculature, which may be indicative of chemo-sensitivity.

SIGNIFICANCE

ULM provides additional utility for tumor therapy response evaluation by offering a myriad of morphological and functional quantitative indices for gauging treatment effect(s).

Dietary tomato inhibits angiogenesis in TRAMP prostate cancer but is not protective with a Western-style diet in this pilot study

Prostate cancer (PCa) remains the second most diagnosed cancer worldwide. Higher body weight is associated with chronic inflammation, increased angiogenesis, and treatment-resistant tumor phenotypes. Dietary tomato reduces PCa risk, which may be due to tomato inhibition of angiogenesis and disruption of androgen signaling. This pilot study investigated the interplay between tomato powder (TP), incorporated into control (CON) and obesogenic (OB) diets, and PCa tumor growth and blood perfusion over time in a transgenic model of PCa (TRAMP). Ultrasound microvessel imaging (UMI) results showed good agreement with gold-standard immunohistochemistry quantification of endothelial cell density, indicating that this technique can be applied to non-invasively monitor tumor blood perfusion in vivo. Greater body weight was positively associated with tumor growth. We also found that TP significantly inhibited prostate tumor angiogenesis but that this inhibition differentially affected measured outcomes depending on CON or OB diets. TP led to reduced tumor growth, intratumoral inflammation, and intratumoral androgen-regulated gene expression (srd5a1, srd5a2) when incorporated with the CON diet but greater tumor growth and intratumoral gene expression when incorporated with the OB diet. Results from this study show that protective benefits from dietary tomato are lost, or may become deleterious, when combined with a Western-style diet.

Real-Time In Vivo Imaging of Human Liver Vasculature Using Coherent Flow Power Doppler: A Pilot Clinical Study

Power Doppler (PD) is a commonly used technique for flow detection and vessel visualization in radiology clinics. Despite its broad set of applications, PD suffers from multiple noise sources and artifacts, such as thermal noise, clutter, and flash artifacts. In addition, a tradeoff exists between acquisition time and Doppler image quality. These limit the ability of clinical PD imaging in deep-lying and small-vessel detection and visualization, particularly among patients with high body mass indices (BMIs). To improve the Doppler vessel detection, we have previously proposed coherent flow PD (CFPD) imaging and demonstrated its performance on porcine vasculature. In this article, we report on a pilot clinical study of CFPD imaging on healthy human volunteers and patients with high BMI to assess the clinical feasibility of the technique in liver imaging. In this study, we built a real-time CFPD imaging system using a graphical processing unit (GPU)-based software beamformer and a CFPD processing module. Using the real-time CFPD imaging system, the liver vasculature of 15 healthy volunteers with normal BMI below 25 and 15 patients with BMI greater than 25 was imaged. Both PD and CFPD image streams were produced simultaneously. The generalized contrast-to-noise ratio (gCNR) of the PD and CFPD images was measured to provide the quantitative evaluation of image quality and vessel detectability. Comparison of PD and CFPD image shows that gCNR is improved by 35% in healthy volunteers and 28% in high BMI patients with CFPD compared to PD. Example images are provided to show that the improvement in the Doppler image gCNR leads to greater detection of small vessels in the liver. In addition, we show that CFPD can suppress in vivo reverberation clutter in clinical imaging.

Super-Resolution Ultrasound Localization Microscopy on a Rabbit Liver VX2 Tumor Model: An Initial Feasibility Study

Ultrasound localization microscopy can image microvasculature in vivo without sacrificing imaging penetration depth. However, the reliance on super-resolution inference limits the applicability of the technique because subpixel tissue motion can corrupt microvascular reconstruction. Consequently, the majority of previous pre-clinical research on this super-resolution procedure has been restricted to low-motion experimental models with ample motion correction or data rejection, which precludes the imaging of organ sites that exhibit a high degree of respiratory and other motion. In this article, we present a novel anesthesia protocol in rabbits that induces safe, controllable periods of apnea to enable the long image-acquisition times required for ultrasound localization microscopy. We apply this protocol to a VX2 liver tumor model undergoing sorafenib therapy and compare the results to super-resolution images from conventional high-dose isoflurane anesthesia. We find that the apneic protocol was necessary to correctly identify the poorly vascularized tumor cores, as verified by immunohistochemistry, and to reveal the tumoral microvascular architecture.

Multidimensional Clutter Filtering of Aperture Domain Data for Improved Blood Flow Sensitivity

Singular value decomposition (SVD) is a valuable factorization technique used in clutter rejection filtering for power Doppler imaging. Conventionally, SVD is applied to a Casorati matrix of radio frequency data, which enables filtering based on spatial or temporal characteristics. In this article, we propose a clutter filtering method that uses a higher order SVD (HOSVD) applied to a tensor of aperture data, e.g., delayed channel data. We discuss temporal, spatial, and aperture domain features that can be leveraged in filtering and demonstrate that this multidimensional approach improves sensitivity toward blood flow. Further, we show that HOSVD remains more robust to short ensemble lengths than conventional SVD filtering. Validation of this technique is shown using Field II simulations and in vivo data.

Simultaneous Noise Suppression and Incoherent Artifact Reduction in Ultrafast Ultrasound Vascular Imaging

Ultrasound vascular imaging based on ultrafast plane wave imaging and singular value decomposition (SVD) clutter filtering has demonstrated superior sensitivity in blood flow detection. However, ultrafast ultrasound vascular imaging is susceptible to electronic noise due to the weak penetration of unfocused waves, leading to a lower signal-to-noise ratio (SNR) at larger depths. In addition, incoherent clutter artifacts originating from strong and moving tissue scatterers that cannot be completely removed create a strong mask on top of the blood signal that obscures the vessels. Herein, a method that can simultaneously suppress the background noise and incoherent artifacts is proposed. The method divides the tilted plane or diverging waves into two subgroups. Coherent spatial compounding is performed within each subgroup, resulting in two compounded data sets. An SVD-based clutter filter is applied to each data set, followed by a correlation between the two data sets to produce a vascular image. Uncorrelated noise and incoherent artifacts can be effectively suppressed with the correlation process, while the coherent blood signal can be preserved. The method was evaluated in wire-target simulations and phantom, in which around 7-10-dB SNR improvement was shown. Consistent results were found in a flow channel phantom with improved SNR by the proposed method (39.9 ± 0.2 dB) against conventional power Doppler (29.1 ± 0.6 dB). Last, we demonstrated the effectiveness of the method combined with block-wise SVD clutter filtering in a human liver, breast tumor, and inflammatory bowel disease data sets. The improved blood flow visualization may facilitate more reliable small vessel imaging for a wide range of clinical applications, such as cancer and inflammatory diseases.

Improved Ultrasound Microvessel Imaging Using Deconvolution with Total Variation Regularization

Singular value decomposition-based clutter filters can robustly reject tissue clutter, allowing for detection of slow blood flow in imaging microvasculature. However, to identify microvessels, high ultrasound frequency must be used to increase the spatial resolution at the expense of shorter depth of penetration. Deconvolution using Tikhonov regularization is an imaging processing method widely used to improve spatial resolution. The ringing artifact of Tikhonov regularization, though, can produce image artifacts such as non-existent microvessels, which degrade image quality. Therefore, a deconvolution method using total variation is proposed in this study to improve spatial resolution and mitigate the ringing artifact. Performance of the proposed method was evaluated using chicken embryo brain, ex ovo chicken embryo chorioallantoic membrane and tumor data. Results revealed that the reconstructed power Doppler (PD) images are substantially improved in spatial resolution compared with original PD images: the full width half-maximum (FWHM) of the cross-sectional profile of a microvessel was improved from 132 to 83 µm. Two neighboring microvessels that were 154 µm apart were better separated using the proposed method than conventional PD imaging. Additionally, 223 FWHMs measured from the cross-sectional profiles of 223 vessels were used to determine the improvement in FWHM with the proposed method statistically. The mean ± standard deviation of the FWHM without and with the proposed method was 233.19 ± 85.08 and 172.31 ± 75.11 μm, respectively; the maximum FWHM without and with the proposed method was 693.01 and 668.69 μm; and the minimum FWHM without and with the proposed method was 73.92 and 45.74 μm. There were statistically significant differences between FWHMs with and without the proposed method according to the rank-sum test, p < 0.0001. The contrast-to-noise ratio improved from 1.06 to 4.03 dB with use of the proposed method. We also compared the proposed method with Tikhonov regularization using ex ovo chicken embryo chorioallantoic membrane data. We found that the proposed method outperformed Tikhonov regularization as false microvessels appeared using the Tikhonov regularization but not with the proposed method. These results indicate that the proposed method is capable of providing more robust PD images with higher spatial resolution and higher contrast-to-noise ratio.

Ultrafast Power Doppler Imaging Using Frame-Multiply-and-Sum-Based Nonlinear Compounding

Ultrafast power Doppler imaging based on coherent compounding (UPDI-CC) has become a promising technique for microvascular imaging due to its high sensitivity to slow blood flows. However, since this method utilizes a limited number of plane-wave or diverging-wave transmissions for high-frame-rate imaging, it suffers from degraded image quality because of the low contrast resolution. In this article, an ultrafast power Doppler imaging method based on a nonlinear compounding framework, called frame-multiply-and-sum (UPDI-FMAS), is proposed to improve contrast resolution. In UPDI-FMAS, unlike conventional channel-domain delay-multiply-and-sum (DMAS) beamforming, the signal coherence is estimated based on autocorrelation function over plane-wave angle frames. To avoid phase distortion of blood flow signals during the autocorrelation process, clutter filtering is preferentially applied to individual beamformed plane-wave data set. Therefore, only coherent blood flow signals are emphasized, while incoherent background noise is suppressed. The performance of the UPDI-FMAS was evaluated with simulation, phantom, and in vivo studies. For the simulation and phantom studies with a constant laminar flow, the UPDI-FMAS showed improvements in the signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) to those of UPDI-CC, i.e., over 10 and 7 dB for 13 plane waves, respectively, and the performances were improved as the number of plane waves increased. Moreover, the enhancement of the image quality due to the increased SNR and CNR in UPDI-FMAS was more clearly depicted with the in vivo study, in which a human kidney and a tumor-bearing mouse were evaluated. These results indicate that the FMAS compounding can improve the image quality of UPDI for microvascular imaging without loss of temporal resolution.

Blood Flow Imaging in the Neonatal Brain Using Angular Coherence Power Doppler

Using ultrasound to image small vessels in the neonatal brain can be difficult in the presence of strong clutter from the surrounding tissue and with a neonate motion during the scan. We propose a coherence-based beamforming method, namely the short-lag angular coherence (SLAC) beamforming that suppresses incoherent noise and motion artifacts in Ultrafast data, and we demonstrate its applicability to improve detection of blood flow in the neonatal brain. Instead of estimating spatial coherence across the receive elements, SLAC utilizes the principle of acoustic reciprocity to estimate angular coherence from the beamsummed signals from different plane-wave transmits, which makes it computationally efficient and amenable to advanced beamforming techniques, such as f-k migration. The SLAC images of a simulated speckle phantom show similar edge resolution and texture size as the matching B-mode images, and reduced random noise in the background. We apply SLAC power Doppler (PD) to free-hand imaging of neonatal brain vasculature with long Doppler ensembles and show that: 1) it improves visualization of small vessels in the cortex compared to conventional PD and 2) it can be used for tracking of blood flow in the brain over time, meaning it could potentially improve the quality of free-hand functional ultrasound.

Real time SVD-based clutter filtering using randomized singular value decomposition and spatial downsampling for micro-vessel imaging on a Verasonics ultrasound system

Singular value decomposition (SVD)-based clutter filters can robustly reject the tissue clutter as compared with the conventional high pass filter-based clutter filters. However, the computational burden of SVD makes real time SVD-based clutter filtering challenging (e.g. frame rate at least 10-15 Hz with region of interest of about 4 × 4 cm²). Recently, we proposed an acceleration method based on randomized SVD (rSVD) clutter filtering and randomized spatial downsampling, which can significantly reduce the computational complexity without compromising the clutter rejection capability. However, this method has not been implemented on an ultrasound scanner and tested for its performance. In this study, we implement this acceleration method on a Verasonics scanner using a multi-core CPU architecture, and evaluate the selections of the imaging and processing parameters to enable real time micro-vessel imaging. The Blood-to-Clutter Ratio (BCR) performance was evaluated on a Verasonics machine with different settings of parameters such as block size and ensemble size. The demonstration of real time process was implemented on a 12-core CPU (downsampling factor of 12, 12-threads in this study) host computer. The processing time of the rSVD-based clutter filter was less than 30 ms and BCRs were higher than 20 dB as the block size, ensemble size and the rank of tissue clutter subspace were set as 30 × 30, 45 and 26 respectively. We also demonstrate that the micro-vessel imaging frame rate of the proposed architecture can reach approximately 22 Hz when the block size, ensemble size and the rank of tissue clutter subspace were set as 20 × 20 pixels, 45 and 26 respectively (using both images and supplementary videos). The proposed method may be important for real time 2D scanning of tumor microvessels in 3D to select and store the most representative 2D view with most abnormal micro-vessels for better diagnosis.

Short Acquisition Time Super-Resolution Ultrasound Microvessel Imaging via Microbubble Separation

Super-resolution ultrasound localization microscopy (ULM), based on localization and tracking of individual microbubbles (MBs), offers unprecedented microvascular imaging resolution at clinically relevant penetration depths. However, ULM is currently limited by the requirement of dilute MB concentrations to ensure spatially sparse MB events for accurate localization and tracking. The corresponding long imaging acquisition times (tens of seconds or several minutes) to accumulate sufficient isolated MB events for full reconstruction of microvasculature preclude the clinical translation of the technique. To break this fundamental tradeoff between acquisition time and MB concentration, in this paper we propose to separate spatially overlapping MB events into sub-populations, each with sparser MB concentration, based on spatiotemporal differences in the flow dynamics (flow speeds and directions). MB localization and tracking are performed for each sub-population separately, permitting more robust ULM imaging of high-concentration MB injections. The superiority of the proposed MB separation technique over conventional ULM processing is demonstrated in flow channel phantom data, and in the chorioallantoic membrane of chicken embryos with optical imaging as an in vivo reference standard. Substantial improvement of ULM is further demonstrated on a chicken embryo tumor xenograft model and a chicken brain, showing both morphological and functional microvasculature details at super-resolution within a short acquisition time (several seconds). The proposed technique allows more robust MB localization and tracking at relatively high MB concentrations, alleviating the need for dilute MB injections, and thereby shortening the acquisition time of ULM imaging and showing great potential for clinical translation.

Ultrasound localization microscopy of renal tumor xenografts in chicken embryo is correlated to hypoxia

Ultrasound localization microscopy (ULM) permits the reconstruction of super-resolved microvascular images at clinically relevant penetration depths, which can be potentially leveraged to provide non-invasive quantitative measures of tissue hemodynamics and hypoxic status. We demonstrate that ULM microbubble data processing methods, applied to images acquired with a Verasonics Vantage 256 system, can provide a non-invasive imaging surrogate biomarker of tissue oxygenation status. This technique was applied to evaluate the microvascular structure, vascular perfusion, and hypoxia of a renal cell carcinoma xenograft model grown in the chorioallantoic membrane of chicken embryos. Histological microvascular density was significantly correlated to ULM measures of intervessel distance (R = -0.92, CI95 = [-0.99,-0.42], p = 0.01). The Distance Metric, a measure of vascular tortuosity, was found to be significantly correlated to hypoxyprobe quantifications (R = 0.86, CI95 = [0.17, 0.99], p = 0.03). ULM, by providing non-invasive in vivo microvascular structural information, has the potential to be a crucial clinical imaging modality for the diagnosis and therapy monitoring of solid tumors.

Adaptive background noise bias suppression in contrast-free ultrasound microvascular imaging

Non-invasive, contrast-free imaging of small vessel blood flow is diagnostically invaluable for detection, diagnosis and monitoring of disease. Recent advances in ultrafast imaging and tissue clutter-filtering have considerably improved the sensitivity of power Doppler (PD) imaging in detecting small vessel blood flow. However, suppression of tissue clutter exposes the depth-dependent time-gain compensated noise bias that noticeably degrades the PD image. We hypothesized that background suppression of PD images based on noise bias estimated from the entire clutter-filtered singular value spectrum can considerably improve flow signal visualization compared to currently existing techniques. To test our hypothesis, in vivo experiments were conducted on suspicious breast lesions in 10 subjects and deep-seated hepatic and renal microvasculatures in four healthy volunteers. Ultrasound PD images were acquired using a clinical ultrasound scanner, implemented with compounded plane wave imaging. The time gain compensated noise field was computed from the clutter-filtered Doppler ensemble (CFDE) based on its local spatio-temporal correlation, combined with low-rank signal estimation. Subsequently, the background bias in the PD images was suppressed by subtracting the estimated noise field. Background-suppressed PD images obtained using the proposed technique substantially improved visualization of the blood flow signal. The background bias in the noise suppressed PD images varied  <0.6 dB, independent of depth, which otherwise increased up to 13.8 dB. Further, the results demonstrated that the proposed technique efficaciously suppressed the background noise bias associated with smaller Doppler ensembles, which are challenging due to increased overlap between blood flow and noise components in the singular value spectrum. These preliminary results demonstrate the utility of the proposed technique to improve the visualization of small vessel blood flow in contrast-free PD images. The results of this feasibility study were encouraging, and warrant further development and additional in vivo validation.

Wide Field-of-View Ultrafast Curved Array Imaging Using Diverging Waves

Ultrafast ultrasound imaging provides great opportunities for very high frame rate applications, such as shear wave elastography and microvascular imaging. However, ultrafast imaging with curved array transducers remains challenging in terms of element directivity and a limited field-of-view (FOV) for a fully synthetic area. In this paper, a wide FOV ultrafast curved array imaging method based on diverging wave transmissions is presented for high frame rate abdominal ultrasound applications. For this method, a theoretical model for a diverging wave solution based on a virtual point source originating from a circular line is proposed, and the FOV and element directivity are analyzed by this model. Furthermore, an integrated model for plane wave and diverging wave imaging along the location of the virtual point source is derived. The proposed method was evaluated with simulation, phantom, and in vivo studies. In the simulation and phantom studies, the image quality (i.e., spatial resolution, cystic resolution, and contrast-to-noise ratio), and effective FOV were assessed. For the in vivo study, a preliminary result from abdominal microvascular imaging, where diverging wave excitation was utilized to depict the vasculature, was also presented. In the renal cortex microvessels, the diverging wave imaging yielded a higher signal-to-clutter ratio value than the plane wave imaging, i.e., 6.35 vs. 4.26 dB, due to the wider synthetic field. These studies demonstrated that the proposed ultrafast curved array imaging technique based on diverging wave excitation allowed for an extended FOV with high spatiotemporal resolution.

Debiasing-Based Noise Suppression for Ultrafast Ultrasound Microvessel Imaging

Ultrasound microvessel imaging (UMI) based on the combination of singular value decomposition (SVD) clutter filtering and ultrafast plane wave imaging has recently demonstrated significantly improved Doppler sensitivity, especially to small vessels that are invisible to conventional Doppler imaging. Practical implementation of UMI is hindered by the high computational cost associated with SVD and low blood signal-to-noise ratio (SNR) in deep regions of the tissue due to the lack of transmit focusing of plane waves. Concerning the high computational cost, an accelerated SVD clutter filtering method based on randomized SVD (rSVD) and randomized spatial downsampling (rSD) was recently proposed by our group, which showed the feasibility of real-time implementation of UMI. Concerning the low blood flow SNR in deep imaging regions, here we propose a noise suppression method based on noise debiasing that can be easily applied to the accelerated SVD method to bridge the gap between real-time implementation and high imaging quality. The proposed method experimentally measures the noise-induced bias by collecting the noise signal using the identical imaging sequence as regular UMI, but with the ultrasound transmission turned off. The estimated bias can then be subtracted from the original power Doppler (PD) image to obtain effective noise suppression. The feasibility of the proposed method was validated under different ultrasound imaging parameters [including transmitting voltages and time-gain compensation (TGC) settings] with a phantom experiment. The noise-debiased images showed an increase of up to 15.3 and 13.4 dB in SNR as compared to original PD images on the blood flow phantom and an in vivo human kidney data set, respectively. The proposed noise suppression method has negligible computational cost and can be conveniently combined with the previously proposed accelerated SVD clutter filtering technique to achieve high quality, real-time UMI imaging.

Functional Ultrasound Imaging of Spinal Cord Hemodynamic Responses to Epidural Electrical Stimulation: A Feasibility Study

This study presents the first implementation of functional ultrasound (fUS) imaging of the spinal cord to monitor local hemodynamic response to epidural electrical spinal cord stimulation (SCS) on two small and large animal models. SCS has been successfully applied to control chronic refractory pain and recently was evolved to alleviate motor impairment in Parkinson's disease and after spinal cord injury. At present, however, the mechanisms underlying SCS remain unclear, and current methods for monitoring SCS are limited in their capacity to provide the required sensitivity and spatiotemporal resolutions to evaluate functional changes in response to SCS. fUS is an emerging technology that has recently shown promising results in monitoring a variety of neural activities associated with the brain. Here we demonstrated the feasibility of performing fUS on two animal models during SCS. We showed in vivo spinal cord hemodynamic responses measured by fUS evoked by different SCS parameters. We also demonstrated that fUS has a higher sensitivity in monitoring spinal cord response than electromyography. The high spatial and temporal resolutions of fUS were demonstrated by localized measurements of hemodynamic responses at different spinal cord segments, and by reliable tracking of spinal cord responses to patterned electrical stimulations, respectively. Finally, we proposed optimized fUS imaging and post-processing methods for spinal cord. These results support feasibility of fUS imaging of the spinal cord and could pave the way for future systematic studies to investigate spinal cord functional organization and the mechanisms of spinal cord neuromodulation in vivo.

Noninvasive Contrast-Free 3D Evaluation of Tumor Angiogenesis with Ultrasensitive Ultrasound Microvessel Imaging

Ultrasound microvessel imaging (UMI), when applied with ultrafast planewave acquisitions, has demonstrated superior blood signal sensitivity in comparison to conventional Doppler imaging. Here we propose a high spatial resolution and ultra-sensitive UMI that is based on conventional line-by-line high-frequency ultrasound imagers and singular value decomposition (SVD) clutter filtering for the visualization and quantification of tumor microvasculature and perfusion. The technique was applied to a chicken embryo tumor model of renal cell carcinoma that was treated with two FDA-approved anti-angiogenic agents at clinically relevant dosages. We demonstrate the feasibility of 3D evaluation with UMI to achieve highly sensitive detection of microvasculature using conventional line-by-line ultrasound imaging on a preclinical and commercially available high-frequency ultrasound device without software or hardware modifications. Quantitative parameters (vascularization index and fractional moving blood volume) derived from UMI images provide significantly improved evaluation of anti-angiogenic therapy response as compared with conventional power Doppler imaging, using histological analysis and immunohistochemistry as the reference standard. This proof-of-concept study demonstrates that high-frequency UMI is a low-cost, contrast-agent-free, easily applicable, accessible, and quantitative imaging tool for tumor characterization, which may be very useful for preclinical evaluation and longitudinal monitoring of anti-cancer treatment.

Non-Local Based Denoising Framework for In Vivo Contrast-Free Ultrasound Microvessel Imaging

Vascular networks can provide invaluable information about tumor angiogenesis. Ultrafast Doppler imaging enables ultrasound to image microvessels by applying tissue clutter filtering methods on the spatio-temporal data obtained from plane-wave imaging. However, the resultant vessel images suffer from background noise that degrades image quality and restricts vessel visibilities. In this paper, we addressed microvessel visualization and the associated noise problem in the power Doppler images with the goal of achieving enhanced vessel-background separation. We proposed a combination of patch-based non-local mean filtering and top-hat morphological filtering to improve vessel outline and background noise suppression. We tested the proposed method on a flow phantom, as well as in vivo breast lesions, thyroid nodules, and pathologic liver from human subjects. The proposed non-local-based framework provided a remarkable gain of more than 15 dB, on average, in terms of contrast-to-noise and signal-to-noise ratios. In addition to improving visualization of microvessels, the proposed method provided high quality images suitable for microvessel morphology quantification that may be used for diagnostic applications.