3-D Velocity and Volume Flow Measurement In Vivo Using Speckle Decorrelation and 2-D High-Frame-Rate Contrast-Enhanced Ultrasound

Ultrasound Image Velocimetry for High Spatiotemporal Resolution Blood Flow Velocity Field Mapping in Mice

OBJECTIVE

Abnormal hemodynamics is thought to play an essential role in the development of cardiovascular diseases. Mouse models are widely used for elucidating the underlying mechanisms; however, their small size and high heart rates make it difficult to perform quantitative flow velocity field mapping with sufficient temporal resolution. Our objective was to develop a noninvasive method for quantitative flow field mapping in mice based on speckle-tracking from high-frequency ultrasound B-mode imaging.

METHODS

Ultrasound ECG-gated kilohertz visualization (EKV) was performed on a mouse-aorta-sized tubular flow phantom at frame rates up to 10,000 fps. Unexpected velocity underestimations were elucidated by simulating EKV reconstruction and performing ultrasound image velocimetry (UIV) in silico. A technique for error correction was developed and validated in vitro, and demonstrated in vivo.

RESULTS

In flow phantoms, EKV-UIV underestimated velocity in the beam lateral direction by 50%-70%. This was attributed to loss of speckle contiguity owing to EKV's retrospective strip-based reconstruction of the two-dimensional B-mode image. The proposed correction technique reduced the errors to <10% by accounting only for speckle movement within each image strip. A preliminary in vivo study showed that vortex shapes and near-wall expansion movement inside a mouse left ventricle were more aligned with physical expectations after correction.

CONCLUSION

A novel technique was developed to quantitatively map blood flow with high spatiotemporal resolution. Further optimization will enable longitudinal studies in mice to gain insights on the role of local hemodynamic forces in the development of cardiovascular diseases.

Angle-Independent Blood Flow Velocity Measurement With Ultrasound Speckle Decorrelation Analysis

Precise measurement of the blood flow velocity in major arteries is important for the assessment of circulation dysfunction but challenging when using a one-dimensional (1D) ultrasound transducer array. Current available ultrasound velocimetry methods are susceptible to the probe-to-vessel angle and require the vessels to be well-aligned within the imaging plane of the 1D transducer array. In this study, a novel angle-independent velocimetry (VT-vUS) based on the ultrasound speckle decorrelation analysis of the ultrasound field signal is proposed to measure the blood flow velocity using a conventional 1D ultrasound transducer array. We first introduced the principle and evaluated this technique with numerical simulation and phantom experiments, which demonstrated that VT-vUS can accurately reconstruct the velocity magnitude of blood flow at arbitrary probe-to-vessel angles for different preset flow speeds (up to ~2.5 m/s). Further, we applied VT-vUS to measure the pulsatile flow of the radial artery and carotid artery in a healthy volunteer. Results show that the absolute velocity profiles obtained with VT-vUS at different probe-to-vessel angles have high consistency and agree well with the absolute speed obtained with the color Doppler-corrected velocimetry throughout the cardiac cycle. With the ability to alleviate the dependency on probe-to-vessel angle, VT-vUS has the potential for circulation-related disease screening in clinical practices.

Measurement of Wall Shear Rate Across the Entire Vascular Wall Using Ultrasound Speckle Decorrelation

OBJECTIVE

The accurate measurement of the wall shear rate (WSR) plays a crucial role in the early diagnosis of cardiovascular disease progression and acute events such as aneurysms and atherosclerotic plaque ruptures. To address this need, the speckle decorrelation (SDC) technique has been used to measure WSR based on the 2-D out-of-plane blood flow speed. This technique is particularly advantageous because it enables the use of a 1-D array transducer to measure WSR over the entire luminal area. This study aims to develop a region-based singular value decomposition (SVD) filtering technique that selectively suppresses clutter noise in the vascular region to measure WSR using SDC.

METHOD

Ultrasound simulations, in-vitro flow experiments, and an in-vivo human study were conducted to evaluate the feasibility of this method's clinical application.

RESULTS

The results demonstrated that WSR can be effectively measured across entire vascular walls using a conventional 1-D array transducer along with the proposed methodology.

CONCLUSION

This study successfully demonstrates a noninvasive and accurate SDC-based method for measuring vital vascular WSR. This approach holds significant promise for assessing vascular WSR in both healthy individuals and high-risk cardiovascular disease patients.

Ultrasound Speckle Decorrelation Analysis-Based Velocimetry for 3D-Velocity-Components Measurement Using a 1D Transducer Array

Ultrasound velocimetry has been widely used for blood flow imaging. However, the flow measurements are constrained to resolve the in-plane 2D flow components when using a 1D transducer array. In this work, an ultrasound speckle decorrelation analysis-based velocimetry (3C-vUS) is proposed for 3D velocity components measurement using a 1D transducer array. The 3C-vUS theory is first derived and validated with numerical simulations and phantom experiments. The in vivo testing results show that 3C-vUS can accurately measure the blood flow 3D-velocity-components of the human carotid artery at arbitrary probe-to-vessel angles throughout the cardiac cycle. With such capability, the 3C-vUS will alleviate the requirement of operators and promote disease screening for blood flow-related disorders.

Time-Lagged Functional Ultrasound for Multi-Parametric Cerebral Hemodynamic Imaging

We introduce an ultrasound speckle decorrelation-based time-lagged functional ultrasound technique (tl-fUS) for the quantification of the relative changes in cerebral blood flow speed (rCBF [Formula: see text]), cerebral blood volume (rCBV) and cerebral blood flow (rCBF) during functional stimulations. Numerical simulations, phantom validations, and in vivo mouse brain experiments were performed to test the capability of tl-fUS to parse out and quantify the ratio change of these hemodynamic parameters. The blood volume change was found to be more prominent in arterioles compared to venules and the peak blood flow changes were around 2.5 times the peak blood volume change during brain activation, agreeing with previous observations in the literature. The tl-fUS shows the ability of distinguishing the relative changes of rCBFspeed, rCBV, and rCBF, which can inform specific physiological interpretations of the fUS measurements.

Transverse oscillation tensor velocity imaging using a row-column addressed array: Experimental validation

Tensor velocity imaging (TVI) performance with a row-column probe was assessed for constant flow in a straight vessel phantom and pulsatile flow in a carotid artery phantom. TVI, i.e., estimating the 3-D velocity vector as a function of time and spatial position, was performed using the transverse oscillation cross-correlation estimator, and the flow was acquired with a Vermon 128+128 row-column array probe connected to a Verasonics 256 research scanner. The emission sequence used 16 emissions per image, and a TVI volume rate of 234 Hz was obtained for a pulse repetition frequency (f_prf) of 15 kHz. The TVI was validated by comparing estimates of the flow rate through several cross-sections with the flow rate set by the pump. For the constant 8 mL/s flow in the straight vessel phantom with relative estimator bias (RB) and standards deviation (RSD) was found in the range of -2.18% to 0.55% and 4.58% to 2.48% in measurements performed with an f_prf of 15, 10, 8, and 5 kHz. The pulsatile flow in the carotid artery phantom the was set to an average flow rate of 2.44 mL/s, and the flow was acquired with an f_prf of 15, 10, and 8 kHz. The pulsatile flow was estimated from two measurement sites: one at a straight section of the artery and one at the bifurcation. In the straight section, the estimator predicted the average flow rate with an RB value ranging from -7.99% to 0.10% and an RSD value ranging from 10.76% to 6.97%. At the bifurcation, RB and RSD values were between -7.47% to 2.02% and 14.46% to 8.89%. This demonstrates that an RCA with 128 receive elements can accurately capture the flow rate through any cross-section at a high sampling rate.

Ultrasound Speckle Decorrelation-Based Blood Flow Measurements

Ultrasound imaging is the preferred noninvasive technique to measure blood flow to diagnose cardiovascular disease such as heart failure, carotid stenosis, and renal failure. Conventional ultrasound techniques such as Doppler ultrasound, ultrasound imaging velocimetry, vector Doppler and transverse oscillation beamforming have been used for blood flow velocity profile measurement. However, these techniques were limited to measuring blood flow velocities within the 2-D lateral (across the ultrasound beam) plane of a vessel, and the blood flow velocity profile was derived by assuming that blood vessels have a circular cross-section with axis symmetry. This assumption is incorrect because most vessels have complex geometries, such as tortuosity and branches, and an asymmetric flow profile in the presence of vascular plaque. Consequently, ultrasound speckle decorrelation has been proposed to measure blood flow from transverse views of blood vessels wherein the ultrasound beam is perpendicular to the vessel axis. In this review, we present a summary of recent progress in ultrasound speckle decorrelation-based blood flow measurement techniques.

Velocity field estimation in transcranial small vessel using super-resolution ultrasound imaging velocimetry

Based on the diameter and position information of small vessels obtained by transcranial super-resolution imaging using 3 MHz low-frequency chirp plane waves, a Gaussian-like non-linear compression was adopted to compress the blood flow signals in spatiotemporal filtering (STF) data to a precise region, and then estimate the blood flow velocity field inside the region over the adjacent time intervals using ultrasound imaging velocimetry (UIV). Imaging parameters, such as the mechanical index (MI), frame rate, and microbubble (MB) concentration, are critical during the estimation of velocity fields over a short time at high MB contrast agent concentrations. These were optimized through experiments and algorithms, in which dividing the connected domain was proposed to calculate MB cluster spot centroid spacing (SCS) and the spot-to-flow area ratio (SFAR) to determine the suitable MB concentration. The results of the in vitro experiments showed that the estimation of the small vessel flow velocity field was consistent with the theoretical results; the velocity field resolution for vessels with diameters of 0.5 mm and 0.3 mm was 36 μm and 21 μm, and the error between the mean velocity and the theoretical value was 0.7 % and 0.67 %, respectively.

Fast and Selective Super-Resolution Ultrasound In Vivo With Acoustically Activated Nanodroplets

Perfusion by the microcirculation is key to the development, maintenance and pathology of tissue. Its measurement with high spatiotemporal resolution is consequently valuable but remains a challenge in deep tissue. Ultrasound Localization Microscopy (ULM) provides very high spatiotemporal resolution but the use of microbubbles requires low contrast agent concentrations, a long acquisition time, and gives little control over the spatial and temporal distribution of the microbubbles. The present study is the first to demonstrate Acoustic Wave Sparsely-Activated Localization Microscopy (AWSALM) and fast-AWSALM for in vivo super-resolution ultrasound imaging, offering contrast on demand and vascular selectivity. Three different formulations of acoustically activatable contrast agents were used. We demonstrate their use with ultrasound mechanical indices well within recommended safety limits to enable fast on-demand sparse activation and destruction at very high agent concentrations. We produce super-localization maps of the rabbit renal vasculature with acquisition times between 5.5 s and 0.25 s, and a 4-fold improvement in spatial resolution. We present the unique selectivity of AWSALM in visualizing specific vascular branches and downstream microvasculature, and we show super-localized kidney structures in systole (0.25 s) and diastole (0.25 s) with fast-AWSALM outperforming microbubble based ULM. In conclusion, we demonstrate the feasibility of fast and selective imaging of microvascular dynamics in vivo with subwavelength resolution using ultrasound and acoustically activatable nanodroplet contrast agents.

Evaluation of aortic stenosis: From Bernoulli and Doppler to Navier-Stokes

Uni-dimensional Doppler echocardiography data provide the mainstay of quantative assessment of aortic stenosis, with the transvalvular pressure drop a key indicator of haemodynamic burden. Sophisticated methods of obtaining velocity data, combined with improved computational analysis, are facilitating increasingly robust and reproducible measurement. Imaging modalities which permit acquisition of three-dimensional blood velocity vector fields enable angle-independent valve interrogation and calculation of enhanced measures of the transvalvular pressure drop. This manuscript clarifies the fundamental principles of physics that underpin the evaluation of aortic stenosis and explores modern techniques that may provide more accurate means to grade aortic stenosis and inform appropriate management.

Recent advances in vascular ultrasound imaging technology and their clinical implications

In this paper recent advances in vascular ultrasound imaging technology are discussed, including three-dimensional ultrasound (3DUS), contrast-enhanced ultrasound (CEUS) and strain- (SE) and shear-wave-elastography (SWE). 3DUS imaging allows visualisation of the actual 3D anatomy and more recently of flow, and assessment of geometrical, morphological and mechanical features in the carotid artery and the aorta. CEUS involves the use of microbubble contrast agents to estimate sensitive blood flow and neovascularisation (formation of new microvessels). Recent developments include the implementation of computerised tools for automated analysis and quantification of CEUS images, and the possibility to measure blood flow velocity in the aorta. SE, which yields anatomical maps of tissue strain, is increasingly being used to investigate the vulnerability of the carotid plaque, but is also promising for the coronary artery and the aorta. SWE relies on the generation of a shear wave by remote acoustic palpation and its acquisition by ultrafast imaging, and is useful for measuring arterial stiffness. Such advances in vascular ultrasound technology, with appropriate validation in clinical trials, could positively change current management of patients with vascular disease, and improve stratification of cardiovascular risk.

Volumetric Flow Estimation in a Coronary Artery Phantom Using High-Frame-Rate Contrast-Enhanced Ultrasound, Speckle Decorrelation, and Doppler Flow Direction Detection

The coronary flow reserve (CFR), relating to the volumetric flow rate, is an effective functional parameter to assess the stenosis in the left anterior descending (LAD) coronary artery. We have recently proposed to use high-frame-rate (HFR) contrast-enhanced ultrasound (CEUS) to estimate the volumetric flow rate using ultrasound (US) speckle decorrelation (SDC) without any assumptions about the velocity profile. However, this method still has challenges in imaging deep and small vessels, such as LAD. In this study, we proposed to address the challenges and demonstrate the feasibility of volumetric flow rate measurement in a coronary mimicking phantom with pulsatile flow using a 1-D array cardiac probe, vector Doppler, and an optimal probe rotation/tilting for flow direction detection. Both simulations and in vitro experiments were conducted to validate the proposed method. It is shown that in-plane velocities estimated by vector Doppler under a 10° probe tilting resulted in smaller percentage error (+5.2%) in flow rate estimates than that in US imaging velocimetry (-20.2%) although their relative standard deviations were very close, being 2.6 and 2.8 ml/min, respectively. The flow rate estimated by SDC without direction detection had an error higher than 70%. A 10° tilting of the probe had the best results in flow rate estimation compared to the 5° or 15° tilting. Realistic global motions in the LAD increased the flow rate estimation error from 5.2% to 14.2%. It is concluded that it is feasible to measure the volumetric flow rate in a coronary artery flow phantom with a conventional cardiac probe, using HFR acquisition, Doppler, and SDC analysis. Potentially, this technique could also be applied to investigate the volumetric flow rate in other small vessels similar to the LAD.

Optimum Speckle Tracking Based on Ultrafast Ultrasound for Improving Blood Flow Velocimetry

Speckle tracking using optimum comparison frames (STO) is proposed to improve the blood flow velocity profile (BFVP) estimation based on ultrafast ultrasound with coherent plane-wave compounding. The optimum comparison frames are as far as possible from the reference frame image while possessing a speckle correlation above a given threshold. The correlation thresholds for different kernel sizes are determined via an experiment based on a vascular-mimicking phantom. In in vitro experiments with different peak velocities of the flow ranging from 0.38 to 1.18 m/s, the proposed STO method with three kernel sizes ( 0.46 × 0.46 , 0.31 × 0.69 , and 0.92 × 0.92 mm²) is used for the BFVP estimations. The normalized root mean square errors (NRMSEs) between the estimated and theoretical BFVPs are calculated and compared with the results based on the speckle tracking using adjacent-frame images. For the three kernel sizes, the mean relative decrements in the STO-based NRMSEs are 46.6%, 44.7%, and 52.9%, and the standard deviations are 36.8%, 37.6%, and 35.9%, respectively. The STO method is also validated by in vivo experiments using rabbit iliac arteries with contrast agents. With parabolic curves fitting to the mean velocity estimates, the average relative increments for the STO-based R² (coefficients of determination) are 7.22% and 6.25% for kernel sizes of 0.46 × 0.46 and 0.31 × 0.69 mm², respectively. In conclusion, the STO method improves the BFVP measurement accuracy, whereby accurate diagnosis information can be acquired for clinical applications.

Contrast-Enhanced High-Frame-Rate Ultrasound Imaging of Flow Patterns in Cardiac Chambers and Deep Vessels

Cardiac function and vascular function are closely related to the flow of blood within. The flow velocities in these larger cavities easily reach 1 m/s, and generally complex spatiotemporal flow patterns are involved, especially in a non-physiologic state. Visualization of such flow patterns using ultrasound can be greatly enhanced by administration of contrast agents. Tracking the high-velocity complex flows is challenging with current clinical echographic tools, mostly because of limitations in signal-to-noise ratio; estimation of lateral velocities; and/or frame rate of the contrast-enhanced imaging mode. This review addresses the state of the art in 2-D high-frame-rate contrast-enhanced echography of ventricular and deep-vessel flow, from both technological and clinical perspectives. It concludes that current advanced ultrasound equipment is technologically ready for use in human contrast-enhanced studies, thus potentially leading to identification of the most clinically relevant flow parameters for quantifying cardiac and vascular function.

Functional Ultrasound Speckle Decorrelation-Based Velocimetry of the Brain

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

Determining Haemodynamic Wall Shear Stress in the Rabbit Aorta In Vivo Using Contrast-Enhanced Ultrasound Image Velocimetry

Abnormal blood flow and wall shear stress (WSS) can cause and be caused by cardiovascular disease. To date, however, no standard method has been established for mapping WSS in vivo. Here we demonstrate wide-field assessment of WSS in the rabbit abdominal aorta using contrast-enhanced ultrasound image velocimetry (UIV). Flow and WSS measurements were made independent of beam angle, curvature or branching. Measurements were validated in an in silico model of the rabbit thoracic aorta with moving walls and pulsatile flow. Mean errors over a cardiac cycle for velocity and WSS were 0.34 and 1.69%, respectively. In vivo time average WSS in a straight segment of the suprarenal aorta correlated highly with simulations (PC = 0.99) with a mean deviation of 0.29 Pa or 5.16%. To assess fundamental plausibility of the measurement, UIV WSS was compared to an analytic approximation derived from the Poiseuille equation; the discrepancy was 17%. Mapping of WSS was also demonstrated in regions of arterial branching. High time average WSS (TAWSS_xz = 3.4 Pa) and oscillatory flow (OSI_xz = 0.3) were observed near the origin of conduit arteries. In conclusion, we have demonstrated that contrast-enhanced UIV is capable of measuring spatiotemporal variation in flow velocity, arterial wall location and hence WSS in vivo with high accuracy over a large field of view.

Optimization of 3-D Divergence-Free Flow Field Reconstruction Using 2-D Ultrasound Vector Flow Imaging

3-D blood vector flow imaging is of great value in understanding and detecting cardiovascular diseases. Currently, 3-D ultrasound vector flow imaging requires 2-D matrix probes, which are expensive and suffer from suboptimal image quality. Our recent study proposed an interpolation algorithm to obtain a divergence-free reconstruction of the 3-D flow field from 2-D velocities obtained by high-frame-rate ultrasound particle imaging velocimetry (High Frame Rate echo-Particle Imaging Velocimetry, also known as HFR Ultrasound Imaging Velocimetry (UIV)), using a 1-D array transducer. The aim of this work was to significantly improve the accuracy and reduce the time-to-solution of our previous approach, thereby paving the way for clinical translation. More specifically, accuracy was improved by optimising the divergence-free basis to reduce Runge phenomena near domain boundaries, and time-to-solution was reduced by demonstrating that under certain conditions, the resulting system could be solved using widely available and highly optimised generalised minimum residual algorithms. To initially illustrate the utility of the approach, coarse 2-D subsamplings of an analytical unsteady Womersely flow solution and a steady helical flow solution obtained using computational fluid dynamics were used successfully to reconstruct full flow solutions, with 0.82% and 4.8% average relative errors in the velocity field, respectively. Subsequently, multiplane 2-D velocity fields were obtained through HFR UIV for a straight-tube phantom and a carotid bifurcation phantom, from which full 3-D flow fields were reconstructed. These were then compared with flow fields obtained via computational fluid dynamics in each of the two configurations, and average relative errors of 6.01% and 12.8% in the velocity field were obtained. These results reflect 15%-75% improvements in accuracy and 53- to 874-fold acceleration of reconstruction speeds for the four cases, compared with the previous divergence-free flow reconstruction method. In conclusion, the proposed method provides an effective and fast method to reconstruct 3-D flow in arteries using a 1-D array transducer.

Measurement of Flow Volume in the Presence of Reverse Flow with Ultrasound Speckle Decorrelation

Direct measurement of volumetric flow rate in the cardiovascular system with ultrasound is valuable but has been a challenge because most current 2-D flow imaging techniques are only able to estimate the flow velocity in the scanning plane (in-plane). Our recent study demonstrated that high frame rate contrast ultrasound and speckle decorrelation (SDC) can be used to accurately measure the speed of flow going through the scanning plane (through-plane). The volumetric flow could then be calculated by integrating over the luminal area, when the blood vessel was scanned from the transverse view. However, a key disadvantage of this SDC method is that it cannot distinguish the direction of the through-plane flow, which limited its applications to blood vessels with unidirectional flow. Physiologic flow in the cardiovascular system could be bidirectional due to its pulsatility, geometric features, or under pathologic situations. In this study, we proposed a method to distinguish the through-plane flow direction by inspecting the flow within the scanning plane from a tilted transverse view. This method was tested on computer simulations and experimental flow phantoms. It was found that the proposed method could detect flow direction and improved the estimation of the flow volume, reducing the overestimation from over 100% to less than 15% when there was flow reversal. This method showed significant improvement over the current SDC method in volume flow estimation and can be applied to a wider range of clinical applications where bidirectional flow exists.