3dem / relion

Image-processing software for cryo-electron microscopy
https://relion.readthedocs.io/en/latest/
GNU General Public License v2.0
444 stars 197 forks source link

Reconstructing with small number of frames to assess damage #945

Open cgisriel opened 1 year ago

cgisriel commented 1 year ago

Hello, I have an ~2 angstrom resolution structure processed using Relion 3.1. We would like to create the reconstruction with 2 frames, 3, 4, etc., so that we can assess radiation damage to the metallocofactor. What is the best way to go about this? I have previously started over entirely, selecting only the first 3 frames in motion correction and proceeded to perform refinement using the coordinates from the previous full dose workflow. This required subsequent CTF refinement and polishing. But, I would prefer to simply reconstruct the "full dose" map using only the first few frames.

One paper states that after they produced the full dose map, they performed an additional polishing step, setting the frames used to only the first 2 or 3, then used relion_reconstruct to reconstruct without refinement of positions and orientations. I assume they somehow selected which frames to use in the reconstruction using relion_reconstruct? I don't know how to use relion_reconstruct for this purpose.

I would be extremely grateful for any assistance! Thanks

biochem-fan commented 1 year ago

they performed an additional polishing step, setting the frames used to only the first 2 or 3, then used relion_reconstruct to reconstruct without refinement of positions and orientations. I assume they somehow selected which frames to use in the reconstruction using relion_reconstruct?

After Polish using the selected frames, the outputs (shiny particles) obviously contain contributions from the selected frames only. So you don't need to select frames in relion_reconstruct.

You can also use relion_movie_reconstruct to make per-frame reconstructions. See the following discussions in CCPEM.

cgisriel commented 1 year ago

Ah, so it is as simple as selecting only the number of frames I would like to reconstruct with in polishing. Thanks very much and I will give it a try with relion_movie_reconstruct, too. I will report back.