High Performance Computing in Electron Microscope Tomography of Complex Biological Structures

Series expansion reconstruction methods using smooth basis functions are evaluated in the framework of electron tomography of complex biological specimens, with a special emphasis upon the computational perspective. Efficient iterative algorithms that are characterized by a fast convergence rate have been used to tackle the image reconstruction problem. The use of smooth basis functions provides the reconstruction algorithms with an implicit regularization mechanism, very appropriate for noisy conditions. High Performance Computing (HPC) techniques have been applied so as to face the computational requirements demanded by the reconstruction of large volumes. An efficient domain decomposition scheme has been devised that leads to a parallel approach with capabilities of interprocessor communication latency hiding. Comparisons with Weighted BackProjection (WBP), the standard method in the field, are presented in terms of computational demands as well as in reconstruction quality under highly noisy conditions. The combination of efficient iterative algorithms and HPC techniques have proved to be well suited for the reconstruction of large biological specimens in electron tomography, yielding solutions in reasonable computational times.

[1]  Yair Censor,et al.  BICAV: An Inherently Parallel Algorithm for Sparse Systems with Pixel-Dependent Weighting , 2001, IEEE Trans. Medical Imaging.

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

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

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

[5]  Robert M. Lewitt,et al.  A comparison of transform and iterative reconstruction techniques for a volume-imaging PET scanner with a large axial acceptance angle , 1995 .

[6]  Barry Wilkinson,et al.  Parallel programming , 1998 .

[7]  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.

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

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

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

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

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