Wednesday, July 20, 2022

Model bias in cryo-EM - Einstein from noise

 
 
 
In the article of Maxim Shatsky et al.  entitled "A Method for the Alignment of Heterogeneous Macromolecules from Electron Microscopy" (J Struct Biol. 2009, 166(1): 67–78. doi:10.1016/j.jsb.2008.12.008.), the authors show that pure noise images aligned to a reference image give an average image resembling the reference. This is a perfect example of the model bias occurring during the 2D classification step. Let's do it!!!

1- Pre-processing

In the article, they chose the iconic image of Einstein by Arthur Sasse. I just crop and resize this image in 128x128. Additionally, I convert it in 32-bits (Image > Type > 32-bits) and rescale the image in order to have a zero mean and a sigma equal to 1.0 (Process > Math > Macro with a formula v=(v-mean)/stdev ). To get the values of mean and stdev, the simplest way is to calculate the histogram (Analyze > Histogram) (Fig.1).

Fig.1: Reference image inspired by those used in the article of Shatsky et al. This image is copyrighted and I did not granted permission, so this image will be removed upon request.

2. 2D alignment by cross-correlation

3. The source code

The script below is composed of three parts:
  • (i) Initialization 
  • (ii) Creation of the stack of gaussian noise slices
  • (iii) 2D alignment of noise to the reference by cross-correlation
 
The computation of the average image from the "aligned" stack is not included in the script. You can easily do that in (Image > Stacks > Z-Project and choose Average Intensity).
 
+++ IJ JavaScript snippet +++ +++End of JavaScript snippet +++
 

4. Results 

 

Fig. ???: Correlations maps

 
Fig. ???: Aligned noisy images


Fig. ???: Average Images calculated from 25,100, 500, and 1000 noisy images, respectively.

 

 Thank you for reading.

5. Other crazybiocomputing posts

Further readings are available in ...
  • CryoEM Series  [Link]
  • Image Processing TOC [Link]

 

No comments:

Post a Comment