Regularized reconstruction in electrical impedance tomography using a variance uniformization constraint.

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

Image reconstruction in electrical impedance tomography using the sensitivity theorem is generally based on the assumption that the initial conductivity distribution of the body being imaged is uniform. The technique of image reconstruction using this method is described and reconstructed images are presented. Improvements in image quality and accuracy are demonstrated when accurate a priori `anatomical' information, in the form of a model of the distribution of conductivity within the region to be imaged, is used to construct the sensitivity matrix. In practice correct a priori information is not available, for example the conductivity values within the various anatomical regions will not be known. An iterative algorithm is presented which allows the conductivity parameters of the a priori model to be determined during reconstruction. Keywords: electrical impedance tomography, image reconstruction, a priori information

Image reconstruction in electrical impedance tomography is, generally, an ill-posed nonlinear inverse problem. Regularization methods are widely used to ensure a stable solution. In this paper we present a case study, which uses a novel electrical impedance tomography method for reconstruction of layered biological tissues with piece-wise continuous plane-stratified profiles. The algorithm implements the recently proposed reconstruction scheme for piecewise constant conductivity profiles, based on an improved Prony method in conjunction with Bessel transform. It is shown that the proposed algorithm is capable of successfully reconstructing piece-wise continuous conductivity profiles with moderate (WKB) slop. This reconstruction procedure, which calculate both the locations and the conductivities, repetitively provides inhomogeneous depth discretization, i.e., the depths grid is not equispaced. Incorporation of this specific inhomogeneous grid, the main result of our investigation, in the widely used MLS reconstruction procedure results in a stable reconstruction, whereas, utilization of equispaced depth grid leads to inaccurate reconstruction. The effectiveness and the significance of inhomogeneous discretization in EIT reconstruction procedure is further demonstrated and verified via numerical simulations.

Resistivity reconstruction of practical phantoms is studied with Projection Error Propagation-based Regularization (PEPR) method to improve the image quality in Electrical Impedance Tomography (EIT). PEPR method calculates the regularization parameter as a function of the projection error produced due to the mismatch between experimental measurements and the calculated data. The regularization parameter in the reconstruction algorithm automatically adjusts its magnitude depending on the noise level present in measured data as well as the ill-posedness of the Hessian matrix. Resistivity images are reconstructed from the practical phantom data using the Electrical Impedance Diffuse Optical Reconstruction Software (EIDORS). The resistivity images reconstructed with PEPR method are compared with the Single step Tikhonov Regularization (STR) and Modified Levenberg Regularization (LMR) techniques. The results show that, the PEPR technique reduces the reconstruction errors in each iteration and improves the reconstructed images with better contrast to noise ratio (CNR), percentage of contrast recovery (PCR), coefficient of contrast (COC) and diametric resistivity profile (DRP).KeywordsElectrical Impedance Tomography (EIT)Projection Error Propagation-based Regularization (PEPR)Practical PhantomsBoundary potentialsEIDORSResistivity reconstructionSTRLMR

Impact of the Radial Basis Function Spread Factor onto Image Reconstruction in Electrical Impedance Tomography

A real-time EIT imaging system based on the split augmented Lagrangian shrinkage algorithm

Image reconstruction in electrical impedance tomography (EIT) deals with an ill-posed and non-linear inverse problem. It has the objective of minimizing the difference between simulated (virtual) object data and electric voltage measurements performed on a non-simulated (real) object. In this paper, a new approach to the simulated annealing method applied to the reconstruction of EIT images is described. The main advantage of this approach is that all conductivity parameters are updated simultaneously. Most methods that employ simulated annealing to the problem of EIT usually evaluate each conductivity parameter individually resulting in high computational cost. The algorithm was tested both with computationally generated data and with measurements performed on a physical tank. In both cases, the method was able to make data inversion, determining the position, the dimensions and the conductivity of materials in an opaque object plane.

Comparison of statistical inversion with iteratively regularized Gauss Newton method for image reconstruction in electrical impedance tomography

Image reconstruction in electrical impedance tomography (EIT) is a non-linear inverse problem. The linear model is always used in most of the reconstruction algorithm in dynamic EIT, which causes large errors of image reconstruction. In this paper, we proposed a new image reconstruction method in dynamic EIT. In the method the artificial neural network based on the error back propagation algorithm (BP ANN) is used to express the non-linear relation between the impedance change position inside the measured object and the voltage change value measured on the surface of the object. Thus, the location ofthe impedance change can be decided by the measured voltage change on the surface, and then the impedance change image will be reconstructed with linear approximated method. The reconstructed error will be decreased largely, because the impedance change position can be detected precisely by our proposed method. The experimental results indicate that the precision of the reconstructed image with our method is greatly higher than that with the back projection method.

Shape-based solutions have recently received attention for certain ill-posed inverse problems. Their advantages include implicit imposition of relevant constraints and reduction in the number of unknowns, especially important for nonlinear ill-posed problems. We apply the shape-based approach to current-injection electrical impedance tomography (EIT) reconstructions. We employ a boundary element method (BEM) based solution for EIT. We introduce two shape models, one based on modified B-splines, and the other based on spherical harmonics, for BEM modeling of shapes. These methods allow us to parameterize the geometry of conductivity inhomogeneities in the interior of the volume. We assume the general shape of piecewise constant inhomogeneities is known but their conductivities and their exact location and shape is not. We also assume the internal conductivity profile is piecewise constant, meaning that each region has a constant conductivity. We propose and test three different regularization techniques to be used with either of the shape models. The performance of our methods is illustrated via both simulations in a digital torso model and phantom experiments when there is a single internal object. We observe that in the noisy environment, either simulated noise or real sources of noise in the experimental study, we get reasonable reconstructions. Since the signal-to-noise ratio (SNR) expected in modern EIT instruments is higher than that used in this study, these reconstruction methods may prove useful in practice.

Magnetic induction tomography (MIT) of biological tissue is used for the reconstruction of the complex conductivity distribution κ inside the object under investigation. It is based on the perturbation of an alternating magnetic field caused by the object and can be used in all applications of electrical impedance tomography (EIT) such as functional lung monitoring and assessment of tissue fluids. In contrast to EIT, MIT does not require electrodes and magnetic fields can also penetrate non-conducting barriers such as the skull. As in EIT, the reconstruction of absolute conductivity values is very difficult because of the method's sensitivity to numerical errors and noise. To overcome this problem, image reconstruction in EIT is often done differentially. Analogously, this concept has been adopted for MIT. Two different methods for differential imaging are applicable. The first one is state-differential, for example when the conductivity change between inspiration and expiration in the lung regions is being detected. The second one is frequency-differential, which is of high interest in motionless organs like the brain, where a state-differential method cannot be applied. An equation for frequency-differential MIT was derived taking into consideration the frequency dependence of the sensitivity matrix. This formula is valid if we can assume that only small conductivity changes occur. In this way, the non-linear inverse problem of MIT can be approximated by a linear one (depending only on the frequency), similar to in EIT. Keeping this limitation in mind, the conductivity changes between one or more reference frequencies and several measurement frequencies were reconstructed, yielding normalized conductivity spectra. Due to the differential character of the method, these spectra do not provide absolute conductivities but preserve the shape of the spectrum. The validity of the method was tested with artificial data generated with a spherical perturbation within a conducting cylinder as well as for real measurement data. The measurement data were obtained from a potato immersed in saline. The resulting spectra were compared with reference measurements and the preservation of the shape of the spectra was analyzed.

This paper introduces a new method of simultaneous acquisition of electrocardiogram (ECG) signals and electrical impedance tomography (EIT) measurements on active electrodes, i.e. electrodes that are applying current, in an EIT system. Howland noise reduction followed by an integrate-and-dump filter and adaptive filtering are used to extract the ECG signal without loss of amplitude. Results of simultaneous ECG signals and EIT reconstructions on human subjects are shown. It is also demonstrated that the recovered ECG waveforms from 32 electrodes can be used to reconstruct the heart current source vectors during a cardiac cycle, yielding an approximate solution to the inverse problem of electrocardiography. Acquiring ECG signals from all electrodes attached to the body simultaneously with EIT data acquisition eliminates the need for a separate system to collect ECG signals. Time-aligned ECG signals are helpful in interpreting pulsatile perfusion images by providing exact timing of each image in relation to the cardiac cycle. Additionally, collecting ECG signals from many or all of the EIT electrodes enables reconstruction of the heart's moving dipole moment at the same time as the EIT images.

Image reconstruction for Electrical Impedance Tomography (EIT) is a highly nonlinear and ill-posed inverse problem. It requires the design and employment of feasible reconstruction methods capable to guarantee trustworthy image generation. Deep Neural Networks (DNN) have a powerful ability to express complex nonlinear functions. This research paper introduces a novel framework based on DNN aiming to achieve EIT image reconstruction. The proposed DNN model, comprises of the following two layers, namely: The Stacked Autoencoder (SAE) and the Logistic Regression (LR). It is trained using the large lab samples which are obtained by the COMSOL simulation software (a cross platform finite elements analysis solver). The relationship between the voltage measurement and the internal conductivity distribution is determined. The untrained voltage measurement samples are used as input to the trained DNN, and the output is an estimate for image reconstruction of the internal conductivity distribution. The results show that the proposed model can achieve reliable shape and size reconstruction. When white Gaussian noise with a signal-to-noise ratio of 30, 40 and 50 were added to test set, the proposed DNN structure still has good imaging results, which proved the anti-noise capability of the network. Furthermore, the network that was trained using simulation data sets, would be applied for the EIT image reconstruction based on the experimental data that were produced after preprocessing.

Sparsity reconstruction in electrical impedance tomography: An experimental evaluation

Electrical impedance tomography (EIT) is an imaging modality based on electrical currents injected through the imaging domain and voltages measured on the boundary. It has a wide variety of applications ranging e.g. from medical imaging to industrial process monitoring. Computing the EIT reconstruction from the measurements is a severely ill-posed inverse problem, and hence it is not a straightforward task. Therefore, the EIT reconstruction algorithms remain an active area of research. To facilitate this research, we present in this article object-oriented EIT (OOEIT), an open source software package for MATLAB. It is designed for ease of use for those with no deep understanding of EIT, while also featuring a modular structure that allows for straightforward customization for those interested in developing their own algorithms. OOEIT includes implementations for both 2D and 3D computations, for both static and dynamic reconstructions, and supports user defined meshes. As examples, we present reconstructions computed from the Kuopio Tomography Challenge (KTC2023) data set.