The CTF Model(s) used in EMAN2
UNDER CONSTRUCTION
Like other aspects of EMAN2, the CTF model infrastructure is extensible, though less so than the Processors, Reconstructors, etc. At present, there are 2 types of CTF instances in EMAN2: EMAN1Ctf and EMAN2Ctf. The EMAN1Ctf is provided for backwards compatibility, so if particles processed in EMAN1 are read into EMAN2, they will have a valid CTF definition in their header.
The new EMAN2 CTF model is defined as follows:
- float defocus; // defocus in microns, positive underfocus float dfdiff; // defocus difference for astigmatism, defocus is the major elliptical axis float dfang; // angle of the major elliptical axis in degrees measured counterclockwise from x
float bfactor; // B-factor in 1/A2 expressed using the x-ray convention (e-B/4 s2 in amplitude space) EMAN1 used E-B s^2 float ampcont; // amplitude contrast as a percentage ie- this should be 10 for 10% amp contrast float voltage; // microscope voltage in kV float cs; // Cs in mm float apix; // A/pix value used when generating 2D results float dsbg; // ds value for background and SNR vector<float> background; // background intensity, 1 value per radial pixel (NX/2, corners ignored) vector<float> snr; // SNR, 1 value per radial pixel (NX/2, corners assumed 0)
Astigmatism is, at present, disabled/unused, but it is on my list of things to finalize just AFTER the workshop. You can certainly store the values for it there, but the generated curves won't contain it.
The parameters can be initialized as:
ctf=EMAN2Ctf() ctf.from_dict("defocus":2.5,"voltage":300.0", ...)
- or -
ctf.defocus=2.5 ctf.voltage=300.0
In an EMData() object:
img=EMData(256,256) ctf=EMAN2Ctf() ctf.from_dict("defocus":2.5,"voltage":300.0", ...) img.set_attr("ctf",ctf)
This image attribute can be archived in HDF and BDB format, but will be lost in other formats.
While these parameters can obviously be extracted from the CTF object itself, this is not the normal way of using these objects. Each Ctf object (other subclasses could be added) has two functions:
vector < float >compute_1d(int size,float ds, CtfType type, XYData * struct_factor = 0); void compute_2d_complex(EMData * image, CtfType type, XYData * struct_factor = 0);
These functions return computed CTF-derived curves in 1-D or 2-D. In 1-D, specify the size of the returned array, and ds (the change in spatial frequency for each element, starting at 0 for the first). For 2-D you pass in a complex image of the desired dimensions, and its contents are replaced by the requested curve (in this case, the internal apix value is used to determine an appropriate ds). The struct_fact is not required for any current operations. CtfType is one of:
Ctf.CtfType. -->
- CTF_AMP, // ctf ampltidue only CTF_SIGN, // ctf sign (+-1) CTF_BACKGROUND, // Background, no ctf oscillation CTF_SNR, // Signal to noise ratio CTF_WIENER_FILTER, // Weiner Filter = 1/(1+1/snr) CTF_TOTAL // AMP*AMP+NOISE
For example, to pad (oversample) an image, and perform phase flipping, one could:
- ctf=im1["ctf"] im1.clip_inplace(Region(-(ys2*(oversamp-1)/2),-(ys2*(oversamp-1)/2),ys,ys)) # pad in real-space fft1=im1.do_fft() flipim=fft1.copy()
ctf.compute_2d_complex(flipim,Ctf.CtfType.CTF_SIGN) # get a phase flipping map (this could/should be cached for multiple images) fft1.mult(flipim) # flip phases out=fft1.do_ift() out["ctf"]=ctf out.clip_inplace(Region(int(ys2*(oversamp-1)/2.0),int(ys2*(oversamp-1)/2.0),ys2,ys2)) # unpad out.write_image(phaseflip,i)