Computation of Magnetic Field Distortions and Impact on T₂∗-weighted MRI, with Applications to Magnetic Susceptibility Parameter Estimation

A two-step numerical computation of T₂∗ signal weighting maps in gradient echo magnetic resonance imaging in the presence of an object with varied susceptibility property is presented. In the first step, the magnetic scalar potential is computed for an arbitrary 2D magnetic susceptibility distribution using an algebraic solver. The corresponding magnetic field disturbance is computed from the magnetic scalar potential. In the second step, nonlinear operations are used to compute T₂∗ from the magnetic field disturbance and then to generate a map of T₂∗ signal weighting. The linearity of the first step of the solution process is used to implement a superposition of basis solutions approach that increases computational efficiency. Superposition of basis solutions, computed from a system composed of a single node of differing magnetic susceptibility from the surround, herein referred to as the base system, is found to provide an accurate estimation of the scalar potential for arbitrary susceptibility distributions. Afterwards, nonlinear computation of the T₂∗ signal weighting maps can be performed. The properties of the algebraic magnetic scalar potential solver are discussed in this work. Finally, the linearity of the magnetic scalar potential solver is used to estimate the magnetic susceptibility of various objects from in vitro MR-imaging data acquired at 9.4 T.