Fast Plane Wave 2-D Vector Flow Imaging Using Transverse Oscillation and Directional Beamforming

A Pseudo-Spectral Method for Wall Shear Stress Estimation from Doppler Ultrasound Imaging in Coronary Arteries

PURPOSE

The Wall Shear Stress (WSS) is the component tangential to the boundary of the normal stress tensor in an incompressible fluid, and it has been recognized as a quantity of primary importance in predicting possible adverse events in cardiovascular diseases, in general, and in coronary diseases, in particular. The quantification of the WSS in patient-specific settings can be achieved by performing a Computational Fluid Dynamics (CFD) analysis based on patient geometry, or it can be retrieved by a numerical approximation based on blood flow velocity data, e.g., ultrasound (US) Doppler measurements. This paper presents a novel method for WSS quantification from 2D vector Doppler measurements.

METHODS

Images were obtained through unfocused plane waves and transverse oscillation to acquire both in-plane velocity components. These velocity components were processed using pseudo-spectral differentiation techniques based on Fourier approximations of the derivatives to compute the WSS.

RESULTS

Our Pseudo-Spectral Method (PSM) is tested in two vessel phantoms, straight and stenotic, where a steady flow of 15 mL/min is applied. The method is successfully validated against CFD simulations and compared against current techniques based on the assumption of a parabolic velocity profile. The PSM accurately detected Wall Shear Stress (WSS) variations in geometries differing from straight cylinders, and is less sensitive to measurement noise. In particular, when using synthetic data (noise free, e.g., generated by CFD) on cylindrical geometries, the Poiseuille-based methods and PSM have comparable accuracy; on the contrary, when using the data retrieved from US measures, the average error of the WSS obtained with the PSM turned out to be 3 to 9 times smaller than that obtained by state-of-the-art methods.

CONCLUSION

The pseudo-spectral approach allows controlling the approximation errors in the presence of noisy data. This gives a more accurate alternative to the present standard and a less computationally expensive choice compared to CFD, which also requires high-quality data to reconstruct the vessel geometry.

Dynamic Coronary Blood Flow Velocity and Wall Shear Stress Estimation Using Ultrasound in an Ex Vivo Porcine Heart

PURPOSE

Wall shear stress (WSS) is a critically important physical factor contributing to atherosclerosis. Mapping the spatial distribution of local, oscillatory WSS can identify important mechanisms underlying the progression of coronary artery disease.

METHODS

In this study, blood flow velocity and time-varying WSS were estimated in the left anterior descending (LAD) coronary artery of an ex vivo beating porcine heart using ultrasound with an 18 MHz linear array transducer aligned with the LAD in a forward-viewing orientation. A pulsatile heart loop with physiologically-accurate flow was created using a pulsatile pump. The coronary artery wall motion was compensated using a local block matching technique. Next, 2D and 3D velocity magnitude and WSS maps in the LAD coronary artery were estimated at different time points in the cardiac cycle using an ultrafast Doppler approach. The blood flow velocity estimated using the presented approach was compared with a commercially-available, calibrated single element blood flow velocity measurement system.

RESULTS

The resulting root mean square error (RMSE) of 2D velocity magnitude acquired from a high frequency, linear array transducer was less than 8% of the maximum velocity estimated by the commercial system.

CONCLUSION

When implemented in a forward-viewing intravascular ultrasound device, the presented approach will enable dynamic estimation of WSS, an indicator of plaque vulnerability in coronary arteries.

Thermal analysis and SNR comparison of CMUT and PZT transducers using coded excitation

Coded excitation (CE) has the ability to enhance image quality and penetration depth by improving the signal-to-noise ratio (SNR). Their usefulness has been extensively proven in the literature, however, there are very few publications that have discussed the practicality of using CEs, as they can increase the operating temperature of a transducer beyond the safety limits. In this paper, the potential for capacitive micromachined ultrasonic transducers (CMUTs) to handle CEs is investigated and compared to a geometrically similar Lead Zirconate Titanate (PZT) probe. It is hypothesized that CMUTs are comparatively advantageous for CE and can generate CE signals more effectively while operating within the safety limits of temperature specified by the International Electrotechnical Commission (IEC). Simple chirp signals were used for excitation and the temperature of the transducer assembly was measured under two different test conditions following the IEC standards. The still air test showed that at the pulse repetition frequency of 4 kHz and a signal duration of 8μs, the PZT probe reached the safe limit of 27 °C within 3 min of acquisition, whereas the CMUT probe showed an increase of only 4-5 °C over 30 min of scan time. The phantom experiments showed that for the maximum temperature rise allowed at the object-transducer interface by the IEC, the gain in SNR possible for the CMUT probe was 13-14 dB and the increase in penetration depth was 30 mm, whereas no gain in SNR and penetration depth was achieved for the PZT probe. The peak mechanical index and the derated spatial peak temporal averaged intensity were found to be less than half the FDA limits for both probes in all cases, proving temperature to be the first limiting factor when using CE. Therefore, the CMUT has shown to be thermally more efficient than the PZT probe and a good candidate for CE imaging.

Frame rate improvement in ultrafast coherent plane wave compounding

Coherent plane wave compounding (CPWC), as an ultrafast ultrasound imaging technique, makes a significant breakthrough in frame rate enhancement. However, there exists a compromise between the quality of the final image and the frame rate in CPWC. In this paper, we propose an efficient method to minimize the number of required emissions, and consequently, improve the frame rate, while maintaining the image quality. To this end, we down-sample the angle interval using two specific sampling factors. More precisely, we construct two different subsets, each of which consists of a few numbers of emissions. The optimal values of the angle intervals are achieved based on the beampattern that corresponds to the reference case (that is, the case where all plane waves are used). Finally, in order to keep the image quality comparable with the reference case, we apply some modifications to the image reconstruction procedure. In the proposed algorithm, the Delay-and-Sum beamformed images of two considered subsets are convolved to achieve the final reconstructed image. The obtained results confirm the efficiency of the proposed method in terms of frame rate improvement compared to the reference case. In particular, by using the proposed method, the required emissions in PICMUS data reduce to 16, which is 4.6 times smaller compared to the reference case. Also, the gCNR values of the proposed method and the reference case are obtained as 0.98 and 0.97, respectively, for in-vivo dataset. This demonstrates that the proposed method successfully preserves the quality of the reconstructed image by using much fewer emissions.

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.

Fast volumetric ultrasound facilitates high-resolution 3D mapping of tissue compartments

Volumetric ultrasound imaging has the potential for operator-independent acquisition and enhanced field of view. Panoramic acquisition has many applications across ultrasound; spanning musculoskeletal, liver, breast, and pediatric imaging; and image-guided therapy. Challenges in high-resolution human imaging, such as subtle motion and the presence of bone or gas, have limited such acquisition. These issues can be addressed with a large transducer aperture and fast acquisition and processing. Programmable, ultrafast ultrasound scanners with a high channel count provide an unprecedented opportunity to optimize volumetric acquisition. In this work, we implement nonlinear processing and develop distributed beamformation to achieve fast acquisition over a 47-centimeter aperture. As a result, we achieve a 50-micrometer -6-decibel point spread function at 5 megahertz and resolve in-plane targets. A large volume scan of a human limb is completed in a few seconds, and in a 2-millimeter dorsal vein, the image intensity difference between the vessel center and surrounding tissue was ~50 decibels, facilitating three-dimensional reconstruction of the vasculature.

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

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

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.

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

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

Forward-viewing estimation of 3D blood flow velocity fields by intravascular ultrasound: Influence of the catheter on velocity estimation in stenoses

Coronary artery disease is the most common type of cardiovascular disease, affecting > 18 million adults, and is responsible for > 365 k deaths per year in the U.S. alone. Wall shear stress (WSS) is an emerging indicator of likelihood of plaque rupture in coronary artery disease, however, non-invasive estimation of 3-D blood flow velocity and WSS is challenging due to the requirement for high spatial resolution at deep penetration depths in the presence of significant cardiac motion. Thus we propose minimally-invasive imaging with a catheter-based, 3-D intravascular forward-viewing ultrasound (FV US) transducer and present experiments to quantify the effect of the catheter on flow disturbance in stenotic vessel phantoms with realistic velocities and luminal diameters for both peripheral (6.33 mm) and coronary (4.74 mm) arteries. An external linear array ultrasound transducer was used to quantify 2-D velocity fields in vessel phantoms under various conditions of catheter geometry, luminal diameter, and position of the catheter relative to the stenosis at a frame rate of 5000 frames per second via a particle imaging velocimetry (PIV) approach. While a solid catheter introduced an underestimation of velocity measurement by > 20% relative to the case without a catheter, the hollow catheter introduced < 10% velocity overestimation, indicating that a hollow catheter design allowing internal blood flow reduces hemodynamic disturbance. In addition, for both peripheral and coronary arteries, the hollow catheter introduced < 3% deviation in flow velocity at the minimum luminal area compared to the control case. Finally, an initial comparison was made between velocity measurements acquired using a low frequency, catheter-based, 3-D intravascular FV US transducer and external linear array measurements, with relative error < 12% throughout the region of interest for a flow rate of 150 mL/min. While further system development is required, results suggest intravascular ultrasound characterization of blood flow velocity fields in stenotic vessels could be feasible with appropriate catheter design.

Pulse Wave Imaging Coupled With Vector Flow Mapping: A Phantom, Simulation, and In Vivo Study

Pulse wave imaging (PWI) is an ultrasound imaging modality that estimates the wall stiffness of an imaged arterial segment by tracking the pulse wave propagation. The aim of the present study is to integrate PWI with vector flow imaging, enabling simultaneous and co-localized mapping of vessel wall mechanical properties and 2-D flow patterns. Two vector flow imaging techniques were implemented using the PWI acquisition sequence: 1) multiangle vector Doppler and 2) a cross-correlation-based vector flow imaging (CC VFI) method. The two vector flow imaging techniques were evaluated in vitro using a vessel phantom with an embedded plaque, along with spatially registered fluid structure interaction (FSI) simulations with the same geometry and inlet flow as the phantom setup. The flow magnitude and vector direction obtained through simulations and phantom experiments were compared in a prestenotic and stenotic segment of the phantom and at five different time frames. In most comparisons, CC VFI provided significantly lower bias or precision than the vector Doppler method ( ) indicating better performance. In addition, the proposed technique was applied to the carotid arteries of nonatherosclerotic subjects of different ages to investigate the relationship between PWI-derived compliance of the arterial wall and flow velocity in vivo. Spearman's rank-order test revealed positive correlation between compliance and peak flow velocity magnitude ( rs = 0.90 and ), while significantly lower compliance ( ) and lower peak flow velocity magnitude ( ) were determined in older (54-73 y.o.) compared with young (24-32 y.o.) subjects. Finally, initial feasibility was shown in an atherosclerotic common carotid artery in vivo. The proposed imaging modality successfully provided information on blood flow patterns and arterial wall stiffness and is expected to provide additional insight in studying carotid artery biomechanics, as well as aid in carotid artery disease diagnosis and monitoring.

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.

Vortex dynamics and transport phenomena in stenotic aortic models using Echo-PIV

Atherosclerosis is the most fatal cardiovascular disease. As disease progresses, stenoses grow inside the arteries blocking their lumen and altering blood flow. Analysing flow dynamics can provide a deeper insight on the stenosis evolution. In this work we combined Eulerian and Lagrangian descriptors to analyze blood flow dynamics and fluid transport in stenotic aortic models with morphology, mechanical and optical properties close to those of real arteries. To this end, vorticity, particle residence time (PRT), particle's final position (FP) and finite time Lyapunov's exponents (FTLE) were computed from the experimental fluid velocity fields acquired using ultrasonic particle imaging velocimetry (Echo-PIV). For the experiments, CT-images were used to create morphological realistic models of the descending aorta with 0%, 35% and 50% occlusion degree with same mechanical properties as real arteries. Each model was connected to a circuit with a pulsatile programmable pump which mimics physiological flow and pressure conditions. The pulsatile frequency was set to ≈0.9 Hz (55 bpm) and the upstream peak Reynolds number (Re) was changed from 1100 to 2000. Flow in the post-stenotic region was composed of two main structures: a high velocity jet over the stenosis throat and a recirculation region behind the stenosis where vortex form and shed. We characterized vortex kinematics showing that vortex propagation velocity increases withRe. Moreover, from the FTLE field we identified Lagrangian coherent structures (i.e. material barriers) that dictate transport behind the stenosis. The size and strength of those barriers increased withReand the occlusion degree. Finally, from the PRT and FP maps, we showed that independently ofRe, the same amount of fluid remains on the stenosis over more than a pulsatile period.

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.

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

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

Estimation of High Velocities in Synthetic Aperture Imaging: II: Experimental Investigation

The paper describes the performance of a new pulse sequence design and estimation approach for increasing the maximum detectable velocity in synthetic aperture (SA) velocity imaging. Measurements are conducted for conventional imaging for comparing the velocity range detectable by a directional Transverse Oscillation (TO) autocorrelation estimator to a new cross-correlation estimator. For conventional focused emissions a 192-elements, 3 MHz convex array transducer is used together with the SARUS experimental scanner on a flow rig at beam-to-flow angles of 60°, 70° and 90°. Here the new estimator always yields a higher precision, and the aliasing limit is increased by a factor 3. The new SA inter-spaced scheme was investigated using Field II simulations and SARUS measurements. A 3 MHz, 128-elements phased array was employed with a 5 virtual source emissions scheme for flow estimation and 15 emissions for B-mode imaging. The scheme was interleaved three times for a positive, negative, and positive transmission, so that non-linear pulse inversion also could be made. The experiments were conducted at three angles and for 4 different pulse repetition frequencies. A peak transverse velocity of 0.51 m/s could be estimated at fprf=450 Hz, translating to 5.6 m/s at fprf=5 kHz showing the theoretical increase of a factor 10 predicted in the accompanying theory paper.

Accuracy and Precision of a Plane Wave Vector Flow Imaging Method in the Healthy Carotid Artery

The objective of the study described here was to investigate the accuracy and precision of a plane wave 2-D vector flow imaging (VFI) method in laminar and complex blood flow conditions in the healthy carotid artery. The approach was to study (i) the accuracy for complex flow by comparing the velocity field from a computational fluid dynamics (CFD) simulation to VFI estimates obtained from the scan of an anthropomorphic flow phantom and from an in vivo scan; (ii) the accuracy for laminar unidirectional flow in vivo by comparing peak systolic velocities from VFI with magnetic resonance angiography (MRA); (iii) the precision of VFI estimation in vivo at several evaluation points in the vessels. The carotid artery at the bifurcation was scanned using both fast plane wave ultrasound and MRA in 10 healthy volunteers. The MRA geometry acquired from one of the volunteers was used to fabricate an anthropomorphic flow phantom, which was also scanned using the fast plane wave sequence. The same geometry was used in a CFD simulation to calculate the velocity field. Results indicated that similar flow patterns and vortices were estimated with CFD and VFI in the phantom for the carotid bifurcation. The root-mean-square difference between CFD and VFI was within 0.12 m/s for velocity estimates in the common carotid artery and the internal branch. The root-mean-square difference was 0.17 m/s in the external branch. For the 10 volunteers, the mean difference between VFI and MRA was -0.17 m/s for peak systolic velocities of laminar flow in vivo. The precision in vivo was calculated as the mean standard deviation (SD) of estimates aligned to the heart cycle and was highest in the center of the common carotid artery (SD = 3.6% for velocity magnitudes and 4.5° for angles) and lowest in the external branch and for vortices (SD = 10.2% for velocity magnitudes and 39° for angles). The results indicate that plane wave VFI measures flow precisely and that estimates are in good agreement with a CFD simulation and MRA.

Ultrasound Open Platforms for Next-Generation Imaging Technique Development

Open platform (OP) ultrasound systems are aimed primarily at the research community. They have been at the forefront of the development of synthetic aperture, plane wave, shear wave elastography, and vector flow imaging. Such platforms are driven by a need for broad flexibility of parameters that are normally preset or fixed within clinical scanners. OP ultrasound scanners are defined to have three key features including customization of the transmit waveform, access to the prebeamformed receive data, and the ability to implement real-time imaging. In this paper, a formative discussion is given on the development of OPs from both the research community and the commercial sector. Both software- and hardware-based architectures are considered, and their specifications are compared in terms of resources and programmability. Software-based platforms capable of real-time beamforming generally make use of scalable graphics processing unit architectures, whereas a common feature of hardware-based platforms is the use of field-programmable gate array and digital signal processor devices to provide additional on-board processing capacity. OPs with extended number of channels (>256) are also discussed in relation to their role in supporting 3-D imaging technique development. With the increasing maturity of OP ultrasound scanners, the pace of advancement in ultrasound imaging algorithms is poised to be accelerated.

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.

Spiral Flow Phantom for Ultrasound Flow Imaging Experimentation

As new ultrasound flow imaging methods are being developed, there is a growing need to devise appropriate flow phantoms that can holistically assess the accuracy of the derived flow estimates. In this paper, we present a novel spiral flow phantom design whose Archimedean spiral lumen naturally gives rise to multi-directional flow over all possible angles (i.e., from 0° to 360°). Developed using lost-core casting principles, the phantom geometry comprised a three-loop spiral (4-mm diameter and 5-mm pitch), and it was set to operate in steady flow mode (3 mL/s flow rate). After characterizing the flow pattern within the spiral vessel using computational fluid dynamics (CFD) simulations, the phantom was applied to evaluate the performance of color flow imaging (CFI) and high-frame-rate vector flow imaging. Significant spurious coloring artifacts were found when using CFI to visualize flow in the spiral phantom. In contrast, using vector flow imaging (least-squares multi-angle Doppler based on a three-transmit and three-receive configuration), we observed consistent depiction of flow velocity magnitude and direction within the spiral vessel lumen. The spiral flow phantom was also found to be a useful tool in facilitating demonstration of dynamic flow visualization based on vector projectile imaging. Overall, these results demonstrate the spiral flow phantom's practical value in analyzing the efficacy of ultrasound flow estimation methods.

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.