Suitability of poroelastic and viscoelastic mechanical models for high and low frequency MR elastography.

PURPOSE: Descriptions of the structure of brain tissue as a porous cellular matrix support application of a poroelastic (PE) mechanical model which includes both solid and fluid phases. However, the majority of brain magnetic resonance elastography (MRE) studies use a single phase viscoelastic (VE) model to describe brain tissue behavior, in part due to availability of relatively simple direct inversion strategies for mechanical property estimation. A notable exception is low frequency intrinsic actuation MRE, where PE mechanical properties are imaged with a nonlinear inversion algorithm. METHODS: This paper investigates the effect of model choice at each end of the spectrum of in vivo human brain actuation frequencies. Repeat MRE examinations of the brains of healthy volunteers were used to compare image quality and repeatability for each inversion model for both 50 Hz externally produced motion and ≈1 Hz intrinsic motions. Additionally, realistic simulated MRE data were generated with both VE and PE finite element solvers to investigate the effect of inappropriate model choice for ideal VE and PE materials. RESULTS: In vivo, MRE data revealed that VE inversions appear more representative of anatomical structure and quantitatively repeatable for 50 Hz induced motions, whereas PE inversion produces better results at 1 Hz. Reasonable VE approximations of PE materials can be derived by equating the equivalent wave velocities for the two models, provided that the timescale of fluid equilibration is not similar to the period of actuation. An approximation of the equilibration time for human brain reveals that this condition is violated at 1 Hz but not at 50 Hz. Additionally, simulation experiments when using the "wrong" model for the inversion demonstrated reasonable shear modulus reconstructions at 50 Hz, whereas cross-model inversions at 1 Hz were poor quality. Attenuation parameters were sensitive to changes in the forward model at both frequencies, however, no spatial information was recovered because the mechanisms of VE and PE attenuation are different. CONCLUSIONS: VE inversions are simpler with fewer unknown properties and may be sufficient to capture the mechanical behavior of PE brain tissue at higher actuation frequencies. However, accurate modeling of the fluid phase is required to produce useful mechanical property images at the lower frequencies of intrinsic brain motions.

Multiresolution MR elastography using nonlinear inversion.

PURPOSE: Nonlinear inversion (NLI) in MR elastography requires discretization of the displacement field for a finite element (FE) solution of the "forward problem", and discretization of the unknown mechanical property field for the iterative solution of the "inverse problem". The resolution requirements for these two discretizations are different: the forward problem requires sufficient resolution of the displacement FE mesh to ensure convergence, whereas lowering the mechanical property resolution in the inverse problem stabilizes the mechanical property estimates in the presence of measurement noise. Previous NLI implementations use the same FE mesh to support the displacement and property fields, requiring a trade-off between the competing resolution requirements. METHODS: This work implements and evaluates multiresolution FE meshes for NLI elastography, allowing independent discretizations of the displacements and each mechanical property parameter to be estimated. The displacement resolution can then be selected to ensure mesh convergence, and the resolution of the property meshes can be independently manipulated to control the stability of the inversion. RESULTS: Phantom experiments indicate that eight nodes per wavelength (NPW) are sufficient for accurate mechanical property recovery, whereas mechanical property estimation from 50 Hz in vivo brain data stabilizes once the displacement resolution reaches 1.7 mm (approximately 19 NPW). Viscoelastic mechanical property estimates of in vivo brain tissue show that subsampling the loss modulus while holding the storage modulus resolution constant does not substantially alter the storage modulus images. Controlling the ratio of the number of measurements to unknown mechanical properties by subsampling the mechanical property distributions (relative to the data resolution) improves the repeatability of the property estimates, at a cost of modestly decreased spatial resolution. CONCLUSIONS: Multiresolution NLI elastography provides a more flexible framework for mechanical property estimation compared to previous single mesh implementations.

Retrospective image correction in the presence of nonlinear temporal magnetic field changes using multichannel navigator echoes.

Spatio-temporal magnetic field changes in the brain caused by breathing or body movements can lead to image artifacts. This is especially a problem in T(2)(*)-weighted sequences. With the acquisition of an extra echo (navigator), it is possible to measure the magnetic field change induced frequency offset for a given slice during image acquisition. However, substantial local variation across a slice can occur. This work describes an extension of the conventional navigator technique that improves the estimation of the magnetic field distribution in the brain during strong field fluctuations. This is done using the combination of signals from multiple coil elements, the coil sensitivity profiles, and frequency encoding: termed sensitivity-encoded navigator echoes. In vivo validation was performed in subjects who performed normal breathing, nose touching, and deep breathing during scanning. The sensitivity-encoded navigator technique leads to an error reduction in estimating the field distribution in the brain of 73% ± 16% compared with 56% ± 14% for conventional estimation. Image quality can be improved via incorporating this navigator information appropriately into the image reconstruction. When the sensitivity-encoded navigator technique was applied to a T(2)(*)-weighted sequence at 7 T, a ghosting reduction of 47% ± 13% was measured during nose touching experiments compared with no correction.

Joint estimation of water/fat images and field inhomogeneity map.

Water/fat separation in the presence of B 0 field inhomogeneity is a problem of considerable practical importance in MRI. This article describes two complementary methods for estimating the water/fat images and the field inhomogeneity map from Dixon-type acquisitions. One is based on variable projection (VARPRO) and the other on linear prediction (LP). The VARPRO method is very robust and can be used in low signal-to-noise ratio conditions because of its ability to achieve the maximum-likelihood solution. The LP method is computationally more efficient, and is shown to perform well under moderate levels of noise and field inhomogeneity. These methods have been extended to handle multicoil acquisitions by jointly solving the estimation problem for all the coils. Both methods are analyzed and compared and results from several experiments are included to demonstrate their performance.

Accelerating Advanced MRI Reconstructions on GPUs.

Computational acceleration on graphics processing units (GPUs) can make advanced magnetic resonance imaging (MRI) reconstruction algorithms attractive in clinical settings, thereby improving the quality of MR images across a broad spectrum of applications. This paper describes the acceleration of such an algorithm on NVIDIA's Quadro FX 5600. The reconstruction of a 3D image with 128(3) voxels achieves up to 180 GFLOPS and requires just over one minute on the Quadro, while reconstruction on a quad-core CPU is twenty-one times slower. Furthermore, relative to the true image, the error exhibited by the advanced reconstruction is only 12%, while conventional reconstruction techniques incur error of 42%.