Hi Jenis,
Great question! freud doesn't have anything specifically for this use case, however we can modify the problem at bit so freud can still be helpful.
I think the best way to do your calculation is to project your points onto the desired subspaces and then compute the RDF for the projected 2D space. For example, to get the perpendicular RDF, project all points onto a 2D plane P which is perpendicular to the semimajor axis. Then, rotate P such that all points in P are in the XY plane. At this point, freud can compute the RDF for the 2D system.
Give this a try and let me know if its useful!
- Tommy