Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 154630, 15 pages

doi:10.1155/2012/154630

Research Article

Combined Energy Minimization for Image Reconstruction from Few Views

Wei Wei,¹ Xiao-Lin Yang,² Bin Zhou,³ Jun Feng,² and Pei-Yi Shen⁴

¹School of Computer Science and Engineering, Xi’an University of Technology, Xi’an 710048, China

²College of Management Science, Chengdu University of Technology, Chengdu 610059, China

³College of Science, Xi’an University of Science and Technology, Xi’an 710054, China

⁴National School of Software, Xidian University, Xi’an 710071, China

Received 15 June 2012; Revised 27 August 2012; Accepted 22 September 2012

Academic Editor: Kui Fu Chen

Copyright © 2012 Wei Wei 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.

Reconstruction from few views is an important problem in medical imaging and applied mathematics. In this paper, a combined energy minimization is proposed for image reconstruction. energy of the image gradient is introduced in the lower density region, and it can accelerate the reconstruction speed and improve the results. Total variation of the image is introduced in the higher density region, and the image features can be preserved well. Nonlinear conjugate gradient method is introduced to solve the problem. The efficiency and accuracy of our method are shown in several numerical experiments.

1. Introduction

Computed tomography (CT) is one of the most important advance in diagnostic radiology in recent decades. CT uses multiple X-ray images to build up cross-sectional and 3D pictures of structures inside the human body which enable doctors to view internal organs with unprecedented precision. However, the use of ionizing radiation in CT may induce cancer in the exposed individual after a latent period [1–3]. Cancer induction by ionizing radiation is a probabilistic process. Reduction of radiation dose used in CT will therefore lead to a reduction in the number of induced cancer cases.

Some ways can be used to reduce the radiation dose from CT such as decreasing intensity of X-ray beam, handling scattered radiation, restricting exposure area. Reducing the X-ray exposing time is a simple one. Here we focus on the low-dose X-ray imaging strategy that only a limited number of projection images are taken for the reconstruction, which is called limited-view reconstruction [4–6]. Image reconstruction from few views would enable rapid scanning with a reduced X-ray dose delivered to the patient. But it is possibe to build up blurring artifacts and lose critical spatial resolution in the reconstructed image. So the trade-off needs to be carefully examined for intended use. However, it is beyond the scope of this paper. In a word, we will consider image reconstruction from projection data at few views.

CT reconstruction methods can roughly be categorized as analytic reconstruction methods and iterative reconstruction methods. The analytic reconstruction methods, such as filtered back-projection (FBP) methods [7, 8], require sufficient projection data with low noise level. As the limited-view reconstruction is considered, the analytic methods may induce more noise and produce significant artifacts. The iterative reconstruction methods, such as the algebraic reconstruction technique (ART) [9, 10], require less data than FBP methods and are more robust to the effects of noise, but need much more computation.

Recently, the minimization of the image total variation (TV) has been introduced to divergent-beam CT and an new iterative image reconstruction algorithm was presented [11]. Lustig et al. applied compressed sensing theory to rapid magnetic resonance imaging [12]. Many approaches have been presented based on this [13–16].

In this paper, a novel image reconstruction model is proposed. The image total variation and the energy of the image gradient are combined to a new energy functional. Then the functional is introduced to the constrained optimization problem for the reconstruction from 2D parallel-beam data at few views. Our algorithm are performed with various insufficient data problems in fan-beam CT and the numerical results show the efficiency and accuracy of proposed method. The algorithm can be generalized to fan-beam CT and cone-beam CT as well as other tomographic imaging modalities.

2. Reconstruction from Few Views and Total Variation Minimization

There are many approaches about tomographic reconstruction from limited views projection data [17–19]. Algebraic reconstruction technique (ART) [9, 20, 21] and the expectation-maximization (EM) algorithm [22, 23] have been widely used in this field. As the image is discretized on the grids, each projection is regarded as a linear equation of the discrete density distribution. Then a system of simultaneous equations can be obtained and ART tends to solve it via iterative method. ART algorithm can find the image that is consistent with the projection data and the sum-of-squares of the density values is minimized. The EM algorithm applies to positive integral equations, seeking to minimize the Kullback-Liebler distance between the measured data and the projection of the estimated image [11].

However, it is known that the ray does harm to human body and abundant irradiation may lead to cancer [17, 24]. So researchers begin to study the tomographic reconstruction with projection data as little as possible.

2.1. Reconstruction from Few Views

Tomographic reconstruction from few views projection data is an efficient way to reduce the harm caused by ray irradiation, and there are some approaches about it [11, 12].

As the gray image to be reconstructed can be denoted by below

Here and mean the size of image.

The projection can be denoted as the following equations:

Here is the vector form rearranged form . means the projection data. More exactly, is the product of the number of views and the number of detector’s pixels. is the projection matrix which can be precomputed.   is the same size as . The reconstruction procedure equals to solve (2.2).

Unfortunately, this equations are indeterminate if the reconstruction was based on few views. In other words, the number of the equations are less than the number of variables (). In practice, it is more often that . From the linear algebraic theory, the solution is not unique and the traditional methods cannot be applied. In fact, even if is satisfied, it is still compromised to deal with the consistency of the projection data and lead to artifacts in the reconstructed image.

The ART can be applied to solve this equation and it means to solve the following problem:

Here means the norm of . Because of the serious insufficiency of the projection data, ART algorithm can hardly provide satisfactory result. The same as to EM algorithm. So some other models should be discussed.

2.2. Total Variation Minimization

The total variation (TV) was first introduced by Rudin et al. [25], and it can be utilized in image processing for images denoising while edges preserved. Candès et al. applied it to image reconstruction with insufficient parallel-beam data [26]. More exactly, they considered the following problem:

Here means the norm of which is the rearranged form of matrix . satisfy . can be computed as

Based on these, Sidky et al. developed an iterative image reconstruction algorithm for fan-beam CT in [11].

The TV minimization can efficiently reduce errors and preserve features in the image reconstruction. In next section, we will concentrate on developing a new model to improve the convergence speed and reduce errors based on this TV model.

3. Minimization of a Combined Energy for Image Reconstruction

It is known that the convergence speed will be enhanced when the norm of image gradient is considered as shown in following:

But this result can also blur the image features. So some combined energies can be considered.

3.1. A Combined Energy of Image

The natural idea is to combine the norm and TV directly. The combined energy can be denoted by

However, the new result cannot be improved much more than the TV result though the convergence speed may be accelerated in some sense. In fact, the new result is a weighted sum of the result and TV result.

To improve the TV result of image denoising, Chambolle and Lions [27] proposed a combined functional (CL energy)

where

Here is a fixed positive number and it is a threshold of . In some way, it means an approximation of the critical value which can be used to distinguish image features and noise.

3.2. The Reconstruction Model

Based on the discrete form of the CL energy and the rearranged vector , we can get

Here , . It can be found that the energy of image gradient is considered in the noise part while the TV energy is computed in features part . Then the new model for image reconstruction can be denoted as follows:

With the Lagrange method applied, this constraint optimization problem can be rewritten as an unconstraint optimization problem of following combined Chambolle-Lions (CCL) energy:

3.3. Conjugate Gradient Descend Algorithm

The gradient of can be computed as

where

In practice, the parameter is set to be .

Set a initial value , the conjugate gradient descend algorithm can be given as Algorithm 1. There the time step is denoted by .

% Initialization

maxGrad = ; maxIter = 100; maxTau = ; = 4; = 0;

= 1; = ; = −;

% Iterations

while ( > maxGrad and < maxIter and > maxTau){

        %Linear search

        min = (); opt = 0; = −;

        while ( <= s){

                ;

                if (( + Δ) < min){min = ( + Δ); opt = ; }

                ; }

         = opt;

         = ; = (); = ;

         = − + ;

         = + 1; }

Algorithm 1: Iteration algorithm for CL minimization reconstruction.

From the experiments results, we will find some advantages of the proposed model. These are chiefly due to the different optimization problems and it is related with the algorithm. More exactly, the first term in (3.5) can help to enhance the convergence speed and reduce some artifacts in the smooth region. But the efficiency role of this energy of image gradient is depended on the well define of . It is related with the characters of the image to be reconstructed. In our experiments, it is set to be 0.01 times the range of phantoms. The advanced researches about this will be approached in our next work. The general metric peak signal to noise ratio (PSNR) is introduced to evaluate the results and it can be computed as

where MSE is the mean square error of the gray image.

4. Numerical Experiments

Example 4.1. Reconstruction of Shepp-Logan phantom from 72 views.

The true image is taken to be the Shepp-Logan image shown in Figure 1(a) discretized on a pixel grid. The computational parameters are set as shown in the algorithm, the same to the following experiments. The reconstruction from 72 views is completed after 31 iterations. Figures 1(b) and 1(c) show the ART result and TV result while Figure 1(d) shows our CL result. It can be found that the ART result is enhanced by TV much more. Many artifacts have been removed or slighted. Figures 1(e) and 1(f) show the gray distributions of row 128 and column 128. There are few differences between the reconstructed horizontal gray and the real one. It is similar to the reconstructed vertical gray. The evolutions of PSNR and are shown in Figures 2(a) and 2(b). The PSNR has been improved from 46.4040 (TV result) to 50.5664 (CL result). Though the TV result is almost accurate, our CL result improves it significantly.

(a) Original

(b) Art result

(c) TV result

(d) CL result

(e) Horizontal gray

(f) Vertical gray

Figure 1: Reconstruction of Shepp-Logan phantom from 72 views.

(a) Evolution of PSNR

(b) Evolution of

Figure 2: Result analysis of Example 4.1.

Example 4.2. Reconstruction of Shepp-Logan phantom from 24 views.

The true image is still taken to be the same Shepp-Logan image as Example 4.1 while the reconstruction views are reduced acutely from 72 to 24. This reconstruction is completed after 100 iterations. Figures 3(b), 3(c), and 3(d) show the ART result, TV result and our CL result. The gray distributions of row 128 and column 128 are shown in Figures 3(e) and 3(f). Figures 4(a) and 4(b) show the evolutions of PSNR and . The PSNR has been improved from 21.3451 (TV result) to 34.4123 (CL result). It can be found that there are a lot of serious line artifacts and block artifacts in TV result. Most of these artifacts has been eliminated and the rest has been slighted significantly in our CL result.

(a) Original

(b) Art result

(c) TV result

(d) CL result

(e) Horizontal gray

(f) Vertical gray

Figure 3: Reconstruction of Shepp-Logan phantom from 24 views.

(a) Evolution of PSNR

(b) Evolution of

Figure 4: Result analysis of Example 4.2.

Example 4.3. Reconstruction of fruits image from 72 views.

A fruits image with size is taken to be the true one. The reconstruction views are set to be 72. It costs 25 iterations to finish the reconstruction. Figures 5(b), 5(c), and 5(d) show the ART result, TV result, and our CL result. The gray distributions of row 128 and column 128 are shown in Figures 5(e) and 5(f). Figure 6 shows the evolutions of PSNR and . The PSNR has been improved from 31.8559 (TV result) to 32.2934 (CL result). There are few differences between TV result and our CL result though CL result is more satisfactory than TV result in numerical value.

(a) Original

(b) Art result

(c) TV result

(d) CL result

(e) Horizontal gray

(f) Vertical gray

Figure 5: Reconstruction of fruits from 72 views.

(a) Evolution of PSNR

(b) Evolution of

Figure 6: Result analysis of Example 4.3.

Example 4.4. Reconstruction of fruits image from 30 views.

The same fruits image as Example 4.3 is taken to be the true one while the views are reduced from 72 to 30. This reconstruction is completed after 59 iterations. ART result, TV result, and our CL result are shown in Figures 7(b), 7(c), and 7(d). The gray distributions of row 128 and column 128 are shown in Figures 7(e) and 7(f). The evolutions of PSNR and are shown in Figures 8(a) and 8(b). The PSNR has been improved from 25.9273 (TV result) to 27.3870 (CL result). It can be found that there are many artifacts in TV result and a lot of them have been eliminated or slighted in our CL result.

(a) Original

(b) Art result

(c) TV result

(d) CL result

(e) Horizontal gray

(f) Vertical gray

Figure 7: Reconstruction of fruits from 30 views.

(a) Evolution of PSNR

(b) Evolution of

Figure 8: Result analysis of Example 4.4.

Example 4.5. Reconstruction of a synopsis phantom from 72 views.

A synopsis phantom with size is taken to be the true image. The reconstruction views are set to be 72. After 24 iterations, the reconstruction is finished. Figures 9(b), 9(c), and 9(d) show the ART result, TV result, and our CL result. The gray distributions of row 128 and column 128 are shown in Figures 9(e) and 9(f). Figure 10 shows the evolutions of PSNR and . The PSNR has been improved from 48.5852 (TV result) to 49.6270 (CL result). It can be found that CL result is more satisfactory than TV result in numerical value, but there are few differences between TV result and our CL result from vision terms.

(a) Original

(b) Art result

(c) TV result

(d) CL result

(e) Horizontal gray

(f) Vertical gray

Figure 9: Reconstruction of a synopsis phantom from 72 views.

(a) Evolution of PSNR

(b) Evolution of

Figure 10: Result analysis of Example 4.5.

Example 4.6. Reconstruction of a synopsis phantom from 20 views.

The same synopsis phantom as Example 4.5 is taken to be the true one while the views are reduced from 72 to 20. This reconstruction is completed after 68 iterations. ART result, TV result, and our CL result are shown in Figures 11(c), 11(b), and 11(d). The gray distributions of row 128 and column 128 are shown in Figures 11(e) and 11(f). The evolutions of PSNR and are shown in Figures 12(a) and 12(b). The PSNR has been improved from 31.7194 (TV result) to 36.9629 (CL result). It can be found that there are more artifacts in TV result than in our CL result.

(a) Original

(b) Art result

(c) TV result

(d) CL result

(e) Horizontal gray

(f) Vertical gray

Figure 11: Reconstruction of a synopsis phantom from 20 views.

(a) Evolution of PSNR

(b) Evolution of

Figure 12: Result analysis of Example 4.6.

5. Conclusion

In this paper, a novel model for image reconstruction from few views in parallel-beam data is proposed. First, the energy of the image gradient and the total variation of the image are combined to the CL energy. The energy is applied in the lower density region, and it can accelerate the reconstruction speed. The total variation is applied in the higher density region, and it can preserve the image features well. Contributed to the lagrange method and nonlinear conjugate gradient algorithm, the related optimization problem can be solved.

Acknowledgments

The authors would like to thank the anonymous reviewers for their constructive feedback and valuable input. Due thanks are for the supports to their program from the TI, the XILINX, and the Software School of Xidian University. This program is partially supported by NSFC (Grant nos. 61072105, 61007011) and by the Open Projects Program of National Laboratory of Pattern Recognition. The Project is also partially supported by Natural Science Basic Research Plan in Shaanxi Province of China (Program no. 2010JM8005) and Scientific Research Program Funded by Shaanxi Provincial Education Department (Program no. 11JK0504).

References

D. J. Brenner and E. J. Hall, “Computed tomography—an increasing source of radiation exposure,” New England Journal of Medicine, vol. 357, no. 22, pp. 2277–2284, 2007.

R. E. Van Gelder, H. W. Venema, J. Florie et al., “CT colonography: feasibility of substantial dose reduction—comparison of medium to very low doses in identical patients,” Radiology, vol. 232, no. 2, pp. 611–620, 2004.

N. Savage, “Medical imagers lower the dose,” IEEE Spectrum, vol. 4, no. 3, pp. 14–16, 2010.

R. Rangayyan, A. Prakash Dhawan, and R. Gordon, “Algorithms for limited-view computed tomography: an anotated bibliography and a challenge,” Applied Optics, vol. 24, no. 23, pp. 4000–4012, 1985.

J. Velikina, S. Leng, and G. H. Chen, “Limited view angle tomographic image reconstruction via total variation minimization,” in Medical Imaging: Physics of Medical Imaging, vol. 6510 of Proceedings of SPIE, p. 20, February 2007.

T. Wu, A. Stewart, M. Stanton et al., “Tomographic mammography using a limited number of low-dose cone-beam projection images,” Medical Physics, vol. 30, no. 3, pp. 365–380, 2003.

P. Grangeat, “Mathematical framework of cone beam 3d reconstruction via the first derivative of the radon transform,” in Mathematical Methods in Tomography, pp. 66–97, 1991.

D. Soimu, I. Buliev, and N. Pallikarakis, “Studies on circular isocentric cone-beam trajectories for 3D image reconstructions using FDK algorithm,” Computerized Medical Imaging and Graphics, vol. 32, no. 3, pp. 210–220, 2008.

R. Gordon, R. Bender, and G. T. Herman, “Algebraic Reconstruction Techniques (ART) for three-dimensional electron microscopy and X-ray photography,” Journal of Theoretical Biology, vol. 29, no. 3, pp. 471–481, 1970.

Y. Pan, R. Whitaker, A. Cheryauka, and D. Ferguson, “Tv-regularized iterative image reconstruction on a mobile c-arm ct,” vol. 7622 of Proceedings of SPIE, p. 2L, 2010.

E. Y. Sidky, C. M. Kao, and X. Pan, “Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT,” Journal of X-Ray Science and Technology, vol. 14, no. 2, pp. 119–139, 2006.

M. Lustig, D. Donoho, and J. M. Pauly, “Sparse MRI: the application of compressed sensing for rapid MR imaging,” Magnetic Resonance in Medicine, vol. 58, no. 6, pp. 1182–1195, 2007.

S. Brunner, B. E. Nett, R. Tolakanahalli, and G. H. Chen, “Prior image constrained scatter correction in cone-beam computed tomography image-guided radiation therapy,” Physics in Medicine and Biology, vol. 56, no. 4, pp. 1015–1030, 2011.

G. T. Herman and R. Davidi, “Image reconstruction from a small number of projections,” Inverse Problems, vol. 24, no. 4, pp. 45011–45028, 2008.

J. Song, Q. H. Liu, G. A. Johnson, and C. T. Badea, “Sparseness prior based iterative image reconstruction for retrospectively gated cardiac micro-CT,” Medical Physics, vol. 34, no. 11, pp. 4476–4483, 2007.

J. Tang, B. E. Nett, and G. H. Chen, “Performance comparison between total variation (TV)-based compressed sensing and statistical iterative reconstruction algorithms,” Physics in Medicine and Biology, vol. 54, no. 19, pp. 5781–5804, 2009.

V. Singh, L. Mukherjee, P. M. Dinu, J. Xu, and K. R. Hoffmann, “Limited view CT reconstruction and segmentation via constrained metric labeling,” Computer Vision and Image Understanding, vol. 112, no. 1, pp. 67–80, 2008.

F. Dennerlein, F. Noo, H. Schöndube, G. Lauritsch, and J. Hornegger, “A factorization approach for cone-beam reconstruction on a circular short-scan,” IEEE Transactions on Medical Imaging, vol. 27, no. 7, pp. 887–896, 2008.

B. P. Medoff, W. R. Brody, M. Nassi, and A. Macovski, “Iterative convolution backprojection algorithms for image reconstruction from limited data,” Journal of the Optical Society of America, vol. 73, no. 11, pp. 1493–1500, 1983.

T. G. Herman and B. L. Meyer, “Algebraic reconstruction techniques can be made computationally efficient,” IEEE Transactions on Medical Imaging, vol. 12, no. 3, pp. 600–609, 1993.

F. Nattererm and F. Wübbeling, Mathematical Methods in Image Reconstruction, vol. 5, Society for Industrial Mathematics, 2001.

V. Y. Panin, G. L. Zeng, and G. T. Gullberg, “Total variation regulated EM algorithm,” IEEE Transactions on Nuclear Science, vol. 46, no. 6, pp. 2202–2210, 1999.

H. H. Barrett and K. J. Myers, Foundations of Image Science, vol. 44, Wiley-Interscience, 2004.

B. L. Cohen, “The cancer risk from low level radiation,” in Radiation Dose from Multidetector CT, pp. 61–79, 2012.

L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D, vol. 60, no. 1–4, pp. 259–268, 1992.

E. J. Candès, J. Romberg, and T. Tao, “Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information,” IEEE Transactions on Information Theory, vol. 52, no. 2, pp. 489–509, 2006.

A. Chambolle and P. L. Lions, “Image recovery via total variation minimization and related problems,” Numerische Mathematik, vol. 76, no. 2, pp. 167–188, 1997.