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 study, 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 are presented on experimental EIT data from the ACT4 and KIT4 EIT systems.


I. INTRODUCTION
E LECTRICAL IMPEDANCE TOMOGRAPHY (EIT) im- ages traditionally display the tissue-dependent 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]- [6].The resulting images are of high-contrast 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 ill-posed nonlinear inverse problem, and thus requires a noise-robust regularization strategy for stability.The 'D-bar method', the only proven regularization strategy for the full nonlinear problem [7], provides real-time noise-robust image recovery by using a low-pass 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 real-time 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 (non-iterative) reconstruction algorithms do not need the simulation of the forward operator.One such method is known as the D-bar algorithm which is based on a nonlinear Fourier transformation of the measured surface current/voltage data.The method employs a low-pass filtering of this transformed data as a regularization strategy to stabilize the image reconstruction process against noise in the measured data.This filtering results in reconstructed images that suffer from a significant loss of sharpness, but unlike other modalities (such as CT) the D-bar method does not suffer from undersampling artefacts.Furthermore, unlike iterative methods, the direct D-bar method is robust to incorrect or incomplete knowledge of electrode locations as well as errors in boundary shape [12], therefore removing the need for meticulously specific forward models that make or break most EIT methods for absolute EIT imaging.
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 real-time imaging in many applications.Convolutional Neural Networks (CNN) are especially suitable for post-processing initial reconstructions that come from algorithms based on, or related to, Fourier transforms, as suggested in [13].Such initial reconstructions typically suffer from a loss of spatial resolution, due to some sort of low-pass 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 [13], [14], MRI [15], and PAT [16], [17].Although the EIT problem is nonlinear in nature, the low-pass filtered images from the low-passed D-bar method naturally fit into this setting.
In this study we formulate a real-time capable reconstruction algorithm that produces high quality sharp EIT images by combining the D-bar algorithm with subsequent processing by a CNN.For this task we utilize an established CNN architecture, known as U-net, adjusted to cope with the typical image structures of D-bar 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 D-bar algorithm and is especially important as the need for good training data is often the bottleneck for the success of such network-based approaches for other imaging modalities, [13], [17], [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 D-bar, coined 'Deep D-bar' 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 D-BAR RECONSTRUCTION METHOD
Electrical impedance tomography is a nonlinear inverse problem in which we aim to determine the interior conductivity from current-to-voltage measurements at the boundary.The problem can be formulated as a generalized Laplace equation modeling the electrical potential u inside the domain Ω ⊂ R n for a given conductivity σ, with the Neumann boundary condition describing the boundary voltage occurring from the applied mean-free current ϕ.The measurement data consists of pairs of current and voltage measurements and is modeled by the current-to-voltage map R σ defined by This measurement operator is also known as the Neumannto-Dirichlet (ND) map, and knowledge of it allows one to predict the resulting voltage for any injected current pattern for n = 2, 3.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 D-bar algorithm we use below requires the corresponding Dirichlet-to-Neumann (DN) map, which can be obtained as the inverse of the ND map, Λ σ = (R σ ) −1 , for full (vs.partial) boundary data.In this work we consider the n = 2 case as the D-bar reconstruction framework is further developed in 2D.However, we expect a natural extension to 3D [19].

A. Real-time reconstructions using an approximate D-bar method
By the D-bar method, we refer to the regularized D-bar method [7] based on the theoretical proof given in [20].The approach uses a nonlinear Fourier transform, called a scattering transform, tailor-made 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 ∂ k or 'D-bar' equation, giving the method its name.Note that the derivative operators ∂ z and ∂ z are defined as The conductivity σ is then recovered directly from the solution to the D-bar equation.The D-bar approach [20] is to transform the physical conductivity equation ∇ • σ∇u = 0 into a nonphysical Schrödinger equation, solve that problem instead using the D-bar methods popularized by Beals and Coifman [21], and then transform back to the physical setting.The change of variables ũ = σ 1/2 u and q(z) = σ −1/2 (z)∆σ 1/2 (z) produces the desired Schrödinger equation [−∆ + q(z)]ũ(z) = 0, where z ∈ Ω. Provided that σ(z) is constant in a neighborhood of the boundary, without loss of generality σ = 1 near ∂Ω, the conductivity can be extended from Ω to the entire plane by setting σ(z) ≡ 1 for z ∈ C \Ω.Note that this gives the potential q(z) compact support in Ω.We make use of special solutions ψ(z, k) to the Schrödinger equation where e(z, k) := exp{i(kz + kz)} and t(k) is the nonlinear scattering data defined by Note that this scattering data t can be thought of as nonlinear Fourier data by the following observation.Replacing the CGO solutions µ(z, k) in (4) with the asymptotic behavior 1 yields and thus the 'Born' approximation t exp is essentially a shifted Fourier transform of the potential q.A connection to the measurement data Λ σ can be established via Alessandrini's identity [22] In this work we use this 'Born' approximation t exp to the scattering data, first presented in [23], as it allows the D-bar method to solve the EIT problem fast enough to be considered 'real-time' [24] and is robust against noisy data.The main steps in the algorithm are outlined below: The numbers on top of the blue bars denote the channels for each layer.The resolution for each multilevel decomposition is shown in gray on the left.Each convolutional layer is equipped with a Rectified Linear Unit as nonlinearity, given by ReLU(x) = max(0, x).
Step 1: For each k ∈ C \{0}, evaluate the approximate scattering data Step 2: For each z ∈ Ω, solve the D-bar equation (3) using the integral equation and recover the approximate conductivity

III. DEEP D-BAR
The aim of this study is to formulate a real-time reconstruction algorithm for electrical impedance tomography that produces sharp and robust images.To achieve this we combine the D-bar algorithm, described in Section II-A, with a convolutional neural network (CNN).This idea relies on a network architecture known as U-Net [25], originally developed for image segmentation.It has been shown for several linear inverse problems [13]- [17] 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 [26], where the authors used artificial neural networks to postprocess initial reconstructions from one step of a linear Gauss-Newton algorithm.This approach is fundamentally different, since it requires the simulation of a forward model.
The network architecture we have chosen relies on the established U-Net [25], 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 output of our network is the sharpened image, since the main goal is to identify organ boundaries and deconvolve the reconstruction, similar to the original purpose of U-Net for image segmentation.This is in contrast to the studies in [13]- [17], where the authors learn a residual update to the initial reconstruction.As an additional modification to the original architecture we needed to increase the convolutional filter size to 5 × 5 (compared to 3 × 3), presumably to deal with the nonlinearity of the reconstructions and enforce consistency of the reconstructions.

A. Training of the network
Given the true conductivity σ, we simulate measurement data, as will be described in Section IV-A, and reconstruct the approximate conductivity σ exp with the D-bar method outlined in II-A.Since the reconstruction step (6) in the Dbar algorithm can be done for any z ∈ R 2 we reconstruct σ exp on the square [−1, 1] 2 to obtain a square image as input to the network.The resolution is chosen to be 64 × 64.The ground truth σ is similarly extended to [−1, 1] 2 by extending the background conductivity.
Having obtained the training set {σ i , σ exp i } i , we train the Deep D-bar network, denoted by D θ , for the set of network parameters θ, i.e. the convolutional filters and biases in each convolutional layer.Given the output of the network σ = D θ (σ exp ) we seek to minimize the 2 -error of network output to phantom, given by the loss The network is implemented with the Python library Ten-sorFlow and the optimization is performed for 60 epochs in batches of 16, with TensorFlow's implementation of the Adam algorithm and an initial learning rate of 10 −4 .The training procedure takes 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 D-bar 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 D-bar method using experimental data from two different EIT machines: ACT4 [27], [28] from Rensselaer Polytechnic Institute (RPI) as well as KIT4 [29] 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 SFB-7 bioimpedance meter 1 .Note that the ACT4 system applies voltages and measures currents rather than vice-versa.In these experiments, trigonometric voltage patterns of maximum amplitude 0.5V (and frequency 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.When reconstructing the conductivity via the D-bar method described in Section II-A, a change of basis was performed on the measured ACT4 currents and voltages to synthesize the voltages that would have occurred if 1mA maximum amplitude trigonometric current patterns had been applied.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.

A. Simulation of 2D EIT data
The boundary conditions of ( 1) assume a continuum model for the boundary measurements, completely ignoring 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 [30], 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) [31], however our simplified continuum electrode model proved sufficient for this proof of concept study.
We aim to represent the ND map as matrix approximation R σ with respect to an orthonormal basis on the boundary.Let L be an even number of electrodes, then the basis functions are chosen for The ACT4 system uses L = 32 electrodes and the KIT4 system uses L = 16.The measured voltages are then projected to a continuum trace g n , see [30], [32], and we obtain the ND matrix R σ by evaluating inner products in L 2 (∂Ω) as follows The matrix approximation of the DN map, L σ , is then formed by inverting the ND matrix, i.e.L σ = (R σ ) −1 .If the radius r of the domain is not 1, the DN matrix can be scaled by r to correspond to the data that would be obtained if the radius were 1.Similarly, if σ = σ 0 = 1 near ∂Ω, the DN matrix is scaled by 1 σ0 to produce the DN matrix that would correspond to σ = 1 near the boundary.If an estimate for σ 0 is not available, the best constant conductivity approximation to the data can be formed as described in [23].The scaling is undone at the end of the D-bar algorithm by multiplying the conductivity by σ 0 .

B. 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 'HLSA 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. 2, 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 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.A total of 2048 such simulations were made.
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].Each lung had an independent chance of such an injury (60%).More complex injuries could be simulated but are outside the scope of this study.An additional 2048 of such simulations were performed for this approach and combined with the 2048 above giving a total of 4096 phantoms 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 IV-A.Relative white noise with variance of 10 −4 was added to the measured voltages.The resulting simulated voltages/currents were used to solve the inverse problem using the D-bar method described in Section II-A with a low-pass filtering radius of R = 4.5 in the scattering domain using the procedure outlined in [33]  Training data for the KIT4 experiments was simulated in a similar manner.In this case, one to three circular inclusions of varying radii drawn from the uniform distribution on [0.2, 0.4] were simulated.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 drawn from [0.027, 0.033] S/m.The conductivity ranges for the inclusions were determined from initial higher-pass (larger filtering radius in the scattering domain) D-bar 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 non-negative conductivities bounded away from zero and infinity.Nevertheless, the method provides useful conductor/resistor information.The same k and z grids were used in the D-bar reconstructions of the KIT4 example as in the ACT4 example.However, we reduced the non-uniform 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.

C. Computational Notes for D-bar Reconstructions from Experimental EIT Data
The D-bar reconstructions from the experimental ACT4 and KIT4 data were computed in the same manner as the simulated data case described above in Section IV-B with the exception of the formation of the DN matrices L σ and L 1 , 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 trigonometric current patterns of amplitude 1mA had been applied by using a change of basis.
Using the approach introduced in [23], the (m, n) entry of the ND matrix R σ was formed as the discrete inner product where φ m denotes the m-th normalized current pattern, v n the voltage resulting from applying the m-th current pattern (normalized such that the voltages sum to zero and are scaled by the 2 norms of the applied current patterns), and |e | denotes the area of the -th electrode.This formula holds for a system with L electrodes where L − 1 linearly independent current patterns have been applied.The discrete DN matrix L σ was then formed by L σ = (R σ ) −1 and scaled by r σ0 as described above.As the scattering data t exp (5) requires the difference in DN matrices (L σ − L 1 ), the discrete DN matrices L 1 must be formed for both the ACT4 and KIT4 system with L = 32 and L = 16 electrodes, respectively.To this end, the conductivity equation in (1) was solved with boundary conditions given by the Complete Electrode Model [31] on a FEM mesh with 3,493 triangular elements simulating injected trigonometric current patterns of amplitude 1mA and non-optimized constant contact impedances of 0.00024.
The scattering data t exp (5) was evaluated as a simple Simpson's rule approximation where z is the vector of the positions of the centers of the electrodes, T denotes the traditional matrix transpose, and is the expansion of the asymptotic behavior ψ ∼ e ikz at the center of the -th electrode, z , in the basis of normalized applied current patterns φ.Then, the d-bar equation ( 6) can be solved using Fast Fourier Transforms using Vainikko's method [34].The interested reader is referred to [33] for further details.

V. RESULTS
We now demonstrate the effect of the Deep D-bar method on simulated as well as experimental data.

A. 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 Structural Similarity Indices (SSIMs) were computed for the ACT4 and KIT4 examples, shown in Figures 8 and 9, respectively.Additionally, we evaluated the minimized 2 -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 28.05% to 9.92% and for the KIT4 test data from 16.82% to 9.12% relative 2 -error.

B. Reconstructions from Experimental Data
Now we proceed to the reconstructions from experimental data.Figure 10 depicts the results of the Deep D-bar approach on four experiments with ACT4 data: HEALTHY and INJURIES 1-3 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 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 [17] 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 Low-pass D-bar and Deep D-bar reconstructions are shown on the same scale, which does mask the true dynamic range of the Low-pass D-bar 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 L 1 for the ACT4 and KIT4 experimental data examples was not finely tuned to either EIT device (which is required for iterative minimization-based methods), and 3) the D-bar solver was not optimized for the respective ACT4/KIT4 data.Rather it was used merely to provide the low-pass 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.
Furthermore, evaluation of the CNN is highly efficient on a GPU and takes in average 7.65ms for a single sample, hence we expect Deep D-bar to be real-time capable.This can be done by combining the D-bar reconstruction, as outlined in [24], with the application of the CNN in a unified framework to reduce overhead due to data transmission. 2n 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 to changes in the system.Furthermore, Figures 6 and 10 illustrate that the network can handle reconstructions of phantoms that do not conform with the training data.Additionally the robustness of obtaining the reconstructions as input to the network is crucial, here we point out that the D-bar algorithm is known to be robust with respect to perturbations of electrode positions or domain shapes [12].

VII. CONCLUSION
The D-bar method for 2D EIT provides reliable reconstructions of the conductivity but suffers from a blurring due to a low-pass filtering of the scattering data.Sharp improvements in image quality can be achieved by coupling the D-bar 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 D-bar method provides real-time conductivity reconstructions [24], the post-processing 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 D-bar computational framework has been further developed.

Fig. 1 .
Fig. 1.Deep D-bar network structure.The input is given by the D-bar reconstruction σ exp with a resolution of 64 × 64 and the output is denoted by σ.The numbers on top of the blue bars denote the channels for each layer.The resolution for each multilevel decomposition is shown in gray on the left.Each convolutional layer is equipped with a Rectified Linear Unit as nonlinearity, given by ReLU(x) = max(0, x).

Fig. 2 .
Fig.2.Experimental Setups for test phantoms taken on the ACT4 system from RPI. Agar/graphite targets were used to simulate a chest phantom with a heart, two lungs, aorta, and spine.The first image shows the healthy phantom.Three injuries are explored: 'Injury 1', replaced the cut portion of the right lung with agar/graphite of the same conductivity as the heart target to simulate a potential pleural effusion, 'Injury 2', replaced the cut portion of the right lung replaced with three plastic tubes, and'Injury 3', replaced the cut portion with three copper tubes.

Fig. 3 .
Fig. 3. Experimental Setups with conductive and resistive targets on the KIT4 EIT system from UEF.The white objects are made of solid plastic and are resistive.The hollow circular objects are conductive metal rings.

Fig. 4 .
Fig. 4. Depiction of the simulation of the training data for the ACT4 experiments of Figure 2. The first image shows the healthy phantom from which the 'true boundary' (black dots) and 'approximate boundary' (red stars) were extracted, shown in the second image.The third and fourth images display sample simulated phantoms using in the training data with and without injuries with the true boundaries overlaid in black dots.
and uniformly spaced 64 × 64 k and z-grids on [−4.5, 4.5] 2 with stepsize h k = 0.3234, and [−1, 1] 2 with stepsize h z = 0.0317, respectively.A non-uniform cutoff threshold was enforced on the scattering data for frequencies such that t exp (k) = 0 if either | (t exp (k)| or | (t exp (k)| exceeds 24.Then, the 4096 pairs of data in the form of 'Truth' and 'Low-pass D-bar Reconstruction' were used to train the convolutional neural network described in Section III.

Fig. 5 .
Fig. 5. Depiction of the simulation of the training data for the KIT4 experiments of Figure 3.The images shown are sample simulations of inclusions present in the training data.Zero to three, non-overlapping, circular inclusions were allowed in each simulation.

Fig. 6 .
Fig. 6. Results for simulated test data from the ACT4 geometry.The phantom in the first row conforms with the training data and the phantoms in the second and third row include pathologies not supported by the training data.The initial D-bar image is compared to the Deep D-bar image.The D-bar images, on the full square are used as the 'input' images for the CNN.Images are displayed here on the circular geometry of the tank, for presentation only.Each row is plotted on its own scale.

Fig. 7 .Fig. 8 .Fig. 9 .
Fig. 7. Results for simulated test data from the KIT4 geometry.All phantom are drawn from the same distribution as the training data.The initial D-bar image is compared to the Deep D-bar image.The D-bar images, on the full square are used as the 'input' images for the CNN.Images are displayed here on the circular geometry of the tank, for presentation only.Each row is plotted on its own scale.

Fig. 10 .
Fig. 10.ACT4 Results for the various test scenarios: Healthy, Injuries 1-3 corresponding to conductive agar, plastic tubes, and conductive copper tubes, respectively.The initial D-bar image is compared to the Deep D-bar image.The D-bar images, on the full square are used as the 'input' images for the CNN.Images are displayed here on the circular geometry of the tank, for presentation only.Each row is plotted on its own scale.

Fig. 11 .
Fig. 11.SSIM measurements are compared for the D-bar method and the new 'Deep D-bar' method for the ACT4 experimental data reconstructions shown in Figure 10.Note that meaningful SSIMs could not be computed for 'Injury 3' due to the copper inclusions which have infinite conductivity.

1 Fig. 12 .
Fig. 12. KIT4 Results for the various phantoms with conductive and/or resistive targets, as shown in the first column.The initial D-bar image is compared to the Deep D-bar image.The D-bar images, on the full square are used as the 'input' images for the CNN.Images are displayed here on the circular geometry of the tank, for presentation only.Each row is plotted on its own scale.