Physicist here. I have a dataset. The data is the emission from a molecule that has two dipoles. Molecules can only emit along these dipoles.
As I rotate the molecule, I will selectively excite the dipoles, thus my gaussian emission should change and become squished and more ellipsoid. And it does!
e.g squished gaussian emission
Now, I was hoping to do some PCA (principal component analysis) and I should get two principal components which would've changed as the direction of the emission changes and correspond to the dipoles. (i.e. two component eigenvalues which make up all the emission), see this image from wikipedia
However, when I've done PCA I've always done it on individual data points (think scatter plot)- "coordinate like" data points. But my data is really the total number of counts on each pixel, not individual points on a graph. The data looks like a matrix like:
(0,0,0.5,0,0) (0,0.5,1,0.5,0) (0,0.5,1,0.5,0) (0,0,0.5,0,0)
where the numbers correspond to the level of emission from that point. I could create a number datapoints which correspond to the intensity at pixel.
I feel that it is dodgy to infer the individual counts from a total count. The mapping between sums and their constituents is not bijective.
So, Can I still use PCA in this scenario? Is there a better method?