Hi, community!
I am implementing a Python script to quantify the agreement between an experimental Cryo-EM particle image and a re-projection of the corresponding 3-D reconstruction.
The pipeline is as follows:
-
Reconstruct a consensus 3-D volume
Min cryoSPARC. -
Compute the 3-D Fourier transform of
M(without normalization) to obtainMʹ. -
Center the zero-frequency component of
Mʹ. -
Extract a central 2-D slice
SfromMʹat the posePreported for the particle. -
Transform the experimental particle image
I(corresponding to poseP) into reciprocal space and apply CTF correction to yieldIʹ -
Center the zero-frequency component of
Iʹ. -
Evaluate the correlation between
SandIʹ.
Because S and Iʹ should share the same pose P, their correlation is expected to be high. In practice, however, the value is essentially zero, indicating no detectable correlation.
Are there any steps I might be missing in my current pipeline?