I Introduction
Electrical Impedance Tomography (EIT) images traditionally display the tissuedependent conductivity distribution of a patient in the plane of the attached measurement electrodes allowing, e.g., visualization of heart and lung function as well as injuries [1, 2, 3, 4, 5, 6]. The resulting images are of highcontrast and data acquisition is done by harmless electrical measurements without the need for contrast agents or ionizing radiation. However, the image recovery process of forming the EIT image from the current/voltage measurement data is a severely illposed nonlinear inverse problem, and thus requires a noiserobust regularization strategy for stability. The ‘Dbar method’, the only proven regularization strategy for the full nonlinear problem [7], provides realtime noiserobust image recovery by using a lowpass filter of the associated nonlinear Fourier data. Unfortunately, this results in images that suffer a loss of sharp features often important in medical imaging applications. In this work, we propose combining Dbar with Deep Learning, specifically with a Convolutional Neural Network, to ‘learn’ and undo the image blurring resulting in realtime sharp EIT images.
EIT reconstructions are typically computed with iterative algorithms that are based on minimizing a penalty functional, such as [8, 9]. These methods perform very well in reconstruction quality due to a flexibility of incorporating prior knowledge, but require careful modeling of the boundary shape in the repeated simulation of the forward problem. Possibilities to overcome the boundary sensitivity are proposed in [10, 11]
, but tend to be computationally demanding. On the other hand, direct (noniterative) reconstruction algorithms do not need the repeated simulation of the forward operator. One such method is known as the Dbar algorithm which is based on a nonlinear Fourier transformation of the measured surface current/voltage data. The method employs a lowpass filtering of this transformed data as a regularization strategy to stabilize the image reconstruction process against noise in the measured data. Consequently, this filtering results in reconstructed images that suffer from a significant loss of sharpness. It has been shown that the direct Dbar method is robust to incorrect or incomplete knowledge of electrode locations as well as errors in boundary shape, see for instance
[12] and the discussion in Section IIB. Iterative methods on the other hand are either very sensitive to the correct forward model or are based on sophisticated modelling to cope with uncertainties in the model, such as unknown electrode locations, boundary shape, or contact impedances [11, 13, 10].Recent advances in the larger field of image reconstruction have demonstrated the power of Deep Learning and Neural Networks for improving low quality or corrupted images. In particular, combining fast direct reconstruction procedures with deep neural networks can provide high quality images with low latency, leading to prospective realtime imaging in many applications. Convolutional Neural Networks (CNN) are especially suitable for postprocessing initial reconstructions that come from algorithms based on, or related to, Fourier transforms, as suggested in [14]. Such initial reconstructions typically suffer from a loss of spatial resolution, due to some sort of lowpass filtering, as well as additional undersampling artefacts. Training a CNN to remove these artefacts to improve the information content of the reconstructed image has been studied for several linear inverse problems in medical imaging, including CT [14, 15], MRI [16], and PAT [17, 18]. Although the EIT problem is nonlinear in nature, the lowpass filtered images from the lowpassed Dbar method naturally fit into this setting.
In this study we formulate a realtime capable reconstruction algorithm that produces high quality sharp absolute EIT images by combining the Dbar algorithm with subsequent processing by a CNN. For this task we utilize an established CNN architecture, known as Unet, adjusted to cope with the typical image structures of Dbar EIT reconstructions. We train the network on simulated training data and directly apply the trained network to experimental data with no training on experimental data itself. This successful transition to experimental data underlines the robustness of the Dbar algorithm and is especially important as the need for good training data is often the bottleneck for the success of such networkbased approaches for other imaging modalities, [14, 19, 18].
This paper is organized as follows. Section II presents a brief review of the mathematical problem of EIT and the Dbar solution method. The deep learning CNN for Dbar, coined ‘Deep Dbar’ is introduced in Section III. The experimental setup as well as simulation of training data are described in Section IV and results presented in Section V. A discussion of the results is given in Section VI and conclusions drawn in Section VII. The reader is encouraged to view the manuscript on a computer screen as details in the image contrast may be masked in printed versions.
Ii Electrical Impedance Tomography and the Dbar reconstruction method
Electrical impedance tomography is a nonlinear inverse problem in which we aim to determine the interior conductivity from currenttovoltage measurements at the boundary. The problem can be formulated as a generalized Laplace equation
(1) 
modeling the electrical potential inside the domain for a given conductivity , with the Neumann boundary condition describing the boundary voltage occurring from the applied meanfree current . The measurement data consists of pairs of current and voltage measurements and is modeled by the currenttovoltage map defined by
This measurement operator is also known as the NeumanntoDirichlet (ND) map, and knowledge of it allows one to predict the resulting voltage for any injected current pattern for . In practice, an approximation to the ND map is formed by applying a basis of current patterns and tracking the responses of the voltages. The Dbar algorithm we use below requires the corresponding DirichlettoNeumann (DN) map, which can be obtained as the inverse of the ND map, , for full (vs. partial) boundary data. In this work we consider the case as the Dbar reconstruction framework is further developed in 2D. However, we expect a natural extension to 3D [20].
Iia Realtime reconstructions using an approximate Dbar method
By the Dbar method, we refer to the regularized Dbar method [7] based on the theoretical proof given in [21]. The approach uses a nonlinear Fourier transform, called a scattering transform, tailormade for the EIT problem which is applied to the measured current/voltage data in the form of the DN map
. That scattering data is then used as input data into a partial differential equation, a
or ‘Dbar’ equation, giving the method its name. Note that the derivative operators and are defined as and , where . The conductivity is then recovered directly from the solution to the Dbar equation.The Dbar approach [21] is to transform the physical conductivity equation into a nonphysical Schrödinger equation, solve that problem instead using the Dbar methods popularized by Beals and Coifman [22], and then transform back to the physical setting. The change of variables and produces the desired Schrödinger equation , where . Provided that is constant in a neighborhood of the boundary, without loss of generality near , the conductivity can be extended from to the entire plane by setting for . Note that this gives the potential compact support in . We make use of special solutions to the Schrödinger equation
(2) 
called Complex Geometrical Optics (CGO) solutions, that have a specific asymptotic behavior for large or , . Note that we associate with via the mapping and thus denotes complex multiplication. The CGO solutions solve a Dbar equation in the nonphysical scattering variable
(3) 
where and is the nonlinear scattering data defined by
(4) 
Note that this scattering data can be thought of as nonlinear Fourier data by the following observation. Replacing the CGO solutions in (4) with the asymptotic behavior yields
and thus the ‘Born’ approximation is essentially a shifted Fourier transform of the potential . A connection to the measurement data can be established via Alessandrini’s identity [23]
In this work we use this ‘Born’ approximation to the scattering data, first presented in [24], as it allows the Dbar method to solve the EIT problem fast enough to be considered ‘realtime’ [25] and is robust against noisy data. The main steps in the algorithm are outlined below:
Step 1: For each , evaluate the approximate scattering data
(5) 
Step 2: For each , solve the Dbar equation (3) using the integral equation
(6) 
and recover the approximate conductivity
(7) 
IiB Robustness of Dbar Methods for EIT
Recent studies [12, 26] suggest that Dbar based reconstruction methods for 2D EIT are robust to incorrect electrode locations and boundary shape. This robustness holds for absolute, as well as timedifference, imaging with both images behaving similarly to incorrect boundary shape and electrode locations. This may be due to the fact that incorrect domain modeling leads to EIT data from a DN map that is only possible for an anisotropic conductivity, even when the true conductivity is isotropic. While the anisotropic conductivity cannot be recovered uniquely, one can recover a unique isotropization, , of the matrixvalued anisotropic conductivity, interpreted as a deformation of the true anisotropic conductivity by isothermal coordinates. In [27, 28], it is proved that the equations in the Dbar reconstruction methods are identical for anisotropic and isotropic EIT data, helping to explain why Dbar methods have still produced quality images even on anisotropic conductivities and imprecisely known boundary shapes. Here we focus on absolute images.
Iii Deep Dbar
The aim of this study is to formulate a realtime reconstruction algorithm for electrical impedance tomography that produces sharp and robust absolute EIT images. To achieve this we combine the Dbar algorithm, described in Section IIA, with a convolutional neural network (CNN). This idea relies on a network architecture known as UNet [29], originally developed for image segmentation. It has been shown for several linear inverse problems [14, 15, 16, 17, 18] that this particular network structure can be modified to successfully remove artefacts in medical image reconstructions. The basic recipe is to use a fast and simple reconstruction algorithm to obtain corrupted images and then train the network to remove those artefacts. A related study for electrical impedance tomography is [30], where the authors used artificial neural networks (ANNs) to postprocess initial reconstructions from one step of a linear GaussNewton algorithm for 3D timedifference EIT imaging. Our approach is fundamentally different as it recovers absolute EIT images.
The network architecture we have chosen relies on the established UNet [29], which consists of a multilevel decomposition and several skip connections to avoid singularities in the training procedure, see Figure 1 for an illustration of our specific architecture. The original purpose of UNet was image segmentation. This is very similar to our application, where the main goal is to identify organ boundaries and deconvolve the reconstruction, hence the output of our network is a sharpened image. Therefore, we believe that the UNet architecture is a suitable choice for the purpose of EIT imaging, since the multilevel structure can deal efficiently with the nonlinearity and sharpening over large image areas. Additionally, as discussed in [31], pooling layers leads to translational invariance, which is important to reduce locational bias in the reconstruction process and detect injuries not present in the training set. As a modification to the original architecture we needed to increase the convolutional filter size to (compared to ), presumably to deal with the nonlinearity of the reconstructions and enforce consistency of the reconstructions. We would like to note, that in contrast to the studies in [14, 15, 16, 17, 18], where the authors learn a residual update to the initial reconstruction, we train the network to produce a single sharpened version of the input.
Iiia Training of the network
Given the true conductivity , we simulate measurement data, as will be described in Section IVA, and reconstruct the approximate conductivity with the Dbar method outlined in IIA. Since the reconstruction step (6) in the Dbar algorithm can be done for any we reconstruct on the square to obtain a square image as input to the network. The resolution is chosen to be . The ground truth is similarly extended to by extending the background conductivity.
Having obtained the training set , we train the Deep Dbar network, denoted by , for the set of network parameters , i.e. the convolutional filters and biases in each convolutional layer. Given the output of the network we seek to minimize the error of network output to phantom, given by the loss
The network is implemented with the Python library TensorFlow and the optimization is performed for 1,000 epochs in batches of 16, with TensorFlow’s implementation of the Adam algorithm and an initial learning rate of
. The training procedure takes only 4 hours on a single Titan XP GPU with 12GB memory. As we will discuss in the following section, we do not need to perform a transfer training to apply the trained Deep Dbar network to experimental data, the training on simulated data proved to be sufficient.Iv Experimental Setup and Computational Notes
We will demonstrate the new Deep Dbar method using experimental data from two different EIT machines: ACT4 [32, 33] from Rensselaer Polytechnic Institute (RPI) as well as KIT4 [34] from the University of Eastern Finland (UEF).
The ACT4 data uses agar (4%) based targets with added graphite (10%) to simulate a heart, two lungs, an aorta, and a spine. All images are shown in DICOM orientation, meaning that the right lung corresponds to the viewer’s left, as if we are looking up through the patient’s feet. Injuries were simulated in the right (DICOM) lung away from the heart by removing a portion of the lung and (1) replacing the missing portion with a piece of agar/graphite with the same conductivity as the heart to simulate an injury such as a pleural effusion, (2) placing three plastic tubes in the missing region to simulate an area of very low conductivity such as a pneumothorax, and (3) replacing the missing portion with three metal tubes. The experiments are shown in Figure 2. The approximate conductivities of the targets are displayed in Table I. The admittivity spectrum of the agar/graphite targets were measured on testcells with Impedimed’s SFB7 bioimpedance meter^{1}^{1}1https://www.impedimed.com/products/sfb7forbodycomposition/. Note that the ACT4 system applies voltages and measures currents rather than viceversa. In these experiments, trigonometric voltage patterns of maximum amplitude 0.5V (and frequency 3.3kHz) were applied on a circular tank (radius 15cm), with 32 electrodes (width 2.5cm), filled with saline (0.3 S/m) to a height of 2.25cm.
Measured Values  Simulated Values  
(S/m)  Ranges (S/m)  
Heart/Aorta  0.67781  [0.5, 0.8] 
Lungs/Spine  0.056714  [0.01, 0.2] 
Saline Background  0.3  [0.29, 0.31] 
Injury 1: Agar/Graphite  0.67781  [0.01, 1.5] 
Injury 2: Plastic Tubes  0  [0.01, 1.5] 
Injury 3: Copper Tubes  infinite  [0.01, 1.5] 
The KIT4 data was taken on a circular tank of radius 14cm with 16 electrodes of width 2.5cm and tap water with conductivity 0.03 S/m filled to a height of 7cm. Conductive (metal) and resistive (plastic) targets were placed in the tank, as shown in Figure 3, and adjacent current patterns with amplitude 2mA were applied at 1kHz. We remark that while this data may not satisfy safety standards for human imaging, it is included for illustrative purposes and potential industrial applications.
Iva Simulation of 2D EIT data
The boundary conditions of (1) assume a continuum model for the boundary measurements, completely ignoring the discrete positioning of the electrodes. When simulating the training data, we use a modified version of the continuum model, called the continuum electrode model introduced in [35], which was developed to simulate realistic electrode data in a continuum setting. In essence, the continuum current/voltage traces are optimally projected onto subsets of the boundary corresponding to the electrode locations. The training could be done with a more complicated electrode model, such as the Complete Electrode Model (CEM) [36], however our simplified continuum electrode model proved sufficient for this proof of concept study.
We aim to represent the ND map as matrix approximation with respect to an orthonormal basis on the boundary. Let be an even number of electrodes, then the basis functions are chosen for as
The ACT4 system uses electrodes and the KIT4 system uses . The measured voltages are then projected to a continuum trace , see [35, 37], and we obtain the ND matrix by evaluating inner products in as follows
(8) 
The matrix approximation of the DN map, , is then formed by inverting the ND matrix, i.e. . If the maximal radius of the domain is not 1, the DN matrix can be scaled by to correspond to the data that would be obtained if the radius were 1. Similarly, if near , the DN matrix is scaled by to produce the DN matrix that would correspond to
near the boundary. If an estimate for
is not available, the best constant conductivity approximation to the data can be formed as described in [24]. The scaling is undone at the end of the Dbar algorithm by multiplying the conductivity by . The matrix approximation to the DN map , required to evaluate the scattering data via (5), is simulated using the constant conductivity .IvB Simulation of Training Data
Training data for the neural network was created using solely simulated data: one group for the ACT4 data and another group for the KIT4 data.
The ACT4 training data was created as follows. Using the ‘Healthy’ image, shown in Figure 2 (top left), approximate organ boundaries were extracted by clicking around the targets in the image for the heart, aorta, left lung, right lung, and spine (Fig. 4, top right). Random numbers were generated to decide whether each individual target was included, heart (95%), aorta (95%), left lung (90%), right lung (90%), spine (100%). If a given target was included, white Gaussian noise (25db) was added to the approximate boundary points of the target using the awgn command in Matlab to create ‘noisy’ boundary locations. Figure 4
(bottom) shows the effect of the white noise on the boundary locations. Noise was added to each target/organ independently. Conductivities were assigned for each included target by generating a random number from a uniform distribution in the ranges shown in Table
I, last column. Elementary injuries were simulated by generating a horizontal dividing line in the lung and assigning randomly generated values in each of the two portions of the divided lung from the uniform distribution of values in [0.01, 1.5], see Fig. 4 bottom right. Each lung had an independent chance of such an injury (30%). More complex injuries could be simulated but are outside the scope of this study. A total of 4,096 simulations were performed for the ACT4 training.After each conductivity phantom was constructed, the mathematical forward problem (1) was solved to recover the corresponding theoretical boundary voltages and currents using a FEM mesh with 65,536 triangular elements using the continuum electrode model described in Section IVA
. Relative white noise with variance of
was added to the measured voltages. The resulting simulated voltages/currents were used to solve the inverse problem using the Dbar method described in Section IIA with a lowpass filtering radius of in the scattering domain using the procedure outlined in [38] and uniformly spaced and grids on with stepsize , and with stepsize , respectively. A nonuniform cutoff threshold was enforced on the scattering data for frequencies such that if either or exceeds 24. Then, the 4,096 pairs of data in the form of ‘Truth’ and ‘Lowpass Dbar Reconstruction’ were used to train the convolutional neural network described in Section III.Training data for the KIT4 experiments was simulated in a similar manner. In this case, one to three circular inclusions were simulated, with varying radii drawn from the uniform distribution on [0.2, 0.4], with center in [0, 0.6], and an angle in [0, 2
]. Inclusions were not allowed to overlap and each inclusion had an equal probability of being ‘conductive’ or ‘resistive’, and values were assigned from the Uniform distributions [0.05,0.12] and [0.005,0.015] in S/m, respectively. The conductivity of the background was drawn from [0.027, 0.033] S/m. The conductivity ranges for the inclusions were determined from initial higherpass (larger filtering radius in the scattering domain) Dbar reconstructions of the KIT4 current/voltage data. While we note that the infinite (metal) conductors should have much higher conductivities, this range was observed in the initial testing and proved sufficient for classifying objects as conductors/resistors in our study. Note that such infinite conductors/resistors violate the theory of Dbar which requires inclusions to have nonnegative conductivities bounded away from zero and infinity. Nevertheless, the method provides useful conductor/resistor information. The same
and grids were used in the Dbar reconstructions of the KIT4 example as in the ACT4 example. However, we reduced the nonuniform cutoff threshold of the scattering data from 24 to 8 to reduce oscillatory artefacts that can result from noise in the scattering data for higher frequencies. Figure 5 shows sample phantoms used in the training data for the KIT4 example. A total of 4,096 simulations were performed and pairs of ‘Truth’ vs. ‘Lowpass Dbar Reconstruction’ used to train the network.IvC Computational Notes for Dbar Reconstructions from Experimental EIT Data
The Dbar reconstructions from the experimental ACT4 and KIT4 data were computed in the same manner as the simulated data case described above in Section IVB with the exception of the formation of the DN matrices and , which now come from discrete vs. continuous measurements. For convenience, for both the ACT4 and KIT4 data, we synthesized the current/voltage measurements that would have occurred if orthonormal trigonometric current patterns had been applied, by using a change of basis. Define to be the value of the th normalized trigonometric current pattern on the th electrode, following Isaacson et al. [24],
for . This ensures that the matrix of current patterns are orthonormal allowing the solution method of DeAngelo and Mueller [39]. Alternative methods such as GramSchmidt based approaches could also be used to produce a matrix of orthonormal currents. The corresponding voltages are formed using a change of basis matrix relating the physical currents and the normalized trig currents. As ACT4 applies voltages and measures currents, we must enforce orthonormality of the currents. Similarly, the KIT4 data used adjacent current patterns which are not orthogonal and must be converted.
Using the approach introduced in [24], the entry of the ND matrix was formed as the discrete inner product
(9) 
where denotes the th normalized current pattern, the voltage resulting from applying the th current pattern (normalized such that the voltages sum to zero and are scaled by the norms of the applied current patterns), and denotes the area of the th electrode. This formula holds for a system with electrodes where linearly independent current patterns have been applied.
The discrete DN matrix was then formed by and scaled by as described above. As the scattering data (5) requires the difference in DN matrices , the discrete DN matrices must be formed for both the ACT4 and KIT4 system with and electrodes, respectively. To this end, the conductivity equation in (1) was solved, using , with boundary conditions given by the Complete Electrode Model [36] on a FEM mesh with triangular elements (ACT4: 4,149; KIT4: 3,493) simulating injected trigonometric current patterns of amplitude 1mA and nonoptimized constant contact impedances of 0.00024 mm.
The scattering data (5) was evaluated as a simple Simpson’s rule approximation, following [39],
where
is the vector of the positions of the centers of the electrodes,
denotes the traditional matrix transpose, andis the expansion of the asymptotic behavior at the center of the th electrode, , in the basis of normalized applied current patterns . Then, the dbar equation (6) can be solved using Fast Fourier Transforms using Vainikko’s method [40]. The interested reader is referred to [38] for further details.
V Results
We now demonstrate the effect of the Deep Dbar method on simulated, as well as experimental, data for absolute EIT imaging.
Va Reconstructions from Simulated Data
We begin with purely simulated data for the ACT4 and KIT4 examples. For the ACT4 setting, we consider three scenarios, as shown in Figure 6: one consistent with the training data but not used in the training (top), and two examples violating the training data  one with three horizontally divided regions in the left lung (middle) and the final with a vertical division in the left lung (bottom). Note that the ‘Lowpass Dbar’ and ‘Deep Dbar’ images are shown on the same scale for visual comparison. The complete ‘input’ and ‘output’ of the CNN are on the unit square . Reconstructions here are clipped to the disc for visualization purposes only. Next, Figure 7 shows the results of the new algorithm on simulated data for the KIT4 example. Three scenarios consistent with the training data, but not used in the training, are presented.
Structural Similarity Indices (SSIMs) computed for the simulated ACT4 and KIT4 examples are shown in Figures 8 and 9, respectively. Additionally, we evaluated the minimized loss by computing the mean relative error for a test set of 16 samples drawn from the same distribution as the training data. The For the ACT4 simulations we improved from to and for the KIT4 test data from to relative error.
VB Reconstructions from Experimental Data
Next, we proceed to reconstructions from experimental data. Figure 10 depicts the results of the Deep Dbar approach on four experiments with ACT4 data: Healthy and Injuries 13 as shown in Figure 2. The black dots represent the approximate boundaries of the ‘healthy’ organs, extracted from the photograph. SSIMs (Figure 11) were computed for the experimental reconstructions with the exception of Injury 3, which has the infinite conductors (copper tubes). The SSIM comparisons used approximate ‘truth’ images formed by assigning the measured conductivity values (Table I) in the respective regions.
Lastly, Figure 12 shows results of the method on the four KIT4 scenarios shown in Figure 3. The overlaid black dots depict the approximate ‘true’ locations of the targets as extracted from their corresponding photographs. No SSIMs were computed here since the objects are infinite conductors and resistors. For a comparison to an iterative method with a total variation prior [9], performed on the same KIT4 data, we refer the reader to the documentation [34].
Vi Discussion
The reconstructions shown in Figures 6, 7, 10, and 12 demonstrate that Deep Dbar provides superior reconstructions giving both visual and quantitative improvements. In particular, the SSIMs (Figs. 8, 9, and 11) show significant SSIM increases for the Deep Dbar vs. Lowpass Dbar. Note that for the SSIM computation for ACT4 Injury 2 (plastic tubes), the ‘truth’ image was unrealistically set to zero in the lower portion of the right lung, even though the tubes do not entirely fill that region.
We remind the reader that no experimental (truth, reconstruction) pairs were used in training the network and no adaptation to the experimental system was necessary, apart from the number of electrodes in the system. The training was done purely with simulated data. In most applications, either a transfer training [18] or training with a golden standard from the same system must be performed. This demonstrates the robustness of our approach. Additionally, we expect further improvements in the ACT4 reconstructions if more complicated injuries are included in the training and remind the reader that the Lowpass Dbar and Deep Dbar reconstructions are shown on the same scale, which does mask the true dynamic range of the individual images.
We review additional simplifications used in our process: 1) we used the continuum electrode model for the boundary conditions in the training data, 2) the FEM solver used to form for the ACT4 and KIT4 experimental data examples was not finely tuned to either EIT device (which is required for iterative minimizationbased methods), and 3) the Dbar solver was not optimized for the respective ACT4/KIT4 data. Rather it was used merely to provide the lowpass reconstructions used as inputs in the CNN. These simplifications were used to demonstrate the robustness of the approach to both noise in the data and tolerance to modeling errors at multiple stages of the reconstruction process.
Initial experiments performed with the original UNet architecture, i.e. convolutional filters of size , did not perform satisfactorily leading us to increase the filter size in this study to . This of course leads to an increase in parameters from to
resulting in longer training times. No batch normalization was needed in our training processes and training times were only 4 hours per network, due to the rather small image size.
The evaluation of the CNN is highly efficient on a GPU and took on average ms for a single sample, hence we expect Deep Dbar to be realtime capable. This can be done by combining the Dbar reconstruction, as outlined in [25], with the application of the CNN in a unified framework to reduce overhead due to data transmission.
Via Generalization
An important aspect for medical imaging is the robustness and consistency of reconstructions. The successful transition to experimental data suggests that the proposed Deep Dbar method is robust enough for translational imaging. Furthermore, Figures 6 and 10 (ACT4) illustrate that the network can handle reconstructions of phantoms that do not conform with the training data. However, while we were able to localize the inclusions in Figure 12 (Phantom 2.2, KIT4), the sharp angular boundaries of the triangular target were not recovered when using only circular inclusion training data. Our initial testing suggests that this can be improved upon by including significant training on triangular inclusions. Challenges recovering triangular shapes have also been observed in [41]. In terms of image quality, our Deep Dbar approach is comparable to results from (slower) iterative methods on similar data from the KIT4 system, see [34, 9, 41]. Additionally, as the ACT4 injuries we simulated were elementary (only using a horizontal dividing line rather than the true diagonal cut and incomplete regional replacements), we expect that the reconstructions may improve further if more complex injuries were introduced. For human targets, a larger database of training data could be employed and built from an anatomical atlas or collection of CT/MR scans both including and not including abnormalities/injuries.
Crucial for the success of the postprocessing network is consistency in the input reconstructions. In order to improve flexibility of the network, one can train the network on reconstructions from scattering data with varying cutoff radii. This allows the user to decide on the quality of the measured data at hand and adjust the cutoff radii as needed for the input reconstruction to the network. First tests have shown that this procedure indeed improves consistency and stability of the reconstructions as illustrated in Figure 13, where we have trained the ACT4 network with varying cutoff radii . While the SSIMs remained consistent, the localization and recovered conductivity of the injury did improve with new variable radii network.
While we chose, in this study, to match Dbar and CNNs due to the robustness of Dbar for absolute and timedifference imaging and the convolved nature of the Dbar reconstructions, alternative reconstruction methods for the input images could also be used.
Vii Conclusions
The Dbar method for 2D EIT provides reliable reconstructions of the conductivity but suffers from a blurring due to a lowpass filtering of the scattering data. Sharp improvements in absolute EIT image quality can be achieved by coupling the Dbar reconstruction method with a convolutional neural network. We demonstrated that a CNN can effectively learn the deblurring using only simulated data and still transition to experimental data without including any experimental data in the training itself. As the training can be done offline ahead of time, and the Dbar method provides realtime conductivity reconstructions [25], the postprocessing step by the trained CNN adds minimal time to the overall image recovery process, due to the highly efficient evaluation on a GPU. While this work is shown in 2D, we expect the approach to extend to 3D once the Dbar computational framework has been further developed.
Acknowledgments
We gratefully acknowledge the support of NVIDIA Corporation with the donation of the Titan Xp GPU used for this research. The authors would like to thank the Isaac Newton Institute for Mathematical Sciences, Cambridge, for support and hospitality during the programme ‘Variational methods and effective algorithms for imaging and vision’ where work on this paper was undertaken, EPSRC grant no EP/K032208/1. AH acknowledges support from EPSRC project EP/M020533/1, ‘Medical image computing for nextgeneration healthcare technology’. We additionally thank the EIT group at RPI ^{2}^{2}2https://www.ecse.rpi.edu/homepages/saulnier/eit/eit.html for their assistance and for providing the ACT4 tank data.
References
 [1] Gilda Cinnella, Salvatore Grasso, Pasquale Raimondo, Davide D’Antini, Lucia Mirabella, Michela Rauseo, and Michele Dambrosio. Physiological effects of the open lung approach in patients with early, mild, diffuse acute respiratory distress syndromean electrical impedance tomography study. The Journal of the American Society of Anesthesiologists, 123(5):1113–1121, 2015.
 [2] C.A. Grant, T. Pham, J. Hough, T. Riedel, C. Stocker, and A. Schibler. Measurement of ventilation and cardiac related impedance changes with electrical impedance tomography. Critical Care, 15(1):R37, 2011.
 [3] Christian Karagiannidis, Andreas D. Waldmann, Carlos Ferrando Ortolá, Manuel Muñoz Martinez, Anxela Vidal, Arnoldo Santos, Peter L. Róka, Manuel Perez Márquez, Stephan H. Bohm, and Fernando SuarezSpimann. Positiondependent distribution of ventilation measured with electrical impedance tomography. European Respiratory Journal, 46(suppl 59), 2015.
 [4] Antonio Pesenti, Guido Musch, Daniel Lichtenstein, Francesco Mojoli, Marcelo B. P. Amato, Gilda Cinnella, Luciano Gattinoni, and Michael Quintel. Imaging in acute respiratory distress syndrome. Intensive Care Medicine, 42(5):686–698, 2016.
 [5] H. Reinius, J. B. Borges, F. Fredén, L. Jideus, E. D. L. B. Camargo, M. B. P. Amato, G. Hedenstierna, Larsson A., and F. Lennmyr. Realtime ventilation and perfusion distributions by electrical impedance tomography during onelung ventilation with capnothorax. Acta Anaesthesiologica Scandinavica, 59(3):354–368, 2015.
 [6] A. Schlibler, T.M.Y. Pham, A.A. Moray, and C. Stocker. Ventilation and cardiac related impedance changes in children undergoing corrective open heart surgery. Physiological Measurement, 34:1319–1327, 2013.
 [7] K. Knudsen, M. Lassas, J.L. Mueller, and S. Siltanen. Regularized Dbar method for the inverse conductivity problem. Inverse Problems and Imaging, 3(4):599–624, 2009.
 [8] Zhou Zhou, Gustavo Sato dos Santos, Thomas Dowrick, James Avery, Zhaolin Sun, Hui Xu, and David S Holder. Comparison of total variation algorithms for electrical impedance tomography. Physiological measurement, 36(6):1193, 2015.
 [9] Gerardo González, Ville Kolehmainen, and Aku Seppänen. Isotropic and anisotropic total variation regularization in electrical impedance tomography. Computers & Mathematics with Applications, 2017.
 [10] Jérémi Dardé, N Hyvönen, A Seppänen, and Stratos Staboulis. Simultaneous reconstruction of outer boundary shape and admittivity distribution in electrical impedance tomography. SIAM Journal on Imaging Sciences, 6(1):176–198, 2013.
 [11] A. Nissinen, V. Kolehmainen, and J. P. Kaipio. Compensation of modelling errors due to unknown domain boundary in electrical impedance tomography. IEEE Transaction on Medical Imaging, 30(2):231–242, 2011.
 [12] E. K. Murphy and J. L. Mueller. Effect of domainshape modeling and measurement errors on the 2d Dbar method for electrical impedance tomography. IEEE Transactions on Medical Imaging, 28(10):1576–1584, 2009.
 [13] V. Kolehmainen, M. Lassas, and P. Ola. Electrical impedance tomography problem with inaccurately known boundary and contact impedances. Medical Imaging, IEEE Transactions on, 27(10):1404 –1414, oct. 2008.
 [14] Kyong Hwan Jin, Michael T McCann, Emmanuel Froustey, and Michael Unser. Deep convolutional neural network for inverse problems in imaging. IEEE Transactions on Image Processing, 26(9):4509–4522, 2017.
 [15] Eunhee Kang, Junhong Min, and Jong Chul Ye. A deep convolutional neural network using directional wavelets for lowdose xray ct reconstruction. Medical Physics, 44(10), 2017.
 [16] Christopher M Sandino, Neerav Dixit, Joseph Y Cheng, and Shreyas S Vasanawala. Deep convolutional neural networks for accelerated dynamic magnetic resonance imaging. In NIPS 2017, Medical Imaging Meets NIPS Workshop., 2017.
 [17] Stephan Antholzer, Markus Haltmeier, and Johannes Schwab. Deep learning for photoacoustic tomography from sparse data. arXiv preprint arXiv:1704.04587, 2017.
 [18] A. Hauptmann, F. Lucka, M. Betcke, N. Huynh, J. Adler, B. Cox, P. Beard, S. Ourselin, and S. Arridge. Model based learning for accelerated, limitedview 3d photoacoustic tomography. IEEE Transactions on Medical Imaging, 2018.
 [19] Jonas Adler and Ozan Öktem. Solving illposed inverse problems using iterative deep neural networks. Inverse Problems, 33(12):124007, 2017.
 [20] Fabrice Delbary, Per Christian Hansen, and Kim Knudsen. Electrical impedance tomography: 3d reconstructions using scattering transforms. Applicable Analysis, 91(4):737–755, 2012.
 [21] A. I. Nachman. Global uniqueness for a twodimensional inverse boundary value problem. Annals of Mathematics, 143:71–96, 1996.
 [22] Richard Beals and Ronald R. Coifman. Multidimensional inverse scatterings and nonlinear partial differential equations. In Pseudodifferential operators and applications (Notre Dame, Ind., 1984), pages 45–70. Amer. Math. Soc., Providence, RI, 1985.
 [23] G. Alessandrini. Stable determination of conductivity by boundary measurements. Applicable Analysis, 27:153–172, 1988.
 [24] D. Isaacson, J. L. Mueller, J. C. Newell, and S. Siltanen. Reconstructions of chest phantoms by the Dbar method for electrical impedance tomography. IEEE Transactions on Medical Imaging, 23:821–828, 2004.
 [25] Melody Dodd and Jennifer L Mueller. A realtime Dbar algorithm for 2D electrical impedance tomography data. Inverse problems and imaging, 8(4):1013–1031, 2014.
 [26] S. J. Hamilton, J. L. Mueller, and T. R. Santos. Robust computation of 2d EIT absolute images with dbar methods. (under review, arXiv preprint: 1712.00379).
 [27] Gennadi Henkin and Matteo Santacesaria. On an inverse problem for anisotropic conductivity in the plane. Inverse Problems, 26(9):095011, 2010.
 [28] S. J. Hamilton, M. Lassas, and S. Siltanen. A Direct Reconstruction Method for Anisotropic Electrical Impedance Tomography. Inverse Problems, 30:(075007), 2014.
 [29] Olaf Ronneberger, Philipp Fischer, and Thomas Brox. Unet: Convolutional networks for biomedical image segmentation. In International Conference on Medical Image Computing and ComputerAssisted Intervention, pages 234–241. Springer, 2015.
 [30] Sébastien Martin and Charles TM Choi. A postprocessing method for threedimensional electrical impedance tomography. Scientific Reports, 7, 2017.

[31]
Thomas Wiatowski and Helmut Bölcskei.
A mathematical theory of deep convolutional neural networks for feature extraction.
IEEE Transactions on Information Theory, 2017.  [32] N. Liu, G.J. Saulnier, J.C. Newell, D. Isaacson, and TJ. Kao. Act4: a highprecision, multifrequency electrical impedance tomograph. Presented at 6th Conference on Biomedical Applications of Electrical Impedance Tomography, June 2005. London, U.K.
 [33] G. Saulnier, N. Liu, C. P. Tamma, H. Xia, T.J. Kao, J. Newell, and D. Isaacson. An electrical impedance spectroscopy system for breast cancer detection. Proc. 29th Int. Conf. IEEE Eng. Med. Biol. Soc, (1):4154–4157, Aug. 2007.
 [34] Andreas Hauptmann, Ville Kolehmainen, Nguyet Minh Mach, Tuomo Savolainen, Aku Seppänen, and Samuli Siltanen. Open 2d Electrical Impedance Tomography data archive. arXiv:1704.01178, 2017.
 [35] Andreas Hauptmann. Approximation of fullboundary data from partialboundary electrode measurements. Inverse Problems (accepted), arXiv preprint arXiv:1703.05550, 2017.
 [36] Erkki Somersalo, Margaret Cheney, and David Isaacson. Existence and uniqueness for electrode models for electric current computed tomography. SIAM Journal on Applied Mathematics, 52(4):1023–1040, 1992.
 [37] N. Hyvönen. Approximating idealized boundary data of electric impedance tomography by electrode measurements. Mathematical Models and Methods in Applied Sciences, 19(07):1185–1202, 2009.
 [38] J. Mueller and S. Siltanen. Linear and Nonlinear Inverse Problems with Practical Applications, volume 10 of Computational Science and Engineering. SIAM, 2012.
 [39] M. DeAngelo and J. L. Mueller. 2d Dbar reconstructions of human chest and tank data using an improved approximation to the scattering transform. Physiological Measurement, 31:221–232, 2010.
 [40] Gennadi Vainikko. Fast solvers of the LippmannSchwinger equation. In Direct and inverse problems of mathematical physics (Newark, DE, 1997), volume 5 of Int. Soc. Anal. Appl. Comput., pages 423–440. Kluwer Acad. Publ., Dordrecht, 2000.
 [41] D. Liu, A. K. Khambampati, and J. Du. A parametric level set method for electrical impedance tomography. IEEE Transactions on Medical Imaging, 37(2):451–460, Feb 2018.