rot vs psi

[kimdn]

@Guillawme

As I ran angdist, I see 'rot vs tilt' distribution.

Can I see psi distribution as well such as 'rot vs psi'?

Since you developed working code already, I think that it is trivial to add 'psi'?

Just add
rlnAnglePsi
to
angles = star_data['particles'][['rlnAngleRot', 'rlnAngleTilt']]

and change
hb = ax.hexbin(angles.rlnAngleRot, angles.rlnAngleTilt, bins='log', cmap=colormap, gridsize=gridsize)
similarly (as user specified an option)?

Thanks!

[Guillawme]

rlnAnglePsi is less useful, because it is the in-plane rotation that relates particle images that are projections of the same orientation in the two other directions (so, particle images that have the same or very close rlnAngleRot and rlnAngleTilt). Also, which distribution do you think would make sense to look at? rlnAnglePsi versus rlnAngleRot? Or rlnAnglePsi versus rlnAngleTilt? Or any of the other possible combinations?

I might add this option, but can't promise I will do it any time soon because I have a lot of other work and this tool is complete for my use case. But I will consider a pull request, so if you want to make the change you suggest please go ahead.

[kimdn]

I see, thank you for your explanation. I learned a lot.
I still need to study about Euler angles for cryo-EM to better answer your questions.

I just explored by using part of your code into my jupyter-lab to generate below (if I end up using your code for a publishable work, I will definitely cite!).