Speckle tracking for multi-dimensional flow estimation

High-Frame-Rate Ultrasound Velocimetry in the Healthy Femoral Bifurcation: A Comparative Study Against 4-D Flow Magnetic Resonance Imaging

OBJECTIVE

Local flow dynamics impact atherosclerosis yet are difficult to quantify with conventional ultrasound techniques. This study investigates the performance of ultrasound vector flow imaging (US-VFI) with and without ultrasound contrast agents in the healthy femoral bifurcation.

METHODS

High-frame-rate ultrasound data with incremental acoustic outputs were acquired in the femoral bifurcations of 20 healthy subjects before (50V) and after contrast injection (2V, 5V and 10V). 2-D blood-velocity profiles were obtained through native blood speckle tracking (BST) and contrast tracking (echo particle image velocimetry [echoPIV]). As a reference, 4-D flow magnetic resonance imaging (4-D flow MRI) was acquired. Contrast-to-background ratio and vector correlation were used to assess the quality of the US-VFI acquisitions. Spatiotemporal velocity profiles were extracted, from which peak velocities (PSV) were compared between the modalities. Furthermore, root-mean-square error analysis was performed.

RESULTS

US-VFI was successful in 99% of the cases and optimal VFI quality was established with the 10V echoPIV and BST settings. A good correspondence between 10V echoPIV and BST was found, with a mean PSV difference of -0.5 cm/s (limits of agreement: -14.1-13.2). Both US-VFI techniques compared well with 4-D flow MRI, with a mean PSV difference of 1.4 cm/s (-18.7-21.6) between 10V echoPIV and MRI, and 0.3 cm/s (-23.8-24.4) between BST and MRI. Similar complex flow patterns among all modalities were observed.

CONCLUSION

2-D blood-flow quantification of femoral bifurcation is feasible with echoPIV and BST. Both modalities showed good agreement compared to 4-D flow MRI. For the femoral tract the administration of contrast was not needed to increase the echogenicity of the blood for optimal image quality.

A GPU-Based, Real-Time Dealiasing Framework for High-Frame-Rate Vector Doppler Imaging

Vector Doppler is well regarded as a potential way of deriving flow vectors to intuitively visualize complex flow profiles, especially when it is implemented at high frame rates. However, this technique's performance is known to suffer from aliasing artifacts. There is a dire need to devise real-time dealiasing solutions for vector Doppler. In this article, we present a new methodological framework for achieving aliasing-resistant flow vector estimation at real-time throughput from precalculated Doppler frequencies. Our framework comprises a series of compute kernels that have synergized: 1) an extended least squares vector Doppler (ELS-VD) algorithm; 2) single-instruction, multiple-thread (SIMT) processing principles; and 3) implementation on a graphical processing unit (GPU). Results show that this new framework, when executed on an RTX-2080 GPU, can effectively generate aliasing-free flow vector maps using high-frame-rate imaging datasets acquired from multiple transmit-receive angle pairs in a carotid phantom imaging scenario. Over the entire cardiac cycle, the frame processing time for aliasing-resistant vector estimation was measured to be less than 16 ms, which corresponds to a minimum processing throughput of 62.5 frames/s. In a human femoral bifurcation imaging trial with fast flow (150 cm/s), our framework was found to be effective in resolving two-cycle aliasing artifacts at a minimum throughput of 53 frames/s. The framework's processing throughput was generally in the real-time range for practical combinations of ELS-VD algorithmic parameters. Overall, this work represents the first demonstration of real-time, GPU-based aliasing-resistant vector flow imaging using vector Doppler estimation principles.

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.

Advanced Fourier migration for Plane-Wave vector flow imaging

Ultrafast ultrasound imaging modalities have been studied extensively in the ultrasound community. It breaks the compromise between the frame rate and the region of interest by imaging the whole medium with wide unfocused waves. Continuously available data allow monitoring fast transient dynamics at hundreds to thousands of frames per second. This feature enables a more accurate and robust velocity estimation in vector flow imaging (VFI). On the other hand, the huge amount of data and real-time processing demands are still challenging in VFI. A solution is to provide a more efficient beamforming approach with smaller computation complexity than the conventional time-domain beamformer like delay-and-sum (DAS). Fourier-domain beamformers are shown to be more computationally efficient and can provide equally good image quality as DAS. However, previous studies generally focus on B-mode imaging. In this study, we propose a new framework for VFI which is based on two advanced Fourier migration methods, namely, slant stack migration (SSM) and ultrasound Fourier slice beamform (UFSB). By carefully modifying the beamforming parameters, we successfully apply the cross-beam technique within the Fourier beamformers. The proposed Fourier-based VFI is validated in simulation studies, in vitro, and in vivo experiments. The velocity estimation is evaluated via bias and standard deviation and the results are compared with conventional time-domain VFI using the DAS beamformer. In the simulation, the bias is 6.4%, -6.2%, and 5.7%, and the standard deviation is 4.3%, 2.4%, and 3.9% for DAS, UFSB, and SSM, respectively. In vitro studies reveal a bias of 4.5%, -5.3%, and 4.3% and a standard deviation of 3.5%, 1.3%, and 1.6% from DAS, UFSB, and SSM, respectively. The in vivo imaging of the basilic vein and femoral bifurcation also generate similar results using all three methods. With the proposed Fourier beamformers, the computation time can be shortened by up to 9 times and 14 times using UFSB and SSM.

To investigate the correlation between normal fetal biventricular myocardial function and gestational age using velocity vector imaging

Objective

The aim of this study was to evaluate the left and right ventricular segmental and global myocardial function of normal fetuses using velocity vector imaging and explore the correlation between global myocardial function parameters and gestational age.

Methods

A total of 127 normal fetuses were selected and divided into five groups according to gestational age for the measurement of their left and right ventricular segmental and global velocity, strain, and strain rate. This study also explored the change trend in the global myocardial function parameters at different gestational ages and analyzed its correlation with gestational age.

Results

The peak velocities of the biventricular segments of the normal fetuses showed a decreasing trend from the basal to the middle to the apex segment, and the differences were statistically significant (P < 0.05). However, the strain and peak strain rate between adjacent segments showed no significant differences (P > 0.05). The peak global velocity of both ventricles increased with the gestational age, and it was moderately correlated with gestational age; however, the correlation of strain and peak strain rate with gestational age was not statistically significant (P > 0.05).

Conclusion

In normal fetuses, the peak myocardial velocity of the biventricular segments showed a decreasing trend from the basal to the apical segment. The global peak myocardial velocity was linearly correlated with gestational age; however, the global strain and peak strain rate did not change as gestational age increased, indicating that the myocardial deformability of the fetus' ventricles was constant in the middle and late trimesters.

Quantification of left ventricular strain and torsion by joint analysis of 3D tagging and cine MR images

Cardiovascular magnetic resonance (CMR) imaging is the gold standard for the non-invasive assessment of left-ventricular (LV) function. Prognostic value of deformation metrics extracted directly from regular SSFP CMR images has been shown by numerous studies in the clinical setting, but with some limitations to detect torsion of the myocardium. Tagged CMR introduces trackable features in the myocardium that allow for the assessment of local myocardial deformation, including torsion; it is, however, limited in the quantification of radial strain, which is a decisive metric for assessing the contractility of the heart. In order to improve SSFP-only and tagged-only approaches, we propose to combine the advantages of both image types by fusing global shape motion obtained from SSFP images with the local deformation obtained from tagged images. To this end, tracking is first performed on SSFP images, and subsequently, the resulting motion is utilized to mask and track tagged data. Our implementation is based on a recent finite element-based motion tracking tool with mechanical regularization. Joint SSFP and tagged images registration performance is assessed based on deformation metrics including LV strain and twist using human and in-house porcine datasets. Results show that joint analysis of SSFP and 3DTAG images provides better quantification of LV strain and twist as either data source alone.

Blood speckle imaging compared with conventional Doppler ultrasound for transvalvular pressure drop estimation in an aortic flow phantom

BACKGROUND

Transvalvular pressure drops are assessed using Doppler echocardiography for the diagnosis of heart valve disease. However, this method is highly user-dependent and may overestimate transvalvular pressure drops by up to 54%. This work aimed to assess transvalvular pressure drops using velocity fields derived from blood speckle imaging (BSI), as a potential alternative to Doppler.  METHODS: A silicone 3D-printed aortic valve model, segmented from a healthy CT scan, was placed within a silicone tube. A CardioFlow 5000MR flow pump was used to circulate blood mimicking fluid to create eight different stenotic conditions. Eight PendoTech pressure sensors were embedded along the tube wall to record ground-truth pressures (10 kHz). The simplified Bernoulli equation with measured probe angle correction was used to estimate pressure drop from maximum velocity values acquired across the valve using Doppler and BSI with a GE Vivid E95 ultrasound machine and 6S-D cardiac phased array transducer.

RESULTS

There were no significant differences between pressure drops estimated by Doppler, BSI and ground-truth at the lowest stenotic condition (10.4 ± 1.76, 10.3 ± 1.63 vs. 10.5 ± 1.00 mmHg, respectively; p > 0.05). Significant differences were observed between the pressure drops estimated by the three methods at the greatest stenotic condition (26.4 ± 1.52, 14.5 ± 2.14 vs. 20.9 ± 1.92 mmHg for Doppler, BSI and ground-truth, respectively; p < 0.05). Across all conditions, Doppler overestimated pressure drop (Bias = 3.92 mmHg), while BSI underestimated pressure drop (Bias = -3.31 mmHg).

CONCLUSIONS

BSI accurately estimated pressure drops only up to 10.5 mmHg in controlled phantom conditions of low stenotic burden. Doppler overestimated pressure drops of 20.9 mmHg. Although BSI offers a number of theoretical advantages to conventional Doppler echocardiography, further refinements and clinical studies are required with BSI before it can be used to improve transvalvular pressure drop estimation in the clinical evaluation of aortic stenosis.

Improving image contrast and accuracy in velocity estimation by convolution filters for intracardiac blood flow imaging

In this study, the point spread function (PSF) of an ultrasound imaging system was estimated and used as a reference signal in a filtering method for improvement of image quality. The PSF of the imaging system was estimated from measured echo signals from an imaging target. Convolution filters (including deconvolution) were used for improvement of image contrast and spatial resolution. Furthermore, the accuracy in estimation of velocity vectors was evaluated for investigation of the impact of the proposed filters on velocity estimation. In the phantom experiment, contrast of the B-mode image was improved from 76.4 dB to 81.1 dB and 77.8 dB using the convolution and deconvolution filters, respectively. Also, the two-dimensional (2D) velocity distribution in the phantom was estimated by the block matching method, and the bias error (BE) in the estimated lateral velocity was reduced from -19.7% to 2.16% and 2.29% using the convolution and the deconvolution filters, respectively.

On the Depth-Dependent Accuracy of Plane-Wave-Based Vector Velocity Measurements With Linear Arrays

High-frame-rate vector Doppler methods are used to measure blood velocities over large 2-D regions, but their accuracy is often estimated over a short range of depths. This article thoroughly examines the dependence of velocity measurement accuracy on the target position. Simulations were carried out on flat and parabolic flow profiles, for different Doppler angles, and considering a 2-D vector flow imaging (2-D VFI) method based on plane wave transmission and speckle tracking. The results were also compared with those obtained by the reference spectral Doppler (SD) method. Although, as expected, the bias and standard deviation generally tend to worsen at increasing depths, the measurements also show the following. First, the errors are much lower for the flat profile (from ≈ -4 ± 3% at 20 mm to ≈ -17 ± 4% at 100 mm) than for the parabolic profile (from ≈ -4 ± 3% to ≈ -38 ±%). Second, only part of the relative estimation error is related to the inherent low resolution of the 2-D VFI method. For example, even for SD, the error bias increases (on average) from -0.7% (20 mm) to -17% (60 mm) up to -26% (100 mm). Third, conversely, the beam divergence associated with the linear array acoustic lens was found to have a great impact on the velocity measurements. By simply removing such lens, the average bias for 2-D VFI at 60 and 100 mm dropped down to -9.4% and -19.4%, respectively. In conclusion, the results indicate that the transmission beam broadening on the elevation plane, which is not limited by reception dynamic focusing, is the main cause of velocity underestimation in the presence of high spatial gradients.

Tensor Velocity Imaging With Motion Correction

This article presents a motion compensation procedure that significantly improves the accuracy of synthetic aperture tensor velocity estimates for row-column arrays. The proposed motion compensation scheme reduces motion effects by moving the image coordinates with the velocity field during summation of low-resolution volumes. The velocity field is estimated using a transverse oscillation cross-correlation estimator, and each image coordinate's local tensor velocity is determined by upsampling the field using spline interpolation. The motion compensation procedure is validated using Field II simulations and flow measurements acquired using a 3-MHz row-column addressed probe and the research scanner SARUS. For a peak velocity of 25 cm/s, a pulse repetition frequency of 2 kHz, and a beam-to-flow angle of 60°, the proposed motion compensation procedure was able to reduce the relative bias from -27.0% to -9.4% and the standard deviation from 8.6% to 8.1%. In simulations performed with a pulse repetition frequency of 10 kHz, the proposed method reduces the bias in all cases with beam-to-flow angles of 60° and 75° and peak velocities between 10 and 150 cm/s.

CNN-Based Ultrasound Image Reconstruction for Ultrafast Displacement Tracking

Thanks to its capability of acquiring full-view frames at multiple kilohertz, ultrafast ultrasound imaging unlocked the analysis of rapidly changing physical phenomena in the human body, with pioneering applications such as ultrasensitive flow imaging in the cardiovascular system or shear-wave elastography. The accuracy achievable with these motion estimation techniques is strongly contingent upon two contradictory requirements: a high quality of consecutive frames and a high frame rate. Indeed, the image quality can usually be improved by increasing the number of steered ultrafast acquisitions, but at the expense of a reduced frame rate and possible motion artifacts. To achieve accurate motion estimation at uncompromised frame rates and immune to motion artifacts, the proposed approach relies on single ultrafast acquisitions to reconstruct high-quality frames and on only two consecutive frames to obtain 2-D displacement estimates. To this end, we deployed a convolutional neural network-based image reconstruction method combined with a speckle tracking algorithm based on cross-correlation. Numerical and in vivo experiments, conducted in the context of plane-wave imaging, demonstrate that the proposed approach is capable of estimating displacements in regions where the presence of side lobe and grating lobe artifacts prevents any displacement estimation with a state-of-the-art technique that relies on conventional delay-and-sum beamforming. The proposed approach may therefore unlock the full potential of ultrafast ultrasound, in applications such as ultrasensitive cardiovascular motion and flow analysis or shear-wave elastography.

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.

Use of Speckle Tracking Echocardiography to Detect Induced Regional Strain Changes in the Murine Myocardium by Acoustic Radiation Force

BACKGROUND

It is difficult to simulate the abnormal myocardial strain patterns caused by ischemic coronary artery disease (CAD) which are a precursor to heart failure (HF) within an animal model. Simulation of these strain changes could contribute to better understanding of the early formative stages of HF. This is especially important in investigating the poorly understood pathogenesis of heart failure with preserved ejection fraction (HFpEF). Here, we discuss delivery of high intensity focused ultrasound (HIFU) in a murine model to alter left ventricular (LV) regional longitudinal strain (RLS), and use of speckle tracking echocardiography to detect these changes.

METHODS

HIFU pulses (pressure amplitude 1.7 MPa) were generated by amplifying a sinusoidal waveform from a function generator into a piezoelectric transducer. These pulses were then directed extracorporeally towards the anterior LV surface of C57BI6 mice during three time periods (early, mid, and late diastole). Speckle tracking echocardiography was then used to quantify changes in RLS within six segments of the LV.

RESULTS

We observed an increase in LV RLS with acoustic augmentation during all three time periods. This augmentation was most prominent near the anterior apical region in early diastole and near the posterior basilar region during late diastole.

CONCLUSIONS

Our findings demonstrate the application of HIFU to non-invasively induce changes in RLS within a murine model. Our results also reflect the capability of speckle tracking echocardiography to analyze and quantify these changes. These findings represent the first demonstration of ultrasound-induced augmentation in LV RLS within a small animal model.

In vitro performance of echoPIV for assessment of laminar flow profiles in a carotid artery stent

Purpose: Detailed blood flow studies may contribute to improvements in carotid artery stenting. High-frame-rate contrast-enhanced ultrasound followed by particle image velocimetry (PIV), also called echoPIV, is a technique to study blood flow patterns in detail. The performance of echoPIV in presence of a stent has not yet been studied extensively. We compared the performance of echoPIV in stented and nonstented regions in an in vitro flow setup.

Approach: A carotid artery stent was deployed in a vessel-mimicking phantom. High-frame-rate contrast-enhanced ultrasound images were acquired with various settings. Signal intensities of the contrast agent, velocity values, and flow profiles were calculated.

Results: The results showed decreased signal intensities and correlation coefficients inside the stent, however, PIV analysis in the stent still resulted in plausible flow vectors.

Conclusions: Velocity values and laminar flow profiles can be measured in vitro in stented arteries using echoPIV.

A Deep Learning Approach to Resolve Aliasing Artifacts in Ultrasound Color Flow Imaging

Despite being used clinically as a noninvasive flow visualization tool, color flow imaging (CFI) is known to be prone to aliasing artifacts that arise due to fast blood flow beyond the detectable limit. From a visualization standpoint, these aliasing artifacts obscure proper interpretation of flow patterns in the image view. Current solutions for resolving aliasing artifacts are typically not robust against issues such as double aliasing. In this article, we present a new dealiasing technique based on deep learning principles to resolve CFI aliasing artifacts that arise from single- and double-aliasing scenarios. It works by first using two convolutional neural networks (CNNs) to identify and segment CFI pixel positions with aliasing artifacts, and then it performs phase unwrapping at these aliased pixel positions. The CNN for aliasing identification was devised as a U-net architecture, and it was trained with in vivo CFI frames acquired from the femoral bifurcation that had known presence of single- and double-aliasing artifacts. Results show that the segmentation of aliased CFI pixels was achieved successfully with intersection over union approaching 90%. After resolving these artifacts, the dealiased CFI frames consistently rendered the femoral bifurcation's triphasic flow dynamics over a cardiac cycle. For dealiased CFI pixels, their root-mean-squared difference was 2.51% or less compared with manual dealiasing. Overall, the proposed dealiasing framework can extend the maximum flow detection limit by fivefold, thereby improving CFI's flow visualization performance.

A Novel 2-D Speckle Tracking Method for High-Frame-Rate Echocardiography

Speckle tracking echocardiography (STE) is a clinical tool to noninvasively assess regional myocardial function through the quantification of regional motion and deformation. Even if the time resolution of STE can be improved by high-frame-rate (HFR) imaging, dedicated HFR STE algorithms have to be developed to detect very small interframe motions. Therefore, in this article, we propose a novel 2-D STE method, purposely developed for HFR echocardiography. The 2-D motion estimator consists of a two-step algorithm based on the 1-D cross correlations to separately estimate the axial and lateral displacements. The method was first optimized and validated on simulated data giving an accuracy of ~3.3% and ~10.5% for the axial and lateral estimates, respectively. Then, it was preliminarily tested in vivo on ten healthy volunteers showing its clinical applicability and feasibility. Moreover, the extracted clinical markers were in the same range as those reported in the literature. Also, the estimated peak global longitudinal strain was compared with that measured with a clinical scanner showing good correlation and negligible differences (-20.94% versus -20.31%, p -value = 0.44). In conclusion, a novel algorithm for STE was developed: the radio frequency (RF) signals were preferred for the axial motion estimation, while envelope data were preferred for the lateral motion. Furthermore, using 2-D kernels, even for 1-D cross correlation, makes the method less sensitive to noise.

Blind Source Separation for Clutter and Noise Suppression in Ultrasound Imaging: Review for Different Applications

Blind source separation (BSS) refers to a number of signal processing techniques that decompose a signal into several "source" signals. In recent years, BSS is increasingly employed for the suppression of clutter and noise in ultrasonic imaging. In particular, its ability to separate sources based on measures of independence rather than their temporal or spatial frequency content makes BSS a powerful filtering tool for data in which the desired and undesired signals overlap in the spectral domain. The purpose of this work was to review the existing BSS methods and their potential in ultrasound imaging. Furthermore, we tested and compared the effectiveness of these techniques in the field of contrast-ultrasound super-resolution, contrast quantification, and speckle tracking. For all applications, this was done in silico, in vitro, and in vivo. We found that the critical step in BSS filtering is the identification of components containing the desired signal and highlighted the value of a priori domain knowledge to define effective criteria for signal component selection.

Data-Adaptive 2-D Tracking Doppler for High-Resolution Spectral Estimation

Spectral broadening in pulsed-wave Doppler caused by the transit-time effect deteriorates the frequency resolution and may cause overestimation of maximum velocities in high-velocity blood flow regions and for large beam-to-flow angles. Data-adaptive spectral estimators have been shown to provide improved frequency resolution, especially for small ensemble lengths, but offer little or no improvement when the transit-time effect dominates. In this work, a method is presented that combines a data-adaptive spectral estimation method, the power spectral Capon, and 2-D tracking Doppler to enable improved frequency resolution for both high and low velocities. For each velocity, a time signal is extracted by tracking scatterers over time and space to decrease the transit-time effect, and power spectral Capon is used for spectral estimation. The method is evaluated using simulations, flow phantom recordings, and recordings from healthy and stenotic carotid arteries. Simulation results showed that the spectral width was decreased by 60% compared to 2-D tracking Doppler for velocities around 2.3 m/s using 12 time samples. The reduction was estimated to be 66% using the flow phantom results for 0.85-m/s mean velocity. A 5-dB SNR gain was observed from the in vivo results compared with Welch's method. Computer simulations confirm that in the presence of velocity gradients or out-of-plane motion, the proposed method can be used to reduce spectral broadening by requiring shorter observation windows.

Combining Automatic Angle Correction and 3-D Tracking Doppler for the Assessment of Aortic Stenosis Severity

Aortic valve stenosis (AS) is a narrowing of the aortic valve opening, which causes increased load on the left ventricle. Untreated, this condition can eventually lead to heart failure and death. According to current recommendations, an accurate diagnosis of AS mandates the use of multiple acoustic windows to determine the highest velocity. Furthermore, the optimal positioning of both patient and transducer to reduce the beam-to-flow angle is emphasized. Being operator dependent, the beam alignment is a potential source of uncertainty. In this work, we perform noncompounded 3-D plane wave imaging for retrospective estimation of maximum velocities in aortic jets with automatic angle correction. This is achieved by combining a hybrid 3-D speckle tracking method to estimate the jet direction and 3-D tracking Doppler to generate angle-corrected sonograms, using the direction from speckle tracking as input. Results from simulations of flow through an orifice show that 3-D speckle tracking can estimate the jet orientation with acceptable accuracy for signal-to-noise ratios above 10 dB. Results from 12 subjects show that sonograms recorded from a standard apical view using the proposed method yield a maximum velocity that matches continuous wave (CW) Doppler sonograms recorded from the acoustic window with the lowest angle within a ±10% margin, provided that a high enough pulse repetition frequency could be achieved. These results motivate further validation and optimization studies.

Imaging modalities to diagnose carotid artery stenosis: progress and prospect

In the past few decades, imaging has been developed to a high level of sophistication. Improvements from one-dimension (1D) to 2D images, and from 2D images to 3D models, have revolutionized the field of imaging. This not only helps in diagnosing various critical and fatal diseases in the early stages but also contributes to making informed clinical decisions on the follow-up treatment profile. Carotid artery stenosis (CAS) may potentially cause debilitating stroke, and its accurate early detection is therefore important. In this paper, the technical development of various CAS diagnosis imaging modalities and its impact on the clinical efficacy is thoroughly reviewed. These imaging modalities include duplex ultrasound (DUS), computed tomography angiography (CTA) and magnetic resonance angiography (MRA). For each of the imaging modalities considered, imaging methodology (principle), critical imaging parameters, and the extent of imaging the vulnerable plaque are discussed. DUS is usually the initial recommended CAS diagnostic examination. However, for the therapeutic intervention, either MRA or CTA is recommended for confirmation, and for added information on intracranial cerebral circulation and aortic arch condition for procedural planning. Over the past few decades, the focus of CAS diagnosis has also shifted from pure stenosis quantification to plaque characterization. This has led to further advancement in the existing imaging tools and development of other potential imaging tools like Optical coherence tomography (OCT), photoacoustic tomography (PAT), and infrared (IR) thermography.

Vector Flow Imaging of a Highly Laden Suspension in a Zinc-Air Flow Battery Model

Flow batteries using suspension electrodes, e.g., zinc-air flow batteries (ZABs), have recently gained renewed interest as potential candidates for grid energy storage or mobile applications. The performance of ZABs depends on the local flow conditions of the suspension in the electrochemical cell, which acts as an electrode. Hence, it is crucial to measure and understand the complex flow characteristics of such solid-liquid suspensions. The investigated suspension electrode is an opaque slurry that consists of microscopic zinc particles and an aqueous potassium hydroxide electrolyte. Commonly, ultrasound Doppler velocimetry is used for flow imaging in opaque fluids. However, due to the high particle concentration in the suspension electrode, strong scattering and wavefront distortions of the ultrasound are introduced. In this paper, we show that this results in an increased measurement uncertainty for Doppler-based velocity estimation. Instead, ultrasound image velocimetry is applied to measure the 2-D and two-component flow field in the zinc-electrolyte suspension. This is possible by adapting the measurement system to the suspension with a calibration setup. The total measurement uncertainties of 4.1% and 2.5% for the axial and lateral flow components are derived from the calibration measurements. For the first time, the flow field of such a suspension could be measured in a scaled fluidic model of a ZAB. The comparison of the estimated flow rates from the velocity profiles showed good agreement to a gravimetric reference. A significant difference in the flow characteristics of a macroscopically homogeneous electrolyte and the same electrolyte loaded with 8 vol.-% zinc particles, i.e., the suspension electrode, was found. Along with the demonstration of the measurement technique for opaque, concentrated suspensions, the measurement data will be used to calibrate and validate numerical models for comparable multiphase fluids.

Fast Flow-Line-Based Analysis of Ultrasound Spectral and Vector Velocity Estimators

A new technique, termed FLUST (FlowLine Ultrasound Simulation Tool), is proposed as a computationally cheap alternative to simulations based on randomly positioned scatterers for the simulation of stationary blood velocity fields. In FLUST, the flow field is represented as a collection of flow lines. Point spread functions are first calculated at regularly spaced positions along the flow lines before realizations of single scatterers traversing the flow lines are generated using temporal interpolation. Several flow-line realizations are then generated by convolution with temporal noise filters, and finally, flow-field realizations are obtained by the summation of the individual flow-line realizations. Flow-field realizations produced by FLUST are shown to correspond well with conventional Field II simulations both quantitatively and qualitatively. The added value of FLUST is demonstrated by using the proposed simulation technique to obtain multiple realizations of realistic 3-D flow fields at a significantly reduced computational cost. This information is utilized for a performance assessment of different spectral and vector velocity estimators for carotid and coronary imaging applications. The computational load of FLUST does not increase substantially with the number of realizations or simulated frames, and for the examples shown, it is the fastest alternative when the total number of simulated frames exceeds 48. In the examples, the standard deviation and bias of the velocity estimators are calculated using 100 FLUST realizations, in which case the proposed method is two orders of magnitude faster than simulations based on random scatterer positions.

4-D Intracardiac Ultrasound Vector Flow Imaging-Feasibility and Comparison to Phase-Contrast MRI

In vivo characterization of intracardiac blood velocity vector fields may provide new clinical information but is currently not available for bedside evaluation. In this paper, 4-D vector flow imaging for intracardiac flow assessment is demonstrated using a clinical ultrasound (US) system and a matrix array transducer, without the use of contrast agent. Two acquisition schemes were developed, one for full volumetric coverage of the left ventricle (LA) at 50 vps and a 3-D thick-slice setup with continuous frame acquisition (4000 vps), both utilizing ECG-gating. The 3-D vector velocity estimates were obtained using a novel method combining phase and envelope information. In vitro validation in a rotating tissue-mimicking phantom revealed velocity estimates in compliance with the ground truth, with a linear regression slope of 0.80, 0.77, and 1.03 for the , , and velocity components, and with standard deviations of 2.53, 3.19, and 0.95 cm/s, respectively. In vivo measurements in a healthy LV showed good agreement with PC-MRI. Quantitative analysis of energy loss (EL) and kinetic energy (KE) further showed similar trends, with peak KE at 1.5 and 2.4 mJ during systole and 3.6 and 3.1 mJ for diastole for US and PC-MRI. Similar for EL, 0.15- 0.2 and 0.7 mW was found during systole and 0.6 and 0.7 mW during diastole, for US and PC-MRI, respectively. Overall, a potential for US as a future modality for 4D cardiac vector flow imaging was demonstrated, which will be further evaluated in clinical studies.

High-Frame Rate Vector Flow Imaging of the Carotid Bifurcation in Healthy Adults: Comparison With Color Doppler Imaging

OBJECTIVES

To evaluate the carotid bifurcation in healthy adults using a commercial system equipped with high-frame rate vector flow imaging (VFI) based on the plane wave and to compare VFI with color Doppler imaging.

METHODS

Carotid bifurcation diameters and flow characteristics of 60 vessels in 60 healthy volunteers were evaluated quantitatively and qualitatively to assess complex flow patterns and their extension and duration.

RESULTS

Complex flow in the internal carotid artery (ICA) was associated with a statistically significant difference in the ΔICA sinus-to-common carotid artery (CCA) diameter ratio (the relative change in diameter between the CCA and ICA sinus.) Vector flow imaging and color Doppler imaging were in accordance when detecting complex flow in 96.7% of cases; in 3.3% of cases, only VFI identified small recirculation areas of short duration. Vector flow imaging highlighted a larger extension of the complex flow (mean ± SD, 47.7 ± 28.5 mm² ; median, 45.5 mm² ) compared with color Doppler imaging (mean, 29.2 ± 19.9 mm² ; median, 29.5 mm² ) and better depicted different complex flow patterns; a strong correlation (r = 0.84) was found between the ΔICA sinus-to-CCA diameter ratio and the complex flow extension. Vector flow imaging showed a longer duration of the flow disturbances (mean, 380 ± 218 milliseconds; median, 352.5 milliseconds) compared with color Doppler imaging (mean, 325 ± 206 milliseconds; median, 333 milliseconds), and there was a strong correlation (r = 0.92).

CONCLUSIONS

Vector flow imaging is as effective as color Doppler imaging in the detection of flow disturbances, but it is more powerful in the assessment of complex flow patterns.

Role of the CHADS2 Score in the Evaluation of Carotid Atherosclerosis in Patients with Atrial Fibrillation Undergoing Carotid Artery Ultrasonography

Objective

This study investigated the characteristics of carotid atherosclerosis in patients with atrial fibrillation (AF) and determined the feasibility and significance of the CHADS₂ score in predicting the degree of carotid atherosclerosis.

Methods

Consecutive patients (n = 109) with nonvalvular AF were registered and classified into two groups, the paroxysmal AF group (n = 59) and persistent AF group (n = 50). Fifty healthy patients, matched by sex and age, were considered the control group. All patients were examined using carotid ultrasound and velocity vector imaging (VVI).

Results

Compared with the control group, the mean intimal-medial thickness in the paroxysmal AF group (0.56 ± 0.11 versus 0.61 ± 0.10, respectively, P < 0.05) and the persistent AF group (0.56 ± 0.11 versus 0.64 ± 0.13, respectively, P < 0.001) was significantly increased. The plaque index (PI) in the persistent AF group was significantly higher than that observed in the paroxysmal AF group (1.05 ± 1.33 versus 1.42 ± 1.47, respectively, P < 0.001). Regarding the VVI indices, those reflecting the long-axis longitudinal motion function of carotid arteries were significantly decreased in both AF groups. Compared with the control group, a significantly lower total longitudinal displacement (tLoD) index was observed in the persistent AF group (0.73 ± 0.66 versus 0.31 ± 0.23, respectively, P < 0·0001) and the paroxysmal AF group (0.73 ± 0.66 versus 0.34 ± 0.17, P < 0·0001). The CHADS₂ score was related to indicators reflecting the structure and function of the carotid artery.

Conclusions

Carotid arterial structure and function were significantly altered in patients with AF. The degree of carotid atherosclerosis depended on the duration of AF. The CHADS₂ score may be useful as a predictor of the extent of carotid atherosclerosis in patients with AF.

Real-Time 2-D Phased Array Vector Flow Imaging

Echocardiography examination of the blood flow is currently either restricted to 1-D techniques in real-time or experimental offline 2-D methods. This paper presents an implementation of transverse oscillation for real-time 2-D vector flow imaging (VFI) on a commercial BK Ultrasound scanner. A large field-of-view (FOV) sequence for studying flow dynamics at 11 frames per second (fps) and a sequence for studying peak systolic velocities (PSVs) with a narrow FOV at 36 fps were validated. The VFI sequences were validated in a flow rig with continuous laminar parabolic flow and in a pulsating flow pump system before being tested in vivo, where measurements were obtained on two healthy volunteers. Mean PSV from 11 cycles was 155 cms^-1 with a precision of ±9.0% for the pulsating flow pump. In vivo, PSV estimated in the ascending aorta was 135 cms^-1 ± 16.9% for eight cardiac cycles. Furthermore, in vivo flow dynamics of the left ventricle and in the ascending aorta were visualized. In conclusion, angle independent 2-D VFI on a phased array has been implemented in real time, and it is capable of providing quantitative and qualitative flow evaluations of both the complex and fully transverse flow.

Combined 2-D Vector Velocity Imaging and Tracking Doppler for Improved Vascular Blood Velocity Quantification

Measurement of the maximum blood flow velocity is the primary means for determining the degree of carotid stenosis using ultrasound. The current standard for estimating the maximum velocity is pulsed-wave Doppler with manual angle correction, which is prone to error and interobserver variability. In addition, spectral broadening in the velocity spectra leads to overestimation of maximal velocities. In this paper, we propose to combine two velocity estimation methods to reduce the bias and variability in maximum velocity measurements. First, the direction of the blood flow is estimated using an aliasing-resistant least squares vector Doppler technique. Then, tracking Doppler is performed on the same data, using the direction of the vector Doppler estimate as the tracking direction. Simulations show that the method can estimate a maximum velocity of 2 m/s with accuracy 5% for beam-to-flow angles between 20° and 75°, and that the primary source of error is inaccuracy in the flow direction estimate from vector Doppler. Simulations of complex flow in a carotid bifurcation demonstrated that the combined technique provided spectral velocity profiles corresponding well with the true maximum velocity trace, and that the bias originating from the directional estimate was within 5% for all spatial points. A healthy volunteer and a volunteer with carotid artery stenosis were imaged, showing in vivo feasibility of the method, for high velocities and with beam-to-flow angles varying throughout the cardiac cycle.

Vector flow imaging techniques: An innovative ultrasonographic technique for the study of blood flow

Doppler ultrasonography is routinely used to identify abnormal blood flow. Nevertheless, conventional Doppler can be used to determine only the axial component of blood flow velocity and is angle dependent. A new method of multidimensional angle-independent estimation of flow velocity, called Vector Flow Imaging (VFI), has been proposed. It quantitatively evaluates the true velocity vector's amplitude and direction at any location into a vessel and displays a more intuitive depiction of the flow movements. High frame rate VFI, based on plane wave imaging, allows a detailed dynamic visualization of complex flow by showing even transient events, otherwise undetectable. © 2017 Wiley Periodicals, Inc. J Clin Ultrasound 45:582-588, 2017.

In Vivo Intracardiac Vector Flow Imaging Using Phased Array Transducers for Pediatric Cardiology

Two-dimensional blood speckle tracking (ST) has shown promise for measuring complex flow patterns in neonatal hearts using linear arrays and high-frame-rate plane wave imaging. For general pediatric applications, however, the need for phased array probes emerges due to the limited intercostal acoustic window available. In this paper, a clinically approved real-time duplex imaging setup with phased array probes was used to investigate the potential of blood ST for the 2-D vector flow imaging of children with congenital heart disease. To investigate transmit beam pattern and tracking accuracy, straight tubes with parabolic flow were simulated at three depths (4.5, 7, and 9.5 cm). Due to the small aperture available, diffraction effects could be observed when approaching 10 cm, which limited the number of parallel receive beams that could be utilized. Moving to (slightly) diverging beams was shown to solve this issue at the expense of a loss in signal-to-noise ratio. To achieve consistent estimates, a forward-backward tracking scheme was introduced to avoid measurement bias occurring due to tracking kernel averaging artifacts at flow domain boundaries. Promising results were observed for depths <10 cm in two pediatric patients, where complex cardiac flow patterns could be estimated and visualized. As a loss in penetration compared with color flow imaging is expected, a larger clinical study is needed to establish the clinical feasibility of this approach.

High-frame rate vector flow imaging of the carotid bifurcation

Abstract

Carotid artery atherosclerotic disease is still a significant cause of cerebrovascular morbidity and mortality. A new angle-independent technique, measuring and visualizing blood flow velocities in all directions, called vector flow imaging (VFI) is becoming available from several vendors. VFI can provide more intuitive and quantitative imaging of vortex formation, which is not clearly distinguishable in the color Doppler image. VFI, as quantitative method assessing disturbed flow patterns of the carotid bifurcation, has the potential to allow better understanding of the diagnostic value of complex flow and to enhance risk stratification. This pictorial review article will show which new information VFI adds for the knowledge of hemodynamics in comparison to the conventional ultrasound techniques.

Teaching points

• VFI is an angle-independent technique measuring flow velocities in all directions.

• This kind of VFI is based on a plane wave multidirectional excitation technique.

• VFI allows quantitative assessment of carotid streamlines progression and visualizes vorticity.

• VFI does not allow a precise comprehension of streamlines’ 3D shape.

• VFI allows a better understanding of carotid artery complex flows.

Electronic supplementary material

The online version of this article (doi:10.1007/s13244-017-0554-5) contains supplementary material, which is available to authorized users.

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.

Management and control of isotonic contraction generated stress: evaluation of masseter muscle deformation pattern by means of ecography

PURPOSE

The objective of the following study is to observe the behavior of the six layers of the masseter during an isometric contraction at maximum exertion with the deformation pattern analysis method.

MATERIALS AND METHODS

This study has been conducted by use of an ultrasound machine (MicrUs ext-1H Telemed Medical Systems Milano) and a linear probe (L12-5l40S-3 5-12 MHz 40 mm) which allowed us to record a video (DCM) comprised of 45 frames per second. The probe was fixed to a brace and the patient was asked to clench their teeth as hard as possible, obtain the muscle's maximum exertion, for 5 seconds three times, with 30 seconds intervals in between. Both right and left masseter muscles were analyzed. Then we applied to the resulting video a software (Mudy 1.7.7.2 AMID Sulmona Italy) for the analysis of muscle deformation patterns (contraction, dilatation, cross-plane, vertical strain, horizontal strain, vertical shear, horizontal shear, horizontal displacement, vertical displacement). The number of videos of masseter muscles in contraction at maximum exertion due to dental clenching made during this research is around 12,000. Out of these we chose 1,200 videos which examine 200 patients (100 females, 100 males).

RESULTS

The analysis of the deformation patterns of the masseter allows us to observe how the six layers of the muscle have different and specific functions each, which vary depending on the applied force (application point, magnitude and direction) so that we find it impossible to assign to one of the three sections of the muscle a mechanical predominance. Therefore it appears that the three parts of the muscle have specific and synergistic tasks.

Ultrasound and analysis of the deformation patterns of the masseter muscle: comparing surgical anatomy, ultrasound and functional anatomy

PURPOSE

We have tried to demonstrate whether the analysis of the muscle strain allows us to identify the three distinct functional areas of the architecture of the masseter, as one would see them by performing or viewing an anatomical dissection of said muscle, and whether these sections have behave differently in terms of origin and coping of the strain they face (quantitative analysis).

MATERIALS AND METHODS

This work has been elaborated by the use of an ultrasound machine (MicrUs ext-1H Telemed Medical Systems Milano) and a linear probe (L12-5l40S-3 5-12 MHz 40 mm) which allowed us to record a 45 frame per second video (DCM). Videos has been elaborated by use of an ultrasound machine (MicrUs ext-1H Telemed Medical Systems Milano) and a linear probe (L12-5l40S-3 5-12 MHz 40 mm) which allowed us to record a 45 frame per second video (DCM). We applied to the resulting video a software (Mudy 1.7.7.2 AMID Sulmona Italy) for the analysis of muscle deformation patters (contraction, dilatation, cross-plane, vertical strain, horizontal strain, vertical shear, horizontal shear, horizontal displacement, vertical displacement). The number of videos of masseter muscles in contraction at maximum exertion due to dental clenching made during this research is around 12,000. Out of these we chose 1,200 videos which examine 200 patients (100 females, 100 males).

RESULTS

The deformation pattern analysis of the skeletal muscle on ultrasound basis seems to be an adequate instrument to use during the investigation of the functional structure of the masseter muscle given its ability to highlight the distinct activity of each separate part of the muscle.

CONCLUSIONS

Moreover the strain does not apply to the muscle uniformly; instead it varies according to the observed area.

Ultrasound imaging velocimetry with interleaved images for improved pulsatile arterial flow measurements: a new correction method, experimental and in vivo validation

Blood velocity measurements are important in physiological science and clinical diagnosis. Doppler ultrasound is the most commonly used method but can only measure one velocity component. Ultrasound imaging velocimetry (UIV) is a promising technique capable of measuring two velocity components; however, there is a limit on the maximum velocity that can be measured with conventional hardware which results from the way images are acquired by sweeping the ultrasound beam across the field of view. Interleaved UIV is an extension of UIV in which two image frames are acquired concurrently, allowing the effective interframe separation time to be reduced and therefore increasing the maximum velocity that can be measured. The sweeping of the ultrasound beam across the image results in a systematic error which must be corrected: in this work, we derived and implemented a new velocity correction method which accounts for acceleration of the scatterers. We then, for the first time, assessed the performance of interleaved UIV for measuring pulsatile arterial velocities by measuring flows in phantoms and in vivo and comparing the results with spectral Doppler ultrasound and transit-time flow probe data. The velocity and flow rate in the phantom agreed within 5-10% of peak velocity, and 2-9% of peak flow, respectively, and in vivo the velocity difference was 9% of peak velocity. The maximum velocity measured was 1.8 m s^-1, the highest velocity reported with UIV. This will allow flows in diseased arteries to be investigated and so has the potential to increase diagnostic accuracy and enable new vascular research.

Assessing the Performance of Ultrafast Vector Flow Imaging in the Neonatal Heart via Multiphysics Modeling and In Vitro Experiments

Ultrafast vector flow imaging would benefit newborn patients with congenital heart disorders, but still requires thorough validation before translation to clinical practice. This paper investigates 2-D speckle tracking (ST) of intraventricular blood flow in neonates when transmitting diverging waves at ultrafast frame rate. Computational and in vitro studies enabled us to quantify the performance and identify artifacts related to the flow and the imaging sequence. First, synthetic ultrasound images of a neonate's left ventricular flow pattern were obtained with the ultrasound simulator Field II by propagating point scatterers according to 3-D intraventricular flow fields obtained with computational fluid dynamics (CFD). Noncompounded diverging waves (opening angle of 60°) were transmitted at a pulse repetition frequency of 9 kHz. ST of the B-mode data provided 2-D flow estimates at 180 Hz, which were compared with the CFD flow field. We demonstrated that the diastolic inflow jet showed a strong bias in the lateral velocity estimates at the edges of the jet, as confirmed by additional in vitro tests on a jet flow phantom. Furthermore, ST performance was highly dependent on the cardiac phase with low flows (<5 cm/s), high spatial flow gradients, and out-of-plane flow as deteriorating factors. Despite the observed artifacts, a good overall performance of 2-D ST was obtained with a median magnitude underestimation and angular deviation of, respectively, 28% and 13.5° during systole and 16% and 10.5° during diastole.

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

Conventional color Doppler imaging is limited, since it only provides velocity estimates along the ultrasound beam direction for a restricted field of view at a limited frame rate. High-frame-rate speckle tracking, using plane wave transmits, has shown potential for 2-D blood velocity estimation. However, due to the lack of focusing in transmit, image quality gets reduced, which hampers speckle tracking. Although ultrafast imaging facilitates improved clutter filtering, it still remains a major challenge in blood velocity estimation. Signal dropouts and poor velocity estimates are still present for high beam-to-flow angles and low blood flow velocities. In this paper, ultrafast plane wave imaging was combined with multiscale speckle tracking to assess the 2-D blood velocity vector in a common carotid artery (CCA) flow field. A multiangled plane wave imaging sequence was used to compare the performance of displacement compounding, coherent compounding, and compound speckle tracking. Zero-degree plane wave imaging was also evaluated. The performance of the methods was evaluated before and after clutter filtering for the large range of velocities (0-1.5 m/s) that are normally present in a healthy CCA during the cardiac cycle. An extensive simulation study was performed, based on a sophisticated model of the CCA, to investigate and evaluate the performance of the methods at different pulse repetition frequencies and signal-to-noise levels. In vivo data were acquired of a healthy carotid artery bifurcation to support the simulation results. In general, methods utilizing compounding after speckle tracking, i.e., displacement compounding and compound speckle tracking, were least affected by clutter filtering and provided the most robust and accurate estimates for the entire velocity range. Displacement compounding, which uses solely axial information to estimate the velocity vector, provided most accurate velocity estimates, although it required sufficiently high pulse repetition frequencies in high blood velocity phases and reliable estimates for all acquisition angles. When this latter requirement was not met, compound speckle tracking was most accurate, because it uses the possibility to discard angular velocity estimates corrupted by clutter filtering. Similar effects were observed for in vivo data obtained at the carotid artery bifurcation. Investigating a combination of these two compounding techniques is recommended for future research.

Least-Squares Multi-Angle Doppler Estimators for Plane-Wave Vector Flow Imaging

Designing robust Doppler vector estimation strategies for use in plane-wave imaging schemes based on unfocused transmissions is a topic that has yet to be studied in depth. One potential solution is to use a multi-angle Doppler estimation approach that computes flow vectors via least-squares fitting, but its performance has not been established. Here, we investigated the efficacy of multi-angle Doppler vector estimators by: 1) comparing its performance with respect to the classical dual-angle (cross-beam) Doppler vector estimator and 2) examining the working effects of multi-angle Doppler vector estimators on flow visualization quality in the context of dynamic flow path rendering. Implementing Doppler vector estimators that use different combinations of transmit (Tx) and receive (Rx) steering angles, our analysis has compared the classical dual-angle Doppler method, a 5-Tx version of dual-angle Doppler, and various multi-angle Doppler configurations based on 3 Tx and 5 Tx. Two angle spans (10°, 20°) were examined in forming the steering angles. In imaging scenarios with known flow profiles (rotating disk and straight-tube parabolic flow), the 3-Tx, 3-Rx and 5-Tx, 5-Rx multi-angle configurations produced vector estimates with smaller variability compared with the dual-angle method, and the estimation results were more consistent with the use of a 20° angle span. Flow vectors derived from multi-angle Doppler estimators were also found to be effective in rendering the expected flow paths in both rotating disk and straight-tube imaging scenarios, while the ones derived from the dual-angle estimator yielded flow paths that deviated from the expected course. These results serve to attest that using multi-angle least-squares Doppler vector estimators, flow visualization can be consistently achieved.

Investigation of Crossbeam Multi-receiver Configurations for Accurate 3-D Vector Doppler Velocity Estimation

An accurate estimation of low blood velocities whose Doppler shifts span the wall filter cutoff, such as near the wall in recirculation or disturbed flow regions, is important for accurate mapping of velocities to achieve improved estimations of wall shear stress and turbulence, which are known risk factors for atherosclerosis and stroke. This paper presents the comparative benefit of increasing the number of receiver beams above three for an improved estimation of low 3-D velocities. The 3-D crossbeam vector Doppler ultrasound configurations were studied in terms of the number of receiver beams, interbeam angle, and beam selection method (criterion for discriminating between tissue and blood Doppler signals) for a range of velocity orientations, which may prove useful in the design of a future 2-D array for vascular imaging. For maximum velocity resolution, a shallow gradient of low flow velocities up to 5 cm/s was generated across a large-diameter (2.46 cm) straight vessel. Data were acquired using a linear array rotated around the central transmit beam axis to generate three- to eight-receiver (3R-8R) configurations;the rotation of each configuration relative to the flow axis was used to mimic a broad range of velocity vector orientations. Accuracy and precision for ≥5 receivers were consistently better over all velocity orientations and for all selection methods. For a velocity magnitude of 2 cm/s, the best accuracy and precision in both magnitude and direction (~21% ± 13%, <1° ± 9°, respectively) were seen with a 5R configuration using a weighted least-squares selection method. Asymmetry in the 5R configuration led to an improved accuracy and precision compared with that in symmetrical 6R and 8R configurations. The results demonstrated relatively little to no benefit from more than five receiver beams.

Cardiovascular Magnetic Resonance Myocardial Feature Tracking: Concepts and Clinical Applications

Heart failure-induced cardiovascular morbidity and mortality constitute a major health problem worldwide and result from diverse pathogeneses, including coronary artery disease, nonischemic cardiomyopathies, and arrhythmias. Assessment of cardiovascular performance is important for early diagnosis and accurate management of patients at risk of heart failure. During the past decade, cardiovascular magnetic resonance myocardial feature tracking has emerged as a useful tool for the quantitative evaluation of cardiovascular function. The method allows quantification of biatrial and biventricular mechanics from measures of deformation: strain, torsion, and dyssynchrony. The purpose of this article is to review the basic principles, clinical applications, accuracy, and reproducibility of cardiovascular magnetic resonance myocardial feature tracking, highlighting the prognostic implications. It will also provide an outlook on how this field might evolve in the future.

Performance study of a new time-delay estimation algorithm in ultrasonic echo signals and ultrasound elastography

Time-delay estimation has countless applications in ultrasound medical imaging. Previously, we proposed a new time-delay estimation algorithm, which was based on the summation of the sign function to compute the time-delay estimate (Shaswary et al., 2015). We reported that the proposed algorithm performs similar to normalized cross-correlation (NCC) and sum squared differences (SSD) algorithms, even though it was significantly more computationally efficient. In this paper, we study the performance of the proposed algorithm using statistical analysis and image quality analysis in ultrasound elastography imaging. Field II simulation software was used for generation of ultrasound radio frequency (RF) echo signals for statistical analysis, and a clinical ultrasound scanner (Sonix® RP scanner, Ultrasonix Medical Corp., Richmond, BC, Canada) was used to scan a commercial ultrasound elastography tissue-mimicking phantom for image quality analysis. The statistical analysis results confirmed that, in overall, the proposed algorithm has similar performance compared to NCC and SSD algorithms. The image quality analysis results indicated that the proposed algorithm produces strain images with marginally higher signal-to-noise and contrast-to-noise ratios compared to NCC and SSD algorithms.

Accurate Angle Estimator for High-Frame-Rate 2-D Vector Flow Imaging

This paper presents a novel approach for estimating 2-D flow angles using a high-frame-rate ultrasound method. The angle estimator features high accuracy and low standard deviation (SD) over the full 360° range. The method is validated on Field II simulations and phantom measurements using the experimental ultrasound scanner SARUS and a flow rig before being tested in vivo. An 8-MHz linear array transducer is used with defocused beam emissions. In the simulations of a spinning disk phantom, a 360° uniform behavior on the angle estimation is observed with a median angle bias of 1.01° and a median angle SD of 1.8°. Similar results are obtained on a straight vessel for both simulations and measurements, where the obtained angle biases are below 1.5° with SDs around 1°. Estimated velocity magnitudes are also kept under 10% bias and 5% relative SD in both simulations and measurements. An in vivo measurement is performed on a carotid bifurcation of a healthy individual. A 3-s acquisition during three heart cycles is captured. A consistent and repetitive vortex is observed in the carotid bulb during systoles.

Influence of ultrasound speckle tracking strategies for motion and strain estimation

Speckle Tracking is one of the most prominent techniques used to estimate the regional movement of the heart based on ultrasound acquisitions. Many different approaches have been proposed, proving their suitability to obtain quantitative and qualitative information regarding myocardial deformation, motion and function assessment. New proposals to improve the basic algorithm usually focus on one of these three steps: (1) the similarity measure between images and the speckle model; (2) the transformation model, i.e. the type of motion considered between images; (3) the optimization strategies, such as the use of different optimization techniques in the transformation step or the inclusion of structural information. While many contributions have shown their good performance independently, it is not always clear how they perform when integrated in a whole pipeline. Every step will have a degree of influence over the following and hence over the final result. Thus, a Speckle Tracking pipeline must be analyzed as a whole when developing novel methods, since improvements in a particular step might be undermined by the choices taken in further steps. This work presents two main contributions: (1) We provide a complete analysis of the influence of the different steps in a Speckle Tracking pipeline over the motion and strain estimation accuracy. (2) The study proposes a methodology for the analysis of Speckle Tracking systems specifically designed to provide an easy and systematic way to include other strategies. We close the analysis with some conclusions and recommendations that can be used as an orientation of the degree of influence of the models for speckle, the transformation models, interpolation schemes and optimization strategies over the estimation of motion features. They can be further use to evaluate and design new strategy into a Speckle Tracking system.

Analysis of Systolic Backflow and Secondary Helical Blood Flow in the Ascending Aorta Using Vector Flow Imaging

Secondary rotational flow and systolic backflow are seen in the ascending aorta and, in this study, were analyzed with the vector velocity method transverse oscillation. Twenty-five patients were scanned intra-operatively, and the vector velocities were related to estimates of transesophageal echocardiography and pulmonary artery catheter thermodilution, and associated with gender, age, aortic diameter, atherosclerotic plaques, left ventricular ejection fraction and previous myocardial infarctions. Secondary flow was present for all patients. The duration and rotational frequency (p < 0.001) and the duration and flow direction of the secondary flow (p < 0.002) were associated. Systolic backflow was present in 40% of the patients and associated with systolic velocities (p < 0.002) and the presence of atherosclerotic plaques (p < 0.001). No other significant associations were observed. The study indicates that backflow is injurious and that secondary flow is a normal flow phenomenon. The study also shows that transverse oscillation can provide new information on blood flow in the ascending aorta.

Plane-wave transverse oscillation for high-frame-rate 2-D vector flow imaging

Transverse oscillation (TO) methods introduce oscillations in the pulse-echo field (PEF) along the direction transverse to the ultrasound propagation direction. This may be exploited to extend flow investigations toward multidimensional estimates. In this paper, the TOs are coupled with the transmission of plane waves (PWs) to reconstruct high-framerate RF images with bidirectional oscillations in the pulse-echo field. Such RF images are then processed by a 2-D phase-based displacement estimator to produce 2-D vector flow maps at thousands of frames per second. First, the capability of generating TOs after PW transmissions was thoroughly investigated by varying the lateral wavelength, the burst length, and the transmission frequency. Over the entire region of interest, the generated lateral wavelengths, compared with the designed ones, presented bias and standard deviation of -3.3 ± 5.7% and 10.6 ± 7.4% in simulations and experiments, respectively. The performance of the ultrafast vector flow mapping method was also assessed by evaluating the differences between the estimated velocities and the expected ones. Both simulations and experiments show overall biases lower than 20% when varying the beam-to-flow angle, the peak velocity, and the depth of interest. In vivo applications of the method on the common carotid and the brachial arteries are also presented.