...
- Find and kill the stray light source. But we can't go back in time and do that, so instead just cut off y<500 from the frames, to make things cleaner
- compute the median of all columns and subtract that from each column, to get rid of the horizontal streaking
- Either smooth the image with a Gaussian kernel with sigma~3 pixels to account for residual Bayer pixel artifacts and fragmentation of the PSF
- run a source detection code with a fairly high threshold of significance, and then select out the bright sources that have centroids near where we expect them.
- compute mean of all 7 centroids, call that the 'average position' of all 7. That is the origin of the center-of-spots coordinate system.
- Compute the mean position of each spot in that new coordinate system. This accounts for boresight jitter.
- For each frame and for each spot , subtract that average from each of the 7 positions, which is differential motion of each spotoff that mean-spot-position. This is the instantaneous offset of each spot. It's a 2-d vector
- For each pair of spots (i,j) in each frame, compute the vector difference r_i-r_j of their instantaneous offset vectors. Compute the norm of that quantity. If spot motions are 100% correlated, this quantity would be zero.
- Plot this magnitude-of-offset-difference vs. separation.
- I'd also compute the mean intensity for all 7 spots, and take the ratio of each spot's intensity to that mean. This will let us look at scintillation as well as image motion.
- Put the results into some sensible data structure(s)
For reasons of time I'm going to skip steps 1,2,3 and run Source Extractor on the four frames, with a tweaked configuration file that requires high significance for detection.
...
Those are the lines on the plot above.
A few other comments.
1) the assumption of isotropy does not hold. Motion in y is more than motion in x:
2) If image motions are totally uncorrelated, we expect the RMS of their separation to be sqrt(2) times the RMS of a single spot. So asymptotic value should be around 5.5 pixels or so?
3) We need to span a wider separation in order to reach asymptotic value, where relative motion becomes independent of separation, or else have smaller amounts of turbulence, or both.
MATLAB code DIMSUM.m is here: DIMSUM.m
...