Interpret numpy.fft.fft2 output

TheChymera picture TheChymera · Jan 26, 2014 · Viewed 13.2k times · Source

My goal is to obtain a plot with the spatial frequencies of an image - kind of like doing a fourier transformation on it. I don't care about the position on the image of features with the frequency f (for instance); I'd just like to have a graphic which tells me how much of every frequency I have (the amplitude for a frequency band could be represented by the sum of contrasts with that frequency).

I am trying to do this via the numpy.fft.fft2 function.

Here is a link to a minimal example portraying my use case.

As it turns out I only get distinctly larger values for frequencies[:30,:30], and of these the absolute highest value is frequencies[0,0]. How can I interpret this?

  • What exactly does the amplitude of each value stand for?
  • What does it mean that my highest value is in frequency[0,0] What is a 0 Hz frequency?
  • Can I bin the values somehow so that my frequency spectrum is orientation agnostic?

Answer

unutbu picture unutbu · Jan 26, 2014

freq has a few very large values, and lots of small values. You can see that by plotting

plt.hist(freq.ravel(), bins=100)

(See below.) So, when you use

ax1.imshow(freq, interpolation="none")

Matplotlib uses freq.min() as the lowest value in the color range (which is by default colored blue), and freq.max() as the highest value in the color range (which is by default colored red). Since almost all the values in freq are near the blue end, the plot as a whole looks blue.

You can get a more informative plot by rescaling the values in freq so that the low values are more widely distributed on the color range.

For example, you can get a better distribution of values by taking the log of freq. (You probably don't want to throw away the highest values, since they correspond to frequencies with the highest power.)

import matplotlib as ml
import matplotlib.pyplot as plt
import numpy as np
import Image
file_path = "data"
image = np.asarray(Image.open(file_path).convert('L'))
freq = np.fft.fft2(image)
freq = np.abs(freq)

fig, ax = plt.subplots(nrows=2, ncols=2, figsize=(14, 6))
ax[0,0].hist(freq.ravel(), bins=100)
ax[0,0].set_title('hist(freq)')
ax[0,1].hist(np.log(freq).ravel(), bins=100)
ax[0,1].set_title('hist(log(freq))')
ax[1,0].imshow(np.log(freq), interpolation="none")
ax[1,0].set_title('log(freq)')
ax[1,1].imshow(image, interpolation="none")
plt.show()

enter image description here


From the docs:

The output, analogously to fft, contains the term for zero frequency in the low-order corner of the transformed axes,

Thus, freq[0,0] is the "zero frequency" term. In other words, it is the constant term in the discrete Fourier Transform.