The image.save_fits() function was incorrectly shaping non-square images in Q, U, and V. E.g., an image that was supposed to have a shape of (100,20) would instead be saved with a shape of (100,20) in I but (20,100) in Q, U, and V, which caused issues when reading the file. The fix was simply to swap the xdim and ydim in some np.reshape() calls inside of the function.
The image.save_fits() function was incorrectly shaping non-square images in Q, U, and V. E.g., an image that was supposed to have a shape of (100,20) would instead be saved with a shape of (100,20) in I but (20,100) in Q, U, and V, which caused issues when reading the file. The fix was simply to swap the xdim and ydim in some np.reshape() calls inside of the function.