Lee Junghoon, Zheng Yili, Doerschuk Peter C
Sch. of Electr. & Comput. Eng., Purdue Univ., West Lafayette, IN 47907, USA.
Conf Proc IEEE Eng Med Biol Soc. 2006;2006:2538-41. doi: 10.1109/IEMBS.2006.260210.
In a cryo electron microscopy experiment, the data is noisy 2-D projection images of the 3-D electron scattering intensity where the orientation of the projections is not known. In previous work we have developed a solution for this problem based on a maximum likelihood estimator that is computed by an expectation maximization algorithm. In the expectation maximization algorithm the expensive step is the expectation which requires numerical evaluation of 3- or 5-dimensional integrations of a square matrix of dimension equal to the number of Fourier series coefficients used to describe the 3-D reconstruction. By taking advantage of the rotational properties of spherical harmonics, we can reduce the integrations of a matrix to integrations of a scalar. The key property is that a rotated spherical harmonic can be expressed as a linear combination of the other harmonics of the same order and the weights in the linear combination factor so that each of the three factors is a function of only one of the Euler angles describing the orientation of the projection. Numerical example of the reconstructions is provided based on Nudaurelia Omega Capensis virus.