Model-based elastography: a survey of approaches to the inverse elasticity problem.

Elastography is emerging as an imaging modality that can distinguish normal versus diseased tissues via their biomechanical properties. This paper reviews current approaches to elastography in three areas--quasi-static, harmonic and transient--and describes inversion schemes for each elastographic imaging approach. Approaches include first-order approximation methods; direct and iterative inversion schemes for linear elastic; isotropic materials and advanced reconstruction methods for recovering parameters that characterize complex mechanical behavior. The paper's objective is to document efforts to develop elastography within the framework of solving an inverse problem, so that elastography may provide reliable estimates of shear modulus and other mechanical parameters. We discuss issues that must be addressed if model-based elastography is to become the prevailing approach to quasi-static, harmonic and transient elastography: (1) developing practical techniques to transform the ill-posed problem with a well-posed one; (2) devising better forward models to capture the complex mechanical behavior of soft tissues and (3) developing better test procedures to evaluate the performance of modulus elastograms.

Investigating the impact of spatial priors on the performance of model-based IVUS elastography.

This paper describes methods that provide pre-requisite information for computing circumferential stress in modulus elastograms recovered from vascular tissue-information that could help cardiologists detect life-threatening plaques and predict their propensity to rupture. The modulus recovery process is an ill-posed problem; therefore, additional information is needed to provide useful elastograms. In this work, prior geometrical information was used to impose hard or soft constraints on the reconstruction process. We conducted simulation and phantom studies to evaluate and compare modulus elastograms computed with soft and hard constraints versus those computed without any prior information. The results revealed that (1) the contrast-to-noise ratio of modulus elastograms achieved using the soft prior and hard prior reconstruction methods exceeded those computed without any prior information; (2) the soft prior and hard prior reconstruction methods could tolerate up to 8% measurement noise, and (3) the performance of soft and hard prior modulus elastograms degraded when incomplete spatial priors were employed. This work demonstrates that including spatial priors in the reconstruction process should improve the performance of model-based elastography, and the soft prior approach should enhance the robustness of the reconstruction process to errors in the geometrical information.

Imaging the elastic properties of tissue: the 20 year perspective.

After 20 years of innovation in techniques that specifically image the biomechanical properties of tissue, the evolution of elastographic imaging can be viewed from its infancy, through a proliferation of approaches to the problem to incorporation on research and then clinical imaging platforms. Ultimately this activity has culminated in clinical trials and improved care for patients. This remarkable progression represents a leading example of translational research that begins with fundamentals of science and engineering and progresses to needed improvements in diagnostic and monitoring capabilities applied to major categories of disease, surgery and interventional procedures. This review summarizes the fundamental principles, the timeline of developments in major categories of elastographic imaging, and concludes with recent results from clinical trials and forward-looking issues.

Effects of frequency- and direction-dependent elastic materials on linearly elastic MRE image reconstructions.

The mechanical model commonly used in magnetic resonance elastography (MRE) is linear elasticity. However, soft tissue may exhibit frequency- and direction-dependent (FDD) shear moduli in response to an induced excitation causing a purely linear elastic model to provide an inaccurate image reconstruction of its mechanical properties. The goal of this study was to characterize the effects of reconstructing FDD data using a linear elastic inversion (LEI) algorithm. Linear and FDD phantoms were manufactured and LEI images were obtained from time-harmonic MRE acquisitions with variations in frequency and driving signal amplitude. LEI responses to artificially imposed uniform phase shifts in the displacement data from both purely linear elastic and FDD phantoms were also evaluated. Of the variety of FDD phantoms considered, LEI appeared to tolerate viscoelastic data-model mismatch better than deviations caused by poroelastic and anisotropic mechanical properties in terms of visual image contrast. However, the estimated shear modulus values were substantially incorrect relative to independent mechanical measurements even in the successful viscoelastic cases and the variations in mean values with changes in experimental conditions associated with uniform phase shifts, driving signal frequency and amplitude were unpredictable. Overall, use of LEI to reconstruct data acquired in phantoms with FDD material properties provided biased results under the best conditions and significant artifacts in the worst cases. These findings suggest that the success with which LEI is applied to MRE data in tissue will depend on the underlying mechanical characteristics of the tissues and/or organs systems of clinical interest.

Interstitial fluid pressure in soft tissue as a result of an externally applied contact pressure.

Manipulation of interstitial fluid pressure (IFP) has a clinical potential when used in conjunction with near-infrared spectroscopy for the detection of breast cancer. In order to better interpret how the applied pressure alters the vascular space and interstitial water volumes in breast tissue, a study on tissue-mimicking, gelatin phantoms was carried out to mimic the translation of external force into internal pressures. A complete set of three-dimensional (3D) pressure maps were obtained for the interior volumes of phantoms as an external force of 10 mmHg was applied, using mixtures of elastic moduli 19 and 33 kPa to simulate adipose and fibroglandular values of breast tissue. Corresponding linear elastic finite element analysis (FEA) cases were formulated. Shear stress, nonlinear mechanical properties, gravity and tissue geometry were all observed to contribute to internal pressure distribution, with surface shear stresses increasing internal pressures near the surface to greater than twice the applied external pressure. Average pressures by depth were predicted by the linear elastic FEA models. FEA models were run for cases mimicking a 93 kPa tumor inclusion within regions of adipose, fibroglandular tissue, and a composite of the two tissue types to illustrate the localized high fluid pressures caused by a tumor when an external force is applied. The conclusion was that external contact forces can generate potentially clinically useful fluid pressure magnitudes in regions of sharp effective elastic modulus gradients, such as tumor boundaries.

Morphological and mechanical information of coronary arteries obtained with intravascular elastography; feasibility study in vivo.

AIMS: Plaque composition is a major determinant of coronary related clinical syndromes. In vitro experiments on human coronary and femoral arteries have demonstrated that different plaque types were detectable with intravascular ultrasound elastography. The aim of this study was to investigate the feasibility of applying intravascular elastography during interventional catheterization procedures. METHODS AND RESULTS: Data were acquired in patients (n=12) during PTCA procedures with an EndoSonics InVision echoapparatus equipped with radiofrequency output. The systemic pressure was used to strain the tissue, and the strain was determined using cross-correlation analysis of sequential frames. A likelihood function was determined to obtain the frames with minimal motion of the catheter in the lumen, since motion of the catheter prevents reliable strain estimation. Minimal motion was observed near end-diastole. Reproducible strain estimates were obtained within one pressure cycle and over several pressure cycles. Validation of the results was limited to the information provided by the echogram. Strain in calcified material (0.20%+/-0.07) was lower (P<0.001) than in non-calcified tissue (0.51%+/-0.20). CONCLUSION: In vivo intravascular elastography is feasible. Significantly higher strain values were found in non-calcified plaques than in calcified plaques.

A freehand elastographic imaging approach for clinical breast imaging: system development and performance evaluation.

A prototype freehand elastographic imaging system has been developed for clinical breast imaging. The system consists of a fast data acquisition system, which is able to capture sequences of intermediate frequency echo frames at full frame rate from a commercial ultrasound scanner whilst the breast is deformed using hand-induced transducer motion. Two-dimensional echo tracking was used in combination with global distortion compensation and multi-compression averaging to minimise decorrelation noise incurred when stress is applied using hand-induced transducer motion. Experiments were conducted on gelatine phantoms to evaluate the quality of elastograms produced using the prototype system relative to those produced using mechanically induced transducer motion. The strain sensitivity and contrast-to-noise ratio of freehand elastograms compared favourably with elastograms produced using mechanically induced transducer motion. However, better dynamic range and signal-to-noise ratio was achieved when elastograms were created using mechanically induced transducer motion. Despite the loss in performance incurred when stress is applied using hand-induced transducer motion, it was concluded that the prototype system performed sufficiently well to warrant clinical evaluation.

Advancing intravascular ultrasonic palpation toward clinical applications.

This paper describes the first reported attempt to develop a real-time intravascular ultrasonic palpation system. We also report on our first experience in the catherization laboratory with this new elastographic imaging technique. The prototype system was based on commercially available intravascular ultrasound (US) scanner that was equipped with a 20-MHz array catheter. Digital beam-formed radiofrequency (RF) echo data (i.e., 12 bits, 100 Hz) was captured at full frame rate from the scanner and transferred to personal computer (PC) memory using a fast data-acquisition system. Composite palpograms were created by applying a one-dimensional (1-D) echo tracking technique in combination with global motion compensation and multiframe averaging to several pairs of RF echo frames that were obtained in the diastolic phase of the cardiac cycle. The quality of palpograms was assessed by conducting experiments on vessel phantoms and on patients. The results demonstrated that robust and consistent palpograms could be generated in almost real-time using the proposed system. Good correlation was observed between low strain values and regions of calcification as identified from the intravascular US (IVUS) sonograms. Although the clinical results are clearly preliminary, it was concluded that the prototype system performed sufficiently well to warrant further and more in-depth clinical investigation.

Evaluation of an iterative reconstruction method for quantitative elastography.

This paper describes an inverse reconstruction technique based on a modified Newton Raphson iterative scheme and the finite element method, which has been developed for computing the spatial distribution of Young's modulus from within soft tissues. Computer simulations were conducted to determine the relative merits of reconstructing tissue elasticity using knowledge of (a) known displacement boundary conditions (DBC), and (b) known stress boundary conditions (SBC). The results demonstrated that computing Young's modulus using knowledge of SBC allows accurate quantification of Young's modulus. However, the quality of the images produced using this reconstruction approach was dependent on the Young's modulus distribution assumed at the start of the reconstruction procedure. Computing Young's modulus from known DBC provided relative estimates of tissue elasticity which, despite the disadvantage of not being able to accurately quantify Young's modulus, formed images that were generally superior in quality to those produced using the known SBC, and were not affected by the trial solution. The results of preliminary experiments on phantoms demonstrated that this reconstruction technique is capable in practice of improving the fidelity of tissue elasticity images, reducing the artefacts otherwise present in strain images, and recovering Young's modulus images that possess excellent spatial and contrast resolution.