位相限定相関のnumpy実装について

Question

前提・実現したいこと

opencvなどで実装されている位相限定相関法をnumpyで実装して、画像のズレを計算したいと思っています。
2枚の画像をnumpyで作成し、そのズレを計算する方法を実装したのですが、うまくずれが計算できません。
自分の理解が誤っている部分がどこかにあるはずなのですが、上手く行かない原因がわかりません。
ご存知の方がいらっしゃいましたら、教えていただけますと嬉しいです。
よろしくお願いいたします。

参考にしたURL: http://www.sanko-shoko.net/note.php?id=13hc

発生している問題・エラーメッセージ

以下のようなコードを作成していますが、うまく画像のシフトが計算できません。

x,y = np.meshgrid(np.arange(0,1000),np.arange(0,1000)) 
map1 = np.sqrt((x-300)**2 + (y-600)**2)
map2 = np.sqrt((x-400)**2 + (y-850)**2)
mask1 = (map1<100)
mask2 = (map2<100)
plt.imshow(mask1)
plt.show()
plt.imshow(mask2)
plt.show()

f = mask1
g = mask2

#2次元FFT処理
F = np.fft.fft2(f)
G = np.fft.fft2(g)
G_ = np.conj(f2)

R = F*G_/np.abs(F*G_)
fft_R = np.real(np.fft.ifft2(R))

plt.imshow(fft_R)
np.unravel_index(np.argmax(fft_R), array.shape)

補足情報（FW/ツールのバージョンなど）

python3.6

Answer 1

コードにミスがあるように思います。

G_ = np.conj(f2)

これを

G_ = np.conj(G)

と変更すべきかと思います。
複素共役をとるのはGですよね。