...
2.20
- Ellipticity components
- Make new figure with running mean lines displayed for x/y source separation over 1000 images
Elana sent me a source that breaks down ellipticity and what it means.its components mean.
I have a code block that calculates the ellipticity components e1 and e2 given the approximate coordinates of a star. I calculate ellipticity and its components e1 and e2 (stretching along x/y and y=x/y=-x). Once I calculate the components, the ellipse's major axis is at an angle a = 1/2 arctan(e2/e1) from the positive x-axis, so then I can calculate what the rotation of the focal plane was for this image, and from that back out this angle in that reference frame.
...
2.21
Just realized that although the running mean I calculated in the real notebook yesterday was correct and in the rotated frame, my plots above are NOT in the rotated plane and should be re-made. New plots (for four of the pairs):
...
- make x-axis time
- make y-axis magnitude (unsigned)
- keep same y-scale on both graphs
Want to figure out the simplest way to convert img_num to time. I think when I sample temperature I can also take the time average for each temperature item, and then plot all graphs against those times instead of image numbers. Updated image:
...
The scaling on y is difficult, but I can set the distance between ticks to be the same. However, the actual distance between them is different.
The updated figure looks like this:
...
I also have to subtract off the running mean from the other series of plots (over the course of the night). Doing that now.
That plot now looks like this:
...
Deciding whether or not to keep the final 11:51-11:59 plot in the sequence.
Now calculating the ellipticity components. We have the position angle calculated, and we calculate rotazel in the same way as before (with calculate_rotazel function). Now, we rotate the points before calculating the ellipticity components. Wondering if we can just rotate the angle?