Limited angle reconstruction using stabilized algorithms

The stability of solutions to the limited-angle tomography reconstruction problem obtained by using the projections-onto-convex-sets (POCS) technique are examined. Although POCS techniques provide a feasible solution to the reconstruction problem, the solution is only one sample from the intersection of the closed convex sets that define the solution space. A method for evaluating the ensemble of possible solution waveforms that are in the neighborhood of a solution is presented. The ensemble characteristics are used to construct an inverse filter which is then applied to the computed solution. The results obtained using this method are less sensitive to noise amplification and are less dependent on both starting data and the number of iterations. An estimate of the object-dependent extrapolation that is possible using either linear or nonlinear constraints is provided. >

[1]  Tomaso Poggio,et al.  Computational vision and regularization theory , 1985, Nature.

[2]  H. Trussell,et al.  Convergence criteria for iterative restoration methods , 1983 .

[3]  K. Miller.  Least Squares Methods for Ill-Posed Problems with a Prescribed Bound , 1970 .

[4]  M. Sezan,et al.  Tomographic Image Reconstruction from Incomplete View Data by Convex Projections and Direct Fourier Inversion , 1984, IEEE Transactions on Medical Imaging.

[5]  Gabor T. Herman.  Image Reconstruction from Projections , 1979 .

[6]  H Stark,et al.  One-step image reconstruction from incomplete data in computer tomography. , 1989, IEEE transactions on medical imaging.

[7]  R.W. Schafer,et al.  Constrained iterative restoration algorithms , 1981, Proceedings of the IEEE.

[8]  M. E. Davison,et al.  The Ill-Conditioned Nature of the Limited Angle Tomography Problem , 1983 .

[9]  E. P. Tinnel,et al.  Electron microscope tomography: transcription in three dimensions. , 1983, Science.

[10]  A. Devaney,et al.  Inverse Source and Scattering Problems in Ultrasonics , 1983, IEEE Transactions on Sonics and Ultrasonics.

[11]  O. Strand.  Theory and methods related to the singular-function expansion and Landweber's iteration for integral equations of the first kind , 1974 .

[12]  Carazo Jm,et al.  Restoration of direct Fourier three-dimensional reconstructions of crystalline specimens by the method of convex projections. , 1987 .

[13]  M. Reha Civanlar,et al.  Digital signal restoration using fuzzy sets , 1986, IEEE Trans. Acoust. Speech Signal Process..

[14]  R. Bracewell.  Strip Integration in Radio Astronomy , 1956 .

[15]  Joachim Frank,et al.  Three-Dimensional Reconstruction of Nonperiodic Macromolecular Assemblies from Electron Micrographs , 1986 .

[16]  L. Landweber.  An iteration formula for Fredholm integral equations of the first kind , 1951 .

[17]  H. Trussell,et al.  The feasible solution in signal restoration , 1984 .

[18]  José María Carazo,et al.  Information recovery in missing angular data cases: an approach by the convex projections method in three dimensions , 1987 .

[19]  Carl Wunsch,et al.  Ocean acoustic tomography: a scheme for large scale monitoring , 1979 .

[20]  R M Rangayyan,et al.  Image restoration by Wiener deconvolution in limited-view computed tomography. , 1985, Applied optics.

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

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

[23]  A P Dhawan,et al.  Algorithms for limited-view computed tomography: an annotated bibliography and a challenge. , 1985, Applied optics.

[24]  L. Shepp,et al.  Limited Angle Reconstruction in Tomography via Squashing , 1987, IEEE Transactions on Medical Imaging.

[25]  S. Ivansson,et al.  Seismic borehole tomography—Theory and computational methods , 1986, Proceedings of the IEEE.

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