Differences between revisions 29 and 30
Revision 29 as of 2013-07-23 19:23:41
Size: 2599
Editor: penczek
Comment:
Revision 30 as of 2017-05-07 11:51:17
Size: 2636
Editor: penczek
Comment:
Deletions are marked like this. Additions are marked like this.
Line 17: Line 17:
 * To apply CTF object to an image, use '''[[filt_ctf|filt_ctf(img, ctf)]]''', where ctf is a CTF object.  * To apply CTF object to an image, use '''[[filt_ctf|filt_ctf(img, ctf)]]''', where ctf is a CTF object created using '''[[generate_ctf]]'''.

How to use CTF in SPARX

In SPARX, the EMAN CTF object are used to store and apply all CTF parameters.

  • To generate a CTF object, see function generate_ctf(p), here p is a list of parameters [defocus, cs, voltage, apix, bfactor, amp_contrast, astigmatism amplitude, astigmatism anglet]. Their definitions are:

    defocus

    defocus in `µ`m. Positive number means underfocus.

    cs
    Spherical aberration constant in mm. (Default 2.0)
    voltage

    microscope voltage in `kV`. (Default 300 kV)

    apix

    pixel size in `Å`.

    amp_contrast
    Amplitude contrast ratio, unitless, given as percentage of (0,100). (Default 10%)
    bfactor

    B-factor defining exponential decay of Fourier amplitudes decay `Å^2`. (Default 0.0, i.e., no decay).

    astigmatism amplitude
    in [nm]. (Default 0.0)
    astigmatism angle
    in [degrees]. (Default 0.0)

    Notice: In some other software packages, the defocus value is in Å and amplitude contrast ratio is in the range of (0, 1). This function will automatically determine whether these two CTF parameters are in this convention and convert them into our convention.

  • To apply CTF object to an image, use filt_ctf(img, ctf), where ctf is a CTF object created using generate_ctf.

  • To put a CTF object into the header of an image, it can be done in two ways:
    1. use function set_ctf(img, p), where p is a list of parameters [defocus, cs, voltage, apix, bfactor, ampcont].

    2. if the CTF object already exists, one can use img.set_attr("ctf", ctf) directly.

  • To recover the numerical values of CTF parameters from CTF object stored in the header of an image, use function p=get_ctf(img).

  • To access or set individual parameters in a ctf object, use:
  • bfactor = ctf.bfactor
  • ctf.bfactor = 200.0

Description

  • The Contrast Transfer Function of electron microscope has the following form:
  • `gamma(K,"defocus")=2pi(-1/2"defocus"*lambda*K^2+1/4*"Cs"*lambda^3*K^4)`

  • `CTF(K)=(sqrt(1-A^2)*sin(gamma(K,"defocus")) - A*cos(gamma(K,"defocus"))*e^(-"B_factor"/4*K^2)`

  • where `K=sqrt(K_x^2+K_y^2+K_z^2)` and `A` is amplitude contrast.

Reference

1. J. Frank Optik (1973) 38:519

2. R. Wade & J. Frank Optik (1974) 49:81.

3. Z. Huang, P. R. Baldwin, S.Mullapudi, and P .A. Penczek, Automated determination of parameters describing power spectra of micrograph images in electron microscopy. J. Struct. Biol. 144 (2003), pp. 79-94.

CTF_info (last edited 2017-05-07 11:51:17 by penczek)