Linear Shift Values
I modified some lines for the array subtraction operations to get rid of shape mismatch issues. However, I'm not sure if they are really the values you wanted.
What match_drift_ref is doing currently:
- Match drift_coords (SkyCoord of the VR image coordinates) to ref_coords(SkyCoord of the main catalogue coordinates).
- Subtract each of the drift coordinates from its matched point in the catalogue and get xshift, yshift.
- Take the first values in the xshift and yshift arrays as references (?) and add them to each of the centroid coordinates -> drift_xy.