Projects:QuantitativeSusceptibilityMapping

There is increasing evidence that excessive iron deposition in specific regions of the brain is associated with neurodegenerative disorders such as Alzheimer's and Parkinson's disease [1]. The role of iron in the pathogenesis of these diseases remains unknown and is difficult to determine without a non-invasive method to quantify its concentration in-vivo. Since iron is a ferromagnetic substance, changes in iron concentration result in local changes in the magnetic susceptibility of tissue. In magnetic resonance imaging (MRI) experiments, differences in magnetic susceptibility cause perturbations in the local magnetic field, which can be computed from the phase of the MR signal (in a gradient echo sequence, the observed field is proportional to the MR phase).

Description

In MRI, magnetic susceptibility differences cause measurable perturbations in the local magnetic field that can be modeled as the convolution of a dipole-like kernel with the spatial susceptibility distribution. In the Fourier domain, the kernel exhibits zeros at the magic angle, preventing direct inversion of the field map; also, limited observations make the problem ill-posed. The observed data is also corrupted by confounding fields (ie. those from tissue/air interfaces, mis-set shims, and other non-local sources). Previous work has shown that MR images can be successfully reconstructed from under-sampled observations by exploiting the sparsity of in-vivo data under various transforms using methods from compressed sensing [1]. In susceptibility estimation, the forward model results in under-sampling of the data in the Fourier domain, but accurate estimates can be obtained using the Laplacian and L1 norm, which promote sparse solutions while removing external field artifacts.

Our variational method for susceptibility estimation is described in Fig. 1.

Forward Model (Eq. 1)

The forward model relates the perturbing field to the unknown susceptibility through a local term and convolution of the second z-derivative of the Green’s function of the Laplacian [2].

Bias Field Elimination (Eq. 2)

Applying the Laplacian removes non-local phase effects such as shim fields, which are a solution to the Laplace equation.

Objective Function (Eq. 3)

The first term provides regularization, penalizing solutions with large differences in spatial frequency structure relative to the magnitude image. The second term penalizes departures from Eq. 2, by enforcing agreement of high frequency phase effects while eliminating low order bias fields.

Data Acquisition

Cylindrical and rectangular phantoms were made using Magnevist (gadopentetate dimeglumine) solutions of 0.5, 1.0, 2.0, and 3.0 mM corresponding to susceptibility values of 0.15, 0.31, 0.62, and 0.94 ppm [3,4]. Field maps were obtained using a 3D multi-echo GRE sequence on a 3T Siemens Trio MRI.

Results

Application of the Laplacian removes substantial inhomogeniety in the field map in both phantoms as shown in Fig. 2. Rectangular phantom: mean estimated susceptibility values for water and Magnevist were -9.049 and 0.6273 ppm, with true values of -9.050 and 0.6270 ppm. Cylindrical phantom: the estimated susceptibility map allowed different concentrations of Magnevist to be clearly identified and reasonable estimates were obtained in the presence of significant noise and bias due to external field effects.