Correlation doesn't give the sign of the phase difference though as it commutes.
Taking the synchronous FFT of both signals and then dividing the complex amplitudes of the peak frequency can do this though. Division doesn't commute.
The problem is that I suspect the two signals are not sampled synchronously, but in round robin, adding 1/2 a sample time's skew between them - to be calibrated out?