A Comparison Between Compounding Techniques Using Large Beam-Steered Plane Wave Imaging for Blood Vector Velocity Imaging in a Carotid Artery Model

Reconstructive interpolation for pulse wave estimation to improve local PWV measurement of carotid artery

Ultrasonic transit time (TT)-based local pulse wave velocity (PWV) measurement is defined as the distance between two beam positions on a segment of common carotid artery (CCA) divided by the TT in the pulse wave propagation. However, the arterial wall motions (AWMs) estimated from ultrasonic radio frequency (RF) signals with a limited number of frames using the motion tracking are typically discrete. In this work, we develop a method involving motion tracking combined with reconstructive interpolation (MTRI) to reduce the quantification errors in the estimated PWs, and thereby improve the accuracy of the TT-based local PWV measurement for CCA. For each beam position, normalized cross-correlation functions (NCCFs) between the reference (the first frame) and comparison (the remaining frames) RF signals are calculated. Thereafter, the reconstructive interpolation is performed in the neighborhood of the NCCFs' peak to identify the interpolation-deduced peak locations, which are more exact than the original ones. According to which, the improved AWMs are obtained to calculate their TT along a segment of the CCA. Finally, the local PWV is measured by applying a linear regression fit to the time-distance result. In ultrasound simulations based on the pulse wave propagation models of young, middle-aged, and elderly groups, the MTRI method with different numbers of interpolated samples was used to estimate AWMs and local PWVs. Normalized root mean squared errors (NRMSEs) between the estimated and preset values of the AWMs and local PWVs were calculated and compared with ones without interpolation. The means of the NRMSEs for the AWMs and local PWVs based on the MTRI method with one interpolated sample decrease from 1.14% to 0.60% and 7.48% to 4.61%, respectively. Moreover, Bland-Altman analysis and coefficient of variation were used to validate the performance of the MTRI method based on the measured local PWVs of 30 healthy subjects. In conclusion, the reconstructive interpolation for the pulse wave estimation improves the accuracy and repeatability of the carotid local PWV measurement.

FLUST: A fast, open source framework for ultrasound blood flow simulations

BACKGROUND AND OBJECTIVE

Ultrasound based blood velocity estimation is a continuously developing frontier, where the vast number of possible acquisition setups and velocity estimators makes it challenging to assess which combination is better suited for a given imaging application. FLUST, the Flow-Line based Ultrasound Simulation Tool, may be used to address this challenge, providing a common platform for evaluation of velocity estimation schemes on in silico data. However, the FLUST approach had some limitations in its original form, including reduced robustness for phase sensitive setups and the need for manual selection of integrity parameters. In addition, implementation of the technique and therefore also documentation of signal integrity was left to potential users of the approach.

METHODS

In this work, several improvements to the FLUST technique are proposed and investigated, and a robust, open source simulation framework developed. The software supports several transducer types and acquisition setups, in addition to a range of different flow phantoms. The main goal of this work is to offer a robust, computationally cheap and user-friendly framework to simulate ultrasound data from stationary blood velocity fields and thereby facilitate design and evaluation of estimation schemes, including acquisition design, velocity estimation and other post-processing steps.

RESULTS

The technical improvements proposed in this work resulted in reduced interpolation errors, reduced variability in signal power, and also automatic selection of spatial and temporal discretization parameters. Results are presented illustrating the challenges and the effectiveness of the solutions. The integrity of the improved simulation framework is validated in an extensive study, with results indicating that speckle statistics, spatial and temporal correlation and frequency content all correspond well with theoretical predictions. Finally, an illustrative example shows how FLUST may be used throughout the design and optimization process of a velocity estimator.

CONCLUSIONS

The FLUST framework is available as a part of the UltraSound ToolBox (USTB), and the results in this paper demonstrate that it can be used as an efficient and reliable tool for the development and validation of ultrasound-based velocity estimation schemes.

Two-Dimensional Wavenumber Analysis Implemented in Ultrasonic Vector Doppler Method with Focused Transmit Beams

The multi-angle Doppler method was introduced for the estimation of velocity vectors by measuring axial velocities from multiple directions. We have recently reported that the autocorrelation-based velocity vector estimation could be ameliorated significantly by estimating the wavenumbers in two dimensions. Since two-dimensional wavenumber estimation requires a snapshot of an ultrasonic field, the method was first implemented in plane wave imaging. Although plane wave imaging is predominantly useful for examining blood flows at an extremely high temporal resolution, it was reported that the contrast in a B-mode image obtained with a few plane wave emissions was lower than that obtained with focused beams. In this study, the two-dimensional wavenumber analysis was first implemented in a framework with focused transmit beams. The simulations showed that the proposed method achieved an accuracy in velocity estimation comparable to that of the method with plane wave imaging. Furthermore, the performances of the methods implemented in focused beam and plane wave imaging were compared by measuring human common carotid arteries in vivo. Image contrasts were analyzed in normal and clutter-filtered B-mode images. The method with focused beam imaging achieved a better contrast in normal B-mode imaging, and similar velocity magnitudes and angles were obtained by both the methods with focused beam and plane wave imaging. In contrast, the method with plane wave imaging gave a better contrast in a clutter-filtered B-mode image and smaller variances in velocity magnitudes than those with focused beams.

On the Investigation of Autocorrelation-Based Vector Doppler Method With Plane Wave Imaging

Although color flow imaging is one of the representative applications of the Doppler method, it can estimate only the velocity component in the direction of ultrasonic propagation, that is, the axial velocity component. The vector Doppler method with high-frame-rate plane wave imaging overcomes such a limitation by estimating the blood flow velocity vectors using the axial velocities obtained by emitting plane waves in multiple directions. The autocorrelation technique can be used for the estimation of the axial velocity using the phase shift of an ultrasonic echo signal between two transmit-receive events. The technique also requires the frequency of the received echo signal. Although the center frequency of the emitted ultrasonic signal is commonly used in the estimation of axial velocities, the center frequency should be estimated from the received signals. In this study, a method for the estimation of the center frequency designed particularly for the high-frame-rate plane wave imaging was developed. The proposed method estimates the wavenumbers of the received signal in lateral and vertical directions to estimate the wavenumber in the axial direction, from which the center frequency was estimated. The beam steering angle was also estimated from the wavenumbers in the two directions. The effect of the proposed method was validated in simulations. The absolute bias error (ABE) and root-mean squared error in estimated velocity vectors obtained by plane wave imaging with three beam steering angles (-15°, 0°, and 15°) were reduced from 9.27% and 14.80% to 1.15% and 8.75%, respectively, by the proposed method. The applicability of the proposed method to in vivo measurements was also demonstrated using the in vivo recordings of human common carotid arteries. Physiologically consistent blood flow velocity distributions were obtained with respect to three subjects using the proposed method.

Acoustic-Field Beamforming-Based Generalized Coherence Factor for Handheld Ultrasound

Handheld ultrasound devices have been widely used for diagnostic applications. The use of the acoustic-field beamforming (AFB) method has been proposed for handheld ultrasound to reduce electricity consumption and avoid battery and unwanted heat issues. However, the image quality, such as the contrast ratio and contrast-to-noise-ratio, are poorer with this technique than with the conventional delay-and-sum method. To address the problems associated with the worse image quality in AFB imaging, in this paper we propose the use of an AFB-based generalized coherence factor (GCF) technique, in which the GCF weighting developed for adaptive beamforming is extended to AFB. Simulation data, experimental results, and in vivo testing verified the efficacy of our proposed AFB-based GCF technique.

Emerging Technology: Ultrasound Vector Flow Imaging—A Novel Approach to Arterial Hemodynamic Quantification

Accurate, reliable, and easily obtainable quantification of peripheral arterial hemodynamic states has long been a holy grail of vascular ultrasound. While conventional Doppler modalities have been relied upon for decades to provide velocity, directionality, and flow volume data for integration into patient management schema, they carry limitations in accurately and reproducibly quantifying complex arterial hemodynamic patterns. Advances in ultrasound imaging architecture, such as virtual beamforming, integration of “big data” capabilities, and the use of enhanced digital signal processing methods have opened the door for a novel approach to arterial hemodynamic mapping and quantification—ultrasound vector flow imaging (VFI). This article presents an overview of the technological underpinnings of VFI, compares it with conventional pulsed wave and color Doppler methods, and describes the potential clinical benefits of this emerging vascular ultrasound modality.

Acoustic-Field Beamforming for Low-Power Portable Ultrasound

Portable ultrasound has been extensively used for diagnostic applications in health monitoring, emergency rooms, and ambulances. However, these handheld ultrasound systems may suffer from heat and battery issues attributed to the large power consumption of the transmitter. Additionally, the largest portion of the direct current (DC) power consumption can be attributed to the amplifier in the digital-to-analog converter (DAC) of the transmitter and to the analog-to-digital converter (ADC) of the receiver. Therefore, the number of transmit/receive channels in a portable ultrasound instrument is one of the crucial design factors regarding heat and battery related issues. To address these problems, we propose an acoustic-field beamforming (AFB) technique for low-power portable ultrasound systems with a single receive and five transmit channels. Finally, the simulation, experimental, and in vivo results verified the feasibility of this approach.

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.

Flow Complexity Estimation in Dysfunctional Arteriovenous Dialysis Fistulas using Vector Flow Imaging

Non-invasive assessment is preferred for monitoring arteriovenous dialysis fistulas (AVFs). Vector concentration assesses flow complexity, which may correlate with stenosis severity. We determined whether vector concentration could assess stenosis severity in dysfunctional AVFs. Vector concentration was estimated in four stenotic phantoms at different pulse repetition frequencies. Spectral Doppler peak velocity and vector concentration were measured in 12 patients with dysfunctional AVFs. Additionally, 5 patients underwent digital subtraction angiography (DSA). In phantoms, vector concentration exhibited an inverse relationship with stenosis severity and was less affected by aliasing in severe stenoses. In nine stenoses of 5 patients undergoing DSA, vector concentration correlated strongly with stenosis severity (first stenosis: r = -0.73, p = 0.04; other stenoses; r = -0.69, p = 0.02) and mid-stenotic diameter (first stenosis: r = 0.87, p = 0.006; other stenoses: r = 0.70, p = 0.02) as opposed to peak velocities (p > 0.05). Vector concentration is less affected by aliasing in severe stenoses and correlates with DSA in patients with dysfunctional AVF.

Optimization of transmission and reconstruction parameters in angular displacement compounding using plane wave ultrasound

In ultrasound elastography, plane-wave acquisitions and angular displacement compounding (ADC) are often used and combined to allow high frame rates and to improve accuracy of lateral displacement estimates, respectively. This study investigates the performance of displacement and strain estimation for ADC as a function of; the main-to-grating-lobe-amplitude ratio which decreases as a function of steering angle; plane-wave acquisition and Delay-and-Sum (DaS)-related parameters; and grating-lobe filter cut-off frequency. Three experiments were conducted with a block phantom to test ADC performance for displacement fields of varying complexity: a lateral transducer shift, phantom rotation and phantom deformation. Experiments were repeated for four linear array transducers (pitch-to-lambda ratios between 0.6 and 1.4). Best ADC performance was found for steering angles that resulted in a theoretically derived main-to-grating-lobe-amplitude ratio of 1.7 dB for pure lateral translation and 6 dB for predominately lateral strain or rotation. Temporal filtering to reduce grating lobe signal or shifting of the receive aperture to receive angles below or above the optimal angle, as dictated by the main-to-grating-lobe-amplitude ratio, did not improve results. The accuracy of lateral displacement and strain estimates was improved by apodization in transmission and a dedicated F-number in DaS (0.75) allowing incidence angles within ± 33° in the active aperture. ADC with the optimized settings as found in this study improves the accuracy of displacements and strain estimates up to 80.7% compared to non-ADC. Compared to ADC settings described in current literature, our optimization improved the accuracy by 11.9% to 75.3% for lateral displacement and strain, and by 89.3% to 96.2% for rotation. The accuracy of ADC in rotation seemed to depend highly on plane-wave and DaS-related parameters which may explain the major improvement compared to settings in current literature. The overall improvement by optimized ADC was statistically significant compared to non-ADC (p = 0.003) and literature (p = 0.002).

In Vivo Blood Velocity Vector Imaging Using Adaptive Velocity Compounding in the Carotid Artery Bifurcation

Visualization and quantification of blood flow are considered important for early detection of atherosclerosis and patient-specific diagnosis and intervention. As conventional Doppler imaging is limited to 1-D velocity estimates, 2-D and 3-D techniques are being developed. We introduce an adaptive velocity compounding technique that estimates the 2-D velocity vector field using predominantly axial displacements estimated by speckle tracking from dual-angle plane wave acquisitions. Straight-vessel experiments with a 7.8-MHz linear array transducer connected to a Verasonics Vantage ultrasound system revealed that the technique performed with a maximum velocity magnitude bias and angle bias of -3.7% (2.8% standard deviation) and -0.16° (0.41° standard deviation), respectively. In vivo, complex flow patterns were visualized in two healthy and three diseased carotid arteries and quantified using a vector complexity measure that increased with increasing wall irregularity. This measure could potentially be a relevant clinical parameter which might aid in early detection of atherosclerosis.

Spatial Angular Compounding With Affine-Model-Based Optical Flow for Improvement of Motion Estimation

Tissue motion estimation is an essential step for ultrasound elastography. Our previous study has shown that the affine-model-based optical flow (OF) method outperforms the normalized cross-correlation-based block matching (BM) method in motion estimation. However, the quality of lateral estimation using OF is still low due to inherent limitation of ultrasound imaging. BM-based spatial angular compounding (SAC) has been developed to obtain better motion estimation. In this paper, OF-based SAC (OF-SAC) is proposed to further improve the performance of lateral (and axial) estimation, and it is compared with BM-based SAC (BM-SAC). Plane wave as well as focused wave is transmitted in both simulations and phantom experiments on a linear array. In order to compare the performance quantitatively, the root-mean-square error (RMSE) of axial/lateral displacement and strain, and signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) of axial/lateral strain are used as the evaluation criteria in the simulations. In the phantom experiments, the SNR and CNR are used to assess the quality of axial/lateral strain. The results show that for both OF and BM, SAC improves the performance of motion estimation, regardless of using plane or focused wave transmission. More importantly, OF-SAC is shown to outperform BM-SAC with lower RMSE, higher SNR, and higher CNR. In addition, preliminary in vivo experiments on the carotid artery of a healthy human subject also prove the superiority of OF-SAC. These results suggest that OF-SAC is preferred for both axial and lateral motion estimation to BM-SAC.

High-Frame-Rate Contrast-Enhanced Echocardiography Using Diverging Waves: 2-D Motion Estimation and Compensation

Combining diverging ultrasound waves and microbubbles could improve contrast-enhanced echocardiography (CEE), by providing enhanced temporal resolution for cardiac function assessment over a large imaging field of view. However, current image formation techniques using coherent summation of echoes from multiple steered diverging waves (DWs) are susceptible to tissue and microbubble motion artifacts, resulting in poor image quality. In this study, we used correlation-based 2-D motion estimation to perform motion compensation for CEE using DWs. The accuracy of this motion estimation method was evaluated with Field II simulations. The root-mean-square velocity errors were 5.9% ± 0.2% and 19.5% ± 0.4% in the axial and lateral directions, when normalized to the maximum value of 62.8 cm/s which is comparable to the highest speed of blood flow in the left ventricle (LV). The effects of this method on image contrast ratio (CR) and contrast-to-noise ratio (CNR) were tested in vitro using a tissue mimicking rotating disk with a diameter of 10 cm. Compared against the control without motion compensation, a mean increase of 12 dB in CR and 7 dB in CNR were demonstrated when using this motion compensation method. The motion correction algorithm was tested in vivo on a CEE data set acquired with the Ultrasound Array Research Platform II performing coherent DW imaging. Improvement of the B-mode and contrast-mode image quality with cardiac motion and blood flow-induced microbubble motion was achieved. The results of motion estimation were further processed to interpret blood flow in the LV. This allowed for a triplex cardiac imaging technique, consisting of B mode, contrast mode, and 2-D vector flow imaging with a high frame rate of 250 Hz.

Simultaneous Vascular Strain and Blood Vector Velocity Imaging Using High-Frequency Versus Conventional-Frequency Plane Wave Ultrasound: A Phantom Study

Plaque strain and blood vector velocity imaging of stenosed arteries are expected to aid in diagnosis and prevention of cerebrovascular disease. Ultrafast plane wave imaging enables simultaneous strain and velocity estimation. Multiple ultrasound vendors are introducing high-frequency ultrasound probes and systems. This paper investigates whether the use of high-frequency ultrafast ultrasound is beneficial for assessing blood velocities and strain in arteries. The performance of strain and blood flow velocity estimation was compared between a high-frequency transducer (MS250, fc = 21 MHz) and a clinically utilized transducer (L12-5, fc = 9 MHz). Quantitative analysis based on straight tube phantom experiments revealed that the MS250 outperformed the L12-5 in the superficial region: low velocities near the wall were more accurately estimated and wall strains were better resolved. At greater than 2-cm echo depth, the L12-5 performed better due to the high attenuation of the MS250 probe. Qualitative comparison using a perfused patient-specific carotid bifurcation phantom confirmed these findings. Thus, in conclusion, for strain and blood velocity estimation for depths up to ~2 cm, a high-frequency probe is recommended.

Low-Cost 3-D Flow Estimation of Blood With Clutter

Volumetric flow rate estimation is an important ultrasound medical imaging modality that is used for diagnosing cardiovascular diseases. Flow rates are obtained by integrating velocity estimates over a cross-sectional plane. Speckle tracking is a promising approach that overcomes the angle dependency of traditional Doppler methods, but suffers from poor lateral resolution. Recent work improves lateral velocity estimation accuracy by reconstructing a synthetic lateral phase (SLP) signal. However, the estimation accuracy of such approaches is compromised by the presence of clutter. Eigen-based clutter filtering has been shown to be effective in removing the clutter signal; but it is computationally expensive, precluding its use at high volume rates. In this paper, we propose low-complexity schemes for both velocity estimation and clutter filtering. We use a two-tiered motion estimation scheme to combine the low complexity sum-of-absolute-difference and SLP methods to achieve subpixel lateral accuracy. We reduce the complexity of eigen-based clutter filtering by processing in subgroups and replacing singular value decomposition with less compute-intensive power iteration and subspace iteration methods. Finally, to improve flow rate estimation accuracy, we use kernel power weighting when integrating the velocity estimates. We evaluate our method for fast- and slow-moving clutter for beam-to-flow angles of 90° and 60° using Field II simulations, demonstrating high estimation accuracy across scenarios. For instance, for a beam-to-flow angle of 90° and fast-moving clutter, our estimation method provides a bias of -8.8% and standard deviation of 3.1% relative to the actual flow rate.