Min Rao

Spatial-angular compounding for elastography using beam steering on linear array transducers

Spatial-angular compounding is a new technique that enables the reduction of noise artifacts in ultrasound elastography. Under this method, compounded elastograms are obtained from a spatially weighted average of local strain estimated from radio frequency (rf) echo signals acquired at different insonification angles. In previous work, the acquisition of the rf signals was performed through the lateral translation of a phased-array transducer. Clinical applications of angular compounding would, however, require the utilization of beam steering on linear-array transducers to obtain angular data sets, which is more efficient than translating phased-array transducers. In this article, we investigate the performance of angular compounding for elastography by using beam steering on a linear-array transducer. Quantitative experimental results demonstrate that spatial angular compounding provides significant improvement in both the elastographic signal-to-noise ratio and the contrast-to-noise ratio. For the linear array transducer used in this study, the optimum angular increment is around 1.5 degrees-3.75 degrees, and the maximum angle that can be used in angular compounding should not exceed 10 degrees.

Axial-Shear Strain Imaging for Differentiating Benign and Malignant Breast Masses

Axial strain imaging has been utilized for the characterization of breast masses for over a decade; however, another important feature namely the shear strain distribution around breast masses has only recently been used. In this article, we examine the feasibility of utilizing in vivo axial-shear strain imaging for differentiating benign from malignant breast masses. Radio-frequency data was acquired using a VFX 13-5 linear array transducer on 41 patients using a Siemens SONOLINE Antares real-time clinical scanner at the University of Wisconsin Breast Cancer Center. Free-hand palpation using deformations of up to 10% was utilized to generate axial strain and axial-shear strain images using a two-dimensional cross-correlation algorithm from the radio-frequency data loops. Axial-shear strain areas normalized to the lesion size, applied strain and lesion strain contrast was utilized as a feature for differentiating benign from malignant masses. The normalized axial-shear strain area feature estimated on eight patients with malignant tumors and 33 patients with fibroadenomas was utilized to demonstrate its potential for lesion differentiation. Biopsy results were considered the diagnostic standard for comparison. Our results indicate that the normalized axial-shear strain area is significantly larger for malignant tumors compared with benign masses such as fibroadenomas. Axial-shear strain pixel values greater than a specified threshold, including only those with correlation coefficient values greater than 0.75, were overlaid on the corresponding B-mode image to aid in diagnosis. A scatter plot of the normalized area feature demonstrates the feasibility of developing a linear classifier to differentiate benign from malignant masses. The area under the receiver operator characteristic curve utilizing the normalized axial-shear strain area feature was 0.996, demonstrating the potential of this feature to noninvasively differentiate between benign and malignant breast masses.

Spatial Angular Compounding for Elastography without the Incompressibility Assumption

Spatial-angular compounding is a new technique that enables the reduction of noise artifacts in ultrasound elastography. Previous results using spatial angular compounding, however, were based on the use of the tissue incompressibility assumption. Compounded elastograms were obtained from a spatially-weighted average of local strain estimated from radiofrequency echo signals acquired at different insonification angles. In this paper, we present a new method for reducing the noise artifacts in the axial strain elastogram utilizing a least-squares approach on the angular displacement estimates that does not use the incompressibility assumption. This method produces axial strain elastograms with higher image quality, compared to noncompounded axial strain elastograms, and is referred to as the least-squares angular-compounding approach for elastography. To distinguish between these two angular compounding methods, the spatial-angular compounding with angular weighting based on the tissue incompressibility assumption is referred to as weighted compounding. In this paper, we compare the performance of the two angular-compounding techniques for elastography using beam steering on a linear-array transducer. Quantitative experimental results demonstrate that least-squares compounding provides comparable but smaller improvements in both the elastographic signal-to-noise ratio and the contrast-to-noise ratio, as compared to the weighted-compounding method. Ultrasound simulation results suggest that the least-squares compounding method performs better and provide accurate and robust results when compared to the weighted compounding method, in the case where the incompressibility assumption does not hold.

Correlation analysis of three-dimensional strain imaging using ultrasound two-dimensional array transducers

Two-dimensional (2D) transducer arrays represent a promising solution for implementing real time three-dimensional (3D) ultrasound elastography. 2D arrays enable electronic steering and focusing of ultrasound beams throughout a 3D volume along with improved slice thickness performance when compared to one-dimensional (1D) transducer arrays. Therefore, signal decorrelation caused by tissue motion in the elevational (out-of-plane) direction needs to be considered. In this paper, a closed form expression is derived for the correlation coefficient between pre- and postdeformation ultrasonic radio frequency signals. Signal decorrelation due to 3D motion of scatterers within the ultrasonic beam has been considered. Computer simulations are performed to corroborate the theoretical results. Strain images of a spherical inclusion phantom generated using 1D and 2D array transducers are obtained using a frequency domain simulation model. Quantitative image quality parameters, such as the signal-to-noise and contrast-to-noise ratios obtained using 1D, 2D, and 3D motion tracking algorithms, are compared to evaluate the performance with the 3D strain imaging system. The effect of the aperture size for 2D arrays and other factors that affect signal decorrelation are also discussed.

Shear strain imaging using shear deformations

In this article we investigate the generation of shear strain elastograms induced using a lateral shear deformation. Ultrasound simulation and experimental results demonstrate that the shear strain elastograms obtained under shear deformation exhibit significant differences between bound and unbound inclusions in phantoms, when compared to shear strain images induced upon an axial compression. A theoretical model that estimates the decorrelation between pre- and postdeformation radio frequency signals, as a function of extent of shear deformation, is also developed. Signal-to-noise ratios of shear strain elastograms obtained at different shear angles are investigated theoretically and verified using ultrasound simulations on a uniformly elastic phantom. For the simulation and experiment, a two-dimensinal block-matching-based algorithm is used to estimate the axial and lateral displacement. Shear strains are obtained from the displacement vectors using a least-squares strain estimator. Our results indicate that the signal-to-noise ratio (SNR) of shear strain images increases to reach a maximum and saturates, and then decreases with increasing shear angle. Using typical system parameters, the maximum achievable SNR for shear strain elastography is around 8 (18 dB), which is comparable to conventional axial strain elastography induced by axial compression. Shear strain elastograms obtained experimentally using single inclusion tissue-mimicking phantoms with both bound and unbound inclusions (mimicking cancerous masses and benign fibroadenomas, respectively) demonstrate the characteristic differences in the depiction of these inclusions on the shear strain elastograms.

Correlation analysis of the beam angle dependence for elastography

Signal decorrelation is a major source of error in the displacements estimated using correlation techniques for elastographic imaging. Previous papers have addressed the variation in the correlation coefficient as a function of the applied compression for a finite window size and an insonification angle of zero degrees. The recent use of angular beam-steered radio-frequency echo signals for spatial angular compounding and shear strain estimation have demonstrated the need for understanding signal decorrelation artifacts for data acquired at different beam angles. In this paper, we provide both numerical and closed form theoretical solutions of the correlation between pre- and post-compression radio-frequency echo signals acquired at a specified beam angle. The expression for the correlation coefficient obtained is a function of the beam angle and the applied compression for a finite duration window. Accuracy of the theoretical results is verified using tissue-mimicking phantom experiments on a uniformly elastic phantom using beam-steered data acquisitions on a linear array transducer. The theory predicts a faster decorrelation with changes in the beam or insonification angle for longer radio-frequency echo signal segments and at deeper locations in the medium. Theoretical results provide useful information for improving angular compounding and shear strain estimation techniques for elastography.

Estimation of the Optimal Maximum Beam Angle and Angular Increment for Normal and Shear Strain Estimation

In the current practice of ultrasound elastography, only the axial component of the displacement vector is estimated and used to produce strain images. A method was recently proposed by our group to estimate both the axial and lateral components of a displacement vector using RF echo signal data acquired along multiple angular insonification directions of the ultrasound beam. Previous work has demonstrated that it is important to choose appropriate values for the maximum beam angle and angular increment to achieve optimal performance with this technique. In this paper, we present error propagation analysis using the least-square fitting process for the optimization of the angular increment and the maximum beam steered angle. Ultrasound simulations are performed to corroborate the theoretical prediction of the optimal values for the maximum beam angle and angular increment. Selection of the optimal parameters depends on system parameters, such as center frequency and aperture size. For typical system parameters, the optimal maximum beam angle is around 10deg for axial strain estimation and around 15deg for lateral strain estimation. The optimal angular increment is around 4deg -6deg, which indicates that only five to seven beam angles are required for this strain-tensor estimation technique.

Correlation Analysis for Angular Compounding in Strain Imaging

Spatial angular compounding for elastography is a new technique that enables the reduction of noise artifacts in elastograms. This technique is most effective when the angular strain estimates to be averaged or compounded are uncorrelated. In this paper, we present a theoretical analysis of the correlation between pre- and postcompression radio-frequency echo signals acquired from the same location but at different beam insonification angles. The accuracy of the theoretical results is verified using radio- frequency pre- and postcompression echo signals acquired using a real-time clinical scanner on tissue-mimicking uniformly elastic and homogenous phantoms. The theory predicts an increased signal decorrelation with an increase in the beam-steered insoniflcation angle as the applied strain increases and for increasing depths in the medium. Theoretical results provide useful information regarding the correlation of the angular strain estimates obtained from different beam angles that helps in finding optimum compounding schemes for elastography.

Simulation of ultrasound two-dimensional array transducers using a frequency domain model.

Ultrasound imaging with two-dimensional (2D) arrays has garnered broad interest from scanner manufacturers and researchers for real time three-dimensional (3D) applications. Previously the authors described a frequency domain B-mode imaging model applicable for linear and phased array transducers. In this paper, the authors extend this model to incorporate 2D array transducers. Further approximations can be made based on the fact that the dimensions of the 2D array element are small. The model is compared with the widely used ultrasound simulation program FIELD II, which utilizes an approximate form of the time domain impulse response function. In a typical application, errors in simulated RF waveforms are less than 4% regardless of the steering angle for distances greater than 2 cm, yet computation times are on the order of 1/35 of those incurred using FIELD II. The 2D model takes into account the effects of frequency-dependent attenuation, backscattering, and dispersion. Modern beam-forming techniques such as apodization, dynamic aperture, dynamic receive focusing, and 3D beam steering can also be simulated.

TU‐E‐201C‐02: Breast Mass Differentiation Using Axial Shear Strain Imaging

Purpose: Since cancers infiltrate surrounding normal tissue, evoke a desmoplastic, scirrhous reaction and become firmly attached to background tissue, they tend to be less mobile than benign masses like fibroadenomas. We test the feasibility of using in‐vivo axial‐shear strain features to determine if the bonding of masses to background tissue can differentiate benign from malignant. Method and Materials: Radiofrequency data was acquired using a VFX 13‐5 linear array transducer using a Siemens SONOLINE Antares real‐time clinical scanner. Data were acquired for 41 biopsy‐proven breast masses (8 malignant tumors and 33 benign fibroadenomas). Free‐hand palpation using the transducer and deformations of up to 10% was utilized for data acquisition. A two‐dimensional cross‐correlation algorithm was used to generate axial strain and axial‐shear strain images. Axial‐shear strain values normalized to the breast mass dimensions, applied strain and strain contrast were utilized to calculate the feature called the “normalized axial‐shear strain area”, for differentiating benign from malignant masses. Results: The normalized axial‐shear strain area is significantly larger for malignant masses when compared to benign fibroadenomas. Scatter plots of the normalized axial‐shear strain area demonstrates the feasibility for differentiating benign from malignant masses. Receiver operator characteristic analysis demonstrates the improvement in the classification obtained using the normalized axial‐shear strain area. The area under the ROC curve of 0.996 suggests that this feature can effectively differentiate malignant tumors from the benign masses. Conclusions: Axial‐shear strain images may provide important additional information which along with currently utilized axial strain and B‐mode images may improve differentiation of benign and malignant breast masses. This work was supported by Komen grant BCTR0601153 and NTH‐NCI grants R01CA112192‐S103, R01CA100373 and R01 CA111289.