摘要

Magnetic resonance electrical impedance tomography (MREIT) takes advantage of internal information to solve its nonlinear inverse problem of recovering a conductivity distribution inside an imaging object. When we inject current into the imaging object, there occurs a distribution of internal magnetic flux density B = (B-x, B-y, B-z). In MREIT we utilize a magnetic resonance imaging scanner with its main magnetic field in the z direction to acquire B-z data. The harmonic B-z algorithm was invented in 2001 to reconstruct cross-sectional conductivity images from B-z data sets subject to multiple injection currents. Utilizing internal B-z data, it overcomes the inherent ill posedness in electrical impedance tomography. We can set up the inverse problem in MREIT as a coefficient identification problem of finding sigma appearing in del . (sigma del u) = 0 from acquired data of the z component of del x (sigma del u). The harmonic Bz algorithm has shown an excellent performance in numerical simulations and phantom experiments. Experimental MREIT studies have now reached the stage of in vivo animal and human imaging experiments. However, there is not much work on rigorous mathematical theories of error estimate and convergence analysis yet. The purpose of this paper is to provide a posteriori error estimate in MREIT conductivity image reconstructions. This enables us to evaluate a difference between a reconstructed conductivity image and the unknown true conductivity image. We also describe a convergence analysis of the harmonic B-z algorithm, which improves the previous result of [J. J. Liu, J. K. Seo, M. Sini, and E. J. Woo, SIAM J. Appl. Math., 67 (2007), pp. 1259-1282] in the sense that assumptions on the conductivity are much relaxed.