A fast reconstruction algorithm for electron microscope tomography.

We have implemented a Fast Fourier Summation algorithm for tomographic reconstruction of three-dimensional biological data sets obtained via transmission electron microscopy. We designed the fast algorithm to reproduce results obtained by the direct summation algorithm (also known as filtered or R-weighted backprojection). For two-dimensional images, the new algorithm scales as O(N(theta)M log M)+O(MN log N) operations, where N(theta) is the number of projection angles and M x N is the size of the reconstructed image. Three-dimensional reconstructions are constructed from sequences of two-dimensional reconstructions. We demonstrate the algorithm on real data sets. For typical sizes of data sets, the new algorithm is 1.5-2.5 times faster than using direct summation in the space domain. The speed advantage is even greater as the size of the data sets grows. The new algorithm allows us to use higher order spline interpolation of the data without additional computational cost. The algorithm has been incorporated into a commonly used package for tomographic reconstruction.

[1]  G. Herman,et al.  Linograms in Image Reconstruction from Projections , 1987, IEEE Transactions on Medical Imaging.

[2]  Yoram Bresler,et al.  O(N2log2N) filtered backprojection reconstruction algorithm for tomography , 2000, IEEE Trans. Image Process..

[3]  S. Deans.  The Radon Transform and Some of Its Applications , 1983 .

[4]  A. Cormack.  Representation of a Function by Its Line Integrals, with Some Radiological Applications , 1963 .

[5]  P. Gilbert.  The reconstruction of three-dimensional structure from projections and its application to electron microscopy II. Direct methods , 1972, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[6]  Johan Waldén,et al.  Analysis of the direct Fourier method for computer tomography , 2000, IEEE Transactions on Medical Imaging.

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

[8]  I. J. Schoenberg.  Cardinal Spline Interpolation , 1987 .

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

[10]  R. Veselska,et al.  Reaction of the skin fibroblast cytoskeleton to micromanipulation interventions. , 2001, Journal of structural biology.

[11]  Salvatore Lanzavecchia,et al.  Fast computation of 3D radon transform via a direct Fourier method , 1998, Bioinform..

[12]  J. Frank.  Electron Tomography , 1992, Springer US.

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

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

[15]  D. Mastronarde.  Dual-axis tomography: an approach with alignment methods that preserve resolution. , 1997, Journal of structural biology.

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

[17]  F. Natterer.  The Mathematics of Computerized Tomography , 1986 .

[18]  D. J. De Rosier,et al.  Reconstruction of Three Dimensional Structures from Electron Micrographs , 1968, Nature.

[19]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[20]  Achi Brandt,et al.  A Fast and Accurate Multilevel Inversion of the Radon Transform , 1999, SIAM J. Appl. Math..

[21]  J R Kremer,et al.  Computer visualization of three-dimensional image data using IMOD. , 1996, Journal of structural biology.

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

[23]  Gabriele Steidl,et al.  A new linogram algorithm for computerized tomography , 2001 .

[24]  Frank Natterer,et al.  Mathematical methods in image reconstruction , 2001, SIAM monographs on mathematical modeling and computation.