Job: Local Motion Correction

Local motion correction.

At a Glance

Track and correct the motion of individual particles in movies.

Description

Local Motion Correction job corrects beam-induced motion on a per-particle basis, taking full-frame-motion corrected movies and particle pick locations as inputs. A square patch around each particle is extracted from each frame of the movie. The trajectory of these patches is modeled and then used to extract a motion-corrected patch for each particle and each frame. Local motion correction also applies Dose Weighting. Local motion correction is based on alignparts_lbfgs (Rubinstein et al. 2015).

Two versions of Local Motion Correction are available - Local Motion Correction runs on a single GPU, while Local Motion Correction (Multi-GPU) uses multiple available GPUs to parallelize processing.

This job does not require a high-quality 3D reference volume (unlike Reference Based Motion Correction), but it does require well-centered particle picks. A typical workflow would be to first use patch motion correction on the raw movies, perform picking, and then re-perform motion correction with Local motion correction.

See Motion Correction for more details about the different types of motion correction available in CryoSPARC.

Inputs

Movies

Local motion correction requires raw movies in .mrc , .tif , .eer, or .mrc.bz2 format, typically from an Import Movies job.

Particles

Particles should be well-centered and relatively clean of junk or empty picks, most often from a 2D Classification or any 3D refinement job.

Commonly Adjusted Parameters

Only process this many movies

This parameter selects the first n movies to process, rather than the entire input stack. It is most useful when working with a subset of movies to assess data quality before committing to a full processing pipeline.

If collection parameters changed during the collection, it may be better to use Exposure Sets to select a random subset instead.

Save results in 16-bit floating point

Saves the output micrographs in 16-bit floating point. We recommend that this option is turned on to save space at minimal loss of accuracy. See the 16-bit floating point article for more information.

Override e/A^2/frame

When Apply exposure weighting is enabled, this parameter can be used to override the exposure level of the movie data, in case the level was not correctly set at import time. In contrast to the Import Movies job's Total exposure dose (e/A^2) parameter, Local Motion Correction's parameter must be specified per-frame. Exposure weighting is based on the reference curves from Grant and Grigorieff.

Extraction box size

The side length of the extraction box, in pixels. The general rule of thumb is twice the width of the particle, but more precise guidance is given by Rosenthal and Henderson as

Box Size (A˚)=d + (2λ×δr)\mathrm{Box\ Size\ (\AA{})} = d\ +\ (2\lambda{}\times{}\frac{\delta{}}{r})

where

  • d is the particle diameter in Å,

  • λ is the electron wavelength (approximately 0.020 Å and 0.025 Å for 300 kV and 200 kV microscopes, respectively),

  • δ is the defocus (in Å, not μm), and

  • r is the expected resolution of the final reconstruction (the average CTF fit would be a good starting value if data quality is not known)

The first term of this equation is simply the particle diameter, an obvious lower limit on box size, since the entire particle must fit in the box. The second term models the displacement of signal by the CTF. With higher defocus, higher frequency information is displaced by a greater distance, hence the direct dependence on defocus and inverse dependence on final resolution, since higher resolution decreases r.

Fourier crop to box size (pix)

Extracted, motion-corrected particle images will be Fourier cropped to this box size.

Outputs

Particles extracted

Since particle locations are already known and motion of the entire movie is not modeled, this job outputs motion corrected particle images rather than whole micrographs.

Diagnostic plots

Local Motion Correction also produces a plot of the modeled movement of each particle:

The local motion plot shows the motion modeled for each particle patch. The trajectories are scaled by 40x to be visible on the micrograph scale.

The grey line shows the raw trajectory. Because single-particle data typically has a very poor signal-to-noise ratio, these trajectories are very noisy and likely not a good representation of the true movement of the particle. The per-frame trajectory is therefore smoothed out to produce the final result. This smooth trajectory is in red.

Common Next Steps

Particle images from Local Motion Correction are suitable for use in a variety of jobs, including 3D Refinements and heterogeneity analysis techniques.

If a high-quality reconstruction is available for the dataset, Reference Based Motion Correction will likely perform better than Local Motion Correction. Like Local Motion Correction, it directly models particle trajectories. However, Reference Based Motion Correction also calculates empirical dose weights (rather than use the reference curves), which can lead to improved final reconstruction quality.

References

  1. Rubinstein JL, Brubaker MA. Alignment of cryo-EM movies of individual particles by optimization of image translations. J Struct Biol (2015).

  2. Rosenthal, P. B. & Henderson, R. Optimal Determination of Particle Orientation, Absolute Hand, and Contrast Loss in Single-particle Electron Cryomicroscopy. Journal of Molecular Biology 333, 721–745 (2003).

  3. A table of electron wavelengths for a given accelerating voltage can be found here: https://www.jeol.com/words/emterms/20121023.071258.php#gsc.tab=0

  4. Grant, T. & Grigorieff, N. Measuring the optimal exposure for single particle cryo-EM using a 2.6 Å reconstruction of rotavirus VP6. eLife 4, e06980 (2015).

Last updated