Mathematical Problems in Engineering

Volume 2013, Article ID 128172, 8 pages

http://dx.doi.org/10.1155/2013/128172

## A Superresolution Image Reconstruction Algorithm Based on Landweber in Electrical Capacitance Tomography

Higher Educational Key Laboratory for Measuring and Control Technology and Instrumentations of Heilongjiang Province, Harbin University of Science and Technology, Harbin, China

Received 25 April 2013; Accepted 19 August 2013

Academic Editor: Zhengguang Wu

Copyright © 2013 Chen Deyun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

According to the image reconstruction accuracy influenced by the “soft field” nature and ill-conditioned problems in electrical capacitance tomography, a superresolution image reconstruction algorithm based on Landweber is proposed in the paper, which is based on the working principle of the electrical capacitance tomography system. The method uses the algorithm which is derived by regularization of solutions derived and derives closed solution by fast Fourier transform of the convolution kernel. So, it ensures the certainty of the solution and improves the stability and quality of image reconstruction results. Simulation results show that the imaging precision and real-time imaging of the algorithm are better than Landweber algorithm, and this algorithm proposes a new method for the electrical capacitance tomography image reconstruction algorithm.

#### 1. Introduction

Electrical capacitance tomography (ECT) is a high-tech cross-section detection technology in recent years, and it can achieve the visualization measurement of the multiphase media component parameter in the closed pipe and process container equipment. The ECT technology could be able to obtain the permittivity distribution inside the pipeline through measuring the capacitance values between the different plates of the capacitive sensing electrode array around the pipeline. It has the advantages of noninvasive, simple structure, low cost, fast response, good safety performance, wide application range, and so forth. So, it can be used for measuring the two-phase flow in the petroleum, chemical and electric power industries, and so forth.

There are three principal factors with image reconstruction for the electrical capacitance tomography system. Firstly, it is the “soft field” effect; the sensitivity performance is in homogeneity affected by the medium distribution. Secondly, the measurements are less than the pixels of the image, and the underdetermined equation has numerous solutions. Thirdly, the solution stability is poor; even small perturbation in the measured capacitance data can cause great error. Due to the nonlinear and the “soft field” effect, the stability of the ECT system is poor and seriously pathological. So, it is difficult for image reconstruction [1]. At present, there are many methods about the ECT image reconstruction, such as the linear back projection (LBP), regularization method, Landweber iterative method [2], genetic algorithms [3], projection Landweber iterative method, conjugate gradient method [4], improved Gauss-Newton method [5], and acceleration algorithm based on polynomial [6]. The characteristic of the LBP is simple and fast reconstruction. But because its imaging quality is relatively poor, strictly speaking, the algorithm is only a qualitative algorithm. Regular methods commonly use empirical values, because the regular selection of parameters influences the imaging quality. It is specially worth mentioning that the online preiterative method (OIOR) based on Landweber iterative method has the same imaging quality with several Landweber iterations under the same speed with LBP, but its spatial resolution has yet to exceed Landweber method.

#### 2. Basic Principle of ECT System

Electrical capacitance tomography system mainly consists of three parts: capacitive sensor, data acquisition and signal processing, and image reconstruction, as shown in Figure 1. Pairs of the metal electrode plates are mounted uniformly on the outside of the insulating pipe, and exterior uses shielding cover. This constitutes the capacitive sensor. Each phase separation medium of the multiphase fluid in the pipeline has a different dielectric constant. When the phase concentration and distribution of the fluid change, the equivalent permittivity of the multiphase flow mixture will change. So, the measuring capacitance values between the plates will be changing too. The size of the capacitor values reflects the size and distribution of the multiphase flow medium concentration. The data acquisition and signal-processing section samples the capacitance values between each pair of the plates and obtains the data in the different angular positions after superposition. These data will be sent to the computer after processing. These measurement values reflect the distribution of the dielectric constant in the pipeline. The distribution information of the fluid medium in the pipeline is obtained through some image reconstruction algorithm by the computer.

is the total number of independent electrode pairs, which is available for the ECT system with the -electrode plate. Consider In this paper, the typical 8-electrode capacitance sensor is the research object, we can get 28 independent measured values as . The capacitance value of degree of influence on the same plate pair has different effects on the medium of the different locations in the pipeline. That is the soft field characteristic. On the other hand, the same point of the medium within the pipeline has different effects on the capacitance values of different plate pairs. The capacitance values between a pair of capacitor plates can actually be approximately regarded as the varying degrees superposition by the capacitance of all the points within the pipeline. The capacitance measurement values can be seen as the projection data of the multiphase fluid on a direction within the pipeline. The electrical capacitance tomography obtains each phase parameter of the multiphase flow from the projection data; then reconstruct the image. According to the electrical principles, ignoring the impact of the distribution changes of each subphase of the pipe radial multiphase fluid and the shield, the measure value which is for the capacitance between any two electrodes in the array type of the 8-electrode capacitance sensor is where , is the pipeline cross-sections, represents the distribution function of the dielectric in the pipeline cross-section, and is the sensitivity distribution function of the measuring capacitance , namely, the sensitivity of the capacitance to the media on point .

The sensitivity distribution data is the priori data required for the image reconstruction. The sensitivity distribution which is also known as the sensitive field can reduce the soft field error of the sensitive field and improve the quality of the image reconstruction. It is extremely necessary to carry out in-depth analysis on the sensitivity distribution.

The calculation formula of the sensitivity function is where is empty-pipe capacitance, that is, the capacitance when the pipeline is filled with the discrete phase (oil), and it is a constant, where . When the dielectric constant of the th microelement is in the pipeline and the dielectric constant of other microelements is , is the capacitance value.

#### 3. The Principle of the Landweber Iterative Algorithms

The ECT image reconstruction is a typical ill-posed problem, and its solution is unstable. In order to obtain meaningful results of reconstruction, we should use the method that can improve the quality and stability of the solution. Regularization technique is an effective method to deal with pathological problems, and Landweber iteration method is a typical regularization techniques. The Landweber iterative format can be expressed as [7]

The iterative method actually uses the steepest decent method with the step size as to solve the minimum of the quadratic functional . Actually, if

is the steepest descent iterative sequence with the step size as . On the other hand, constant in the iterative scheme can also be seen as a relaxation factor, and it is a modified form from .

The is expressed as ; the operator is

Use the singular system of compact operator to expand. The expression of is as follows: where the filtering function may be expressed as follows:

When constructing the regularization solution sequence by Landweber iteration, the number of iterations should be given, and it is as the iteration stopping criteria. One common method is that iteration stops when the first which meets formula (9) is disappeared:

The so-called semiconvergence phenomenon will disappear when solves the inverse problem with Landweber iterative method. That is, in the early stages of the iteration approximate solution can be stably improved, and the effect of self-regularization is shown, but the iterations will tend to diverge when iterations exceed a certain threshold [8]. So, the key question is to find a suitable termination principle in order to match the number of iterations and the error level of original data. However, in practical applications, the error level of the original data is difficult to obtain. So, there is a common method to determine the number of iterations by empirical methods.

Landweber iterative method has various modifications and most improvements concentrate on accelerating the problem. From the numerical optimization point, Landweber iterative method belongs to the most descent methods, and its convergence rate is relatively slow. -method is an accelerated method of Landweber algorithm. The method by Brakhage is a polynomial acceleration method whose applications are broader than the method by Chebyshev. -method can be expressed as the following iterative scheme [9]: where ; when , iterative scheme is as follows:

In fact, standard Landweber iterative method can be seen as a special kind of steepest descent method. Standard landweber iterative method steepest descent method have the same iteration format, and the only difference is the selected aspect of the iteration step. The steepest descent method requires searching from each of iteration by the 1-D to obtain the optimal step size. Landweber iteration method makes the step size in the interval .

#### 4. Superresolution Image Reconstruction Algorithm

##### 4.1. Superresolution Reconstruction Model

In the ECT image reconstruction, small changes of capacitance value which is measured can cause the fluctuation of the solution. This can cause the deterioration of the imaging quality. Sequence images about the same object with different angles can make up such issues to a certain extent. The destination of image restoration is to recover image that has been degraded and make sure that the processed image is close to the original as possible Figure 2 simulated the degradation processing of sequence images.

In order to establish mathematical model of super-resolution image reconstruction, suppose that is th low-resolution image, is high-resolution image which needs to be reconstruct, is the geometry transformation matrix, is the fuzzy matrix, is the downsampling matrix, and is the additional noise. We can obtain superresolution reconstruction model which is as follows [10]: where is the vector, which is super-resolution image vectors permutation that we need to solve. shows a low-resolution image with a total of frames. is the vector, which is said to get the low-resolution image whose size is from the -frame. is the multiple of the corresponding super-resolution reconstruction. is the motion compensation matrix of which does pixel registration. And it is obtained by general multilevel block-matching registration. is blurring function matrix . is downsampling matrix . is the Gaussian noise obtained by the low-resolution image of the th frame of the image. Array formula (12) by a column vector of the image and as follows: abbreviated as

##### 4.2. Regularization of Solutions

Due to the arbitrariness of the displacement among the low-resolution images and the number of the reconstruction low-resolution images with the effect of the drop clarity matrix, the previous reconstruction problem is usually a pathological problem. That is, the above reconstruction requires introducing the regularization techniques for solving the ill-posed problem. The minimized formula solved by formula (14) is shown as follows: where is the smooth operator for solving ill-posed problem and it usually uses the Laplacian in super-resolution reconstruction. is the regularization parameter and it can be obtained by experiences in engineering approximation or other methods. Transformation is shown as follows:

Formula (15) becomes

Minimizing the previous equation is equivalent to solving the following equations:

Since the matrix cannot be diagonalized, the solution of the previously mentioned equation is computationally intensive part in the super-resolution image reconstruction. Approximate the registration results among the images by iterative solution algorithm. And solve the previous equation by preconditioned conjugate gradient method. The computational complexity of the algorithm is high, and the computing complexity with the number of low-resolution images increases exponentially. Then, analyze the previous equation, and derive a direct calculation algorithm, so that the complexity of computation increases linearly with the number of low-resolution image.

##### 4.3. The Export of Closed Solution

The solutions of (19) can be directly expressed as follows:

According to the matrix inversion lemma, The equation obtained by formula (20) is shown as follows:

The formula obtained by solution is shown as follows:

According to formula, (14) and (17), is a block matrix, and the block is whose size is . The calculation uses the calculation model of cyclic matrix, so , and are all the block cyclic matrixes. According to the nature of the block cyclic matrix, still is a block cyclic matrix whose size is . Since is the next sampling matrix, the super resolution multiples on the horizontal and vertical directions are the same. is still a block cyclic matrix whose size is . Since block cyclic matrix plus the unit matrix is still a faster cycle matrix, matrix is a partitioned matrix; each block of the matrix is still the block cyclic matrix whose size is . According to the nature that the block cyclic matrix can be diagonalized, let be the block cyclic matrix of : where is the diagonal matrix which is diagonalized by the corresponding block circulant matrix. is the corresponding orthogonal matrix, so the following formula can be obtained:

The corresponding calculation mainly converts to seek the inverse of the corresponding block diagonalization matrix:

Let be an element permutation matrix of block matrix [11]. . The function is to substitute the element of the first block matrix for the element of the block matrix.

So, where every is a matrix. The inverse of formula (27) is shown as follows:

According to the nature and inverse of the block diagonal matrix in formula (29), its inverse is also the block diagonalization of the corresponding matrix inverse. That is,

Thus, the inversion of the matrix whose size is transforms into the inversion of matrix. According to the nature of the permutation matrix, formula (31) is shown as follows: where every is diagonal matrix and substituted into formula (26):

According to formula (14), (16), (17), (23), and (32), formula (33) can be obtained.

Diagonalize the corresponding block circulant matrix:

Substitute formulae (14) and (35) into formula (34):

The simplified formula is shown as follows:

So, the superresolution image reconstruction computable closed formula is derived. In the derivation process, there are many methods to be used, such as the diagonalization of the block circulation diagonal matrix, fast Fourier transform, downsampling matrix transpose multiplying, small matrix inversion, and vector product. The block circulant matrix can be diagonalized by fast Fourier transform of the convolution kernel. Down sampling matrix transpose multiplication is equivalent to add zero in the corresponding position. For the small size of the matrix inversion, it can be calculated directly.

#### 5. Experimental Results and Analysis

In order to verify the effectiveness of the algorithm, do numerical simulation test by the aforementioned algorithm, and compare the reconstruction results with the Landweber iterative algorithm. Matlab simulation is used on a computer with Core (TM) 2 Duo CPU, 2 GHz, and 1 G memory. The ECT system uses the 8-electrode; display uses subdivision, and the imaging region is divided into 1024 pixels; there are 856 imaging units on the effective area of the pipeline cross-section. The research object, respectively, uses the typical stratified flow, core flow, and trickle flow. The image reconstruction results are as shown in Table 1 (Black area is water, and white area is oil).

In order to evaluate the imaging performance of the super resolution image reconstruction algorithm based on Landweber in Electrical Capacitance Tomography, Landweber algorithm is used for comparative evaluation. This paper uses three objective evaluation methods: mean square error, signal to noise ratio, and peak signal to noise ratio. The details are as follows.(1)*Mean Square Error.*

Mean square error indicates the proximity of two images of the same size. Usually, in the super resolution reconstruction experiment, the reconstructed images are compared with the known images after obtaining the low-resolution image sequences using the original images: where are the length and width of the image. represents the pixels of the original image. represents the pixels of the reconstructed image.

The exact figures of the evaluations simulated by Matlab are shown in Table 2.(2)*Signal to Noise Ratio.*

Consider where are the length and width of the image. represents the pixels of the original image.

Signal to noise ratio is used to measure the ratio between the reconstructed image and the original image. The larger the ratio is, the more information the image contains. On the contrary, the smaller the ratio is, the less information the image contains.

The exact figures of the evaluations simulated by Matlab are shown in Table 3.(3)*Peak Signal to Noise Ratio.*

Peak signal to noise ratio indicates the ratio between the maximum possible power signal and the destructive noise power which affects the representation accuracy. The larger the ratio is, the less the distortion is:

The exact figures of the evaluations simulated by Matlab are shown in Table 4.

In Table 1, the quality of the reconstructed images for the flow pattern (a), (b), (c) has improved by using the algorithm proposed in this paper. It can be seen that the resolution in the periphery and center of the imaging area is better than Landweber reconstruction algorithm. The data listed in Tables 2, 3, and 4 are consistent with the results of the corresponding images listed in Table 1, from the MSE, SNR, and PSNR.

#### 6. Conclusions

This paper presents the super-resolution image reconstruction algorithm based on Landweber iteration and gives the mathematical theory calculation steps of the algorithm based on the basic principle of the electrical capacitance tomography. Through the experimental simulation analysis, the algorithm can overcome the “soft field” nature of the electrical capacitance tomography and has a higher imaging accuracy. Experimental results show that the super resolution image reconstruction algorithm based on Landweber iterative can significantly improve imaging quality, stability, and so on. So, it proposes a new method for ECT image reconstruction algorithm.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (60572153, 60972127), Doctoral Program of Higher Education of the Specialized Research Fund (200802140001), The Chunhui Plan of the Education Ministry (Z2007-1-15013), Natural Science Foundation of Heilongjiang Province (F200609, QC2012C059), Project of Higher Educational Key Laboratory for Measuring and Control Technology and Instrumentations of Heilongjiang Province, and Project of Education Department of Heilongjiang Province (11541040, 12511097, 12531094).

#### References

- D. Y. Chen, X. Y. Yin, and L. Q. Sun, “Optimal design and simulation of electrical capacitance sensor for electrical tomography system,”
*Journal of Electronic Measurement and Instrument*, vol. 20, no. 1, pp. 22–27, 2006. View at Google Scholar - W. Q. Yang and L. H. Peng, “Image reconstruction algorithms for electrical capacitance tomography,”
*Measurement Science and Technology*, vol. 14, no. 1, pp. R1–R13, 2003. View at Publisher · View at Google Scholar · View at Scopus - D. Chen, C. J. Yang, G. B. Zheng, X. Yu, and L. Sun, “Novel method of plotting sensor sensitivity field and image reconstruction algorithm for electrical capacitance tomography system,”
*Journal of Scientific Instrument*, vol. 26, no. 3, pp. 229–234, 2005. View at Google Scholar · View at Scopus - H.-X. Wang, X.-M. Zhu, and L.-F. Zhang, “Conjugate gradient algorithm for electrical capacitance tomography,”
*Journal of Tianjin University Science and Technology*, vol. 38, no. 1, pp. 1–4, 2005. View at Google Scholar · View at Scopus - D.-Y. Chen, Y. Chen, L.-L. Wang, and X.-Y. Yu, “A novel Gauss-Newton image reconstruction algorithm for Electrical Capacitance Tomography System,”
*Acta Electronica Sinica*, vol. 4, no. 37, pp. 739–743, 2009. View at Google Scholar · View at Scopus - Y. Chen, D. Y. Chen, and L. L. Wang, “Image reconstruction algorithm accelerated by polynomial for electrical capacitance tomography system,”
*Chinese Journal of Scientific Instrument*, vol. 29, no. 12, pp. 2538–2542, 2008. View at Google Scholar - K. M. Huang and X. Zhao,
*Inverse Problems and Applications of Electromagnetic Fields*, Beijing Science Press, 2005. - T. Y. Xiao, S. G. Yu, and Y. F. Wang,
*Numerical Solution of the Inverse Problem*, Beijing Science Press, 2003. - Y. F. Wang,
*Calculation Method and Application of the Inversion Problems*, Beijing Higher Education Press, 2007. - N. Nguyen, P. Milanfar, and G. H. Golub, “A computationally efficient superresolution image reconstruction algorithm,”
*IEEE Transactions on Image Processing*, vol. 10, no. 4, pp. 573–583, 2001. View at Publisher · View at Google Scholar · View at Scopus - X. D. Zhang,
*Linear Algebra in Signal Processing*, Science Press, Beijing, China, 1997, Chinese.