plt.imshow(np.abs(np.fft.fftshift(fft)), interpolation='nearest')
plt.show()
plt.imshow(np.abs(np.fft.fftshift(fft))[224:288, 224:288], interpolation=
    'nearest')
plt.show()
Grayscale image to NumPy array for Fourier transform