CryoSPARC Guide
  • About CryoSPARC
  • Current Version
  • Licensing
    • Non-commercial license agreement
  • Setup, Configuration and Management
    • CryoSPARC Architecture and System Requirements
    • CryoSPARC Installation Prerequisites
    • How to Download, Install and Configure
      • Obtaining A License ID
      • Downloading and Installing CryoSPARC
      • CryoSPARC Cluster Integration Script Examples
      • Accessing the CryoSPARC User Interface
    • Deploying CryoSPARC on AWS
      • Performance Benchmarks
    • Using CryoSPARC with Cluster Management Software
    • Software Updates and Patches
    • Management and Monitoring
      • Environment variables
      • (Optional) Hosting CryoSPARC Through a Reverse Proxy
      • cryosparcm reference
      • cryosparcm cli reference
      • cryosparcw reference
    • Software System Guides
      • Guide: Updating to CryoSPARC v4
      • Guide: Installation Testing with cryosparcm test
      • Guide: Verify CryoSPARC Installation with the Extensive Validation Job (v4.3+)
      • Guide: Verify CryoSPARC Installation with the Extensive Workflow (≤v4.2)
      • Guide: Performance Benchmarking (v4.3+)
      • Guide: Download Error Reports
      • Guide: Maintenance Mode and Configurable User Facing Messages
      • Guide: User Management
      • Guide: Multi-user Unix Permissions and Data Access Control
      • Guide: Lane Assignments and Restrictions
      • Guide: Queuing Directly to a GPU
      • Guide: Priority Job Queuing
      • Guide: Configuring Custom Variables for Cluster Job Submission Scripts
      • Guide: SSD Particle Caching in CryoSPARC
      • Guide: Data Management in CryoSPARC (v4.0+)
      • Guide: Data Cleanup (v4.3+)
      • Guide: Reduce Database Size (v4.3+)
      • Guide: Data Management in CryoSPARC (≤v3.3)
      • Guide: CryoSPARC Live Session Data Management
      • Guide: Manipulating .cs Files Created By CryoSPARC
      • Guide: Migrating your CryoSPARC Instance
      • Guide: EMDB-friendly XML file for FSC plots
    • Troubleshooting
  • Application Guide (v4.0+)
    • A Tour of the CryoSPARC Interface
    • Browsing the CryoSPARC Instance
    • Projects, Workspaces and Live Sessions
    • Jobs
    • Job Views: Cards, Tree, and Table
    • Creating and Running Jobs
    • Low Level Results Interface
    • Filters and Sorting
    • View Options
    • Tags
    • Flat vs Hierarchical Navigation
    • File Browser
    • Blueprints
    • Workflows
    • Inspecting Data
    • Managing Jobs
    • Interactive Jobs
    • Upload Local Files
    • Managing Data
    • Downloading and Exporting Data
    • Instance Management
    • Admin Panel
  • Cryo-EM Foundations
    • Image Formation
      • Contrast in Cryo-EM
      • Waves as Vectors
      • Aliasing
  • Expectation Maximization in Cryo-EM
  • Processing Data in cryoSPARC
    • Get Started with CryoSPARC: Introductory Tutorial (v4.0+)
    • Tutorial Videos
    • All Job Types in CryoSPARC
      • Import
        • Job: Import Movies
        • Job: Import Micrographs
        • Job: Import Particle Stack
        • Job: Import 3D Volumes
        • Job: Import Templates
        • Job: Import Result Group
        • Job: Import Beam Shift
      • Motion Correction
        • Job: Patch Motion Correction
        • Job: Full-Frame Motion Correction
        • Job: Local Motion Correction
        • Job: MotionCor2 (Wrapper) (BETA)
        • Job: Reference Based Motion Correction (BETA)
      • CTF Estimation
        • Job: Patch CTF Estimation
        • Job: Patch CTF Extraction
        • Job: CTFFIND4 (Wrapper)
        • Job: Gctf (Wrapper) (Legacy)
      • Exposure Curation
        • Job: Micrograph Denoiser (BETA)
        • Job: Micrograph Junk Detector (BETA)
        • Interactive Job: Manually Curate Exposures
      • Particle Picking
        • Interactive Job: Manual Picker
        • Job: Blob Picker
        • Job: Template Picker
        • Job: Filament Tracer
        • Job: Blob Picker Tuner
        • Interactive Job: Inspect Particle Picks
        • Job: Create Templates
      • Extraction
        • Job: Extract from Micrographs
        • Job: Downsample Particles
        • Job: Restack Particles
      • Deep Picking
        • Guideline for Supervised Particle Picking using Deep Learning Models
        • Deep Network Particle Picker
          • T20S Proteasome: Deep Particle Picking Tutorial
          • Job: Deep Picker Train and Job: Deep Picker Inference
        • Topaz (Bepler, et al)
          • T20S Proteasome: Topaz Particle Picking Tutorial
          • T20S Proteasome: Topaz Micrograph Denoising Tutorial
          • Job: Topaz Train and Job: Topaz Cross Validation
          • Job: Topaz Extract
          • Job: Topaz Denoise
      • Particle Curation
        • Job: 2D Classification
        • Interactive Job: Select 2D Classes
        • Job: Reference Based Auto Select 2D (BETA)
        • Job: Reconstruct 2D Classes
        • Job: Rebalance 2D Classes
        • Job: Class Probability Filter (Legacy)
        • Job: Rebalance Orientations
        • Job: Subset Particles by Statistic
      • 3D Reconstruction
        • Job: Ab-Initio Reconstruction
      • 3D Refinement
        • Job: Homogeneous Refinement
        • Job: Heterogeneous Refinement
        • Job: Non-Uniform Refinement
        • Job: Homogeneous Reconstruction Only
        • Job: Heterogeneous Reconstruction Only
        • Job: Homogeneous Refinement (Legacy)
        • Job: Non-uniform Refinement (Legacy)
      • CTF Refinement
        • Job: Global CTF Refinement
        • Job: Local CTF Refinement
        • Job: Exposure Group Utilities
      • Conformational Variability
        • Job: 3D Variability
        • Job: 3D Variability Display
        • Job: 3D Classification
        • Job: Regroup 3D Classes
        • Job: Reference Based Auto Select 3D (BETA)
        • Job: 3D Flexible Refinement (3DFlex) (BETA)
      • Postprocessing
        • Job: Sharpening Tools
        • Job: DeepEMhancer (Wrapper)
        • Job: Validation (FSC)
        • Job: Local Resolution Estimation
        • Job: Local Filtering
        • Job: ResLog Analysis
        • Job: ThreeDFSC (Wrapper) (Legacy)
      • Local Refinement
        • Job: Local Refinement
        • Job: Particle Subtraction
        • Job: Local Refinement (Legacy)
      • Helical Reconstruction
        • Helical symmetry in CryoSPARC
        • Job: Helical Refinement
        • Job: Symmetry search utility
        • Job: Average Power Spectra
      • Utilities
        • Job: Exposure Sets Tool
        • Job: Exposure Tools
        • Job: Generate Micrograph Thumbnails
        • Job: Cache Particles on SSD
        • Job: Check for Corrupt Particles
        • Job: Particle Sets Tool
        • Job: Reassign Particles to Micrographs
        • Job: Remove Duplicate Particles
        • Job: Symmetry Expansion
        • Job: Volume Tools
        • Job: Volume Alignment Tools
        • Job: Align 3D maps
        • Job: Split Volumes Group
        • Job: Orientation Diagnostics
      • Simulations
        • Job: Simulate Data (GPU)
        • Job: Simulate Data (Legacy)
    • CryoSPARC Tools
    • Data Processing Tutorials
      • Case study: End-to-end processing of a ligand-bound GPCR (EMPIAR-10853)
      • Case Study: DkTx-bound TRPV1 (EMPIAR-10059)
      • Case Study: Pseudosymmetry in TRPV5 and Calmodulin (EMPIAR-10256)
      • Case Study: End-to-end processing of an inactive GPCR (EMPIAR-10668)
      • Case Study: End-to-end processing of encapsulated ferritin (EMPIAR-10716)
      • Case Study: Exploratory data processing by Oliver Clarke
      • Tutorial: Tips for Membrane Protein Structures
      • Tutorial: Common CryoSPARC Plots
      • Tutorial: Negative Stain Data
      • Tutorial: Phase Plate Data
      • Tutorial: EER File Support
      • Tutorial: EPU AFIS Beam Shift Import
      • Tutorial: Patch Motion and Patch CTF
      • Tutorial: Float16 Support
      • Tutorial: Particle Picking Calibration
      • Tutorial: Blob Picker Tuner
      • Tutorial: Helical Processing using EMPIAR-10031 (MAVS)
      • Tutorial: Maximum Box Sizes for Refinement
      • Tutorial: CTF Refinement
      • Tutorial: Ewald Sphere Correction
      • Tutorial: Symmetry Relaxation
      • Tutorial: Orientation Diagnostics
      • Tutorial: BILD files in CryoSPARC v4.4+
      • Tutorial: Mask Creation
      • Case Study: Yeast U4/U6.U5 tri-snRNP
      • Tutorial: 3D Classification
      • Tutorial: 3D Variability Analysis (Part One)
      • Tutorial: 3D Variability Analysis (Part Two)
      • Tutorial: 3D Flexible Refinement
        • Installing 3DFlex Dependencies (v4.1–v4.3)
      • Tutorial: 3D Flex Mesh Preparation
    • Webinar Recordings
  • Real-time processing in cryoSPARC Live
    • About CryoSPARC Live
    • Prerequisites and Compute Resources Setup
    • How to Access cryoSPARC Live
    • UI Overview
    • New Live Session: Start to Finish Guide
    • CryoSPARC Live Tutorial Videos
    • Live Jobs and Session-Level Functions
    • Performance Metrics
    • Managing a CryoSPARC Live Session from the CLI
    • FAQs and Troubleshooting
  • Guides for v3
    • v3 User Interface Guide
      • Dashboard
      • Project and Workspace Management
      • Create and Build Jobs
      • Queue Job, Inspect Job and Other Job Actions
      • View and Download Results
      • Job Relationships
      • Resource Manager
      • User Management
    • Tutorial: Job Builder
    • Get Started with CryoSPARC: Introductory Tutorial (v3)
    • Tutorial: Manually Curate Exposures (v3)
  • Resources
    • Questions and Support
Powered by GitBook
On this page
  • Aliasing
  • Nyquist Frequency
  • CTF Aliasing
  1. Cryo-EM Foundations
  2. Image Formation

Aliasing

PreviousWaves as VectorsNextExpectation Maximization in Cryo-EM

Last updated 5 months ago

Summary: When a signal is sampled at some sampling rate, the maximum frequency that sample can represent is twice the sampling rate. In images, the sampling rate is the pixel size. Thus, the pixel size sets an upper bound on the best resolution achievable for a given set of images.

Aliasing

In the physical world, signals are continuous. The pressure waves which make up sound vary smoothly between high and low pressure; light waves which form images vary from light to dark, etc. However, computers must represent these signals with discrete samples. When we sample a signal, we take measurements of that signal at evenly spaced positions in space or time.

For a concrete example, consider the samples below:

By eye, it seems obvious that these samples come from a simple sine wave which oscillates once:

However, these samples are explained just as well by a wave which oscillates seven times:

or thirteen times:

Indeed, there are in fact an infinite family of waves which perfectly fit any finite set of samples. There is no way of knowing which of these infinite waves truly gave rise to the samples we observe, so by convention we select the lowest-frequency wave which fits the observed data. However, what happens when there really is high-frequency information in our images which our samples cannot capture?

Nyquist Frequency

Say you collect an image of some object using a camera sensor which is 12 Å wide and has 6 pixels. Your pixel size is therefore 2 Å. Put another way, you are sampling the incoming image with a sampling rate of 12A˚\frac{1}{2 \AA}2A˚1​ , or one sample for every two Å. If your object is a sine wave which oscillates with a frequency of 112A˚\frac{1}{12 \AA{}}12A˚1​, the image is unambiguous:

There are plenty of samples along the entire wave to accurate capture its shape. What happens when the wave’s frequency increases?

When you image an object with a wavelength of 6 Å, the result is still correct. The exact position of the object relative to the pixels has imposed some asymmetry in the image, but if the object shifted left or right the result would again be symmetric.

What about an object with even higher frequency?

When the object has a wavelength of 3 Å, the samples line up perfectly with samples from a 6 Å wave as well. Because we always take the lowest frequency wave that explain our data, we incorrectly interpret our image as resulting from a 6 Å wave. This incorrect result is called “aliasing” and is an important effect in SPA.

The frequency beyond which aliasing occurs is called the Nyquist frequency, and is half the sampling frequency.

fNyquist=12×fsampling=2×1pixel sizef_{\mathrm{Nyquist}} = \frac{1}{2} \times f_{\mathrm{sampling}} = 2 \times \frac{1}{\mathrm{pixel\ size}}fNyquist​=21​×fsampling​=2×pixel size1​

Put another way, any spatial features in the object which are smaller than twice the pixel size will be aliased in the image. This results in the commonly quoted maximum resolution for a particle stack of twice the image pixel size.

Frequencies faster than the Nyquist frequency alias to the frequency below Nyquist and the same distance away from it. This operation is commonly described as “folding over” the Nyquist frequency. Note also that in the 3 Å wavelength example above, the aliased 6 Å wave has opposite phase. This is another property of aliasing: the “folding” operation also flips the phase of the aliased wave.

The phase flipping here is a property of all discrete representations of continuous signals. It is distinct from the specific mechanism of phase inversion that occurs in the electron microscope, which is modeled by the CTF.

In the animation below, the moment the true frequency crosses the Nyquist frequency the aliased wave begins to decrease in frequency, rather than increase, since it is folding over the Nyquist frequency.

CTF Aliasing

Just as particle images are represented with pixels in real space, their Fourier transform is represented with pixels in Fourier space. The box in Fourier space is the same size as the box in real space, but each pixel represents a certain range of frequencies rather than a certain region of space. Specifically,

1 Fourier pixel=1N×pixel size\mathrm{1\ Fourier\ pixel} = \frac{1}{N \times \mathrm{pixel\ size}}1 Fourier pixel=N×pixel size1​

where N is the box size in pixels and the pixel size is the size of a pixel in the real space image. For instance, if a particle is represented with a box of 128 pixels which each represent 2 Å, the Fourier space particle image has a Nyquist frequency of 12×128A˚−1=1256A˚−1\frac{1}{2 \times 128} \mathrm{\AA{}^{-1}} = \frac{1}{256} \mathrm{\AA{}^{-1}}2×1281​A˚−1=2561​A˚−1. If the CTF oscillates with a frequency greater than this, the contrast transfer function will be aliased.

Consider the true contrast transfer function of an image collected in a typical electron microscope used in single-particle analysis (300 kV, 2.7 Cs) with 2 µm defocus.

However, we cannot capture the full, continuous contrast transfer function. We must model it from our images, which are sampled using discrete Fourier pixels. If the Fourier pixels are small, the resulting image of the contrast transfer function may look like this:

At the higher frequencies there is some disruption of the amplitude, but the representation is mostly accurate. However, if the Fourier pixels are too large, the contrast transfer function will be significantly aliased.

In this rather extreme example, the contrast transfer function modeled from the image (black) is significantly different from the true contrast transfer function (grey), especially past approximately 4 Å. This happens because in this region the contrast transfer function starts oscillating more quickly than the Nyquist limit of the Fourier space image. The high-frequency contrast transfer function oscillations are therefore aliased to incorrect lower frequencies.

This effect is more pronounced with:

  • higher defocus values, which increase the true oscillation rate of the CTF,

  • smaller pixel sizes, which increase the real space Nyquist limit, requiring the Fourier image to represent a larger region of the CTF, and

  • smaller box sizes, which reduce the sampling rate in Fourier space.

Since defocus cannot be changed after acquisition, if significant CTF aliasing is observed, the remaining means of removing it are therefore

  • downsampling, which “zooms in” on a subregion of the CTF by removing higher frequencies, or

  • using a larger box which gives more pixels in Fourier space to represent the CTF.

Of course, the problem of contrast transfer function aliasing is often rendered moot by the poor signal-to-noise ratio of cryo-EM data. For instance, even the severely-aliased example above more-or-less correctly models data up to the resolution to which the CTF fits the data well.

Seven dots, clearly arranged along a single cycle of a sine wave
The seven points are connected by a sine wave which oscillates once
The same seven points are now connected with a sine wave which oscillates seven times
The same seven points are connected by a wave which oscillates thirteen times
The same set of points are shown twice. At top, they are connected by a sine wave which oscillates once. On the bottom, they are connected by steps at the height of each sample, and a dotted wave which oscillates once.
The wave at the top now oscillates twice, and the steps at the bottom also oscillate twice. The dotted wave is correct, also oscillating twice.
The wave at the top now oscillates four times. At the bottom, the steps only oscillate twice, so the dotted wave only oscillates twice as well. The phase of the steps and dotted wave are also flipped.
An arrow pointing up at a frequency of 0.6 times the sampling frequency is "folde" across the Nyquist frequency to become an arrow pointing down (representing a phase shift) at 0.4 times the sampling frequency
A graph of the contrast transfer function
The sane CTF as above, but discretely sampled. A the highest frequencies the amplitude is slightly off, but the overall signal is correct.
The same CTF as above, but now the sampling rate is too low to accurately capture the CTF.
Four graphs of the CTF are shown. At the top is the continuous CTF. Below that, the aliased CTF shown previously. Third from the top, the aliased CTF is shown for a downsampled image. This produces a CTF which only extends to around 4 Å, where aliasing is not significant. Finally, the CTF is shown from a bigger box. The larger box means the CTF has more samples, reducing aliasing.