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

[Echocardiography with high frame rates in the clinical practice : Principles, applications and perspectives]

Continuous developments in cardiovascular imaging, software and hardware have led to technological advancements that open new ways for assessing myocardial mechanics, hemodynamics, and function. Through new scan modalities, echocardiographic scanners can nowadays achieve very high frame rates up to 5000 frames s-1 which enables a wide variety of new applications, including shear wave elastography, ultrafast speckle tracking, the visualization of intracardiac blood flow and myocardial perfusion imaging. This review provides an overview of these advances and demonstrates possible applications and their potential added value in the clinical practice.

Intra-Cardiac Flow from Geometry Prescribed Computational Fluid Dynamics: Comparison with Ultrasound Vector Flow Imaging

PURPOSE

This paper investigates the accuracy of blood flow velocities simulated from a geometry prescribed computational fluid dynamics (CFD) pipeline by applying it to a dynamic heart phantom. The CFD flow patterns are compared to a direct flow measurement by ultrasound vector flow imaging (VFI). The hypothesis is that the simulated velocity magnitudes are within one standard deviation of the measured velocities.

METHODS

The CFD pipeline uses computed tomography angiography (CTA) images with 20 volumes per cardiac cycle as geometry input. Fluid domain movement is prescribed from volumetric image registration using the CTA image data. Inlet and outlet conditions are defined by the experimental setup. VFI is systematically measured in parallel planes, and compared to the corresponding planes in the simulated time dependent three dimensional fluid velocity field.

RESULTS

The measured VFI and simulated CFD have similar flow patterns when compared qualitatively. A quantitative comparison of the velocity magnitude is also performed at specific regions of interest. These are evaluated at 11 non-overlapping time bins and compared by linear regression giving R2 = 0.809, SD = 0.060 m/s, intercept = - 0.039 m/s, and slope = 1.09. Excluding an outlier at the inlet, the correspondence between CFD and VFI improves to: R2 = 0.823, SD = 0.048 m/s, intercept = -0.030 m/s, and slope = 1.01.

CONCLUSION

The direct comparison of flow patterns shows that the proposed CFD pipeline provide realistic flow patterns in a well-controlled experimental setup. The demanded accuracy is obtained close to the inlet and outlet, but not in locations far from these.

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

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

Retrospective Ultrasound Doppler Quantification Using a Single Acquisition in Healthy Adults

OBJECTIVE

Using an experimental tool for retrospective ultrasound Doppler quantification-with high temporal resolution and large spatial coverage-simultaneous flow and tissue measurements were obtained. We compared and validated these experimental values against conventional measurements to determine if the experimental acquisition produced trustworthy tissue and flow velocities.

METHODS

We included 21 healthy volunteers. The only exclusion criterion was the presence of an irregular heartbeat. Two ultrasound examinations were performed for each participant, one using conventional and one using experimental acquisition. The experimental acquisition used multiple plane wave emissions combined with electrocardiography stitching to obtain continuous data with over 3500 frames per second. With two recordings covering a biplane apical view of the left ventricle, we retrospectively extracted selected flow and tissue velocities.

RESULTS

Flow and tissue velocities were compared between the two acquisitions. Statistical testing showed a small but significant difference. We also exemplified the possibility of extracting spectral tissue Doppler from different sample volumes in the myocardium within the imaging sector, showing a decrease in the velocities from the base to the apex.

CONCLUSION

This study demonstrates the feasibility of simultaneous, retrospective spectral and color Doppler of both tissue and flow from an experimental acquisition covering a full sector width. The measurements were significantly different between the two acquisitions but were still comparable, as the biases were small compared to clinical practice, and the two acquisitions were not done simultaneously. The experimental acquisition also enabled the study of deformation by simultaneous spectral velocity traces from all regions of the image sector.

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.

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.

Intracardiac energy inefficiency during atrial fibrillation rhythm with heart failure: Assessment by echo-vector flow mapping

AIMS

Intracardiac dynamics during atrial fibrillation(AF) complicated by heart failure(HF) are not fully understood. The aim of this study was to evaluate the impact of intracardiac dynamics assessed by echo-vector flow mapping on AF complicated by HF.

METHODS AND RESULTS

We enrolled 76 AF patients receiving sinus restoration therapy and assessed energy loss(EL) by echo-vector flow mapping during both AF rhythm and sinus rhythm. Patients were divided into 2 groups according to serum NT-proBNP level: high NT-proBNP group(≥1800 pg/ml during AF rhythm: N = 19), and low NT-proBNP group(N = 57). Outcome measures were average ELs/stroke volume(SV) in left ventricle(LV) and left atrium(LA). Average EL/SVs during AF rhythm in the LV and LA were significantly larger in the high NT-proBNP group than the low NT-proBNP group(54.2 mE/m·L versus 41.2 mE/m·L, P = 0.02; 3.2 mE/m·L versus 1.9 mE/m·L, P = 0.01). The significantly larger EL/SV in the high NT-proBNP group was recorded for maximum EL/SVs. In patients with high NT-proBNP, large vortex formation with extreme EL was detected in the LV and LA during the diastolic phase. After sinus restoration, the average decrease of EL/SV in the LV and LA was larger in the high NT-proBNP group(-21.4 mE/m·L versus +2.6 mE/m·L, P = 0.04; -1.6 mE/m·L versus -0.3 mE/m·L, P = 0.02). Average EL/SV during sinus rhythm did not significantly differ between the high and low NT-proBNP groups in the LV and LA.

CONCLUSIONS

High EL during AF rhythm as intracardiac energy inefficiency was associated with high serum NT-proBNP levels and improved after sinus restoration.

Assessment of Early Diastolic Intraventricular Pressure Difference in Children by Blood Speckle-Tracking Echocardiography

BACKGROUND

The lack of reliable echocardiographic techniques to assess diastolic function in children is a major clinical limitation. Our aim was to develop and validate the intraventricular pressure difference (IVPD) calculation using blood speckle-tracking (BST) and investigate the method's potential role in the assessment of diastolic function in children.

METHODS

Blood speckle-tracking allows two-dimensional angle-independent blood flow velocity estimation. Blood speckle-tracking images of left ventricular (LV) inflow from the apical 4-chamber view in 138 controls, 10 patients with dilated cardiomyopathies (DCMs), and 21 patients with hypertrophic cardiomyopathies (HCMs) <18 years of age were analyzed to study LV IVPD during early diastole. Reproducibility of the IVPD analysis was assessed, IVPD estimates from BST and color M mode were compared, and the validity of the BST-based IVPD calculations was tested in a computer flow model.

RESULTS

Mean IVPD was significantly higher in controls (-2.28 ± 0.62 mm Hg) compared with in DCM (-1.21 ± 0.39 mm Hg, P < .001) and HCM (-1.57 ± 0.47 mm Hg, P < .001) patients. Feasibility was 88.3% in controls, 80% in DCM patients, and 90.4% in HCM patients. The peak relative negative pressure occurred earlier at the apex than at the base and preceded the peak E-wave LV filling velocity, indicating that it represents diastolic suction. Intraclass correlation coefficients for intra- and interobserver variability were 0.908 and 0.702, respectively. There was a nonsignificant mean difference of 0.15 mm Hg between IVPD from BST and color M mode. Estimation from two-dimensional velocities revealed a difference in peak IVPD of 0.12 mm Hg (6.6%) when simulated in a three-dimensional fluid mechanics model.

CONCLUSIONS

Intraventricular pressure difference calculation from BST is highly feasible and provides information on diastolic suction and early filling in children with heart disease. Intraventricular pressure difference was significantly reduced in children with DCM and HCM compared with controls, indicating reduced early diastolic suction in these patient groups.

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

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

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.

Left ventricular high frame rate echo-particle image velocimetry: clinical application and comparison with conventional imaging

Background

Echo-Particle Image Velocimetry (echoPIV) tracks speckle patterns from ultrasound contrast agent(UCA), being less angle-sensitive than colour Doppler. High frame rate (HFR) echoPIV enables tracking of high velocity flow in the left ventricle (LV). We aimed to demonstrate the potential clinical use of HFR echoPIV and investigate the feasibility and accuracy in patients.

Methods

Nineteen patients admitted for heart failure were included. HFR contrast images were acquired from an apical long axis view (ALAX), using a fully-programmable ultrasound system. A clinical UCA was continuously infused with a dedicated pump. Additionally, echocardiographic images were obtained using a clinical system, including LV contrast-enhanced images and pulsed-wave (PW) Doppler of the LV inflow and outflow in ALAX. 11 patients underwent CMR and 4 cardiac CT as clinically indicated. These CMR and CT images were used as reference. In 10 patients with good echoPIV tracking and reference imaging, the intracavitary flow was compared between echoPIV, conventional and UCA echocardiography.

Results

EchoPIV tracking quality was good in 12/19 (63%), moderate in 2/19 (10%) and poor in 5/19 (26%) subjects. EchoPIV could determine inflow velocity in 17/19 (89%), and outflow in 14/19 (74%) patients. The correlation of echoPIV and PW Doppler was good for the inflow (R² = 0.77 to PW peak; R² = 0.80 PW mean velocity) and moderate for the outflow (R² = 0.54 to PW peak; R² = 0.44 to PW mean velocity), with a tendency for echoPIV to underestimate PW velocities. In selected patients, echoPIV was able in a single acquisition to demonstrate flow patterns which required multiple interrogations with classical echocardiography. Those flow patterns could also be linked to anatomical abnormalities as seen in CMR or CT.

Conclusion

HFR echoPIV tracks multidirectional and complex flow patterns which are unapparent with conventional echocardiography, while having comparable feasibility. EchoPIV tends to underestimate flow velocities as compared to PW Doppler. It has the potential to provide in one acquisition all the functional information obtained by conventional imaging, overcoming the angle dependency of Doppler and low frame rate of classical contrast imaging.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12947-022-00283-4.

Performance Assessment of Row-Column Transverse Oscillation Tensor Velocity Imaging Using Computational Fluid Dynamics Simulation of Carotid Bifurcation Flow

In this work, the accuracy of row-column tensor velocity imaging (TVI), i.e., 3-D vector flow imaging (VFI) in 3-D space over time, is quantified on a complex, clinically relevant flow. The quantification is achieved by transferring the flow simulated using computational fluid dynamics (CFD) to a Field II simulation environment, and this allows for a direct comparison between the actual and estimated velocities. The carotid bifurcation flow simulations were performed with a peak inlet velocity of 80 cm/s, nonrigid vessel walls, and a flow cycle duration of 1.2 s. The flow was simulated from two observation angles, and it was acquired using a 3-MHz 62+62 row-column addressed array (RCA) at a pulse repetition frequency ( f_prf ) of 10 and 20 kHz. The tensor velocities were obtained at a frame rate of 208.3 Hz, at f_prf = 10 kHz , and the results from two velocity estimators were compared. The two estimators were the directional transverse oscillation (TO) cross correlation estimator and the proposed autocorrelation estimator. Linear regression between the actual and estimated velocity components yielded, for the cross correlation estimator, an R ² value in the range of 0.89-0.91, 0.46-0.77, and 0.91-0.97 for the x -, y -, and z -components, and 0.87-0.89, 0.40-0.83, and 0.91-0.96 when using the autocorrelation estimator. The results demonstrate that an RCA can, with just 62 receive channels, measure complex 3-D flow fields at a high volume rate.

Contrast Agent-Free Assessment of Blood Flow and Wall Shear Stress in the Rabbit Aorta using Ultrasound Image Velocimetry

Blood flow velocity and wall shear stress (WSS) influence and are influenced by vascular disease. Their measurement is consequently useful in the laboratory and clinic. Contrast-enhanced ultrasound image velocimetry (UIV) can estimate them accurately but the need to inject contrast agents limits utility. Singular value decomposition and high-frame-rate imaging may render contrast agents dispensable. Here we determined whether contrast agent-free UIV can measure flow and WSS. In simulation, accurate measurements were achieved with a signal-to-noise ratio of 13.5 dB or higher. Signal intensity in the rabbit aorta increased monotonically with mechanical index; it was lowest during stagnant flow and uneven across the vessel. In vivo measurements with contrast-free and contrast-enhanced UIV differed by 4.4% and 1.9% for velocity magnitude and angle and by 9.47% for WSS. Bland-Altman analysis of waveforms revealed good agreement between contrast-free and contrast-enhanced UIV. In five rabbits, the root-mean-square errors were as low as 0.022 m/s (0.81%) and 0.11 Pa (1.7%). This study indicates that with an optimised protocol, UIV can assess flow and WSS without contrast agents. Unlike contrast-enhanced UIV, contrast-free UIV could be routinely employed.

Intraventricular Vector Flow Imaging with Blood Speckle Tracking in Adults: Feasibility, Normal Physiology and Mechanisms in Healthy Volunteers

This study examines the feasibility of blood speckle tracking for vector flow imaging in healthy adults and describes the physiologic flow pattern and vortex formation in relation to the wall motion in the left ventricle. The study included 21 healthy volunteers and quantified and visualized flow patterns with high temporal resolution down to a depth of 10-12 cm without the use of contrast agents. Intraventricular flow seems to originate during the isovolumetric relaxation with a propagation of blood from base to apex. With the E-wave, rapid inflow and vortex formation occurred on both sides of the valve basally. During diastasis the flow gathers in a large vortex before the pattern from the E-wave repeats during the A-wave. In isovolumetric contraction, the flow again gathers in a large vortex that seems to facilitate the flow out in the aorta during systole. No signs of a persistent systolic vortex were visualized. The geometry of the left ventricle and the movement of the AV-plane is important in creating vortices that are favorable for the blood flow and facilitate outflow. The quantitative measurements are in concordance with these findings, but the clinical interpretation must be evaluated in future clinical studies.

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.

Optimization of Microbubble Concentration and Acoustic Pressure for Left Ventricular High-Frame-Rate EchoPIV in Patients

High-frame-rate (HFR) echo-particle image velocimetry (echoPIV) is a promising tool for measuring intracardiac blood flow dynamics. In this study, we investigate the optimal ultrasound contrast agent (UCA: SonoVue) infusion rate and acoustic output to use for HFR echoPIV (PRF = 4900 Hz) in the left ventricle (LV) of patients. Three infusion rates (0.3, 0.6, and 1.2 ml/min) and five acoustic output amplitudes (by varying transmit voltage: 5, 10, 15, 20, and 30 V-corresponding to mechanical indices of 0.01, 0.02, 0.03, 0.04, and 0.06 at 60-mm depth) were tested in 20 patients admitted for symptoms of heart failure. We assess the accuracy of HFR echoPIV against pulsed-wave Doppler acquisitions obtained for mitral inflow and aortic outflow. In terms of image quality, the 1.2-ml/min infusion rate provided the highest contrast-to-background ratio (CBR) (3-dB improvement over 0.3 ml/min). The highest acoustic output tested resulted in the lowest CBR. Increased acoustic output also resulted in increased microbubble disruption. For the echoPIV results, the 1.2-ml/min infusion rate provided the best vector quality and accuracy; mid-range acoustic outputs (corresponding to 15-20-V transmit voltages) provided the best agreement with the pulsed-wave Doppler. Overall, the highest infusion rate (1.2 ml/min) and mid-range acoustic output amplitudes provided the best image quality and echoPIV results.

Right Ventricular Flow Dynamics in Dilated Right Ventricles: Energy Loss Estimation Based on Blood Speckle Tracking Echocardiography-A Pilot Study in Children

Using blood speckle tracking (BST) based on high-frame-rate echocardiography (HFRE), we compared right ventricle (RV) flow dynamics in children with atrial septal defects (ASDs) and repaired tetralogy of Fallot (rTOF). Fifty-seven children with rTOF with severe pulmonary insufficiency (PI) (n = 21), large ASDs (n = 11) and healthy controls (CTL, n = 25) were included. Using a flow phantom, we studied the effects of imaging plane and smoothing parameters on 2-D energy loss (EL). RV diastolic EL was similar in ASD and rTOF, but both were greater than in CTL. Locations of high EL were similar in all groups in systole, occurring in the RV outflow tract and around the tricuspid valve leaflets in early diastole. An additional apical early diastolic area of EL was noted in rTOF, corresponding to colliding tricuspid inflow and PI. The flow phantom revealed that EL varied with imaging plane and smoothing settings but that the EL trend was preserved if kept consistent.

Real-Time 3-D Spectral Doppler Analysis With a Sparse Spiral Array

2-D sparse arrays may push the development of low-cost 3-D systems, not needing to control thousands of elements by expensive application-specific integrated circuits (ASICs). However, there is still some concern about their suitability in applications, such as Doppler investigation, which inherently involve poor signal-to-noise ratios (SNRs). In this article, a novel real-time 3-D pulsed-wave (PW) Doppler system, based on a 256-element 2-D spiral array, is presented. Coded transmission (TX) and matched filtering were implemented to improve the system SNR. Standard sonograms as well as multigate spectral Doppler (MSD) profiles, along lines that can be arbitrarily located in different planes, are presented. The performance of the system was assessed quantitatively on experimental data obtained from a straight tube flow phantom. An SNR increase of 11.4 dB was measured by transmitting linear chirps instead of standard sinusoidal bursts. For a qualitative assessment of the system performance in more realistic conditions, an anthropomorphic phantom of the carotid arteries was used. Finally, real-time B-mode and MSD images were obtained from healthy volunteers.

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.

Tapered Vector Doppler for Improved Quantification of Low Velocity Blood Flow

A new vector velocity estimation scheme is developed, termed tapered vector Doppler (TVD), aiming to improve the accuracy of low velocity flow estimation. This is done by assessing the effects of singular value decomposition (SVD) and finite impulse response (FIR) filters and designing an estimator which accounts for signal loss due to filtering. Synthetic data created using a combination of in vivo recordings and flow simulations were used to investigate scenarios with low blood flow, in combination with true clutter motion. Using this approach, the accuracy and precision of TVD was investigated for a range of clutter-to-blood and signal-to-noise ratios. The results indicated that for the investigated carotid application and setup, the SVD filter performed as a frequency-based filter. For both SVD and FIR filters, suppression of the clutter signal resulted in large bias and variance in the estimated blood velocity magnitude and direction close to the vessel walls. Application of the proposed tapering technique yielded significant improvement in the accuracy and precision of near-wall vector velocity measurements, compared to non-TVD and weighted least squares approaches. In synthetic data, for a blood SNR of 5 dB, and in a near-wall region where the average blood velocity was 9 cm/s, the use of tapering reduced the average velocity magnitude bias from 26.3 to 1.4 cm/s. Complex flow in a carotid bifurcation was used to demonstrate the in vivo performance of TVD, and it was shown that tapering enables vector velocity estimation less affected by clutter and clutter filtering than what could be obtained by adaptive filter design only.

Characterization of Vortex Flow in a Mouse Model of Ventricular Dyssynchrony by Plane-Wave Ultrasound Using Hexplex Processing

The rodent heart is frequently used to study human cardiovascular disease (CVD). Although advanced cardiovascular ultrasound imaging methods are available for human clinical practice, application of these techniques to small animals remains limited due to the temporal and spatial-resolution demands. Here, an ultrasound vector-flow workflow is demonstrated that enables visualization and quantification of the complex hemodynamics within the mouse heart. Wild type (WT) and fibroblast growth factor homologous factor 2 (FHF2)-deficient mice (Fhf2 ^KO/Y ), which present with hyperthermia-induced ECG abnormalities highly reminiscent of Brugada syndrome, were used as a mouse model of human CVD. An 18-MHz linear array was used to acquire high-speed (30 kHz), plane-wave data of the left ventricle (LV) while increasing core body temperature up to 41.5 °C. Hexplex (i.e., six output) processing of the raw data sets produced the output of vector-flow estimates (magnitude and phase); B-mode and color-Doppler images; Doppler spectrograms; and local time histories of vorticity and pericardium motion. Fhf2 ^WT/Y mice had repeatable beat-to-beat cardiac function, including vortex formation during diastole, at all temperatures. In contrast, Fhf2 ^KO/Y mice displayed dyssynchronous contractile motion that disrupted normal inflow vortex formation and impaired LV filling as temperature rose. The hexplex processing approach demonstrates the ability to visualize and quantify the interplay between hemodynamic and mechanical function in a mouse model of human CVD.

Fast 3-D Velocity Estimation in 4-D Using a 62 + 62 Row-Column Addressed Array

This article presents an imaging scheme capable of estimating the full 3-D velocity vector field in a volume using row-column addressed arrays (RCAs) at a high volume rate. A 62 + 62 RCA array is employed with an interleaved synthetic aperture sequence. It contains repeated emissions with rows and columns interleaved with B-mode emissions. The sequence contains 80 emissions in total and can provide continuous volumetric data at a volume rate above 125 Hz. A transverse oscillation cross correlation estimator determines all three velocity components. The approach is investigated using Field II simulations and measurements using a specially built 3-MHz 62 + 62 RCA array connected to the SARUS experimental scanner. Both the B-mode and flow sequences have a penetration depth of 14 cm when measured on a tissue-mimicking phantom (0.5-dB/[ [Formula: see text]] attenuation). Simulations of a parabolic flow in a 12-mm-diameter vessel at a depth of 30 mm, beam-to-flow angle of 90°, and xy-rotation of 45° gave a standard deviation (SD) of (3.3, 3.4, 0.4)% and bias of (-3.3, -3.9, -0.1)%, for ( v_x , v_y , and v_z ). Decreasing the beam-to-flow angle to 60° gave an SD of (8.9, 9.1, 0.8)% and bias of (-7.6, -9.5, -7.2)%, showing a slight increase. Measurements were carried out using a similar setup, and pulsing at 2 kHz yielded comparable results at 90° with an SD of (5.8, 5.5, 1.1)% and bias of (1.4, -6.4, 2.4)%. At 60°, the SD was (5.2, 4.7 1.2)% and bias (-4.6, 6.9, -7.4)%. Results from measurements across all tested settings showed a maximum SD of 6.8% and a maximum bias of 15.8% for a peak velocity of 10 cm/s. A tissue-mimicking phantom with a straight vessel was used to introduce clutter, tissue motion, and pulsating flow. The pulsating velocity magnitude was estimated across ten pulse periods and yielded an SD of 10.9%. The method was capable of estimating transverse flow components precisely but underestimated the flow with small beam-to-flow angles. The sequence provided continuous data in both time and space throughout the volume, allowing for retrospective analysis of the flow. Moreover, B-mode planes can be selected retrospectively anywhere in the volume. This shows that tensor velocity imaging (full 3-D volumetric vector flow imaging) can be estimated in 4-D ( x, y, z, and t ) using only 62 channels in receive, making 4-D volumetric imaging implementable on current scanner hardware.

Helical toroid phantom for 3D flow imaging investigations

The medical physics community has hitherto lacked an effective calibration phantom to holistically evaluate the performance of three-dimensional (3D) flow imaging techniques. Here, we present the design of a new omnidirectional, three-component (3-C) flow phantom whose lumen is consisted of a helical toroid structure (4 mm lumen diameter; helically winded for 5 revolutions over a torus with 10 mm radius; 5 mm helix radius). This phantom's intraluminal flow trajectory embraces all combinations of x, y, and z directional components, as confirmed using computational fluid dynamics (CFD) simulations. The phantom was physically fabricated via lost-core casting with polyvinyl alcohol cryogel (PVA) as the tissue mimic. 3D ultrasound confirmed that the phantom lumen expectedly resembled a helical toroid geometry. Pulsed Doppler measurements showed that the phantom, when operating under steady flow conditions (3 ml s^-1 flow rate), yielded flow velocity magnitudes that agreed well with those derived from CFD at both the inner torus (-47.6 ± 5.7 versus -52.0 ± 2.2 cm s^-1; mean ± 1 S.D.) and the outer torus (49.5 ± 4.2 versus 48.0 ± 1.7 cm s^-1). Additionally, 3-C velocity vectors acquired from multi-angle pulsed Doppler showed good agreement with CFD-derived velocity vectors (<7% and 10° difference in magnitude and flow angle, respectively). Ultrasound color flow imaging further revealed that the phantom's axial flow pattern was aligned with the CFD-derived flow profile. Overall, the helical toroid phantom has strong potential as an investigative tool in 3D flow imaging innovation endeavors, such as the development of flow vector estimators and visualization algorithms.

Concepts and applications of ultrafast cardiac ultrasound imaging

The concept of ultrafast echocardiographic imaging has been around for decades. However, only recent progress in ultrasound machine hardware and computer technology allowed to apply this concept to echocardiography. High frame rate echocardiography can visualize phenomena that have never been captured before. It enables a wide variety of potential new applications, including shear wave imaging, speckle tracking, ultrafast Doppler imaging, and myocardial perfusion imaging. The principles of these applications and their potential clinical use will be presented in this manuscript.

Non-invasive estimation of relative pressure for intracardiac flows using virtual work-energy

Intracardiac blood flow is driven by differences in relative pressure, and assessing these is critical in understanding cardiac disease. Non-invasive image-based methods exist to assess relative pressure, however, the complex flow and dynamically moving fluid domain of the intracardiac space limits assessment. Recently, we proposed a method, νWERP, utilizing an auxiliary virtual field to probe relative pressure through complex, and previously inaccessible flow domains. Here we present an extension of νWERP for intracardiac flow assessments, solving the virtual field over sub-domains to effectively handle the dynamically shifting flow domain. The extended νWERP is validated in an in-silico benchmark problem, as well as in a patient-specific simulation model of the left heart, proving accurate over ranges of realistic image resolutions and noise levels, as well as superior to alternative approaches. Lastly, the extended νWERP is applied on clinically acquired 4D Flow MRI data, exhibiting realistic ventricular relative pressure patterns, as well as indicating signs of diastolic dysfunction in an exemplifying patient case. Summarized, the extended νWERP approach represents a directly applicable implementation for intracardiac flow assessments.

Myelin detection in fluorescence microscopy images using machine learning

BACKGROUND

The myelin sheath produced by glial cells insulates the axons, and supports the function of the nervous system. Myelin sheath degeneration causes neurodegenerative disorders, such as multiple sclerosis (MS). There are no therapies for MS that promote remyelination. Drug discovery frequently involves screening thousands of compounds. However, this is not feasible for remyelination drugs, since myelin quantification is a manual labor-intensive endeavor. Therefore, the development of assistive software for expedited myelin detection is instrumental for MS drug discovery by enabling high-content image-based drug screens.

NEW METHOD

In this study, we developed a machine learning based expedited myelin detection approach in fluorescence microscopy images. Multi-channel three-dimensional microscopy images of a mouse stem cell-based myelination assay were labeled by experts. A spectro-spatial feature extraction method was introduced to represent local dependencies of voxels both in spatial and spectral domains. Feature extraction yielded two data set of over forty-seven thousand annotated images in total.

RESULTS

Myelin detection performances of 23 different supervised machine learning techniques including a customized-convolutional neural network (CNN), were assessed using various train/test split ratios of the data sets. The highest accuracy values of 98.84±0.09% and 98.46±0.11% were achieved by Boosted Trees and customized-CNN, respectively.

COMPARISON WITH EXISTING METHODS

Our approach can detect myelin in a common experimental setup. Myelin extending in any orientation in 3 dimensions is segmented from 3 channel z-stack fluorescence images.

CONCLUSIONS

Our results suggest that the proposed expedited myelin detection approach is a feasible and robust method for remyelination drug screening.

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

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

3-D Intravascular Characterization of Blood Flow Velocity Fields with a Forward-Viewing 2-D Array

Risk stratification in coronary artery disease is an ongoing challenge for which few tools are available for quantifying physiology within coronary arteries. Recently, anatomy-driven computational fluid dynamic modeling has enabled the mapping of local flow dynamics in coronary stenoses, with derived parameters such as WSS exhibiting a strong capability for predicting adverse clinical events on a patient-specific basis. As cardiac catheterization is common in patients with coronary artery disease, minimally invasive technologies capable of identifying pathologic flow in situ in real time could have a significant impact on clinical decision- making. As a step toward in vivo quantification of slow flow near the arterial wall, proof-of-concept for 3-D intravascular imaging of blood flow dynamics is provided using a 118-element forward-viewing ring array transducer and a research ultrasound system. Blood flow velocity components are estimated in the direction of primary flow using an unfocused wave Doppler approach, and in the lateral and elevation directions, using a transverse oscillation approach. This intravascular 3-D vector velocity system is illustrated by acquiring real-time 3-D data sets in phantom experiments and in vivo in the femoral artery of a pig. The effect of the catheter on blood flow dynamics is also experimentally assessed in flow phantoms with both straight and stenotic vessels. Results indicate that 3-D flow dynamics can be measured using a small form factor device and that a hollow catheter design may provide minimal disturbance to flow measurements in a stenosis (peak velocity: 54.97 ± 2.13 cm/s without catheter vs. 51.37 ± 1.08 cm/s with hollow catheter, 6.5% error). In the future, such technologies could enable estimation of 3-D flow dynamics near the wall in patients already undergoing catheterization.

High-Frame-Rate Color Doppler Echocardiography: A Quantitative Comparison of Different Approaches

Ultrasound color Doppler imaging (CDI) provides a map of the axial blood flow velocities in a 2-D/3-D region of interest. While CDI is clinically effective for a qualitative analysis of abnormal blood flows, e.g., for valvular disease in cardiology, it is in limited use for quantitative measures, mainly hampered by low frame rate and measurement bias. These limitations can be reduced by different approaches toward high-frame-rate (HFR) imaging at the expense of reduced image quality and penetration depth. The aim of this study was to compare the impact of different HFR sequences on CDI quantitatively. Different cardiac scan sequences, including diverging waves and multiline transmission, were designed, implemented on a research system, and compared in terms of patient safety parameters, image quality, and penetration depth. Furthermore, in vivo images were acquired and compared for healthy volunteers. Results showed that the HFR techniques spread artifacts on larger areas than the standard single-line scans (> +50%). In addition, due to patient safety limitations, they reduce the penetration depth up to -5 cm. On the other hand, the HFR techniques provide comparable velocity estimates (relative difference <6%) and enhance the time resolution of the color Doppler images, achieving frame rates up to 625 Hz in continuous acquisition.

4-D Echo-Particle Image Velocimetry in a Left Ventricular Phantom

Left ventricular (LV) blood flow is an inherently complex time-varying 3-D phenomenon, where 2-D quantification often ignores the effect of out-of-plane motion. In this study, we describe high frame rate 4-D echocardiographic particle image velocimetry (echo-PIV) using a prototype matrix transesophageal transducer and a dynamic LV phantom for testing the accuracy of echo-PIV in the presence of complex flow patterns. Optical time-resolved tomographic PIV (tomo-PIV) was used as a reference standard for comparison. Echo-PIV and tomo-PIV agreed on the general profile of the LV flow patterns, but echo-PIV smoothed out the smaller flow structures. Echo-PIV also underestimated the flow rates at greater imaging depths, where the PIV kernel size and transducer point spread function were large relative to the velocity gradients. We demonstrate that 4-D echo-PIV could be performed in just four heart cycles, which would require only a short breath-hold, providing promising results. However, methods for resolving high velocity gradients in regions of poor spatial resolution are required before clinical translation.

Spectral Doppler Measurements With 2-D Sparse Arrays

The 2-D sparse arrays, in which a few hundreds of elements are distributed on the probe surface according to an optimization procedure, represent an alternative to full 2-D arrays, including thousands of elements usually organized in a grid. Sparse arrays have already been used in B-mode imaging tests, but their application to Doppler investigations has not been reported yet. Since the sparsity of the elements influences the acoustic field, a corresponding influence on the mean frequency (Fm), bandwidth (BW), and signal-to-noise ratio (SNR) of the Doppler spectra is expected. This article aims to assess, by simulations and experiments, to what extent the use of a sparse rather than a full gridded 2-D array has an impact on spectral Doppler measurements. Parabolic flows were investigated by a 3 MHz, 1024-element gridded array and by a sparse array; the latter was obtained by properly selecting a subgroup of 256 elements from the full array. Simulations show that the mean Doppler frequency does not change between the sparse and the full array while there are significant differences on the BW (average reduction of 17.2% for the sparse array, due to different apertures of the two probes) and on the signal power (Ps) (22 dB, due to the different number of active elements). These results are confirmed by flow phantom experiments, which also highlight that the most critical difference between sparse and full gridded array in Doppler measurements is in terms of SNR (-16.8 dB).

Blood Speckle-Tracking Based on High-Frame Rate Ultrasound Imaging in Pediatric Cardiology

BACKGROUND

Flow properties play an important role in cardiac function, remodeling, and morphogenesis but cannot be displayed in detail with today's echocardiographic techniques. The authors hypothesized that blood speckle-tracking (BST) could visualize and quantify flow patterns. The aim of this study was to determine the feasibility, accuracy, and potential clinical applications of BST in pediatric cardiology.

METHODS

BST is based on high-frame rate ultrasound, using a combination of plane-wave imaging and parallel receive beamforming. Pattern-matching techniques are used to quantify blood speckle motion. Accuracy of BST velocity measurements was validated using a rotating phantom and by comparing BST-derived inflow velocities with pulsed-wave Doppler obtained in the left ventricles of healthy control subjects. To test clinical feasibility, 102 subjects (21 weeks to 11.5 years of age) were prospectively enrolled, including healthy fetuses (n = 4), healthy control subjects (n = 51), and patients with different cardiac diseases (n = 47).

RESULTS

The phantom data showed a good correlation (r = 0.95, with a tracking quality threshold of 0.4) between estimated BST velocities and reference velocities down to a depth of 8 cm. There was a good correlation (r = 0.76) between left ventricular inflow velocity measured using BST and pulsed-wave Doppler. BST displayed lower velocities (mean ± SD, 0.59 ± 0.14 vs 0.82 ± 0.21 m/sec for pulsed-wave Doppler). However, the velocity amplitude in BST increases with reduced smoothing. The clinical feasibility of BST was high, as flow patterns in the area of interest could be visualized in all but one case (>99%).

CONCLUSIONS

BST is highly feasible in fetal and pediatric echocardiography and provides a novel approach for visualizing blood flow patterns. BST provides accurate velocity measurements down to 8 cm, but compared with pulsed-wave Doppler, BST displays lower velocities. Studying blood flow properties may provide novel insights into the pathophysiology of pediatric heart disease and could become an important diagnostic tool.

High-Frame-Rate Tri-Plane Echocardiography With Spiral Arrays: From Simulation to Real-Time Implementation

Major cardiovascular diseases (CVDs) are associated with (regional) dysfunction of the left ventricle. Despite the 3-D nature of the heart and its dynamics, the assessment of myocardial function is still largely based on 2-D ultrasound imaging, thereby making diagnosis heavily susceptible to the operator's expertise. Unfortunately, to date, 3-D echocardiography cannot provide adequate spatiotemporal resolution in real-time. Hence, tri-plane imaging has been introduced as a compromise between 2-D and true volumetric ultrasound imaging. However, tri-plane imaging typically requires high-end ultrasound systems equipped with fully populated matrix array probes embedded with expensive and little flexible electronics for two-stage beamforming. This article presents an advanced ultrasound system for real-time, high frame rate (HFR), and tri-plane echocardiography based on low element count sparse arrays, i.e., the so-called spiral arrays. The system was simulated, experimentally validated, and implemented for real-time operation on the ULA-OP 256 system. Five different array configurations were tested together with four different scan sequences, including multi-line and planar diverging wave transmission. In particular, the former can be exploited to achieve, in tri-plane imaging, the same temporal resolution currently used in clinical 2-D echocardiography, at the expenses of contrast (-3.5 dB) and signal-to-noise ratio (SNR) (-8.7 dB). On the other hand, the transmission of planar diverging waves boosts the frame rate up to 250 Hz, but further compromises contrast (-10.5 dB), SNR (-9.7 dB), and lateral resolution (+46%). In conclusion, despite an unavoidable loss in image quality and sensitivity due to the limited number of elements, HFR tri-plane imaging with spiral arrays is shown to be feasible in real-time and may enable real-time functional analysis of all left ventricular segments of the heart.

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.