Next Article in Journal
Ambient Light Rejection Integrated Circuit for Autonomous Adaptation on a Sub-Retinal Prosthetic System
Next Article in Special Issue
SAFIR-I: Design and Performance of a High-Rate Preclinical PET Insert for MRI
Previous Article in Journal
High-Throughput Phenotyping Approach for the Evaluation of Heat Stress in Korean Ginseng (Panax ginseng Meyer) Using a Hyperspectral Reflectance Image
Previous Article in Special Issue
Feasibility Study of a Time-of-Flight Brain Positron Emission Tomography Employing Individual Channel Readout Electronics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Kwancheol Shin
1,* and
Jennifer L. Mueller
2
1
Department of Mathematics, Chungbuk National University, Cheongju 28644, Korea
2
Department of Mathematics and School of Biomedical Engineering, Colorado State University, Fort Collins, CO 80523, USA
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(16), 5635; https://doi.org/10.3390/s21165635
Submission received: 16 June 2021 / Revised: 31 July 2021 / Accepted: 13 August 2021 / Published: 21 August 2021
(This article belongs to the Collection Biomedical Imaging & Instrumentation)

Abstract

:
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.

1. Introduction

Electrical impedance tomography has been demonstrated to be a promising technique for pulmonary imaging at the bedsides of critically ill patients and for those with chronic lung disease. Recent survey articles on clinical applications of EIT include [1,2] and the reader is also referred to the articles [3,4,5,6,7,8,9] for an introduction.
Calderón’s method is a direct reconstruction method based on complex geometrical optics (CGO) solutions to Laplace’s equation. The method was outlined by Calderón in his 1980 paper On an inverse boundary value problem [10], in which he proved that the linearized problem has a unique solution and proposed a direct method of reconstruction. The method can be viewed as a linearization of the D-bar method [11,12] based on the 1996 global uniqueness proof by Nachman [13]. While the D-bar method makes use of a certain type of nonlinear Fourier transform known as the scattering transform as an intermediate function in the solution of a ¯ (D-bar) equation to directly compute the conductivity, Calderón’s method relies on computing the inverse Fourier transform of a bilinear form involving the data and CGO solutions to obtain the conductivity. Both methods require knowledge of the Dirichlet-to-Neumann or voltage-to-current density map. Calderón’s method was first implemented numerically in [14] for experimental EIT data collected on a saline-filled tank and on a healthy human subject reconstructed on a circular domain. In [15], this was extended to noncircular but symmetric domains. A real-time implementation in 2-D on subject-specific domains was presented in [16,17]. A “higher-order” Calderón’s method and a method of including a spatial prior was presented in [12] with demonstrated results from data collected on a circular saline-filled tank with agar and vegetable targets. A 3-D Calderón-based method was provided in [18] in the planar geometry and [19] for a cylinder.
In this paper, an implementation in which the exact locations of the electrodes are known is provided. Previous work on a subject-specific domain [16,17] assumed that the electrodes were equally spaced in their angular positions around the boundary, leading to some simplification in the computational formulas. A further novel contribution of this paper is a method of including a static prior derived from CT scans and a preliminary reconstruction to compute dynamic images of ventilation and perfusion. The relevant medical application for this scenario would be a case when a CT scan and EIT scan are performed upon hospitalization, and then the patient is monitored with EIT for pulmonary complications, such as a pneumothorax or pleural effusion, both of which are conditions that may develop, for example, in patients with acute respiratory distress syndrome (ARDS) receiving mechanical ventilation.
The use of a spatial prior to improve the resolution of reconstructed EIT images has proved to be effective in iterative methods [20,21,22,23,24,25,26,27,28,29] and in the D-bar method [30,31,32,33]. In the D-bar method, the scattering transform computed from the data is appended with a scattering transform computed from a conductivity prior. In [34], a statistical approach using the Schur complement was used to construct a prior for the D-bar method. In [12], the Fourier transform of the prior conductivity distribution is appended to the bilinear form before inverting. Dynamic images of the human chest prove more challenging in the construction of a prior conductivity distribution. In [35] a dynamic prior was constructed for the D-bar method applied to ventilation data from inspiratory and expiratory CT scans of the patient. In this work, a static prior is constructed by segmenting a CT scan of a patient with cystic fibrosis, and the conductivity values for the lung, heart, and background are assigned from a preliminary reconstruction with Calderón’s method.
Since Calderon’s method has attracted attention as a direct real-time reconstruction method, we aim to show its potential for practical implementation, which has both advantages and shortcomings. Other methods may result in higher resolution images, some with significant trade-offs, such as computation time. Gauss–Newton methods require a forward model for the iterative updates, which to date precludes real-time reconstruction. However, careful modeling also contributes to improvements in resolution. Including a prior further improves performance [20,21,22,23,24,25,26,28,29]. The D-bar method is a direct reconstruction method for which spatial priors have been included by appending the scattering transform computed from the measured data with the scattering transform of the prior. The D-bar method also can be implemented in real time, although Calderón’s method produces reconstructions even faster [16]. Shape reconstruction methods [36,37,38,39] have the advantage that they can incorporate geometry and prior information directly and preserve sharp edges, also reducing the computational burden of the full reconstruction problem, which also shows promise for lung imaging.
The paper is organized as follows. In Section 2, the mathematical model of EIT is provided. In Section 3, Calderón’s method and our numerical implementation with exact electrode modeling and the inclusion of a spatial prior are presented. The results of the method on experimental data collected on a chest-shaped tank, simulated data with several simulated pulmonary pathologies, and human data collected on a patient are provided in Section 4. The conclusion is presented in Section 5.

2. Background

Modeling of EIT

Let γ ( x ) γ 0 > 0 be the electrical conductivity with a positive constant lower bound γ 0 , u ( x ) be the electrical potential and Ω R 2 be a bounded domain. The governing equation of EIT is
· γ ( x ) u ( x ) = 0 , x Ω .
The applied current density j on the boundary corresponds to the Neumann boundary condition
γ ( x ) u ν ( x ) = j ( x ) , x Ω ,
where ν is the outward normal vector to Ω . We denote the voltage distribution on the boundary by f so that
u ( x ) = f ( x ) , x Ω ,
is the corresponding Dirichlet boundary condition for (1). If γ ( x ) and one of the boundary conditions (2) or (3) are given, the forward problem is to solve for u ( x ) in Ω . The inverse problem is to find the unknown γ ( x ) from knowledge of the Neumann-to-Dirichlet (ND) map,
R γ : γ ( x ) u ν ( x ) u ( x ) , x Ω ,
which is also called the current density-to-voltage map. However, in most of the mathematical literature, the theory is developed using the Dirichlet-to-Neumann (DN) map
Λ γ : u ( x ) γ ( x ) u ν ( x ) , x Ω ,
or voltage-to-current density map.

3. Methods

3.1. Calderon’s Method

For completeness, we summarize Calderón’s method from [10]. Assume that γ ( x ) is a small perturbation δ ( x ) L ( Ω ) from a background conductivity of 1 so that γ ( x ) = 1 + δ ( x ) . Calderón’s method uses a special type of harmonic function known as complex geometrical optics (CGO) solutions defined by
f ( x ; k , a ) = e π i ( k · x ) + π ( a · x ) , g ( x ; k , a ) = e π i ( k · x ) π ( a · x ) ,
where k , a R 2 are nonphysical frequency variables with | k | = | a | , k · a = 0 . Next, let ω i be defined on Ω as
ω i = u i + v i in Ω , i = 1 , 2 ,
where ω 1 = f and ω 2 = g on Ω , u i H 1 ( Ω ) , and v i H 0 1 ( Ω ) for i = 1 , 2 . Here, H 1 ( Ω ) denotes the Sobolev space of functions with one weak derivative in L 2 ( Ω ) and H 0 1 ( Ω ) the Sobolev space H 1 ( Ω ) with trace zero [40]. Then, by (1) and integration by parts,
Ω f ( x , k ) Λ γ g ( x , k ) d s ( x ) = 2 π 2 | k | 2 Ω γ ( x ) exp 2 π i x · k d x + R ˜ ( k ) ,
where R ˜ ( k ) is considered to be a remainder term. Rearranging terms results in
Ω γ ( x ) exp 2 π i x · k d x = 1 2 π 2 | k | 2 Ω f ( x , k ) Λ γ g ( x , k ) d s ( x ) + R ^ ( k ) ,
where R ^ ( k ) = R ˜ ( k ) 2 π 2 | k | 2 . Assuming γ ( x ) is constant outside of Ω , the left-hand side is the Fourier transform of γ χ Ω , which we denote by γ ^ ( k ) . We denote the first term on the right-hand side by
F ^ γ abs ( k ) 1 2 π 2 | k | 2 Ω f ( x , k ) Λ γ g ( x , k ) d s ( x ) .
Calderón showed that when | | δ ( x ) | | L ( Ω ) is small, | R ^ ( k ) | is small for small values of | k | . Since γ ( x ) χ Ω is zero outside of Ω , | γ ^ ( k ) | has to decrease to zero as | k | becomes larger. Therefore, by multiplying a mollifying function η ^ ( k ) such that η ^ ( 0 ) = 1 , η ^ C and that decreases fast to zero as | k | becomes larger, F ^ γ abs η ^ approximates γ ^ , and it is shown in [10] that the inverse Fourier transform of R ^ ( k ) η ^ ( k ) is negligible when | | δ | | L is small. The use of the mollification is equivalent to applying a low-pass filter to F ^ γ abs . This low-pass filtering can be carried out alternatively by truncating F ^ γ abs ( k ) into a disk of radius R 1 . We call R 1 the truncation radius. Therefore, we obtain an approximation of γ ( x ) by
γ abs ( x ) | k | R 1 F ^ γ abs ( k ) η ^ ( k / t ) e 2 π i x · k d k ,
where η ^ ( k / t ) is a mollifying function with η ^ ( 0 ) = 1 and t is the mollfication parameter. We call the reconstruction from (7) the absolute image of γ.
Let Λ 1 denote the Dirichlet-to-Neumann map corresponding to a constant conductivity distribution of γ = 1 . By subtracting F ^ 1 abs ( k ) from F ^ γ abs ( k ) , defined in Equation (6), we define
F ^ hom diff ( k ) 1 2 π 2 | k | 2 Ω f ( x , k ) ( Λ γ Λ 1 ) g ( x , k ) d s ( x ) ,
Inverting (8) on the truncated domain | k | R 1 results in an approximation of the perturbation δ ( x ) , which we denote by
δ hom diff ( x ) | k | R 1 F ^ hom diff ( k ) η ^ ( k / t ) e 2 π i x · k d k ,
and call the difference image of δ with respect the homogeneous data.
Note, however, that the computation of the difference image requires knowledge of Λ 1 , which in most practical applications such as in medical imaging, cannot be measured. On the other hand, the action of Λ 1 on g ( x , k ) can be computed analytically,
Ω f ( x , k ) Λ 1 g ( x , k ) d x = 2 π 2 | k | 2 Ω e 2 π i k · x d x ,
and Equation (8) can be replaced by
F ^ δ abs ( k ) 1 2 π 2 | k | 2 Ω f ( x , k ) Λ γ g ( x , k ) d s ( x ) Ω e 2 π i k · x d x .
We define
δ abs ( x ) | k | R 1 F ^ δ abs ( k ) η ^ ( k / t ) e 2 π i x · k d k ,
and call it the the absolute image of δ.
While one can compute γ ( x ) directly by using (7), the use of (9) or (11) avoids the Gibbs phenomena near Ω , which is introduced by the sharp discontinuity in γ ( x ) χ Ω along Ω . Under the assumption that δ ( x ) is zero and therefore flat near Ω , there are no Gibbs phenomena in the reconstructed images. Gibbs phenomena in the absolute images of γ can be observed in [12,14].
Another possibility that is particularly relevant to medical imaging is to choose one frame in a dynamic sequence of images as a reference image. For ventilatory images, this is often the frame corresponding to maximal exhalation and for perfusion images, this could correspond to peak systole, peak diastole, or a time point in between. Denoting the Dirichlet-to-Neumann map for the reference frame by Λ γ ref , subtracting F ^ γ ref abs for γ ref from Equation (6) results in
F ^ ref diff ( k ) 1 2 π 2 | k | 2 Ω f ( x , k ) ( Λ γ Λ γ ref ) g ( x , k ) d s ( x ) ,
Taking the Fourier transform to (12) in the region | k | R 1 , we obtain an approximation to γ γ ref which we will denote by δ ref diff ( x ) and call it the difference image of δ with respect to the reference frame. We will make clear whether we are using a homogeneous data as in (8) and (9) or a specific frame in a dynamic sequence of data as in (12) when we refer to the difference image.

3.2. Numerical Implementation

Let { T i } be a set of linearly independent current patterns, i = 1 , , L 1 , applied to L electrodes. For i = 1 , , L 1 , l = 1 , , L , let T l i denote the ith current pattern on the lth electrode and V l i the measured voltage. By Kirchhoff’s law, the elements of T i must satisfy l = 1 L T l i = 0 , and for a unique choice of ground we require l = 1 L V l i = 0 . Let t i denote the normalized current t i = ( T i ) / ( | | T i | | 2 ) and v i the normalized voltage v i = ( V i ) / ( | | T i | | 2 ) , where | | T i | | 2 = l = 1 L ( T l i ) 2 . In this work, the adjacent current patterns were applied on the electrodes and are given by
T l i = M , l = i , M , l = i + 1 mod L , 0 , otherwise .
We model the current density j ( x ) on the boundary by the gap model,
j ( x ) = I l A , x e l 0 , otherwise ,
where e l denotes the lth electrode and A is the area of an electrode, which is assumed to be the same for all electrodes.
Assume the locations of the centers of each of the electrodes is known, and let x | Ω = ( x 1 ( θ ) , x 2 ( θ ) ) = r ( θ ) ( cos θ , sin θ ) be a parameterization of the boundary of Ω by θ , and the line element d s = ( x 1 ) 2 + ( x 2 ) 2 d θ = r 2 + r 2 d θ . Denote the angle that the center of the lth electrode makes with a reference point chosen to be 0 degrees by ( θ l ) l = 1 l = L , and the angle between the center of the lth and l + 1 st electrode by Δ θ l = θ l + 1 θ l . For functions r , s of θ such that r , s : R L R , let ( r ( · ) , s ( · ) ) L denote the discrete inner product defined by ( r ( · ) , s ( · ) ) L = l = 1 L r ( θ l ) s ( θ l ) .
Expand f ( x l , k ) and g ( x l , k ) with respect to the normalized current patterns and the measured voltages as follows:
f ( x l , k ) = i = 1 L 1 f k i t l i ,
g ( x l , k ) = j = 1 L 1 g k , γ j v l j , γ ,
and let f k = f k i i = 1 L 1 and g k , γ = g k , γ j j = 1 L 1 be the coefficient vectors, where the subscripts k and γ indicate the dependence of the coefficients on the variable k and the conductivity γ . Since the applied normalized current patterns and the measured voltages are not orthogonal in general, in order to compute f k and g k , γ , we need to solve systems of linear equations as follow. Taking the inner product with v i , γ ( i = 1 , 2 , . . . , L 1 ) on both sides of (14), we obtain V γ g k , γ = c k , γ , where V γ ( i , j ) = ( v i , γ , v j , γ ) L and c k , γ ( i ) = g ( x l , k ) , v l i , γ L . Therefore, g k , γ = V γ 1 c k , γ . Similarly, f k = T 1 d k , where T ( i , j ) = ( t i , t j ) L and d k ( i ) = f ( x l , k ) , t l i L . Now, with (13) and (14),
Ω f ( x , k ) Λ γ g ( x , k ) d s ( x ) = 0 2 π i = 1 L 1 f k i t i ( θ ) Λ γ j = 1 L 1 g k , γ j v j , γ ( · ) d s ( θ ) = i = 1 L 1 j = 1 L 1 f k i g k , γ j 0 2 π t i ( θ ) Λ γ v j , γ ( · ) d s ( θ ) = 1 A i = 1 L 1 j = 1 L 1 f k i g k , γ j l = 1 L t l i t l j e w , = e w A f k T T g k , γ ,
where e w is the width of one electrode, where all are assumed to be equal, and f k T is the transpose of f k . Notice that the coefficients f k and g k , γ encode the information about the location of electrodes and therefore the boundary shape provided that we use the exact values of f ( x l , k ) and g ( x l , k ) in (13) and (14). From (8), (10), and (12), and by the similar calculation,
F ^ homo diff ( k ) = e w 2 π 2 A | k | 2 f k T T ( g k , γ g k , 1 ) ,
F ^ δ abs ( k ) = e w 2 π 2 A | k | 2 f k T T g k , γ Ω e 2 π i k · x d x ,
and
F ^ ref diff ( k ) = e w 2 π 2 A | k | 2 f k T T ( g k , γ g k , ref ) .
The approximations δ homo diff , δ abs and δ ref diff are computed by taking the inverse Fourier transform of (15)–(17), respectively. The inverse Fourier transform is computed for | k | R 1 , where R 1 is the truncation radius. This radius is unitless, and it depends on the scaling of the domain. If, for example, the maximum value of | x | is 180 mm, and we scale the domain to have a maximum of | x | = 1 ; then, we rescale the truncation radius to R 1 180 Simpson’s quadrature rule is used to compute the integral of the inverse Fourier transform. Once we compute δ homo diff , δ abs and δ ref diff , we add the background conductivity to those to get the conductivity distributions γ homo diff , γ abs and γ ref diff , respectively.

3.3. Calderón’s Method with a Spatial Prior

In [12], a method of including a spatial prior to Calderón’s method is introduced. We implement the method in this paper. We assume that we have a priori knowledge for the boundary of organs and the approximate regional conductivity for each organ. We denote this a priori conductivity distribution by δ pr ( x ) . By adding a high-pass filtered prior R 1 | k | R 2 δ ^ pr ( k ) η ^ ( k / t ) e 2 π i k · x d k to the reconstructions obtained from (11) and (12), we define
δ pr abs ( x ) = δ abs ( x ) + pri δ ^ pr ( k ) η ^ ( k / t ) e 2 π i k · x d k ,
and
δ pr diff ( x ) = δ diff ( x ) + R 1 | k | R 2 δ ^ pr ( k ) η ^ ( k / t ) e 2 π i k · x d k ,
for some constant R 2 , where δ diff is either δ homo diff or δ ref diff depending on the applications.We define γ pr abs and γ pr diff by adding the background conductivity to (18) and (19), respectively.

4. Experimental Results

4.1. Chest Shaped Tank Data

In this section, the improvement in resolution when correct electrode positions are included in the algorithm is demonstrated. In this section, no spatial prior is used in the reconstruction algorithm. Data was collected on a chest-shaped tank of with circumference 1.016 m, simulating the shape of a human subject. The tank was filled with saline of conductivity of 0.2 S/m to a height of 0.0204 m, and three inclusions featuring two lungs of conductivity 0.09 S/m and the heart of conductivity 0.45 S/m. The width of the electrodes was 0.254 m. From the center of the tank, the electrodes are spaced nonuniformly in angle. For the homogeneous data for the difference images of δ , a data set was collected with only saline in the tank. The data was taken with the Active Complex Electrode (ACE1) system (see, [41,42]) in the EIT lab at Colorado State University. The frequency of the system was 125 kHz and the current amplitude was 3.3 mA.
See Figure 1 for a photo of the experimental configuration and reconstructions of a difference image of δ homo diff with the algorithm suggested in [16] in which the location of electrodes is assumed to be spaced uniformly in angle (upper right). Without correct electrode modeling, absolute images were unattainable. The truncation radius is 1.2 for all image. The lower left image in Figure 1 is a reconstruction of δ homo diff with correct electrode modeling suggested in this paper. The lower right image is an absolute image, δ abs with the same truncation radius. With correct electrodes modeling, the absolute image of δ is comparable to the difference image The L2 norms of the differences between those reconstructions and the ground truth are 0.1128, 0.0973 and 0.0981 in order.

4.2. Simulated Data on a 2D Chest-Shaped Tank with a Spatial Prior

The purpose of including a spatial prior in the reconstruction algorithm in the clinical setting is to improve resolution and detection of pulmonary pathology without biasing the reconstruction towards any particular pathology. With this goal, we take the point of view that the spatial prior should represent a normal chest with healthy lungs. To study the effectiveness of this approach, reconstructions from simulated data with 0.1 % added noise representing healthy lungs and four different pathologies were computed using a healthy lung prior. In all cases, the correct electrode locations were used in the reconstructions. Four simulated pathologies were studied: a small pneumothorax and a large pneumothorax in the left lung, a contusion in the right lung, and a pleural effusion in the left lung. Reconstructions are found in Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, respectively. The prior was constructing from the healthy lung case. The conductivity values for each of the organs and pathologies are found in Table 1, and were chosen in agreement with the references [43,44]. In the simulation, data with no inclusions are computed and used for the difference images.
In each example with a pathology, the pathology is clearly visible in the difference image created by subtracting the data from the simulated healthy heart and lung case from the data corresponding to the case with the pathology (the lower right image in each of these figures). This simulates detection of a pathology that occurs while the patient is being monitored, and pre-pathology data is available. The images of γ abs with no prior, γ diff with no prior, γ abs with the spatial prior, γ diff with the spatial prior have more subtle differences from one another. In the case of the small and large pneumothoraces, the left lung is visibly larger, and resistive region near the top of the left lung expands. In the simulation of the pulmonary contusion, the contusion is arguably invisible in the reconstructions, but clearly visible in the difference image in the lower right. In the simulation of the pleural effusion, the effusion is clearly visible in each of the reconstructions, although not well-localized.
The results suggest that the method is best-suited for monitoring for developing pathologies as opposed to detection of a pathology already present with no reference image in which the pathology is absent. We are not assuming that we can collect EIT data on both a healthy and injured lung at the same time. In the case of a contusion, which is caused by blunt-force trauma to the lung, the clinical scenario most relevant here would be monitoring the resolution or healing of the contusion, using the contusion image as reference.

4.3. Human Subject Data

We demonstrate the effectiveness of Calderón’s method with accurate shape and electrode modeling and a spatial prior on archival data collected as part of a larger study conducted in accordance with the amended Declaration of Helsinki—Ethical Principles for Medical Research Involving Human Subjects. Data were collected at Children’s Hospital Colorado, Aurora, CO, under the approval of the Colorado Multiple Institutional Review Board (COMIRB) (approval number COMIRB 14-0652). Informed written parental consent and children’s informed assent was obtained from the subject. Data from one 10-year-old male cystic fibrosis patient collected during a routine annual exam is considered here.
EIT data were collected during tidal breathing and during breathholding on 22 pediatric Philips EKG electrodes of height 3.33 cm and width 2.22 cm placed around the circumference of the chest. When the electrodes were removed, fiducial markers were placed at the center of each electrode, and a CT scan was performed as part of the subject’s standard care. Bipolar adjacent current patterns were applied with a current amplitude of 2.9 mA on the injecting electrodes. Electrode centers were identified from the inspiratory CT scans with the slices overlaid. The CT scan slide in Figure 7 was used to approximate the boundary shape for the domain and the prior as well as the shape and location of the lungs and heart to create the spatial prior. For the ventilation image sequence, a reference frame was chosen by finding the minimum value of the first component in a principal component analysis (PCA) plot of the data. This frame corresponds to maximum exhalation over the sequence and was the 200th frame. For the breath-holding sequence, the voltage data on each electrode for a fixed current pattern was averaged over all of the frames to create a data set representing a data set averaged over all frames.
After running a preliminary difference reconstruction of δ diff with no prior, the maximum conductivity value in the region of the heart was assigned to the heart region in the spatial prior, the minimum conductivity value in the region of the lungs was assigned to the lung regions in the spatial prior, and the background conductivity was set to zero, since these are difference images. Difference images of ventilation data with no prior and truncation radius 1.2 are found in Figure 8. While ventilatory changes are clearly visible, the spatial resolution is very poor. Difference images with a spatial prior with an inner truncation radius of 1.2 and outer truncation radius of 3 are found in Figure 9. The spatial resolution has improved when compared to the CT scan and the ventilatory changes are still clearly visible. In Figure 10, time traces of a pixel from the lung region and from the heart region are superimposed both for the no-prior reconstructions and the reconstructions with the static spatial prior.
Difference images from the breath holding data with no prior and truncation radius 1.2 are found in Figure 11, and reconstructions with the prior are found in Figure 12. In this case, a reference DN map was computed by averaging the voltages over all frames to compute an average reference frame. That is, for each fixed current pattern and each electrode, an average voltage over all frames was assigned to the voltage component in the computed reference data. The plotted images correspond to systole and diastole, as selected manually from the images. A movie of the full reconstruction sequence can be found in the Supplementary Data. The BIOPAC 3-lead EKG data are plotted in Figure 13. Since the ACE 1 system with 22 electrodes collects data at 36.28 frames/s, each frame corresponds to 0.276 s of data. Since start time of the BIOPAC data collection was delayed by 1.06 s, frame 100, the first frame found in Figure 11 and Figure 12, correspond to 3.82 s, and we display the relevant portion of the EKG output. For reference, each small rectangle in the EKG figure represents 0.04 s or 1.45 frames. Frames from Figure 11 and Figure 12 are marked with arrows, and the frames corresponding to systole have red box around them. One sees that systole in the EIT images closely coincides with the QRS complex in the EKG data, as would be expected. The average heart rate can be seen in the lower half of the figure. From the images, one also sees that the heart region is less red when the heart ejects blood during systole, and becomes deeper red during diastole. The lung regions shrink and become less blue during systole as blood enters the lungs and they become more conductive, and they become bluer as the heart fills during diastole, as the conductive blood leaves the lungs and returns to the heart. In Figure 10, we plot the time trace of two pixels which are indicated with squares in the first image of Figure 12. The left plot is the time trace of the ventilation data and the right one is for the breath-holding data. In both plots, each graph represents, from the top to the bottom, the heart region with and without a prior and the lung region with and without a prior. The two plots show the breath and the beating of the heart clearly. The left plot shows frames 190 to 450, demonstrating about 15.5 heart beats and five breath cycles. The right plot shows frames 100 to 196, depicting about 5.5 heart beats and the perfusion in the lung region. In the plot from the breath-holding data, the time trace of the lung pixel and the time trace of the heart pixel are out-of-phase, as is expected.
The full image sequences of reconstructions for Figure 8, Figure 9, Figure 11 and Figure 12 can be found as movies on https://youtu.be/0qxE4Aosmtw, https://youtu.be/oHTjZGnWLKQ, https://youtu.be/1PP7c9E17sM, and https://youtu.be/q_yugxxhdHM all accessed on 13 August 2021, in order.
Reconstructions of ventilation computed on a circular domain of the same perimeter as the true domain are included in Figure 14 and Figure 15. These images show separation between the lungs, as may be familiar to the reader for EIT reconstructions, but it is an artifact of the circular domain, which has the effect of pulling the targets toward the boundary, as has been previously reported in Figure 14 and Figure 15 for Calderon’s method. Due to the very close proximity of the lungs to each other near the heart (see the CT scan in Figure 7), the lungs should not show a large separation. In fact, the reconstructions from the prior show an indentation near the bottom (at the location of the spine), which is to be expected in a difference image, and this is further evidence of the improved resolution with the prior.

5. Conclusions

In this work, Calderón’s method in two dimensions was implemented with correct electrode positions for experimental and simulated data and shown to yield improved difference and absolute images compared to when the electrodes are assumed to be equally spaced. In addition, a spatial prior was included in the reconstruction algorithm and demonstrated to improve spatial resolution and accuracy of conductivity values. Examples from human data demonstrate the method’s ability to image conductivity changes due to ventilation and pulsatile pulmonary perfusion.

Supplementary Materials

The original files for Figure 8, Figure 9, Figure 10 and Figure 11 are available online at https://www.mdpi.com/article/10.3390/s21165635/s1.

Author Contributions

Formal analysis, K.S. and J.L.M.; investigation, K.S. and J.L.M.; methodology, K.S. and J.L.M.; validation, K.S. and J.L.M.; writing—original draft, K.S. and J.L.M.; writing—review and editing, K.S. and J.L.M. Both authors have read and agreed to the published version of the manuscript.

Funding

Research reported in this publication was supported by the National Institute Of Biomedical Imaging And Bioengineering of the National Institutes of Health under Award Number R01EB026710. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Colorado Multiple Institutional Review Board (COMIRB) (approval number COMIRB 14-0652).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Costa, E.L.V.; Lima, R.G.; Amato, M.B.P. Electrical impedance tomography. In Intensive Care Medicine; Vincent, J.L., Ed.; Springer: New York, NY, USA, 2009; pp. 394–404. [Google Scholar]
  2. Nguyen, D.T.; Thiagalingam, J.C.; McEwan, A.L.A. A review on electrical impedance tomography for pulmonary perfusion imaging. Physiol. Meas. 2012, 33, 695–706. [Google Scholar] [CrossRef]
  3. Arad, M.; Zlochiver, S.; Davidson, T.; Shoenfeld, Y.; Adunsky, A.; Abboud, A. The detection of pleural effusion using a parametric eit technique. Physiol. Meas. 2009, 30, 421–428. [Google Scholar] [CrossRef]
  4. Costa, E.L.; Chaves, C.N.; Gomes, S.; Beraldo, M.A.; Volpe, M.S.; Tucci, M.R.; Schettino, I.A.; Bohm, S.H.; Carvalho, C.R.; Tanaka, H.; et al. Real-time detection of pneumothorax using electrical impedance tomography. Crit. Care Med. 2008, 36, 1230–1238. [Google Scholar] [CrossRef]
  5. Frerichs, I.; Pulletz, S.; Elke, G.; Reifferscheid, F.; Schädler, D.; Scholz, J.; Weiler, N. Assessment of changes in distribution of lung perfusion by electrical impedance tomography. Respiration 2009, 77, 282–291. [Google Scholar] [CrossRef] [PubMed]
  6. Lowhagen, K.; Lundin, S.; Stenqvist, O. Regional intratidal gas distribution in acute lung injury and acute respiratory distress syndrome—Assessed by electric impedance tomography. Minerva Anestesiol. 2010, 76, 1024–1035. [Google Scholar]
  7. Muders, T.; Luepschen, H.; Putensen, C. Impedance tomography as a new monitoring technique. Curr. Opin. Crit. Care 2010, 16, 269–275. [Google Scholar] [CrossRef]
  8. Reinius, H.; Borges, J.B.; Fredén, F.; Jideus, L.; Camargo, E.D.; Amato, M.B.; Hedenstierna, G.; Larsson, A.; Lennmyr, F. Real-time ventilation and perfusion distributions by electrical impedance tomography during one-lung ventilation with capnothorax. Acta Anaesthesiol. Scand. 2015, 59, 354–368. [Google Scholar] [CrossRef] [PubMed]
  9. Victorino, J.A.; Borges, J.B.; Okamoto, V.N.; Matos, G.F.J.; Tucci, M.R.; Caramez, M.P.R.; Tanaka, H.; Sipmann, F.S.; Santos, D.C.B.; Barbas, C.S.V.; et al. Imbalances in regional lung ventilation: A validation study on electrical impedance tomography. Am. J. Respir. Crit. Care Med. 2004, 169, 791–800. [Google Scholar] [CrossRef] [PubMed]
  10. Calderón, A.P. On an inverse boundary value problem. In Seminar on Numerical Analysis and Its Applications to Continuum Physics; Sociedade Brasileira de Matemàtica: Rio de Janeiro, Brazil, 1980; pp. 65–73. [Google Scholar]
  11. Knudsen, K.; Lassas, M.; Mueller, J.L.; Siltanen, S. D-Bar method for electrical impedance tomography with discontinuous conductivities. SIAM J. Appl. Math. 2007, 7, 893–913. [Google Scholar] [CrossRef] [Green Version]
  12. Shin, K.; Mueller, J. A second order Calderón’s method with a correction term and a priori information. Inverse Probl. 2020, 32, 124005. [Google Scholar] [CrossRef]
  13. Nachman, A.I. Global uniqueness of a two-dimensional inverse boundary value problem. Ann. Math. 1966, 2, 71–96. [Google Scholar] [CrossRef]
  14. Bikowski, J.; Mueller, J.L. 2D EIT reconstructions using Calderón’s method. Inverse Probl. Imaging 2008, 2, 43. [Google Scholar] [CrossRef]
  15. Muller, P.A.; Isaacson, D.; Newell, J.C.; Saulnier, G.J. Calderón’s method on an elliptical domain. Physiol. Meas. 2013, 34, 609–622. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Muller, P.A.; Mueller, J.L.; Mellenthin, M.M. Real-Time Implementation of Calderón’s Method on Subject-Specific Domains. IEEE Trans. Med. Imaging 2017, 36, 1868–1875. [Google Scholar] [CrossRef]
  17. Muller, P.A.; Mueller, J.L. Reconstruction of complex conductivities by calderon’s method on subject-specific domains. In Proceedings of the 2018 International Applied Computational Electromagnetics Society Symposium (ACES), Denver, CO, USA, 25–29 March 2018. [Google Scholar]
  18. Choi, M.H.; Kao, T.; Isaacson, D.; Saulnier, G.J.; Newell, J.C. A reconstruction Algorithm for breast cancer imaging with electrical impedance tomography in mammography geometry. IEEE Trans. Biomed. Eng. 2007, 54, 700–710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Shin, K.; Ahmad, S.; Mueller, J.L. Three dimensional Calderón’s method for EIT on the cylindrical geometry. IEEE Trans. Biomed. Eng. 2021, 68, 1487–1495. [Google Scholar] [CrossRef]
  20. Dobson, D.C.; Santosa, F. An image-enhancement technique for electrical impedance tomography. Inverse Probl. 1994, 10, 317–334. [Google Scholar] [CrossRef]
  21. Camargo, E.D.L.B. Development of an Absolute Electrical Impedance Imaging Algorithm for Clinical Use; University of São Paulo: São Paulo, Brazil, 2013. [Google Scholar]
  22. Ferrario, D.; Grychtol, B.; Adler, A.; Sola, J.; Bohm, S.H.; Bodenstein, M. Toward Morphological Thoracic EIT: Major Signal Sources Correspond to Respective Organ Locations in CT. IEEE Trans. Med. Imaging 2012, 59, 3000–3008. [Google Scholar] [CrossRef]
  23. Flores-Tapia, D.; Pistorius, S. Electrical impedance tomography reconstruction using a monotonicity approach based on a priori knowledge. In Proceedings of the 2010 Annual International Conference of the IEEE, Buenos Aires, Argentina, 31 August–4 September 2010; pp. 4996–4999. [Google Scholar]
  24. Dehghani, H.; Barber, D.C.; Basarab-Horwath, I. Incorporating a priori anatomical information into image reconstruction in electrical impedance tomography. Physiol. Meas. 1999, 20, 87–102. [Google Scholar] [CrossRef] [Green Version]
  25. Kaipio, J.P.; Kolehmainen, V.; Vauhkonen, M.; Somersalo, E. Inverse problems with structural prior information. Inverse Probl. 1999, 15, 713. [Google Scholar] [CrossRef]
  26. Vauhkonen, M.; Vadasz, D.; Karjalainen, P.A.; Somersalo, E.; Kaipio, J.P. Tikhonov regularization and prior information in electrical impedance tomography. IEEE Trans. Med. Imaging 1998, 17, 285–293. [Google Scholar] [CrossRef] [PubMed]
  27. Avis, N.J.; Barber, D.C. Incorporating a priori information into the Sheffield filtered backprojection algorithm. Physiol. Meas. 1995, 16, A111–A122. [Google Scholar] [CrossRef] [PubMed]
  28. Soleimani, M. Electrical impedance tomography imaging using a priori ultrasound data. BioMed. Eng. Online 2006, 5, 1–8. [Google Scholar] [CrossRef] [PubMed]
  29. Baysal, U.; Eyüboglu, B.M. Use of a priori information in estimating tissue resistivities—A simulation study. Phys. Med. Biol. 1998, 43, 3589–3606. [Google Scholar] [CrossRef] [PubMed]
  30. Alsaker, M. Computational Advancements in the D-Bar Reconstruction Method for 2-D Electrical Impedance Tomography. Ph.D. Dissertation, Colorado State University, Fort Collins, CO, USA, 2016. [Google Scholar]
  31. Alsaker, M.; Hamilton, S.J.; Hauptmann, A. A direct D-bar method for partial boundary data electrical impedance tomography with a prior information. Inverse Probl. Imaging 2017, 11, 427–454. [Google Scholar] [CrossRef]
  32. Alsaker, M.; Mueller, J.L. A D-bar algorithm with a priori information for 2 dimensional electrical impedance tomography. SIAM J. Imaging Sci. 2016, 9, 1619–1654. [Google Scholar] [CrossRef]
  33. Alsaker, M.; Mueller, J.L. Use of an optimized spatial prior in D-bar reconstructions of EIT tank data. Inverse Probl. Imaging 2018, 12, 883–901. [Google Scholar] [CrossRef] [Green Version]
  34. Santos, T.; Nakanishi, R.M.; Kaipio, J.P.; Mueller, J.L.; Lima, R.G. Introduction of sample based prior into the D-Bar method through a Schur complement property. IEEE Trans. Med. Imaging 2020, 39, 4085–4093. [Google Scholar] [CrossRef]
  35. Alsaker, M.; Mueller, J.L.; Murthy, R. Dynamic optimized priors for D-bar reconstructions of human ventilation using electrical impedance tomography. J. Comput. Appl. Math. 2019, 362, 276–294. [Google Scholar] [CrossRef]
  36. Ren, S.; Sun, K.; Liu, D.; Dong, F. A Statistical Shape-Constrained Reconstruction Framework for Electrical Impedance Tomography. IEEE Trans. Med. Imaging 2019, 38, 2400–2410. [Google Scholar] [CrossRef]
  37. Liu, D.; Gu, D.; Smyl, D.; Khambampati, A.K.; Deng, J.; Du, J. Shape-Driven EIT Reconstruction Using Fourier Representations. IEEE Trans. Med. Imaging 2021, 40, 481–490. [Google Scholar] [CrossRef]
  38. Liu, D.; Gu, D.; Smyl, D.; Deng, J.; Du, J. B-Spline Level Set Method for Shape Reconstruction in Electrical Impedance Tomography. IEEE Trans. Med. Imaging 2020, 39, 1917–1929. [Google Scholar] [CrossRef]
  39. Liu, D.; Du, J. A Moving Morphable Components Based Shape Reconstruction Framework for Electrical Impedance Tomography. IEEE Trans. Med. Imaging 2019, 38, 2937–2948. [Google Scholar] [CrossRef] [PubMed]
  40. Evans, L.C. Partial Differential Equations; American Mathematical Society: Providence, RI, USA, 2010. [Google Scholar]
  41. Mellenthin, M.M.; Mueller, J.L.; de Camargo, E.D.L.B.; Moura, F.S.D.; Santos, T.B.R.; Lima, R.G.; Hamilton, S.J.; Muller, P.A.; Alsaker, M. The ACE1 Electrical Impedance Tomography System for Thoracic Imaging. IEEE Trans. Instrum. Meas. 2019, 68, 3137–3150. [Google Scholar] [CrossRef]
  42. Mellenthin, M.M.; Meuller, J.L.; de Camargo, E.D.L.B.; Moura, F.S.D.; Himilton, S.J.; Lima, R.G. The ACE1 thoracic Electrical Impedance Tomography system for ventilation and perfusion. In Proceedings of the 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Milan, Italy, 25–29 August 2015. [Google Scholar]
  43. Alsaker, M.; Cárdenas, D.A.C.; Furuie, S.S.; Mueller, J.L. Complementary use of priors for pulmonary imaging with electrical impedance and ultrasound computed tomography. J. Comput. Appl. Math. 2021, 395, 113591. [Google Scholar] [CrossRef] [PubMed]
  44. Hamilton, S.J.; Mueller, J.L. Direct EIT reconstructions of complex admittivities on a chest-shaped domain in 2-D. IEEE Trans. Med. Imaging 2013, 32, 757–769, Epub 9 January 2013. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Top left: The chest shaped tank filled with saline bath and inclusions simulating two low conductive lungs and the high conductive heart. Top right: A difference image without modeling the location of electrodes. Bottom left: A difference image with modeling the location of electrodes. Bottom right: An absolute image with modeling the location of electrodes. The L2 norms of the differences between those reconstructions and the ground truth are 0.1128, 0.0973 and 0.0981 in order.
Figure 1. Top left: The chest shaped tank filled with saline bath and inclusions simulating two low conductive lungs and the high conductive heart. Top right: A difference image without modeling the location of electrodes. Bottom left: A difference image with modeling the location of electrodes. Bottom right: An absolute image with modeling the location of electrodes. The L2 norms of the differences between those reconstructions and the ground truth are 0.1128, 0.0973 and 0.0981 in order.
Sensors 21 05635 g001
Figure 2. Healthy heart and lungs. From left: Reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ diff with the spatial prior. The L2 norms of the difference between those reconstruction and the ground truth, the first figure, are 0.2952, 0.3119, 0.2920 and 0.3085 in order.
Figure 2. Healthy heart and lungs. From left: Reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ diff with the spatial prior. The L2 norms of the difference between those reconstruction and the ground truth, the first figure, are 0.2952, 0.3119, 0.2920 and 0.3085 in order.
Sensors 21 05635 g002
Figure 3. Small pneumothorax. From top left: The true conductivity distribution, reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ pr diff with the spatial prior, a difference image in which the data for the example of healthy heart and lungs in Figure 2 with no pathology is chosen for γ ref . The L2 norms of the differences between those reconstructions and the ground truth are 0.3038, 0.3137, 0.3013 and 0.3112 in order.
Figure 3. Small pneumothorax. From top left: The true conductivity distribution, reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ pr diff with the spatial prior, a difference image in which the data for the example of healthy heart and lungs in Figure 2 with no pathology is chosen for γ ref . The L2 norms of the differences between those reconstructions and the ground truth are 0.3038, 0.3137, 0.3013 and 0.3112 in order.
Sensors 21 05635 g003
Figure 4. Large pneumothorax. From top left: The true conductivity distribution, reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ pr diff with the spatial prior, a difference image in which the data for the example of healthy heart and lungs in Figure 2 with no pathology is chosen for γ ref . The L2 norms of the differences between those reconstructions and the ground truth are 0.3076, 0.3261, 0.3047 and 0.3232 in order.
Figure 4. Large pneumothorax. From top left: The true conductivity distribution, reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ pr diff with the spatial prior, a difference image in which the data for the example of healthy heart and lungs in Figure 2 with no pathology is chosen for γ ref . The L2 norms of the differences between those reconstructions and the ground truth are 0.3076, 0.3261, 0.3047 and 0.3232 in order.
Sensors 21 05635 g004
Figure 5. Simulated pulmonary contusion. From left: The true conductivity distribution, reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ pr diff with the spatial prior, a difference image in which the data for the example of healthy heart and lungs in Figure 2 with no pathology is chosen for γ ref . The L2 norms of the differences between those reconstructions and the ground truth are 0.2958, 0.3071, 0.2927 and 0.3040 in order.
Figure 5. Simulated pulmonary contusion. From left: The true conductivity distribution, reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ pr diff with the spatial prior, a difference image in which the data for the example of healthy heart and lungs in Figure 2 with no pathology is chosen for γ ref . The L2 norms of the differences between those reconstructions and the ground truth are 0.2958, 0.3071, 0.2927 and 0.3040 in order.
Sensors 21 05635 g005
Figure 6. Simulated pleural effusion. From left: The true conductivity distribution, reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ pr diff with the spatial prior, a difference image in which the data for the example of healthy heart and lungs in Figure 2 with no pathology is chosen for γ ref . The L2 norms of the differences between those reconstructions and the ground truth are 0.3000, 0.3065, 0.2966 and 0.3029 in order.
Figure 6. Simulated pleural effusion. From left: The true conductivity distribution, reconstruction of γ abs with no prior, γ diff with no prior, γ pr abs with the spatial prior, γ pr diff with the spatial prior, a difference image in which the data for the example of healthy heart and lungs in Figure 2 with no pathology is chosen for γ ref . The L2 norms of the differences between those reconstructions and the ground truth are 0.3000, 0.3065, 0.2966 and 0.3029 in order.
Sensors 21 05635 g006
Figure 7. One slice of the CT scan of the subject and the static prior constructed from the CT image.
Figure 7. One slice of the CT scan of the subject and the static prior constructed from the CT image.
Sensors 21 05635 g007
Figure 8. Difference images of ventilation data with no spatial prior included in the algorithm. Frame 200 is used for the reference image. The truncation radius is 1.2.
Figure 8. Difference images of ventilation data with no spatial prior included in the algorithm. Frame 200 is used for the reference image. The truncation radius is 1.2.
Sensors 21 05635 g008
Figure 9. Difference images of ventilation data with the spatial prior included in the algorithm. Frame 200 is used for the reference image. The truncation radius is 1.2 and the outer truncation radius is 3.
Figure 9. Difference images of ventilation data with the spatial prior included in the algorithm. Frame 200 is used for the reference image. The truncation radius is 1.2 and the outer truncation radius is 3.
Sensors 21 05635 g009
Figure 10. The time trace plots of a pixel in the right lung region and a pixel in the heart region from the reconstructions. Two pixels are marked with squares in the first image of Figure 12. Left is the plot from the ventilation data and the right one is for the breath-holding data. In both plots, from the top to the bottom, each graph represents the time trace of the pixel in the heart with and without a prior, and the time trace of the right lung pixel with and without a prior. The plots of the lung pixels in the perfusion time trace plot are shifted upward by 0.02 in order to better demonstrate the out-of-phase nature.
Figure 10. The time trace plots of a pixel in the right lung region and a pixel in the heart region from the reconstructions. Two pixels are marked with squares in the first image of Figure 12. Left is the plot from the ventilation data and the right one is for the breath-holding data. In both plots, from the top to the bottom, each graph represents the time trace of the pixel in the heart with and without a prior, and the time trace of the right lung pixel with and without a prior. The plots of the lung pixels in the perfusion time trace plot are shifted upward by 0.02 in order to better demonstrate the out-of-phase nature.
Sensors 21 05635 g010
Figure 11. Difference images of breath holding data with no spatial prior included in the algorithm. Averaged voltage data were used for the reference data. The truncation radius is 1.2.
Figure 11. Difference images of breath holding data with no spatial prior included in the algorithm. Averaged voltage data were used for the reference data. The truncation radius is 1.2.
Sensors 21 05635 g011
Figure 12. Difference images of breath holding data with the spatial prior included in the algorithm. Averaged voltage data were used for the reference data. The inner truncation radius is 1.2 and the outer truncation radius is 3. The small squares superimposed on the upper left figure indicate the pixels chosen for the time traces in Figure 10.
Figure 12. Difference images of breath holding data with the spatial prior included in the algorithm. Averaged voltage data were used for the reference data. The inner truncation radius is 1.2 and the outer truncation radius is 3. The small squares superimposed on the upper left figure indicate the pixels chosen for the time traces in Figure 10.
Sensors 21 05635 g012
Figure 13. EKG data collected with BIOPAC during breath-holding. Frames from Figure 11 and Figure 12 are annotated on the EKG trace in the upper half of the image. One sees that systole in the EIT images closely coincides with the QRS complex in the EKG data, as would be expected. The average heart rate is found in the lower half of the figure.
Figure 13. EKG data collected with BIOPAC during breath-holding. Frames from Figure 11 and Figure 12 are annotated on the EKG trace in the upper half of the image. One sees that systole in the EIT images closely coincides with the QRS complex in the EKG data, as would be expected. The average heart rate is found in the lower half of the figure.
Sensors 21 05635 g013
Figure 14. Difference images of ventilation data with modeling the domain to be a circle. No prior is added. Frame 200 is used for the reference image. The truncation radius is 1.2.
Figure 14. Difference images of ventilation data with modeling the domain to be a circle. No prior is added. Frame 200 is used for the reference image. The truncation radius is 1.2.
Sensors 21 05635 g014
Figure 15. Difference images of ventilation data with modeling the domain to be a circle. A spatial prior is added. Frame 200 is used for the reference image. The truncation radius is 1.2.
Figure 15. Difference images of ventilation data with modeling the domain to be a circle. A spatial prior is added. Frame 200 is used for the reference image. The truncation radius is 1.2.
Sensors 21 05635 g015
Table 1. Conductivity values for the simulated data by organ.
Table 1. Conductivity values for the simulated data by organ.
Organ/RegionConductivity (S/m)True δ abs Value
Background0.80
Heart1.10.3
Lung0.5−0.3
Spine0.15 −0.65
Pneumothorax0.25 −0.55
Contusion 0.6 −0.2
Effusion 1.0 0.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shin, K.; Mueller, J.L. Calderón’s Method with a Spatial Prior for 2-D EIT Imaging of Ventilation and Perfusion. Sensors 2021, 21, 5635. https://doi.org/10.3390/s21165635

AMA Style

Shin K, Mueller JL. Calderón’s Method with a Spatial Prior for 2-D EIT Imaging of Ventilation and Perfusion. Sensors. 2021; 21(16):5635. https://doi.org/10.3390/s21165635

Chicago/Turabian Style

Shin, Kwancheol, and Jennifer L. Mueller. 2021. "Calderón’s Method with a Spatial Prior for 2-D EIT Imaging of Ventilation and Perfusion" Sensors 21, no. 16: 5635. https://doi.org/10.3390/s21165635

APA Style

Shin, K., & Mueller, J. L. (2021). Calderón’s Method with a Spatial Prior for 2-D EIT Imaging of Ventilation and Perfusion. Sensors, 21(16), 5635. https://doi.org/10.3390/s21165635

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop