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.

Robust imaging using electrical impedance tomography: review of current tools

Electrical impedance tomography (EIT) is a medical imaging technique with many advantages and great potential for development in the coming years. Currently, some limitations of EIT are related to the ill-posed nature of the problem. These limitations are translated on a practical level by a lack of genericity of the developed tools. In this paper, the main robust data acquisition and processing tools for EIT proposed in the scientific literature are presented. Their relevance and potential to improve the robustness of EIT are analysed, in order to conclude on the feasibility of a robust EIT tool capable of providing resistivity or difference of resistivity mapping in a wide range of applications. In particular, it is shown that certain measurement acquisition tools and algorithms, such as faulty electrode detection algorithm or particular electrode designs, can ensure the quality of the acquisition in many circumstances. Many algorithms, aiming at processing acquired data, are also described and allow to overcome certain difficulties such as an error in the knowledge of the position of the boundaries or the poor conditioning of the inverse problem. They have a strong potential to faithfully reconstruct a quality image in the presence of disturbances such as noise or boundary modelling error.

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.

Calderón's Method with a Spatial Prior for 2-D EIT Imaging of Ventilation and Perfusion

Bedside imaging of ventilation and perfusion is a leading application of 2-D medical electrical impedance tomography (EIT), in which dynamic cross-sectional images of the torso are created by numerically solving the inverse problem of computing the conductivity from voltage measurements arising on electrodes due to currents applied on electrodes on the surface. Methods of reconstruction may be direct or iterative. Calderón’s method is a direct reconstruction method based on complex geometrical optics solutions to Laplace’s equation capable of providing real-time reconstructions in a region of interest. In this paper, the importance of accurate modeling of the electrode location on the body is demonstrated on simulated and experimental data, and a method of including a priori spatial information in dynamic human subject data is presented. The results of accurate electrode modeling and a spatial prior are shown to improve detection of inhomogeneities not included in the prior and to improve the resolution of ventilation and perfusion images in a human subject.

Fusing electrical and elasticity imaging

Electrical and elasticity imaging are promising modalities for a suite of different applications, including medical tomography, non-destructive testing and structural health monitoring. These emerging modalities are capable of providing remote, non-invasive and low-cost opportunities. Unfortunately, both modalities are severely ill-posed nonlinear inverse problems, susceptive to noise and modelling errors. Nevertheless, the ability to incorporate complimentary datasets obtained simultaneously offers mutually beneficial information. By fusing electrical and elastic modalities as a joint problem, we are afforded the possibility to stabilize the inversion process via the utilization of auxiliary information from both modalities as well as joint structural operators. In this study, we will discuss a possible approach to combine electrical and elasticity imaging in a joint reconstruction problem giving rise to novel multi-modality applications for use in both medical and structural engineering. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.

Shape-Driven EIT Reconstruction Using Fourier Representations

Shape-driven approaches have been proposed as an effective strategy for the electrical impedance tomography (EIT) reconstruction problem in recent years. In order to augment the shape-driven approaches, we propose a new method that transforms the shape to be reconstructed as basic primitives directly modeled by using Fourier representations. To allow automatic topological changes between the basic primitives and surrounding objects simultaneously, Boolean operations are employed. The Boolean operations with direct representation of primitives can be utilized for dimensionality and ill-posedness reduction, enabling feasible shape and topology optimization with shape-driven approaches. As a proof of principle, we leverage the proposed method for two dimensional shape reconstruction in EIT with various conductivity distributions. We demonstrate that our method is able to improve EIT reconstructions by enabling accurate shape and topology optimization.

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.

Shape-Driven Difference Electrical Impedance Tomography

This work proposes a novel shape-driven reconstruction approach for difference electrical impedance tomography (EIT). In the proposed approach, the reconstruction problem is formulated as a shape reconstruction problem and solved via an explicit and geometrical methodology, where the geometry of the embedded inclusions is represented by a shape and topology description function (STDF). To incorporate more geometry and prior information directly into the reconstruction and to provide better flexibility in the solution process, the concept of a moving morphable component (MMC) is applied here implying that MMC is treated as the basic building block of the embedded inclusions. Simulations, phantom studies, and in vivo pig data are used to test the proposed approach for the most popular biomedical application of EIT - lung imaging - and the performance is compared with the conventional linear approach. In addition, the modality's robustness is studied in cases where (i) modeling errors are caused by inhomogeneity in the background conductivity, and (ii) uncertainties in the contact impedances and reference state are present. The results of this work indicate that the proposed approach is tolerant to modeling errors and is fairly robust to typical EIT uncertainties, producing greatly improved image quality compared to the conventional linear approach.

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.

Reconstruction of Organ Boundaries With Deep Learning in the D-Bar Method for Electrical Impedance Tomography

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.

B-Spline Level Set Method for Shape Reconstruction in Electrical Impedance Tomography

A B-spline level set (BLS) based method is proposed for shape reconstruction in electrical impedance tomography (EIT). We assume that the conductivity distribution to be reconstructed is piecewise constant, transforming the image reconstruction problem into a shape reconstruction problem. The shape/interface of inclusions is implicitly represented by a level set function (LSF), which is modeled as a continuous parametric function expressed using B-spline functions. Starting from modeling the conductivity distribution with the B-spline based LSF, we show that the shape modeling allows us to compute the solution by restricting the minimization problem to the space spanned by the B-splines. As a consequence, the solution to the minimization problem is obtained in terms of the B-spline coefficients. We illustrate the behavior of this method using simulated as well as water tank data. In addition, robustness studies considering varying initial guesses, differing numbers of control points, and modeling errors caused by inhomogeneity are performed. Both simulation and experimental results show that the BLS-based approach offers clear improvements in preserving the sharp features of the inclusions in comparison to the recently published parametric level set method.

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.

A Moving Morphable Components Based Shape Reconstruction Framework for Electrical Impedance Tomography

This paper presents a new computational framework in electrical impedance tomography (EIT) for shape reconstruction based on the concept of moving morphable components (MMC). In the proposed framework, the shape reconstruction problem is solved in an explicit and geometrical way. Compared with the traditional pixel or shape-based solution framework, the proposed framework can incorporate more geometry and prior information into shape and topology optimization directly and therefore render the solution process more flexibility. It also has the afford potential to substantially reduce the computational burden associated with shape and topology optimization. The effectiveness of the proposed approach is tested with noisy synthetic data and experimental data, which demonstrates the most popular biomedical application of EIT: lung imaging. In addition, robustness studies of the proposed approach considering modeling errors caused by non-homogeneous background, varying initial guesses, differing numbers of candidate shape components, and differing exponent in the shape and topology description function are performed. The simulation and experimental results show that the proposed approach is tolerant to modeling errors and is fairly robust to these parameter choices, offering significant improvements in image quality in comparison to the conventional absolute reconstructions using smoothness prior regularization and total variation regularization.

A Statistical Shape-Constrained Reconstruction Framework for Electrical Impedance Tomography

A statistical shape-constrained reconstruction (SSCR) framework is presented to incorporate the statistical prior information of human lung shapes for lung electrical impedance tomography. The prior information is extracted from 8000 chest-computed tomography scans across 800 patients. The reconstruction framework is implemented with two approaches-a one-step SSCR and an iterative SSCR in lung imaging. The one-step SSCR provides fast and high accurate reconstructions of healthy lungs, whereas the iterative SSCR allows to simultaneously estimate the pre-injured lung and the injury lung part. The approaches are evaluated with the simulated examples of thorax imaging and also with the experimental data from a laboratory setting, with difference imaging considered in both the approaches. It is demonstrated that the accuracy of lung shape reconstruction is significantly improved. In addition, the proposed approaches are proved to be robust against measurement noise, modeling error caused by inaccurately known domain boundary, and the selection of the regularization parameters.

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.

Comparison of Selected Machine Learning Algorithms for Industrial Electrical Tomography

The main goal of this work was to compare the selected machine learning methods with the classic deterministic method in the industrial field of electrical impedance tomography. The research focused on the development and comparison of algorithms and models for the analysis and reconstruction of data using electrical tomography. The novelty was the use of original machine learning algorithms. Their characteristic feature is the use of many separately trained subsystems, each of which generates a single pixel of the output image. Artificial Neural Network (ANN), LARS and Elastic net methods were used to solve the inverse problem. These algorithms have been modified by a corresponding increase in equations (multiply) for electrical impedance tomography using the finite element method grid. The Gauss-Newton method was used as a reference to machine learning methods. The algorithms were trained using learning data obtained through computer simulation based on real models. The results of the experiments showed that in the considered cases the best quality of reconstructions was achieved by ANN. At the same time, ANN was the slowest in terms of both the training process and the speed of image generation. Other machine learning methods were comparable with the deterministic Gauss-Newton method and with each other.

A Parametric Level Set-Based Approach to Difference Imaging in Electrical Impedance Tomography

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.