Tomographic Image Reconstruction from Incomplete View Data by Convex Projections and Direct Fourier Inversion

We consider the problem of reconstructing CAT imagery by the direct Fourier method (DFM) when not all view data are available. To restore the missing information we use the method of projections onto convex sets (POCS). POCS is a recursive image restoration technique that finds a solution consistent with the measured data and a priori known constraints in both the space and Fourier domain. Because DFM reconstruction is a frequency-domain technique it is ideally matched to POCS restoration when, for one reason or another, we are forced to generate an image from a less than complete set of view data. We design and apply an algorithm (PRDF) which interpolates/extrapolates the missing Fourier domain information by POCS and reconstructs an image by DFM. A simulated human thorax cross section is restored and reconstructed. The restorations using POCS are compared with the Gerchberg-Papoulis extrapolation method and shown to be superior. Applications of PRDF to other types of medical imaging modalities are discussed.

[1]  Z. Opial.  Weak convergence of the sequence of successive approximations for nonexpansive mappings , 1967 .

[2]  Boris Polyak,et al.  The method of projections for finding the common point of convex sets , 1967 .

[3]  R. M. Mersereau,et al.  Digital reconstruction of multidimensional signals from their projections , 1974 .

[4]  A. Papoulis.  A new algorithm in spectral analysis and band-limited extrapolation. , 1975 .

[5]  P. Dreike,et al.  Convolution reconstruction of fan beam projections , 1976 .

[6]  H. Stark.  Sampling theorems in polar coordinates , 1979 .

[7]  M. Kaveh,et al.  Reconstructive tomography and applications to ultrasonics , 1979, Proceedings of the IEEE.

[8]  R. Hingorani,et al.  An Investigation of Computerized Tomography by Direct Fourier Inversion and Optimum Interpolation , 1981, IEEE Transactions on Biomedical Engineering.

[9]  K. Tam,et al.  Tomographical imaging with limited-angle input , 1981 .

[10]  R. Hingorani,et al.  An Investigation of Computerized Tomography by Direct Fourier Inversion and Optimum Interpolation , 1981 .

[11]  A. Lent,et al.  An iterative method for the extrapolation of band-limited functions , 1981 .

[12]  M. Nakajima,et al.  A Reconstruction Algorithm of Body Sections with Opaque Obstructions , 1982, IEEE Transactions on Sonics and Ultrasonics.

[13]  M. Sezan,et al.  Image Restoration by the Method of Convex Projections: Part 2-Applications and Numerical Results , 1982, IEEE Transactions on Medical Imaging.

[14]  D. Youla,et al.  Image Restoration by the Method of Convex Projections: Part 1ߞTheory , 1982, IEEE Transactions on Medical Imaging.

[15]  J. Greenleaf,et al.  Limited Angle Multifrequency Diffraction Tomography , 1982, IEEE Transactions on Sonics and Ultrasonics.

[16]  Z. Cho,et al.  Fourier transform nuclear magnetic resonance tomographic imaging , 1982, Proceedings of the IEEE.

[17]  Henry Stark,et al.  Comments and corrections on the use of polar sampling theorems in CT , 1983 .

[18]  William R. Brody,et al.  Iterative convolution backprojection algorithms for image reconstruction from limited data , 1983 .

[19]  H Stark,et al.  Image restoration by convex projections in the presence of noise. , 1983, Applied optics.

[20]  J. Greenleaf,et al.  Computerized tomography with ultrasound , 1983, Proceedings of the IEEE.

[21]  R. Lewitt.  Reconstruction algorithms: Transform methods , 1983, Proceedings of the IEEE.

[22]  A. Lent,et al.  An introduction to NMR imaging: From the Bloch equation to the imaging equation , 1983, Proceedings of the IEEE.

[23]  Direct Fourier Reconstruction Techniques in NMR Tomography , 1984 .