Gridding-based direct Fourier inversion of the three-dimensional ray transform.

We describe a fast and accurate direct Fourier method for reconstructing a function f of three variables from a number of its parallel beam projections. The main application of our method is in single particle analysis, where the goal is to reconstruct the mass density of a biological macromolecule. Typically, the number of projections is extremely large, and each projection is extremely noisy. The projection directions are random and initially unknown. However, it is possible to determine both the directions and f by an iterative procedure; during each stage of the iteration, one has to solve a reconstruction problem of the type considered here. Our reconstruction algorithm is distinguished from other direct Fourier methods by the use of gridding techniques that provide an efficient means to compute a uniformly sampled version of a function g from a nonuniformly sampled version of Fg, the Fourier transform of g, or vice versa. We apply the two-dimensional reverse gridding method to each available projection of f, the function to be reconstructed, in order to obtain Ff on a special spherical grid. Then we use the three-dimensional gridding method to reconstruct f from this sampled version of Ff. This stage requires a proper weighting of the samples of Ff to compensate for their nonuniform distribution. We use a fast method for computing appropriate weights that exploits the special properties of the spherical sampling grid for Ff and involves the computation of a Voronoi diagram on the unit sphere. We demonstrate the excellent speed and accuracy of our method by using simulated data.

[1]  A Leith,et al.  SPIDER and WEB: processing and visualization of images in 3D electron microscopy and related fields. , 1996, Journal of structural biology.

[2]  N Grigorieff,et al.  Three-dimensional structure of bovine NADH:ubiquinone oxidoreductase (complex I) at 22 A in ice. , 1998, Journal of molecular biology.

[3]  R. H. Wade A brief look at imaging and contrast transfer , 1992 .

[4]  L. I. Matveenko,et al.  The aperture synthesis. , 1983 .

[5]  J Frank,et al.  A NEW 3‐D RECONSTRUCTION SCHEME APPLIED TO THE 50S RIBOSOMAL SUBUNIT OF E. COLI , 1986, Journal of microscopy.

[6]  Lanzavecchia,et al.  Fast and accurate three-dimensional reconstruction from projections with random orientations via radon transforms , 1999, Journal of structural biology.

[7]  Vladimir Rokhlin,et al.  Fast Fourier Transforms for Nonequispaced Data , 1993, SIAM J. Sci. Comput..

[8]  Pawel A Penczek,et al.  Three-dimensional spectral signal-to-noise ratio for a class of reconstruction algorithms. , 2002, Journal of structural biology.

[9]  J. D. O'Sullivan,et al.  A Fast Sinc Function Gridding Algorithm for Fourier Inversion in Computer Tomography , 1985, IEEE Transactions on Medical Imaging.

[10]  W. O. Saxton,et al.  The correlation averaging of a regularly arranged bacterial cell envelope protein , 1982, Journal of microscopy.

[11]  Salvatore Lanzavecchia,et al.  Electron tomography in conical tilt geometry. The accuracy of a direct Fourier method (DFM) and the suppression of non-tomographic noise , 1996 .

[12]  J. Frank Three-Dimensional Electron Microscopy of Macromolecular Assemblies , 2006 .

[13]  J M Carazo,et al.  3D reconstruction in electron microscopy using ART with smooth spherically symmetric volume elements (blobs). , 1998, Ultramicroscopy.

[14]  Jan Timmer,et al.  The gridding method for image reconstruction by Fourier transformation , 1995, IEEE Trans. Medical Imaging.

[15]  Volker Rasche,et al.  Resampling of data between arbitrary grids using convolution interpolation , 1999, IEEE Transactions on Medical Imaging.

[16]  J. Frank,et al.  The ribosome at improved resolution: new techniques for merging and orientation refinement in 3D cryo-electron microscopy of biological particles. , 1994, Ultramicroscopy.

[17]  Robert J. Renka,et al.  Algorithm 772: STRIPACK: Delaunay triangulation and Voronoi diagram on the surface of a sphere , 1997, TOMS.

[18]  J. Gillis,et al.  Methods in Computational Physics , 1964 .

[19]  Ya. A. Ilyushin,et al.  On the theory of three-dimensional reconstruction , 1997 .

[20]  J. Frank,et al.  Solution Structure of the E. coli 70S Ribosome at 11.5 Å Resolution , 2000, Cell.

[21]  G. Beylkin On the Fast Fourier Transform of Functions with Singularities , 1995 .

[22]  A. Macovski,et al.  Selection of a convolution function for Fourier inversion using gridding [computerised tomography application]. , 1991, IEEE transactions on medical imaging.

[23]  Daniel Zwillinger,et al.  CRC standard mathematical tables and formulae; 30th edition , 1995 .

[24]  Salvatore Lanzavecchia,et al.  SPARK, a kernel of software programs for spatial reconstruction in electron microscopy , 1993 .

[25]  Gabriele Steidl,et al.  New Fourier reconstruction algorithms for computerized tomography , 2000, SPIE Optics + Photonics.

[26]  M. Heel,et al.  Exact filters for general geometry three dimensional reconstruction , 1986 .

[27]  J Frank,et al.  Three-dimensional reconstruction of single particles embedded in ice. , 1992, Ultramicroscopy.

[28]  Hermann Schomberg,et al.  Notes on direct and gridding-based Fourier inversion methods , 2002, Proceedings IEEE International Symposium on Biomedical Imaging.