  Magnetic resonance electrical property tomography is a recent medical imaging
modality for visualizing the electrical tissue properties of the human body
using radio-frequency magnetic fields. It uses the fact that in magnetic
resonance imaging systems the eddy currents induced by the radio-frequency
magnetic fields reflect the conductivity ($\sigma$) and permittivity
($\epsilon$) distributions inside the tissues through Maxwell's equations. The
corresponding inverse problem consists of reconstructing the admittivity
distribution ($\gamma=\sigma+i\omega\epsilon$) at the Larmor frequency
($\omega/2\pi=$128 MHz for a 3 tesla MRI machine) from the positive circularly
polarized component of the magnetic field ${\bf H}=(H_x,H_y,H_z)$. Previous
methods are usually based on an assumption of local homogeneity
($\nabla\gamma\approx 0$) which simplifies the governing equation. However,
previous methods that include the assumption of homogeneity are prone to
artifacts in the region where $\gamma$ varies. Hence, recent work has sought a
reconstruction method that does not assume local-homogeneity. This paper
presents a new magnetic resonance electrical property tomography reconstruction
method which does not require any local homogeneity assumption on $\gamma$. We
find that $\gamma$ is a solution of a semi-elliptic partial differential
equation with its coefficients depending only on the measured data $H^+$, which
enable us to compute a blurred version of $\gamma$. To improve the resolution
of the reconstructed image, we developed a new optimization algorithm that
minimizes the mismatch between the data and the model data as a highly
nonlinear function of $\gamma$. Numerical simulations are presented to
illustrate the potential of the proposed reconstruction method.
