PROJECTIONS with EULER ANGLES 

The problem addressed here is given any set of Euler angles to know what is the projection that Xmipp (and Spider) will associate to them. The Euler angles are defined as follow:
 
 
rotational or azimuthal: First rotation.  Around Z axis
tilting: Second rotation.  Around new Y axis 
psi:  Third rotation.  Around new Z axis

All angles are measured using the left hand thumb rule (when aligning the left hand thumb with the rotation axis, positive angles are measured according to the direction of the rest of the fingers).

There is an approach which consists of rotating the volume according to the Euler angles and then projecting from Z-axis. An alternative one, and the one I have followed and found more practical is to move the view point (or better said the projection plane) according to the Euler Angles. Both ways will be explained here.

First, the way Xmipp define the coordinate axes is like this.


MOVING THE PROJECTION PLANE



The idea is to define a projection plane and move it around the volume until we get the desired point of view. It would seem that we had to make the rotations in the same order as they are defined in the Euler angles but it is not totally true. We can interchange tilt and rotation to make an easiest method to get manually the projection. I will describe a chain of steps in order to get the projection.

Getting Started

1)Construct a projection plane: take a little piece of paper and mark on it the image coordinate system will define the orientation of the projection. We will look the projection in the orientation just as it is shown.

2) Our point of view is defined just now for a little point on the bottom of the sphere, while the projection plane is just on the opposite side.
 
 






We begin now to read eulerian angles
3) TILT:
move the point of view as many degrees as necessary towards right (positive tilting angles) or left (negative tilting angles). Projection plane must move as if the system point_of_view-plane_of_projection were a whole.For example, 90º of tilting would be






4) ROTATIONAL (AZIMUTHAL) ANGLE:
move the point of view towards you (positive azimuthal angles) or away from you (negative azimuthal  angles) (Notice that this directions are simplifications for when tilting angles are positive, you must keep in mind the direction of turn and not the description of towards or away from you. If tilt angle is negative this directions are reversed in respect to you but not in the direction of turn).  For example from the above situation of tilt=90, now let's take rotational=90. The new situation would be
 
 

5) PSI:
for taking psi into account we must turn the projection plane together with the point of view. Looking from the point of view positive psi angles turns counterclockwise. For example, psi=90, would give,






PROJECTING

Now we have to project the volume on the projection plane. As an example I have defined a volume which has got a sphere along y-axis, a larger ellipsoid along z-axis and a cylinder along x-axis as you can see in the next figure. If we project for (tilt=90, rot=90, psi=90) we would obtain the following by direct projection

But as we know that our projection plane always has got its reference with X axis facing to the right and the Y axis looking down, we should rotate the projection obtained until it is well-referenced






 In the las figures we can check the consistency of this reasoning. The projection with tilt=90, rot=90, psi=90 has been taken with SPIDER. As we can see the results are identical.

PROJECTING in a more Euler fashion

So far we have devised a method to manually compute projections using a piece of paper and a fixed phantom. But this is not the way really they are taken. In this part we will move the phantom and letting the projection plane fixed. We start again from the same initial point as before. Now we will move the internal reference axes in the reverse order (click here for further information for why this is like that).

Starting Point:





3) PSI ANGLE

The rotational angle is the first Euler angle to measure, around Z. It is a left hand angle, so if the left hand thumb is aligned with Z then positive angles are measured in the direction of the rest of the fingers. We will develop the same example as before, so the new internal reference axes are





4) TILTING ANGLE

The tilting angle is a left handed turn around Y, but the new Y!!!. Applying tilt=90 in our example we get





5) ROTATIONAL ANGLE

Finally the third angle is again a left handed turn around Z, again the new one!!!. Using rot=90 we finally get
 
 

FACTS ABOUT THE EULER MATRIX and SPIDER