Next Article in Journal
Influence of the Porosity of Polymer Foams on the Performances of Capacitive Flexible Pressure Sensors
Next Article in Special Issue
Adaptive Selection of Truncation Radius in Calderon’s Method for Direct Image Reconstruction in Electrical Capacitance Tomography
Previous Article in Journal
A Unified Multiple-Target Positioning Framework for Intelligent Connected Vehicles
Previous Article in Special Issue
Comparison of Selected Machine Learning Algorithms for Industrial Electrical Tomography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Lagrange-Newton Method for EIT/UT Dual-Modality Image Reconstruction †

1
Tianjin Key Laboratory of Process Measurement and Control, School of Electrical and Information Engineering, Tianjin University, Tianjin 300072, China
2
Institute of Biomedical Engineering, CAMS & PUMC (Chinese Academy of Medical Sciences and Peking Union Medical College), Tianjin 300192, China
*
Authors to whom correspondence should be addressed.
This paper is an extended version of our paper published in Liang, G.; Ren, S.; Dong, F. An Inclusion Boundary Reconstruction Method Using Electrical Impedance and Ultrasound Reflection Dual-Modality Tomography. In Proceedings of the 9th World Congress on Industrial Process Tomography, Bath, UK, 2–6 September 2018.
Sensors 2019, 19(9), 1966; https://doi.org/10.3390/s19091966
Submission received: 30 March 2019 / Revised: 23 April 2019 / Accepted: 23 April 2019 / Published: 26 April 2019

Abstract

:
An image reconstruction method is proposed based on Lagrange-Newton method for electrical impedance tomography (EIT) and ultrasound tomography (UT) dual-modality imaging. Since the change in conductivity distribution is usually accompanied with the change in acoustic impedance distribution, the reconstruction targets of EIT and UT are unified to the conductivity difference using the same mesh model. Some background medium distribution information obtained from ultrasound transmission and reflection measurements can be used to construct a hard constraint about the conductivity difference distribution. Then, the EIT/UT dual-modality inverse problem is constructed by an equality constraint equation, and the Lagrange multiplier method combining Newton-Raphson iteration is used to solve the EIT/UT dual-modality inverse problem. The numerical and experimental results show that the proposed dual-modality image reconstruction method has a better performance than the single-modality EIT method and is more robust to the measurement noise.

1. Introduction

Electrical impedance tomography (EIT) is a non-invasive technique to reconstruct the distribution of the media within a closed vessel based on the differences of their conductivities. Due to its advantages of high speed, non-radiation, and low cost, EIT has been widely researched and used in many fields, such as, industrial process imaging [1,2], medical imaging [3,4], geological exploration [5,6], structural health monitoring [7,8], and nondestructive evaluation [9,10]. However, as a soft-field imaging technique, EIT usually suffers from the low spatial resolution and imaging accuracy for its inherent nonlinearity and ill-posedness [11]. It has hindered the further application of EIT in many fields, especially for medical imaging, where a higher imaging accuracy is usually expected [12]. Thus, the study of improving the EIT spatial resolution and imaging accuracy can always attract the attention of the researchers.
Many researchers have devoted to the study of improving the spatial resolution and imaging accuracy of EIT and many methods have been developed. For the nonlinearity of EIT, the directly linearization [13], iterative linearization [14], and directly nonlinear methods [15] have been proposed to formulate the EIT inverse problem. For the ill-posedness of EIT, the regularization methods, such as, Tikhonov regularization [14], Total Variation (TV) regularization [16], and sparse regularization [17], have been proposed to constrain the solution space of EIT inverse problem by introducing the prior information about the targets to be reconstructed. Thus, the construction of the prior information about the target is important for the EIT inverse problem. The prior information is usually predetermined and added into the EIT inverse problem. A commonly used method for constructing the prior information is to use the characteristics of the target itself. For example, in the imaging problems of small targets, the sparse prior constraints are commonly used [17]. For specific target imaging problems, the statistical modeling methods can be used to extract the similarity information of the targets, such as structure or shape, and construct the statistical prior information about the targets. For example, in EIT lung function monitoring, the statistical shape prior extracted from computed tomography (CT) images can be used to alleviate the ill-posedness of EIT inverse problem and improve the EIT imaging accuracy and measurement noise robustness [18].
Although the spatial resolution and imaging accuracy have been improved to some extent by the regularization constraint from the prior information in EIT inverse problem, the acquisition of the prior information about the target usually falls into the problems of lack of the data sources and feature extraction difficulties. In addition, the individual difference of the imaging targets can also invalidate the prior information in some practical applications. Thus, relying on a single-modality imaging technique is difficult to substantially improve the EIT spatial resolution and imaging accuracy because the available prior information is usually limited. Considering the potential information complementarity between different imaging modalities, the dual-modality imaging techniques have attracted more and more attention in the past decade. For example, the positron emission computed tomography (PET)/CT [19], magnetic resonance imaging (MRI)/CT [20] and Ultrasound/MRI [21] have been successfully used in medical diagnosis. The ultrasound tomography (UT)/electrical capacitance tomography (ECT) [22], CT/ECT [23], Gamma Densitometry Tomography (GDT)/EIT [24] have been widely used in multiphase flow process parameters detection. Thus, the dual-modality imaging technique is promising to improve the EIT imaging accuracy by adding the prior information from the other imaging modality.
UT is also a non-invasive technique to reconstruct the media distribution based on the differences of their acoustic impedances and has been widely used in industrial imaging [25] and medical imaging [26]. Different from the EIT, UT is a ‘hard-field’ imaging technique and the spatial resolution and imaging accuracy of UT is mainly depended on the amount of measured data available. The more measured data means that more ultrasonic transducers are needed in UT [27]. However, the number of the ultrasonic transducer is usually not enough compared to the need for the imaging accuracy in practical applications. Fortunately, UT can obtain some accurate local media distribution information even the number of the ultrasonic transducer is limited [28,29]. Thus, the prior information about the media distribution from UT is promising to improve the imaging accuracy of EIT.
Some studies have been reported about the EIT/UT dual-modality imaging in the past decade. In 2008, Steiner et al. proposed an EIT/ultrasound reflection tomography (URT) dual-modality imaging method for small object detection in biomedical applications, where the measurement information from URT was used to construct the weighted regularization matrix in EIT inverse problem [30]. In 2010, Wan et al. presented a transrectal ultrasound coupled EIT imaging method for transrectal cancer detection, where the contour of the target was captured by ultrasound measurement, and then used to guide the EIT image reconstruction [31]. In 2015, Teniou et al. proposed an electrical resistance tomography (ERT)/URT dual-modality imaging method to detect the small lesions in biological soft tissue, where the measurements from URT was used as hard constraints during the EIT image reconstruction process [32]. In 2017, Pusppanathan et al. proposed an ECT/ultrasound transmission tomography (UTT) dual-modality imaging method for multiphase flow imaging, where the images from the ECT and UTT are fused by the pixel-based fuzzy logic method, and then used to distinguish oil/gas/water three-phase media [33]. In 2018, Liang et al. proposed a shape-based EIT/URT dual-modality imaging method, where some accurate boundaries information about the inclusions detected from URT are used to improve the EIT inclusion boundary reconstruction accuracy [34]. In addition, the emerging impediography can improve the EIT image accuracy by coupling the mechanical effects driven from ultrasound waves [35,36,37]. Although the EIT/UT dual-modality imaging has made some progresses, the way of information fusion between EIT and UT is still an open issue in EIT/UT dual-modality imaging researches.
In this paper, a dual-modality image reconstruction method combining EIT and UT is proposed. Some prior information about the background medium distribution is determined by the ultrasound reflection and transmission measurements, and then used as a hard constraint added into the EIT inverse problem through an equality constraint equation. Then, the ultrasound constraint EIT image reconstruction problem is reformed by the Lagrange multiplier method and iteratively solved with Newton-Raphson method. To test the performance of the proposed dual-modality imaging method, the numerical and experimental tests are conducted.

2. Mathematical Model of EIT and UT

2.1. EIT Mathematical Model

EIT mathematical model is derived from Maxwell’s electromagnetic theory under the condition of low frequencies. In the quasi-static approximation, the electrical fields can be described in terms of a scalar voltage potential u satisfying the following equation:
· σ u = 0 in Ω ,
where σ represents the conductivity in the imaging domain Ω.
Considering the contact impedance between the background medium and electrodes, the complete electrode model (CEM) is used to define the boundary conditions [38]:
{ σ u · n = 0 , on Ω \ l = 1 L e l e l σ u · n = I l , u + z l σ u · n = U l , on e l ,
where ∂Ω is the boundary of the imaging domain Ω, e l Ω represents the location covered by the l-th electrode, L is the total number of the electrodes, n is the outward normal vector on ∂Ω, z l is the contact impedance between the l-th electrode and background medium, I l and U l are the electric current and electric potential at the l-th electrode, respectively. Subject to the law of the conservation of charge, the electric current and voltage in the electrodes satisfy the following equation:
l = 1 L I l = 0 and l = 1 L U l = 0 ,
Thus, with the definition of EIT mathematical model and the corresponding boundary conditions, the electric potential, conductivity, electric current, and contact impedance can be formulated as the following nonlinear equation:
U = u ( σ , z , I ) ,
Considering a small change on the conductivity, the corresponding change on the boundary voltage can be linearly represented using Equation (4) as follows:
Δ U = u ( σ , z , I ) u ( σ 0 , z , I ) = d u ( σ , z , I ) d σ | σ = σ 0 ( σ σ 0 ) + O ( ( σ σ 0 ) 2 ) = u ( σ 0 , z , I ) Δ σ + O ( ( Δ σ ) 2 ) ,
If the change of conductivity Δ σ is small enough, the high order terms O ( ( Δ σ ) 2 ) can be neglected, and Equation (5) will be simplified as follows:
Δ U = u ( σ 0 , z , I ) Δ σ ,
further, the Equation (6) is usually discretized as follows:
y = J x ,
where y R M × 1 represents the boundary voltage difference, x R N × 1 represents the unknown conductivity difference, J R M × N represents the Jacobian matrix, which can be calculated using the Geselowitz’s sensitivity theorem [39]:
J m , n = Ω n u ( I i ) I i u ( I j ) I j d S ,
where J m , n is the element at the m-th row and n column of the Jacobian matrix, u ( I i ) is the electric potential for the i-th electrode pair with electric current I i , u ( I j ) is the electric potential for the j-th electrode pair with electric current I j , Ω n is the region covered by n-th pixel.
For a given conductivity difference x, the boundary voltage y can be determined by the Equations (1)–(3). This process is usually mentioned as the EIT forward problem, where the Boundary Element Method (BEM) solver developed by Ren et al. is used to calculate the EIT forward problem and obtain the sensitivity matrix by the perturbation method [40,41]. In contrast, given the boundary voltage y, the EIT inverse problem is to calculate the conductivity difference x by the Equation (7). Since the sensitivity matrix J is usually singular, the unknown conductivity difference x cannot be directly obtained from the Equation (7) in EIT inverse problem. The regularization methods are commonly used to solve it, which will be introduced in the latter section.

2.2. UT Mathematical Model

When the influence of medium heterogeneity on acoustic wave propagation is not considered, the propagation of ultrasonic wave in medium can be approximated as linear propagation model. This model is usually referred as the geometrical/ray acoustic model and widely used in UT image reconstruction. Thus, the mathematical model of UT can be referred to the geometrical optics or ray optics model, and the most commonly considered modes of the ultrasonic wave propagation are transmission and reflection modes, as shown in Figure 1. In addition, the piston transducer is used in the paper, where the energy of the ultrasonic wave is mainly concentrated in a narrow main lobe area [42]. The power of the transmission and reflection wave is mainly affected by the acoustic impedance of the media where the ultrasonic wave pass through, and it can be described as follows:
α R = ( P r P i ) 2 = ( Z 2 Z 1 Z 2 + Z 1 ) 2 ,
α T = ( P t P i ) 2 = ( 2 Z 2 Z 2 + Z 1 ) 2 ,
where P i ,   P r ,   P t represent incident, reflection, and transmission sound pressures, respectively. α R and α T represent the reflection and transmission coefficients, respectively. Z 1 and Z 2 are the acoustic impedance of the media and equal to the product of medium density ρ and sound velocity c.

2.2.1. Transmission Mode Mathematical Model

Under the assumption that the ultrasonic wave is propagating along the straight line, the mathematical model of transmission mode UT is similar to the x-ray CT. The ultrasonic waves sending from the exciting transducer will pass through the media and arrive at the measuring transducers. Consider that the acoustic attenuation effect of the medium, the transmission mode UT usually utilizes the energy attenuation of the ultrasonic waves to invert the medium distribution in the ultrasonic wave propagation paths. Inspired from the Equation (10), we reformulate the sending acoustic energy I e from the exciting transducer and the receiving acoustic energy I m from the measuring transducer as the following equation:
Ω k μ d S = ln I m I e ,
where μ is the medium’s attenuation coefficient, Ω k is the region where the k-th ultrasonic wave pass through.
Consider a scenario that medium located in the region Ω k changes, the corresponding attenuation coefficient μ in Ω k will be also changing to μ + δ μ . Then, the energy attenuation of the ultrasonic wave in the region Ω k will change and the Equation (11) can be rewritten as follows:
Ω k ( δ μ ) d S = Ω k ( μ + δ μ ) d S Ω k μ d S = ln I m a I e ( ln I m b I e ) = ln I m b I m a ,
where δ μ is the variation of the attenuation coefficient in Ω k , I m b and I m a represent the acoustic energy before and after the attenuation coefficient changes, respectively.

2.2.2. Reflection Mode Mathematical Model

Similar to the propagation of optics ray in the medium, the mathematical model of the reflection mode UT can be derived from the Fermat’s principle, where the ultrasonic waves will reflect at the media interface and continue to propagate to the measuring transducer along the shortest propagation path. Thus, the ultrasonic waves sending from the exciting transducer will be reflected at the interface where the acoustic impedances change and continue to propagate in the background medium to the measuring transducer. When the exciting and measuring transducer are determined, the detected points will fall on an elliptical arc. This elliptical arc can be determined by the coordinates of the exciting and measuring transducer and the time of flight (TOF) of the ultrasonic wave, which can be described as follows:
1 c · [ ( x d x e ) 2 + ( y d y e ) 2 + ( x d x m ) 2 + ( y d y m ) 2 ] = t f ,
where c is the sound speed in the background medium, t f is the TOF of the ultrasonic wave in the background medium, ( x d , y d ) is the coordinate of the detected point, which obeys the arc distribution, ( x e , y e ) and ( x m , y m ) are the coordinates of the exciting and measuring transducer, respectively.

2.3. Complementarity of EIT and URT Sensitive Field

The boundary measurements of EIT and UT have different sensitivities about the changes of the sensitive parameters in different spaces of the imaging domain, and higher sensitivity means higher detection resolution and accuracy. In EIT, the electric current is injected into the imaging domain through the electrodes attached at the boundary. The electric field lines near the edge of the imaging domain are dense than the central region of the imaging domain. Therefore, the EIT measurements have a higher sensitivity to the changes in conductivity near the edge of the imaging domain. According to the reciprocity principle, the sensitivity of EIT measurements about the changes of conductivity is proportional to the integration of the electric field intensity [43]. As shown in Figure 2a, EIT has higher sensitivity near the edge of the imaging domain. In UT, the ultrasonic waves emitted by the transducers are mainly concentrated in a narrow main lobe range along the radial of the transducer. Only the changes of acoustic impedance occurring at the propagation path of ultrasonic waves will cause the change on reflection or transmission measurements. In the imaging domain, all the ultrasonic waves propagation paths of the transducers are intersected in the central of the imaging domain. Therefore, the UT measurements have a higher sensitivity about the changes in acoustic impedance near the central of the imaging domain [44]. According to the radon transform theory, the sensitivity of UT measurements about the change of acoustic impedance can be formed by the superimposition of all single transducer ultrasonic waves propagation paths. As shown in Figure 2b, UT has higher sensitivity near the central of the imaging domain. Based on the complementarity of EIT and UT sensitivity space distribution, the combination of EIT and UT is expected to improve the imaging resolution and accuracy.

3. EIT/UT Dual-Modality Image Reconstruction

3.1. Unification of Mesh Model and Reconstruction Target

In the paper, the purpose of both EIT and UT is to reconstruct the media distribution based on the change of conductivity or acoustic impedance distribution. In some applications, when the media distribution changes, the conductivity and acoustic impedance will also change. For example, in gas/liquid two-phase flow cross-section reconstruction problem, the presence of bubble will change the distribution of conductivity and acoustic impedance, the targets of EIT and UT are all the reconstruction of the bubble distribution. Thus, the reconstruction target of EIT and UT can be uniformed to the changes of the media distribution, which means that the background medium has the same representation in EIT and UT. In addition, another factor to consider is the mesh model. Based on this demand of reconstructing the same unknown media distribution using EIT and UT, the mesh model for EIT and UT should be also uniformed in the same imaging region. To achieve the purpose, the same mesh model is used in EIT and UT, as shown in Figure 3.

3.2. Total Variation Regularization for EIT

It’s known that EIT inverse problem is a typical ill-posed inverse problem due to the large condition number of the Jacobian matrix J, which means that it’s unrealistic to directly solve the unknown x from the Equation (7). To overcome this problem, the regularization methods are usually adopted. First, we convert the EIT inverse problem to find the minimum from the least-square equation J x y 2 2 . Then, a penalty term, also referred as regularization term, is added to stabilize the solution process by constraining the solution space. This process can be formulated in the following equations:
min x R N f ( x ) ,
f ( x ) = 1 2 J x y 2 2 + β R ( x ) ,
where the regularization parameter β is a positive constant and balances the residual term J x y 2 2 and the penalty term R ( x ) in the inverse problem solving process. There are many forms of the regularization term based on the difference of the solution space constraints. Here, the Total Variation (TV) regularization term is adopted, and the regularization term can be written as:
R ( x ) = Ω | x | d S ,
However, the Equations (14)–(16) cannot be directly solved due to the non-differentiability of | x | . To overcome this problem, a smooth approximation is usually adopted [45]. Then, the regularization term R ( x ) can be rewritten as:
R ( x ) = Ω | x | 2 + θ d S ,
where θ > 0 is a predetermined small constant and is selected as 1 × 10 6 here.

3.3. Constraints Information from UT Measurements

As mentioned above, the purpose of EIT and UT is to reconstruct the media distribution based on the change of conductivity or acoustic impedance from the boundary measurements. Since there is no change in the conductivity and acoustic impedance in the background medium region, the image reconstruction of the background medium satisfies the following relationship in EIT and UT:
Ω b = Ω E I T b ( δ σ = 0 ) = Ω U T b ( δ Z = 0 , δ μ = 0 ) ,
where Ω b represents the background medium region, Ω E I T b and Ω U T b are the background medium region reconstructed from EIT and UT, respectively.
Thus, the background medium region reconstructed from UT can be used as the prior information in EIT image reconstruction. Consider two scenarios whether there is a target on the ultrasonic wave’s propagation path, as shown in Figure 4a,b. If there is no target on the ultrasonic wave propagation path, the ultrasonic waves sending from the exciting transducer will continue to propagate to the measuring transducer and the attenuation of the acoustic energy is only caused by the background medium. And then, we can mark the pixels on the transmission ultrasonic waves propagation path as background medium region Ω U T b . In contrast, if there is a target on the ultrasonic wave propagation path, the ultrasonic waves sending from the exciting transducer will reflect on the target interface, and then continue to propagate to the measuring transducer through the background medium. The time of flight t f can be extracted from the reflection waves. And then, using the Equation (13) and the main lobe area of the transducer, we can calculate the region covered by the reflection ultrasonic waves propagation paths. After that, we can also mark the pixels on the reflection ultrasonic waves propagation path as background medium region Ω U T b .
Since the conductivity of the background medium region does not change, we can derive the following equation:
h ( x ) = C   x = 0 ,
where h ( x ) represents the constraint function of UT measurements with respect to the conductivity difference x, C is the corresponding extend constraint matrix, and it’s calculated by the analytical method:
C ( i , j ) = { 1 , i = 1 , , M U ,   j = n ( i ) 0 , otherwise ,
where C R M U × N , M U is the total number of pixels in region of Ω U T b , n ( i ) is the serial number of the pixels located at Ω U T b in the set of all pixels.

3.4. Lagrange-Newton Solution for Dual-Modality Inverse Problem

In summary, we can convert the media distribution image reconstruction problem into the conductivity change image reconstruction problem for the EIT/UT dual-modality imaging. The EIT/UT dual-modality inverse problem can be described as follows:
min x R N f ( x ) s .   t . h ( x ) = 0
where f(x) and h(x) are defined by Equations (15) and (19), respectively.
Obviously, the problem (21) is a typical equality constraint optimization problem. The Lagrange function combining with Newton method is commonly used to solve it. First, we formulate the following Lagrange function:
Ψ ( x , λ ) = f ( x ) λ T h ( x ) = 1 2 J x y 2 2 + β R ( x ) λ T h ( x ) ,
where λ is the Lagrange multiplier. The stability point of the above Lagrange function corresponds to the optimal solution of the unknown x, and satisfies the following conditions:
{ Ψ x = f ( x ) h ( x ) T λ = 0 Ψ λ = h ( x ) = 0 ,
Then, we solve the problem (22) using the Newton-Raphson method, and the iteration solution can be described as x k + 1 = x k + ( δ x ) k , where the iteration step length ( δ x ) k and the Lagrange multiplier ( δ λ ) k satisfy the following equation:
( W ( x k , λ k ) h ( x k ) T h ( x k ) 0 )   ( ( δ x ) k ( δ λ ) k ) = ( f ( x k ) h ( x k ) T λ k h ( x k ) ) ,
where
{ h ( x ) = ( C x ) = C f ( x ) = ( 1 2 J x y 2 2 + β R ( x ) ) = J T ( J x y ) + β L θ ( x ) x W ( x , λ ) = 2 ( Ψ ( x , λ ) ) = 2 ( f ( x ) λ T h ( x ) ) = J T J + β L θ ( x ) L θ ( x ) = L a T ( x ) H 1 ( x ) L a ( x ) , ,
where J T is the transpose of J, L a R Q × N is a sparse matrix, which can be determined by the topological relationship between the pixels, referring to the literature [45], H R Q × Q is a diagonal matrix whose diagonal entry can be described as follows:
H q q ( x ) = L a , q ( x ) x 2 + θ , q = 1 , 2 , , Q ,
where L a , q R 1 × N is the q-th row element of L a , Q is the total number of adjacent element pairs along the horizontal or vertical direction, referring to the literature [45].
To analyze the iteration convergence and determine the iteration termination condition of the problem (21), we define the following penalty function:
G ( x , λ ) = f ( x ) h ( x ) T λ 2 2 + h ( x ) 2 2 ,
In summary, the implement of the proposed EIT/UT dual-modality image reconstruction method can be described as Algorithm 1.
Algorithm 1: Lagrange-Newton Method for EIT/UT Dual-Modality Image Reconstruction
Step 1: Initialization.
  Calculate the initial value of x using linear back-projection (LBP) algorithm [13]: x 0 = J T y ,
  Given the value of the iteration termination parameter: ε 0 ,   ε = 1 × 10 6 ,
  Given the value of the intermediate parameter in iteration: κ ( 0 ,   1 ) ,   κ = 0.5 , η = 1 .
Step 2: Termination Condition Judgment.
  Calculate G ( x k ,   λ k ) by Equation (27), and judge:
  If G ( x k ,   λ k ) ε , stop iteration and return the x k as the optimal estimation value,
  Otherwise, calculate ( δ x ) k and ( δ λ ) k by the Equation (24).
Step 3: Compute Step Length.
  Calculate G ( x k + η ( δ x ) k , λ k + η ( δ λ ) k ) and G ( x k ,   λ k ) by Equation (27), and judge:
  If G ( x k + η ( δ x ) k , λ k + η ( δ λ ) k ) ( 1 κ η ) G ( x k , λ k ) , go to Step 4,
  Otherwise, set η = η / 4 , continue to the Step 3.
Step 4: Iteration Update.
  Update the conductivity difference x: x k + 1 = x k + η ( δ x ) k ,
  Update the Lagrange multiplier λ: λ k + 1 = λ k + η ( δ λ ) k ,
  Set k=k+1 and go to Step 2.

4. Results

To evaluate the performance of the proposed EIT/UT dual-modality imaging method, a series of numerical and experimental tests are carried out and compared with single-modality EIT method. In the tests, three typical EIT image reconstruction methods, Newton’s One-Step Error Reconstructor (NOSER), L1 regularization, and TV regularization, are used in single-modality EIT image reconstruction. Considering that the selection of regularization parameters has a great influence on the imaging results, the regularization parameters in all methods are empirically selected in the range of 10 6 to 10 3 , and the one that reached the most similar result to the real model was used in the comparison [45]. In EIT, the adjacent electric current excitation and adjacent voltage measurement model is used to obtain the boundary measurements. In UT, the piston transducers are used in numerical and experimental tests, and the acoustic energy can be better concentrated in a narrow beam range. In transmission mode, three transducers directly opposite to the exciting transducer are used to obtain the transmission measurements. In reflection mode, three transducers on each adjacent side of the exciting transducer are used to obtain the reflection measurements. The ultrasound frequency is 1MHz in the numerical and experimental tests. In addition, the reconstruction results from EIT/UT dual-modality method is simplified as UET in the paper.

4.1. Numerical Results and Discussion

The numerical tests are conducted by the commercial multi-physical coupling software COMSOL and MATLAB. The simulation model is designed by simulating the discrete bubbles in a conductive liquid. The conductivity of the bubble is 1 × 10 6 Sm−1 and the background medium conductivity is 1 × 10 2 Sm−1. The density of the bubble is 1.29 kgm−3 and sound velocity is 340 ms−1. The density of the background medium is 1000 kgm−3 and sound velocity is 1400 ms−1. The radius of the imaging domain is 75 mm, and 16 electrodes and 16 ultrasound transducers are staggered evenly and arranged on the periphery of the observation domain, as shown in Figure 3. The piston ultrasound transducers sequentially send and receive the three periods sinusoidal pulse signal in the narrow field range.

4.1.1. Quantitative Index

To evaluate the reconstruction results, four quantitative indexes of relative image error (RE), correlation coefficient (CC), position error (PE) and area difference (AD), are used in the paper. The RE and CC are usually used to describe the degree of difference of the true conductivity difference and reconstructed conductivity difference. The definition of RE and CC can be expressed as follows:
R E = x x 2 x 2 ,
C C = n = 1 N ( x n x ¯ ) ( x n x ¯ ) n = 1 N ( x n x ¯ ) 2 n = 1 N ( x n x ¯ ) 2 ,
where x is the calculated conductivity difference vector, x is the real conductivity difference vector, x n and x n are the n-th element of x and x , respectively, x ¯ and x ¯ are the mean values of x and x , respectively. A lower value of RE and a higher value of CC mean a good image reconstruction result.
The PE and AD are respectively used to evaluate the target position and shape reconstructed accuracy by comparing the segmented reconstructed images and the real distribution images. The segmented threshold is selected as 25% of the minimum pixel amplitudes, as follows:
[ x q ] n = { 1 , if [ x ] n 1 4 min ( x ) 0 , otherwise ,
where x is the calculated image pixel amplitudes (conductivity difference), x q is the segmented image pixel amplitudes. The PE and AD are defined as follows:
P E = | c t c q | ,
A D = sum ( x b x q ) sum ( x b x q ) N ,
where c t is the center of gravity of the target in real conductivity difference distribution image, c q is the center of gravity of the target in reconstructed conductivity difference distribution image, x b is the binary image of the real conductivity difference distribution, x q is the segmented image of the reconstructed conductivity difference distribution image, N is the total number of the pixels.

4.1.2. Image Reconstruction Results Analysis

In numerical tests, five different models are used to compare the proposed EIT/UT dual-modality imaging method (UET) and the single-modality EIT imaging method using NOSER, L1, TV. The image reconstruction results are shown in Figure 5, Figure 6, Figure 7 and Figure 8 when adding no noise, 40 dB, 25 dB, 20 dB Gaussian noise in EIT measurements, respectively. The real targets boundaries are indicated by the solid black line in the reconstructed images. It can be seen from the Figure 5, Figure 6, Figure 7 and Figure 8 that the UET has a cleaner background and better imaging performance compared with the other methods, which can be also verified by the quantitative indexes, as shown in Figure 9. For different quantitative indexes, we respectively calculate their mean values in five media distribution models under a certain noise condition. When the noise level increases, the imaging results from single-modality EIT using NOSER, L1 and TV are gradually getting worse, while the proposed UET method can always better reconstruct the media distribution in different noise levels. Thus, the proposed UET method is more robust to the EIT measurement noise. In addition, the expected time of UET is about 0.32 s for per image reconstruction. The expected time of EIT using NOSER, L1, TV are about 0.02 s, 0.2 s, 0.3 s for per image reconstruction, respectively. Thus, the imaging speed of UET is comparable to that of single-modality EIT using TV method, but the imaging accuracy of UET have an obvious improvement.

4.1.3. Comparison of Edge Reconstruction Ability

It can be seen from the numerical results that the UET has a cleaner edge than other methods. To verify this conclusion, we further analyze the numerical results of the model with one circle and one square in Figure 5. First, the pixel values of the reconstructed images are normalized to the interval [0, 1]. Then, the contour of the pixel values along the center horizontal line of the reconstructed images is drawn in Figure 10. It’s obvious that the contour of UET reconstructed image is closer to the real contour. Thus, the proposed UET method has a better edge reconstruction ability than the single-modality EIT methods.

4.2. Experimental Results and Discussion

To further verify the performance of the proposed method, a set of experimental tests was carried out. The experimental system is shown in Figure 11, which is developed by Tianjin University, China [46]. It mainly contains three parts, the sensor array, the data acquisition and control unit, and the image reconstruction and visualization unit. The radius of the sensor is 75 mm and 16 electrodes and 16 ultrasonic transducers are evenly placed on the periphery of the sensor. The contact resistance between the electrodes and background medium is 37.1 Ω/cm2. The background medium consists of Na2SO4 solution with the conductivity of 1.41 mS/cm. The data acquisition speed of EIT is about 1000-frames/s, and the data acquisition speed of UT is about 17-frames/s. The expected time of the image reconstruction in experimental tests is similar to numerical tests.
The experimental results are shown in Figure 12, where five different distribution models (A, B, C, D, E) of the nylon rod ( 10 11 mS/cm) are used. It can be seen from the results that the proposed UET method has a better imaging performance compared with the single-modality EIT methods. The quantitative analysis of the experimental results is shown in Figure 13, which can also verify the conclusion that the proposed UET method shows better reconstruction performance.

5. Conclusions

A novel EIT/UT dual-modality imaging method is proposed, where the measurement information obtained from UT can be used as prior constraint and improve the EIT image reconstruction performance. The dual-modality inverse problem is constructed by the equality constraint equation, and then, solved by the Lagrange-Newton method. The numerical and experimental tests show that the proposed EIT/UT dual-modality imaging method has a better image reconstruction performance than the traditional single-modality EIT methods, where the reconstructed images of UET have clearer targets edges and less artifacts in the background region. In the experimental tests, the average RE of UET have reduced 29.7%, 31.1%, 14.9% than single-modality EIT using NOSER, L1, TV, respectively. The average CC of UET have improved 7.7%, 16.9%, 7.4% than single-modality EIT using NOSER, L1, TV, respectively. The average PE of UET have reduced 32.9%, 68.2%, 71.3% than single-modality EIT using NOSER, L1, TV, respectively. The average AD of UET have reduced 55.1%, 50.6%, 36.5% than single-modality EIT using NOSER, L1, TV, respectively. In addition, the UET shows better performance of edge reconstruction, where the contour of UET reconstructed image is closer to the real contour than single-modality EIT methods. Furthermore, the proposed UET method is more robust to the EIT measurement noise, where the quantitative indexes of UET have no significant change with the increase of the measurement noise. The future work will focus on the research of enriching the UT measurements by adding the number of the transducers or increasing the diffusion angle of the transducers, and further improve the image reconstruction performance of EIT/UT dual-modality imaging.

Author Contributions

G.L. developed the algorithm and ran the testes; S.R. designed the study; S.Z. interpreted the data and images; F.D. designed the conception; all the authors were involved in writing and editing the manuscript.

Funding

This research was funded in part by the National Natural Science Foundation of China, grant number 61571321, and in part by the Natural Science Foundation of Tianjin City, grant number 17JCQNJC03500.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Dong, F.; Jiang, Z.X.; Qiao, X.T.; Xu, L.A. Application of electrical resistance tomography to two-phase pipe flow parameters measurement. Flow Meas. Instrum. 2003, 14, 183–192. [Google Scholar] [CrossRef]
  2. Yao, J.F.; Takei, M. Application of process tomography to multiphase flow measurement in industrial and biomedical fields—A review. IEEE Sens. J. 2017, 17, 8196–8205. [Google Scholar] [CrossRef]
  3. Liu, D.; Khambampati, A.K.; Du, J.F. A parametric level set method for electrical impedance tomography. IEEE Trans. Med. Imaging 2018, 37, 451–460. [Google Scholar] [CrossRef]
  4. Nissinen, A.; Kaipio, J.P.; Vauhkonen, M.; Kolehmainen, V. Contrast enhancement in EIT imaging of the brain. Physiol. Meas. 2015, 37, 1–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Storz, H.; Storz, W.; Jacobs, F. Electrical resistivity tomography to investigate geological structures of the earth’s upper crust. Geophys. Prospect. 2000, 48, 455–471. [Google Scholar] [CrossRef]
  6. Loke, M.H.; Chambers, J.E.; Rucker, D.F.; Kuras, O.; Wilkinson, P.B. Recent developments in the direct-current geoelectrical imaging method. J. Appl. Geophys. 2013, 95, 135–156. [Google Scholar] [CrossRef]
  7. Tallman, T.N.; Gungor, S.; Wang, K.W.; Bakis, C.E. Damage detection and conductivity evolution in carbon nanofiber epoxy via electrical impedance tomography. Smart Mater. Struct. 2014, 23, 045034. [Google Scholar] [CrossRef]
  8. Tallman, T.N.; Semperlotti, F.; Wang, K.W. Enhanced delamination detection in multifunctional composites through nanofiller tailoring. J. Intell. Mater. Syst. Struct. 2015, 26, 2565–2576. [Google Scholar] [CrossRef]
  9. Hallaji, M.; Seppänen, A.; Pour-Ghaz, M. Electrical resistance tomography to monitor unsaturated moisture flow in cementitious materials. Cem. Concr. Res. 2015, 69, 10–18. [Google Scholar] [CrossRef]
  10. Hallaji, M.; Seppänen, A.; Pour-Ghaz, M. Electrical impedance tomography-based sensing skin for quantitative imaging of damage in concrete. Smart Mater. Struct. 2014, 23, 085001. [Google Scholar] [CrossRef]
  11. Holder, D.S. Electrical Impedance Tomography: Methods, History and Applications; Institute of Physics Publishing: London, UK, 2005. [Google Scholar]
  12. Murphy, E.K.; Mahara, A.; Halter, R.J. Absolute reconstructions using rotational electrical impedance tomography for breast cancer imaging. IEEE Trans. Med. Imaging 2017, 36, 892–903. [Google Scholar] [CrossRef]
  13. Barber, D.C.; Brown, B.H. Applied potential tomography. J. Phys E Sci. Instrum. 1984, 17, 723. [Google Scholar] [CrossRef]
  14. 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]
  15. Liu, D.; Kolehmainen, M.; Siltanen, S.; Laukkanen, A.M.; Seppanen, A. Estimation of conductivity changes in a region of interest with electrical impedance tomography. Inverse Probl. Imag. 20l5, 9, 211–229. [Google Scholar] [CrossRef]
  16. Chen, B.; Abascal, J.; Soleimani, M. Electrical resistance tomography for visualization of moving objects using a spatiotemporal total variation regularization algorithm. Sensors 2018, 18, 1704. [Google Scholar] [CrossRef]
  17. Liu, S.H.; Jia, J.B.; Zhang, Y.D.; Yang, Y.J. Image reconstruction in electrical impedance tomography based on structure-aware sparse bayesian learning. IEEE Trans. Med. Imaging 2018, 37, 2090–2102. [Google Scholar] [CrossRef]
  18. Ren, S.J.; Sun, K.; Liu, D.; Dong, F. A statistical shape constrained reconstruction framework for electrical impedance tomography. IEEE Trans. Med. Imaging 2019, in press. [Google Scholar] [CrossRef]
  19. Cherry, S.R. Multimodality imaging: beyond PET/CT and SPECT/CT. Semin. Nucl. Med. 2009, 39, 348–353. [Google Scholar] [CrossRef]
  20. Veninga, T.; Huisman, H.; Maazen, R.W.; Huizenga, H. Clinical validation of the normalized mutual information method for registration of CT and MR images in radiotherapy of brain tumours. J. Appl. Clin. Med. Phys. 2004, 5, 66–79. [Google Scholar] [CrossRef]
  21. Kaplan, I.; Oldenburg, N.E.; Meskell, P.; Blake, M.; Church, P.; Holupka, E.J. Real time MRI-ultrasound image guided stereotactic prostate biopsy. Magn. Reson. Imaging 2002, 20, 295–299. [Google Scholar] [CrossRef]
  22. Steiner, G. Sequential fusion of ultrasound and electrical capacitance tomography. Int. J. Inf. Syst. Sci. 2006, 2, 484–497. [Google Scholar]
  23. Zhang, R.H.; Wang, Q.; Wang, H.X.; Zhang, M.; Li, H.L. Data fusion in dual-mode tomography for imaging oil–gas two-phase flow. Flow Meas. Instrum. 2014, 37, 1–11. [Google Scholar] [CrossRef]
  24. Johansen, G.A.; Abro, E. A new CdZnTe detector system for low-energy gamma-ray measurement. Sens. Actuators A Phys. 1996, 54, 493–498. [Google Scholar] [CrossRef]
  25. Liu, S.; Liu, S.; Ren, T. Acoustic tomography reconstruction method for the temperature distribution measurement. IEEE Trans. Instrum. Meas. 2017, 66, 1936–1945. [Google Scholar] [CrossRef]
  26. Duric, N.; Boyd, N.; Littrup, P.; Sak, M.; Myc, L.; Li, C.; West, E.; Minkin, S.; Martin, L.; Yaffe, M.; et al. Breast density measurements with ultrasound tomography: A comparison with film and digital mammography. Med. Phys. 2013, 40, 013501. [Google Scholar] [CrossRef] [Green Version]
  27. Langener, S.; Vogt, M.; Ermert, H.; Musch, T. A real-time ultrasound process tomography system using a reflection-mode reconstruction technique. Flow Meas. Instrum. 2017, 53, 107–115. [Google Scholar] [CrossRef]
  28. Liang, G.H.; Ren, S.J.; Dong, F. An inclusion boundary reconstruction method using electrical impedance and ultrasound reflection dual-modality tomography. In Proceedings of the 9th World Congress on Industrial Process Tomography (WCIPT), Bath, UK, 2–6 September 2018. [Google Scholar]
  29. Liang, G.H.; Ren, S.J.; Dong, F. Ultrasound guided electrical impedance tomography for 2D free-interface reconstruction. Meas. Sci. Technol. 2017, 28, 074003. [Google Scholar] [CrossRef]
  30. Steiner, G.; Soleimani, M.; Watzenig, D. A bio-electromechanical imaging technique with combined electrical impedance and ultrasound tomography. Physiol. Meas. 2008, 29, 63–75. [Google Scholar] [CrossRef] [PubMed]
  31. Wan, Y.; Halter, R.; Borsic, A.; Manwaring, P.; Hartov, A.; Paulsen, K. Sensitivity study of an ultrasound coupled transrectal electrical Impedance Tomography system for prostate imaging. Physiol. Meas. 2010, 31, 17–29. [Google Scholar] [CrossRef] [PubMed]
  32. Teniou, S.; Meribout, M. A multimodal image reconstruction method using ultrasonic waves and electrical resistance tomography. IEEE Trans. Image Process. 2015, 24, 3512–3521. [Google Scholar] [CrossRef] [PubMed]
  33. Pusppanathan, J.; Rahim, R.A.; Phang, F.A.; Mohamad, E.J.; Seong, C.K. Single-plane dual-modality tomography for multiphase flow imaging by integrating electrical capacitance and ultrasonic sensors. IEEE Sens. J. 2017, 17, 6368–6377. [Google Scholar] [CrossRef]
  34. Liang, G.H.; Ren, S.J.; Dong, F. An augmented Lagrangian trust region method for inclusion boundary reconstruction using ultrasound/electrical dual-modality tomography. Meas. Sci. Technol. 2018, 29, 074008. [Google Scholar] [CrossRef] [Green Version]
  35. Hassan, H.; Semperlotti, F.; Wang, K.W.; Tallman, T.N. Enhanced imaging of piezoresistive nanocomposites through the incorporation of nonlocal conductivity changes in electrical impedance tomography. J. Intell. Mater. Syst. Struct. 2018, 29, 1850–1861. [Google Scholar] [CrossRef]
  36. Zhao, L.; Yang, J.; Wang, K.W.; Semperlotti, F. An application of impediography to the high sensitivity and high resolution identification of structural damage. Smart Mater. Struct. 2015, 24, 065044. [Google Scholar] [CrossRef] [Green Version]
  37. Ammari, H.; Bonnetier, E.; Capdeboscq, Y.; Tanter, M.; Fink, M. Electrical impedance tomography by elastic deformation. SIAM J. Appl. Math. 2008, 68, 1557–1573. [Google Scholar] [CrossRef]
  38. Somersalo, E.; Cheney, M.; Isaacson, D. Existence and uniqueness for electrode models for electric current computed tomography. SIAM J. Appl. Math. 1992, 52, 1023–1040. [Google Scholar] [CrossRef]
  39. Geselowitz, D. An application of electrocardiographic lead theory to impedance plethysmography. IEEE Trans. Biomed. Eng. 1971, 18, 38–41. [Google Scholar] [CrossRef] [PubMed]
  40. Ren, S.J.; Soleimani, M.; Dong, F. Inclusion boundary reconstruction and sensitivity analysis in electrical impedance tomography. Inverse Probl. Sci. Eng. 2018, 26, 1037–1061. [Google Scholar] [CrossRef]
  41. Ren, S.J.; Wang, Y.; Liang, G.H.; Dong, F. A robust inclusion boundary reconstructor for electrical impedance tomography with geometric constraints. IEEE Trans. Instrum. Meas. 2019, 68, 762–773. [Google Scholar] [CrossRef]
  42. Li, N.; Cao, M.C.; Xu, K.; Jia, J.B.; Du, H.B. Ultrasonic transmission tomography sensor design for bubble identification in gas-liquid bubble column reactors. Sensors 2018, 18, 4256. [Google Scholar] [CrossRef] [PubMed]
  43. Brandstatter, B. Jacobian calculation for electrical impedance tomography based on the reciprocity principle. IEEE Trans. Magn. 2003, 39, 1309–1312. [Google Scholar] [CrossRef]
  44. Rahiman, M.H.F.; Rahim, R.A.; Rahim, H.A.; Zakaria, Z.; Pusppanathan, J. A study on forward and inverse problems for an ultrasonic tomography. J. Teknol. 2014, 70, 113–117. [Google Scholar] [CrossRef]
  45. Song, X.Z.; Xu, Y.B.; Dong, F. A spatially adaptive total variation regularization method for electrical resistance tomography. Meas. Sci. Technol. 2015, 26, 125401. [Google Scholar] [CrossRef]
  46. Dong, F.; Xu, C.; Zhang, Z.Q.; Ren, S.J. Design of parallel electrical resistance tomography system for measuring multiphase flow. Chin. J. Chem. Eng. 2012, 20, 368–379. [Google Scholar] [CrossRef]
Figure 1. Ultrasound transmission and reflection physical model.
Figure 1. Ultrasound transmission and reflection physical model.
Sensors 19 01966 g001
Figure 2. 3-D sensitivity distribution: (a) Electrical impedance tomography (EIT), (b) ultrasound tomography (UT).
Figure 2. 3-D sensitivity distribution: (a) Electrical impedance tomography (EIT), (b) ultrasound tomography (UT).
Sensors 19 01966 g002
Figure 3. Mesh model of EIT/UT dual-modality imaging.
Figure 3. Mesh model of EIT/UT dual-modality imaging.
Sensors 19 01966 g003
Figure 4. Schematic diagram of ultrasonic wave propagation path: (a) there is no target on the ultrasonic waves propagation path, (b) there is a target on the ultrasonic waves propagation path.
Figure 4. Schematic diagram of ultrasonic wave propagation path: (a) there is no target on the ultrasonic waves propagation path, (b) there is a target on the ultrasonic waves propagation path.
Sensors 19 01966 g004
Figure 5. Reconstruction results with no noise.
Figure 5. Reconstruction results with no noise.
Sensors 19 01966 g005
Figure 6. Reconstruction results with 40 dB noise.
Figure 6. Reconstruction results with 40 dB noise.
Sensors 19 01966 g006
Figure 7. Reconstruction results with 25 dB noise.
Figure 7. Reconstruction results with 25 dB noise.
Sensors 19 01966 g007
Figure 8. Reconstruction results with 20 dB noise.
Figure 8. Reconstruction results with 20 dB noise.
Sensors 19 01966 g008
Figure 9. Mean value of the quantitative indexes.
Figure 9. Mean value of the quantitative indexes.
Sensors 19 01966 g009
Figure 10. Contour of the reconstructed images along the center horizontal line.
Figure 10. Contour of the reconstructed images along the center horizontal line.
Sensors 19 01966 g010
Figure 11. EIT/UT dual-modality data acquisition system.
Figure 11. EIT/UT dual-modality data acquisition system.
Sensors 19 01966 g011
Figure 12. Reconstruction results of the experimental tests.
Figure 12. Reconstruction results of the experimental tests.
Sensors 19 01966 g012
Figure 13. Quantitative analysis of the experimental results.
Figure 13. Quantitative analysis of the experimental results.
Sensors 19 01966 g013

Share and Cite

MDPI and ACS Style

Liang, G.; Ren, S.; Zhao, S.; Dong, F. A Lagrange-Newton Method for EIT/UT Dual-Modality Image Reconstruction. Sensors 2019, 19, 1966. https://doi.org/10.3390/s19091966

AMA Style

Liang G, Ren S, Zhao S, Dong F. A Lagrange-Newton Method for EIT/UT Dual-Modality Image Reconstruction. Sensors. 2019; 19(9):1966. https://doi.org/10.3390/s19091966

Chicago/Turabian Style

Liang, Guanghui, Shangjie Ren, Shu Zhao, and Feng Dong. 2019. "A Lagrange-Newton Method for EIT/UT Dual-Modality Image Reconstruction" Sensors 19, no. 9: 1966. https://doi.org/10.3390/s19091966

APA Style

Liang, G., Ren, S., Zhao, S., & Dong, F. (2019). A Lagrange-Newton Method for EIT/UT Dual-Modality Image Reconstruction. Sensors, 19(9), 1966. https://doi.org/10.3390/s19091966

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