Precise two-dimensional D-bar reconstructions of human chest and phantom tank via sinc-convolution algorithm

BioMedical Engineering OnLine

Mahdi Abbasi¹  Ahmad-Reza Naghsh-Nilchi¹ 

[1] Department of Computer Engineering, Engineering Faculty, University of Isfahan, Isfahan, Iran

关键词: Human chest;    Chest phantom;    Accuracy measures;    Sinc-convolution;    D-bar;    EIT;   

Others  :  798072
DOI  :  10.1186/1475-925X-11-34

 received in 2012-01-06, accepted in 2012-05-07,  发布年份 2012

【 摘 要 】

Background

Electrical Impedance Tomography (EIT) is used as a fast clinical imaging technique for monitoring the health of the human organs such as lungs, heart, brain and breast. Each practical EIT reconstruction algorithm should be efficient enough in terms of convergence rate, and accuracy. The main objective of this study is to investigate the feasibility of precise empirical conductivity imaging using a sinc-convolution algorithm in D-bar framework.

Methods

At the first step, synthetic and experimental data were used to compute an intermediate object named scattering transform. Next, this object was used in a two-dimensional integral equation which was precisely and rapidly solved via sinc-convolution algorithm to find the square root of the conductivity for each pixel of image. For the purpose of comparison, multigrid and NOSER algorithms were implemented under a similar setting. Quality of reconstructions of synthetic models was tested against GREIT approved quality measures. To validate the simulation results, reconstructions of a phantom chest and a human lung were used.

Results

Evaluation of synthetic reconstructions shows that the quality of sinc-convolution reconstructions is considerably better than that of each of its competitors in terms of amplitude response, position error, ringing, resolution and shape-deformation. In addition, the results confirm near-exponential and linear convergence rates for sinc-convolution and multigrid, respectively. Moreover, the least degree of relative errors and the most degree of truth were found in sinc-convolution reconstructions from experimental phantom data. Reconstructions of clinical lung data show that the related physiological effect is well recovered by sinc-convolution algorithm.

Conclusions

Parametric evaluation demonstrates the efficiency of sinc-convolution to reconstruct accurate conductivity images from experimental data. Excellent results in phantom and clinical reconstructions using sinc-convolution support parametric assessment results and suggest the sinc-convolution to be used for precise clinical EIT applications.

【 授权许可】

   
2012 Abbasi and Naghsh-Nilchi; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

Figure 7.

【 参考文献 】

• [1]Webster JG: Electrical impedance tomography. 1st edition. New York: Adam Hilger; 1990.
• [2]Adler A, Arnold JH, Bayford R, Borsic A, Brown B, Dixon P, Faes TJC, Frerichs I, Gagnon H, Gärber Y, et al.: GREIT: a unified approach to 2D linear EIT reconstruction of lung images. Physiol Meas 2009, 30(6):S35.
• [3]Heinrich S, Schiffmann H, Frerichs A, Klockgether-Radke A, Frerichs I: Body and head position effects on regional lung ventilation in infants: an electrical impedance tomography study. Intensive Care Med 2006, 32(9):1392-1398.
• [4]Flores-Tapia D, O'Halloran M, Pistorius S: A bimodal reconstruction method for breast cancer imaging. Prog Electromagn Res 2011, 118:461-486.
• [5]Elke G, Pulletz S, Schädler D, Zick G, Gawelczyk B, Frerichs I, Weiler N: Measurement of regional pulmonary oxygen uptake—a novel approach using electrical impedance tomography. Physiol Meas 2011, 32(7):877.
• [6]Calderon AP: On inverse boundary value problem. Math App Comput 2006, 25(2):133-138.
• [7]Kotre CJ: EIT image reconstruction using sensitivity weighted filtered backprojection. Physiol Meas 1994, 15(2A):A125.
• [8]Wang H, Xu G, Zhang S, Yin N, Yan W: Implementation of generalized back projection algorithm in 3-D EIT. IEEE Trans Magn 2011, 47(5):1466-1469.
• [9]Cheney M, Isaacson D, Newell JC, Goble J, Simske S: Noser: an algorithm for solving the inverse conductivity problem. Int J Imaging Syst Technol 1990, 2:66-75.
• [10]Lionheart WRB: EIT reconstruction algorithms: Pitfalls, challenges and recent developments. Physiol Meas 2004, 25(1):125-142.
• [11]Hua P, Woo EJ, Webster JG, Tompkins WJ: Iterative reconstruction methods using regularization and optimal current patterns in electrical impedance tomography. IEEE Trans Med Imaging 1991, 10(4):621-628.
• [12]Kallman JS, Berryman JG: Weighted least-squares criteria for electrical impedance tomography. IEEE Trans Med Imaging 1992, 11(2):284-292.
• [13]Kaipio JP, Kolehmainen V, Somersalo E, Vauhkonen M: Statistical inversion and Monte Carlo sampling methods in electrical impedance tomography. Inverse Problems 2000, 16(5):1487-1522.
• [14]Rao L, He R, Wang Y, Yan W, Bai J, Ye D: An efficient improvement of modified Newton–Raphson algorithm for electrical impedance tomography. IEEE Trans Magn 1999, 35(3 PART 1):1562-1565.
• [15]Somersalo E, Cheney M, Isaacson D, Isaacson E: Layer stripping: A direct numerical method for impedance imaging. Inverse Problems 1991, 7(6):899-926.
• [16]Nachman AI: Global uniqueness for a two-dimensional inverse boundary value problem. Ann Math 1996, 143(1):71-96.
• [17]Abbasi A, Vosoughi-Vahdat B, Ebrahimi-Fakhim G: An inverse solution for 2D electrical impedance tomography based on electrical properties of materials. J App Sci 2009, 9(10):1962-1967.
• [18]Abbasi A, Vosoghi Vahdat B: Improving forward solution for 2D block electrical impedance tomography using modified equations. Sci Res Essays 2010, 5(11):1260-1263.
• [19]Abbasi A, Vosoghi Vahdat B: A non-itertive linear inverse solution for block approach in EIT. Comput Sci 2010, 1:190-196.
• [20]Nachman AI, Ablowitz MJ: Multidimensional inverse-scattering method. Stud Appl Math 1984, 71(3):243-250.
• [21]Isaacson D, Mueller JL, Newell JC, Siltanen S: Reconstructions of chest phantoms by the D-bar method for electrical impedance tomography. IEEE Trans Med Imaging 2004, 23(7):821-828.
• [22]Knudsen K, Mueller J, Siltanen S: Numerical solution method for the dbar-equation in the plane. J Comput Phys 2004, 198(2):500-517.
• [23]Siltanen S, Mueller J, Isaacson D: An implementation of the reconstruction algorithm of A Nachman for the 2D inverse conductivity problem. Inverse Problems 2000, 16(3):681-699.
• [24]Mueller JL, Siltanen S, Isaacson D: A direct reconstruction algorithm for electrical impedance tomography. IEEE Trans Med Imaging 2002, 21(6):555-559.
• [25]Vainikko G: Fast solvers of the Lippmann-Schwinger equation. Direct and Inverse problems of Mathical Physics 2000, 5:423-440.
• [26]DeAngelo M, Mueller JL: 2D D-bar reconstruction of the human chest and tank data using an improved approximation to the scattering transform. Physiol Meas 2010, 31:221-232.
• [27]Abbasi M, Naghsh-Nilchi AR: Iterative sinc-convolution method for solving planar D-bar equation with application to EIT. Biomed Eng: Int J Numer Meth 2012, 28(8):838-860.
• [28]Brown RM, Uhlman G: Uniqueness in the inverse conductivity problem for non-smooth conductivities in two dimensions. Commun Partial Differ Eq 1997, 22(5):1009-1027.
• [29]Astala K, Mueller JL, Päivärinta L, Siltanen S: Numerical computation of complex geometrical optics solutions to the conductivity equation. Appl Comput Harmon Anal 2010, 29(1):2-17.
• [30]Beals R, Coifman RR: The D-bar approach to inverse scattering and nonlinear evolutions. Physica D: Nonlinear Phenomena 1986, 18(1–3):242-249.
• [31]Faddeev LD: Increasing solutions of the schrodinger equation. Sov Phys Dokl 1965, 10:1033-1035.
• [32]Knudsen K, Lassas M, Mueller JL, Siltanen S: Regularized d-bar method for the inverse conductivity problem. Inverse Problems and Imaging 2009, 3(4):599-624.
• [33]Mueller JL, Siltanen S: Direct reconstructions of conductivities from boundary measurements. SIAM J Sci Comput 2003, 24(4):1232-1263.
• [34]Isaacson D, Mueller JL, Newell JC, Siltanen S: Imaging cardiac activity by the D-bar methd for electrical impedance tomography. Physiol Meas 2006, 27:S43-S50.
• [35]Du K, Nevanlinna O: A note on ℝ-linear GMRES for solving a class of ℝ-linear systems. Numer Linear Algebra Appl 2011. BioMed Central Full Text
• [36]Murphy EK, Mueller JL, Newell JC: Reconstructions of conductive and insulating targets using the D-bar method on an elliptical domain. Physiol Meas 2007, 28(7):S101-S114.
• [37]Murphy EK, Mueller JL: Effect of domain shape modeling and measurement errors on the 2-D D-bar method for EIT. IEEE Trans Med Imaging 2009, 28(10):1576-1584.
• [38]Cheng KS, Isaacson D, Newell JC, Gisser DG: Electrode models for electric current computed tomography. IEEE Trans Biomed Eng 1989, 36(9):918-924.
• [39]Paulson K, Breckon W, Pidcock M: Electrode modelling in electrical impedance tomography. SIAM J Appl Math 1992, 52(4):1012-1022.
• [40]Schroberl J: NETGEN: an advancing front 2D/3D-mesh generator based on abstract rules. Comput Vis Sci 1997, 1:41-52.
• [41]EIDORS: Calibrated chest shaped targets in 2D circular tank with ACT 3 EIT system. [http://eidors3d.sourceforge.net/data_contrib/ds_tank_phantom/ds_tank_phantom.shtml webcite]
• [42]Edic PM, Saulnier GJ, Newell JC, Isaacson D: A real-time electrical impedance tomograph. IEEE Trans Biomed Eng 1995, 42(9):849-859.
• [43]Isaacson D: Distinguishability of conductivities by electric current computed tomography. IEEE Trans Med Imaging 1986, MI-5(2):91-95.
• [44]EIDORS: EIT images of spontaneously breathing neonate. [http://eidors3d.sourceforge.net/data_contrib/if-neonate-spontaneous/index.shtml webcite]
• [45]Stenger F: Numerical methods based on sinc and analytic functions. 1st edition. New York: Springer; 1993.