The first goal is to produce diagnostic information about the offset (displacement) of the spectrum in both directions relative to a reference spectrum.
A secondary goal is to produce small Quicklook plots showing the cross-correlation function.
Try using the upsampled DFT in scipy
The analysis method should be part of the Analyze2D class.
Notes from discussion with Sam, Andrew, and BJ: