Contrast enhancement of adaptive ultrasound imaging using eigenspace-based minimum variance beamfoming

Constraint matrix projection minimum variance beamformer based on adaptive decomposition threshold for ultrasound imaging

In order to improve the efficiency of the eigenspace-based minimum variance (ESBMV) algorithm, a low-complexity minimum variance combined with eigenvalue decomposition (IBMV) algorithm is proposed. Firstly, the echo signals are transformed into beam domain by the discrete cosine transformation. Then, the eigenvalue decomposition of the sample covariance matrix is used to extract the signal subspace. The largest eigenvalue and its corresponding eigenvector can be extracted and other eigenvalues are taken the same value when the trace of sample covariance matrix remains unchanged. By these ways, the inverse operation of matrix can be simplified to the multiplication of vectors. In order to verify the validity of the proposed algorithm, Field II is introduced to test point targets and the cyst phantom. Experimental results indicate that the proposed IBMV with robustness to noises is obviously more efficient than ESBMV. Besides, the imaging quality of IBMV is better than delay and sum (DAS) method, minimum variance (MV) algorithm and beam domain minimum variance (BMV) method in terms of resolution and contrast ratio.

Adaptive beamforming methods for ultrasound imaging have been studied to improve image resolution and contrast. The most common approach is the minimum variance (MV) beamformer which minimizes the power of the beamformed output while maintaining the response from the direction of interest constant. The method achieves higher resolution and better contrast than the delay-and-sum (DAS) beamformer, but it suffers from high computational cost. This cost is mainly due to the computation of the spatial covariance matrix and its inverse, which requires O(L(3)) computations, where L denotes the subarray size. In this study, we propose a computationally efficient MV beamformer based on QR decomposition. The idea behind our approach is to transform the spatial covariance matrix to be a scalar matrix σI and we subsequently obtain the apodization weights and the beamformed output without computing the matrix inverse. To do that, QR decomposition algorithm is used and also can be executed at low cost, and therefore, the computational complexity is reduced to O(L(2)). In addition, our approach is mathematically equivalent to the conventional MV beamformer, thereby showing the equivalent performances. The simulation and experimental results support the validity of our approach.

In this paper, the adaptive, minimum variance (MV) beam-former is applied to medical ultrasound imaging. The significant resolution and contrast gain provided by the adaptive, minimum variance (MV) beamformer, introduces the possibility of plane wave (PW) ultrasound imaging. Data is obtained using Field II and a 7 MHz, 128-elements, linear array transducer with lambda/2-spacing. MV is compared to the conventional delay-and-sum (DS) beamformer with Boxcar and Hanning weights. Furthermore, the PW images are compared to the a conventional ultrasound image, obtained from a linear scan sequence. The four approaches, {Linear Scan, DS Boxcar, DS Hanning, MV}, have full width at half maximum of {0.82, 0.71, 1.28, 0.12} mm and peak side-lobe levels of {-40.1,-16.8,-34.4,-57.0} dB.

Although the minimum variance (MV) beamformer can provide enhancement in both resolution and contrast of ultrasound images when compared with conventional delay-and-sum (DAS) beamforming, its clinical application is limited by its sensitivity to phase aberrations. Several robust MV beamformers have been proposed, but present limitations when faced with second order phase aberration. Additionally, the information that can be obtained from phase aberration correction (PAC) methods is not utilized and the gain in robustness is usually accompanied by a loss in resolution or contrast. In the current work, we propose a new robust MV beamformer that uses a locally adaptive PAC (LAPAC-MV) in order to cope with the aformentioned difficulties. The method is based on a PAC estimation at every point in the image domain and a profile decision logic is included in the algorithm to take into account regions with low signal power. Preliminary results in full wave simulations and experimental phantom data are reported. In the resulting images, LAPAC-MV achieved enhancement in resolution of 0.03 mm, and reduction in sidelobe levels of around 8 dB when compared with precorrected DAS (DAS-PC) or precorrected MV (MV-PC). Contrast of LAPAC-MV is equivalent to MV-PC in cyst target simulations but the former showed improvements of approximately 1 dB in contrast in experiments with a 20 mm thick tissue aberrator.

In the recent years, various studies have been performed on adaptive beamforming especially for medical ultrasound imaging. In this paper, at first we applied the adaptive Minimum Variance (MV) beamforming to the simulated data derived from Field II. In contrast to conventional beamforming (DAS), MV adaptive beamforming can significantly improve the imaging resolution while retaining the contrast of the image. Then we combined MV with Coherence Factor (CF), by which the resolution and contrast of the medical ultrasound images are simultaneously increased. We implemented these beamformers in the MATLAB and survey the authenticity of them, then fulfil the beamformers with the Visual Studio codes. Those had the same results but execution times of Visual Studio codes are about 13 times lower than other. The major parameters which affect on the image quality are number of elements, height of elements, and sampling frequency. We investigated the effect of these parameters on the image quality by some simulated examples. The results illustrated that the reduction of the number of elements from M=96 to 24 and sampling frequency from 100MHz to 20MHz will degrade quality of the images. Nevertheless, the results of MV beamformer with M=24 and Fs=20MHz are comparable to DAS with M=96 and Fs=100MHz. We illustrated that if height of elements is equal or greater than 5mm, the off-plane targets have a negligible effect on the image quality but if the height of elements is in the same order of the width the off-plane targets have a significant role on the image quality. It should be taken into consideration in using of planar probs.

Minimum variance (MV) beamformers have been introduced in medical ultrasound imaging to improve image quality. In most cases, the MV beamformers have been investigated in terms of resolution improvement. However, the contrast-to-noise ratio (CNR) is also a clinically important metrics and gathers attention recently. In this study, we examined the diagonal loading parameter σ in MV beamforming and determined its appropriate value by evaluating image quality evaluation metrics including CNR. In order to further improve the image quality, a method for determining the value of σ based on the difference in statistical properties of received ultrasonic echo signals was also investigated. The phantom experimental results showed that the proposed method achieved a better CNR than the conventional MV beamformer while keeping resolution significantly better than that in delay-and-sum beamforming.

Delay-and-sum (DAS) beamformers, when applied to photoacoustic (PA) image reconstruction, produce strong sidelobes due to the absence of transmit focusing. Consequently, DAS PA images are often severely degraded by strong off-axis clutter. For preclinical in vivo cardiac PA imaging, the presence of these noise artifacts hampers the detectability and interpretation of PA signals from the myocardial wall, crucial for studying blood-dominated cardiac pathological information and to complement functional information derived from ultrasound imaging. In this article, we present PA subaperture processing (PSAP), an adaptive beamforming method, to mitigate these image degrading effects. In PSAP, a pair of DAS reconstructed images is formed by splitting the received channel data into two complementary nonoverlapping subapertures. Then, a weighting matrix is derived by analyzing the correlation between subaperture beamformed images and multiplied with the full-aperture DAS PA image to reduce sidelobes and incoherent clutter. We validated PSAP using numerical simulation studies using point target, diffuse inclusion and microvasculature imaging, and in vivo feasibility studies on five healthy murine models. Qualitative and quantitative analysis demonstrate improvements in PAI image quality with PSAP compared to DAS and coherence factor weighted DAS (DAS CF ). PSAP demonstrated improved target detectability with a higher generalized contrast-to-noise (gCNR) ratio in vasculature simulations where PSAP produces 19.61% and 19.53% higher gCNRs than DAS and DAS CF , respectively. Furthermore, PSAP provided higher image contrast quantified using contrast ratio (CR) (e.g., PSAP produces 89.26% and 11.90% higher CR than DAS and DAS CF in vasculature simulations) and improved clutter suppression.

In the article "Effect of element directivity on adaptive beamforming applied to high-frame-rate ultrasound" [ibid., vol. 62, no. 3, pp. 511-523, 2015], the contrast C was correctly defined. The authors note that therefore, 20 log <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">10</sub> 4 = 12.04 dB should be added to all the contrast values in the paper to obtain the correct contrast values.

One of the challenges in Super-Resolution Ultrasound Imaging (SRUI) based on Ultrasound Localization Microscopy (ULM) is that we need a large number of localized points to form the final image, so a significant number of frames is needed. In this study, we propose a method that significantly reduces the number of frames while eventually having an acceptable output image. Also, in this method, we no longer need to localize the center of the microbubbles, and we will create images only by the accumulation of binary images. This requires using a suitable adaptive beamformer for each frame to achieve the best possible quality of work. Investigation of beamformers in our other work proved that the Eigenspace-Based Minimum Variance (EIB-MV) beamformer with Coherence Factor (CF) coefficient was better than other beamformers in creating suitable binary frames. To compare the proposed method with the ULM method, parameters of SNR, PSNR and MSE were investigated. First, 100 frames were used to create the final images. Obtained results showed that the proposed method improves SNR and PSNR by 2 dB. Also, to investigate the effect of the number of frames on imaging parameters, only 8 frames were used to create the final image in the proposed method. By comparing the results of our proposed method with 8 frames and the ULM method with 100 frames, it was observed that the proposed method improves the SNR and PSNR by 0.8 dB relative to the ULM method.

In order to enhance the robustness and contrast in the minimum variance (MV) beamformer, adaptive diagonal loading method was proposed. The conventional diagonal loading technique has already been used in the MV beamformer, but has the drawback that its level is specified by predefined parameter and without consideration of input-data. To alleviate this problem, the level of diagonal loading was computed appropriately and automatically from the given data by shrinkage method in the proposed adaptive diagonal loaded beamformer. The performance of the proposed beamformer was tested on the simulated point target and cyst phantom was obtained using Field II. In the point target simulation, it is shown that the proposed method has higher lateral resolution than the conventional delay-and-sum beamformer and could be more robust in estimating the amplitude peak than the MV beamformer when acoustic velocity error exists. In the cyst phantom simulation, the proposed beamformer has shown that it achieves an improvement in contrast ratio and without distorting the edges of cyst.

A Low-complexity Minimum-variance Beamformer Based on Orthogonal Decomposition of the Compounded Subspace.

In photoacoustic imaging, delay-and-sum (DAS) beamformer is a common beamforming algorithm having a simple implementation. However, it results in a poor resolution and high sidelobes. To address these challenges, a new algorithm namely delay-multiply-and-sum (DMAS) was introduced having lower sidelobes compared to DAS. To improve the resolution of DMAS, a beamformer is introduced using minimum variance (MV) adaptive beamforming combined with DMAS, so-called minimum variance-based DMAS (MVB-DMAS). It is shown that expanding the DMAS equation results in multiple terms representing a DAS algebra. It is proposed to use the MV adaptive beamformer instead of the existing DAS. MVB-DMAS is evaluated numerically and experimentally. In particular, at the depth of 45mm MVB-DMAS results in about 31, 18, and 8dB sidelobes reduction compared to DAS, MV, and DMAS, respectively. The quantitative results of the simulations show that MVB-DMAS leads to improvement in full-width-half-maximum about 96%, 94%, and 45% and signal-to-noise ratio about 89%, 15%, and 35% compared to DAS, DMAS, MV, respectively. In particular, at the depth of 33mm of the experimental images, MVB-DMAS results in about 20dB sidelobes reduction in comparison with other beamformers.

Minimum variance beamforming (MVBF) is an adaptive beamforming technique, which aims to improve the lateral resolution by computing and applying signal-dependent apodization rather than predetermined apodization as typically done in conventional delay-and-sum (DAS) beamforming. Although studies have shown that the improvement in lateral resolution associated with MVBF is significant, the axial resolution remains unaffected. In this work, we combine MVBF and spiking deconvolution to improve both lateral and axial resolutions in synthetic aperture ultrasound imaging. We implement our new method and evaluate its performance using experimental datasets from a tissue-mimicking phantom. Our results show that our new method yields improved axial and lateral resolutions as well as image contrast.

In contrast with the common delay-and-sum (DAS) beamformer, adaptive beamformers can be used to increase the resolution of acoustic measurement systems without sacrificing sidelobe levels. The non-linear nature of adaptive beamforming methods may cause changes to measured scattered amplitude distributions. For example, large differences in speckle statistics have been observed in some cases when using the adaptive minimum variance distortionless response (MVDR) beamformer in medical ultrasound, unless preventative measures were implemented. Low complexity adaptive (LCA) beamforming techniques have shown similar, but smaller, differences. Recently, experiments have been conducted in the Oslofjord, Norway, to determine whether low-complexity adaptive (LCA) beamforming techniques can improve the performance of bathymetric echosounders. When using LCA with these types of sonar systems, it is important to understand any changes to the backscatter statistics, as these changes may affect, e.g., seafloor classification methods. In this talk we will present comparisons between backscatter statistics estimated using a weighted delay-and-sum beamformer and those estimated using a low complexity adaptive beamformer. We will discuss any differences seen, possible mechanisms causing these differences and potential remedies. [Research supported by The Norwegian Research Council and Kongsberg Maritime. The first author is employed in Kongsberg Maritime.]