Reconstructions of chest phantoms by the D-bar method for electrical impedance tomography

ACT5 Electrical Impedance Tomography System

OBJECTIVE

This article introduces the Adaptive Current Tomograph 5 (ACT5) Electrical Impedance Tomography (EIT) system. ACT5 is a 32 electrode applied-current multiple-source EIT system that can display real-time images of conductivity and susceptivity at 27 frames per second. The adaptive current sources in ACT5 can apply fully programmable current patterns with frequencies varying from 5 kHz to 500 kHz. The system also displays real-time ECG readings during the EIT imaging process.

METHODS

The hardware and software design and specifications are presented, including the current source design, FPGA hardware, safety features, calibration, and shunt impedance measurement.

RESULTS

Images of conductivity and susceptivity are presented from ACT5 data collected on tank phantoms and a human subject illustrating the system's ability to provide real-time images of pulsatile perfusion and ECG traces.

SIGNIFICANCE

The portability, high signal-to-noise ratio, and flexibility of applied currents over a wide range of frequencies enable this instrument to be used to obtain useful human subject data with relative clinical ease.

Domain independent post-processing with graph U-nets: applications to electrical impedance tomographic imaging⋆

Objective.To extend the highly successful U-Net Convolutional Neural Network architecture, which is limited to rectangular pixel/voxel domains, to a graph-based equivalent that works flexibly on irregular meshes; and demonstrate the effectiveness on electrical impedance tomography (EIT).Approach.By interpreting the irregular mesh as a graph, we develop a graph U-Net with new cluster pooling and unpooling layers that mimic the classic neighborhood based max-pooling important for imaging applications.Mainresults.The proposed graph U-Net is shown to be flexible and effective for improving early iterate total variation (TV) reconstructions from EIT measurements, using as little as the first iteration. The performance is evaluated for simulated data, and on experimental data from three measurement devices with different measurement geometries and instrumentations. We successfully show that such networks can be trained with a simple two-dimensional simulated training set, and generalize to very different domains, including measurements from a three-dimensional device and subsequent 3D reconstructions.Significance.As many inverse problems are solved on irregular (e.g. finite element) meshes, the proposed graph U-Net and pooling layers provide the added flexibility to process directly on the computational mesh. Post-processing an early iterate reconstruction greatly reduces the computational cost which can become prohibitive in higher dimensions with dense meshes. As the graph structure is independent of 'dimension', the flexibility to extend networks trained on 2D domains to 3D domains offers a possibility to further reduce computational cost in training.

Image Reconstruction Using Supervised Learning in Wearable Electrical Impedance Tomography of the Thorax

Electrical impedance tomography (EIT) is a non-invasive technique for visualizing the internal structure of a human body. Capacitively coupled electrical impedance tomography (CCEIT) is a new contactless EIT technique that can potentially be used as a wearable device. Recent studies have shown that a machine learning-based approach is very promising for EIT image reconstruction. Most of the studies concern models containing up to 22 electrodes and focus on using different artificial neural network models, from simple shallow networks to complex convolutional networks. However, the use of convolutional networks in image reconstruction with a higher number of electrodes requires further investigation. In this work, two different architectures of artificial networks were used for CCEIT image reconstruction: a fully connected deep neural network and a conditional generative adversarial network (cGAN). The training dataset was generated by the numerical simulation of a thorax phantom with healthy and illness-affected lungs. Three kinds of illnesses, pneumothorax, pleural effusion, and hydropneumothorax, were modeled using the electrical properties of the tissues. The thorax phantom included the heart, aorta, spine, and lungs. The sensor with 32 area electrodes was used in the numerical model. The ECTsim custom-designed toolbox for Matlab was used to solve the forward problem and measurement simulation. Two artificial neural networks were trained with supervision for image reconstruction. Reconstruction quality was compared between those networks and one-step algebraic reconstruction methods such as linear back projection and pseudoinverse with Tikhonov regularization. This evaluation was based on pixel-to-pixel metrics such as root-mean-square error, structural similarity index, 2D correlation coefficient, and peak signal-to-noise ratio. Additionally, the diagnostic value measured by the ROC AUC metric was used to assess the image quality. The results showed that obtaining information about regional lung function (regions affected by pneumothorax or pleural effusion) is possible using image reconstruction based on supervised learning and deep neural networks in EIT. The results obtained using cGAN are strongly better than those obtained using a fully connected network, especially in the case of noisy measurement data. However, diagnostic value estimation showed that even algebraic methods allow us to obtain satisfactory results.

Implicit Solutions of the Electrical Impedance Tomography Inverse Problem in the Continuous Domain with Deep Neural Networks

Electrical impedance tomography (EIT) is a non-invasive imaging modality used for estimating the conductivity of an object Ω from boundary electrode measurements. In recent years, researchers achieved substantial progress in analytical and numerical methods for the EIT inverse problem. Despite the success, numerical instability is still a major hurdle due to many factors, including the discretization error of the problem. Furthermore, most algorithms with good performance are relatively time consuming and do not allow real-time applications. In our approach, the goal is to separate the unknown conductivity into two regions, namely the region of homogeneous background conductivity and the region of non-homogeneous conductivity. Therefore, we pose and solve the problem of shape reconstruction using machine learning. We propose a novel and simple jet intriguing neural network architecture capable of solving the EIT inverse problem. It addresses previous difficulties, including instability, and is easily adaptable to other ill-posed coefficient inverse problems. That is, the proposed model estimates the probability for a point of whether the conductivity belongs to the background region or to the non-homogeneous region on the continuous space Rd∩Ω with d∈{2,3}. The proposed model does not make assumptions about the forward model and allows for solving the inverse problem in real time. The proposed machine learning approach for shape reconstruction is also used to improve gradient-based methods for estimating the unknown conductivity. In this paper, we propose a piece-wise constant reconstruction method that is novel in the inverse problem setting but inspired by recent approaches from the 3D vision community. We also extend this method into a novel constrained reconstruction method. We present extensive numerical experiments to show the performance of the architecture and compare the proposed method with previous analytic algorithms, mainly the monotonicity-based shape reconstruction algorithm and iteratively regularized Gauss-Newton method.

Fast absolute 3D CGO-based electrical impedance tomography on experimental tank data

Objective.To present the first 3D CGO-based absolute EIT reconstructions from experimental tank data.Approach.CGO-based methods for absolute EIT imaging are compared to traditional TV regularized non-linear least squares reconstruction methods. Additional robustness testing is performed by considering incorrect modeling of domain shape.Main Results.The CGO-based methods are fast, and show strong robustness to incorrect domain modeling comparable to classic difference EIT imaging and fewer boundary artefacts than the TV regularized non-linear least squares reference reconstructions.Significance.This work is the first to demonstrate fully 3D CGO-based absolute EIT reconstruction on experimental data and also compares to TV-regularized absolute reconstruction. The speed (1-5 s) and quality of the reconstructions is encouraging for future work in absolute EIT.

Improved resolution of D-bar images of ventilation using a Schur complement property and an anatomical atlas

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.

Adaptively Regularized Bases-Expansion Subspace Optimization Methods for Electrical Impedance Tomography

OBJECTIVE

In this work, to deal with the difficulties in choosing regularization weighting parameters and low spatial resolution problems in difference electrical impedance tomography (EIT), we propose two adaptively regularized bases-expansion subspace optimization methods (AR-BE-SOMs).

METHODS

Firstly, an adaptive L1-norm based total variation (TV) regularization is introduced under the framework of BE-SOM. Secondly, besides the additive regularization, an adaptive weighted L2-norm multiplicative regularization is further proposed. The regularized objective functions are optimized by conjugate gradient method, where the unknowns in both methods are update alternatively between induced contrast current (ICC) and conductivity domain.

CONCLUSION

Both numerical and experimental tests are conducted to validate the proposed methods, where AR-BE-SOMs show better edge-preserving, anti-noise performance, lower relative errors, and higher structure similarity indexes than BE-SOM.

SIGNIFICANCE

Unlike the common regularization techniques in EIT, the proposed regularization factors can be obtained adaptively during the optimization process. More importantly, ARBE-SOMs perform well in reconstructions of some challenging geometry with sharp corners such as the heart and lung phantoms, deformation phantoms, triangles and even rectangles. It is expected that the proposed AR-BE-SOMs will find their applications for high-quality lung health monitoring and other clinical applications.

Improvement of performance and sensitivity of 2D and 3D image reconstruction in EIT using EFG forward model

This paper presents a pure element-free Galerkin method (EFGM) forward model for image reconstruction in 2D and 3D electrical impedance tomography (EIT) using an adaptive current injection method. In EIT systems with the adapting current injection method, both static and dynamic images can be reconstructed; however, determination of electrode contact impedances in the complete electrode model is difficult and the Gap model is used. In this paper, in the EIT forward problem a weak form functional based on the Gap model and a pure EFGM approach are developed, and in the EIT inverse problem, Jacobian matrix is computed by the EFGM, and a fast integration technique is introduced to calculate the entries of the Jacobian matrix within an adequate computation time. The influence of increasing the density of nodes at and near the electrodes with steep electric potential gradients on the accuracy of FEM and EFGM forward solutions is investigated, and the performance of the image reconstruction algorithm with the proposed fast integration technique is examined. The numerical results reveal that the proposed EFGM forward model with the fast integration technique has an efficient performance both in terms of mean relative imaging errors and computational time.

A Rotational Invariant Neural Network for Electrical Impedance Tomography Imaging without Reference Voltage: RF-REIM-NET

Background: Electrical Impedance Tomography (EIT) is a radiation-free technique for image reconstruction. However, as the inverse problem of EIT is non-linear and ill-posed, the reconstruction of sharp conductivity images poses a major problem. With the emergence of artificial neural networks (ANN), their application in EIT has recently gained interest. Methodology: We propose an ANN that can solve the inverse problem without the presence of a reference voltage. At the end of the ANN, we reused the dense layers multiple times, considering that the EIT exhibits rotational symmetries in a circular domain. To avoid bias in training data, the conductivity range used in the simulations was greater than expected in measurements. We also propose a new method that creates new data samples from existing training data. Results: We show that our ANN is more robust with respect to noise compared with the analytical Gauss–Newton approach. The reconstruction results for EIT phantom tank measurements are also clearer, as ringing artefacts are less pronounced. To evaluate the performance of the ANN under real-world conditions, we perform reconstructions on an experimental pig study with computed tomography for comparison. Conclusions: Our proposed ANN can reconstruct EIT images without the need of a reference voltage.

Graph Convolutional Networks for Model-Based Learning in Nonlinear Inverse Problems

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.

3D ELECTRICAL IMPEDANCE TOMOGRAPHY RECONSTRUCTIONS FROM SIMULATED ELECTRODE DATA USING DIRECT INVERSION texp AND CALDERÓN METHODS

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.

Model Selection Based Algorithm in Neonatal Chest EIT

This paper presents a new method for selecting a patient specific forward model to compensate for anatomical variations in electrical impedance tomography (EIT) monitoring of neonates. The method uses a combination of shape sensors and absolute reconstruction. It takes advantage of a probabilistic approach which automatically selects the best estimated forward model fit from pre-stored library models. Absolute/static image reconstruction is performed as the core of the posterior probability calculations. The validity and reliability of the algorithm in detecting a suitable model in the presence of measurement noise is studied with simulated and measured data from 11 patients. The paper also demonstrates the potential improvements on the clinical parameters extracted from EIT images by considering a unique case study with a neonate patient undergoing computed tomography imaging as clinical indication prior to EIT monitoring. Two well-known image reconstruction techniques, namely GREIT and tSVD, are implemented to create the final tidal images. The impacts of appropriate model selection on the clinical extracted parameters such as center of ventilation and silent spaces are investigated. The results show significant improvements to the final reconstructed images and more importantly to the clinical EIT parameters extracted from the images that are crucial for decision-making and further interventions.

Introduction of Sample Based Prior into the D-Bar Method Through a Schur Complement Property

Electrical impedance tomography (EIT) is a non-invasive medical imaging technique in which images of the conductivity in a region of interest in the body are computed from measurements of voltages on electrodes arising from low-frequency, low-amplitude applied currents. Mathematically, the inverse conductivity problem is nonlinear and ill-posed, and the reconstructions have characteristically low spatial resolution. One approach to improve the spatial resolution of EIT images is to include anatomically and physiologically-based prior information in the reconstruction algorithm. Statistical inversion theory provides a means of including prior information from a representative sample population. In this paper, a method is proposed to introduce statistical prior information into the D-bar method based on Schur complement properties. The method presents an improvement of the image obtained by the D-bar method by maximizing the conditional probability density function of an image that is consistent with a prior information and the model, given a D-bar image computed from the voltage measurements. Experimental phantoms show an improved spatial resolution by the use of the proposed method for the D-bar image reconstructions.

A second order Calderón's method with a correction term and a priori information

Calderón's method is a direct linearized reconstruction method for the inverse conductivity problem with the attribute that it can provide absolute images of both conductivity and permittivity with no need for forward modeling. In this work, an explicit relationship between Calderón's method and the D-bar method is provided, facilitating a "higher-order" Calderón's method in which a correction term is included, derived from the relationship to the D-bar method. Furthermore, a method of including a spatial prior is provided. These advances are demonstrated on simulated data and on tank data collected with the ACE1 EIT system.

The D-bar method for electrical impedance tomography-demystified

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.

Induced-Current Learning Method for Nonlinear Reconstructions in Electrical Impedance Tomography

Electrical impedance tomography (EIT) is an attractive technique that aims to reconstruct the unknown electrical property in a domain from the surface electrical measurements. In this work, the induced-current learning method (ICLM) is proposed to solve nonlinear electrical impedance tomography (EIT) problems. Specifically, the cascaded end-to-end convolutional neural network (CEE-CNN) architecture is designed to implement the ICLM. The CEE-CNN greatly decreases the nonlinearities in EIT problems by designing a combined objective function and introducing multiple labels. A noticeable characteristic of the proposed CNN scheme is that the input parameters are chosen as both induced contrast current (ICC) and the updated electrical field from a spectral analysis and the output is chosen as ICC, which is fundamentally different from prevailing CNN schemes. Further, several skip connections are introduced to focus on learning only the unknown part of ICC. ICLM is verified with both numerical and experimental tests on typical EIT problems, and it is found that ICLM is able to solve typical EIT problems in less than 1 second with high image qualities. More importantly, it is also highly robust to measurement noises and modeling errors, such as inaccurate boundary data.

Improved Imaging Resolution of Electrical Impedance Tomography Using Artificial Neural Networks for Image Reconstruction

Electrical impedance tomography (EIT) is a noninvasive and non-radiative medical imaging technique based on detecting the inhomogeneous electrical properties of the tissue. The inverse problem of EIT is a highly nonlinear ill-posed problem, which is the main reason that affects image quality. Our goal is to solve the EIT inverse problem using the nonlinear mapping properties of artificial neural networks (ANNs) and convolutional neural networks (CNNs). In this paper, the adaptive moment estimation (ADAM) optimization method and mean-square-error (MSE) function are used to train an ANN to solve the inverse problem and a CNN to process the ANN image. The networks are trained on datasets of simulated data, and tested on datasets of simulated data and experimental data. Results for time-difference EIT (td-EIT) images are presented using simulated EIT data from EIDORS and experimental EIT data from our EIT systems. The results are used to compare the proposed method with the one-step Gauss-Newton linear method and RBFNN method. The proposed method offers improved resolution (RES), low position error (PE) and excellent artefact removal compared to the existing methods. The experimental results show that our method can improve the RES by 50 to 70 percent and reduce the PE by 60 to 70 percent. The improvements in RES and processing speed are essential for clinical EIT measurement of dynamic physiological processes.

Towards a thoracic conductive phantom for EIT

Phantom experiments are a crucial step for testing new hardware or imaging algorithms for electrical impedance tomography (EIT) studies. However, constructing an accurate phantom for EIT research remains critical; some studies have attempted to model the skull and breasts, and even fewer, as yet, have considered the thorax. In this study, a critical comparison between the electrical properties (impedance) of three materials is undertaken: a polyurethane foam, a silicone mixture and a thermoplastic polyurethane filament. The latter was identified as the most promising material and adopted for the development of a flexible neonatal torso. The validation is performed by the EIT image reconstruction of the air filled cavities, which mimic the lung regions. The methodology is reproducible for the creation of any phantom that requires a slight flexibility.

DYNAMIC OPTIMIZED PRIORS FOR D-BAR RECONSTRUCTIONS OF HUMAN VENTILATION USING ELECTRICAL IMPEDANCE TOMOGRAPHY

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.

The ACE1 Electrical Impedance Tomography System for Thoracic Imaging

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.

Beltrami-net: domain-independent deep D-bar learning for absolute imaging with electrical impedance tomography (a-EIT)

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.

A Review of Electrical Impedance Tomography in Lung Applications: Theory and Algorithms for Absolute Images

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.

Comparing D-bar and common regularization-based methods for electrical impedance tomography

Objective:

To compare D-bar difference reconstruction with regularized linear reconstruction in electrical impedance tomography.

Approach:

A standard regularized linear approach using a Laplacian penalty and the GREIT method for comparison to the D-bar difference images. Simulated data was generated using a circular phantom with small objects, as well as a ‘Pac-Man’ shaped conductivity target. An L-curve method was used for parameter selection in both D-bar and the regularized methods.

Main results:

We found that the D-bar method had a more position independent point spread function, was less sensitive to errors in electrode position and behaved differently with respect to additive noise than the regularized methods.

Significance:

The results allow a novel pathway between traditional and D-bar algorithm comparison.

Dominant-Current Deep Learning Scheme for Electrical Impedance Tomography

OBJECTIVE

Deep learning has recently been applied to electrical impedance tomography (EIT) imaging. Nevertheless, there are still many challenges that this approach has to face, e.g., targets with sharp corners or edges cannot be well recovered when using circular inclusion training data. This paper proposes an iterative-based inversion method and a convolutional neural network (CNN) based inversion method to recover some challenging inclusions such as triangular, rectangular, or lung shapes, where the CNN-based method uses only random circle or ellipse training data.

METHODS

First, the iterative method, i.e., bases-expansion subspace optimization method (BE-SOM), is proposed based on a concept of induced contrast current (ICC) with total variation regularization. Second, the theoretical analysis of BE-SOM and the physical concepts introduced there motivate us to propose a dominant-current deep learning scheme for EIT imaging, in which dominant parts of ICC are utilized to generate multi-channel inputs of CNN.

RESULTS

The proposed methods are tested with both numerical and experimental data, where several realistic phantoms including simulated pneumothorax and pleural effusion pathologies are also considered.

CONCLUSIONS AND SIGNIFICANCE

Significant performance improvements of the proposed methods are shown in reconstructing targets with sharp corners or edges. It is also demonstrated that the proposed methods are capable of fast, stable, and high-quality EIT imaging, which is promising in providing quantitative images for potential clinical applications.

Deep D-Bar: Real-Time Electrical Impedance Tomography Imaging With Deep Neural Networks

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.

Robust computation in 2D absolute EIT (a-EIT) using D-bar methods with the 'exp' approximation

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.

Real-Time Implementation of Calderón's Method on Subject-Specific Domains

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.

EIT Imaging of admittivities with a D-bar method and spatial prior: experimental results for absolute and difference imaging

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.

Incorporating a Spatial Prior into Nonlinear D-Bar EIT Imaging for Complex Admittivities

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.

Solution of the nonlinear inverse scattering problem by T-matrix completion. I. Theory

We propose a conceptually different method for solving nonlinear inverse scattering problems (ISPs) such as are commonly encountered in tomographic ultrasound imaging, seismology, and other applications. The method is inspired by the theory of nonlocality of physical interactions and utilizes the relevant formalism. We formulate the ISP as a problem whose goal is to determine an unknown interaction potential V from external scattering data. Although we seek a local (diagonally dominated) V as the solution to the posed problem, we allow V to be nonlocal at the intermediate stages of iterations. This allows us to utilize the one-to-one correspondence between V and the T matrix of the problem. Here it is important to realize that not every T corresponds to a diagonal V and we, therefore, relax the usual condition of strict diagonality (locality) of V. An iterative algorithm is proposed in which we seek T that is (i) compatible with the measured scattering data and (ii) corresponds to an interaction potential V that is as diagonally dominated as possible. We refer to this algorithm as to the data-compatible T-matrix completion. This paper is Part I in a two-part series and contains theory only. Numerical examples of image reconstruction in a strongly nonlinear regime are given in Part II [H. W. Levinson and V. A. Markel, Phys. Rev. E 94, 043318 (2016)10.1103/PhysRevE.94.043318]. The method described in this paper is particularly well suited for very large data sets that become increasingly available with the use of modern measurement techniques and instrumentation.

Interface and permittivity simultaneous reconstruction in electrical capacitance tomography based on boundary and finite-elements coupling method

Electrical capacitance tomography (ECT) is a non-destructive detection technique for imaging the permittivity distributions inside an observed domain from the capacitances measurements on its boundary. Owing to its advantages of non-contact, non-radiation, high speed and low cost, ECT is promising in the measurements of many industrial or biological processes. However, in the practical industrial or biological systems, a deposit is normally seen in the inner wall of its pipe or vessel. As the actual region of interest (ROI) of ECT is surrounded by the deposit layer, the capacitance measurements become weakly sensitive to the permittivity perturbation occurring at the ROI. When there is a major permittivity difference between the deposit and the ROI, this kind of shielding effect is significant, and the permittivity reconstruction becomes challenging. To deal with the issue, an interface and permittivity simultaneous reconstruction approach is proposed. Both the permittivity at the ROI and the geometry of the deposit layer are recovered using the block coordinate descent method. The boundary and finite-elements coupling method is employed to improve the computational efficiency. The performance of the proposed method is evaluated with the simulation tests. This article is part of the themed issue 'Supersensing through industrial process tomography'.

EIT image reconstruction based on a hybrid FE-EFG forward method and the complete-electrode model

This paper presents the application of the hybrid finite element-element free Galerkin (FE-EFG) method for the forward and inverse problems of electrical impedance tomography (EIT). The proposed method is based on the complete electrode model. Finite element (FE) and element-free Galerkin (EFG) methods are accurate numerical techniques. However, the FE technique has meshing task problems and the EFG method is computationally expensive. In this paper, the hybrid FE-EFG method is applied to take both advantages of FE and EFG methods, the complete electrode model of the forward problem is solved, and an iterative regularized Gauss-Newton method is adopted to solve the inverse problem. The proposed method is applied to compute Jacobian in the inverse problem. Utilizing 2D circular homogenous models, the numerical results are validated with analytical and experimental results and the performance of the hybrid FE-EFG method compared with the FE method is illustrated. Results of image reconstruction are presented for a human chest experimental phantom.

A shape-based quality evaluation and reconstruction method for electrical impedance tomography

Linear methods of reconstruction play an important role in medical electrical impedance tomography (EIT) and there is a wide variety of algorithms based on several assumptions. With the Graz consensus reconstruction algorithm for EIT (GREIT), a novel linear reconstruction algorithm as well as a standardized framework for evaluating and comparing methods of reconstruction were introduced that found widespread acceptance in the community. In this paper, we propose a two-sided extension of this concept by first introducing a novel method of evaluation. Instead of being based on point-shaped resistivity distributions, we use 2759 pairs of real lung shapes for evaluation that were automatically segmented from human CT data. Necessarily, the figures of merit defined in GREIT were adjusted. Second, a linear method of reconstruction that uses orthonormal eigenimages as training data and a tunable desired point spread function are proposed. Using our novel method of evaluation, this approach is compared to the classical point-shaped approach. Results show that most figures of merit improve with the use of eigenimages as training data. Moreover, the possibility of tuning the reconstruction by modifying the desired point spread function is shown. Finally, the reconstruction of real EIT data shows that higher contrasts and fewer artifacts can be achieved in ventilation- and perfusion-related images.

Linearly constrained minimum variance spatial filtering for localization of conductivity changes in electrical impedance tomography

We localize dynamic electrical conductivity changes and reconstruct their time evolution introducing the spatial filtering technique to electrical impedance tomography (EIT). More precisely, we use the unit-noise-gain constrained variation of the distortionless-response linearly constrained minimum variance spatial filter. We address the effects of interference and the use of zero gain constraints. The approach is successfully tested in simulated and real tank phantoms. We compute the position error and resolution to compare the localization performance of the proposed method with the one-step Gauss-Newton reconstruction with Laplacian prior. We also study the effects of sensor position errors. Our results show that EIT spatial filtering is useful for localizing conductivity changes of relatively small size and for estimating their time-courses. Some potential dynamic EIT applications such as acute ischemic stroke detection and neuronal activity localization may benefit from the higher resolution of spatial filters as compared to conventional tomographic reconstruction algorithms.

Direct 2-D reconstructions of conductivity and permittivity from EIT data on a human chest

A novel direct D-bar reconstruction algorithm is presented for reconstructing a complex conductivity distribution from 2-D EIT data. The method is applied to simulated data and archival human chest data. Permittivity reconstructions with the aforementioned method and conductivity reconstructions with the previously existing nonlinear D-bar method for real-valued conductivities depicting ventilation and perfusion in the human chest are presented. This constitutes the first fully nonlinear D-bar reconstructions of human chest data and the first D-bar permittivity reconstructions of experimental data. The results of the human chest data reconstructions are compared on a circular domain versus a chest-shaped domain.

Impedance imaging with first-order TV regularization

EIT problem is a typical inverse problem with serious ill-posedness. In general, regularization techniques are necessary for such ill-posed inverse problems. To overcome ill-posedness, the total variation (TV) regularization is widely used and it is also successfully applied to EIT. For realtime monitoring, a fast and robust image reconstruction algorithm is required. By exploiting recent advances in optimization, we propose a first-order TV algorithm for EIT, which simply consists of matrix-vector multiplications and in which the sparse structure of the system can be easily exploited. Furthermore, a typical smoothing parameter to overcome nondifferentibility of the TV term is not needed and a closed form solution can be applied in part using soft thresholding. It shows a fast reconstruction in the beginning. Numerical experiments using simulated data and real experimental data support our claim.

A Real-time D-bar Algorithm for 2-D Electrical Impedance Tomography Data

The aim of this paper is to show the feasibility of the D-bar method for real-time 2-D EIT reconstructions. A fast implementation of the D-bar method for reconstructing conductivity changes on a 2-D chest-shaped domain is described. Cross-sectional difference images from the chest of a healthy human subject are presented, demonstrating what can be achieved in real time. The images constitute the first D-bar images from EIT data on a human subject collected on a pairwise current injection system.

Optical imaging of phantoms from real data by an approximately globally convergent inverse algorithm

A numerical method for an inverse problem for an elliptic equation with the running source at multiple positions is presented. This algorithm does not rely on a good first guess for the solution. The so-called "approximate global convergence" property of this method is shown here. The performance of the algorithm is verified on real data for Diffusion Optical Tomography. Direct applications are in near-infrared laser imaging technology for stroke detection in brains of small animals.

Calderón's method on an elliptical domain

One possible application for electrical impedance tomography is in medical imaging where lung and heart function may be monitored. One drawback of current algorithms is that they are implemented for use in a circular domain, but a human thorax is more elliptical than circular. In this paper, a reconstruction algorithm based on the work of Calderón (1980 Seminar on Numerical Analysis and its Applications to Continuum Physics (Rio de Janeiro) pp 65-75) on the inverse conductivity problem is derived for an elliptical domain. It is explained how this reconstruction algorithm uses a transformed Dirichlet-to-Neumann map. Experimental results from an elliptical tank are given to show how correct domain modelling reduces the artefacts produced by this version of Calderón's reconstruction algorithm.

Direct EIT reconstructions of complex admittivities on a chest-shaped domain in 2-D

Electrical impedance tomography (EIT) is a medical imaging technique in which current is applied on electrodes on the surface of the body, the resulting voltage is measured, and an inverse problem is solved to recover the conductivity and/or permittivity in the interior. Images are then formed from the reconstructed conductivity and permittivity distributions. In the 2-D geometry, EIT is clinically useful for chest imaging. In this work, an implementation of a D-bar method for complex admittivities on a general 2-D domain is presented. In particular, reconstructions are computed on a chest-shaped domain for several realistic phantoms including a simulated pneumothorax, hyperinflation, and pleural effusion. The method demonstrates robustness in the presence of noise. Reconstructions from trigonometric and pairwise current injection patterns are included.

A direct D-bar reconstruction algorithm for recovering a complex conductivity in 2-D

A direct reconstruction algorithm for complex conductivities in W^2,∞ (Ω), where Ω is a bounded, simply connected Lipschitz domain in ℝ², is presented. The framework is based on the uniqueness proof by Francini [Inverse Problems 20 2000], but equations relating the Dirichlet-to-Neumann to the scattering transform and the exponentially growing solutions are not present in that work, and are derived here. The algorithm constitutes the first D-bar method for the reconstruction of conductivities and permittivities in two dimensions. Reconstructions of numerically simulated chest phantoms with discontinuities at the organ boundaries are included.

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

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.

Factorization method and its physical justification in frequency-difference electrical impedance tomography

Time-difference electrical impedance tomography (tdEIT) requires two data sets measured at two different times. The difference between them is utilized to produce images of time-dependent changes in a complex conductivity distribution inside the human body. Frequency-difference EIT (fdEIT) was proposed to image frequency-dependent changes of a complex conductivity distribution. It has potential applications in tumor and stroke imaging since it can visualize an anomaly without requiring any time-reference data obtained in the absence of an anomaly. In this paper, we provide a rigorous analysis for the detectability of an anomaly based on a constructive and quantitative physical correlation between a measured fdEIT data set and an anomaly. From this, we propose a new noniterative frequency-difference anomaly detection method called the factorization method (FM) and elaborate its physical justification. To demonstrate its practical applicability, we performed fdEIT phantom imaging experiments using a multifrequency EIT system. Applying the FM to measured frequency-difference boundary voltage data sets, we could quantitatively evaluate indicator functions inside the imaging domain, of which values at each position reveal presence or absence of an anomaly. We found that the FM successfully localizes anomalies inside an imaging domain with a frequency-dependent complex conductivity distribution. We propose the new FM as an anomaly detection algorithm in fdEIT for potential applications in tumor and stroke imaging.

Reconstruction convergence and speed enhancement in electrical impedance tomography for domains with known internal boundaries

An improved approach for electrical impedance tomography (EIT) image reconstruction, based on modifying the forward and inverse solutions, is proposed. In this approach, the EIT forward problem is solved via the finite element method (FEM) using two types of elements. The inverse problem is solved by the modified Newton-Raphson method, whereas the condition number of the Hessian matrix is being monitored. At the early stage of the reconstruction, first-order elements are used, and if the condition number exceeds the allowable limit, the algorithm restarts. Otherwise, if the reconstruction error becomes lower than a predefined threshold, second-order elements are employed in the forward solution in order to preserve the precision of the final results. The latter stage converges in very few iterations. Since the solution speed with the first-order FEM is considerably higher than the second-order FEM, the reconstruction speed improves considerably by this approach, whereas the accuracy of the results is guaranteed by the well-conditioned Hessian matrix. Numerical simulations and experiments are followed by comparisons with other reconstruction methods which demonstrate the reliability and high solution speed of this approach. According to the results, the convergence of the proposed method is significantly improved, and its speed is 2-200 times higher than the previously developed methods with the same level of precision.

Effect of domain shape modeling and measurement errors on the 2-D D-bar method for EIT

The D-bar algorithm based on Nachman's 2-D global uniqueness proof for the inverse conductivity problem (Nachman, 1996) is implemented on a chest-shaped domain. The scattering transform is computed on this chest-shaped domain using trigonometric and adjacent current patterns and the complete electrode model for the forward problem is computed with the finite element method in order to obtain simulated voltage measurements. The robustness and effectiveness of the method is demonstrated on a simulated chest with errors in input currents, output voltages, electrode placement, and domain modeling.