1
|
Chen R, Krueger-Ziolek S, Lovas A, Benyó B, Rupitsch SJ, Moeller K. Structural priors represented by discrete cosine transform improve EIT functional imaging. PLoS One 2023; 18:e0285619. [PMID: 37167237 PMCID: PMC10174522 DOI: 10.1371/journal.pone.0285619] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Accepted: 04/26/2023] [Indexed: 05/13/2023] Open
Abstract
Structural prior information can improve electrical impedance tomography (EIT) reconstruction. In this contribution, we introduce a discrete cosine transformation-based (DCT-based) EIT reconstruction algorithm to demonstrate a way to incorporate the structural prior with the EIT reconstruction process. Structural prior information is obtained from other available imaging methods, e.g., thorax-CT. The DCT-based approach creates a functional EIT image of regional lung ventilation while preserving the introduced structural information. This leads to an easier interpretation in clinical settings while maintaining the advantages of EIT in terms of bedside monitoring during mechanical ventilation. Structural priors introduced in the DCT-based approach are of two categories in terms of different levels of information included: a contour prior only differentiates lung and non-lung region, while a detail prior includes information, such as atelectasis, within the lung area. To demonstrate the increased interpretability of the EIT image through structural prior in the DCT-based approach, the DCT-based reconstructions were compared with reconstructions from a widely applied one-step Gauss-Newton solver with background prior and from the advanced GREIT algorithm. The comparisons were conducted both on simulation data and retrospective patient data. In the simulation, we used two sets of forward models to simulate different lung conditions. A contour prior and a detail prior were derived from simulation ground truth. With these two structural priors, the reconstructions from the DCT-based approach were compared with the reconstructions from both the one-step Gauss-Newton solver and the GREIT. The difference between the reconstructions and the simulation ground truth is calculated by the ℓ2-norm image difference. In retrospective patient data analysis, datasets from six lung disease patients were included. For each patient, a detail prior was derived from the patient's CT, respectively. The detail prior was used for the reconstructions using the DCT-based approach, which was compared with the reconstructions from the GREIT. The reconstructions from the DCT-based approach are more comprehensive and interpretable in terms of preserving the structure specified by the priors, both in simulation and retrospective patient data analysis. In simulation analysis, the ℓ2-norm image difference of the DCT-based approach with a contour prior decreased on average by 34% from GREIT and 49% from the Gauss-Newton solver with background prior; for reconstructions of the DCT-based approach with detail prior, on average the ℓ2-norm image difference is 53% less than GREIT and 63% less than the reconstruction with background prior. In retrospective patient data analysis, the reconstructions from both the DCT-based approach and GREIT can indicate the current patient status, but the DCT-based approach yields more interpretable results. However, it is worth noting that the preserved structure in the DCT-based approach is derived from another imaging method, not from the EIT measurement. If the structural prior is outdated or wrong, the result might be misleadingly interpreted, which induces false clinical conclusions. Further research in terms of evaluating the validity of the structural prior and detecting the outdated prior is necessary.
Collapse
Affiliation(s)
- Rongqing Chen
- Institute of Technical Medicine (ITeM), Furtwangen University, Villingen-Schwenningen, Germany
- Faculty of Engineering, University of Freiburg, Freiburg, Germany
| | - Sabine Krueger-Ziolek
- Institute of Technical Medicine (ITeM), Furtwangen University, Villingen-Schwenningen, Germany
| | - András Lovas
- Department of Anaesthesiology and Intensive Therapy, Kiskunhalas Semmelweis Hospital, Kiskunhalas, Hungary
| | - Balázs Benyó
- Department of Control Engineering and Information Technology, Faculty of Electrical Engineering and Informatics, Budapest University of Technology and Economics, Budapest, Hungary
| | | | - Knut Moeller
- Institute of Technical Medicine (ITeM), Furtwangen University, Villingen-Schwenningen, Germany
| |
Collapse
|
2
|
Santos TBR, Nakanishi RM, de Camargo EDLB, Amato MBP, Kaipio JP, Lima RG, Mueller JL. Improved resolution of D-bar images of ventilation using a Schur complement property and an anatomical atlas. Med Phys 2022; 49:4653-4670. [PMID: 35411573 PMCID: PMC9544658 DOI: 10.1002/mp.15669] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 03/31/2022] [Accepted: 03/31/2022] [Indexed: 11/18/2022] Open
Abstract
BACKGROUND Electrical impedance tomography (EIT) is a nonionizing imaging technique for real-time imaging of ventilation of patients with respiratory distress. Cross-sectional dynamic images are formed by reconstructing the conductivity distribution from measured voltage data arising from applied alternating currents on electrodes placed circumferentially around the chest. Since the conductivity of lung tissue depends on air content, blood flow, and the presence of pathology, the dynamic images provide regional information about ventilation, pulsatile perfusion, and abnormalities. However, due to the ill-posedness of the inverse conductivity problem, EIT images have a coarse spatial resolution. One method of improving the resolution is to include prior information in the reconstruction. PURPOSE In this work, we propose a technique in which a statistical prior built from an anatomical atlas is used to postprocess EIT reconstructions of human chest data. The effectiveness of the method is demonstrated on data from two patients with cystic fibrosis. METHODS A direct reconstruction algorithm known as the D-bar method was used to compute a two-dimensional reconstruction of the conductivity distribution in the plane of the electrodes. Reconstructions using one step in an iterative (regularized) Newton's method were also computed for comparison. An anatomical atlas consisting of 1589 synthetic EIT images computed from X-ray computed tomography (CT) scans of 74 adult male subjects was computed for use as a statistical prior. The resolution of the D-bar images was then improved by maximizing the conditional probability density function of an image that is consistent with the a priori information and the statistical model. A new method to evaluate the accuracy of the EIT images using CT scans of the imaged patient as ground truth is presented. The novel approach is tested on data from two patients with cystic fibrosis. RESULTS AND CONCLUSIONS The D-bar images resulted in better structural similarity index measures (SSIM) and multiscale (MS) SSIM measures for both subjects using the mask or amplitude evaluation approach than the one-step (regularized) Newton's method. Further improvement was achieved using the Schur complement (SC) approach, with MS-SSIM values of 0.718 and 0.682 using SC evaluated with the mask and amplitude approach, respectively, for Patient 1, and MS-SSIM values of 0.726 and 0.692 using SC evaluated with the mask and amplitude approach, respectively, for Patient 2. The results from applying an anatomical atlas and statistical prior to EIT data from two patients with cystic fibrosis suggest that the spatial resolution of the EIT image can be improved to reveal pathology that may be difficult to discern in the original EIT image. The novel metric of evaluation is consistent with the appearance of improved spatial resolution and provides a new way to evaluate the accuracy of EIT reconstructions when a CT scan is available.
Collapse
Affiliation(s)
| | - Rafael Mikio Nakanishi
- Mechanical Engineering DepartmentPolytechnic School of the University of São PauloSão PauloSPBrazil
| | | | | | - Jari P. Kaipio
- Department of MathematicsUniversity of AucklandNew Zealand
- Department of Applied PhysicsUniversity of Eastern FinlandKuopioFinland
| | - Raul Gonzalez Lima
- Mechanical Engineering DepartmentPolytechnic School of the University of São PauloSão PauloSPBrazil
| | - Jennifer L. Mueller
- Department of Mathematics and School of Biomedical Engineering and the Department of Electrical and Computer EngineeringColorado State UniversityColoradoUSA
| |
Collapse
|
3
|
Mueller JL. Evaluation of Pulmonary Structure and Function in Patients with Cystic Fibrosis from Electrical Impedance Tomography Data. Methods Mol Biol 2022; 2393:733-750. [PMID: 34837209 DOI: 10.1007/978-1-0716-1803-5_39] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Electrical impedance tomography (EIT) is a medical imaging technique in which low frequency, low amplitude electromagnetic fields applied through electrodes on the skin are used to compute the conductivity and/or permittivity inside the body and form functional images from the reconstructed values. This work describes methods of computing EIT-derived surrogate measures of pulmonary function and identifying regions of air trapping and consolidation from functional EIT images. These methods were developed for pediatric patients with cystic fibrosis, for whom a real-time non-ionizing imaging modality can be of great benefit for monitoring disease progression or a pulmonary exacerbation.
Collapse
Affiliation(s)
- Jennifer L Mueller
- Department of Mathematics and School of Biomedical Engineering, Colorado State University, Fort Collins, CO, USA.
| |
Collapse
|
4
|
Herzberg W, Rowe DB, Hauptmann A, Hamilton SJ. Graph Convolutional Networks for Model-Based Learning in Nonlinear Inverse Problems. IEEE TRANSACTIONS ON COMPUTATIONAL IMAGING 2021; 7:1341-1353. [PMID: 35873096 PMCID: PMC9307146 DOI: 10.1109/tci.2021.3132190] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The majority of model-based learned image reconstruction methods in medical imaging have been limited to uniform domains, such as pixelated images. If the underlying model is solved on nonuniform meshes, arising from a finite element method typical for nonlinear inverse problems, interpolation and embeddings are needed. To overcome this, we present a flexible framework to extend model-based learning directly to nonuniform meshes, by interpreting the mesh as a graph and formulating our network architectures using graph convolutional neural networks. This gives rise to the proposed iterative Graph Convolutional Newton-type Method (GCNM), which includes the forward model in the solution of the inverse problem, while all updates are directly computed by the network on the problem specific mesh. We present results for Electrical Impedance Tomography, a severely ill-posed nonlinear inverse problem that is frequently solved via optimization-based methods, where the forward problem is solved by finite element methods. Results for absolute EIT imaging are compared to standard iterative methods as well as a graph residual network. We show that the GCNM has good generalizability to different domain shapes and meshes, out of distribution data as well as experimental data, from purely simulated training data and without transfer training.
Collapse
Affiliation(s)
- William Herzberg
- Department of Mathematical and Statistical Sciences; Marquette University, Milwaukee, WI 53233 USA
| | - Daniel B Rowe
- Department of Mathematical and Statistical Sciences; Marquette University, Milwaukee, WI 53233 USA
| | - Andreas Hauptmann
- Research Unit of Mathematical Sciences; University of Oulu, Oulu, Finland and with the Department of Computer Science; University College London, London, United Kingdom
| | - Sarah J Hamilton
- Department of Mathematical and Statistical Sciences; Marquette University, Milwaukee, WI 53233 USA
| |
Collapse
|
5
|
Alsaker M, Cárdenas DAC, Furuie SS, Mueller JL. Complementary use of priors for pulmonary imaging with electrical impedance and ultrasound computed tomography. JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS 2021; 395:113591. [PMID: 34092904 PMCID: PMC8177074 DOI: 10.1016/j.cam.2021.113591] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
For medical professionals caring for patients undergoing mechanical ventilation due to respiratory failure, the ability to quickly and safely obtain images of pulmonary function at the patient's bedside would be highly desirable. Such images could be used to provide early warnings of developing pulmonary pathologies in real time, thereby reducing the incidence of complications and improving patient outcomes. Electrical impedance tomography (EIT) and low-frequency ultrasound computed tomography (USCT) are two imaging techniques with the potential to provide real-time non-ionizing pulmonary monitoring in the ICU setting, and each method has its own unique advantages as well as drawbacks. In this work, we describe a new algorithm for a system in which the strengths of the two modalities are combined in a complementary fashion. Specifically, preliminary reconstructions from each modality are used as priors to stabilize subsequent reconstructions, providing improved spatial resolution, sharper organ boundaries, and enhanced appearance of pathologies and other features. Results are validated using three numerically simulated thoracic phantoms representing pulmonary pathologies.
Collapse
Affiliation(s)
- Melody Alsaker
- Department of Mathematics; Gonzaga University, Spokane, WA 99258 USA
| | | | | | - Jennifer L. Mueller
- Department of Mathematics and School of Biomedical Engineering, Colorado State University, CO 80523 USA
| |
Collapse
|
6
|
Hamilton SJ, Isaacson D, Kolehmainen V, Muller PA, Toivanen J, Bray PF. 3D ELECTRICAL IMPEDANCE TOMOGRAPHY RECONSTRUCTIONS FROM SIMULATED ELECTRODE DATA USING DIRECT INVERSION t exp AND CALDERÓN METHODS. INVERSE PROBLEMS AND IMAGING (SPRINGFIELD, MO.) 2021; 15:1135-1169. [PMID: 35173824 PMCID: PMC8846426 DOI: 10.3934/ipi.2021032] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The first numerical implementation of a t exp method in 3D using simulated electrode data is presented. Results are compared to Calderón's method as well as more common TV and smoothness regularization-based methods. The t exp method for EIT is based on tailor-made non-linear Fourier transforms involving the measured current and voltage data. Low-pass filtering in the non-linear Fourier domain is used to stabilize the reconstruction process. In 2D, t exp methods have shown great promise for providing robust real-time absolute and time-difference conductivity reconstructions but have yet to be used on practical electrode data in 3D, until now. Results are presented for simulated data for conductivity and permittivity with disjoint non-radially symmetric targets on spherical domains and noisy voltage data. The 3D t exp and Calderón methods are demonstrated to provide comparable quality to their 2D counterparts, and hold promise for real-time reconstructions due to their fast, non-optimized, computational cost.
Collapse
Affiliation(s)
- S J Hamilton
- Department of Mathematical and Statistical Sciences; Marquette University, Milwaukee, WI 53233 USA
| | - D Isaacson
- Department of Mathematical Sciences, Rensselaer Polytechnic Institute, Troy, NY, 12180, USA
| | - V Kolehmainen
- Department of Applied Physics, University of Eastern Finland, FI-70210 Kuopio, Finland
| | - P A Muller
- Department of Mathematics & Statistics; Villanova University, Villanova, PA 19085 USA
| | - J Toivanen
- Department of Applied Physics, University of Eastern Finland, FI-70210 Kuopio, Finland
| | - P F Bray
- Department of Mathematics; Drexel University, Philadelphia, PA 19104 USA
| |
Collapse
|
7
|
Abdulla UG, Bukshtynov V, Seif S. Cancer detection through Electrical Impedance Tomography and optimal control theory: theoretical and computational analysis. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2021; 18:4834-4859. [PMID: 34198468 DOI: 10.3934/mbe.2021246] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The Inverse Electrical Impedance Tomography (EIT) problem on recovering electrical conductivity tensor and potential in the body based on the measurement of the boundary voltages on the $ m $ electrodes for a given electrode current is analyzed. A PDE constrained optimal control framework in Besov space is developed, where the electrical conductivity tensor and boundary voltages are control parameters, and the cost functional is the norm difference of the boundary electrode current from the given current pattern and boundary electrode voltages from the measurements. The novelty of the control-theoretic model is its adaptation to the clinical situation when additional "voltage-to-current" measurements can increase the size of the input data from $ m $ up to $ m! $ while keeping the size of the unknown parameters fixed. The existence of the optimal control and Fréchet differentiability in the Besov space along with optimality condition is proved. Numerical analysis of the simulated model example in the 2D case demonstrates that by increasing the number of input boundary electrode currents from $ m $ to $ m^2 $ through additional "voltage-to-current" measurements the resolution of the electrical conductivity of the body identified via gradient method in Besov space framework is significantly improved.
Collapse
Affiliation(s)
- Ugur G Abdulla
- Department of Mathematical Sciences, Florida Institute of Technology, Melbourne, FL 32901, USA
| | - Vladislav Bukshtynov
- Department of Mathematical Sciences, Florida Institute of Technology, Melbourne, FL 32901, USA
| | - Saleheh Seif
- Department of Mathematical Sciences, Florida Institute of Technology, Melbourne, FL 32901, USA
| |
Collapse
|
8
|
Mueller JL, Siltanen S. The D-bar method for electrical impedance tomography-demystified. INVERSE PROBLEMS 2020; 36:093001. [PMID: 33380765 PMCID: PMC7771826 DOI: 10.1088/1361-6420/aba2f5] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Electrical impedance tomography (EIT) is an imaging modality where a patient or object is probed using harmless electric currents. The currents are fed through electrodes placed on the surface of the target, and the data consists of voltages measured at the electrodes resulting from a linearly independent set of current injection patterns. EIT aims to recover the internal distribution of electrical conductivity inside the target. The inverse problem underlying the EIT image formation task is nonlinear and severely ill-posed, and hence sensitive to modeling errors and measurement noise. Therefore, the inversion process needs to be regularized. However, traditional variational regularization methods, based on optimization, often suffer from local minima because of nonlinearity. This is what makes regularized direct (non-iterative) methods attractive for EIT. The most developed direct EIT algorithm is the D-bar method, based on Complex Geometric Optics solutions and a nonlinear Fourier transform. Variants and recent developments of D-bar methods are reviewed, and their practical numerical implementation is explained.
Collapse
Affiliation(s)
- J L Mueller
- Department of Mathematics and School of Biomedical Engineering, Colorado State University, Fort Collins, CO 80525
- Department of Mathematics and Statistics, University of Helsinki, Finland
| | - S Siltanen
- Department of Mathematics and School of Biomedical Engineering, Colorado State University, Fort Collins, CO 80525
- Department of Mathematics and Statistics, University of Helsinki, Finland
| |
Collapse
|
9
|
Capps M, Mueller JL. Reconstruction of Organ Boundaries With Deep Learning in the D-Bar Method for Electrical Impedance Tomography. IEEE Trans Biomed Eng 2020; 68:826-833. [PMID: 32746047 DOI: 10.1109/tbme.2020.3006175] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
OBJECTIVE Medical electrical impedance tomography is a non-ionizing imaging modality in which low-amplitude, low-frequency currents are applied on electrodes on the body, the resulting voltages are measured, and an inverse problem is solved to determine the conductivity distribution in the region of interest. Due the ill-posedness of the inverse problem, the boundaries of internal organs are typically blurred in the reconstructed image. METHODS A deep learning approach is introduced in the D-bar method for reconstructing a 2-D slice of the thorax to recover the boundaries of organs. This is accomplished by training a deep neural network on labeled pairs of scattering transforms and the boundaries of the organs in the data from which the transforms were computed. This allows the network to "learn" the nonlinear mapping between them by minimizing the error between the output of the network and known actual boundaries. Further, a "sparse" reconstruction is computed by fusing the results of the standard D-bar reconstruction with reconstructed organ boundaries from the neural network. RESULTS Results are shown on simulated and experimental data collected on a saline-filled tank with agar targets simulating the conductivity of the heart and lungs. CONCLUSIONS AND SIGNIFICANCE The results demonstrate that deep neural networks can successfully learn the mapping between scattering transforms and the internal boundaries of structures.
Collapse
|
10
|
Alsaker M, Mueller JL, Murthy R. DYNAMIC OPTIMIZED PRIORS FOR D-BAR RECONSTRUCTIONS OF HUMAN VENTILATION USING ELECTRICAL IMPEDANCE TOMOGRAPHY. JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS 2019; 362:276-294. [PMID: 31379404 PMCID: PMC6677406 DOI: 10.1016/j.cam.2018.07.039] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
A method of including dynamic spatial priors in the 2-D D-bar reconstruction algorithm is presented for use on time-difference reconstructions of human subject thoracic data. The conductivity values for the prior are updated at each frame in the reconstruction using an optimization method applied to the scattering transform. The updates of the dynamic spatial priors are guided by a principle component analysis of the data to determine the timepoint in the ventilatory (or cardiac) cycle. The effectiveness of the method is demonstrated on human subject ventilatory data.
Collapse
Affiliation(s)
- Melody Alsaker
- Department of Mathematics; Gonzaga University, Spokane, WA 99258 USA,
| | - Jennifer L Mueller
- Department of Mathematics and School of Biomedical Engineering, Colorado State University, CO 80523 USA,
| | - Rashmi Murthy
- Department of Mathematics, Colorado State University, CO 80523 USA,
| |
Collapse
|
11
|
Mellenthin MM, Mueller JL, de Camargo EDLB, de Moura FS, Santos TBR, Lima RG, Hamilton SJ, Muller PA, Alsaker M. The ACE1 Electrical Impedance Tomography System for Thoracic Imaging. IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT 2019; 68:3137-3150. [PMID: 33223563 PMCID: PMC7678726 DOI: 10.1109/tim.2018.2874127] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
The design and performance of the ACE1 (Active Complex Electrode) electrical impedance tomography system for single-ended phasic voltage measurements is presented. The design of the hardware and calibration procedures allow for reconstruction of conductivity and permittivity images. Phase measurement is achieved with the ACE1 active electrode circuit which measures the amplitude and phase of the voltage and the applied current at the location at which current is injected into the body. An evaluation of the system performance under typical operating conditions includes details of demodulation and calibration and an in-depth look at insightful metrics, such as signal-to-noise ratio variations during a single current pattern. Static and dynamic images of conductivity and permittivity are presented from ACE1 data collected on tank phantoms and human subjects to illustrate the system's utility.
Collapse
Affiliation(s)
| | - Jennifer L Mueller
- Department of Mathematics and School of Biomedical Engineering and the Department of Electrical and Computer Engineering, Colorado State University, CO 80523 USA
| | | | - Fernando Silva de Moura
- Center for Engineering, Modeling and Applied Social Sciences, Federal University of ABC, São Paulo, Brazil
| | | | - Raul Gonzalez Lima
- Mechanical Engineering Department, University of São Paulo, São Paulo, Brazil
| | - Sarah J Hamilton
- Department of Mathematics, Statistics, and Computer Science; Marquette University, Milwaukee, WI, 53233 USA
| | - Peter A Muller
- Department of Mathematics, Colorado State University, CO 80523 USA
| | - Melody Alsaker
- Department of Mathematics, Colorado State University, Fort Collins, C0, 80523 USA
| |
Collapse
|
12
|
Hamilton SJ, Hänninen A, Hauptmann A, Kolehmainen V. Beltrami-net: domain-independent deep D-bar learning for absolute imaging with electrical impedance tomography (a-EIT). Physiol Meas 2019; 40:074002. [PMID: 31091516 PMCID: PMC6816539 DOI: 10.1088/1361-6579/ab21b2] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
OBJECTIVE To develop, and demonstrate the feasibility of, a novel image reconstruction method for absolute electrical impedance tomography (a-EIT) that pairs deep learning techniques with real-time robust D-bar methods and examine the influence of prior information on the reconstruction. APPROACH A D-bar method is paired with a trained convolutional neural network (CNN) as a post-processing step. Training data is simulated for the network using no knowledge of the boundary shape by using an associated nonphysical Beltrami equation rather than simulating the traditional current and voltage data specific to a given domain. This allows the training data to be boundary shape independent. The method is tested on experimental data from two EIT systems (ACT4 and KIT4) with separate training sets of varying prior information. MAIN RESULTS Post-processing the D-bar images with a CNN produces significant improvements in image quality measured by structural SIMilarity indices (SSIMs) as well as relative [Formula: see text] and [Formula: see text] image errors. SIGNIFICANCE This work demonstrates that more general networks can be trained without being specific about boundary shape, a key challenge in EIT image reconstruction. The work is promising for future studies involving databases of anatomical atlases.
Collapse
Affiliation(s)
- S J Hamilton
- Department of Mathematics, Statistics, and Computer Science, Marquette University, Milwaukee, WI 53233, United States of America. Authors to whom any correspondence should be addressed
| | | | | | | |
Collapse
|
13
|
de Castro Martins T, Sato AK, de Moura FS, de Camargo EDLB, Silva OL, Santos TBR, Zhao Z, Möeller K, Amato MBP, Mueller JL, Lima RG, de Sales Guerra Tsuzuki M. A Review of Electrical Impedance Tomography in Lung Applications: Theory and Algorithms for Absolute Images. ANNUAL REVIEWS IN CONTROL 2019; 48:442-471. [PMID: 31983885 PMCID: PMC6980523 DOI: 10.1016/j.arcontrol.2019.05.002] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
Electrical Impedance Tomography (EIT) is under fast development, the present paper is a review of some procedures that are contributing to improve spatial resolution and material properties accuracy, admitivitty or impeditivity accuracy. A review of EIT medical applications is presented and they were classified into three broad categories: ARDS patients, obstructive lung diseases and perioperative patients. The use of absolute EIT image may enable the assessment of absolute lung volume, which may significantly improve the clinical acceptance of EIT. The Control Theory, the State Observers more specifically, have a developed theory that can be used for the design and operation of EIT devices. Electrode placement, current injection strategy and electrode electric potential measurements strategy should maximize the number of observable and controllable directions of the state vector space. A non-linear stochastic state observer, the Unscented Kalman Filter, is used directly for the reconstruction of absolute EIT images. Historically, difference images were explored first since they are more stable in the presence of modelling errors. Absolute images require more detailed models of contact impedance, stray capacitance and properly refined finite element mesh where the electric potential gradient is high. Parallelization of the forward program computation is necessary since the solution of the inverse problem often requires frequent solutions of the forward problem. Several reconstruction algorithms benefit by the Bayesian inverse problem approach and the concept of prior information. Anatomic and physiologic information are used to form the prior information. An already tested methodology is presented to build the prior probability density function using an ensemble of CT scans and in vivo impedance measurements. Eight absolute EIT image algorithms are presented.
Collapse
Affiliation(s)
| | - André Kubagawa Sato
- Computational Geometry Laboratory, Escola Politécnica da Universidade de São Paulo, Brazil
| | - Fernando Silva de Moura
- Universidade Federal do ABC, Center of Engineering, Modeling and Applied Social Sciences, Brazil
| | | | - Olavo Luppi Silva
- Universidade Federal do ABC, Center of Engineering, Modeling and Applied Social Sciences, Brazil
| | | | - Zhanqi Zhao
- Institute of Technical Medicine, Furtwangen University, Germany
- Department of Biomedical Engineering, Fourth Military Medical University, Xi'an, China
| | - Knut Möeller
- Institute of Technical Medicine, Furtwangen University, Germany
| | - Marcelo Brito Passos Amato
- Respiratory Intensive Care Unit, Pulmonary Division, Hospital das Clínicas, Universidade de São Paulo, Brazil
| | - Jennifer L Mueller
- Department of Mathematics, and School of Biomedical Engineering, Colorado State University, United States of America
| | - Raul Gonzalez Lima
- Department of Mechanical Engineering, Escola Politécnica da Universidade de São Paulo, Brazil
| | | |
Collapse
|
14
|
Liu D, Smyl D, Du J. A Parametric Level Set-Based Approach to Difference Imaging in Electrical Impedance Tomography. IEEE TRANSACTIONS ON MEDICAL IMAGING 2019; 38:145-155. [PMID: 30040633 DOI: 10.1109/tmi.2018.2857839] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
This paper presents a novel difference imaging approach based on the recently developed parametric level set (PLS) method for estimating the change in a target conductivity from electrical impedance tomography measurements. As in conventional difference imaging, the reconstruction of conductivity change is based on data sets measured from the surface of a body before and after the change. The key feature of the proposed approach is that the conductivity change to be reconstructed is assumed to be piecewise constant, while the geometry of the anomaly is represented by a shape-based PLS function employing Gaussian radial basis functions (GRBFs). The representation of the PLS function by using GRBF provides flexibility in describing a large class of shapes with fewer unknowns. This feature is advantageous, as it may significantly reduce the overall number of unknowns, improve the condition number of the inverse problem, and enhance the computational efficiency of the technique. To evaluate the proposed PLS-based difference imaging approach, results obtained via simulation, phantom study, and in vivo pig data are studied. We find that the proposed approach tolerates more modeling errors and leads to a significant improvement in image quality compared with the conventional linear approach.
Collapse
|
15
|
Hamilton SJ, Hauptmann A. Deep D-Bar: Real-Time Electrical Impedance Tomography Imaging With Deep Neural Networks. IEEE TRANSACTIONS ON MEDICAL IMAGING 2018; 37:2367-2377. [PMID: 29994023 DOI: 10.1109/tmi.2018.2828303] [Citation(s) in RCA: 71] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The mathematical problem for electrical impedance tomography (EIT) is a highly nonlinear ill-posed inverse problem requiring carefully designed reconstruction procedures to ensure reliable image generation. D-bar methods are based on a rigorous mathematical analysis and provide robust direct reconstructions by using a low-pass filtering of the associated nonlinear Fourier data. Similarly to low-pass filtering of linear Fourier data, only using low frequencies in the image recovery process results in blurred images lacking sharp features, such as clear organ boundaries. Convolutional neural networks provide a powerful framework for post-processing such convolved direct reconstructions. In this paper, we demonstrate that these CNN techniques lead to sharp and reliable reconstructions even for the highly nonlinear inverse problem of EIT. The network is trained on data sets of simulated examples and then applied to experimental data without the need to perform an additional transfer training. Results for absolute EIT images are presented using experimental EIT data from the ACT4 and KIT4 EIT systems.
Collapse
|
16
|
Hamilton SJ, Mueller JL, Santos TR. Robust computation in 2D absolute EIT (a-EIT) using D-bar methods with the 'exp' approximation. Physiol Meas 2018; 39:064005. [PMID: 29846182 DOI: 10.1088/1361-6579/aac8b1] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
OBJECTIVE Absolute images have important applications in medical electrical impedance tomography (EIT) imaging, but the traditional minimization and statistical based computations are very sensitive to modeling errors and noise. In this paper, it is demonstrated that D-bar reconstruction methods for absolute EIT are robust to such errors. APPROACH The effects of errors in domain shape and electrode placement on absolute images computed with 2D D-bar reconstruction algorithms are studied on experimental data. MAIN RESULTS It is demonstrated with tank data from several EIT systems that these methods are quite robust to such modeling errors, and furthermore the artefacts arising from such modeling errors are similar to those occurring in classic time-difference EIT imaging. SIGNIFICANCE This study is promising for clinical applications where absolute EIT images are desirable but previously thought impossible.
Collapse
Affiliation(s)
- S J Hamilton
- Department of Mathematics, Statistics, and Computer Science, Marquette University, Milwaukee, WI 53233, United States of America
| | | | | |
Collapse
|
17
|
Mueller JL, Muller P, Mellenthin M, Murthy R, Capps M, Alsaker M, Deterding R, Sagel SD, DeBoer E. Estimating regions of air trapping from electrical impedance tomography data. Physiol Meas 2018; 39:05NT01. [PMID: 29726838 PMCID: PMC6015736 DOI: 10.1088/1361-6579/aac295] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
OBJECTIVE Electrical impedance tomography (EIT) has been shown to be a viable non-invasive, bedside imaging modality to monitor lung function. This paper introduces a method for identifying regions of air trapping from EIT data collected during tidal breathing and breath-holding maneuvers. APPROACH Ventilation-perfusion index maps are computed from dynamic EIT images. These maps are then used to identify regions of air trapping in the area of the lung as regions that are poorly ventilated but well perfused throughout the breathing and cardiac cycles. These EIT-identified regions are then compared with independently identified regions of low attenuation, or air trapping, on chest CT. Results of this method are demonstrated in two children with cystic fibrosis and on a healthy control subject. MAIN RESULTS In both CF children, the EIT-identified regions of air trapping matched the regions indicated from the chest CT. The EIT-based method is only validated with CT scans within 4 cm of the chest cross-section defined by the electrode plane. SIGNIFICANCE The results indicate the potential use of EIT-derived ventilation-perfusion index maps as a non-invasive method for identifying regions of air trapping.
Collapse
Affiliation(s)
- Jennifer L Mueller
- Department of Mathematics, Colorado State University, Fort Collins, CO, United States of America. School of Biomedical Engineering, Colorado State University, Fort Collins, CO, United States of America
| | | | | | | | | | | | | | | | | |
Collapse
|
18
|
Muller PA, Mueller JL, Mellenthin M, Murthy R, Capps M, Wagner BD, Alsaker M, Deterding R, Sagel SD, Hoppe J. Evaluation of surrogate measures of pulmonary function derived from electrical impedance tomography data in children with cystic fibrosis. Physiol Meas 2018; 39:045008. [PMID: 29565263 DOI: 10.1088/1361-6579/aab8c4] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
OBJECTIVE Lung function monitoring by spirometry plays a critical role in the clinical care of pediatric cystic fibrosis (CF) patients, but many young children are unable to perform spirometry, and the outputs are often normal even in the presence of lung disease. Measures derived from electrical impedance tomography (EIT) images were studied for their utility as potential surrogates for spirometry in CF patients and to assess response to intravenous antibiotic treatment for acute pulmonary exacerbations (PEx) in a subset of patients. APPROACH EIT data were collected on 35 subjects (21 with CF, 14 healthy controls, 8 CF patients pre- and post-treatment for an acute PEx) ages 2 to 20 years during tidal breathing and also concurrently with spirometry on subjects over age 8. EIT-derived measures of FEV1, FVC, and FEV1/FVC were computed globally and regionally from dynamic EIT images. MAIN RESULTS Global EIT-derived FEV1/FVC showed good correlation with spirometry FEV1/FVC values (r = 0.54, p = 0.01), and were able to distinguish between the groups (p = 0.01). Lung heterogeneity was assessed through the spatial coefficient of variation (CV) of EIT difference images between key time points, and the CVs for EIT-derived FEV1 and FVC showed significant correlation with the CV for tidal breathing (r = 0.47, p = 0.01 and r = 0.50, p = 0.01, respectively). Global EIT-derived FEV1/FVC was better able to distinguish between groups than spirometry FEV1 (F-values 776.5 and 146.3, respectively, p < 0.01.) The same held true for the CVs for EIT-derived FEV1, FVC, and tidal breathing (F-values 215.93, 193.89, 204.57, respectively, p < 0.01). SIGNIFICANCE The strong correlation between the CVs for tidal breathing, FEV1, and FVC, and the statistically significant ability of CV for tidal breathing to distinguish between healthy subjects and CF patients, and between the studied CF disease states suggests that the CV may be useful for measuring the extent and severity of structural lung disease.
Collapse
Affiliation(s)
- Peter A Muller
- Department of Mathematics & Statistics, Villanova University, PA, United States of America. was at Department of Mathematics, Colorado State University, CO, United States of America
| | | | | | | | | | | | | | | | | | | |
Collapse
|
19
|
Muller PA, Mueller JL, Mellenthin MM. Real-Time Implementation of Calderón's Method on Subject-Specific Domains. IEEE TRANSACTIONS ON MEDICAL IMAGING 2017; 36:1868-1875. [PMID: 28436855 DOI: 10.1109/tmi.2017.2695893] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
A real-time implementation of Calderón's method for the reconstruction of a 2-D conductivity from electrical impedance tomography data is presented, in which domain-specific modeling is taken into account. This is the first implementation of Calderón's method that accounts for correct modeling of non-symmetric domain boundaries in image reconstruction. The domain-specific Calderón's method is derived and reconstructions from experimental tank data are presented, quantifying the distortion when correct modeling is not included in the reconstruction algorithm. Reconstructions from human subject volunteers are presented, demonstrating the method's effectiveness for imaging changes due to ventilation and perfusion in the human thorax.
Collapse
|
20
|
Hamilton SJ. EIT Imaging of admittivities with a D-bar method and spatial prior: experimental results for absolute and difference imaging. Physiol Meas 2017; 38:1176-1192. [PMID: 28530208 DOI: 10.1088/1361-6579/aa63d7] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Electrical impedance tomography (EIT) is an emerging imaging modality that uses harmless electrical measurements taken on electrodes at a body's surface to recover information about the internal electrical conductivity and or permittivity. The image reconstruction task of EIT is a highly nonlinear inverse problem that is sensitive to noise and modeling errors making the image reconstruction task challenging. D-bar methods solve the nonlinear problem directly, bypassing the need for detailed and time-intensive forward models, to provide absolute (static) as well as time-difference EIT images. Coupling the D-bar methodology with the inclusion of high confidence a priori data results in a noise-robust regularized image reconstruction method. In this work, the a priori D-bar method for complex admittivities is demonstrated effective on experimental tank data for absolute imaging for the first time. Additionally, the method is adjusted for, and tested on, time-difference imaging scenarios. The ability of the method to be used for conductivity, permittivity, absolute as well as time-difference imaging provides the user with great flexibility without a high computational cost.
Collapse
Affiliation(s)
- S J Hamilton
- Department of Mathematics, Statistics, and Computer Science, Marquette University, Milwaukee, WI 53233, United States of America
| |
Collapse
|
21
|
Hamilton SJ, Mueller JL, Alsaker M. Incorporating a Spatial Prior into Nonlinear D-Bar EIT Imaging for Complex Admittivities. IEEE TRANSACTIONS ON MEDICAL IMAGING 2017; 36:457-466. [PMID: 28114061 PMCID: PMC5384275 DOI: 10.1109/tmi.2016.2613511] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
Electrical Impedance Tomography (EIT) aims to recover the internal conductivity and permittivity distributions of a body from electrical measurements taken on electrodes on the surface of the body. The reconstruction task is a severely ill-posed nonlinear inverse problem that is highly sensitive to measurement noise and modeling errors. Regularized D-bar methods have shown great promise in producing noise-robust algorithms by employing a low-pass filtering of nonlinear (nonphysical) Fourier transform data specific to the EIT problem. Including prior data with the approximate locations of major organ boundaries in the scattering transform provides a means of extending the radius of the low-pass filter to include higher frequency components in the reconstruction, in particular, features that are known with high confidence. This information is additionally included in the system of D-bar equations with an independent regularization parameter from that of the extended scattering transform. In this paper, this approach is used in the 2-D D-bar method for admittivity (conductivity as well as permittivity) EIT imaging. Noise-robust reconstructions are presented for simulated EIT data on chest-shaped phantoms with a simulated pneumothorax and pleural effusion. No assumption of the pathology is used in the construction of the prior, yet the method still produces significant enhancements of the underlying pathology (pneumothorax or pleural effusion) even in the presence of strong noise.
Collapse
|
22
|
Alsaker M, Jane Hamilton S, Hauptmann A. A direct D-bar method for partial boundary data electrical impedance tomography with a priori information. ACTA ACUST UNITED AC 2017. [DOI: 10.3934/ipi.2017020] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
23
|
Mellenthin MM, Mueller JL, Bueno de Camargo EDL, Silva de Moura F, Hamilton SJ, Gonzalez Lima R. The ACE1 thoracic Electrical Impedance Tomography system for ventilation and perfusion. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2016; 2015:4073-6. [PMID: 26737189 DOI: 10.1109/embc.2015.7319289] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Electrical Impedance Tomography (EIT) is a technique which can image the varying electrical properties of biological tissues. For clinical use of EIT, it can be advantageous to know both tissue conductivity and permittivity. Presented is the hardware design for the pairwise current injection active complex electrode (ACE1) EIT system which measures phasic voltages for conductivity and permittivity image reconstruction. In this system, alternating current is injected on electrodes on the boundary of a domain and single-ended voltage measurements are used in image reconstructions of the domain's interior and in calculating the current applied at the electrodes.
Collapse
|
24
|
Structural-functional lung imaging using a combined CT-EIT and a Discrete Cosine Transformation reconstruction method. Sci Rep 2016; 6:25951. [PMID: 27181695 PMCID: PMC4867600 DOI: 10.1038/srep25951] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2015] [Accepted: 04/20/2016] [Indexed: 12/14/2022] Open
Abstract
Lung EIT is a functional imaging method that utilizes electrical currents to reconstruct images of conductivity changes inside the thorax. This technique is radiation free and applicable at the bedside, but lacks of spatial resolution compared to morphological imaging methods such as X-ray computed tomography (CT). In this article we describe an approach for EIT image reconstruction using morphologic information obtained from other structural imaging modalities. This leads to recon- structed images of lung ventilation that can easily be superimposed with structural CT or MRI images, which facilitates image interpretation. The approach is based on a Discrete Cosine Transformation (DCT) of an image of the considered transversal thorax slice. The use of DCT enables reduction of the dimensionality of the reconstruction and ensures that only conductivity changes of the lungs are reconstructed and displayed. The DCT based approach is well suited to fuse morphological image information with functional lung imaging at low computational costs. Results on simulated data indicate that this approach preserves the morphological structures of the lungs and avoids blurring of the solution. Images from patient measurements reveal the capabilities of the method and demonstrate benefits in possible applications.
Collapse
|