Ingvild Kinn Ekroll

Coherent Plane Wave Compounding for Very High Frame Rate Ultrasonography of Rapidly Moving Targets

Coherent plane wave compounding is a promising technique for achieving very high frame rate imaging without compromising image quality or penetration. However, this approach relies on the hypothesis that the imaged object is not moving during the compounded scan sequence, which is not the case in cardiovascular imaging. This work investigates the effect of tissue motion on retrospective transmit focusing in coherent compounded plane wave imaging (PWI). Two compound scan sequences were studied based on a linear and alternating sequence of tilted plane waves, with different timing characteristics. Simulation studies revealed potentially severe degradations in the retrospective focusing process, where both radial and lateral resolution was reduced, lateral shifts of the imaged medium were introduced, and losses in signal-to-noise ratio (SNR) were inferred. For myocardial imaging, physiological tissue displacements were on the order of half a wavelength, leading to SNR losses up to 35 dB, and reductions of contrast by 40 dB. No significant difference was observed between the different tilt sequences. A motion compensation technique based on cross-correlation was introduced, which significantly recovered the losses in SNR and contrast for physiological tissue velocities. Worst case losses in SNR and contrast were recovered by 35 dB and 27-35 dB, respectively. The effects of motion were demonstrated in vivo when imaging a rat heart. Using PWI, very high frame rates up to 463 fps were achieved at high image quality, but a motion correction scheme was then required.

Simultaneous quantification of flow and tissue velocities based on multi-angle plane wave imaging

A quantitative angle-independent 2-D modality for flow and tissue imaging based on multi-angle plane wave acquisition was evaluated. Simulations of realistic flow in a carotid artery bifurcation were used to assess the accuracy of the vector Doppler (VD) technique. Reduction in root mean square deviation from 27 cm/s to 6 cm/s and 7 cm/s to 2 cm/s was found for the lateral (vx) and axial (vz) velocity components, respectively, when the ensemble size was increased from 8 to 50. Simulations of a Couette flow phantom (vmax = 2.7 cm/s) gave promising results for imaging of slowly moving tissue, with root mean square deviation of 4.4 mm/s and 1.6 mm/s for the x- and z-components, respectively. A packet acquisition scheme providing both B-mode and vector Doppler RF data was implemented on a research scanner, and beamforming and further post-processing was done offline. In vivo results of healthy volunteers were in accordance with simulations and gave promising results for flow and tissue vector velocity imaging. The technique was also tested in patients with carotid artery disease. Using the high ensemble vector Doppler technique, blood flow through stenoses and secondary flow patterns were better visualized than in ordinary color Doppler. Additionally, the full velocity spectrum could be obtained retrospectively for arbitrary points in the image.

Robust angle-independent blood velocity estimation based on dual-angle plane wave imaging

Two-dimensional blood velocity estimation has shown potential to solve the angle-dependency of conventional ultrasound flow imaging. Clutter filtering, however, remains a major challenge for large beam-to-flow angles, leading to signal drop-outs and corrupted velocity estimates. This work presents and evaluates a compounding speckle tracking (ST) algorithm to obtain robust angle-independent 2-D blood velocity estimates for all beam-to-flow angles. A dual-angle plane wave imaging setup with full parallel receive beamforming is utilized to achieve high-frame-rate speckle tracking estimates from two scan angles, which may be compounded to obtain velocity estimates of increased robustness. The acquisition also allows direct comparison with vector Doppler (VD) imaging. Absolute velocity bias and root-mean-square (RMS) error of the compounding ST estimations were investigated using simulations of a rotating flow phantom with low velocities ranging from 0 to 20 cm/s. In a challenging region where the estimates were influenced by clutter filtering, the bias and RMS error for the compounding ST estimates were 11% and 2 cm/s, a significant reduction compared with conventional single-angle ST (22% and 4 cm/s) and VD (36% and 6 cm/s). The method was also tested in vivo for vascular and neonatal cardiac imaging. In a carotid artery bifurcation, the obtained blood velocity estimates showed that the compounded ST method was less influenced by clutter filtering than conventional ST and VD methods. In the cardiac case, it was observed that ST velocity estimation is more affected by low signal-to-noise (SNR) than VD. However, with sufficient SNR the in vivo results indicated that a more robust angle-independent blood velocity estimator is obtained using compounded speckle tracking compared with conventional ST and VD methods.

Combined vector velocity and spectral Doppler imaging for improved imaging of complex blood flow in the carotid arteries.

Color flow imaging and pulsed wave (PW) Doppler are important diagnostic tools in the examination of patients with carotid artery disease. However, measurement of the true peak systolic velocity is dependent on sample volume placement and the operators ability to provide an educated guess of the flow direction. Using plane wave transmissions and a duplex imaging scheme, we present an all-in-one modality that provides both vector velocity and spectral Doppler imaging from one acquisition, in addition to separate B-mode images of sufficient quality. The vector Doppler information was used to provide automatically calibrated (angle-corrected) PW Doppler spectra at every image point. It was demonstrated that the combined information can be used to generate spatial maps of the peak systolic velocity, highlighting regions of high velocity and the extent of the stenotic region, which could be used to automate work flow as well as improve the accuracy of measurement of true peak systolic velocity. The modality was tested in a small group (N = 12) of patients with carotid artery disease. PW Doppler, vector velocity and B-mode images could successfully be obtained from a single recording for all patients with a body mass index ranging from 21 to 31 and a carotid depth ranging from 16 to 28 mm.

Coherent compounding in doppler imaging

Coherent compounding can provide high frame rates and wide regions of interest for imaging of blood flow. However, motion will cause out-of-phase summation, potentially causing image degradation. In this work the impact of blood motion on SNR and the accuracy of Doppler velocity estimates are investigated. A simplified model for the compounded Doppler signal is proposed. The model is used to show that coherent compounding acts as a low-pass filter on the coherent compounding Doppler signal, resulting in negatively biased velocity estimates. Simulations and flow phantom experiments are used to quantify the bias and Doppler SNR for different velocities and beam-to-flow (BTF) angles. It is shown that the bias in the mean velocity increases with increasing beam-to-flow angle and/or blood velocity, whereas the SNR decreases; losses up to 4 dB were observed in the investigated scenarios. Further, a 2-D motion correction scheme is proposed based on multi-angle vector Doppler velocity estimates. For a velocity of 1.1 vNyq and a BTF angle of 75°, the bias was reduced from 30% to less than 4% in simulations. The motion correction scheme was also applied to flow phantom and in vivo recordings, in both cases resulting in a substantially reduced mean velocity bias and an SNR less dependent on blood velocity and direction.

2-D tracking doppler: a new method to limit spectral broadening in pulsed wave doppler

Transit time broadening is a major limitation in pulsed wave (PW) Doppler, especially when the angle between the flow direction and the ultrasound beam is large. The associated loss in frequency resolution may give severe overestimation of blood velocities, and finer details in the spectral display are lost. By using plane wave transmissions and parallel receive beamforming, multiple PW Doppler signals can be acquired simultaneously in a 2-D region. This enables tracking of the moving blood scatterers over a longer spatial distance to limit transit time broadening. In this work, the new method was tested using in vitro ultrasound recordings from a flow phantom, and in vivo recordings from a human carotid artery. The resulting 2-D tracking Doppler spectra showed significantly reduced spectral broadening compared with Doppler spectra generated by Welchs method. The reduction in spectral broadening was 4-fold when the velocity was 0.82 m/s and the beam-to-flow angle was 62°. A signal model was derived and the expected Doppler power spectra were calculated, showing good agreement with experimental data. Improved spectral resolution was shown for beam-to-flow angles between 40° and 82°.

Spectral doppler estimation utilizing 2-D spatial information and adaptive signal processing

The trade-off between temporal and spectral resolution in conventional pulsed wave (PW) Doppler may limit duplex/triplex quality and the depiction of rapid flow events. It is therefore desirable to reduce the required observation window (OW) of the Doppler signal while preserving the frequency resolution. This work investigates how the required observation time can be reduced by adaptive spectral estimation utilizing 2-D spatial information obtained by parallel receive beamforming. Four adaptive estimation techniques were investigated, the power spectral Capon (PSC) method, the amplitude and phase estimation (APES) technique, multiple signal classification (MUSIC), and a projection-based version of the Capon technique. By averaging radially and laterally, the required covariance matrix could successfully be estimated without temporal averaging. Useful PW spectra of high resolution and contrast could be generated from ensembles corresponding to those used in color flow imaging (CFI; OW = 10). For a given OW, the frequency resolution could be increased compared with the Welch approach, in cases in which the transit time was higher or comparable to the observation time. In such cases, using short or long pulses with unfocused or focused transmit, an increase in temporal resolution of up to 4 to 6 times could be obtained in in vivo examples. It was further shown that by using adaptive signal processing, velocity spectra may be generated without high-pass filtering the Doppler signal. With the proposed approach, spectra retrospectively calculated from CFI may become useful for unfocused as well as focused imaging. This application may provide new clinical information by inspection of velocity spectra simultaneously from several spatial locations.

An Extended Least Squares Method for Aliasing-Resistant Vector Velocity Estimation

An extended least squares method for robust, angle-independent 2-D vector velocity estimation using plane-wave ultrasound imaging is presented. The method utilizes a combination of least squares regression of Doppler autocorrelation estimates and block matching to obtain aliasing-resistant vector velocity estimates. It is shown that the aliasing resistance of the technique may be predicted using a single parameter, which is dependent on the selected transmit and receive steering angles. This parameter can therefore be used to design the aliasing-resistant transmit-receive setups. Furthermore, it is demonstrated that careful design of the transmit-receive steering pattern is more effective than increasing the number of Doppler measurements to obtain robust vector velocity estimates, especially in the presence of higher order aliasing. The accuracy and robustness of the method are investigated using the realistic simulations of blood flow in the carotid artery bifurcation, with velocities up to five times the Nyquist limit. Normalized root-mean-square (rms) errors are used to assess the performance of the technique. At -5 dB channel data blood SNR, rms errors in the vertical and horizontal velocity components were approximately 5% and 15% of the maximum absolute velocity, respectively. Finally, the in vivo feasibility of the technique is shown by imaging the carotid arteries of healthy volunteers.

Reconstruction of specular reflectors by iterative image source localization

A method is presented to reconstruct the geometry of specular reflectors with an ultrasonic array based on the image source principle. The ultrasonic beam is focused at a point in space emulating a point source within the medium. The transmitted wave interacts with the specular reflector and propagates back to the array as if it were generated by an image source. The reflected wave is analyzed with a sound source localization algorithm to estimate the image source location, and the reflector geometry is extracted using the mirror equation for spherical reflectors. The method is validated experimentally and its accuracy is studied. Under ideal conditions the method provides an accurate reconstruction of the position, orientation, and radius of curvature of specular reflectors, with errors Δr < 0.2 mm, Δα < 3°, and ΔR/R0 < 0.2, respectively. The method performs very well in the presence of high levels of thermal and speckle noise, with no degradation of the reconstruction as long as SNR(th) > -3 dB (signal-to-thermal-noise ratio) and SNR(sp) > 7 dB (signal-to-speckle-noise ratio). An iterative scheme based on the proposed method is presented to reconstruct the geometry of arbitrary reflectors by subdividing the reflector boundary into smaller segments. The iterative scheme is demonstrated both numerically and experimentally.

Investigations of spectral resolution and angle dependency in a 2-D tracking doppler method

An important source of error in velocity measurements from conventional pulsed wave (PW) Doppler is the angle used for velocity calibration. Because there are great uncertainties and interobserver variability in the methods used for Doppler angle correction in the clinic today, it is desirable to develop new and more robust methods. In this work, we have investigated how a previously presented method, 2-D tracking Doppler, depends on the tracking angle. A signal model was further developed to include tracking along any angle, providing velocity spectra which showed good agreement with both experimental data and simulations. The full-width at half-maximum (FWHM) bandwidth and the peak value of predicted power spectra were calculated for varying tracking angles. It was shown that the spectra have lowest bandwidth and maximum power when the tracking angle is equal to the beam-to-flow angle. This may facilitate new techniques for velocity calibration, e.g., by manually adjusting the tracking angle, while observing the effect on the spectral display. An in vitro study was performed in which the Doppler angles were predicted by the minimum FWHM and the maximum power of the 2-D tracking Doppler spectra for 3 different flow angles. The estimated Doppler angles had an overall error of 0.24° ± 0.75° when using the minimum FWHM. With an in vivo example, it was demonstrated that the 2-D tracking Doppler method is suited for measurements in a patient with carotid stenosis.