Pixel-based approach to delay multiply and sum beamforming in combination with Wiener filter for improving ultrasound image quality.

Unified pixel-based (PB) beamforming has been implemented for ultrasound imaging, offering significant enhancements in lateral resolution compared to the conventional dynamic focusing. However, it still suffers from clutter and off-axis artifacts, limiting the contrast resolution. This paper proposes an efficient method to improve image quality by integrating filtered delay multiply and sum (F-DMAS) into the framework. This hybrid strategy incorporates the spatial coherence of the received data into the beamforming process to improve contrast resolution and clutter rejection in the generated image. We also integrate a Wiener filter to suppress the spatiotemporal spreading using signals echoed from a single scatterer at the transmit focus as a kernel for the deconvolution. The Wiener filter is applied to the received waveforms before performing the hybrid strategy. The Wiener filter is shown to reduce interference due to the interaction between the excitation pulse and the transfer functions of the transducer elements, thus benefiting the axial resolution of the generated images. We validate the proposed method and compare it with other beamforming strategies through a series of experiments, including simulation, phantom, and in vivo studies. The results show that our approach can substantially improve both spatial resolution and contrast over the unified PB algorithm, while still maintaining the good features of this beamformer. The simplicity and good performance of our method show its potential for use in clinical applications.

Improved ultrasound image quality with pixel-based beamforming using a Wiener-filter and a SNR-dependent coherence factor.

Pixel-based beamforming generates focused data by assuming that the waveforms received on a linear transducer array are composed of spherical pulses. It does not take into account the spatiotemporal spread in the data from the length of the excitation pulse or from the transfer functions of the transducer elements. As a result, these beamformers primarily have impacts on lateral, rather than axial, resolution. This paper proposes an efficient method to improve the axial resolution for pixel-based beamforming. We extend our field pattern analysis and show that the received waveforms should be passed through a Wiener filter before being used in the coherent pixel-based beamformer. This filter is designed based on signals echoed from a single scatterer at the transmit focus. The beamformer output is then combined with a coherence factor, that is adaptive to the signal-to-noise ratio, to improve the image contrast and suppress artifacts that have arisen during the filtering process. We validate the proposed method and compare it with other beamforming strategies using a series of experiments, including simulation, phantom and in vivo studies. It is shown to offer significant improvements in axial resolution and contrast over coherent pixel-based beamforming, as well as other spatial filters derived from synthetic aperture imaging. The method also demonstrates robustness to modeling errors in the experimental data. Overall, the imaging results show that the proposed approach has the potential to be of value in clinical applications.

Photoacoustic-ultrasonic dual-mode microscopy with local speed-of-sound estimation.

Synthetic aperture imaging and virtual point detection have been exploited to extend the depth of view of photoacoustic microscopy. The approach is commonly based on a constant assumed sound speed, which reduces image quality. We propose a new, to the best of our knowledge, self-adaptive technique to estimate the speed of sound when integrated with this hybrid strategy. It is accomplished through linear regression between the square of time of flight detected at individual virtual detectors and the square of their horizontal distances on the focal plane. The imaging results show our proposed method can significantly improve the lateral resolution, imaging intensity, and spatial precision for inhomogeneous tissue.

A Spatial Coherence Approach to Minimum Variance Beamforming for Plane-Wave Compounding.

A new approach to implement minimum variance distortionless response (MVDR) beamforming is introduced for coherent plane-wave compounding (CPWC). MVDR requires the covariance matrix of the incoming signal to be estimated and a spatial smoothing approximation is usually adopted to prevent this calculation from being underconstrained. In the new approach, we analyze MVDR as a spatial filter that decorrelates signals received at individual channels before summation. Based on the analysis, we develop two MVDR beamformers without using any spatial smoothing. First, MVDR weights are applied to the received signals after accumulating the data over transmits at different angles, while the second involves weighting the data collected in individual transmits and compounding over the transducer elements. In both cases, the covariance matrix is estimated using a set of slightly different combinations of the echo data. We show the sufficient statistic for this estimation that can be described by approximating the correlation among the backscattered ultrasound signals to their spatial coherence. Using the van Cittert-Zernike theorem, their statistical similarity is assessed by relating the spatial coherence to the profile of the source intensity. Both spatial-coherence-based MVDR beamformers are evaluated on data sets acquired from simulation, phantom, and in vivo studies. Imaging results show that they offer improvements over simple coherent compounding in terms of spatial and contrast resolutions. They also outperform other existing MVDR-based methods in the literature that are applied to CPWC.

Improvements to ultrasonic beamformer design and implementation derived from the task-based analytical framework.

The task-based framework, previously developed for beamformer comparison [Nguyen, Prager, and Insana, J. Acoust. Soc. Am. 140, 1048-1059 (2016)], is extended to design a new beamformer with potential applications in breast cancer diagnosis. The beamformer is based on a better approximation of the Bayesian strategy. It is a combination of the Wiener-filtered beamformer and an iterative process that adapts the generated image to specific features of the object. Through numerical studies, the new method is shown to outperform other beamformers drawn from the framework, but at an increase in computational cost. It requires a preprocessing step where the scattering field is segmented into regions with distinct statistical properties. Segmentation errors become a major limitation to the beamformer performance. All the beamformers under investigation are tested using data obtained from an instrumented ultrasound machine. They are implemented using a new time delay calculation, recently developed in the pixel-based beamforming studies presented here, which helps to overcome the challenge posed by the shift-variant nature of the imaging system. The efficacy of each beamformer is evaluated based on the quality of generated images in the context of the task-based framework. The in vitro results confirm the conclusions drawn from the simulations.

Ultrasound Pixel-Based Beamforming With Phase Alignments of Focused Beams.

We previously developed unified pixel-based (PB) beamforming to generate high-resolution sonograms, based on field pattern analysis. In this framework, we found that the transmit waveshape away from the focus could be characterized by two spherical pulses. These correspond to the maximal and minimal distances from the imaging point to the active aperture. The beamformer uses this model to select the highest energy signals from backscattered data. A spatiotemporal interpolation formula is used to provide a smooth transition in regions near the focal depth where there is no dominant reflected pulse. In this paper, we show that the unified PB approach is less robust at lower center frequencies. The interpolated data is suboptimal for a longer transmit waveshape. As a result, the spatial resolution at the focal depth is lower than that in other regions. By further exploring the field pattern, we propose a beamformer that is more robust to variations in beamwidth. The new method, named coherent PB beamforming, aligns and compounds the pulse data directly in the transition regions. In simulation and phantom studies, the coherent PB approach is shown to outperform the unified PB approach in spatial resolution. It helps regain optimal resolution at the focal depth while still maintaining good image quality in other regions. We also demonstrate the new method on in vivo data where its improvements over the unified PB method are demonstrated on scanned objects with a more complicated structure.

Minimum Variance Approaches to Ultrasound Pixel-Based Beamforming.

We analyze the principles underlying minimum variance distortionless response (MVDR) beamforming in order to integrate it into a pixel-based algorithm. There is a challenge posed by the low echo signal-to-noise ratio (eSNR) when calculating beamformer contributions at pixels far away from the beam centreline. Together with the well-known scarcity of samples for covariance matrix estimation, this reduces the beamformer performance and degrades the image quality. To address this challenge, we implement the MVDR algorithm in two different ways. First, we develop the conventional minimum variance pixel-based (MVPB) beamformer that performs the MVDR after the pixel-based superposition step. This involves a combination of methods in the literature, extended over multiple transmits to increase the eSNR. Then we propose the coherent MVPB beamformer, where the MVDR is applied to data within individual transmits. Based on pressure field analysis, we develop new algorithms to improve the data alignment and matrix estimation, and hence overcome the low-eSNR issue. The methods are demonstrated on data acquired with an ultrasound open platform. The results show the coherent MVPB beamformer substantially outperforms the conventional MVPB in a series of experiments, including phantom and in vivo studies. Compared to the unified pixel-based beamformer, the newest delay-and-sum algorithm in [1], the coherent MVPB performs well on regions that conform to the diffuse scattering assumptions on which the minimum variance principles are based. It produces less good results for parts of the image that are dominated by specular reflections.

A task-based analytical framework for ultrasonic beamformer comparison.

A task-based approach is employed to develop an analytical framework for ultrasound beamformer design and evaluation. In this approach, a Bayesian ideal-observer provides an idealized starting point and a way to measure information loss in practical beamformer designs. Different approximations of this ideal strategy are shown to lead to popular beamformers in the literature, including the matched filter, minimum variance (MV), and Wiener filter (WF) beamformers. Analysis of the approximations indicates that the WF beamformer should outperform the MV approach, especially in low echo signal-to-noise conditions. The beamformers are applied to five typical tasks from the BIRADS lexicon. Their performance is evaluated based on ability to discriminate idealized malignant and benign features. The numerical results show the advantages of the WF over the MV technique in general; although performance varies predictably in some contrast-limited tasks because of the model modifications required for the MV algorithm to avoid ill-conditioning.

High-Resolution Ultrasound Imaging With Unified Pixel-Based Beamforming.

This paper describes the development and evaluation of a new beamforming strategy based on pixel-based focusing for ultrasound linear array systems. We first implement conventional pixel-based beamforming in which the transmitted wave is assumed as spherical and diverging from the centre of the transmit subaperture. This assumed wave-shape is only valid within a limited angle on each side of the beam and this restricts the number of different subaperture positions from which data can be combined to improve image quality. By analyzing the field patterns, we propose a new unified pixel-based beamforming algorithm that better adapts to the non-spherical wave-shape of the transmit beam. This approach enables us to select the best-possible signal from each transducer waveform for data superposition. In simulations and a phantom study, we show that the unified pixel-based beamformer offers significant improvements in image quality compared to other delay-and-sum methods but at a higher computational cost. The new algorithm also demonstrates robust performance in a limited in vivo study. Overall, the results show that it is potentially of value in clinical applications.

Multidirectional scattering models for 3-dimensional ultrasound imaging.

An ultrasound image is created from backscattered echoes originating from both diffuse and directional scattering. It is potentially useful to separate these two components for the purpose of tissue characterization. This article presents several models for visualization of scattering fields on 3-dimensional (3D) ultrasound imaging. By scanning the same anatomy from multiple directions, we can observe the variation of specular intensity as a function of the viewing angle. This article considers two models for estimating the diffuse and specular components of the backscattered intensity: a modification of the well-known Phong reflection model and an existing exponential model. We examine 2-dimensional implementations and also propose novel 3D extensions of these models in which the probe is not constrained to rotate within a plane. Both simulation and experimental results show that improved performance can be achieved with 3D models.

Ultrasonic imaging of 3D displacement vectors using a simulated 2D array and beamsteering.

Most quasi-static ultrasound elastography methods image only the axial strain, derived from displacements measured in the direction of ultrasound propagation. In other directions, the beam lacks high resolution phase information and displacement estimation is therefore less precise. However, these estimates can be improved by steering the ultrasound beam through multiple angles and combining displacements measured along the different beam directions. Previously, beamsteering has only considered the 2D case to improve the lateral displacement estimates. In this paper, we extend this to 3D using a simulated 2D array to steer both laterally and elevationally in order to estimate the full 3D displacement vector over a volume. The method is tested on simulated and phantom data using a simulated 6-10MHz array, and the precision of displacement estimation is measured with and without beamsteering. In simulations, we found a statistically significant improvement in the precision of lateral and elevational displacement estimates: lateral precision 35.69µm unsteered, 3.70µm steered; elevational precision 38.67µm unsteered, 3.64µm steered. Similar results were found in the phantom data: lateral precision 26.51µm unsteered, 5.78µm steered; elevational precision 28.92µm unsteered, 11.87µm steered. We conclude that volumetric 3D beamsteering improves the precision of lateral and elevational displacement estimates.

A new method for the acquisition of ultrasonic strain image volumes.

This article presents a new method for acquiring three-dimensional (3-D) volumes of ultrasonic axial strain data. The method uses a mechanically-swept probe to sweep out a single volume while applying a continuously varying axial compression. Acquisition of a volume takes 15-20 s. A strain volume is then calculated by comparing frame pairs throughout the sequence. The method uses strain quality estimates to automatically pick out high quality frame pairs, and so does not require careful control of the axial compression. In a series of in vitro and in vivo experiments, we quantify the image quality of the new method and also assess its ease of use. Results are compared with those for the current best alternative, which calculates strain between two complete volumes. The volume pair approach can produce high quality data, but skillful scanning is required to acquire two volumes with appropriate relative strain. In the new method, the automatic quality-weighted selection of image pairs overcomes this difficulty and the method produces superior quality images with a relatively relaxed scanning technique.

A normalization method for axial-shear strain elastography.

The axial-shear strain distribution of soft tissue under load contains information useful for differentiating benign and malignant tumors. This paper describes a novel axial-shear strain normalization method. The algorithm builds on an existing normalization procedure for axial strain to map the shear strain values to the range [ -π/2, π/2]. The normalized shear data do not change sign with the direction of axial probe motion, and therefore can be time averaged without loss of information. Experiments in simulation, in vitro, and in vivo confirm the advantages of normalization. The proposed method is well suited to freehand strain imaging and enables the visualization of subtle slip patterns around inclusions.

Uniform precision ultrasound strain imaging.

Ultrasound strain imaging is becoming increasingly popular as a way to measure stiffness variation in soft tissue. Almost all techniques involve the estimation of a field of relative displacements between measurements of tissue undergoing different deformations. These estimates are often high resolution, but some form of smoothing is required to increase the precision, either by direct filtering or as part of the gradient estimation process. Such methods generate uniform resolution images, but strain quality typically varies considerably within each image, hence a trade-off is necessary between increasing precision in the low-quality regions and reducing resolution in the high-quality regions. We introduce a smoothing technique, developed from the nonparametric regression literature, which can avoid this trade-off by generating uniform precision images. In such an image, high resolution is retained in areas of high strain quality but sacrificed for the sake of increased precision in low-quality areas. We contrast the algorithm with other methods on simulated, phantom, and clinical data, for both 2-D and 3-D strain imaging. We also show how the technique can be efficiently implemented at real-time rates with realistic parameters on modest hardware. Uniform precision nonparametric regression promises to be a useful tool in ultrasound strain imaging.

A quality-guided displacement tracking algorithm for ultrasonic elasticity imaging.

Displacement estimation is a key step in the evaluation of tissue elasticity by quasistatic strain imaging. An efficient approach may incorporate a tracking strategy whereby each estimate is initially obtained from its neighbours' displacements and then refined through a localized search. This increases the accuracy and reduces the computational expense compared with exhaustive search. However, simple tracking strategies fail when the target displacement map exhibits complex structure. For example, there may be discontinuities and regions of indeterminate displacement caused by decorrelation between the pre- and post-deformation radio frequency (RF) echo signals. This paper introduces a novel displacement tracking algorithm, with a search strategy guided by a data quality indicator. Comparisons with existing methods show that the proposed algorithm is more robust when the displacement distribution is challenging.

The general properties including accuracy and resolution of linear filtering methods for strain estimation.

The vast majority of strain imaging systems applies linear filtering to estimate strain from displacement data. Methods such as piecewise-linear least squares regression and staggered strain estimation have come to be widely known and applied, but the properties of these estimators have rarely (or never) been compared quantitatively. Given their tractable properties, careful analysis of linear filters allows us to make numerous observations that are simple, yet valuable. We consider accuracy and resolving power, which raises the question of whether any particular filter offers the best possible accuracy at a given resolution. Our surprising results provide insight at two levels: They highlight general considerations affecting the type of filter that is appropriate for practical applications, and indicate promising avenues for further research.

Dynamic resolution selection in ultrasonic strain imaging.

Ultrasonic strain imaging promises to be a valuable tool in medical diagnostics. Reliability and ease-of-use have become important considerations. These depend on selection of appropriate imaging parameters. Two tasks are undertaken here. The tradeoff between resolution and estimation precision is examined closely to establish models for the relationships with imaging parameters and data properties. These models are then applied in a system that automatically sets the imaging parameters responsive to the data quality and the required estimation precision, so as to produce more meaningful images under varying scan conditions. The new system is applied to simulation, in vitro and in vivo data for validation. It reduces the complexity of the sonographer's role in strain imaging, and produces images of reliable quality even when the level of signal decorrelation varies throughout the ultrasound data.

Rotational motion in sensorless freehand three-dimensional ultrasound.

Freehand three-dimensional ultrasound is usually acquired with a position sensor attached to the ultrasound probe. However, position sensors can be expensive, obtrusive and difficult to calibrate. For this reason, there has been much research on alternative, image-based techniques, with in-plane motion tracked using conventional image registration methods, and out-of-plane motion inferred from the decorrelation between nearby B-scans. However, since out-of-plane motion is not the only source of decorrelation, image-based positions determined in this way suffer from cumulative drift errors. In this paper, we consider the effect of probe rotation on correlation and how this affects the position estimates. We then present a novel technique to compensate for out-of-plane rotations, by making use of orientation measurements from an unobtrusive sensor. Using simulations and in vitro experiments, we demonstrate that the technique is able to reduce the drift error in elevational positioning by 57% on average.

Comparison of freehand 3-D ultrasound calibration techniques using a stylus.

In a freehand 3-D ultrasound system, a probe calibration is required to find the rigid body transformation from the corner of the B-scan to the electrical center of the position sensor. The most intuitive way to perform such a calibration is by locating fiducial points in the scan plane directly with a stylus. The main problem of this approach is the difficulty in aligning the tip of the stylus with the scan plane. The thick beamwidth makes the tip of the stylus visible in the B-scan, even if the tip is not exactly at the elevational center of the scan plane. We present a novel stylus and phantom that simplify the alignment process for more accurate probe calibration. We also compare our calibration techniques with a range of styli. We show that our stylus and cone phantom are both simple in design and can achieve a point reconstruction accuracy of 2.2 mm and 1.8 mm, respectively, an improvement from 3.2 mm and 3.6 mm with the sharp and spherical stylus. The performance of our cone stylus and phantom lie between the state-of-the-art Z-phantom and Cambridge phantom, where accuracies of 2.5 mm and 1.7 mm are achieved.

An intelligent interface for freehand strain imaging.

We present a new, intelligent interface for freehand strain imaging, which has been designed to support clinical trials investigating the potential of ultrasonic strain imaging for diagnostic purposes across a broad range of target pathologies. The aim with this interface is to make scanning easier and to help clinicians learn the necessary scanning technique quickly, by providing real time feedback indicating the quality of the strain data as they are produced. The methods require a pixel-level indicator of estimation precision, which can be calculated in-line with strain estimation. This is exploited in novel approaches to normalisation, persistence and display. The effect of each component is indicated in the results with examples from in vitro and in vivo scanning. As well as providing real-time feedback, the images are easier to interpret because data at unacceptably low signal-to-noise ratios do not reach the display. Additionally, the level of noise in the displayed images is actually reduced compared with other methods that use the same strain estimates with the same level of persistence. The interface also considerably reduces the difficulty in producing volumes of strain data from freehand three-dimensional scans.