Forward projection of grain maps to eta/omega/dty volumes.

[jonwright]

Example / work in progress. Hopefully useful/interesting for someone (@haixing0a ?). Not ready to merge yet, but the idea is to have a plot/figure of "observed intensity" that we can compare to "computed from model". The problem is having a 4D dataset. By looking only at hkl rings, we get to 3D (azimuth, omega, dty), and then project over one of those.

Here is code to compute eta/omega/dty data for different hkl rings from a grain map. These are plotted integrated across dty (eta/omega map is like a pole figure) or integrated over eta (omega/dty is a sinogram) :

See: fwdproj_pbp_hpc761.pdf, original at:
/data/visitor/ihma460/id11/20240621/PROCESSED_DATA/SS316/fwdproj_pbp_hpc761.ipynb

The code for plotting pixels and/or peaks should be in here:

Create_PoleFig_Volumes-ss316.pdf
/visitor/ihma460/id11/20240302/PROCESSED_DATA/JW_AZ91/Create_PoleFig_Volumes-ss316.ipynb

This "just" needs to loop over all the rotation scans. There is a version of the code using full 2D frames rather than sparse peaks (/data/visitor/ma5839/id11/20240118/SCRIPTS/Create_PoleFig_Volumes.ipynb)

Something to be cleaned up after #372 is done, as it relies on matching model to data, it needs to find both.