Reconstruction

5. Reconstruction

Two requirements are needed for 3D structure reconstruction from projections. One is projection (or in our case, particle averages) and the other is projection angle for each projection (see Fig. 1). With these two inputs, we can reconstruct the 3D structures using techniques as seen in medical image reconstruction, which will be briefly described in the following.

Figure 1

Reconstruction

Pipeline of reconstruction, given particle averages and their corresponding orientations.

Particle averages have already been discussed in Classification. The determination of projection orientation can be accomplished using common line method or reference-based method. The first method is to determine the relative orientation between two projections by detecting their common line. The second method assumes an initial 3D model available and each of the particle is cross-correlated with each of the projections of the 3D model. For each particle image, the projection corresponding to the highest cross-correlation coefficient is assumed to have the same orientation as that particle image. Clearly the second method can be made iterative to improve the results and can also be used to solve the particle classification. However, the initial 3D model must be carefully chosen to guarantee the convergence of this iterative process.

CT/MRI Reconstruction

The CT/MRI reconstruction is based on a well-known theorem called the Fourier Slice Theorem, which says that the Fourier transform of any projection of the original structure must be equal to a certain zero-crossing slice of the Fourier transform of the original structure. Fig.2(a) shows an example of a 2D structure reconstruction from its 1D projections. With enough number of projections, we can compute their Fourier transforms, each of which is “embedded” in the Fourier space and then an interpolation technique is used to “reconstruct” the Fourier transform of the original structure f(x, y) in the entire Fourier space (see fig.2 (c)). After we obtain F(u, v), the inverse Fourier transform of F(u, v) gives the original structure f(x,y). Fig.2 (b) illustrates the overall pipeline of reconstruction. It is straightforward to extend this procedure into 3D structure reconstruction.

Figure 2

Reconstruction 1 Reconstruction 2 Reconstruction 3
(a) Different projections (b) Fourier slices and interpolation (c) Illustration of reconstruction

Results

In this section we show one example of reconstruction on P97 dataset (from Dr. Bridget Carragher at The Scripps Research Institute). We used EMAN tool (see [1,2,3]) developed at Baylor College of Medicine, Houston to reconstruct the 3D structures on the above mentioned datasets. However, instead of using the manual or semi-automatic techniques seen in EMNA for picking particles, we used our automatic particle picking algorithm and took the detected particles as input to EMAN.

Figure 3

EMAN 1 EMAN 2 EMAN 3
Reconstruction results on P97 dataset using EMAN based on the particles boxed out by our algorithm.