User talk:Tnakane: Difference between revisions

From Relion
Jump to navigation Jump to search
No edit summary
(Replaced content with "This page is for drafting. <span style="color: red"><b>DO NOT TRUST</b> information on this page.</span> <div style="color: gray"> </div>")
Tag: Replaced
 
Line 2: Line 2:


<div style="color: gray">
<div style="color: gray">
= Orientations =
Orientations (<code>rlnAngleRot</code>, <code>rlnAngleTilt</code>, <code>rlnAnglePsi</code>) in a STAR file rotate the reference into observations (i.e. particle image), while translations (<code>rlnOriginXAngstrom</code> and <code>rlnOriginYAngstrom</code>) shifts observations into the reference projection. For developers, a good starting point for code reading is <code>ObservationModel::predictObservation()</code> in [https://github.com/3dem/relion/blob/ver3.1/src/jaz/obs_model.cpp 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 <code>((int)xdim/2, (int)(ydim/2))</code> (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 (<code>rlnOriginX</code> and <code>rlnOriginY</code>) but changed to ångström (<code>rlnOriginXAngstrom</code> and <code>rlnOriginYAngstrom</code>) in RELION 3.1.
== Particle coordinates ==
The unit of particle coordinates in a micrograph (<code>rlnCoordinateX</code> and <code>rlnCoordinateY</code>) 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 [http://emlab.rose2.brandeis.edu/ctf CTFFIND3], also see the publication by [http://dx.doi.org/10.1016/S1047-8477(03)00069-8 Mindell et al (2003)].
== Higher order aberrations ==
<code>rlnOddZernike</code> contains coefficients for asymmetric (antisymmetric) Zernike polynomials Z<sub>1</sub><sup>-1</sup>, Z<sub>1</sub><sup>1</sup>, Z<sub>3</sub><sup>-3</sup>, Z<sub>3</sub><sup>-1</sup>, Z<sub>3</sub><sup>1</sup>, Z<sub>3</sub><sup>3</sup>, etc in this order. <code>rlnEvenZernike</code> contains coefficients for symmetric Zernike polynomials Z<sub>0</sub><sup>0</sup>, Z<sub>2</sub><sup>-2</sup>, Z<sub>2</sub><sup>0</sup>, Z<sub>2</sub><sup>2</sup>, Z<sub>4</sub><sup>-4</sup>, Z<sub>4</sub><sup>-2</sup>, Z<sub>4</sub><sup>0</sup>, Z<sub>4</sub><sup>2</sup>, Z<sub>4</sub><sup>4</sup>, etc in this order. Thus, the 7-th item in the <code>rlnEvenZernike</code>, Z<sub>4</sub><sup>0</sup>, 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, Z<sub>3</sub><sup>-1</sup> = (3r<sup>3</sup> - 2r) sin θ = 3 (k<sub>x</sub><sup>2</sup> + k<sub>y</sub><sup>2</sup>) k<sub>y</sub> - 2 k<sub>y</sub>, where k<sub>x</sub> and k<sub>y</sub> 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.
<code>rlnMagMatrix_00</code> to <code>rlnMagMatrix_11</code> represent the matrix <code>M</code> in the section 2.4 of [https://www.biorxiv.org/content/10.1101/798066v2 our preprint]. The values become larger when the observed particle in the real space <i>looks larger</i> than the reference projection at the nominal pixel size. This also means that the true pixel size is <i>actually smaller</i> than the nominal pixel size.


</div>
</div>

Latest revision as of 13:55, 25 November 2019

This page is for drafting. DO NOT TRUST information on this page.