Three-dimensional reconstruction of helical structures with fast inversion of very large Fourier transforms

A single projection of a helical distribution of matter allows one to obtain the complete three-dimensional reconstruction of the structure. This task is usually performed by a Fourier-Bessel algorithm, which is more efficient than a customary fast Fourier transform inversion. This article describes how to achieve such a result by a direct Fourier method in a reasonable time. Once the two-dimensional transform of the projection is obtained from the source image, it is possible to build up the three-dimensional transform array, in Cartesian coordinates, that yields the reconstruction by a straightforward Fourier inversion. Images of projected helices should be studied with high sampling rates to enhance the resolution, and the segments of helix should be long enough to give a satisfactory signal-to-noise ratio. These conditions result in three-dimensional transform arrays that would require one or more gigabytes of storage. The strategy proposed here requires much less storage and is fast enough to allow the reconstruction to be performed with different parameters and filters in a very short time without any sacrifice in resolution.