Improved iterative algorithm for sparse object reconstruction and its performance evaluation with micro-CT data

This paper proposes a new row-action type iterative algorithm which is appropriate to reconstruct sparse objects from a limited number of projections. The main idea is to use the L/sub 1/ norm to pick up a sparse solution from a set of feasible solutions to the measurement equation. By perturbing the linear program to a quadratic program, we use the duality of nonlinear programming to construct a row-action type iterative algorithm to find the solution. We also prove that the algorithm converges to the solution under mild assumptions. We show that this method works well in the 3-D blood-vessel reconstruction and its computation time is shorter than those of our previous method and MART method. Furthermore, we apply the method to real data measured with the Micro-CT device developed at Marquette University.

[1]  Erik L Ritman,et al.  The use of microcomputed tomography to study microvasculature in small rodents. , 2002, American journal of physiology. Regulatory, integrative and comparative physiology.

[2]  D W Holdsworth,et al.  Techniques to alleviate the effects of view aliasing artifacts in computed tomography. , 1999, Medical physics.

[3]  Dimitri P. Bertsekas,et al.  Nonlinear Programming , 1997 .

[4]  Françoise J. Prêteux,et al.  Generalized support constraint for three-dimensional reconstruction from incomplete Fourier spectra , 1997, J. Electronic Imaging.

[5]  Nouna Kettaneh,et al.  Statistical Modeling by Wavelets , 1999, Technometrics.

[6]  Peter M. Joseph,et al.  View sampling requirements in fan beam computed tomography. , 1980 .

[7]  Anita H. Dhyani,et al.  Pulmonary arterial morphometry from microfocal X-ray computed tomography. , 2001, American journal of physiology. Heart and circulatory physiology.

[8]  Anita H. Dhyani,et al.  Microfocal X-ray CT imaging and pulmonary arterial distensibility in excised rat lungs. , 2001, American journal of physiology. Heart and circulatory physiology.

[9]  Paul Tseng,et al.  Relaxation methods for problems with strictly convex separable costs and linear constraints , 1987, Math. Program..

[10]  Gabor T. Herman,et al.  Algebraic reconstruction techniques can be made computationally efficient [positron emission tomography application] , 1993, IEEE Trans. Medical Imaging.

[11]  C Picard,et al.  In vivo evaluation of a new system for 3D computerized angiography. , 1994, Physics in medicine and biology.

[12]  Y. Censor.  Finite series-expansion reconstruction methods , 1983, Proceedings of the IEEE.

[13]  Habib Benali,et al.  Hyperfixation point-source reconstruction by the maximum entropy on the mean method , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[14]  O. Mangasarian,et al.  NONLINEAR PERTURBATION OF LINEAR PROGRAMS , 1979 .

[15]  Hiroyuki Kudo,et al.  An accurate iterative reconstruction algorithm for sparse objects: application to 3D blood vessel reconstruction from a limited number of projections. , 2002, Physics in medicine and biology.

[16]  O. Mangasarian.  Iterative Solution of Linear Programs , 1981 .