An Inverse Problems Approach to MR-EPT Image Reconstruction.

Magnetic Resonance-Electrical Properties Tomography (MR-EPT) is an imaging modality that maps the spatial distribution of the electrical conductivity and permittivity using standard MRI systems. The presence of a body within the scanner alters the RF field, and by mapping these alterations it is possible to recover the electrical properties. The field is time-harmonic, and can be described by the Helmholtz equation. Approximations to this equation have been previously used to estimate conductivity and permittivity in terms of first or second derivatives of RF field data. Using these same approximations, an inverse approach to solving the MR-EPT problem is presented here that leverages a forward model for describing the magnitude and phase of the field within the imaging domain, and a fitting approach for estimating the electrical properties distribution. The advantages of this approach are that 1) differentiation of the measured data is not required, thus reducing noise sensitivity, and 2) different regularization schemes can be adopted, depending on prior knowledge of the distribution of conductivity or permittivity, leading to improved image quality. To demonstrate the developed approach, both Quadratic (QR) and Total Variation (TV) regularization methods were implemented and evaluated through numerical simulation and experimentally acquired data. The proposed inverse approach to MR-EPT reconstruction correctly identifies contrasts and accurately reconstructs the geometry in both simulations and experiments. The TV regularized scheme reconstructs sharp spatial transitions, which are difficult to reconstruct with other, more traditional approaches.

Temperature of the magnetic nanoparticle microenvironment: estimation from relaxation times.

Accurate temperature measurements are essential to safe and effective thermal therapies for cancer and other diseases. However, conventional thermometry is challenging so using the heating agents themselves as probes allows for ideal local measurements. Here, we present a new noninvasive method for measuring the temperature of the microenvironment surrounding magnetic nanoparticles from the Brownian relaxation time of nanoparticles. Experimentally, the relaxation time can be determined from the nanoparticle magnetization induced by an alternating magnetic field at various applied frequencies. A previously described method for nanoparticle temperature estimation used a low frequency Langevin function description of magnetic dipoles and varied the excitation field amplitude to estimate the energy state distribution and the corresponding temperature. We show that the new method is more accurate than the previous method at higher applied field frequencies that push the system farther from equilibrium.

Subzone based magnetic resonance elastography using a Rayleigh damped material model.

PURPOSE: Recently, the attenuating behavior of soft tissue has been addressed in magnetic resonance elastography by the inclusion of a damping mechanism in the methods used to reconstruct the resulting mechanical property image. To date, this mechanism has been based on a viscoelastic model for material behavior. Rayleigh, or proportional, damping provides a more generalized model for elastic energy attenuation that uses two parameters to characterize contributions proportional to elastic and inertial forces. In the case of time-harmonic vibration, these two parameters lead to both the elastic modulus and the density being complex valued (as opposed to the case of pure viscoelasticity, where only the elastic modulus is complex valued). METHODS: This article presents a description of Rayleigh damping in the time-harmonic case, discussing the differences between this model and the viscoelastic damping models. In addition, the results from a subzone based Rayleigh damped elastography study of gelatin and tofu phantoms are discussed, along with preliminary results from in vivo breast data. RESULTS: Both the phantom and the tissue studies presented here indicate a change in the Rayleigh damping structure, described as Rayleigh composition, between different material types, with tofu and healthy tissue showing lower Rayleigh composition values than gelatin or cancerous tissue. CONCLUSIONS: It is possible that Rayleigh damping elastography and the concomitant Rayleigh composition images provide a mechanism for differentiating tissue structure in addition to measuring elastic stiffness and attenuation. Such information could be valuable in the use of Rayleigh damped magnetic resonance elastography as a diagnostic imaging tool.

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.

Time-harmonic magnetic resonance elastography of the normal feline brain.

Imaging of the mechanical properties of in vivo brain tissue could eventually lead to non-invasive diagnosis of hydrocephalus, Alzheimer's disease and other pathologies known to alter the intracranial environment. The purpose of this work is to (1) use time-harmonic magnetic resonance elastography (MRE) to estimate the mechanical property distribution of cerebral tissue in the normal feline brain and (2) compare the recovered properties of grey and white matter. Various in vivo and ex vivo brain tissue property measurement strategies have led to the highly variable results that have been reported in the literature. MR elastography is an imaging technique that can estimate mechanical properties of tissue non-invasively and in vivo. Data was acquired in 14 felines and elastic parameters were estimated using a globo-regional nonlinear image reconstruction algorithm. Results fell within the range of values reported in the literature and showed a mean shear modulus across the subject group of 7-8 kPa with all but one animal falling within 5-15 kPa. White matter was statistically stiffer (p<0.01) than grey matter by about 1 kPa on a per subject basis. To the best of our knowledge, the results reported represent the most extensive set of estimates in the in vivo brain which have been based on MRE acquisition of the three-dimensional displacement field coupled to volumetric shear modulus image reconstruction achieved through nonlinear parameter estimation. However, the inter-subject variation in mean shear modulus indicates the need for further study, including the possibility of applying more advanced models to estimate the relevant tissue mechanical properties from the data.