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

J. J. Liu, Jin Keun Seo, M. Sinl, E. J. Woo

Research output: Contribution to journalArticle

26 Citations (Scopus)

Abstract

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
Volume67
Issue number5
DOIs
Publication statusPublished - 2007 Oct 26

Fingerprint

Electrical Impedance Tomography
Acoustic impedance
Magnetic Resonance
Magnetic resonance
Conductivity
Tomography
Harmonic
Computer simulation
Mathematical Analysis
Numerical Simulation
Imaging techniques
Medical imaging
Magnetic flux
Image reconstruction
Uniform Estimates
Exponential Convergence
Electrical Conductivity
Medical Imaging
Guess
Image Reconstruction

All Science Journal Classification (ASJC) codes

  • Applied Mathematics

Cite this

@article{3b0074bd975b4954ae64caad404de78d,
title = "On the convergence of the harmonic BZ algorithm in magnetic resonance electrical impedance tomography",
abstract = "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.",
author = "Liu, {J. J.} and Seo, {Jin Keun} and M. Sinl and Woo, {E. J.}",
year = "2007",
month = "10",
day = "26",
doi = "10.1137/060661892",
language = "English",
volume = "67",
pages = "1259--1282",
journal = "SIAM Journal on Applied Mathematics",
issn = "0036-1399",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "5",

}

On the convergence of the harmonic BZ algorithm in magnetic resonance electrical impedance tomography. / Liu, J. J.; Seo, Jin Keun; Sinl, M.; Woo, E. J.

In: SIAM Journal on Applied Mathematics, Vol. 67, No. 5, 26.10.2007, p. 1259-1282.

Research output: Contribution to journalArticle

TY - JOUR

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

AU - Liu, J. J.

AU - Seo, Jin Keun

AU - Sinl, M.

AU - Woo, E. J.

PY - 2007/10/26

Y1 - 2007/10/26

N2 - 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.

AB - 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.

UR - http://www.scopus.com/inward/record.url?scp=35348986496&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=35348986496&partnerID=8YFLogxK

U2 - 10.1137/060661892

DO - 10.1137/060661892

M3 - Article

VL - 67

SP - 1259

EP - 1282

JO - SIAM Journal on Applied Mathematics

JF - SIAM Journal on Applied Mathematics

SN - 0036-1399

IS - 5

ER -