School of Electrical and Electronic Engineering, University of Manchester, Manchester, M60 1QD, UK.
IEEE Trans Biomed Eng. 2010 Apr;57(4):914-21. doi: 10.1109/TBME.2009.2036733. Epub 2009 Nov 20.
Magnetic induction tomography (MIT) is a noninvasive modality for imaging the complex conductivity (kappa = sigma + jomegaepsilon) or the magnetic permeability (mu) of a target under investigation. Because MIT employs noncontact coils for excitation and detection, MIT may be suitable for imaging biological tissues. In medical applications where high resolutions are sought, image reconstruction is a time and memory consuming task because the associated inverse problem is nonlinear and ill-posed. The time and memory constraints are mainly imposed by the solution of the forward problem within the iterative image reconstruction procedure. This paper investigates the application of a weakly coupled approximation to the solution of the forward problem and examines the accuracy against the computation time and memory gained in adopting this approximation. Initially, an analytical solution for mutual impedance change of a coil pair due to a large planar conductive object is presented based on a full wave theory and used to demonstrate a 10 MHz frequency excitation as an acceptable upper frequency limit under which the approximation is valid. Subsequently, a numerical impedance method adopting the approximation is presented. Here the impedance method is used to solve the forward problem, which employs electrical circuit analogues to mesh the target into a network that can be solved using circuit analysis and sparse matrix technique. The error due to the approximation is further estimated numerically with the impedance method against a commercial finite-element package (commercial FE solver, COMSOL) and results show at 10 MHz excitation a 0.4% of tolerance is achieved for conductivities in the range <0.5 S/m. The results also show the method can be applied for low conductivity medical applications and is computationally efficient compared to equivalent finite-element methods.
磁感应断层成像(MIT)是一种用于对目标的复杂电导率(kappa=sigma+jomegaepsilon)或磁导率(mu)进行成像的非侵入式方法。由于 MIT 采用非接触式线圈进行激励和检测,因此它可能适用于生物组织成像。在需要高分辨率的医学应用中,图像重建是一项耗时且耗内存的任务,因为相关的逆问题是非线性和不适定的。时间和内存的限制主要是由迭代图像重建过程中的正问题求解造成的。本文研究了将弱耦合近似应用于正问题求解,并考察了在采用该近似时计算时间和内存方面获得的准确性。首先,基于全波理论,提出了一种用于大平面导电物体对线圈对互阻抗变化的解析解,并将其用于演示 10 MHz 频率激励作为该近似有效的可接受的上限频率。随后,提出了一种采用该近似的数值阻抗法。这里,阻抗法用于求解正问题,该方法采用电路模拟将目标网格化为一个网络,然后可以使用电路分析和稀疏矩阵技术来求解该网络。进一步通过阻抗法对商业有限元包(商业 FE 求解器,COMSOL)进行数值估计,得到了近似误差,并在 10 MHz 激励下,在 0.5 S/m 范围内的电导率达到 0.4%的容差。结果还表明,该方法可应用于低电导率的医学应用,并且与等效有限元方法相比,计算效率更高。