High-performance electron tomography of complex biological specimens.

We have evaluated reconstruction methods using smooth basis functions in the electron tomography of complex biological specimens. In particular, we have investigated series expansion methods, with special emphasis on parallel computation. Among the methods investigated, the component averaging techniques have proven to be most efficient and have generally shown fast convergence rates. The use of smooth basis functions provides the reconstruction algorithms with an implicit regularization mechanism, very appropriate for noisy conditions. Furthermore, we have applied high-performance computing (HPC) techniques to address the computational requirements demanded by the reconstruction of large volumes. One of the standard techniques in parallel computing, domain decomposition, has yielded an effective computational algorithm which hides the latencies due to interprocessor communication. We present comparisons with weighted back-projection (WBP), one of the standard reconstruction methods in the areas of computational demand and reconstruction quality under noisy conditions. These techniques yield better results, according to objective measures of quality, than the weighted backprojection techniques after a very few iterations. As a consequence, the combination of efficient iterative algorithms and HPC techniques has proven to be well suited to the reconstruction of large biological specimens in electron tomography, yielding solutions in reasonable computation times.

[1]  Willard Miller,et al.  The IMA volumes in mathematics and its applications , 1986 .

[2]  R. Lewitt.  Alternatives to voxels for image representation in iterative reconstruction algorithms , 1992, Physics in medicine and biology.

[3]  W. Baumeister,et al.  Perspectives of molecular and cellular electron tomography. , 1997, Journal of structural biology.

[4]  Joachim Frank,et al.  Electron Tomography , 1992, Springer US.

[5]  T K Narayan,et al.  A methodology for testing for statistically significant differences between fully 3D PET reconstruction algorithms. , 1994, Physics in medicine and biology.

[6]  B. Wilkinson,et al.  A state-wide senior parallel programming course , 1999 .

[7]  Y. Censor,et al.  BICAV: a block-iterative parallel algorithm for sparse systems with pixel-related weighting , 2001, IEEE Transactions on Medical Imaging.

[8]  Yair Censor,et al.  Component averaging: An efficient iterative parallel algorithm for large and sparse unstructured problems , 2001, Parallel Comput..

[9]  Gabor T. Herman,et al.  ALGEBRAIC RECONSTRUCTION TECHNIQUES IN MEDICAL IMAGING , 1998 .

[10]  Samuel Matej,et al.  A comparison of transform and iterative reconstruction techniques for a volume-imaging PET scanner with a large axial acceptance angle , 1994, Proceedings of 1994 IEEE Nuclear Science Symposium - NSS'94.

[11]  J M Carazo,et al.  The effect of overabundant projection directions on 3D reconstruction algorithms. , 2001, Journal of structural biology.

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

[13]  Sérgio Shiguemi Furuie,et al.  Relevance of statistically significant differences between reconstruction algorithms , 1996, IEEE Trans. Image Process..

[14]  Frank Natterer,et al.  Computational Radiology and Imaging , 1999 .

[15]  William Gropp,et al.  Skjellum using mpi: portable parallel programming with the message-passing interface , 1994 .

[16]  R M Lewitt,et al.  Multidimensional digital image representations using generalized Kaiser-Bessel window functions. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[17]  Andrzej Stachurski,et al.  Parallel Optimization: Theory, Algorithms and Applications , 2000, Parallel Distributed Comput. Pract..

[18]  Robert M. Lewitt,et al.  Practical considerations for 3-D image reconstruction using spherically symmetric volume elements , 1996, IEEE Trans. Medical Imaging.

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

[20]  Samuel Matej,et al.  2.5-D simultaneous multislice reconstruction by series expansion methods from Fourier-rebinned PET data , 2000, IEEE Transactions on Medical Imaging.

[21]  M Marko,et al.  The Emergence of Electron Tomography as an Important Tool for Investigating Cellular Ultrastructure , 2001, The journal of histochemistry and cytochemistry : official journal of the Histochemistry Society.

[22]  Samuel Matej,et al.  Efficient 3D grids for image reconstruction using spherically-symmetric volume elements , 1995 .

[23]  G A Perkins,et al.  Electron tomography of large, multicomponent biological structures. , 1997, Journal of structural biology.

[24]  J. Frank,et al.  Double-tilt electron tomography. , 1995, Ultramicroscopy.

[25]  Rajkumar Buyya,et al.  High Performance Cluster Computing , 1999 .

[26]  T. K. Narayan,et al.  Evaluation of task-oriented performance of several fully 3D PET reconstruction algorithms. , 1994, Physics in medicine and biology.

[27]  J M Carazo,et al.  Three-dimensional reconstruction from reduced sets of very noisy images acquired following a single-axis tilt schema: application of a new three-dimensional reconstruction algorithm and objective comparison with weighted backprojection. , 1997, Journal of structural biology.

[28]  José María Carazo,et al.  Fully Three-Dimensional Reconstruction in Electron Microscopy , 1999 .