Elasticity imaging is a medical imaging modality that measures tissue elasticity to aid in diagnosis of certain diseases. Shear wave-based methods have been developed to perform elasticity measurements in soft tissue. Hybrid beamforming applies both the conventional spherical and axicon focusing to produce a beam for generating a shear wave with increased depth-of-field so that measurements can be made with a plane-like shear wave. We present our aperture design and beam optimization performed using Field II simulations. We varied the number of elements devoted to spherical and axicon focusing as well as the opening angle used for axicon focusing. We tested hybrid beamforming in three elastic phantoms and an excised kidney. We evaluated the shear wave speed measurements accuracy in the phantoms as well as the depth-of-field for each hybrid beam. We compared our results with those from using beams generated using spherical and axicon focusing. Our results show that hybrid beamforming is capable of producing a long narrow beam that performs well when among the 128 elements of transducer, 48 elements are allocated to each axicon portion and 32 elements for spherical aperture while the angle of the axicon aperture is set less than 20.
An energy functional approach for inverse characterization of material properties in elastodynamics
We present an inverse problem methodology based on the Error in Constitutive Equations (ECE) approach for the identification of material properties in the context of frequency-domain elastodynamics. In the ECE approach, we define a cost functional based on an energy norm that connects a set of kinematically admissible displacements and a set of dynamically admissible stresses. The set of kinematically admissible displacements is composed of fields that satisfy essential boundary conditions and possess sufficient regularity (i.e. smoothness). The set of dynamically admissible stresses is composed of fields that satisfy conservation of linear momentum and natural (i.e. traction) boundary conditions. The inverse problem is solved by finding material properties along with admissible displacement and stress fields such that the ECE functional is minimized. Experimental data is introduced in the formulation as a quadratic penalty term added to the ECE functional. The talk will focus on the reconstruction of elastic and viscoelastic properties in heterogeneous materials in the context of frequency-domain dynamics. Our findings indicate that ECE methods provide faster and more accurate results than conventional least-squares minimization. We will show numerical and experimental results that demonstrate the salient features of the method.
Quantifying viscoelasticity of boundary sensitive tissues using mechanical wave dispersion ultrasound vibrometry
The cardiovascular diseases such as atherosclerosis, coronary artery disease, hypertension and diastolic dysfunction have been associated with arterial stiffening and decreased ventricular compliance. Noninvasive techniques capable of quantifying elasticity and viscosity of cardiovascular tissues could facilitate early diagnosis and improve treatment. Here, we present a technique that uses focused ultrasound radiation force to excite mechanical waves in the tissue of interest and pulse echo to track the tissue deformation. Analysis of tissue deformation as a function of time using Fourier transforms allows calculation of the phase wave velocity dispersion (change of velocity as a function of frequency) of various modes of deformation. Continuum mechanics equations governing the motion of a viscoelastic plate and a tube are used to model the myocardial wall and arteries, respectively. Dispersion equations are derived for the two geometries and fit to the measured velocity dispersion to estimate tissue elasticity and viscosity. ECG-gated in vivo measurements of porcine myocardial and arterial elasticity and viscosity through a heart cycle are reported. The results show that both elasticity and viscosity increase during systole and decrease during diastole in the myocardium and arteries, consistent with underlying physiology.
Abstract 4047: Introducing a novel acoustic imaging tool for breast cancer detection
Abstract Objective: This paper introduces an alternative non-invasive and non-ionizing breast imaging tool called vibro-acoustography (VA). This method is based on low-frequency harmonic vibrations induced in the object by the radiation force of ultrasound. Contrast in Vibro-acoustography images represent several properties of the object, including the dynamic characteristics at the audio frequency and scattering and absorption properties at the ultrasound frequency. X-ray mammography and ultrasound are the main imaging modality used today for detection of breast lesions. X-ray mammography, however, is ionizing and its sensitivity is greatly reduced in dense breasts. Ultrasound images suffer from speckles or other or artifact that prevents the physician from delineating the mass lesion from the background tissue. Hence, VA may be used as a stand-alone imaging device or as a tool to provide supplementary information to other existing imaging modalities. Materials & Methods: Matching vibro-acoustography and mammography images of human breast are acquired using a combined mammography/vibro-acoustography machine. Images of various types of breast masses are investigated and their characteristics features are compared. Image features, such as lesion borders, interior texture, and speculation pattern, are studied. Results: Resulting images show soft tissue structures, breast lesions and microcalcifications within breast with high contrast, high resolution, and no speckles. The results have been verified using x-ray mammograms of the breast. Vibro-acoustography images demonstrate that calcified masses normally show sharp borders, benign non-calcified masses have smooth texture, and, in some cases, malignant masses are characterized with spiculation. Conclusion: The encouraging results which are based on the data collected from an ongoing study indicate that the above mentioned image features may be used for differentiation of breast masses and suggest that further development of vibro-acoustography technology may lead to a new clinical tool for various imaging applications. Work supported in part by Grant CA121579-01A1from NIH and Grant BCTR#0504550 from Susan G. Komen Breast Cancer Foundation. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 4047. doi:1538-7445.AM2012-4047
New vibration pulses are proposed to increase the power of shear waves induced by ultrasound radiation force in a tissue region with a preferred spectral distribution. The new pulses are sparsely sampled from an orthogonal frequency wave composed of several sinusoidal signals. Those sinusoidal signals have different frequencies and are orthogonal to each other. The phase and amplitude of each sinusoidal signal are adjusted to control the shape of the orthogonal frequency wave. Amplitude of the sinusoidal signal is increased as its frequency increases to compensate for higher loss at higher frequency in the tissue region. The new vibration pulses and detection pulses can be interleaved for array transducer applications. The experimental results show that the new vibration pulses significantly increases induced tissue vibration with the same peak ultrasound intensity, compared with the binary vibration pulses.
Phase velocities and attenuations of shear, Lamb, and Rayleigh waves in plate-like tissues submerged in a fluid (L)
In the past several decades, the fields of ultrasound and magnetic resonance elastography have shown promising results in noninvasive estimates of mechanical properties of soft tissues. These techniques often rely on measuring shear wave velocity due to an external or internal source of force and relating the velocity to viscoelasticity of the tissue. The mathematical relationship between the measured velocity and material properties of the myocardial wall, arteries, and other organs with non-negligible boundary conditions is often complicated and computationally expensive. A simple relationship between the Lamb–Rayleigh dispersion and the shear wave dispersion is derived for both the velocity and attenuation. The relationship shows that the shear wave velocity is around 20% higher than the Lamb–Rayleigh velocity and that the shear wave attenuation is about 20% lower than the Lamb–Rayleigh attenuation. Results of numerical simulations in the frequency range 0–500 Hz are presented.
Transthoracic in vivo measurements of porcine myocardial viscoelasticity using Lamb wave dispersion ultrasound vibrometry
Diastolic dysfunction is characterized by the stiffening of the left-ventricular myocardium. An imaging technique capable of quantifying viscoelasticity of the left-ventricle would aid clinicians in evaluating diastolic function. Lamb wave dispersion ultrasound vibrometry (LDUV) is an ultrasound-based method for measuring viscoelasticity of plate-like materials by exciting Lamb waves and measuring change of wave velocity as a function of frequency (dispersion). Lamb wave dispersion equation is fit to the dispersion data to estimate mechanical properties. We report in vivo transthoracic studies aimed at quantifying elasticity and viscosity of a beating porcine left-ventricle using LDUV. A single ultrasound probe operating at 5 MHz was used for impulse excitation in the myocardium and to track the wave motion. A Fourier space-time analysis of the motion was used to obtain Lamb wave velocity dispersion. ECG-gated transthoracic in vivo measurements of group velocity, elasticity, and viscosity throughout a single heart cycle were obtained. Group velocity, elasticity, and viscosity in the frequency range from 50 to 500 Hz increased by at least two-fold from diastole to systole, consistent with contraction and relaxation of the myocardium.
Two-dimensional shear elasticity imaging using comb-push acoustic radiation force and algebraic direct inversion of the motion differential equation
Tissue mechanical properties can be obtained by algebraic direct inversion (ADI) of the shear-wave motion differential equation, which is insensitive to reflection and geometry of the pushing beam. Shear waves normally generated by a focused ultrasound pushing beam have limited spatial extent in depth and are transient in time, leading to noisy and unstable ADI results. A comb-shape acoustic radiation force distribution can generate shear waves with longer spatial extent and time duration, facilitating more robust ADI. For 400 μs, a linear array transducer simultaneously transmitted four unfocused pushing beams (12 on elements for each pushing beam, 8 off elements between beams) into a tissue-mimicking phantom (∼1.5 kPa) with a cylindrical inclusion (∼4 kPa). The ultrasound system (Verasonics Inc.) then immediately switched to flash imaging mode (frame rate = 2 kHz, spatial resolution = 0.31 mm) to measure shear-wave displacements in a 38 mm by 30 mm domain that was used for ADI. The reconstructed 2D shear elasticity map provides accurate shear elasticity estimates (background = 1.5 kPa with 7% variance; inclusion = 3.9 kPa with 15% variance) and excellent contrast between the background and inclusion.
Measurement of elasticity of thin elastic layers with radiation force and wave propagation methods
In the field of tissue engineering, monitoring the elasticity of developing tissue cultures is important for evaluating their growth and maturation. In this work, we used wave propagation methods to measure the elasticity of thin tissue-mimicking layers to test the feasibility of measuring the elasticity of developing tissue cultures. We developed finite element and analytical models for investigating this problem and found good agreement between the results from both models. We used ultrasound radiation force to induce propagating waves, and the waves are measured using high frame rate ultrasound imaging at 10 kHz. The wave modes were identified and compared with an analytic model for Rayleigh wave propagation in a thin elastic layer bound to a solid substrate. The same gelatin mixture was used to make phantoms 1, 4, and 50 mm thick, where the phantom with 50 mm thickness was used as control and evaluated with shear wave imaging methods. The analytic model was used to fit data from experiments in the large block and thin layers of 1 and 4 mm thick, and the measured shear moduli were 22.7, 21.8, and 22.5 kPa, respectively. These results provide a validation for measurement of elasticity in thin elastic layers.
This paper describes the principles and initial experimental results of a new Doppler method called vibro-acoustic Doppler (VAD). VAD uses acoustic response of a moving object to a highly localized dynamic radiation force of ultrasound to measure the velocity of the object. The low frequency radiation force is exerted by two co-focused ultrasound beams with slightly different frequencies. The acoustic response of the object is detected by a hydrophone. A formula that describes the relation between the Doppler frequency shift of the emitted acoustic field and the velocity of the moving object is reported. To verify the theory, experiments are conducted on a moving string and a fluid flow phantom. Results show that the error in velocity measurement is less than 9.1% for either phantom. An advantage of this method over the traditional ultrasound Doppler is that velocity measurement with VAD is almost independent of the angle between the ultrasound beam and motion direction. It is shows that in the worst case, the error is 10.3% for a ±30 deg angle variation. Potential biomedical applications of VAD will be discussed.