some optimization tests: turns out, at least on my laptop, np.dot is a serious bottleneck.  The fourier transform steps are, unfortunately, the other bottleneck - only thing to do is add more checks in correlate2d to reduce # of FTs, though I'm not sure that will help this code any.
some optimization tests: turns out, at least on my laptop, np.dot is a
serious bottleneck.  The fourier transform steps are, unfortunately, the
other bottleneck - only thing to do is add more checks in correlate2d to
reduce # of FTs, though I'm not sure that will help this code any.

Maybe np.dot -> scipy.flbas.dgemm?

alexlib/image_registration,jueqingsizhe66/image_registration,keflavich/image_registration