User talk:Tnakane
This page is for drafting. DO NOT TRUST information on this page.
Orientations
Orientations (rlnAngleRot
, rlnAngleTilt
, rlnAnglePsi
) in a STAR file rotate the reference into observations (i.e. particle image), while translations (rlnOriginXAngstrom
and rlnOriginYAngstrom
) shifts observations into the reference projection. For developers, a good starting point for code reading is ObservationModel::predictObservation()
in src/jaz/obs_model.cpp.
Coordinate system
In compliance with the Heymann, Chagoyen and Belnap (2005) standard RELION uses a right-handed coordinate system with orthogonal axes X, Y and Z. Right-handed rotations are called positive.
Image center
The center of rotation of a 2D image of dimensions xdim x ydim is defined by ((int)xdim/2, (int)(ydim/2))
(with the first pixel in the upper left being (0,0). Note that for both xdim=ydim=65 and for xdim=ydim=64, the center will be at (32,32). This is the same convention as used in SPIDER and XMIPP. Origin offsets reported for individual images translate the image to its center and are to be applied BEFORE rotations.
The unit of particle translations was pixel (rlnOriginX
and rlnOriginY
) but changed to ångström (rlnOriginXAngstrom
and rlnOriginYAngstrom
) in RELION 3.1.
Particle coordinates
The unit of particle coordinates in a micrograph (rlnCoordinateX
and rlnCoordinateY
) is pixel in the aligned and summed micrograph (possibly binned from super-resolution movies). The origin is the first element in the 2D array of an MRC file. The origin is displayed at the upper-left corner in RELION (other programs might display in other ways).
For compatibility with other EM programs (e.g. UCSF MotionCor2, IMOD), TIFF images are flipped along the slow axis when being read into memory or written to a file. This happens regardless of the TIFFTAG_ORIENTATION value in the header.
Euler angle definitions
Euler angle definitions are according to the Heymann, Chagoyen and Belnap (2005) standard:
* The first rotation is denoted by phi or rot and is around the Z-axis. * The second rotation is called theta or tilt and is around the new Y-axis. * The third rotation is denoted by psi and is around the new Z axis
As such, RELION uses the same Euler angles as XMIPP, SPIDER and FREALIGN.
Contrast Transfer Function
CTF parameters are defined as in CTFFIND3, also see the publication by Mindell et al (2003).
Higher order aberrations
rlnOddZernike
contains coefficients for asymmetric (antisymmetric) Zernike polynomials Z1-1, Z11, Z3-3, Z3-1, Z31, Z33, etc in this order. rlnEvenZernike
contains coefficients for symmetric Zernike polynomials Z00, Z2-2, Z20, Z22, Z4-4, Z4-2, Z40, Z42, Z44, etc in this order. Thus, the 7-th item in the rlnEvenZernike
, Z40, is related to an error in the spherical aberration coefficient.
Look at the table in Wikipedia https://en.wikipedia.org/wiki/Zernike_polynomials#Zernike_polynomials but ignore square root terms, as the coefficients are not normalised in RELION. For example, Z3-1 = (3r3 - 2r) sin θ = 3 (kx2 + ky2) ky - 2 ky, where kx and ky are wave-numbers in the reciprocal space (1 / Å).
Anisotropic magnification corrections
Transformation by anisotropic magnification brings the reference into observations (i.e. particle images) in real space. Note that stretching in real space is shrinking in reciprocal space and vice versa.
rlnMagMatrix_00
to rlnMagMatrix_11
represent the matrix M
in the section 2.4 of our preprint. The values become larger when the observed particle in the real space looks larger than the reference projection at the nominal pixel size. This also means that the true pixel size is actually smaller than the nominal pixel size.