On the convergence of the harmonic BZ algorithm in magnetic resonance electrical impedance tomography

J. J. Liu, J. K. Seo, M. Sinl, E. J. Woo

Research output: Contribution to journalArticlepeer-review

8 Citations (Scopus)


Magnetic resonance electrical impedance tomography (MREIT) is a new medical imaging technique that aims to provide electrical conductivity images with sufficiently high spatial resolution and accuracy. A new MREIT image reconstruction method called the harmonic Bz algorithm was proposed in 2002, and it is based on the measurement of Bz that is a single component of an induced magnetic flux density B = (Bx, By, Bz) subject to an injection current. Since then, MREIT imaging techniques have made significant progress, and recent published numerical simulations and phantom experiments show that we can produce high-quality conductivity images when the conductivity contrast is not very high. Though numerical simulations can explain why we could successfully distinguish different tissues with small conductivity differences, a rigorous mathematical analysis is required to better understand the underlying physical and mathematical principle. The purpose of this paper is to provide such a mathematical analysis of those numerical simulations and experimental results. By using a uniform a priori estimate for the solution of the elliptic equation in the divergent form and an induction argument, we show that, for a relatively small contrast of the target conductivity, the iterative harmonic Bz algorithm with a good initial guess is stable and exponentially convergent in the continuous norm. Both two- and three-dimensional versions of the algorithm are considered, and the difference in the convergence property of these two cases is analyzed. Some numerical results are also given to show the expected exponential convergence behavior.

Original languageEnglish
Pages (from-to)1259-1282
Number of pages24
JournalSIAM Journal on Applied Mathematics
Issue number5
Publication statusPublished - 2007

All Science Journal Classification (ASJC) codes

  • Applied Mathematics


Dive into the research topics of 'On the convergence of the harmonic BZ algorithm in magnetic resonance electrical impedance tomography'. Together they form a unique fingerprint.

Cite this