relion_tomo_align

This program refines the projections that map 3D space onto the images of the tilt series, as well as (optionally) also the beam-induced motion trajectories of the particles. Each projection is optimised using the full 5 degrees of freedom: the assumption of a common tilt axis is abandoned. Even if no beam-induced motion is estimated, the (in that case static) 3D particle-positions are also optimised by the program. This is because those 3D positions cannot be assumed to be known in advance, since they have only been inferred from the observed 2D particle-positions in the individual tilt images. Therefore, this program always looks for the optimal set of 3D positions and projections.

If the particle motion is also being estimated (using the --motion argument), then the 3D positions of the particles are allowed to differ between tilt images, albeit only in a limited and spatially smooth manner. This is done analogously to Bayesian polishing in 2D (implemented in relion_motion_refine), except that the trajectories are estimated in 3D, and that no acceleration constraint is imposed. Because part of the particle motion is rigid, the motion estimation also always includes a rigid alignment of the frames.

Please note that estimating beam-induced motion significantly increases the runtime of the program, so it should only be done in the final rounds of optimisation.

The estimated particle trajectories are written into the output trajectory set, the projection matrices (containing the rigid part of the motion as well as the alignment of the tilt series - these cannot be distinguished) into the output tomogram set and the estimated (initial) particle positions into the output particle set which replace the given trajectory, tomogram and particle sets in the output optimisation set.

General program arguments:

  • --i and/or --p, --t, --mot, --ref<1/2>, --mask and --fsc: input optimisation-set and/or its components (see optimisation set).

  • --b: box size to be used for the estimation. Note that this can be larger than the box size of the reference map. A sufficiently large box size allows more of the high-frequency signal to be captured that has been delocalised by the CTF.

  • --r: maximal assumed error in the initial 2D particle-positions (distances between the projected 3D positions and their true positions in the images), given in pixels.

  • --j: number of threads to be used. This should be set to the number of CPU cores available.

  • --o: name of the output directory (that will be created).

Program output:

After running the program, the output directory (--o argument) will contain the following items:

  • motion.star: a new trajectory set

  • tomograms.star: a new tomogram set containing the newly estimated projection matrices (rlnTomoProj<X/Y/Z/W>)

  • particles.star: a new particle set containing the newly estimated particle positions in 3D (rlnCoordinate<X/Y/Z>)

  • optimisation_set.star: a optimisation set pointing to the new trajectory, tomogram and particle sets.

  • Trajectories/< tomogram_name >_x<1/8>.ply: meshes of 3D trajectories in ply format, scaled by a factor of 1 or 8, respectively. Note that only the trajectories themselves are scaled - their positions are those of the particles. They can be viewed in relation to (low-magnification) tomograms as produced by relion_tomo_reconstruct_tomogram in a viewer that supports both meshes and 3D images (such as e.g. Paraview).

  • note.txt: a text file containing the command issued to run the program.