Next Article in Journal
Comparative Study for Sentiment Analysis of Financial Tweets with Deep Learning Methods
Previous Article in Journal
Integration of Active Clothing with a Personal Cooling System within the NGIoT Architecture for the Improved Comfort of Construction Workers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Computational Imaging Algorithm for Electrical Capacitance Tomography

1
Aircraft Optical Imaging Monitoring and Measurement Technology Laboratory, Xi’an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences, Xi’an 710119, China
2
School of Control and Computer Engineering, North China Electric Power University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(2), 587; https://doi.org/10.3390/app14020587
Submission received: 21 December 2023 / Revised: 4 January 2024 / Accepted: 5 January 2024 / Published: 10 January 2024

Abstract

:
High-precision images enable electrical capacitance tomography (ECT) to obtain more reliable measurement results, meaning that the reconstruction algorithm is particularly important. Some excellent numerical algorithms have successfully solved the inverse problem for ECT imaging, but their imaging quality is relatively low. To solve this problem, this paper proposes a new reconstruction algorithm based on regularized extreme learning machines (RELMs). The implementation of the algorithm is mainly divided into two steps: (1) according to a large number of training samples, the RELM model can be obtained by the iterative split Bregman (ISB) algorithm, which can describe the mapping relationship between the capacitance correlation coefficient and the imaging target well, and (2) the capacitance correlation coefficient is calculated, which is then used as input to the RELM model to predict the final imaging. Both simulation and experimental results show that the RELM algorithm achieves greater improvement in imaging quality and robustness, and provides new development ideas for the ECT.

1. Introduction

Electrical capacitance tomography is a promising process imaging technique. It enables visualization of multiphase flow by reconstructing the distribution of multiphase media in the pipe by measuring the capacitance vector between arrays of electrodes outside the pipe. ECT is a growing research topic, with advantages such as its wide range of applications, non-invasiveness, fast response time, high security performance, low cost, etc. It is widely used in the multiphase flow process detection of insulating media [1,2,3,4], such as monitoring flames and combustion conditions in engine combustion chambers, defect detection of thermal protection materials, measuring solid concentration in cyclone separator, pneumatic measurement of high pressure pulverized coal, etc.
Image reconstruction is an important process for visualizing the detection results of ECT systems. The imaging quality of the ECT determines whether the ECT technique can be practically implemented in engineering applications. However, ECT image reconstruction process is severely ill-conditioned, so it is necessary to study image reconstruction algorithms to weaken this property. At present, some excellent numerical algorithms have successfully solved the inverse problem for ECT imaging. Two major categories of ECT reconstruction methods exist—non-iterative and iterative image reconstruction algorithms [5]. The advantage of the non-iterative imaging algorithm is that as a “one-step” imaging algorithm, it has faster imaging speed and can complete simple imaging tasks. When faced with complex imaging targets, non-iterative imaging algorithms cannot reconstruct the measured objects well and the image quality is poor. Therefore, non-iterative imaging algorithms are mostly used as online imaging algorithms for qualitative analysis or to provide initial values for iterative imaging algorithms. Common iterative algorithms include the gene algorithm (GA) [6], the Landweber algorithm [7], the Newton–Raphson iterative algorithm (NRIA) [8], the Kalman filtering algorithm (KFA) [9], and the extreme learning machine (ELM) algorithm [10]. Compared to non-iterative imaging algorithms, iterative imaging algorithms improve the image reconstruction quality, but this improvement is limited by the high time cost due to many iterations and the difficulty in choosing the parameters of the algorithm. In the past, the aforementioned algorithms have been validated in a number of simple practical applications, which has motivated the development of the ECT techniques. However, as the demand for image reconstruction quality continues to increase, the above algorithms cannot meet the needs of the practical applications of the ECT. To achieve higher image reconstruction quality, this paper proposes a more efficient image reconstruction algorithm.
The stability and robustness of the conventional ELM algorithm itself are not very good, and the image quality is low when used directly for image reconstruction. To address this issue, combining the advantages of the ELM algorithm and the Tikhonov regularization (TR) method, a new image reconstruction algorithm aiming to enhance reconstruction accuracy is proposed. In this paper, the RELM method provides several appealing attributes enumerated in the following:
(1)
The implementation of the algorithm mainly consists of two stages: at the learning stage, according to the training sample set, the RELM model between capacitance correlation coefficient and imaging target is extracted; at the prediction stage, by taking the calculated capacitance correlation coefficient as the input of the RELM model, the corresponding reconstructed image can be predicted.
(2)
The solution method of the capacitance correlation coefficient is analyzed in detail. Based on the TR method, the objective function of the RELM algorithm is designed, in which the L1-norm as data fidelity effectively reduces the influence of outliers in the input data on the reconstruction results, the second-order total variation (STV) norm and the L1-norm as the regularizer effectively improves the stability of the numerical solution. The objective function is solved by the ISB method, which greatly improves the solution efficiency.
(3)
Numerical simulation and experimental results prove that the RELM algorithm improves the stability of the numerical solution, and the problem of poor imaging quality caused by inaccurate capacitance data is effectively alleviated.
The rest of this paper is organized as follows: in Section 2, we propose an imaging methodology based on the RELM. In Section 3, we detail the method for solving the capacitance correlation coefficient. In Section 4, the RELM model between the capacitance correlation coefficient and the imaging target is extracted from a large number of training samples. The numerical simulations and implementation details are described in Section 5 and Section 6. Finally, the conclusions are drawn in Section 7.

2. Imaging Methodology Based on Regularized Extreme Learning Machines

Imaging quality is crucial for ECT applications, so we must find image reconstruction algorithms with strong imaging capabilities. The idea of the RELM algorithm mainly comes from the actual ECT experimental procedure. Through many experiments, it is found that there is a certain relationship between the capacitance correlation coefficient and the imaging target. Therefore, it is particularly important to find the mapping relationship between the two and use this mapping relation to solve the problem of low reconstruction accuracy for the ECT. In this paper, based on the existing numerical methods, we propose a RELM-based computational image reconstruction algorithm that fully exploits the mapping relationship between capacitance correlation coefficients and imaging targets, and achieves high reconstruction accuracy. The specific implementation process can be carried out as follows:
(1)
The training sample set between the capacitance correlation coefficient and the imaging target is obtained through numerical simulation or experiment.
(2)
The mapping model between the capacitance correlation coefficient and the imaging target is obtained by the RELM method.
(3)
The capacitance correlation coefficient corresponding to the measured capacitance value is calculated. The calculation method is the same as that of the training sample set.
(4)
The calculated capacitance correlation coefficient is used as the input of the mapping model obtained in the second step, and the final reconstruction image is predicted.
The proposed RELM-based image reconstruction algorithm is different from the regular ELM algorithm, and their imaging process block diagrams are presented in Figure 1 and Figure 2, respectively. The RELM algorithm has the following advantages: (1) the algorithm effectively uses the prior information obtained from numerical simulations or experimental results, which greatly improves the quality of imaging results; (2) when the mapping model is obtained from the training sample set, the algorithm does not directly establish the imaging model between the measured capacitance and the imaging target, but rather establishes the mapping relationship model between the capacitance correlation coefficient and the true images of targets; and (3) based on the TR technology, the conventional ELM method is improved, which effectively enhances the stability of the numerical solution, and the solution efficiency is greatly improved by the ISB algorithm.
The RELM model has some limitations, such as resolution and number of sensors. When the number of sensors is different and the measured area is divided into different numbers of microelements, the capacitance correlation coefficients obtained are different, which means that the input of the model is different. Therefore, it is necessary to retrain the mapping model.

3. Solving for Capacitance Correlation Coefficients

For two random variables X = x 1 , x 2 , x m and Y = y 1 , y 2 , y m , the correlation coefficient between them is as follows [11]:
γ X , Y = i x i x i ¯ y i y i ¯ i x i x i ¯ 2 i y i y i ¯ 2
where x i ¯ and y i ¯ are the mean values of X and Y , respectively. γ X , Y reflects the correlation degree of the two vectors ( X and Y ), and its value is between −1 and 1.
According to the research conclusions of relevant scholars for decades, it has been a consensus in the industry that the image reconstruction model of the ECT can be expressed using the following mathematical formula when considering capacitive noise [12]:
C + α = S G
where C is the capacitance vector, S is sensitivity matrix, G is the permittivity vector, and α is noise capacitance vector. The image reconstruction process of the ECT is essentially a process of solving Equation (2). However, Equation (2) is only a basic model, and it is very difficult to solve it directly. Therefore, we must find effective numerical methods to solve it.
When the number of electrodes is n and the measured area is divided into t units, the high permittivity medium ( ε 2 ) is placed in the first micro unit and the low permittivity medium ( ε 1 ) is placed in the other units. All capacitances between electrodes are measured and normalized to obtain capacitance vector C s 1 . Similarly, the normalized capacitance matrix can be obtained as follows:
C s = ( C s 1 , C s 2 , , C s i , , C s t )
where C s i = ( c 1 , c 2 , , c m ) T , C s i R m × 1 , C s R N × t , m = n ( n 1 ) 2 . C s is essentially a sensitivity field matrix.
When there is a medium to be measured in the measured area, the measured capacitance value is C d ( C d R m × 1 ). According to Equation (1), the capacitance correlation coefficient γ i = γ C s i , C d can be obtained, and the following capacitance correlation coefficient vector can be obtained in the same way.
γ = γ 1 , γ 2 , , γ t
where γ R 1 × t is essentially a dielectric constant matrix and can be directly used for imaging. Reference [13] has confirmed the advantages of this method. In this paper, we take this further and use it as the input of the RELM model to enhance the imaging quality.

4. Regularized Extreme Learning Machine

In the process of numerical simulation and experiment, we found that there is a certain correspondence between the capacitance correlation coefficient and the imaging target. How to extract and utilize this relationship to improve the stability and reliability of reconstruction results has become the biggest challenge. In view of the excellent numerical performance of the ELM algorithm [14], it is used to extract nonlinear mapping relationships between the capacitance correlation coefficient and the imaging target.

4.1. Extreme Learning Machine

Machine learning algorithms represent nonlinear models of data relationships in a data-driven manner, enabling the effective prediction of new data. As a representative method of machine learning, the ELM algorithms are suitable for image reconstruction due to their advantages of theoretical analytical solutions and fast convergence [15]. The extreme learning machine belongs to the single hidden layer feedforward neural network (SLFN), and it can automatically adjust the structure based on the prediction results of current data segments without manually adjusting parameters and retraining the learning model [16].
For a group of training samples, x i , y i , i = 1,2 , , N , x i ( x i = x i 1 , x i 2 , x i n T R n ) is the i-th sample, y i ( y i = y i 1 , y i 2 , y i m T R m ) is its corresponding label, where n and m are the number of features and categories, respectively [16]. The basic structure of the ELM algorithm is shown in Figure 3. The output of the hidden layer and output layer can be expressed by the following equations [15]:
h = g ( a , b , x )
h x i V = y i , i = 1,2 , , N
where g ( x ) is the activation function.
Equation (6) can be written in matrix form as follows:
H V = Y
where
H = g ( a 1 , b 1 , x 1 ) g ( a 1 , b 1 , x 2 ) g ( a 1 , b 1 , x N ) g ( a 2 , b 2 , x 1 ) g ( a 2 , b 2 , x 2 ) g ( a 2 , b 2 , x N ) g ( a L , b L , x 1 ) g ( a L , b L , x 2 ) g ( a L , b L , x N ) T
V = v 1 T v 2 T v L T , Y = y 1 T y 2 T y N T
where H is the output matrix of the hidden layer, V is the output weight matrix connecting the hidden layer and the output layer, Y is the output matrix of the output layer, a is the input weight, and b is the offset of the hidden layer [17].
In the ELM algorithm, the input weight, a, and the offset, b, are random values without manual setting. When they are determined, the unique output weight matrix, V, can be obtained according to Equation (7). The calculation flow of output weight matrix is shown in Figure 4.

4.2. Regularized Extreme Learning Machine

As a powerful machine learning algorithm, the ELM algorithm has been successfully applied to the ECT imaging, but the stability and robustness of the ELM algorithm itself are not very good. When the conventional ELM algorithm is directly used for image reconstruction, it leads to the following two problems: (1) when directly calculating the least square solution of the cost function, the user cannot change the parameters of the ELM algorithm according to the characteristics of the training sample set, which makes the controllability of the conventional ELM algorithm worse; and (2) when the data set contains some disturbances the conventional ELM algorithm will be seriously affected, making the imaging quality worse. To solve these intractable problems, this paper proposes the RELM algorithm based on the TR technique. The TR theory usually transforms ill-posed problems into optimization problems to solve. Its objective function consists of two parts: the data fidelity term and the regularizer [18]. The generalized expression is as follows:
min V P V + i = 1 n λ i Q i V
where P V is the data fidelity term, which measures the difference between H V and Y , Q i V is the regularizer, and λ i is the regularization parameter, which adjusts the proportion of the data fidelity term and the regularizer, and λ i > 0 [19].
When designing the objective function based on the TR technology, introducing the prior information and the characteristics of the imaging targets into the inverse problem-solving process is helpful to improve the accuracy and stability of the numerical solution. Theoretically, a variety of estimation functions can be used as the data fidelity term and the regularizer, and different choices determine the different stability and accuracy of the numerical solution, so their selection is particularly important. Because reconstructed images in ECT generally have sparse features, the L1 norm has the advantage of generating sparse solutions compared to other norms. Therefore, we deploy the L1-norm as the data fidelity term to enhance the accuracy of numerical solution [20], whose mathematical expression is as follows:
P V = H V Y 1
In addition to the data fidelity term of the solution, the design of the regularizer is also very crucial. The first-order total variation (FTV) regularization has the ability to preserve the discontinuity characteristics of the measured region [21], but this method commonly leads to the ladder effect, so we deploy the STV regularization as the regularizer to get a more stable numerical solution. The STV norm has incomparable advantages in enhancing the sparse characteristics and improving imaging ability for the edge of the imaging target compared to the FTV norm [22]. However, when the STV regularization is employed alone, the imaging results may still have the ladder effect. Therefore, a new penalty function must be added to eliminate the imaging artifacts caused by the STV regularization. Because most imaging targets have sparse characteristics, we add additional sparse constraints through the L1 regularization. The mathematical expression of the regularization term is as follows:
Q V = λ 1 V 1 + λ 2 D x x V + D x y V + D y x V + D y y V
where D x x , D x y , D y x , and D y y are different difference operators, and λ 1 and λ 2 are regularization parameters.
As a result of submitting Equations (11) and (12) into Expression (10), a new objective function can be obtained:
min V H V Y 1 + λ 1 V 1 + λ 2 D x x V + D x y V + D y x V + D y y V
Suppose that W V 1 = D x x V + D x y V + D y x V + D y y V , then Expression (13) can be rewritten as:
min V H V Y 1 + λ 1 V 1 + λ 2 W V 1
where W = [ D x x , D x y , D y x , D y y ] T .
To solve Expression (14), the ISB algorithm is preferred. Expression (14) can be converted to the following form [23]:
min V d 1 1 + λ 1 V 1 + λ 2 d 2 1 s .   t .   d 1 = H V Y , d 2 = W V
According to the ISB algorithm theory, Expression (15) can be rewritten as follows:
V k + 1 , d 1 , k + 1 , d 2 , k + 1 = min V , d 1 , d 2 d 1 1 + λ 1 V 1 + λ 3 2 d 1 H V + Y b 1 , k 2 2 + λ 2 d 2 1 + λ 4 2 d 2 W V b 2 , k 2 2
b 1 , k + 1 = b 1 , k + H V k + 1 Y d 1 , k + 1
b 2 , k + 1 = b 2 , k + W V k + 1 d 2 , k + 1
where λ 1 and λ 2 are regularization parameters.
Equation (16) can be decoupled into three simple sub equations for solution, as shown in Equations (19)–(21) [24]:
V k + 1 = min V λ 1 V 1 + λ 3 2 H V Y + b 1 , k d 1 , k 2 2 + λ 4 2 W V + b 2 , k d 2 , k 2 2
d 1 , k + 1 = min d 1 d 1 1 + λ 3 2 d 1 H V k + 1 + Y b 1 , k 2 2
d 2 , k + 1 = min d 2 λ 2 d 2 1 + λ 4 2 d 2 W V k + 1 b 2 , k 2 2
Solving Equation (16) is equivalent to solving Equations (19)–(21). The solutions of the above three equations can be obtained by iteratively shrinking the threshold (IST) algorithm [25], and the iterative formats of their solutions are as follows:
V k + 1 = s o f t ( V k λ 3 H T H V k Y + b 1 , k d 1 , k λ 4 W T W V k + b 2 , k d 2 , k , λ 1 )
d 1 , k + 1 = s o f t ( H V k + 1 Y + b 1 , k , 1 λ 3 )
d 2 , k + 1 = s o f t ( W V k + 1 + b 2 , k , λ 2 λ 3 )
According to the above discussion, we have obtained a new image reconstruction algorithm, whose execution steps can be summarized as Algorithm 1.
Algorithm 1 The RELM algorithm
Input: parameter initialization
If the iterative condition is met, then
Update V k + 1 : V k + 1 = s o f t ( V k λ 3 H T H V k Y + b 1 , k d 1 , k λ 4 W T W V k + b 2 , k d 2 , k , λ 1 )
Update d 1 , k + 1 : d 1 , k + 1 = s o f t ( H V k + 1 Y + b 1 , k , 1 λ 3 )
Update d 2 , k + 1 : d 2 , k + 1 = s o f t ( W V k + 1 + b 2 , k , λ 2 λ 3 )
Update b 1 , k + 1 : b 1 , k + 1 = b 1 , k + H V k + 1 Y d 1 , k + 1
Update b 2 , k + 1 : b 2 , k + 1 = b 2 , k + W V k + 1 d 2 , k + 1
k++
Else
Output: Y ( H V k + 1 = Y )
end
The RELM algorithm is an improved and optimized version of the ELM technique, integrating the advantages of the ELM technique and TR technology, effectively overcoming the ill-conditioned nature of the ECT inverse problem and the adverse effects of capacitive inaccuracy. The RELM algorithm provides a significant improvement in reliability and reconstruction precision compared to the ELM technique. Its process is shown in Figure 5.

5. Numerical Simulation and Analysis

To validate the reconstruction effect of the regularized extreme learning machine (RELM) algorithm, we carried out a series of numerical simulation experiments. As a result of comparing the reconstruction results of the linear back projection (LBP) technique, the Tikhonov regularization (TR) technique, the truncated singular value decomposition (TSVD) technique, the conjugate gradient (CG) technique, the Landweber technique, the simultaneous iterative reconstruction technique (SIRT), the Gaussian process regression (GPR) method, the extreme learning machine (ELM) method, and the RELM technique, this section qualitatively and quantitatively analyzes the excellent numerical characteristics of the RELM algorithm. The square sensor with 12 electrodes, as shown in Figure 6, is adopted in this simulation study. The parameters of the sensor are listed in Table 1. Through numerical simulation and experimental measurements, we obtained 1500 sets of training samples in this study. The test sample set data is not included in the training sample set. In the sample generation process for the simulation experiments as well as the actual experiments, the inputs are the capacitance correlation coefficients and the outputs are the corresponding true images. In Table 2, we show several training sample pairs. We use relative error (RE) to evaluate the quality of reconstruction results, which is defined as follows:
R E = Y ^ i Y i 2 Y i 2
where Y is the gray value of the real image, and Y ^ is the gray value of the predicted image. The smaller the RE, the higher the reconstructed image quality.

5.1. Case 1

To test the performance of the RELM algorithm, seven groups of simulation experiments are carried out in this section, and the imaging quality of various algorithms is compared. The seven kinds of dielectric constant distributions shown in Figure 7 are taken as the phantom objects for numerical simulation, in which the high permittivity is the yellow area and the low permittivity is the blue area. In Figure 7a, the radius of the circle is 10 mm. For Figure 7b, the radius of the circle is 8 mm. The height of laminar flow in Figure 7c is 20 mm. In Figure 7d and Figure 7g, the sides of the square are 16 mm and 15 mm, respectively. For Figure 7e, there are two rectangles with a length of 60 mm and a width of 15 mm. In Figure 7f, the side length of the square is 10 mm, and the radius of the circle is 5 mm. The TR algorithm with the parameter value of 0.01 provides the initial value for the iterative algorithm mentioned in this paper. The parameters of the CG algorithm, the Landweber algorithm, and the SIRT are listed in Table 3, Table 4 and Table 5, respectively. The number of hidden layer nodes in both ELM and RELM is 1000. The reconstruction results of various algorithms are shown in Figure 8 and Figure 9. To quantitatively evaluate the reconstruction accuracy of the various algorithms, we calculate the REs, shown in Table 6, between the reconstructed and true permittivity.
The imaging results of the LBP technique, the TR technique, and the TSVD method are shown in Figure 8. Although they can roughly recognize the contour of the phantom objects, their imaging results are poor when faced with complex imaging targets. In Figure 9, we can see that the reconstruction results of iterative algorithm, i.e., the CG technique, the Landweber technique, and the SIRT, are significantly superior over non-iterative algorithms. However, the imaging results of all three algorithms still have distortion, which cannot reconstruct the phantom objects well.
In Figure 9(a4–g4) and Figure 9(a5–g5), the GPR algorithm and the ELM algorithm take full use of the advantages of the machine learning algorithm, which considerably improves the final imaging quality compared to the aforementioned traditional imaging algorithms in this paper. When the imaging target is complex or has a right angle, the imaging quality of the classical imaging algorithms is poor. Due to the learning of a large number of datasets, the GPR algorithm and the ELM algorithm can perform better on imaging tasks. However, from Figure 9(a4–g4) and Figure 9(a5–g5), we see that for seven kinds of imaging targets, the imaging results of the ELM algorithm and the GPR algorithm still have different degrees of artifacts. From Figure 9(a6–g6), we observe that the RELM algorithm has less artifacts compared to the ELM algorithm, and the imaging results are closest to the phantom objects. The reason why the RELM can improve the image reconstruction accuracy in engineering applications is that it not only has strong generalization ability and classification capabilities, but also improves the robustness of the numerical solution. From the above analysis, the RELM algorithm can be divided into two steps: (1) the RELM mapping model between the capacitance correlation coefficient and the imaging target is extracted based on the training samples, and (2) the capacitance correlation coefficient is calculated and then it is used as the input of the RELM model to predict the final imaging. The implementation process is relatively simple, and there is no need to manually set the weights and offsets of the hidden layers, which is in line with the practical needs of the project.
High-quality reconstructed images are an important guarantee for reliable measurement results. From the results listed in Table 6, we see that, compared with other algorithms mentioned, the RELM algorithm achieves the highest imaging accuracy for different simulation targets. For phantom objects in Figure 7, the REs of the RELM algorithm are 2.13%, 2.06%, 1.09%, 1.68%, 1.49%, 1.84%, and 1.31%, respectively. The quantitative error comparison results successfully verify the RELM algorithm has high reliability in solving the inverse problem of ECT and greatly improves the reconstruction accuracy.

5.2. Case 2

The robustness of image reconstruction algorithm is another important index to evaluate the algorithm. In practical application, the ECT system will be affected by noise when collecting capacitance data. Small input errors will cause the output result to deviate significantly from the true result. Consequently, to make the simulation results of the RELM algorithm more convincing, we add 5% and 10% random noise to the original capacitance to assess the RELM algorithm. The phantom objects to be reconstructed and the imaging results are presented in Figure 7 and Figure 10 and Figure 11. We executed 100 computations, and the image error mean is displayed in Table 7. The noise level is defined by the following equation:
η = C ^ C C × 100 %
where η is the capacitance noise level, C ^ is the noisy capacitance data, and C is the raw capacitance data.
From the reconstruction results in Figure 10 and Figure 11, we find that under the noise levels of 5% and 10%, the imaging results of phantom objects are still satisfactory, which shows that the RELM algorithm has strong robustness and stability under different noise levels. As can be seen quantitatively from the image error mean values in Table 7, the image errors still keep a small value under 10% noise level. The image error means of phantom object (a)–phantom object (g) are 4.51%, 4.16%, 4.79%, 3.61%, 4.27%, 3.79%, and 3.43%, respectively. The conclusion of the qualitative analysis and quantitative analysis shows that the RELM algorithm has a strong anti-noise performance under different levels of noise for various complex conditions.

5.3. Case 3

When acquiring the RELM mapping model between the capacitance correlation coefficient and the imaging target, it is first necessary to confirm the number of hidden layer nodes (NHLNs). In the actual numerical simulation and experimental process, it is found that the NHLN has an impact on the quality of imaging results. To explore the mechanism of interaction between image accuracy and the NHLN, a series of simulation experiments are carried out in this section, focusing on the relationship between the NHLN and the image errors of the imaging results. The relationship between image errors and NHLN is intuitively illustrated by numerical methods. The reconstructed targets are shown in Figure 7, and the image error as a function of the NHLN is shown in Figure 12.
From Figure 12, we observe that the image errors of the seven kinds of phantom objects decrease rapidly with the increase of the NHLN. When the NHLN is greater than 700, the image errors change very little. However, it is worth noting that although increasing the NHLN can increase the imaging quality of the RELM algorithm, it also increases the computational burden. Especially when the training sample set is large, it increases the time cost significantly. Therefore, in the actual ECT application, the relationship between image errors and the NHLN should be balanced, and the appropriate NHLN should be selected according to the actual needs.

6. Experiment Results and Analysis

A series of experiments are carried out to validate the performance of the proposed method in practical applications. As shown in Figure 13, the ECT imaging system is built by our research team and mainly consists of the following parts: (1) the sensor device. This part is a 12-electrode square sensor with a copper foil wrapped cylindrical shield on the outside, whose cross section is shown in Figure 7. The size parameters of the sensor are presented in Table 1. (2) Data acquisition system. The data acquisition system can switch up to 12 channels, and its resolution is 0.5 Ff. (3) Imaging system. In this experiment, quartz sand is used as the measured material with high dielectric constant, and air is used as the background material with low dielectric constant. The measured objects are plastic polycarbonate tubes filled with quartz sand of different diameters, as shown in (4) in Figure 13. We implement three sets of experiments, and their two-dimensional imaging targets are shown as (a)–(c) in Figure 14. These targets include a concentric rectangular ring (the side length of the large square = 40 mm, and the thickness = 15 mm), big circle (diameter = 16 mm), and small circle (diameter = 6 mm). The empty field calibration and the full field calibration are performed by filling the square tube with air and quartz sand, respectively. For the imaging targets (a)–(c) in Figure 14, the NHLN in the ELM algorithm and the RELM algorithm is 1000, and the parameters of various iterative algorithms are listed in Table 8, Table 9 and Table 10. Table 11 and Table 12 are the images reconstructed of various imaging algorithms, and we quantitatively compare the imaging quality using REs, as presented in Figure 15.
From the image results in Table 11 and Table 12, we can see that the imaging quality of the iterative imaging algorithm is higher than that of the non-iterative imaging algorithm, but the distortion and deformation are still large. From Table 12, we can observe that the GPR algorithm and the ELM algorithm greatly improve the imaging quality, but their imaging results have a certain degree of artifacts due to its poor stability. The results in Table 12 indicate that the RELM algorithm provides a significant improvement for the four kinds of imaging targets, and they are better than other conventional imaging algorithms. In Figure 15, we can clearly see that the image errors of the RELM algorithm are the smallest for the four kinds of experiments, which are 1.43%, 2.54%, and 2.39%, respectively. From experimental results, we conclude that the RELM algorithm effectively solves the image reconstruction problem of the ECT and obtains satisfactory imaging results.

7. Conclusions

To alleviate the ill-conditioned characteristics of the ECT inverse problem, an image reconstruction algorithm based on the RELM is proposed. The main conclusions of this paper are summarized as follows:
(1)
Combined with the prior information of the imaging targets, the RELM algorithm improves the imaging quality. The implementation of the algorithm mainly includes two stages: at the first stage, according to a large number of training samples, the RELM model solved by the ISB algorithm between the capacitance correlation coefficient and the imaging target is extracted; at the second stage, the correlation coefficient of the measured capacitance of ECT system is calculated, which is then used as input for the RELM model to predict the final imaging.
(2)
The solution method of capacitance correlation coefficient is introduced. The newly designed objective function based on L1-norm and STV norm can better model the inverse problem of the ECT, which alleviates the adverse effects caused by the inaccuracy of capacitance data and enhances the sparse characteristics of imaging targets. The ISB method greatly improves the solution efficiency for the objective function.
(3)
Numerical simulation results verify the effectiveness of the RELM algorithm. For the seven kinds of complex flow patterns shown in Figure 7, the RELM algorithm achieves the best imaging quality and the smallest image errors compared to the other above-mentioned algorithms. Under different noise levels, the RELM algorithm has strong robustness and noise resistance. When the number of hidden layer nodes of the RELM is greater than 700, the image error remains at a low level and the attenuation degree is small.
(4)
In this study, a set of the ECT system is built to realize on-line imaging. The experimental results proved that the RELM algorithm has feasibility and effectiveness in the actual ECT image reconstruction process.
This paper proposes a new ECT imaging algorithm. Its advantages are verified by numerical simulations and experiments. The high-precision imaging quality of the RELM algorithm provides a reliable basis for studying the internal mechanism of the imaging target, which should be further studied in the future.

Author Contributions

Conceptualization, Q.Z. and S.L.; methodology, Q.Z.; software, Q.Z.; validation, Q.Z. and S.L.; resources, S.L.; writing—original draft preparation, Q.Z.; writing—review and editing, Q.Z. and S.L.; visualization, W.C.; supervision, W.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant No. 61871181.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article. Further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Gut, Z. Selected application of electrical capacitance tomography in monitoring of combustion process. Trans. Inst. Aviat. 2016, 4, 26–36. [Google Scholar] [CrossRef]
  2. Wen, Y.; Yao, J.; Zhang, Y.; Luo, X.; Wang, H. Defect Detection of Adhesive Layer of Thermal Insulation Materials Based on Improved Particle Swarm Optimization of ECT. Sensor 2017, 17, 2440. [Google Scholar] [CrossRef] [PubMed]
  3. Sun, M.; Liu, S.; Li, Z.H.; Lei, J. Application of Electrical Capacitance Tomography to the Concentration Measurement in a Cyclone Dipleg. Chin. J. Chem. Eng. 2008, 16, 635–639. [Google Scholar] [CrossRef]
  4. Pu, W.; Zhao, C.; Xiong, Y.; Cai, L.; Chen, X.; Peng, L.; Fan, C. Numerical simulation on dense phase pneumatic conveying of pulverized coal in horizontal pipe at high pressure. J. Chem. Ind. Eng. 2008, 65, 2500–2512. [Google Scholar] [CrossRef]
  5. Liu, Z.; Wang, H.; Yang, W.; Ge, R. Optimization of large scale concentric-annulus electrical capacitance tomography sensor with internal electrodes. Measurement 2022, 188, 110604. [Google Scholar] [CrossRef]
  6. Ripon, K.; Kwong, S.; Man, K.F. A real-coding jumping gene genetic algorithm (RJGGA) for multiobjective optimization. Inf. Sci. 2007, 177, 632–654. [Google Scholar] [CrossRef]
  7. Yan, C.; Zhang, D.; Lu, G.; Dong, J. An Improved Algorithm Based on Landweber-Tikhonov Alternating Iteration for ECT Image Reconstruction. J. Phys. Conf. Ser. 2018, 1069, 012178. [Google Scholar] [CrossRef]
  8. Nagares, N.M.; Sison, A.M. A High-Performance Computing of Internal Rate of Return using a Centroid-based Newton-Raphson Iterative Algorithm. Int. J. Adv. Trends Comput. Sci. Eng. 2020, 9, 1888–1892. [Google Scholar] [CrossRef]
  9. Deng, H.B.; Zhang, L.; Ying, W.U.; Zhou, J.; Liu, F. Research on track estimation based on Kalman filtering algorithm. Transducer. Microsyst. Technol. 2012, 31, 4–7. [Google Scholar]
  10. Jun, G.U.; Zou, Q.; Deng, C.; Wang, X. A Novel Robust Online Extreme Learning Machine for the Non-Gaussian Noise. Chin. J. Electron. 2023, 32, 130–139. [Google Scholar]
  11. Sun, H.; Qiang, W. Least square regression with indefinite kernels and coefficient regularization. Appl. Comput. Harmon. A 2011, 30, 96–109. [Google Scholar] [CrossRef]
  12. Zhao, Q.; Li, J.; Liu, S.; Liu, G.; Liu, J. The Sensitivity Optimization Guided Imaging Method for Electrical Capacitance Tomography. IEEE Trans. Instrum. Meas. 2021, 70, 4506215. [Google Scholar] [CrossRef]
  13. Kang, Y.; Liu, S.; Liu, J. Image reconstruction algorithm for electrical capacitance tomography based on data correlation analysis. Flow. Meas. Instrum. 2018, 62, 113–122. [Google Scholar] [CrossRef]
  14. Rui, Z.; Yuan, L.; Huang, G.B.; Xu, Z.B. Universal Approximation of Extreme Learning Machine with Adaptive Growth of Hidden Nodes. IEEE Trans. Neural Netw. Learn. Syst. 2012, 23, 365. [Google Scholar]
  15. Huang, G.B.; Ding, X.; Zhou, H. Optimization method based extreme learning machine for classification. Neurocomputing 2010, 74, 155–163. [Google Scholar] [CrossRef]
  16. Cao, J.; Lin, Z.; Huang, G.B.; Nan, L. Voting based extreme learning machine. Basic Sci. J. Text. Univ. 2012, 185, 66–77. [Google Scholar] [CrossRef]
  17. Liu, X.; Wang, X.; Hu, H.; Li, L.; Yang, X.J. An extreme learning machine combined with Landweber iteration algorithm for the inverse problem of electrical capacitance tomography. Flow Meas. Instrum. 2015, 45, 348–356. [Google Scholar] [CrossRef]
  18. Zhao, Q.; Liu, S.; Chai, X.; Guo, H. A novel computational imaging algorithm based on split Bregman iterative for electrical capacitance tomography. Meas. Sci. Technol. 2021, 32, 125401. [Google Scholar] [CrossRef]
  19. Suo, P.; Sun, J.; Tian, W.; Suna, S.; Xu, L. 3D Image Reconstruction in Planar Array ECT by Combining Depth Estimation and Sparse Representation. IEEE Trans. Instrum. Meas. 2021, 70, 4504909. [Google Scholar] [CrossRef]
  20. Pitao, W.; Huaxiang, W.; Benyuan, S.J. l1-norm-based Image Reconstruction Algorithm for Electrical Capacitance Tomography. Proc. CSEE 2015, 18, 018. [Google Scholar]
  21. Papafitsoros, K.; Schönlieb, C.B. A combined first and second order variational approach for image reconstruction. J. Math. Imaging Vis. 2014, 48, 308–338. [Google Scholar] [CrossRef]
  22. Li, K.; Chandrasekera, T.C.; Li, Y.; Holland, D.J.J. A Non-Linear Reweighted Total Variation Image Reconstruction Algorithm for Electrical Capacitance Tomography. IEEE Sens. J. 2018, 18, 5049–5057. [Google Scholar] [CrossRef]
  23. Jing, W.; Ma, J.; Bo, H.; Qin, L. Split Bregman iterative algorithm for sparse reconstruction of electrical impedance tomography. Signal. Process. 2012, 92, 2952–2961. [Google Scholar]
  24. Liu, X.; Huang, L. Split Bregman iteration algorithm for total bounded variation regularization based image deblurring. J. Math. Anal. Appl. 2010, 372, 486–495. [Google Scholar] [CrossRef]
  25. Beck, A.; Teboulle, M. A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems. SIAM J. Imaging Sci. 2009, 2, 183–202. [Google Scholar] [CrossRef]
Figure 1. Imaging procedure of general extreme learning machine.
Figure 1. Imaging procedure of general extreme learning machine.
Applsci 14 00587 g001
Figure 2. Imaging procedure of regularized extreme learning machine.
Figure 2. Imaging procedure of regularized extreme learning machine.
Applsci 14 00587 g002
Figure 3. Network structure diagram of the ELM.
Figure 3. Network structure diagram of the ELM.
Applsci 14 00587 g003
Figure 4. Calculation flow chart of output weight matrix V .
Figure 4. Calculation flow chart of output weight matrix V .
Applsci 14 00587 g004
Figure 5. Program flowchart of the RELM algorithm.
Figure 5. Program flowchart of the RELM algorithm.
Applsci 14 00587 g005
Figure 6. ECT sensor structure.
Figure 6. ECT sensor structure.
Applsci 14 00587 g006
Figure 7. Phantom objects to be reconstructed.
Figure 7. Phantom objects to be reconstructed.
Applsci 14 00587 g007
Figure 8. Reconstruction results of the non-iterative algorithms. ((a1g1) by LBP algorithm, (a2g2) by TRM algorithm, and (a3g3) by TSVD algorithm for the phantom objects in Figure 7a–g).
Figure 8. Reconstruction results of the non-iterative algorithms. ((a1g1) by LBP algorithm, (a2g2) by TRM algorithm, and (a3g3) by TSVD algorithm for the phantom objects in Figure 7a–g).
Applsci 14 00587 g008
Figure 9. Reconstruction results of the iterative algorithms. ((a1g1) by the CG algorithm, (a2g2) by the Lanweber algorithm, (a3g3) by the SIRT algorithm, (a4g4) by the GPR algorithm, (a5g5) by the ELM algorithm, and (a6g6) by the RELM algorithm for the phantom objects in Figure 7a–g).
Figure 9. Reconstruction results of the iterative algorithms. ((a1g1) by the CG algorithm, (a2g2) by the Lanweber algorithm, (a3g3) by the SIRT algorithm, (a4g4) by the GPR algorithm, (a5g5) by the ELM algorithm, and (a6g6) by the RELM algorithm for the phantom objects in Figure 7a–g).
Applsci 14 00587 g009
Figure 10. Results reconstructed by the RELM algorithm when noise level is 5%.
Figure 10. Results reconstructed by the RELM algorithm when noise level is 5%.
Applsci 14 00587 g010
Figure 11. Results reconstructed by the RELM algorithm when noise level is 10%.
Figure 11. Results reconstructed by the RELM algorithm when noise level is 10%.
Applsci 14 00587 g011
Figure 12. Imaging errors under different numbers of hidden layer neuron (a–g correspond to the imaging errors for the imaging targets in Figure 7a–g).
Figure 12. Imaging errors under different numbers of hidden layer neuron (a–g correspond to the imaging errors for the imaging targets in Figure 7a–g).
Applsci 14 00587 g012
Figure 13. Experimental set-up. (1) Sensor device, (2) data acquisition system, (3) imaging system, (4) measured objects.
Figure 13. Experimental set-up. (1) Sensor device, (2) data acquisition system, (3) imaging system, (4) measured objects.
Applsci 14 00587 g013
Figure 14. Imaging targets.
Figure 14. Imaging targets.
Applsci 14 00587 g014
Figure 15. Comparison of image errors ((ac) correspond to the image errors for the phantom objects in Figure 14a–c).
Figure 15. Comparison of image errors ((ac) correspond to the image errors for the phantom objects in Figure 14a–c).
Applsci 14 00587 g015
Table 1. Structural parameters of the ECT sensor.
Table 1. Structural parameters of the ECT sensor.
ParameterSize
R180 mm
d112 mm
d29 mm
d360 mm
d463 mm
Table 2. Training sample pairs.
Table 2. Training sample pairs.
(a)(b)(c)(d)(e)(f)(g)
input γ a γ b γ c γ d γ e γ f γ g Applsci 14 00587 i001
outputApplsci 14 00587 i002Applsci 14 00587 i003Applsci 14 00587 i004Applsci 14 00587 i005Applsci 14 00587 i006Applsci 14 00587 i007Applsci 14 00587 i008
Table 3. Parameters of the CG algorithm.
Table 3. Parameters of the CG algorithm.
Test Objects(a)(b)(c)(d)(e)(f)(g)
Number of iterations7892159243298413295
Relax factor0.030.030.030.030.030.030.03
Table 4. Parameters of the Landweber algorithm.
Table 4. Parameters of the Landweber algorithm.
Test Objects(a)(b)(c)(d)(e)(f)(g)
Number of iterations406475182390549343
Relax factor1111111
Table 5. Parameters of the SIRT.
Table 5. Parameters of the SIRT.
Test Objects(a)(b)(c)(d)(e)(f)(g)
Number of iterations80105223416786424448
Relax factor0.040.040.040.040.040.040.04
Table 6. Imaging errors (%).
Table 6. Imaging errors (%).
(a)(b)(c)(d)(e)(f)(g)
LBP40.4456.6041.4745.7038.4546.8857.41
TR32.2151.9971.7892.0959.9242.6757.19
TSVD23.6447.1172.15121.4944.3137.7044.20
CG24.6725.5029.1330.6932.1525.3920.61
Landweber22.5421.7827.9930.6635.1128.9926.24
SIRT20.0521.0827.5129.8532.3124.5428.07
GPR9.068.799.648.6910.317.627.01
ELM9.817.638.457.618.787.937.19
RELM2.132.061.091.681.491.841.31
Table 7. Image error mean (%).
Table 7. Image error mean (%).
(a)(b)(c)(d)(e)(f)(g)
0% noise2.132.061.091.681.491.841.31
5% noise3.783.562.692.162.642.932.27
10% noise4.514.164.793.614.273.793.43
Table 8. Parameters of the CG algorithm.
Table 8. Parameters of the CG algorithm.
Test Objects(a)(b)(c)
Number of iterations27486108
Relax factor0.030.030.03
Table 9. Parameters of the Landweber algorithm.
Table 9. Parameters of the Landweber algorithm.
Test Objects(a)(b)(c)
Number of iterations21679127
Relax factor111
Table 10. Parameters of the SIRT.
Table 10. Parameters of the SIRT.
Test Objects(a)(b)(c)
Number of iterations52782149
Relax factor0.040.040.04
Table 11. Experimental results of non-iterative algorithms.
Table 11. Experimental results of non-iterative algorithms.
(a)(b)(c)
LBPApplsci 14 00587 i009Applsci 14 00587 i010Applsci 14 00587 i011Applsci 14 00587 i012
TRApplsci 14 00587 i013Applsci 14 00587 i014Applsci 14 00587 i015
TSVDApplsci 14 00587 i016Applsci 14 00587 i017Applsci 14 00587 i018
Table 12. Experimental results of iterative algorithms.
Table 12. Experimental results of iterative algorithms.
CGLandweberSIRTGPRELMRELM
(a)Applsci 14 00587 i019Applsci 14 00587 i020Applsci 14 00587 i021Applsci 14 00587 i022Applsci 14 00587 i023Applsci 14 00587 i024Applsci 14 00587 i025
(b)Applsci 14 00587 i026Applsci 14 00587 i027Applsci 14 00587 i028Applsci 14 00587 i029Applsci 14 00587 i030Applsci 14 00587 i031
(c)Applsci 14 00587 i032Applsci 14 00587 i033Applsci 14 00587 i034Applsci 14 00587 i035Applsci 14 00587 i036Applsci 14 00587 i037
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, Q.; Liu, S.; Chen, W. A Novel Computational Imaging Algorithm for Electrical Capacitance Tomography. Appl. Sci. 2024, 14, 587. https://doi.org/10.3390/app14020587

AMA Style

Zhao Q, Liu S, Chen W. A Novel Computational Imaging Algorithm for Electrical Capacitance Tomography. Applied Sciences. 2024; 14(2):587. https://doi.org/10.3390/app14020587

Chicago/Turabian Style

Zhao, Qing, Shi Liu, and Weining Chen. 2024. "A Novel Computational Imaging Algorithm for Electrical Capacitance Tomography" Applied Sciences 14, no. 2: 587. https://doi.org/10.3390/app14020587

APA Style

Zhao, Q., Liu, S., & Chen, W. (2024). A Novel Computational Imaging Algorithm for Electrical Capacitance Tomography. Applied Sciences, 14(2), 587. https://doi.org/10.3390/app14020587

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