Investigation of parametric spectral estimation techniques for elasticity imaging

Medical ultrasound: imaging of soft tissue strain and elasticity

After X-radiography, ultrasound is now the most common of all the medical imaging technologies. For millennia, manual palpation has been used to assist in diagnosis, but it is subjective and restricted to larger and more superficial structures. Following an introduction to the subject of elasticity, the elasticity of biological soft tissues is discussed and published data are presented. The basic physical principles of pulse-echo and Doppler ultrasonic techniques are explained. The history of ultrasonic imaging of soft tissue strain and elasticity is summarized, together with a brief critique of previously published reviews. The relevant techniques-low-frequency vibration, step, freehand and physiological displacement, and radiation force (displacement, impulse, shear wave and acoustic emission)-are described. Tissue-mimicking materials are indispensible for the assessment of these techniques and their characteristics are reported. Emerging clinical applications in breast disease, cardiology, dermatology, gastroenterology, gynaecology, minimally invasive surgery, musculoskeletal studies, radiotherapy, tissue engineering, urology and vascular disease are critically discussed. It is concluded that ultrasonic imaging of soft tissue strain and elasticity is now sufficiently well developed to have clinical utility. The potential for further research is examined and it is anticipated that the technology will become a powerful mainstream investigative tool.

Theoretical limitations of the elastic wave equation inversion for tissue elastography

This article examines the theoretical limitations of the local inversion techniques for the measurement of the tissue elasticity. Most of these techniques are based on the estimation of the phase speed or the algebraic inversion of a one-dimensional wave equation. To analyze these techniques, the wave equation in an elastic continuum is revisited. It is proven that in an infinite medium, harmonic shear waves can travel at any phase speed greater than the classically known shear wave speed, mu/rho, by demonstrating this for a special case with cylindrical symmetry. Hence in addition to the mechanical properties of the tissue, the phase speed depends on the geometry of the wave as well. The elastic waves in an infinite cylindrical rod are studied. It is proven that multiple phase speeds can coexist for a harmonic wave at a single frequency. This shows that the phase speed depends not only on the mechanical properties of the tissue but also on its shape. The final conclusion is that the only way to avoid theoretical artifacts in the elastograms obtained by the local inversion techniques is to use the shear wave equation as expressed in the curl of the displacements, i.e., the rotations, for the inversion.

Angular strain estimation method for elastography

In the conventional cross-correlation-based strain estimation, there is a trade-off between the interpolation accuracy and the computational requirement. On the other hand, the autocorrelation-based method does not need interpolation, but it cannot estimate the wide range of displacements for elastography. We have developed a new strain estimator, called the angular strain estimation method, which does not need any interpolation and can estimate strain without restricting the range of displacements. The new method estimates strain utilizing complex correlation between correlated ultrasound signals from pre-and post-compression frames. From simulation and experiments, we found that the angular strain estimation method improves the accuracy and strain image quality compared to the conventional strain estimator using cross correlation with interpolation. Furthermore, it is computationally efficient and can be readily incorporated in ultrasound machines for rea -time elastography.

Comparison of shift estimation strategies in spectral elastography

This paper compares the performance of various spectral shift estimators for use in spectral elastography, namely, the normalized cross-correlation (NCC), sum squared difference (SSD) and sum absolute difference (SAD). Simulation and experimental results demonstrate that the spectral SSD-based elastographic method exhibits no marked difference in performance compared to the more computationally costly NCC-based approach, which has conventionally been the preferred estimator in spectral elastography. The spectral SAD-based strain estimator, despite being computationally less burdening, failed to exhibit performance comparable to that of the NCC- and SSD-based techniques. Furthermore, though spectral subsample estimation techniques using a cosine-fit interpolation method outperformed that of the parabolic-fit method in terms of both reduced bias errors and standard deviations, the latter was analyzed in this study due to computational simplicity. The role of spectral density was evaluated without and with parabolic-based subsample interpolation. Based on minimizing computational complexity, it is concluded that a (low density) spectral SSD strain estimator coupled with parabolic-based subsample estimation is the preferred choice for spectral elastography.

Analysis of a hybrid spectral strain estimation technique in elastography

Conventional spectral elastographic techniques estimate strain using cross-correlation methods. Despite promising results, decorrelation effects compromise the accuracy of these techniques and, subsequently, the tissue strain estimates. Since tissue compression in the time-domain corresponds to upscaling in the frequency-domain, decorrelation effects become more pronounced as tissue strains increase and are a fundamental concern in spectral cross-correlation elastography. In this paper, a two-stage hybrid spectral elastographic technique is introduced. For the first stage, an approximated spectral scaling factor (i.e. initial strain estimate) is employed to compensate for bandwidth broadening (due to tissue compression) between pre- and post-compression power spectra pairs. The second stage then estimates any residual strain information using spectral cross-correlation methods due to improper scaling factor selection in the first stage. This novel hybrid spectral elastographic technique was compared to both conventional spectral and adaptive temporal elastographic methods in simulation and experimentation. In addition to demonstrating enhancement in performance over the conventional spectral elastographic technique, the hybrid spectral-based method introduced in this paper is shown to outperform the adaptive temporal-based elastographic approach.